Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation codes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Operational Stability Limits in Rotating Detonation Engine Numerical Simulations

An instability is described which arises in computational fluid dynamic (CFD) simulations of semi-idealized rotating detonation engines (RDE) configured with a throat at the exit. Its existence is verified by examining output from two independently developed CFD codes simulating the same configuration and producing solutions that agree well. The instability is shown to be thermo-acoustic in that a spatial integral of the product of pressure and heat release fluctuations develops a regular oscillation which grows in time. The instability can become severe enough to cause detonation failure. Its onset is shown to be closely linked to the size of the exit throat and the size of the inlet restriction; both parameters that strongly influence RDE performance. It is shown that the instability places a cap on ideal RDEperformance, but that an optimized exhaust throat and inlet restriction combination still yields substantial pressure gain. Other parametric sensitivities are also examined in terms ofinstability growth. These include axial length, inlet manifold pressure, and air-fuel ratio.

Paxson, Daniel E.↗
DOE OSTI2024

Charge lifetime improvement of the Continuous Electron Beam Accelerator Facility photogun with a biased anode

GaAs photocathodes in dc high-voltage photoguns are highly susceptible to ion back-bombardment, which reduces the photocathode quantum efficiency and limits the useful operating lifetime for producing polarized electron beams. This paper demonstrates that applying a small positive bias to the photogun anode can significantly suppress ion back-bombardment and increase charge lifetime. This technique was studied extensively using the Continuous Electron Beam Accelerator Facility photogun, where highly polarized electron beams created using a strained-superlattice GaAs / GaAsP photocathode were used and charge lifetimes improved by almost a factor of 2. A new simulation code was developed to model ion production and tracking in order to better understand and explain the factors that led to the performance improvement. Results of the experiments and simulations are discussed in detail. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
NASA NTRS1990

Reduced order component models for flexible multibody dynamics simulations

Many flexible multibody dynamics simulation codes require some form of component description that properly characterizes the dynamic behavior of the system. A model reduction procedure for producing low order component models for flexible multibody simulation is described. Referred to as projection and assembly, the method is a Rayleigh-Ritz approach that uses partitions of the system modal matrix as component Ritz transformation matrices. It is shown that the projection and assembly method yields a reduced system model that preserves a specified set of the full order system modes. Unlike classical component mode synthesis methods, the exactness of the method described is obtained at the expense of having to compute the full order system modes. The paper provides a comprehensive description of the method, a proof of exactness, and numerical results demonstrating the method's effectiveness.

Tsuha, Walter S.↗
DOE OSTI2026

Core plasma fueling by fast inward particle transport after hydrogen pellet injection in Wendelstein 7-X

A large database of more than 1000 individual cryogenic hydrogen pellets injected into Wendelstein 7-X for plasma fueling was analyzed to improve the understanding of the three phases of the process: the ablation, deposition and transport of the pellet material. Kilohertz-sampled electron density and temperature measurements revealed a more complex drift behavior than predicted by numerical code simulation. It could be explained by the poloidal plasma E r x B- drift rotation, which plays a significant role in stellarators, but was not previously considered in pellet injection codes like HPI2. The drift results in a fast poloidal rotation of the pellet material around the plasma core, leading to an almost homogeneous deposition over the involved flux surfaces regardless of magnetic high and low field side injection geometry. Additionally, a novel fast inward directed transport mechanism (‘FIT-effect’) was observed. The effect occurs on timescales of tens of milliseconds and cannot be explained by neoclassical transport or diffusion. It might be linked to the turbulence pinch recently found in Wendelstein 7-X. When the FIT-effect occurs, the pellet particles are rapidly transferred from the deposition flux surfaces to the plasma core, causing the plasma density profile to peak, which is beneficial for confinement in Wendelstein 7-X. The large pellet injection database was statistical analyzed with regard to pellet and plasma parameters, which delivered some starting points towards developing an understanding of the physics behind the FIT-effect. The results indicate, that plasma core fueling via pellet injection is largely independent of the injection geometry in stellarators under certain conditions, reducing the technical complexity of the injection system.

Wendelstein 7-X↗
DOE OSTI2024

Semiclassical treatment of photon cascades in nuclei

Here, we present a simple semiclassical treatment of photon cascades, suitable for use in nuclear fission simulation codes. The approximation is here developed for E⁢1 and E⁢2 transitions and its quality is illustrated for a variety of two-photon cascades. Implementation of the treatment into Monte Carlo simulations would make it possible to address photon correlation observables quantitatively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multiphysics modeling of accelerators through code integration

This work aims to improve the ability of particle accelerator researchers to develop high-performance accelerator cavity designs by creating an overall multiphysics framework that integrates and couples existing application codes. This framework will allow accelerator researchers to build multiphysics models that will optimize cavity design, improve understanding of whole-device performance, and reduce the development and fabrication costs of accelerator research. We utilize the open-source VizSchema data standard as an intermediate data structure interface layer to standardize interfaces between individual application codes. VizScema is extensively documented online, and plugins for VizSchema are available for popular visualization packages, including VisIt and ParaView. Currently, the work focuses on coupling the EM field solver COMSOL and the electron gun code MICHELLE to allow COMSOL field-solve results to be seamlessly used by MICHELLE for particle-solve. Later work will extend this integration to include other fields, particles, and thermodynamics simulation codes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI1995

Nuclear Safety [Vol. 36, No. 2, July-December 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 195 The Chornobyl Accident Revisited, Part III: Chernobyl Source Term Release Dynamics and Reconstruction of Events During the Active Phase, A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 218 Second ANS Workshop on the Safety of Soviet-Designed Nuclear Power Plants, R. A. Bari; 234 Elements of a Nuclear Criticality Safety Program, C. M. Hopper; 243 Rickover, Excellence, and Criticality Safety Programs, R. E. Wilson; ACCIDENT ANALYSIS: 249 Transient Analysis of the PIUS Advanced Reactor Design with the TRAC-PF1/MOD2 Code, B. E. Boyack, J. L Steiner, S. C. Harmony, H. J. Stumpf, and J. F. Lime; 278 The Hierarchy-By-Interval Approach to Identifying Important Models that Need Improvement in Severe-Accident Simulation Codes, T. J. Heames, M. Khatib-Rahbar, J. E. Kelly, R. P. Jenks-Johnson, and Y.-S. Chen; 290 RELAP5/MOD3 Code Coupling Model, R. P. Martin; 299 Missiles Caused by Severe Pressurized-Water Reactor Accidents, R. Krieg; DESIGN FEATURES: 310 Validation of COMMIX with Westinghouse AP-600 PCCS Test Data, J. G. Sun, T. H. Chien, J. Ding, and W T. Sha; ENVIRONMENTAL EFFECTS: 321 Spent Nuclear Fuel Characterization for a Bounding Reference Assembly for the Receiving Basin for Off-Site Fuel, S. D. Kahook, R. L. Garrett, L. R. Canas, and M J. Beckum; OPERATING EXPERIENCES: 332 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 335 Reactor Coolant System Blowdown at Wolf Creek on September 17, 1994, J. V. Kauffman and S. L. Israel; RECENT DEVELOPMENTS: 344 Reports, Standards, and Safety Guides, D. S. Oueener; 349 Proposed Rule Changes as of June 30, 1995; ANNOUNCEMENTS: 320 Symposium on Acceptability of Risk From Radiation—Application to Manned Space Flight; 320 24th DOE/NRC Nuclear Air Cleaning and Treatment Conference; 361 Radiation Biology and Radiation Protection— Modern Developments and Tendencies in Radiation Biology; 362 1997 IEEE Sixth Conference on Human Factors and Power Plants; 354 The Authors; 360 Reviewers of Nuclear Safety, Vol. 36.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Prediction of neutron production and energy spectrum by the inverse kinematic reaction between an incident 7 Li 3+ beam and a proton target in PHITS

A neutron source using the inverse kinematic reaction between lithium and proton, p( 7 Li, n) 7 Be, achieves forward-directed neutrons, potentially enhancing neutron yield in the forward direction. Despite the advantage, no evaluated-cross-section data for this reaction can be used in Monte Carlo simulation codes, such as PHITS. To solve this problem, this study aims to evaluate the applicability of the user-defined cross-section data, Frag data, for p ( 7 Li, n) 7 Be in PHITS. The simulations reproduced collisions between 7 Li 3+ ions and polypropylene targets. The Frag data was edited based on the JENDL-5 by utilizing the two-body collision kinematics. The neutron yield and angular distribution were investigated in the simulation. As a result, the forward neutron convergence with a reasonable neutron yield and energy spectrum was observed. The expected neutron yield in the forward 1-steradian area is 2.46 × 10 10 n/s when lithium-ion energy and current are 16.45 MeV and 0.1 mA.

43 PARTICLE ACCELERATORS↗
NASA NTRS1989

Simulation of electrostatic ion instabilities in the presence of parallel currents and transverse electric fields

A spatially two-dimensional electrostatic PIC simulation code was used to study the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a field-aligned drift for L is much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field aligned drift and the E x B drift are individually subcritical. The simulation results show that the growing ion waves are associated with small vortices in the linear stage, which evolve to the nonlinear stage dominated by larger vortices with lower frequencies.

Nishikawa, K.-I.↗
NASA NTRS1990

Simulation of electrostatic turbulence due to sheared flows parallel and transverse to the magnetic field

A spatially two-dimensional electrostatic particle simulation code is used to examine the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a magnetic-field-aligned electron drift for L much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field-aligned drift and the E x B drift are individually subcritical. The simulation results indicate that a broadband turbulence is associated with such an equilibrium.

Nishikawa, K.-I.↗
NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
NASA NTRS1991

Simulating The Performance Of A Scramjet

SCRAM computer program determines, according to one-dimensional mathematical model, cycle performance for airframe-integrated subsonic or supersonic combustion Ramjets. Reliable, efficient, and speedy ramjet designer's software tool usable on all standard computers down to those compatible with IBM PC-AT's. Scramjet cycle analysis code simulates, from nose to tail, hydrogen-fueled, airframe-integrated scramjet in real gas flow with equilibrium thermodynamic properties. Enables rapid scramjet performance estimates. Written in FORTRAN 77.

Walton, James T.↗
NASA NTRS1994

Research on Orbital Plasma-Electrodynamics (ROPE)

Since the development of probe theory by Langmuir and Blodgett, the problem of current collection by a charged spherically or cylindrically symmetric body has been investigated by a number of authors. This paper overviews the development of a fully three-dimensional particle simulation code which can be used to understand the physics of current collection in three dimensions and can be used to analyze data resulting from the future tethered satellite system (TSS). According to the TSS configurations, two types of particle simulation models were constructed: a simple particle simulation (SIPS) and a super particle simulation (SUPS). The models study the electron transient response and its asymptotic behavior around a three dimensional, highly biased satellite. The potential distribution surrounding the satellite is determined by solving Laplace's equation in the SIPS model and by solving Poisson's equation in the SUPS model. Thus, the potential distribution in space is independent of the density distribution of the particles in the SUPS model but it does depend on the density distribution of the particles in the SUPS model. The evolution of the potential distribution in the SUPS model is described. When the spherical satellite is charged to a highly positive potential and immersed in a plasma with a uniform magnetic field, the formation of an electron torus in the equatorial plane (the plane in perpendicular to the magnetic field) and elongation of the torus along the magnetic field are found in both the SIPS and the SUPS models but the shape of the torus is different. The areas of high potential that exist in the polar regions in the SUPS model exaggerate the elongation of the electron torus along the magnetic field. The current collected by the satellite for different magentic field strengths is investigated in both models. Due to the nonlinear effects present in SUPS, the oscillating phenomenon of the current collection curve during the first 10 plasma periods can be seen (this does not appear in SIPS). From the parametric studies, it appears that the oscillating phenomenon of the current collection curve occurs only when the magnetic field strength is less than 0.2 gauss for the present model.

Wu, S. T.↗
NASA NTRS2001

A Comparison of Three-Dimensional Simulations of Traveling-Wave Tube Cold-Test Characteristics Using CST MICROWAVE STUDIO and MAFIA

Previously, it was shown that MAFIA (solutions of Maxwell's equations by the Finite Integration Algorithm), a three-dimensional simulation code, can be used to produce accurate cold-test characteristics including frequency-phase dispersion, interaction impedance, and attenuation for traveling-wave tube (TWT) slow-wave structures. In an effort to improve user-friendliness and simulation time, a model was developed to compute the cold-test parameters using the electromagnetic field simulation software package CST MICROWAVE STUDIO (MWS). Cold-test parameters were calculated for several slow-wave circuits including a ferruled coupled-cavity, a folded waveguide, and a novel finned-ladder circuit using both MWS and MAFIA. Comparisons indicate that MWS provides more accurate cold-test data with significantly reduced simulation times. Both MAFIA and MWS are based on the finite integration (FI) method; however, MWS has several advantages over MAFIA. First, it has a Windows based interface for PC operation, making it very user-friendly, whereas MAFIA is UNIX based. MWS uses a new Perfect Boundary Approximation (PBA), which increases the accuracy of the simulations by avoiding stair step approximations associated with MAFIA's representation of structures. Finally, MWS includes a Visual Basic for Applications (VBA) compatible macro language that enables the simulation process to be automated and allows for the optimization of user-defined goal functions, such as interaction impedance.

Chevalier, C. T.↗
NASA NTRS1994

Axial focusing of impact energy in the Earth's interior: Proof-of-principle tests of a new hypothesis

A causal link between major impact events and global processes would probably require a significant change in the thermal state of the Earth's interior, presumably brought about by coupling of impact energy. One possible mechanism for such energy coupling from the surface to the deep interior would be through focusing due to axial symmetry. Antipodal focusing of surface and body waves from earthquakes is a well-known phenomenon which has previously been exploited by seismologists in studies of the Earth's deep interior. Antipodal focusing from impacts on the Moon, Mercury, and icy satellites has also been invoked by planetary scientists to explain unusual surface features opposite some of the large impact structures on these bodies. For example, 'disrupted' terrains have been observed antipodal to the Caloris impact basis on Mercury and Imbrium Basin on the Moon. Very recently there have been speculations that antipodal focusing of impact energy within the mantle may lead to flood basalt and hotspot activity, but there has not yet been an attempt at a rigorous model. A new hypothesis was proposed and preliminary proof-of-principle tests for the coupling of energy from major impacts to the mantle by axial focusing of seismic waves was performed. Because of the axial symmetry of the explosive source, the phases and amplitudes are dependent only on ray parameter (or takeoff angle) and are independent of azimuthal angle. For a symmetric and homogeneous Earth, all the seismic energy radiated by the impact at a given takeoff angle will be refocused (minus attenuation) on the axis of symmetry, regardless of the number of reflections and refractions it has experienced. Mantle material near the axis of symmetry will experience more strain cycles with much greater amplitude than elsewhere and will therefore experience more irreversible heating. The situation is very different than for a giant earthquake, which in addition to having less energy, has an asymmetric focal mechanism and a larger area. Two independent proof-of-principle approaches were used. The first makes use of seismic simulations, which are being performed with a realistic Earth model to determine the degree of focusing along the axis and to estimate the volume of material, if any, that experiences significant irreversible heating. The second involves two-dimensional hydrodynamic code simulations to determine the stress history, internal energy, and temperature rise as a function of radius along the axis.

Boslough, M. B.↗
NASA NTRS1999

A Mixed-Mode I/II Fracture Criterion and Its Application in Crack Growth Predictions

A crack tip opening displacement (CTOD)-based, mixed mode fracture criterion is developed for predicting the onset and direction of crack growth. The criterion postulates that crack growth occurs in either the Mode I or Mode II direction, depending on whether the maximum in either the opening or the shear component of CTOD, measured at a specified distance behind the crack tip, attains a critical value. For crack growth direction prediction, the proposed CTOD criterion is shown to be equivalent to seven commonly used crack growth criteria under linearly elastic and asymptotic conditions. Under elastic-plastic conditions the CTOD criterion's prediction of the dependence of the crack growth direction on the crack-up mode mixity is in excellent agreement with the Arcan test results. Furthermore, the CTOD criterion correctly predicts the existence of a crack growth transition from mode I to mode II as the mode mixity approaches the mode II loading condition. The proposed CTOD criterion has been implemented in finite element crack growth simulation codes Z1P2DL and FRANC2DL to predict the crack growth paths in (a) a modified Arcan test specimen and fixture made of AL 2024-T34 and (b) a double cantilever beam (DCB) specimen made of AL 7050. A series of crack growth simulations have been carried out for the crack growth tests in the Arcan and DCB specimens and the results further demonstrate the applicability of the mixed mode CTOD fracture criterion crack growth predictions and residual strength analyses for airframe materials.

Sutton, Michael A.↗