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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

ON THE ORIGIN OF HIGH-ALTITUDE OPEN CLUSTERS IN THE MILKY WAY

We present a dynamical study of the effect of the bar and spiral arms on the simulated orbits of open clusters in the Galaxy. Specifically, this work is devoted to the puzzling presence of high-altitude open clusters in the Galaxy. For this purpose we employ a very detailed observationally motivated potential model for the Milky Way and a careful set of initial conditions representing the newly born open clusters in the thin disk. We find that the spiral arms are able to raise an important percentage of open clusters (about one-sixth of the total employed in our simulations, depending on the structural parameters of the arms) above the Galactic plane to heights beyond 200 pc, producing a bulge-shaped structure toward the center of the Galaxy. Contrary to what was expected, the spiral arms produce a much greater vertical effect on the clusters than the bar, both in quantity and height; this is due to the sharper concentration of the mass on the spiral arms, when compared to the bar. When a bar and spiral arms are included, spiral arms are still capable of raising an important percentage of the simulated open clusters through chaotic diffusion (as tested from classification analysis of the resultant high-z orbits), but the bar seems to restrain them, diminishing the elevation above the plane by a factor of about two.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ENERGY DISSIPATION AND LANDAU DAMPING IN TWO- AND THREE-DIMENSIONAL PLASMA TURBULENCE

Plasma turbulence is ubiquitous in space and astrophysical plasmas, playing an important role in plasma energization, but the physical mechanisms leading to dissipation of the turbulent energy remain to be definitively identified. Kinetic simulations in two dimensions (2D) have been extensively used to study the dissipation process. How the limitation to 2D affects energy dissipation remains unclear. This work provides a model of comparison between two- and three-dimensional (3D) plasma turbulence using gyrokinetic simulations; it also explores the dynamics of distribution functions during the dissipation process. It is found that both 2D and 3D nonlinear gyrokinetic simulations of a low-beta plasma generate electron velocity-space structures with the same characteristics as that of the linear Landau damping of Alfvén waves in a 3D linear simulation. The continual occurrence of the velocity-space structures throughout the turbulence simulations suggests that the action of Landau damping may be responsible for the turbulent energy transfer to electrons in both 2D and 3D, and makes possible the subsequent irreversible heating of the plasma through collisional smoothing of the velocity-space fluctuations. Although, in the 2D case where variation along the equilibrium magnetic field is absent, it may be expected that Landau damping is not possible, a common trigonometric factor appears in the 2D resonant denominator, leaving the resonance condition unchanged from the 3D case. The evolution of the 2D and 3D cases is qualitatively similar. However, quantitatively, the nonlinear energy cascade and subsequent dissipation is significantly slower in the 2D case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

gevolution: a cosmological N-body code based on General Relativity

We present a new N-body code, gevolution , for the evolution of large scale structure in the Universe. Our code is based on a weak field expansion of General Relativity and calculates all six metric degrees of freedom in Poisson gauge. N-body particles are evolved by solving the geodesic equation which we write in terms of a canonical momentum such that it remains valid also for relativistic particles. We validate the code by considering the Schwarzschild solution and, in the Newtonian limit, by comparing with the Newtonian N-body codes Gadget-2 and RAMSES . We then proceed with a simulation of large scale structure in a Universe with massive neutrinos where we study the gravitational slip induced by the neutrino shear stress. The code can be extended to include different kinds of dark energy or modified gravity models and going beyond the usually adopted quasi-static approximation. Our code is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantum dot polarisation converter in an optomechanical cavity

We propose a scheme of a quantum photon polarisation converter, which is based on controlled electron – photon – phonon transitions in a hybrid semiconductor nanostructure. This structure consists of a GaAs/InAs quantum dot (QD) that has a parallelepiped shape and contains a single electron, and an optomechanical microcavity (MC) based on a photonic crystal (PC) that supports two orthogonally polarised photonic modes and one mechanical (phonon) mode. Within the framework of the microscopic theory, the QD and MC performance characteristics are found. Populations of states of the system as functions of time and its parameters are calculated. The principal possibility of photon polarisation conversion using transitions in a five-level resonance scheme for coherent (single-photon) and steady-state (subphoton) regimes is shown. The MC optical and mechanical spectra are simulated, and the PC structure parameters are selected to ensure the efficient operation of the converter. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

High-Entropy Metallic Glasses

This chapter applies the concept of high entropy to metallic glasses (MGs), in particular, to those in a bulk shape: bulk metallic glasses (BMGs). The resultant target materials in this chapter are mainly high-entropy bulk metallic glasses (HE-BMGs), which have recently been developed as alloys with characteristics of both high-entropy alloys (HEAs) and BMGs. Here, the contents in this chapter start by introducing historic background of HE-BMGs and by summarizing the differences between HEAs and BMGs. Then, the fundamental properties of representative HE-BMGs found to date are described mainly in terms of thermodynamic and mechanical behaviors. Besides the experiments, the latest computational approach for clarifying the features of HE-BMGs is described based on the results using ab initio molecular dynamics simulations for the atomic structure, chemical interaction, and diffusivity in this unique class of materials. The current status and future prospects of the HE-BMGs by utilizing their unique features are outlined for their future applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Full-wave modeling of RF waves in fusion plasmas with finite element method: Progress in past decades and its future role

This paper reviews the progress in computing radio frequency (RF) wave fields in fusion plasmas, specifically focusing on simulations utilizing the finite element method (FEM) over the past few decades. Computing RF wave fields in fusion plasmas presents unique challenges due to large simulation domains, complex antenna structures, non-local dielectric properties, and wide ranges of spatial scales. It highlights key developments and outlines future directions, primarily addressing waves in the ion cyclotron (IC) to lower hybrid (LH) frequency range. We begin with briefly revisiting earlier developments before the widespread availability of modern computer-aided engineering (CAE) software based on FEM. This historical perspective illuminates early progress and the physics difficulties that motivated ongoing work within the community. Modern wave simulations for RF antennas based on FEM are characterized by the use of detailed 3D antenna model geometry generated from engineering CAD software and localized wave dielectric model. Significant advancements have also been made in improving physics models to include phenomena such as RF sheath rectification and wave scattering. FEM-based RF simulations have also been applied to compute wave propagation in the core region, where the inclusion of non-local dielectric response is crucial. This is a challenging goal, and several promising approaches have been proposed in this area. Additionally, RF simulation development initiatives based on open-source libraries have gained popularity, demonstrating scalability and flexibility in extending physics models. This paper will discuss the advantages and disadvantages of using such a publicly available FEM library.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Multi-sightline observation of narrow absorption lines in lensed quasar SDSS J1029+2623

We exploit the widely separated images of the lensed quasar SDSS J1029+2623 (z{sub em} = 2.197, θ = 22.″5) to observe its outflowing wind through two different sightlines. We present an analysis of three observations, including two with the Subaru telescope in 2010 February and 2014 April, separated by four years, and one with the Very Large Telescope, separated from the second Subaru observation by ∼2 months. We detect 66 narrow absorption lines (NALs), of which 24 are classified as intrinsic NALs that are physically associated with the quasar based on partial coverage analysis. The velocities of intrinsic NALs appear to cluster around values of v{sub ej} ∼ 59,000, 43,000, and 29,000 km s{sup −1}, which is reminiscent of filamentary structures obtained by numerical simulations. There are no common intrinsic NALs at the same redshift along the two sightlines, implying that the transverse size of the NAL absorbers should be smaller than the sightline distance between two lensed images. In addition to the NALs with large ejection velocities of v{sub ej} > 1000 km s{sup −1}, we also detect broader proximity absorption lines (PALs) at z{sub abs} ∼ z{sub em}. The PALs are likely to arise in outflowing gas at a distance of r ≤ 620 pc from the central black hole with an electron density of n {sub e} ≥8.7 × 10{sup 3} cm{sup −3}. These limits are based on the assumption that the variability of the lines is due to recombination. We discuss the implications of these results on the three-dimensional structure of the outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Machine learning bridges local static structure with multiple properties in metallic glasses

A long-standing challenge in the metallic glass (MG) community has been how to quantitatively gauge the influence of the intricate local packing environment on the response (such as the propensity for atomic rearrangement) of the atomic configuration to external stimuli. Here we establish this structure–property relation by representing the complex amorphous structure using a single, flexibility-orientated structural quantity. This structural flexibility (SF) couples to a bona fide structural representation, the pair distribution function (PDF) of individual atoms, through a weighting function that reflects what matters in the static atomic configuration to dynamic responses. Machine learning is used, employing microscopic flexibility volume as the supervisory signal, to establish via direct regression an optimized weighting vector, which is proven robust for all quenching rates, deformation conditions, and different compositions in a given (e.g., Cu x Zr 100-x ) alloy system. Additionally, the SF is evaluated solely from the particle positions (PDF), for any structure variation, from the atomic scale up to sample average. Strong correlations are demonstrated between SF and a broad range of properties, including vibrational, diffusional, as well as elastic and plastic relaxation responses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of 3D-printed surface textures on wear mechanism in 3-body abrasion of soil

This study systematically investigates the enhancement of wear resistance in 3D printed surface textures through both experimental and theoretical approaches. Three distinct surface morphologies (Smooth Surface, Surface with uniformly distributed Pits, and Surface with uniformly distributed Bumps) were fabricated using High-Impact Polystyrene, where the meso-scale textures were precisely controlled through the 3D printing process. Wear behavior was evaluated using a 3-body wear tester in an abrasive particle environment, analyzing the influence of surface textures under various operating conditions. Systematic wear tests revealed that optimally designed surface textures achieved a remarkable 77 % reduction in wear compared to the worst-performing sample. The wear mechanisms were comprehensively characterized through weight loss measurements, Scanning Electron Microscopy (SEM), and Energy Dispersive Spectroscopy (EDS) analyses, elucidating the surface morphology changes and their interaction with wear particles. Notably, the study identified how the geometric characteristics of surface textures influence the movement of wear particles and the distribution of contact stresses. Discrete element method simulations corroborated the experimental findings, providing theoretical validation for the enhanced wear resistance of the optimal structure. The high correlation between simulated wear patterns and experimental results validates the reliability of the proposed design methodology. In conclusion, these results demonstrate that 3D printed surface texturing offers a cost-effective and scalable approach to significantly improve wear resistance in engineering applications, presenting a practical alternative to conventional, high-cost surface engineering methods.

3-body abrasion↗
DOE OSTI2021

Exemption of lattice collapse in Ni–MnO 2 birnessite regulated by the structural water mobility

Lattice collapse and associated mechanical fracture frequently occur in Li-intercalated metal oxide cathodes at the deep charge state upon Li-ion removal, governed by chemical compositions and the resulting electron density of the oxygen atoms. However, similar lattice collapse for metal oxide electrodes in aqueous storage and its mitigation have not been well studied. In this work, we reported the lattice collapse of MnO 2 layered birnessite during the aqueous de-sodiation process at high voltage due to the structural water motion, as evidenced by in situ XRD. Unlike non-aqueous Li-intercalated electrodes, Ni-dopants mitigated the lattice collapse of birnessite at deep charge states. Moreover, density functional theory (DFT) calculations showed that Ni doping induces charge depletion of lattice oxygen due to its higher electronegativity than Mn. This charge reduction, in turn, yields a significant decrease in electrostatic repulsion between the oxygens belonging to lattice and structural water. Classical molecular dynamics simulations based on atomic charges obtained from DFT elucidate that Ni-doping facilitates immobilization of structural water in (Ni)MnO 2 and prevents lattice collapse upon Na-ion removal. (Ni)MnO 2 exempted from lattice collapse shows an improved storage capacity relative to MnO 2 while maintaining similar cycling stability.

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

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Tree tensor network approach to simulating Shor's algorithm

Constructively simulating quantum systems furthers our understanding of qualitative and quantitative features which may be analytically intractable. In this paper, we directly simulate and explore the entanglement structure present in the paradigmatic example for exponential quantum speedups: Shor's algorithm. To perform our simulation, we construct a dynamic tree tensor network which manifestly captures two salient circuit features for modular exponentiation. These are the natural two-register bipartition and the invariance of entanglement with respect to permutations of the top-register qubits. Our construction help identify the entanglement entropy properties, which we summarize by a scaling relation. Further, the tree network is efficiently projected onto a matrix product state from which we efficiently execute the quantum Fourier transform. In conclusion, future simulation of quantum information states with tensor networks exploiting circuit symmetries is discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Effects of Si solute on the glass formation and atomic structure of Pd liquid

Molecular dynamics simulations were performed to study the effects of Si solute on the glass formation and crystallization of Pd liquid. Pure Pd sample prepared by quenching process with cooling rate of 10 13 K/s can be in an amorphous state but the structural analysis indicates there is nearly no glass-forming motif in the sample. However, doping a small amount Si (Si concentration ~ 4%) the sample can be vitrified at a cooling rate of 10 12 K/s. The glass-forming motifs such as Pd-centered Z13, Si-centered Z9-like and Mixed-ICO-Cube clusters with 5-fold local symmetry are found to be the dominant short-range orders in the glassy samples. With the increasing of the Si-doping concentration, these glass-forming motifs tend to aggregate and connect with each other forming a network structure. Our calculated results revealed that Si solutes in liquid Pd can significantly enhance the glass-forming ability.

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

Davis Computational Spectroscopy Workflow—From Structure to Spectra

Here, we describe an automated workflow that connects a series of atomic simulation tools to investigate the relationship between atomic structure, lattice dynamics, materials properties, and inelastic neutron scattering (INS) spectra. Starting from the atomic simulation environment (ASE) as an interface, we demonstrate the use of a selection of calculators, including density functional theory (DFT) and density functional tight binding (DFTB), to optimize the structures and calculate interatomic force constants. We present the use of our workflow to compute the phonon frequencies and eigenvectors, which are required to accurately simulate the INS spectra in crystalline solids like diamond and graphite as well as molecular solids like rubrene. We have also implemented a machine-learning force field based on Chebyshev polynomials called the Chebyshev interaction model for efficient simulation (ChIMES) to improve the accuracy of the DFTB simulations. We then explore the transferability of our DFTB/ChIMES models by comparing simulations derived from different training sets. We show that DFTB/ChIMES demonstrates ~100× reduction in computational expense while retaining most of the accuracy of DFT as well as yielding high accuracy for different materials outside of our training sets. The DFTB/ChIMES method within the workflow expands the possibilities to use simulations to accurately predict materials properties of increasingly complex structures that would be unfeasible with ab initio methods.

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

Effective work function engineering for a TiN/XO(X = La, Zr, Al)/SiO{sub 2} stack structures

In this study, we demonstrated that work function engineering is possible over a wide range (+200 mV to −430 mV) in a TiN/XO (X = La, Zr, or Al)/SiO{sub 2} stack structures. From ab initio simulations, we selected the optimal material for the work function engineering. The work function engineering mechanism was described by metal diffusion into the TiN film and silicate formation in the TiN/SiO{sub 2} interface. The metal doping and the silicate formation were confirmed by transmission electron microscopy and energy dispersive spectroscopy line profiling, respectively. In addition, the amount of doped metal in the TiN film depended on the thickness of the insertion layer XO. From the work function engineering technique, which can control a variety of threshold voltages (Vth), an improvement in transistors with different V{sub th} values in the TiN/XO/SiO{sub 2} stack structures is expected.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Local ordering in Ge/Ge–Sn semiconductor alloy core/shell nanowires revealed by extended x-ray absorption fine structure (EXAFS)

Short-range atomic order in semiconductor alloys is a relatively unexplored topic that may promote design of new materials with unexpected properties. Here, local atomic ordering is investigated in Ge–Sn alloys, a group-IV system that is attractive for its enhanced optoelectronic properties achievable via a direct gap for Sn concentrations exceeding ≈10 at. %. The substantial misfit strain imposed on Ge–Sn thin films during growth on bulk Si or Ge substrates can induce defect formation; however, misfit strain can be accommodated by growing Ge–Sn alloy films on Ge nanowires, which effectively act as elastically compliant substrates. In this work, Ge core/Ge 1–x Sn x (x ≈ 0.1) shell nanowires were characterized with extended x-ray absorption fine structure (EXAFS) to elucidate their local atomic environment. Simultaneous fitting of high-quality EXAFS data collected at both the Ge K-edge and the Sn K-edge reveals a large (≈ 40%) deficiency of Sn in the first coordination shell around a Sn atom relative to a random alloy, thereby providing the first direct experimental evidence of significant short-range order in this semiconductor alloy system. Furthermore, comparison of path length data from the EXAFS measurements with density functional theory simulations provides alloy atomic structures consistent with this conclusion.

36 MATERIALS SCIENCE↗