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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Transfer functions for Q A /Q B international regulatory limits for the safe transport of radioactive materials

This paper presents a proposed revision of the International Atomic Energy Agency transport regulations, related to the A 1 and A 2 limit values used to determine the radioactive transport classification. Based on the 'Q system', a novel methodology was introduced to derive Q A and Q B values related to scenarios involving external exposure from a distant source. These values are key parameters that respectively represent the total effective dose and total equivalent dose to the skin, from all primary and secondary particles contributing to radiation exposure. The International Working Group (WG A 1 /A 2 ) is established and associated with the TRANSSC Technical Expert Group on Radiation Protection. A review of the A 1 and A 2 values is performed in response to identified limitations within the existing Q system. The followed approach is based on Monte Carlo simulations that enabled the development of transfer functions aimed at reducing computational time and increasing the flexibility of dose evaluations for any radionuclide with known particle emission spectra. This method allows updating the Q A and Q B values to account for future data evolutions (decay data, fluence-to-dose conversion coefficients) and standardizing the calculation of regulation limits across all referenced radionuclides and scenarios related to external exposure. The transfer functions are established using three Monte Carlo simulation codes—FLUKA, Geant4, and MCNP—and address the previous limitations of the 'Q system', reflecting the latest International Commission for Radiation Protection recommendations and improvements in calculation techniques. The results of the WG show consistent agreement across the codes, with minor discrepancies observed at low primary energies due to statistical uncertainties and different handling of stopping power for electrons/positrons in the codes. This revised approach aligns with current standards and recommendations, ensuring that the radiological consequences of transport accidents are acceptable for the new A 1 and A 2 limits from a radiological protection perspective.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1995

Topological structure of the magnetotail as a function of interplanetary magnetic field direction

Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Fedder, J. A.↗
DOE OSTI2025

Benchmarking MELCOR's NAC Package to ABCOVE Test AB7

This report presents analyses of the AB7 ABCOVE sodium spray fire experiment with the MELCOR code. This code simulates the progression of accident events for analysis and auditing purposes of nuclear facilities during accident conditions. Historically, the ABCOVE experiments have contributed to the validation of aerosol physics and related phenomena. Given advancements in sodium-cooled reactor designs, characterization of the sodium spray combustion may further the review and validation of newly incorporated sodium properties and physics packages, namely, the sodium equations of state (EOS) and the sodium combustion (NAC) package within MELCOR. Previously, the AB5 and AB6 experiments were analyzed with and without the NAC package. This work builds on the previous analyses with a demonstration of the current code capabilities of MELCOR with a more mild Na spray and pool fire scenario.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Democratizing uncertainty quantification

Uncertainty Quantification (UQ) is vital to safety-critical model-based analyses, but the widespread adoption of sophisticated UQ methods is limited by technical complexity. In this paper, we introduce UM-Bridge (the UQ and Modeling Bridge), a high-level abstraction and software protocol that facilitates universal interoperability of UQ software with simulation codes. It breaks down the technical complexity of advanced UQ applications and enables separation of concerns between experts. UM-Bridge democratizes UQ by allowing effective interdisciplinary collaboration, accelerating the development of advanced UQ methods, and making it easy to perform UQ analyses from prototype to High Performance Computing (HPC) scale. In addition, we present a library of ready-to-run UQ benchmark problems, all easily accessible through UM-Bridge. These benchmarks support UQ methodology research, enabling reproducible performance comparisons. We demonstrate UM-Bridge with several scientific applications, harnessing HPC resources even using UQ codes not designed with HPC support.

Benchmarks↗
NASA NTRS2021

Improvements to the Flight Analysis and Simulation Tool (FAST) and Initial Development of the Genesis Flight Mechanics Simulation for Ascent, Aerocapture, Entry, Descent, and Landing (A2EDL) Trajectory Design

The NASA Engineering and Safety Center (NESC) was requested to provide significant enhancements to the Flight Analysis and Simulation Tool (FAST), a generic, variable-degree-of-freedom, multi-body ascent, aerocapture, entry, descent, and landing (A2EDL) flight simulation code and a key Agency analysis tool. This report details the upgrades to FAST that were completed, describes a change to the assessment’s scope that allowed its goals to be exceeded with a new and innovative approach, and discusses the results of the modified approach.

NASA Engineering and Safety Center; Flight Analysi↗
NASA NTRS2011

Numerical Study Comparing RANS and LES Approaches on a Circulation Control Airfoil

A numerical study over a nominally two-dimensional circulation control airfoil is performed using a large-eddy simulation code and two Reynolds-averaged Navier-Stokes codes. Different Coanda jet blowing conditions are investigated. In addition to investigating the influence of grid density, a comparison is made between incompressible and compressible flow solvers. The incompressible equations are found to yield negligible differences from the compressible equations up to at least a jet exit Mach number of 0.64. The effects of different turbulence models are also studied. Models that do not account for streamline curvature effects tend to predict jet separation from the Coanda surface too late, and can produce non-physical solutions at high blowing rates. Three different turbulence models that account for streamline curvature are compared with each other and with large eddy simulation solutions. All three models are found to predict the Coanda jet separation location reasonably well, but one of the models predicts specific flow field details near the Coanda surface prior to separation much better than the other two. All Reynolds-averaged Navier-Stokes computations produce higher circulation than large eddy simulation computations, with different stagnation point location and greater flow acceleration around the nose onto the upper surface. The precise reasons for the higher circulation are not clear, although it is not solely a function of predicting the jet separation location correctly.

Rumsey, Christopher L.↗
NASA NTRS1993

Hybrid simulations of thin current sheets

A one-dimensional, hybrid simulation code is used to study current sheets with a nonzero normal magnetic field B(sub z) and a dawn-to-dusk electric field E(sub y). Such configurations are dependent upon only two parameters: we use the normalized normal magnetic field B-normalized (sub z) = B(sub z)/(4(pi)(n(sub b)) (v(exp 2 sub T))(exp 1/2) and normalized electric field V-normalized (sub D) = (1/V(sub T)(cE(sub y)/B(sub z)), where V(sub T) is the thermal velocity of ions prior to their interaction with the current sheet and n(sub b) is the number density outside the current sheet (at the simulation boundary). A third parameter that is relevant to the motion of particles in current sheets is kappa(sub A), the value of kappa = (R(sub min)/rho(sub max))(exp 1/2) for particles of average energy. We find that if either B-normalized (sub z) is close to or greater than 1, or if kappa(sub A) is close to 1, a rotational mode develops in which the z = 0 current rotates with the ion sense about the normal magnetic field, while for small values of both B-normalized (sub z) or kappa(sub A), the configuration is quasi-steady. To achieve values of kappa(sub A) of the order of or larger than 1, we decrease the value of V-normalized (sub D) uniformly. We find that the magnetic field fluctuations and particle distribution functions are similar in many respects to what was observed in the day 240, 1986, Active Magnetospheric Particle Tracer Explorer (AMPTE)/CCE current disruption event, an event that appears to be located at the site of initiation of current disruption and related particle energization.

Burkhart, G. R.↗
NASA NTRS1997

Numerical Assessment of Four-Port Through-Flow Wave Rotor Cycles with Passage Height Variation

The potential for improved performance of wave rotor cycles through the use of passage height variation is examined. A Quasi-one-dimensional CFD code with experimentally validated loss models is used to determine the flowfield in the wave rotor passages. Results indicate that a carefully chosen passage height profile can produce substantial performance gains. Numerical performance data are presented for a specific profile, in a four-port, through-flow cycle design which yielded a computed 4.6% increase in design point pressure ratio over a comparably sized rotor with constant passage height. In a small gas turbine topping cycle application, this increased pressure ratio would reduce specific fuel consumption to 22% below the un-topped engine; a significant improvement over the already impressive 18% reductions predicted for the constant passage height rotor. The simulation code is briefly described. The method used to obtain rotor passage height profiles with enhanced performance is presented. Design and off-design results are shown using two different computational techniques. The paper concludes with some recommendations for further work.

Paxson, D. E.↗
DOE OSTI2024

Charged particle transport coefficient challenges in high energy density plasmas

High energy density physics (HEDP) and inertial confinement fusion (ICF) research typically relies on computational modeling using radiation-hydrodynamics codes in order to design experiments and understand their results. These tools, in turn, rely on numerous charged particle transport and relaxation coefficients to account for laser energy absorption, viscous dissipation, mass transport, thermal conduction, electrical conduction, non-local ion (including charged fusion product) transport, non-local electron transport, magnetohydrodynamics, multi-ion-species thermalization, and electron-ion equilibration. In many situations, these coefficients couple to other physics, such as imposed or self-generated magnetic fields. Furthermore, how these coefficients combine are sensitive to plasma conditions as well as how materials are distributed within a computational cell. Uncertainties in these coefficients and how they couple to other physics could explain many of the discrepancies between simulation predictions and experimental results that persist in even the most detailed calculations. This paper reviews the challenges faced by radiation-hydrodynamics in predicting the results of HEDP and ICF experiments with regard to these and other physics models typically included in simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Mathematical description of complex chemical kinetics and application to CFD modeling codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗
NASA NTRS1994

Novel high-gain, improved-bandwidth, finned-ladder V-band Traveling-Wave Tube slow-wave circuit design

The V-band frequency range of 59-64 GHz is a region of the millimeter-wave spectrum that has been designated for inter-satellite communications. As a first effort to develop a high-efficiency V-band Traveling-Wave Tube (TWT), variations on a ring-plane slow-wave circuit were computationally investigated to develop an alternative to the more conventional ferruled coupled-cavity circuit. The ring-plane circuit was chosen because of its high interaction impedance, large beam aperture, and excellent thermal dissipation properties. Despite these advantages, however, low bandwidth and high voltage requirements have, until now, prevented its acceptance outside the laboratory. In this paper, the three-dimensional electrodynamic simulation code MAFIA (solution of MAxwell's Equation by the Finite-Integration-Algorithm) is used to investigate methods of increasing the bandwidth and lowering the operating voltage of the ring-plane circuit. Calculations of frequency-phase dispersion, beam on-axis interaction impedance, attenuation and small-signal gain per wavelength were performed for various geometric variations and loading distributions of the ring-plane TWT slow-wave circuit. Based on the results of the variations, a circuit termed the finned-ladder TWT slow-wave circuit was designed and is compared here to the scaled prototype ring-plane and a conventional ferruled coupled-cavity TWT circuit over the V-band frequency range. The simulation results indicate that this circuit has a much higher gain, significantly wider bandwidth, and a much lower voltage requirement than the scaled ring-plane prototype circuit, while retaining its excellent thermal dissipation properties. The finned-ladder circuit has a much larger small-signal gain per wavelength than the ferruled coupled-cavity circuit, but with a moderate sacrifice in bandwidth.

Kory, Carol L.↗
NASA NTRS1993

Mathematical Description of Complex Chemical Kinetics and Application to CFD Modeling Codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗
NASA NTRS2010

Calculation of Dose Deposition in 3D Voxels by Heavy Ions

The biological response to high-LET radiation is very different from low-LET radiation, and can be partly attributed to the energy deposition by the radiation. Several experiments, notably detection of gamma-H2AX foci by immunofluorescence, has revealed important differences in the nature and in the spatial distribution of double-strand breaks (DSB) induced by low- and high-LET radiations. Many calculations, most of which are based on amorphous track models with radial dose, have been combined with chromosome models to calculate the number and distribution of DSB within nuclei and chromosome aberrations. In this work, the Monte-Carlo track structure simulation code RITRACKS have been used to calculate directly the energy deposition in voxels (3D pixels). A cubic volume of 5 micrometers of side was irradiated by 1) 450 (1)H+ ions of 300 MeV (LET is approximately 0.3 keV/micrometer) and 2) by 1 (56)Fe26+ ion of 1 GeV/amu (LET is approximately 150 keV/micrometer). In both cases, the dose deposited in the volume is approximately 1 Gy. All energy deposition events are recorded and dose is calculated in voxels of 20 micrometers of side. The voxels are then visualized in 3D by using a color scale to represent the intensity of the dose in a voxel. This simple approach has revealed several important points which may help understand experimental observations. In both simulations, voxels which receive low dose are the most numerous, and those corresponding to electron track ends received a dose which is in the higher range. The dose voxels are distributed randomly and scattered uniformly within the volume irradiated by low-LET radiation. The distribution of the voxels shows major differences for the (56)Fe26+ ion. The track structure can still be seen, and voxels with much higher dose are found in the region corresponding to the track "core". These high-dose voxels are not found in the low-LET irradiation simulation and may be responsible for DSB that are more difficult to repair. By applying a threshold on the dose visualization, voxels corresponding to electron track ends are evidenced and the spatial distribution of voxels is very similar to the distribution of DSB observed in gamma H2AX experiments, even if no chromosomes have been included in the simulation. Furthermore, this work has shown that a significant dose is deposited in voxels corresponding to electron track ends. Since some delta-rays from iron ion can travel several millimeters, they may also be of radiobiological importance.

Plante, Ianik↗
NASA NTRS2003

A Nonlinear Model for Fuel Atomization in Spray Combustion

Most gas turbine combustion codes rely on ad-hoc statistical assumptions regarding the outcome of fuel atomization processes. The modeling effort proposed in this project is aimed at developing a realistic model to produce accurate predictions of fuel atomization parameters. The model involves application of the nonlinear stability theory to analyze the instability and subsequent disintegration of the liquid fuel sheet that is produced by fuel injection nozzles in gas turbine combustors. The fuel sheet is atomized into a multiplicity of small drops of large surface area to volume ratio to enhance the evaporation rate and combustion performance. The proposed model will effect predictions of fuel sheet atomization parameters such as drop size, velocity, and orientation as well as sheet penetration depth, breakup time and thickness. These parameters are essential for combustion simulation codes to perform a controlled and optimized design of gas turbine fuel injectors. Optimizing fuel injection processes is crucial to improving combustion efficiency and hence reducing fuel consumption and pollutants emissions.

Liu, Nan-Suey↗
NASA NTRS1991

Evolution of magnetic flux ropes associated with flux transfer events and interplanetary magnetic clouds

Spacecraft observations suggest that flux transfer events and interplanetary magnetic clouds may be associated with magnetic flux ropes which are magnetic flux tubes containing helical magnetic field lines. In the magnetic flux ropes, the azimuthal magnetic field is superposed on the axial field. The time evolution of a localized magnetic flux rope is studied. A two-dimensional compressible MHD simulation code with a cylindrical symmetry is developed to study the wave modes associated with the evolution of flux ropes. It is found that in the initial phase both the fast magnetosonic wave and the Alfven wave are developed in the flux rope. After this initial phase, the Alfven wave becomes the dominant wave mode for the evolution of the magnetic flux rope and the radial expansion velocity of the flux rope is found to be negligible. Numerical results further show that even for a large initial azimuthal component of the magnetic field, the propagation velocity along the axial direction of the flux rope remains the Alfven velocity. It is also found that the localized magnetic flux rope tends to evolve into two separate magnetic ropes propagating in opposite directions. The simulation results are used to study the evolution of magnetic flux ropes associated with flux transfer events observed at the earth's dayside magnetopause and magnetic clouds in the interplanetary space.

Wei, C. Q.↗
NASA NTRS2000

Unsteady Aerodynamic Effects on the Flight Characteristics of an F-16XL Configuration

Unsteady aerodynamic models based on windtunnel forced oscillation test data and analyzed with a fuzzy logic algorithm arc incorporated into an F-16XL flight simulation code. The reduced frequency needed in the unsteady models is numerically calculated by using a limited prior time history of state variables in a least-square sense. Numerical examples arc presented to show the accuracy of the calculated reduced frequency. Oscillatory control inputs are employed to demonstrate the differences in the flight characteristics based on unsteady and quasi-steady aerodynamic models. Application of the unsteady aerodynamic models is also presented and the results are compared with one set of F16XIL longitudinal maneuver flight data. It is shown that the main differences in dynamic response are in the lateral-directional characteristics, with the quasi-steady model being more stable than the flight vehicle, while the unsteady model being more unstable. Similar conclusions can also be made in a simulated rapid sideslipping roll. To improve unsteady aerodynamic modeling, it is recommended to acquire test data with coupled motions in pitch, roll and yaw.

Wang, Zhongjun↗
NASA NTRS2016

A New Code SORD for Simulation of Polarized Light Scattering in the Earth Atmosphere

We report a new publicly available radiative transfer (RT) code for numerical simulation of polarized light scattering in plane-parallel atmosphere of the Earth. Using 44 benchmark tests, we prove high accuracy of the new RT code, SORD (Successive ORDers of scattering). We describe capabilities of SORD and show run time for each test on two different machines. At present, SORD is supposed to work as part of the Aerosol Robotic NETwork (AERONET) inversion algorithm. For natural integration with the AERONET software, SORD is coded in Fortran 90/95. The code is available by email request from the corresponding (first) author or from ftp://climate1.gsfc.nasa.gov/skorkin/SORD/.

AERONET↗
DOE OSTI2026

An Exploration of the Equation of State Dependence of Core-collapse Supernova Explosion Outcomes and Signatures

We explore, using a state-of-the-art simulation code in 3D and late-enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state (EOS). Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear EOSs for a 9 M ⊙ /solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino “light” curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multiprogenitor set of 3D supernova simulations for various EOSs to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest toward a fully quantitative theory of supernova observables.

Rusakov, Aleksandr [Princeton University, NJ (Unit↗