SEARCH · Search NASA

Results for “Structural simulation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗

NASA NTRS1985

Coherent structures in a turbulent mixing layer - A comparison between direct numerical simulations and experiments

An eduction scheme has been developed in an attempt to determine the characteristics of large-scale vortical structures in a turbulent mixing layer. This analysis scheme has been applied to a set of experimental data taken in a new, larger mixing layer facility designed to minimize boundary and resonance effects. A similar scheme has been developed to apply to the results of a direct numerical simulation of a temporally growing mixing layer. A comparison of the two approaches shows important similarities in the coherent structures. The numerical simulations indicate that low levels of coherent forcing can dramatically change the evolution of the mixing layer. In the absence of such forcing, the numerical simulations and experiments show a lack of regularity in the transverse position, spacing, amplitude, shape and spanwise coherence of the large-scale vortical structures.

Metcalfe, R. W.↗

DOE OSTI2025

Quantitative Structure Determination from Experimental Four-Dimensional Scanning Transmission Electron Microscopy via the Scattering Matrix

Considerable inroads have recently been made on algorithms to determine the sample potential from four-dimensional scanning transmission electron microscopy data from thick samples where multiple scattering cannot be neglected. This paper further develops the scattering matrix approach to such structure determination. Through simulation, we demonstrate how this approach can be modified to better handle partial spatial coherence, unknown probe defocus, and information from the dark field region. By combining these developments we reconstruct the electrostatic potential of a monolithic SrTiO 3 crystal showing good quantitative agreement with the expected structure.

4D STEM↗

DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Simulated x-ray absorption spectra calculated with the FEFF code for various models of Cu incorporated into the iron oxide mineral hematite. The structures are drawn from publically available ab initio molecular dynamics simulations.

FEFF simulations of Cu doped hematite.

Prange, Micah P↗

NASA NTRS2002

26th JANNAF Airbreathing Propulsion Subcommittee Meeting

This volume, the first of four volumes, is a collection of 28 unclassified/unlimited-distribution papers which were presented at the Joint Army-Navy-NASA-Air Force (JANNAF) 26th Airbreathing Propulsion Subcommittee (APS) was held jointly with the 38th Combustion Subcommittee (CS), 20th Propulsion Systems Hazards Subcommittee (PSHS), and 2nd Modeling and Simulation Subcommittee. The meeting was held 8-12 April 2002 at the Bayside Inn at The Sandestin Golf & Beach Resort and Eglin Air Force Base, Destin, Florida. Topics covered include: scramjet and ramjet R&D program overviews; tactical propulsion; space access; NASA GTX status; PDE technology; actively cooled engine structures; modeling and simulation of complex hydrocarbon fuels and unsteady processes; and component modeling and simulation.

Fry, Ronald S.↗

NASA NTRS1987

Composite mechanics for engine structures

Recent research activities and accomplishments at Lewis Research Center on composite mechanics for engine structures are summarized. The activities focused mainly on developing procedures for the computational simulation of composite intrinsic and structural behavior. The computational simulation encompasses all aspects of composite mechanics, advanced three-dimensional finite-element methods, damage tolerance, composite structural and dynamic response, and structural tailoring and optimization.

Chamis, Christos C.↗

NASA NTRS1989

Composite mechanics for engine structures

Chamis, Christos C.↗

DOE OSTI2018

Highly parallelisable simulations of time-dependent viscoplastic fluid flow with structured adaptive mesh refinement

We present the extension of an efficient and highly parallelisable framework for incompressible fluid flow simulations to viscoplastic fluids. The system is governed by incompressible conservation of mass, the Cauchy momentum equation, and a generalised Newtonian constitutive law. In order to simulate a wide range of viscoplastic fluids, we employ the Herschel-Bulkley model for yield-stress fluids with nonlinear stress-strain dependency above the yield limit. We utilise Papanastasiou regularisation in our algorithm to deal with the singularity in apparent viscosity. The resulting system of partial differential equations is solved using the IAMR (Incompressible Adaptive Mesh Refinement) code, which uses second-order Godunov methodology for the advective terms and semi-implicit diffusion in the context of an approximate projection method to solve adaptively refined meshes. By augmenting the IAMR code with the ability to simulate regularised Herschel-Bulkley fluids, we obtain efficient numerical software for time-dependent viscoplastic flow in three dimensions, which can be used to investigate systems not considered previously due to computational expense. We validate results from simulations using this new capability against previously published data for Bingham plastics and power-law fluids in the two-dimensional lid-driven cavity. In doing so, we expand the range of Bingham and Reynolds numbers which have been considered in the benchmark tests. Moreover, extensions to time-dependent flow of Herschel-Bulkley fluids and three spatial dimensions offer new insights into the flow of viscoplastic fluids in this test case, and we provide missing benchmark results for these extensions.

42 ENGINEERING↗

DOE OSTI2018

Study of Lasing in Liquid-Crystal Systems with Microgratings

Waveguide lasing in a layer of nematic liquid crystal (NLC) doped with a dye has been experimentally observed. A chromium micrograting with a period of 1.76 μm has been deposited on the surface one of the cell substrates to provide a distributed feedback (DFB) and partially extract laser radiation to a quartz substrate. The DFB provides lasing in the ninth diffraction order for waveguide TM modes. Laser radiation is observed at the output of the substrate end face at an angle of 67.0° ± 1.5° with respect to the normal to the waveguide plane. An increase in voltage across the micrograting electrodes leads to red-shifted multimode lasing. A numerical simulation of a structure imitating the experimental cell has shown good agreement between the calculation results and experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Ship-based Observations and Climate Model Simulations of Cloud Phase over the Southern Ocean

The Southern Ocean (SO) clouds exert a significant influence on the Earth’s radiation budget. Here we analyzed ship-based remote sensing observations of SO clouds over a five-month long DOE ARM Measurements of Aerosols, Radiation and Clouds over the Southern Oceans (MARCUS) field campaign to better understand cloud phase variability. We developed a method to classify eight categories of hydrometeors (ice, liquid, mixed phase, rain, drizzle, snow, aerosols, and clear sky) based on measurements of lidar, radar and radiosondes. Cloud thermodynamic phases (liquid, ice and mixed phase) at coarser scales were further derived to compare with the DOE Energy Exascale Earth System Model version 1 (E3SMv1) simulation and the Earth Model Column Collaboratory (EMC 2 ) instrument simulator. For a scale-aware comparison with climate model simulations, we found that spatially averaging the raw remote sensing data (e.g., backscatter, reflectivity) results in increased cloud cover and cloud liquid, whereas, directly averaging cloud phase from higher to lower resolution maintains clear air regions and is thus recommended for in-cloud frequency comparisons. For cloud thermodynamic phases in stratiform clouds, the E3SM underestimates cloud ice and overestimates cloud liquid at temperatures between -40 and 0°C. When latitudes increase, both observations and simulations show a transition of dominant phase from liquid to ice for cloud tops as well as for the entire cloud columns, but the model underestimation of ice phase is more severe at higher latitudes. Such model bias is unlikely caused by spatial scale differences or lack of heterogeneity in cloud vertical structure in the simulation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Differentiable Predictions for Large Scale Structure with SHAMNet

In simulation-based models of the galaxy-halo connection, theoretical predictions for galaxy clustering and lensing are typically made based on Monte Carlo realizations of a mock universe. In this paper, we use Subhalo Abundance Matching (SHAM) as a toy model to introduce an alternative to stochastic predictions based on mock population, demonstrating how to make simulation-based predictions for clustering and lensing that are both exact and differentiable with respect to the parameters of the model. Conventional implementations of SHAM are based on iterative algorithms such as Richardson-Lucy deconvolution; here we use the JAX library for automatic differentiation to train SHAMNet, a neural network that accurately approximates the stellar-to-halo mass relation (SMHM) defined by abundance matching. In our approach to making differentiable predictions of large-scale structure, we map parameterized PDFs onto each simulated halo, and calculate gradients of summary statistics of the galaxy distribution by using autodiff to propagate the gradients of the SMHM through the statistical estimators used to measure one- and two-point functions. Our techniques are quite general, and we conclude with an overview of how they can be applied in tandem with more complex, higher-dimensional models, creating the capability to make differentiable predictions for the multi-wavelength universe of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Cabana: A Performance Portable Library for Particle-Based Simulations

Particle-based simulations are ubiquitous throughout many fields of computational science and engineering, spanning the atomistic level with molecular dynamics (MD), to mesoscale particle-in-cell (PIC) simulations for solid mechanics, device-scale modeling with PIC methods for plasma physics, and massive N-body cosmology simulations of galaxy structures, with many other methods in between (Hockney & Eastwood, 1989). While these methods use particles to represent significantly different entities with completely different physical models, many low-level details are shared including performant algorithms for short- and/or long-range particle interactions, multi-node particle communication patterns, and other data management tasks such as particle sorting and neighbor list construction. Cabana is a performance portable library for particle-based simulations, developed as part of the Co-Design Center for Particle Applications (CoPA) within the Exascale Computing Project (ECP) (Alexander et al., 2020). The CoPA project and its full development scope, including ECP partner applications, algorithm development, and similar software libraries for quantum MD, is described in (Mniszewski et al., 2021). Cabana uses the Kokkos library for on-node parallelism (Edwards et al., 2014; Trott et al., 2022), enabling simulation on multi-core CPU and GPU architectures, and MPI for GPU-aware, multi-node communication. Cabana provides particle simulation capabilities on almost all current Kokkos backends, including serial execution, OpenMP (including OpenMP-Target for GPUs), CUDA (NVIDIA GPUs), HIP (AMD GPUs), and SYCL (Intel GPUs), providing a clear path for the coming generation of accelerator-based exascale hardware. Cabana builds on Kokkos by providing new particle data structures and particle algorithms resulting in a similar execution policy-based, node-level programming model that is intended to be used in addition to the core Kokkos library within an application. Cabana is designed as an application and physics agnostic, but particle-specific toolkit which can either be used to generate a new application, or to be used as needed in existing applications at various levels of invasiveness including through interfaces that wrap user memory in existing data structures.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Metastable silica high pressure polymorphs as structural proxies of deep Earth silicate melts

Modelling of processes involving deep Earth liquids requires information on their structures and compression mechanisms. However, knowledge of the local structures of silicates and silica (SiO 2 ) melts at deep mantle conditions and of their densification mechanisms is still limited. Here we report the synthesis and characterization of metastable high-pressure silica phases, coesite-IV and coesite-V, using in situ single-crystal X-ray diffraction and ab initio simulations. Their crystal structures are drastically different from any previously considered models, but explain well features of pair-distribution functions of highly densified silica glass and molten basalt at high pressure. Built of four, five-, and six-coordinated silicon, coesite-IV and coesite-V contain SiO 6 octahedra, which, at odds with 3 rd Pauling’s rule, are connected through common faces. Our results suggest that possible silicate liquids in Earth’s lower mantle may have complex structures making them more compressible than previously supposed.

36 MATERIALS SCIENCE↗

NASA NTRS2001

Properties of Air Traffic Conflicts for Free and Structured Routing

This paper analyzes the properties of air traffic conflicts in a future free routing system against those in the current structured routing system. Simulation of en route air traffic operations (above 18,000 ft) over the contiguous United States for a 24-hour period, constructed with initial conditions from actual air traffic data, were conducted using the Future ATM Concepts Evaluation Tool (FACET). Free routes were modeled as great circle (direct) routes from origin to destination, and structured routes were derived from actual flight plans along the current system of air routes. The conflict properties analyzed in this study include: (1) Total number of conflicts; (2) Distributions of key conflict parameters; and, (3) Categorization of conflicts into independent conflicts and two types of interacting conflicts. Preliminary results (for Denver Center traffic) indicate that conflict properties in a free routing system are different from those in the current structured routing system. In particular, a free routing system has significantly fewer conflicts, involving a correspondingly smaller number of aircraft, compared to the current structured routing system. Additionally, the conflict parameter distributions indicate that free routing conflicts are less intrusive than structured routing conflicts, and would therefore require small trajectory deviations for resolution.

Bilimoria, Karl D.↗

NASA NTRS1972

A simulation of saltation dynamics with application to terrestrial and Martian model surfaces.

The feasibility of predicting the saltation threshold from a dynamic model of the saltation process is demonstrated. The model is based on mechanical equations that describe the trajectories of particles (sand grains) near the ground-atmosphere interface. The equations of particle motion are formulated in terms of forces and moments. The forces employed are based on theoretical and experimental results obtained for continuum flow. Surface structuring analytics and simulated fluid turbulence further modify the particle behavior. The model is used to predict the threshold velocities on Mars.

Andres, R. M.↗

DOE OSTI2020

Validation and Sensitivity Analysis of a 1-D Lake Model Across Global Lakes

Lakes have important influence on weather and climate from local to global scales. However, their prediction using numerical models is notoriously difficult because lakes are highly heterogeneous across the globe, but observations are sparse. In this study, we assessed the performance of a 1-D lake model in simulating the thermal structures of 58 lakes with diverse morphometric and geographic characteristics by following the phase 2a local lake protocol of the Inter-sectoral Impact Model Intercomparison Project (ISIMIP2a). After calibration, the root-mean-square errors (RMSE) were below 2 °C for 70% and 75% of the lakes for epilimnion and full-profile temperature simulations, with an average of 1.71 °C and 1.43 °C, respectively. The model performance mainly depended on lake shape rather than location, supporting the possibility of grouping model parameters by lake shape for global applications. Furthermore, through machine-learning based parameter sensitivity tests, we identified turbulent heat fluxes, wind-driven mixing and water transparency as the major processes controlling lake thermal and mixing regimes. Snow density was also important for modeling the ice phenology of high-latitude lakes. The relative influence of the key processes and the corresponding parameters mainly depended on lake latitude and depth. Turbulent heat fluxes showed a decreasing importance in affecting epilimnion temperature with increasing latitude. Wind-driven mixing was less influential to lake stratification for deeper lakes while the impact of light extinction, on the contrary, showed a positive correlation with depth. Our findings may guide improvements in 1-D lake model parameterizations to achieve higher fidelity in simulating global lake thermal dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Investigating the Impact of Thickness, Calendering and Channel Structures of Printed Electrodes on the Energy Density of LIBs - 3D Simulation and Validation

Current lithium ion batteries (LIBs) are expensive and bulky, limited by relatively low charging rates. To increase the rate of charging and reduce weight, thin electrodes with high energy density are required. The increase in energy density can be achieved by several techniques including boosting electrolyte transport, high loading/utilization of active material, employing high conductive electrolytes and electrodes with advanced architectures, and increasing cell temperature. In this paper, a 3D physics-based electrochemical model of LIBs is developed in COMSOL simulation software for different thickness, calendering steps as well as channel structures (conical, cylindrical) to optimize the electrode design and in turn maximize volumetric energy density. The simulation results demonstrated that calendering the electrodes with high initial porosity increases the volumetric energy density of the cell. In addition, cylindrical channel structures with relatively lower edge-to-edge distance also results in increased volumetric energy density. The simulation results of the 3D model was validated by comparing it with experimental results.

improving volumetric energy density↗