Search NASASearch

SEARCH · Search NASA

Results for “numerical simulations”

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.

524 records · Page 6

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics

Coupled, Transient Simulations of the Interaction Between Power and Secondary Flowpaths in Gas Turbines

This paper presents a coupled analysis of the interaction between mainpath and secondary flowpaths in gas turbines using transient simulations. Some of the topics include: 1) Need for Coupled Analysis; 2) Primary-Secondary Coupling Schematic; 3) Secondary Flow Requirement; 4) Objectives of Present Methodology; 5) Current Methodologies Recap; 6) Proposed Coupled Code Methodology; 7) Description of SCISEAL Code; 8) Description of Turbo Code; 9) Code Coupling/Interface Issues; and 10) Current Interface Strategy. This paper is presented in viewgraph form.

M M Athavale

High-Resolution Global Ocean Simulation in the SWOT ERA: Introducing ECCO Llc4320v2

The Surface Water and Ocean Topography (SWOT) mission is opening up new opportunities to evaluate and drive improvement in km-scale global ocean and atmosphere-ocean models via its simultaneous swath observations of sea-surface height, wind speed, and significant wave height down to 2 km scales globally. Motivated by this, we have initiated a collaborative effort at NASA Ames to run a new global ocean-sea-ice simulation in MITgcm at a nominal 1/48-deg resolution of the year 2023 that enables direct comparisons with SWOT observations. New features relative to a prior simulation at a similar resolution of the year 2012 include much more realistic tidal forcing, thermodynamic coupling in ice-shelf cavities, high-frequency hourly atmospheric forcing, much improved river runoff, double the vertical resolution in the deep ocean (173 vs 90 levels), updated bathymetry including improvements in finescale bathymetry derived from SWOT, and updated vertical mixing scheme parameters that improve the equatorial Pacific circulation.

ECCO

Mechanical, Electrochemical & Thermal Modeling of Electric Vehicle Batteries for Crash Simulation (CRADA CRD-19-00811 Final Report)

Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Laboratory of the Rockies (NLR) will cooperate with HMC to develop mathematical models for battery cells and modules for simulating abuse response in batteries subject to the type of mechanical crushing that can occur in a full motor vehicle crash.

33 ADVANCED PROPULSION SYSTEMS

Optical characterization of solar concentrator mirrors using the solar beam of the JPL 25-foot space simulator

The Jet Propulsion Laboratory 25-foot space simulator with its 5.8-m (19-ft) diameter simulated solar beam provides an excellent facility for measuring the optical characteristics of parabolic solar concentrator panels and gores. The virtual source position and size were determined by using a single lamp of the 37 xenon 30-kW source array with only the center lens in the 19-channel optical mixer. This data was used to define the optical test geometry, and it allowed accurate measurement of focal length and surface deviations of the mirror under test. A flux distribution of a typical solar concentrator placed directly on the solar beam gives measurements of performance at the focal point of the parabolic surface.

Maurice J Argoud

Validation of Accurate Cryogenic Fluid Vapor-Liquid Boundary Conditions Via Molecular Simulations

Vapor-liquid interfaces drive many important phenomena in cryogenic fluid management, including heat transfer, evaporation, and capillary flow. Design of cryogenic fluid systems, such as propellant storage, requires accurate predictions of fluid behavior, including evaporation rates. Many models have been proposed for heat and mass transfer at vapor liquid interfaces, but the accuracy of these models in the context of cryogenic fluids has not been performed. We use molecular dynamics simulations, which allow for nanometer scale resolution of fluid phenomena, to evaluate the accuracy of a variety of vapor-liquid boundary conditions at evaporating and condensing interfaces. We find that an anisotropic temperature distribution is a critical ingredient for accurate prediction of intensive evaporation and condensation.

Daniel Vigil

High Temperature Dielectric Properties and Differential Scanning Calorimetry of Lunar Simulants

To guide development of microwave process technology that could be used during in situ construction on the Moon, we measured the high-temperature basic dielectric properties (εʹ and εʺ) of 17 lunar simulants and related materials. In order to confidently use these data one needs to understand the data’s strengths and weaknesses. Therefore, a goal of this publication is to provide insights into the comparative effects of sample composition, pre-treatments, experimental variables, high temperatures, and other factors on the measured response. The dielectric measurements were performed using the cavity perturbation method over a temperature range between room temperature to 1000 °C, or higher, and provided the real and imaginary components of permittivity at six frequencies. The utility of the original values was limited by the varying density of the pellets used in the measurement. Therefore, all of the εʹ and εʺ measurements at the frequency of 2466 MHz have been scaled to a constant density, 1.75 g/cm 3 . Here the data are presented as graphs chosen to aid analysis within and across simulant groups. To gain additional insight into the processes happening at the elevated temperatures in the dielectric measurements, heat capacity data was obtained using differential scanning calorimetry (DSC) on several of the simulant materials. Our data show that over the frequency range 397 MHz – 2985 MHz a material’s behavior does not greatly change, as compared to the scale of differences observed between lunar mare and highland simulants at high temperatures. For example at 1000 °C, the mare simulant JSC-1A absorbs 10 times more power than the highland simulant NUW-LHT-5M. We observe that as melting temperatures are reached both permittivity and dielectric loss rise non-linearly, helping to explain thermal runaway during microware heating. Our data show that even less than a few weight % of many non-lunar minerals, and the use of mixtures in simulants can affect the dielectric behavior at higher temperatures. A comparison of our results with published dielectric data for Apollo samples and with remote sensing of the Moon supports the conclusion the simulants and lunar material at room temperature have very similar dielectric values.

Differential Scanning Calorimetry

Cryogenic Propellant Tank Drain Simulation with a Stratified Ullage Nodal Model

Cryogenic propellant tank pressurization systems maintain tank pressure within specified limits to comply with propellant thermodynamic and tank structural design requirements. Many pump-fed cryogenic liquid propulsion systems control tank pressure using an autogenous pressurization system. Autogenous pressurization systems vaporize and heat propellant to produce the pressurant gas delivered to the tank. While autogenous pressurization systems eliminate the need for a separate pressurant gas storage and delivery system, they couple pressurant mass requirements to propellant mass requirements. Accurate prediction of autogenous pressurant mass during propellant drain is essential in optimizing cryogenic propellant tank designs.

Richard G Hibbs

Atomistic Simulation of Glasses and Amorphous Materials: Challenges and Opportunities for the Next Decade

Atomistic simulations have become indispensable tools for understanding glass structure, dynamics, and properties, yet persistent challenges limit their predictive power. This perspective examines three interconnected issues, namely glass formation procedures, interatomic potential development, and machine learning applications, which emerged from the 5th International Workshop on Challenges of Atomistic Simulations of Glasses and Amorphous Materials. We identify convergent community priorities for (i) standardized validation protocols, (ii) curated benchmark datasets with complete metadata, and (iii) open repositories for glasses. A systematic was forward is provided by a hierarchical validation framework for assessing the structural fidelity, property prediction, and behavioral realism of simulation techniques. Looking ahead, transformative advances are promised by the fusion of classical techniques with machine learning based approaches, for instance, by integrating swap Monte Carlo with machine-learning (ML) potentials, leveraging foundation models through transfer learning, and finetuning ML potentials with experimental data. Progress depends on the community committing to validated models, reproducible protocols, and sustained data sharing.

Krishnan, N. M. Anoop

Effects of Motion Fidelity on Pilot Performance Using Mixed-Reality Visuals

This paper investigates the effects of simulator motion fidelity on pilot performance using mixed-reality visuals. In NASA’s Vertical Motion Simulator, eleven pilots flew a lift-plus-cruise electric vertical takeoff and landing vehicle in a vertiport approach and landing task under four different motion conditions: no motion, small hexapod motion, large hexapod motion, and full motion. The simulator motion condition statistically significantly affected: 1) the severity of cybersickness, 2) motion ratings, 3) final approach and take-off boundary crossing groundspeeds, 4) longitudinal touchdown point, and 5) battery time remaining after landing. These results suggest that the use of motion systems in simulators with virtual- or mixed-reality visuals is warranted.

simulation

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING

Wettability of Two-Dimensional Carbon Allotropes from Molecular Simulations

Force-field Monte Carlo and Molecular Dynamics simulations are used to compare wetting behaviors of model carbon sheets mimicking neat graphene, its saturated derivative, graphane, and related planar allotropes penta-graphene, γ-graphyne, and ψ-graphene in contact with aqueous droplets or an aqueous film confined between parallel carbon sheets. Atomistic and area-integrated surface/water potentials are found to be essentially equivalent in capturing moderate differences between the wetting free energies of tested substrates. Despite notable differences in mechanical and electric properties of distinct allotropes, the predicted allotrope/water contact angles span a narrow window of weakly hydrophilic values. Contact angles in the range of 80 ± 10° indicate modest hydration repulsion incapable of competing with van der Waals attraction between carbon particles. Poor dispersibility in neat water is hence a common feature of studied materials.

Biochemistry & Molecular Biology

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle