Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Materializing rival ground states in the barlowite family of kagome magnets: quantum spin liquid, spin ordered, and valence bond crystal states

The spin- 1 2 kagome antiferromagnet is considered an ideal host for a quantum spin liquid ground state. We find that when the bonds of the kagome lattice are modulated with a periodic pattern, new quantum ground states emerge. Newly synthesized crystalline barlowite (Cu 4 (OH) 6 FBr) and Zn-substituted barlowite demonstrate the delicate interplay between singlet states and spin order on the spin- 1 2 kagome lattice. Comprehensive structural measurements demonstrate that our new variant of barlowite maintains hexagonal symmetry at low temperatures with an arrangement of distorted and undistorted kagome triangles, for which numerical simulations predict a pinwheel va- lence bond crystal (VBC) state instead of a quantum spin liquid (QSL). The presence of interlayer spins eventually leads to an interesting pinwheel q=0 magnetic order. Partially Zn-substituted bar- lowite (Cu 3.44 Zn 0.56 (OH) 6 FBr) has an ideal kagome lattice and shows QSL behavior, indicating a surprising robustness of the QSL against interlayer impurities. The magnetic susceptibility is simi- lar to that of herbertsmithite, even though the Cu 2+ impurities are above the percolation threshold for the interlayer lattice and they couple more strongly to the nearest kagome moment. This system is a unique playground displaying QSL, VBC, and spin order, furthering our understanding of these highly competitive quantum states.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Wind plants can impact long-term local atmospheric conditions

Abstract Long-term weather and climate observatories can be affected by the changing environments in their vicinity, such as the growth of urban areas or changing vegetation. Wind plants can also impact local atmospheric conditions through their wakes, characterized by reduced wind speed and increased turbulence. We explore the extent to which the wind plants near an atmospheric measurement site in the central United States have affected their long-term measurements. Both direct observations and mesoscale numerical weather prediction simulations demonstrate how the wind plants induce a wind deficit aloft, especially in stable conditions, and a wind speed acceleration near the surface, which extend $$\sim 30$$ ∼ 30 km downwind of the wind plant. Turbulence kinetic energy is significantly enhanced within the wind plant wake in stable conditions, with near-surface observations seeing an increase of more than 30% a few kilometers downwind of the plants.

17 WIND ENERGY↗
DOE OSTI2025

Modeling the impact of structure and coverage on the reactivity of realistic heterogeneous catalysts

Adsorbates often cover the surfaces of catalysts densely as they carry out reactions, dynamically altering their structure and reactivity. Understanding adsorbate-induced phenomena and harnessing them in our broader quest for improved catalysts is a substantial challenge that is only beginning to be addressed. Here, in this work, we chart a path toward a deeper understanding of such phenomena by focusing on emerging in silico modeling methodologies, which will increasingly incorporate machine learning techniques. We first examine how adsorption on catalyst surfaces can lead to local and even global structural changes spanning entire nanoparticles, and how this affects their reactivity. We then evaluate current efforts and the remaining challenges in developing robust and predictive simulations for modeling such behavior. Last, we provide our perspectives in four critical areas—integration of artificial intelligence, building robust catalysis informatics infrastructure, synergism with experimental characterization, and adaptive modeling frameworks—that we believe can help surmount the remaining challenges in rationally designing catalysts in light of these complex phenomena.

catalytic mechanisms↗
DOE OSTI2021

Effect of water frustration on water oxidation catalysis in the nanoconfined interlayers of layered manganese oxides birnessite and buserite

The role of geometric frustration of water molecules on the rate of water oxidation in the nanoconfined interlayer of manganese-oxide layered materials (birnessite, buserite) is examined in a well-controlled experiment. Calcium buserite is prepared, and used in a split-batch synthetic protocol to prepare calcium birnessite, sodium buserite, and sodium birnessite, and partially dehydrated sodium birnessite. Thus, four samples are prepared in which features effecting catalytic efficiency (defect density, average manganese oxidation state) are controlled, and the main difference is the degree of hydration of the interlayer (two layers of water in buserites vs one layer of water in birnessite). Here, molecular dynamics simulations predict birnessite samples to exhibit geometric water frustration, which facilitates redox catalysis by lowering the Marcus reorganizaiton energy of electron transfer, while buserite samples exhibit traditional intermolecular hydrogen bonding among the two-layer aqeuous region, leading to slower catalytic behavior akin to redox reactions in bulk water. Water oxdiation activity is investigated using chemical and electrochemical techniques, demonstrating and quantifying the role of water frustration in enhancing catalysis. Calculation and experiment demonstrate dehydrated sodium birnessite to be most effective, and calcium buserite the least effective, with a difference in electrocatlytic overpotential of ~750 mV and a ~20-fold difference in turnover number.

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

Candidate strongly lensed type Ia supernovae in the Zwicky Transient Facility archive

Gravitationally lensed type Ia supernovae (glSNe Ia) are unique astronomical tools that can be used to study cosmological parameters, distributions of dark matter, the astrophysics of the supernovae, and the intervening lensing galaxies themselves. A small number of highly magnified glSNe Ia have been discovered by ground-based telescopes such as the Zwicky Transient Facility (ZTF), but simulations predict that a fainter, undetected population may also exist. We present a systematic search for glSNe Ia in the ZTF archive of alerts distributed from June 1 2019 to September 1 2022. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry (with constraints on the peak absolute magnitude and the distance to a catalogue-matched galaxy, as examples) before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. The candidates were also cross-matched with the DESI spectroscopic catalogue. Seven transients were identified that passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. Although superluminous supernovae (SLSNe) cannot be fully rejected as contaminants, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 ± 3 and 34 ± 1 days for the two events, respectively, which suggests that we have uncovered a population of glSNe Ia with longer time delays. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, such as the Rubin Observatory Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Eulerian finite element simulations of the drop weight impact test with a dislocation Density-based continuum model

During a drop weight impact test, the kinetic energy of the falling weight is transferred into the sample resting on the anvil. The plastic deformation in the sample is an important mechanism for the dissipation of the input kinetic energy. We use Eulerian finite element analysis to simulate the deformation and temperature evolution in a 1,3,5-trinitro-1,3,5-triazine (RDX) sample consisting of multiple crystals. In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, i.e., advection. In an effort to reduce the advection errors, we use a rate form of a dislocation density-based continuum model by Luscher et al. Here, the simulations predict localization of plastic deformation, and plastic dissipation as a significant source of heat generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A 1D coupled transport/cold plasma wave model for parallel ponderomotive density modification near RF actuators

The plasma density profile close to RF actuators plays a vital role in the efficiency of power coupling from the antenna to the core plasma of fusion devices. Nonlinear interactions can be significant when considering the interaction of applied RF power and the tenuous SOL plasma, and the inclusion of this region in computational analysis is essential to providing robust predictive simulations of plasma response to RF application. Of particular interest is the observed density modification during RF operation. Here we present a self consistent description of the nonlinear ponderomotive force in the antenna near field region, via coupling a 1D field aligned transport model with a 1D cold plasma wave model. Results are presented, which show a parallel electric field of 300 Vm -1 drove a density reduction of ~ 20% over 25 RF periods.

Barnett, Rhea↗
DOE OSTI2021

Liquid metallic laser ultrasound transducer for high-temperature applications

This study aims to investigate a laser ultrasound (LUS) transducer for high-temperature (>100 °C) applications. For decades, many researchers have investigated efficient LUS transducers, yet studies on laser ultrasound transducers capable of performing at the high-temperature condition are rarely found in the literature. Most current LUS transducers still utilize a polymer-based composite material, that is, not stable at varying temperature conditions. This study introduces a liquid metallic (LM) LUS transducer that utilizes field's metal, which has a high thermal expansion (~3 × 10 -4 K -1 ). We hypothesized that such a high thermal expansion of the liquid metal can effectively produce laser-generated ultrasound waves, substituting for conventional polymer-based transducers. A numerical simulation predicted that the LM LUS transducer would produce higher LUS intensity (~22 dB) than that without the LUS transducer. Experiment results confirmed that the LM transducer effectively intensifies the ultrasound wave signals, obtaining a signal-to-noise gain over 30 dB. Moreover, the transducer was found capable of transmitting detectable wave packets in relatively high-temperature conditions (~400 °C), while conventional candle soot nanoparticle-polydimethylsiloxane could not perform stably at these elevated temperatures. The investigations introduced in this article are scientifically significant since we demonstrated the engineering feasibility of liquid metallic materials for LUS transducers.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Challenges and transformative opportunities in superconductor vortex physics

In superconductors, the motion of vortices introduces unwanted dissipation that is disruptive to applications. Fortunately, material defects can immobilize vortices, acting as vortex pinning centers, which engenders dramatic improvements in superconductor material properties and device operation. This has motivated decades of research into developing methods of tailoring the disorder landscape in superconductors to increase the strength of vortex pinning. Yet, efficacious materials engineering still eludes us. The electromagnetic properties of real (disordered) superconducting materials cannot yet be reliably predicted, such that designing superconductors for applications remains a largely inefficient process of trial and error. This is ultimately due to large gaps in our knowledge of vortex dynamics: the field is challenged by the extremely complex interplay between vortex elasticity, vortex-vortex interactions, and material disorder. Here, we review obstacles and recent successes in understanding and controlling vortex dynamics in superconducting materials and devices. We further identify major open questions and discuss opportunities for transformative research in the field. This includes improving our understanding of vortex creep, determining and reaching the ceiling for the critical current, advanced microscopy to garner accurate structure-property relationships, frontiers in predictive simulations and the benefits of artificial intelligence, as well as controlling and exploiting vortices in quantum information applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling of small tungsten dust grains in EAST tokamak with NDS-BOUT ++

In order to investigate the transport of small dusts as well as their evolution property along their trajectories, the NDS module is developed under the BOUT++ framework, a highly desirable C++ code package to perform parallel plasma fluid simulations with an arbitrary number of equations in three-dimensional curvilinear coordinates. Due to the severe dust ablation in fusion plasmas, the dust size would decrease from micrometer to nanometer, resulting in impurities. Small dusts in the simulations here are specified as tungsten spheres with the radii on or below the order of submicrometer. The Rayleigh limit is included in the charging process when the dust is ablated to the droplet phase. The simulation results from the NDS module show that a 200 nm radius spherical tungsten dust originated from upper divertor region of EAST Tokamak is ablated completely due to the intense heating from the incoming plasma inside the core region, well consistent with the CCD footage of EAST shot # 81459. Furthermore it is found that the magnetic field dominates the dust transport when the dust radius is below 100 nm during the ablation along the trajectory. Our simulations predict that a 10 nm radius spherical tungsten dust injected from the inner midplane is well constrained by the magnetic field, and it reaches the inner divertor target with a velocity on the order of km/s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Negative capacitance regime in antiferroelectric PbZrO 3

The potential of antiferroelectrics to exhibit a negative capacitance regime has been largely overlooked as all the attention focused on their polar counterparts, ferroelectrics. We use nonequilibrium first-principles-based molecular dynamics to probe a negative capacitance regime in prototypical antiferroelectric PbZrO 3 . Simulations predict that this antiferroelectric can exhibit a negative capacitance/susceptibility regime in response of polarization to an internal electric field, which is a superposition of applied and residual depolarizing fields. Consequently, the regime emerges when the polarization surface charge in the polar phase of antiferroelectric is not fully screened, as is often the case in thin films and nanostructures. The negative capacitance regime occurs below the Curie temperature and disappears in the paraelectric phase. We find that the time the material spends in the negative capacitance regime is proportional to the time needed to complete antipolar–polar (or its reverse) transition and shortens as the frequency of the applied field increases. Furthermore, a negative susceptibility value exhibits strong dependence on the quality of surface charge screening with the largest in magnitude values occurring in the vicinity of the transition into a negative capacitance regime.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Carrier capture and emission by substitutional carbon impurities in GaN vertical diodes

A model was developed for the operation of a GaN pn junction vertical diode which includes rate equations for carrier capture and thermally activated emission by substitutional carbon impurities and carrier generation by ionizing radiation. The model was used to simulate the effect of ionizing radiation on the charge state of carbon. These simulations predict that with no applied bias, carbon is negatively charged in the n-doped layer, thereby compensating n-doping as experimentally observed in diodes grown by metal-organic chemical vapor deposition. With reverse bias, carbon remains negative in the depletion region, i.e., compensation persists in the absence of ionization but is neutralized by exposure to ionizing radiation. This increases charge density in the depletion region, decreases the depletion width, and increases the capacitance. The predicted increase in capacitance was experimentally observed using a pulsed 70 keV electron beam as the source of ionization. In additional confirming experiments, the carbon charge-state conversion was accomplished by photoionization using sub-bandgap light or by the capture of holes under forward bias.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Demonstration of improved laser preheat with a cryogenically cooled magnetized liner inertial fusion platform

We report on progress implementing and testing cryogenically cooled platforms for Magnetized Liner Inertial Fusion (MagLIF) experiments. Two cryogenically cooled experimental platforms were developed: an integrated platform fielded on the Z pulsed power generator that combines magnetization, laser preheat, and pulsed-power-driven fuel compression and a laser-only platform in a separate chamber that enables measurements of the laser preheat energy using shadowgraphy measurements. The laser-only experiments suggest that ~89% ± 10% of the incident energy is coupled to the fuel in cooled targets across the energy range tested, significantly higher than previous warm experiments that achieved at most 67% coupling and in line with simulation predictions. The laser preheat configuration was applied to a cryogenically cooled integrated experiment that used a novel cryostat configuration that cooled the MagLIF liner from both ends. Furthermore, the integrated experiment, z3576, coupled 2.32 ± 0.25 kJ preheat energy to the fuel, the highest to-date, demonstrated excellent temperature control and nominal current delivery, and produced one of the highest pressure stagnations as determined by a Bayesian analysis of the data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

Here, we develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

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

Reduced instability growth and improved radiation trapping with optimized shock timing in double-shell inertial confinement fusion capsules

The double shell is a volume-burn inertial confinement fusion concept consisting of two concentric shells: a low-Z outer shell that collides with and transfers momentum to a high-Z inner shell which compresses and heats the thermonuclear fuel. The increased number of capsule interfaces and severe hydrodynamic instability of the high-density pusher during its acceleration phase provide challenges to the success of the double shell. Two-dimensional radiation-hydrodynamics simulations predict the hydrodynamic instability growth on the outer surface of the pusher can be greatly reduced through appropriate timing of two shocks that cross this interface. One of these shocks, unique to multi-shell designs, arises from radiation-driven ablation of the inner shell ahead of the main shock, the second shock of concern. The shock timing is optimized by increasing the thickness of a low-Z tamper layer exterior to the pusher, resulting in only minimal changes to the implosion timing. Reducing the instability growth on the outer surface of the high-Z pusher can dramatically decrease the modulations that feedthrough to the pusher inner surface, improving the efficacy of radiation trapping in the thermonuclear fuel and increasing the predicted thermonuclear yield by ≳20×.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Numerical simulation of an atmospheric pressure plasma jet with coaxial shielding gas

A two-dimensional ( r, z ) numerical simulation of the discharge characteristics of an atmospheric pressure plasma jet (APPJ), with coaxial shielding gas, was performed. The helium working gas flowed in a central capillary tube, engulfing a needle electrode powered by 13.7 MHz radio frequency sinusoidal voltage. The N 2 shielding gas flowed in the annular space of a coaxial tube. These gases emerged, in laminar flow, in a 78%N 2 -21%O 2 -1%Ar dry air ambient. The characteristics of the APPJ with shielding gas were compared to those of the APPJ without shielding gas. The nitrogen shielding gas hindered the diffusion of oxygen and argon from the ambient air into the helium jet. With the shielding gas present, more nitrogen penetrated into the helium core, causing a shorter plasma ‘plume’. The flow rates of the working and shielding gas, critically affected the gas temperature, and in turn the discharge characteristics. For a He flow of 2 standard liters per minute (slm), switching on the nitrogen shielding gas flow (at 4.5 slm) reduced the on-axis O 2 and Ar mole fractions from 3.9 × 10 - 4 to 6.8 × 10 - 5 and from 1.9 × 10 - 5 to 3.3 × 10 - 6 , respectively, at an axial distance of 3 mm downstream of the nozzle. The radial profiles of the mole fractions of the ambient gases were monotonically and strongly decreasing towards the system axis, for short axial distances from the nozzle (~1 mm), but became progressively flatter at longer distances from the nozzle (3 mm and 5 mm). Simulation predictions captured the salient features of experimental data of ambient species mole fractions in the plasma jet, and the 706 nm optical emission intensity profiles of the He 3 3 S excited state.

Physics↗
DOE OSTI2024

Effect of helium bubbles on the mobility of edge dislocations in copper

Abstract Helium bubbles can form in materials upon exposure to irradiation. It is well known that the presence of helium bubbles can cause changes in the mechanical behavior of materials. To improve the lifetime of nuclear components, it is important to understand deformation mechanisms in helium-containing materials. In this work, we investigate the interactions between edge dislocations and helium bubbles in copper using molecular dynamics (MD) simulations. We focus on the effect of helium bubble pressure (equivalently, the helium-to-vacancy ratio) on the obstacle strength of helium bubbles and their interaction with dislocations. Our simulations predict significant differences in the interaction mechanisms as a function of helium bubble pressure. Specifically, bubbles with high internal pressure are found to exhibit weaker obstacle strength as compared to low-pressure bubbles of the same size due to the formation of super-jogs in the dislocation. Activation energies and rate constants extracted from the MD data confirm this transition in mechanism and enable upscaling of these phenomena to higher length-scale models.

36 MATERIALS SCIENCE↗