Search NASA⌕ Search

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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Non-propagating structures and propagating waves in solar wind turbulence revealed by simulations and observations

Structures and waves are common features of solar wind turbulence at various scales. The interplay between structures and waves is important for processes such as the turbulent energy cascade, plasma heating, and particle scattering. Our understanding of turbulence has been advanced by not only new space missions and numerical simulations, but also techniques that have been developed to interpret the rapidly growing turbulence data. We review basic models of turbulence with a specific focus on the analysis methods for understanding magnetic structures and waves. MHD and kinetic waves in single-spacecraft time series measurements can be identified through mode decomposition or their characteristic polarization signatures. The structures in this paper are considered as zero-frequency, non-propagating or convected modes embedded in the solar wind. The synergy between observations and simulations is most evident in the application of spatial-temporal analysis to multi-spacecraft observation and turbulence simulations. The spatial-temporal analysis has greatly improved our understanding of structures and waves in turbulence. We conclude by discussing prospects for future research.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Data-centric framework for crystal structure identification in atomistic simulations using machine learning

Atomic-level modeling performed at large scales enables the investigation of mesoscale materials properties with atom-by-atom resolution. The spatial complexity of such cross-scale simulations renders them unsuitable for simple human visual inspection. Instead, specialized structure characterization techniques are required to aid interpretation. These have historically been challenging to construct, requiring significant intuition and effort. Here we propose an alternative framework for a fundamental structural characterization task: classifying atoms according to the crystal structure to which they belong. Our approach is data-centric and favors the employment of Machine Learning over heuristic rules of classification. A group of data-science tools and simple local descriptors of atomic structure are employed together with an efficient synthetic training set. We also introduce the first standard and publicly available benchmark data set for evaluation of algorithms for crystal-structure classification. Further, it is demonstrated that our data-centric framework outperforms all of the most popular heuristic methods—especially at high temperatures when lattices are the most distorted—while introducing a systematic route for generalization to new crystal structures. Moreover, through the use of outlier detection algorithms our approach is capable of discerning between amorphous atomic motifs (i.e., noncrystalline phases) and unknown crystal structures, making it uniquely suited for exploratory materials synthesis simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

In crystallo observation of three metal ion promoted DNA polymerase misincorporation

Error-free replication of DNA is essential for life. Despite the proofreading capability of several polymerases, intrinsic polymerase fidelity is in general much higher than what base-pairing energies can provide. Although researchers have investigated this long-standing question with kinetics, structural determination, and computational simulations, the structural factors that dictate polymerase fidelity are not fully resolved. Time-resolved crystallography has elucidated correct nucleotide incorporation and established a three-metal-ion-dependent catalytic mechanism for polymerases. Using X-ray time-resolved crystallography, we visualize the complete DNA misincorporation process catalyzed by DNA polymerase η. The resulting molecular snapshots suggest primer 3´-OH alignment mediated by A-site metal ion binding is the key step in substrate discrimination. Moreover, we observe that C-site metal ion binding preceded the nucleotidyl transfer reaction and demonstrate that the C-site metal ion is strictly required for misincorporation. Our results highlight the essential but separate roles of the three metal ions in DNA synthesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Biomolecular Solvation Structure Revealed by Molecular Dynamics Simulations

To compare ordered water positions from experiment with those from molecular dynamics (MD) simulations, a number of MD models of water structure in crystalline endoglucanase were calculated. The starting MD model was derived from a joint X-ray and neutron diffraction crystal structure, enabling the use of experimentally assigned protonation states. Simulations were performed in the crystalline state, using a periodic 2 × 2 × 2 supercell with explicit solvent. Water X-ray and neutron scattering density maps were computed from MD trajectories using standard macromolecular crystallography methods. In one set of simulations, harmonic restraints were applied to bias the protein structure toward the crystal structure. For these simulations, the recall of crystallographic waters using strong peaks in the MD water electron density was very good, and there also was substantial visual agreement between the boomerang-like wings of the neutron scattering density and the crystalline water hydrogen positions. An unrestrained simulation also was performed. For this simulation, the recall of crystallographic waters was much lower. For both restrained and unrestrained simulations, the strongest water density peaks were associated with crystallographic waters. Overall, the results demonstrate that it is now possible to recover crystallographic water structure using restrained MD simulations but that it is not yet reasonable to expect unrestrained MD simulations to do the same. Further development and generalization of MD water models for force-field development, macromolecular crystallography, and medicinal chemistry applications is now warranted. In particular, the combination of room-temperature crystallography, neutron diffraction, and crystalline MD simulations promises to substantially advance modeling of biomolecular solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Forest structural complexity and ignition pattern influence simulated prescribed fire effects

Background: Forest structural characteristics, the burning environment, and the choice of ignition pattern each influence prescribed fire behaviors and resulting fire effects; however, few studies examine the influences and interactions of these factors. Understanding how interactions among these drivers can influence prescribed fire behavior and effects is crucial for executing prescribed fires that can safely and effectively meet management objectives. To analyze the interactions between the fuels complex and ignition patterns, we used FIRETEC, a three-dimensional computational fluid dynamics fire behavior model, to simulate fire behavior and effects across a range of horizontal and vertical forest structural complexities. For each forest structure, we then simulated three different prescribed fires each with a unique ignition pattern: strip-head, dot, and alternating dot. Results: Forest structural complexity and ignition pattern affected the proportions of simulated crown scorch, consumption, and damage for prescribed fires in a dry, fire-prone ecosystem. Prescribed fires in forests with complex canopy structures resulted in increased crown consumption, scorch, and damage compared to less spatially complex forests. The choice of using a strip-head ignition pattern over either a dot or alternating-dot pattern increased the degree of crown foliage scorched and damaged, though did not affect the proportion of crown consumed. We found no evidence of an interaction between forest structural complexity and ignition pattern on canopy fuel consumption, scorch, or damage. Conclusions: We found that forest structure and ignition pattern, two powerful drivers of fire behavior that forest managers can readily account for or even manipulate, can be leveraged to influence fire behavior and the resultant fire effects of prescribed fire. These simulation findings have critical implications for how managers can plan and perform forest thinning and prescribed burn treatments to meet risk management or ecological objectives.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A fast particle-mesh simulation of non-linear cosmological structure formation with massive neutrinos

Quasi-N-body simulations, such as FastPM, provide a fast way to simulate cosmological structure formation, but have yet to adequately include the effects of massive neutrinos. In this work, we present a method to include neutrino particles in FastPM, enabling computation of the CDM and total matter power spectra to percent-level accuracy in the non-linear regime. The CDM-neutrino cross-power can also be computed at a sufficient accuracy to constrain cosmological observables. To avoid the shot noise that typically plagues neutrino particle simulations, we employ a quasi-random algorithm to sample the relevant Fermi-Dirac distribution when setting the initial neutrino thermal velocities. We additionally develop an effective distribution function to describe a set of non-degenerate neutrinos as a single particle to speed up non-degenerate simulations. The simulation is accurate for the full range of physical interest, M ν ≲ 0.6eV, and applicable to redshifts z ≲ 2. Such accuracy can be achieved by initializing particles with the two-fluid approximation transfer functions (using the REPS package). Convergence can be reached in ~ 25 steps, with a starting redshift of z=99. Probing progressively smaller scales only requires an increase in the number of CDM particles being simulated, while the number of neutrino particles can remain fixed at a value less than or similar to the number of CDM particles. In turn, the percentage increase in runtime-per-step due to neutrino particles is between ~ 5-20% for runs with 1024 3 CDM particles, and decreases as the number of CDM particles is increased. The code has been made publicly available, providing an invaluable resource to produce fast predictions for cosmological surveys and studying reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

Chemistry↗
DOE OSTI2019

Crystallographic structure and molecular dynamics simulations of the major endoglucanase from Xanthomonas campestris pv. campestris shed light on its oligosaccharide products release pattern

Cellulases are critical enzymatic components for the transformation of plant biomass into fuels, renewable materials and green chemicals. In this work, we determined the crystal structure, pattern of hydrolysis products release, and conducted molecular dynamics simulations of the major endoglucanase from the Xanthomonas campestris pv. campestris (XccCel5A). XccCel5A has a TIM barrel fold with the catalytic site centrally placed in a binding groove surrounded by aromatic side chains. Molecular dynamics simulations show that productive position of the substrate is secured by a network of hydrogen bonds in the four main subsites, which differ in details from homologous structures. Capillary zone electrophoresis and computational studies reveal XccCel5A can act both as endoglucanase and licheninase, but there are preferable arrangements of substrate regarding β-1,3 and β-1,4 bonds within the binding cleft which are related to the enzymatic efficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improving scalability of electronic structure code for molecular simulations in the presence of environment

A scalable density functional electronic code with Gaussian basis set, called UTEP-NRLMOL, is developed to perform simulations of molecular systems in the presence of the environment with particular attention to the memory requirements. In the electronic structure calculations, the memory and computation time are proportional to the number of atoms. Memory requirements for density functional calculations scale as N*N, where N is the number of atoms. While the recent advances in HPC offer platforms with large numbers of cores, the limited amount of memory available on a given node and poor scalability of the electronic structure codes hinder their efficient usage of these platforms. We have introduced new scaling and parallelization paradigms using MPI-3 shared-memory functionality combined with usage of sparse algebra and storage of matrices in sparse format. This extends the range of applicability of the UTEP-NRLMOL code to large systems over 10,000 atoms, or using up to 67,000 basis functions, and making use of HPC architectures using over 6,000 processors utilizing all available cores. We have also interfaced code with effective fragment potential and polarizable continuum model libraries. The code was used in simulations of several applications which are published in reputed scientific journals.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Demonstrate multi-turbine simulation with hybrid-structured / unstructured-moving-grid software stack running primarily on GPUs and propose improvements for successful KPP-2

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, capturing the thin boundary layers, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources.

17 WIND ENERGY↗
DOE OSTI2025

Monthly averages of ED2 model simulations initialized with airborne lidar structure, Jan 1981-Dec 2018, Brazilian Amazon

Deforestation and forest degradation (selective logging, fires, fragmentation) have impacted nearly 40% of the original extent of the Brazilian Amazon, and have markedly impacted forest structure across the region. To date, few studies analysed how shifts in forest structure from degradation influence the forest sensitivity to climate extremes, because of the complex interactions between forest structure and micro-environmental conditions. To address this knowledge gap, we carried out a series of simulations across the Brazilian Amazon using the Ecosystem Demography Model (ED2), using observed forest structure derived from 541 airborne lidar transects (375 ha each) and two scenarios representing forest recovery and expansion of degradation to investigate how shifts in forest structure impact ecosystem function under near-average and extreme climate conditions, as part of the manuscript Longo et al 2025 "Degradation and Deforestation Increase the Sensitivity of the Amazon Forest to Climate Extremes". This dataset provides the output results from the ED2 model simulations for the three simulations at monthly time scales, in NetCDF format. For all simulations, we used bias-corrected hourly reanalyses (WFDE5) for most meteorological drivers, except for precipitation, which was obtained from CHIRPS. The meteorological drivers used in the study span 38 years (Jan 1981–Dec 2018). The output results correspond to the last 38 years of simulation (one full cycle of meteorological drivers), in which ED2 simulations used static stand structure (i.e., the forest structure was held constant). The following files are provided:ED2_emean_Global_R004_BrAmaz_s1c0t0l0f0.nc. This corresponds to the Control simulation. The forest structure was obtained from the airborne lidar.ED2_emean_Global_R005_BrAmaz_s1c0t1l1f0.nc. This corresponds to the Degraded simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that expanded deforestation and selective logging across the Amazon.ED2_emean_Global_R006_BrAmaz_s1c0t1l0f0.nc. This corresponds to the Recovery simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that completely halted deforestation and degradation, allowing degraded forests to recover for 38 years.We also provide file ED2_zones_R004_BrAmaz_s1c0t0l0f0.nc, which classifies each grid cell into zones used in the reference manuscript: 1: Southeast. 2: South. 3: West. 4: Central. 5: Northeast. 6: North. 7: Northwest". Index 0 corresponds to grid cells excluded from sub-region analyses because they were dominated by flooded forests, deforestation, and naturally non-forest vegetation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

A computational risk assessment approach to the integration of seismic and flooding hazards with internal hazards

Probabilistic risk assessment (PRA) of nuclear power plants historically focuses on internal events at the plant, rather than external hazards. Although the importance of external hazards risk analysis is now well recognized, the methods for assessing the risk associated with low probability external hazards rely heavily on subjective judgment of specialists, often using conservative elements in the analysis, as is the case with the Conservative Deterministic Failure Margin approach commonly used for seismic PRA. The U.S. Department of Energy Light Water Reactor Sustainability program has been investigating the use of first principles analyses to provide a more realistic assessment of the risk of external hazards, in a process called computational risk assessment (CRA). Efforts undertaken are described in the development of CRA methods to integrate the risk associated with seismic and flooding events into traditional PRA. Further, the results of four case studies are presented.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Simulation of the Electronic Structure of Graphene–Polyvinylidene Fluoride Composite Material

Theoretical models of shockproof composite materials based on two-layer graphenes and multilayer polyvinylidene fluoride C{sub 124}H{sub 40} + n(H–(C{sub 2}H{sub 2}F{sub 2}){sub 5}–H) (n = (1, 8)-bar ) are constructed. The electronic structure is studied using the semiempirical PM3 method that is one of version of the molecular orbital method. The orbital energies, ionization potentials, total electron energies, strength, and other properties of the considered material are calculated based on the theoretical models. The outlooks for application of these materials in the military field for manufacturing superstrong and lightweight flak jackets are considered.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

FY2022 Q4: Demonstrate multi-turbine simulation with hybrid-structured/unstructured-moving-grid software stack running primarily on GPUs and propose improvements for successful KPP-2 [Poster]

Milestone accomplishments staged the ExaWind team for successful completion of KPP-2 challenge problem in FY23, which requires the simulation on Frontier of at least four MW-scale turbines in an atmospheric boundary layer with at least 20B gridpoints. The ExaWind project and software stack is many faceted, with team members working on multiple areas, including linear-system solvers (Trilinos, hypre, AMReX), overset meshes, turbulence modeling, and in situ visualization, all with an aim for high fidelity predictions and performance portability. This milestone marks significant improvements on many fronts and provides the team with a pathway to exascale wind farm simulations in FY23.

17 WIND ENERGY↗
DOE OSTI2024

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Finally, other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗