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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Simulations of cm-wavelength Sunyaev-Zel'dovich galaxy cluster and point source blind sky surveys and predictions for the RT32/OCRA-f and the Hevelius 100-m radio telescope

We investigate the effectiveness of blind surveys for radio sources and galaxy cluster thermal Sunyaev-Zel'dovich effects (TSZEs) using the four-pair, beam-switched OCRA-f radiometer on the 32-m radio telescope in Poland. The predictions are based on mock maps that include the cosmic microwave background, TSZEs from hydrodynamical simulations of large scale structure formation, and unresolved radio sources. We validate the mock maps against observational data, and examine the limitations imposed by simplified physics. We estimate the effects of source clustering towards galaxy clusters from NVSS source counts around Planck-selected cluster candidates, and include appropriate correlations in our mock maps. The study allows us to quantify the effects of halo line-of-sight alignments, source confusion, and telescope angular resolution on the detections of TSZEs. We perform a similar analysis for the planned 100-m Hevelius radio telescope (RTH) equipped with a 49-beam radio camera and operating at frequencies up to 22 GHz.We find that RT32/OCRA-f will be suitable for small-field blind radio source surveys, and will detect 33{sup +17}{sub −11} new radio sources brighter than 0.87 mJy at 30 GHz in a 1 deg{sup 2} field at > 5σ CL during a one-year, non-continuous, observing campaign, taking account of Polish weather conditions. It is unlikely that any galaxy cluster will be detected at 3σ CL in such a survey. A 60-deg{sup 2} survey, with field coverage of 2{sup 2} beams per pixel, at 15 GHz with the RTH, would find <1.5 galaxy clusters per year brighter than 60 μJy (at 3σ CL), and would detect about 3.4 × 10{sup 4} point sources brighter than 1 mJy at 5σ CL, with confusion causing flux density errors ∼< 2% (20%) in 68% (95%) of the detected sources.A primary goal of the planned RTH will be a wide-area (π sr) radio source survey at 15 GHz. This survey will detect nearly 3 × 10{sup 5} radio sources at 5σ CL down to 1.3 mJy, and tens of galaxy clusters, in one year of operation with typical weather conditions. Confusion will affect the measured flux densities by ∼< 1.5% (16%) for 68% (95%) of the point sources. We also gauge the impact of the RTH by investigating its performance if equipped with the existing RT32 receivers, and the performance of the RT32 equipped with the RTH radio camera.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Simulation of the Formation of a Cascade of Displacements and Transient Ionization Processes in Silicon Semiconductor Structures under Neutron Exposure

The formation of a disordered defect region in bulk silicon is simulated using the molecular-dynamics method for various energies of a primary recoil atom. Variations in the volume and number of radiation-induced defects in a cluster during its formation are calculated. The generation rates of nonequilibrium carriers and amplitude-temporal dependences of pulses of ionization currents in test Schottky diodes with hyperhigh frequencies are found theoretically.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Learning to predict the cosmological structure formation

Matter evolved under the influence of gravity from minuscule density fluctuations. Nonperturbative structure formed hierarchically over all scales and developed non-Gaussian features in the Universe, known as the cosmic web. To fully understand the structure formation of the Universe is one of the holy grails of modern astrophysics. Astrophysicists survey large volumes of the Universe and use a large ensemble of computer simulations to compare with the observed data to extract the full information of our own Universe. However, to evolve billions of particles over billions of years, even with the simplest physics, is a daunting task. We build a deep neural network, the Deep Density Displacement Model ( D 3 M ), which learns from a set of prerun numerical simulations, to predict the nonlinear large-scale structure of the Universe with the Zel’dovich Approximation (ZA), an analytical approximation based on perturbation theory, as the input. Our extensive analysis demonstrates that D 3 M outperforms the second-order perturbation theory (2LPT), the commonly used fast-approximate simulation method, in predicting cosmic structure in the nonlinear regime. We also show that D 3 M is able to accurately extrapolate far beyond its training data and predict structure formation for significantly different cosmological parameters. Our study proves that deep learning is a practical and accurate alternative to approximate 3D simulations of the gravitational structure formation of the Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Nanoscale structure of the oil-water interface

X-ray reflectivity (XR) and atomistic molecular dynamics (MD) simulations, carried out to determine the structure of the oil-water interface, provide new insight into the simplest liquid-liquid interface. For several oils (hexane, dodecane, and hexadecane) the XR shows very good agreement with a monotonic interface-normal electron density profile (EDP) broadened only by capillary waves. Similar agreement is also found for an EDP including a sub-Å thick electron depletion layer separating the oil and the water. As a result, the XR and MD derived depletions are much smaller than reported for the interface between solid-supported hydrophobic monolayers and water.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Explicit Filtering and Reconstruction to Reduce Grid Dependence in Convective Boundary Layer Simulations Using WRF-LES

As model grid resolutions move from the mesoscale to the microscale, turbulent structures represented in atmospheric boundary layer simulations change dramatically. At intermediate resolutions, the so-called gray zone, turbulent motions are not resolved accurately, posing a challenge to numerical simulations. The representation of turbulence is also highly sensitive to the choice of closure model. Here, we examine explicit filtering and reconstruction in the gray zone as a technique to better represent atmospheric turbulence. The convective boundary layer is simulated using the Weather Research and Forecasting (WRF) Model with horizontal resolutions ranging from 25 m to 1 km. Four large-eddy simulation (LES) turbulence models are considered: the Smagorinsky model, the TKE-1.5 model, and two versions of the dynamic reconstruction model (DRM). The models are evaluated by their ability to produce consistent mean potential temperature profiles, heat and momentum fluxes, velocity fields, and turbulent kinetic energy spectra as the grids become coarser. The DRM, a mixed model that uses an explicit filtering and reconstruction technique to account for resolvable subfilter-scale (RSFS) stresses, performs very well at resolutions of 500 m and 1 km without any special tuning, whereas the Smagorinsky and TKE-1.5 models produce heavily grid-dependent results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Local structure orders and glass forming ability of Ni-Nb liquids

In this paper, ab initio molecular dynamics simulations were performed to study the atomistic structures of Ni 48 Nb 52 , Ni 59.5 Nb 40.5 and Ni 75 Nb 25 liquids near their melting temperatures. It is found that the total coordination number only changes slightly across the compositions. However, there is significant difference in the chemical ordering. The Cluster Alignment method was used to identify the dominant local structure motifs. It is found that distorted and perfect icosahedral clusters are the most popular motifs around Ni atoms, while Z14, Z15, and Z16 Frank-Kasper polyhedra are dominant around Nb atoms. Among the three compositions, the eutectic composition (Ni 59.5 Nb 40.5 ) has the largest population of these dominant motifs and the highest degree of five-fold local symmetry. Chemical order in the distorted icosahedral motif is studied and the occupation probability of Nb in the distorted region is found to be much different from other sites. The network formed by distorted and perfect icosahedral clusters and that by the interconnection of the perfect icosahedra, Z14, Z15, and Z16 motifs are the strongest at the eutectic composition. Finally, the high percentage of perfect and distorted icosahedral clusters and their strong network structures in this system (especially around the eutectic composition) correlate well with the excellent glass forming ability in this system.

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

Enhancement of nuclear reaction rates in asymmetric binary ionic mixtures

Here, using orbital-free molecular dynamics simulations we study the structure and dynamics of increasingly asymmetric mixtures such as hydrogen–carbon, hydrogen–aluminium, hydrogen–copper, and hydrogen–silver. We show that, whereas the heavy component structure is close to an effective one-component plasma (OCP), the light component appears more structured than the corresponding OCP. This effect is related to the crossover towards a Lorentz-type diffusion triggered by strongly coupled, highly charged heavy ions, and witnessed by the change of temperature scaling laws of diffusion. Finally, this over-correlation translates into an enhancement of nuclear reaction rates much higher than its classical OCP counterpart.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

The Nature of the Technetium Species Formed During the Oxidation of Technetium Dioxide with Oxygen and Water

The oxidation of TcO 2 with O 2 in the presence of H 2 O at 250 °C was investigated; Tc 2 O 7 and a red product were observed after the reaction. UV/Visible spectroscopy and a cerium titration on the “red” product are consistent with the presence of a reduced Tc species. It is proposed that the “red” product was formed by the hydrolysis of Tc 2 O 7 to HTcO 4 followed by the reduction of HTcO 4 . The reaction products of TcO 2 with O 2 /H 2 O at 250 °C were also studied by mass spectrometry. The main dinuclear species were Tc 2 O 7 and Tc 2 O 5 , while the main mononuclear species were TcO 3 and HTcO 4 . The compound HTcO 4 and associated hydrates were simulated to identify candidate structures for the “red” product; however, neither the acid nor any related hydrates exhibited the UV/Visible signature of the “red” product.

Childs, Bradley C.↗
DOE OSTI2025

Discovery of a Ferromagnetic Nickel Chalcogenide Nanocluster Ni 3 S 3 H(PEt 3 ) 5

Atomically precise ligated nanoclusters (NC) are promising cluster-based materials with novel molecular architectures and tunable magnetic properties. Herein, the synthesis and characterization of a nickel sulfide NC Ni 3 S 3 H(PEt 3 ) 5 (PEt 3 = triethylphosphine) with distinct magnetic properties are reported. Magnetization measurements reveal its magnetic moment of 1.5 µ B in the solid phase, consistent with the existence of one unpaired electron predicted by density functional theory (DFT) calculations. Additionally, experimental measurements indicate the presence of ferromagnetic ordering within each Ni 3 S 3 H(PEt 3 ) 5 NC and strong coercivity at temperatures below 20 K. Ion mobility-mass spectrometry is employed in conjunction with DFT calculations and collision cross-section simulations to investigate the structure of the isolated Ni 3 S 3 H(PEt 3 ) 5 . Theoretical studies show that [Ni 3 S 3 H(PEt 3 ) 5 ] + has a planar Ni 3 S 3 core where three Ni atoms are arranged in a triangle with three bridging S atoms residing in the same plane. This structure is preserved in both solution and solid phases, which is confirmed by spectroscopic studies of Ni 3 S 3 H(PEt 3 ) 5 . Additionally, DFT calculations indicate that all spins at the Ni sites are aligned parallel, confirming the presence of ferromagnetic coupling. Overall, this study provides key insights into the structure and magnetic properties of Ni 3 S 3 H(PEt 3 ) 5 , which will facilitate the design of new NC-based magnetic materials.

Nickel sulfide nanocluster↗
DOE OSTI2017

Misfit dislocation patterns of Mg-Nb interfaces

The role of heterogeneous interfaces in improving mechanical properties of polycrystalline aggregates and laminated composites has been well recognized with interface structure being of fundamental importance in designing composites containing multiple interfaces. In this paper, taking the Mg (hexagonal close-packed (hcp))/Nb (body-centered cubic (bcc)) interface as an example, we develop Mg-Nb interatomic potentials for predicting atomic configurations of Mg/Nb interfaces. We systematically characterize interface dislocations of Mg/Nb interfaces with Nishiyama-Wassermann (NW) and Kurdjumov-Sachs (KS) orientation relationships and propose a generalized procedure of characterizing interface structure by combining atomistic simulation and interface dislocation theory, which is applicable for not only hcp/bcc interfaces, but also other systems with complicated interface dislocation configurations.Here, in Mg/Nb, interface dislocation networks of two types of interfaces are significantly different although they originate from partial dislocations of similar character: the NW interface is composed of three sets of partial dislocations, while the KS interface is composed of four sets of interface dislocations - three sets of partial dislocations and one set of full dislocations that forms from the reaction of two close partial dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A computational study of yttria-stabilized zirconia: I. Using crystal chemistry to search for the ground state on a glassy energy landscape

Yttria-stabilized zirconia (YSZ), a ZrO 2 -Y 2 O 3 solid solution that contains a large population of oxygen vacancies, is widely used in energy and industrial applications. Past computational studies correctly predicted the anion diffusivity but not the cation diffusivity, which is important for material processing and stability. One of the challenges lies in identifying a plausible configuration akin to the ground state in a glassy landscape. This is unlikely to come from random sampling of even a very large sample space, but the odds are much improved by incorporating packing preferences revealed by a modest sized configurational library established from empirical potential calculations. Ab initio calculations corroborated these preferences, which prove remarkably robust extending to the fifth cation-oxygen shell about 8 Å away. Yet because of frustration there are still rampant violations of packing preferences and charge neutrality in the ground state, and the approach toward it bears a close analogy to glass relaxations. Fast relaxations proceed by fast oxygen movement around cations, while slow relaxations require slow cation diffusion. The latter is necessarily cooperative because of strong coupling imposed by the long-range packing preferences.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling Concept and Numerical Simulation of Ultrasonic Wave Propagation in a Moving Fluid-Structure Domain based on a Monolithic Approach

Here in the present study, we propose a novel multiphysics model that merges two time-dependent problems – the Fluid-Structure Interaction (FSI) and the ultrasonic wave propagation in a fluid-structure domain with a one directional coupling from the FSI problem to the ultrasonic wave propagation problem. This model is referred to as the “eXtended fluid-structure interaction (eXFSI)” problem. This model comprises isothermal, incompressible Navier-Stokes equations with nonlinear elastodynamics using the Saint-Venant Kirchhoff solid model. The ultrasonic wave propagation problem comprises monolithically coupled acoustic and elastic wave equations. To ensure that the fluid and structure domains are conforming, we use the ALE technique. The solution principle for the coupled problem is to first solve the FSI problem and then to solve the wave propagation problem. Accordingly, the boundary conditions for the wave propagation problem are automatically adopted from the FSI problem at each time step. The overall problem is highly nonlinear, which is tackled via a Newton-like method. The model is verified using several alternative domain configurations. To ensure the credibility of the modeling approach, the numerical solution is contrasted against experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

iTOUGH2: A multiphysics simulation-optimization framework for analyzing subsurface systems

We report iTOUGH2 is a simulation-optimization framework for the TOUGH suite of nonisothermal multiphase flow models and related simulators of geophysical, geochemical, and geomechanical processes. After appropriate parameterization of subsurface structures and their properties, iTOUGH2 runs simulations for multiple parameter sets and analyzes the resulting output for parameter estimation through automatic model calibration, local and global sensitivity analyses, data-worth analyses, and uncertainty propagation analyses. Development of iTOUGH2 is driven by scientific challenges and user needs, with new capabilities continually added to both the forward simulator and the optimization framework. This paper provides a summary description of methods and features implemented in iTOUGH2, and discusses the usefulness and limitations of an integrated simulation-optimization workflow in support of the characterization and analysis of complex multiphysics subsurface systems.

58 GEOSCIENCES↗
DOE OSTI2022

Novel processing route for design and manufacturing of metal toughened nanoceramics: Al–Al 2 O 3 nanocermets

Monolithic alumina is a highly brittle ceramic with a very low fracture toughness, making it an unreliable material for mechanical and structural applications. Ceramic-metal composites (cermets) are the subject of intensive interest to obtain materials which can address the intrinsic brittleness of ceramics without compromising their high strength and light weight. In the present investigation, aluminium-alumina (Al–Al 2 O 3 ) nanocermets were developed. The nanocermets possessed the toughness of Al phase and the strength of reinforced nano Al 2 O 3 particles. Here, the nanocermets green bodies were prepared such that upon sintering each adjacent Al 2 O 3 grain was conformally coated with a continuous Al film. Mechanical milling via a combination of cryo-milling and ball-milling was used to prepare Al–Al 2 O 3 nano-powder mixture. Spark plasma sintering was used to consolidate the green nanocermets. Results show that specimens prepared from 35 h milled powders exhibited superior mechanical properties as compared to other nanocermets, with an improvement in microhardness, elastic modulus, compressive strength, and flexural strength by 88%, 78%, 69% and 40% with respect to sintered pure Al. The indentation toughness of 35 h milled nanocermet showed an increment of 12.5% with respect to sintered pure Al 2 O 3 . Characterization techniques such as XRD, SEM, HR-TEM, X-ray tomography were employed to investigate the phase evolution and structural morphology, while FEM simulations were performed for understanding the crack propagation mechanisms in nanocermets.

36 MATERIALS SCIENCE↗
DOE OSTI2017

The accurate particle tracer code

The Accurate Particle Tracer (APT) code is designed for systematic large-scale applications of geometric algorithms for particle dynamical simulations. Based on a large variety of advanced geometric algorithms, APT possesses long-term numerical accuracy and stability, which are critical for solving multi-scale and nonlinear problems. To provide a flexible and convenient I/O interface, the libraries of Lua and Hdf5 are used. Following a three-step procedure, users can efficiently extend the libraries of electromagnetic configurations, external non-electromagnetic forces, particle pushers, and initialization approaches by use of the extendible module. APT has been used in simulations of key physical problems, such as runaway electrons in tokamaks and energetic particles in Van Allen belt. As an important realization, the APT-SW version has been successfully distributed on the world’s fastest computer, the Sunway TaihuLight supercomputer, by supporting master–slave architecture of Sunway many-core processors. Here, based on large-scale simulations of a runaway beam under parameters of the ITER tokamak, it is revealed that the magnetic ripple field can disperse the pitch-angle distribution significantly and improve the confinement of energetic runaway beam on the same time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

QUIC-fire: A fast-running simulation tool for prescribed fire planning

Coupled fire-atmospheric modeling tools are increasingly used to understand the complex and dynamic behavior of wildland fires. Multiple research tools linking combustion to fluid flow use Navier-Stokes numerical solutions coupled to a thermodynamic model to understand fire-atmospheric feedbacks, but these computational fluid dynamics approaches require high-performance computing resources. We present a new simulation tool called QUIC-Fire to rapidly solve these feedbacks by coupling the mature 3-D rapid wind solver QUIC-URB to a physics-based cellular automata fire spread model Fire-CA. QUIC-Fire uses 3-D fuels inputs similar to the CFD model FIRETEC, allowing this tool to simulate effects of fuel structure on local winds and fire behavior. Results comparing fire behavior metrics to the computational fluid dynamic model FIRETEC show strong agreement. As a result, QUIC-Fire is the first tool intended to provide an opportunity for prescribed fire planners to compare, evaluate and design burn plans, including complex ignition patterns and coupled fire atmospheric feedbacks.

54 ENVIRONMENTAL SCIENCES↗