Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator study”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Facets control charge separation during photoelectrochemical water oxidation with strontium titanate (SrTiO 3 ) single crystals

The photocatalytic overall water splitting reaction provides a pathway to hydrogen fuel from sunlight. Photocatalysts must achieve the reaction without the application of an external bias, which requires an effective charge separation mechanism. Photolabeling studies and electrostatic simulations for the well-known CoOOH/Al:SrTiO 3 /Rh/Cr 2 O 3 photocatalyst suggest that charge separation is driven by work function differences at the (100) and (110) facets of SrTiO 3 , which are electron and hole selective, respectively. Here we use hydrogen annealed SrTiO 3–x single crystals to obtain the first quantitative assessment of the charge separation ability of the (100), or (110), or (111) facets during oxygen evolution. Under UV illumination (60 mW cm –2 ), the crystals exhibit variable water oxidation photocurrents (0.34, 0.82, 1.36 mA cm –2 at 1.23 V versus RHE) and photovoltage values of 1.40, 1.52 and 1.52 V for (100), (110), and (111) SrTiO 3–x , respectively. A surface photovoltage increase in that same order (0.31 V < 0.57 V < 0.67 V) is confirmed independently with vibrating Kelvin probe surface photovoltage spectroscopy (VKP-SPV) under 375 nm (1.91 mW cm –2 ) illumination. Mott Schottky measurements in aqueous K 3/4 [Fe(CN) 6 ] reveal facet-dependent flatband positions of –0.58, –0.71, and –0.74 V RHE for the (100), (110), and (111) crystals respectively. This confirms that the photoelectrochemical water oxidation performance of SrTiO 3–x crystals is controlled by the work function of each facet, which determines the electron transfer barrier height of the respective solid–liquid junctions. After correcting for differences in electron donor concentrations, barriers are found to increase in the order (100) < (111) < (110) and differ by as much as 0.16 eV, similar to an earlier prediction. Altogehter, these results explain the charge separation mechanism in SrTiO 3 photocatalysts and highlight the need for faceted semiconductor crystals as light absorbers in particle-based photocatalysts.

08 HYDROGEN↗
DOE OSTI2025

Solvation directed morphological control in metal oxide nanostructures

The development of structural hierarchy on various length scales during the crystallization process is ubiquitous in biological systems and minerals and is common in synthetic nanomaterials. The driving forces for the formation of complex architectures range from local interfacial interactions, that modify interfacial speciation, local supersaturation, and nucleation barriers, to macroscopic interparticle forces. Although it is enticing to interpret the formation of hierarchical architectures as the assembly of independently nucleated building blocks, crystallization pathways often follow monomer-by-monomer addition with structural complexity arising from interfacial chemical coupling and strongly correlated fluctuation dynamics in the electric double layers. Here, we show that the development of structural hierarchy through heterogeneous nucleation is driven by dipolar and solvation forces. Specifically, coupled simulations and experimental studies revealed that dipole build-up along the slow growth direction can trigger twinning and the development of branched architectures. Enthalpic solvation interactions were shown to either enhance or reduce the dipole moment of the nanoparticles and, thereby, control crystal morphology and architecture. The systematic studies of chemical coupling between different solvents and undercoordinated surface atoms of the growing nanocrystals revealed the mechanism of dimensionality control and the development of structural hierarchy without ligands or structure-directing agents.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Structure and dynamics of ethane confined in silica nanopores in the presence of CO 2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C 2 H 6 ) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C 2 H 6 and CO 2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C 2 H 6 , C 2 H 6 :CO 2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C 2 H 6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C 2 H 6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO 2 molecules are more strongly attracted to the pore surface compared to C 2 H 6 . The C 2 H 6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO 2 molecules play in enhancing C 2 H 6 diffusion (“molecular lubrication”) at high CO 2 :C 2 H 6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

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

Asymmetry between deuterium and tritium turbulent particle flows

The asymmetry between deuterium (D) and tritium (T) turbulent particle fluxes in mixed D–T plasmas is studied with numerical simulations of nonlinear gyrokinetic turbulence in ion temperature gradient-dominated and trapped electron mode-dominated regimes. At 50–50 D–T concentration, the asymmetry, or flow separation, between D and T fluxes is such that the tritium is better confined than the deuterium in both regimes. To supplement the nonlinear simulations, an analytic quasilinear theory of the particle flux symmetry breaking is developed and is valid for general electron dynamics. Furthermore, this theory correctly predicts the numerically computed deviation in the ion density fraction (from 50–50), or ion density gradient (from the electron gradient), required to restore equal deuterium and tritium fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Three-dimensional plasmoid-mediated reconnection and turbulence in Hall magnetohydrodynamics

Plasmoid instability accelerates reconnection in collisional plasmas by transforming a laminar reconnection layer into numerous plasmoids connected by secondary current sheets in two dimensions (2D) and by fostering self-generated turbulent reconnection in three dimensions (3D). In large-scale astrophysical and space systems, plasmoid instability likely initiates in the collisional regime but may transition into the collisionless regime as the fragmentation of the current sheet progresses toward kinetic scales. Hall magnetohydrodynamics (MHD) models are widely regarded as a simplified yet effective representation of the transition from collisional to collisionless reconnection. However, plasmoid instability in 2D Hall MHD simulations often leads to a single-X-line reconnection configuration, which significantly differs from fully kinetic particle-in-cell simulation results. This study shows that single-X-line reconnection is less likely to occur in 3D compared to 2D. Moreover, depending on the Lundquist number and the ratio between the system size and the kinetic scale, Hall MHD can also realize 3D self-generated turbulent reconnection. We analyze the features of the self-generated turbulent state, including the energy power spectra and the scale dependence of turbulent eddy anisotropy.

Physics↗
DOE OSTI2020

Influence of saturable absorber saturation power, modulation depth and relaxation time on pulse parameters of a soliton fibre laser

The influence of the saturable absorber saturation power, modulation depth and relaxation time on the energy and shape of fibre soliton laser pulses has been studied by numerical simulation. The results suggest that the use of a saturable absorber with a lower saturation power in a fibre laser ensures pulsed lasing onset at a lower gain. The fibre laser pulse energy decreases as the saturation power of the absorber or the relaxation time of its open state increases at a constant gain coefficient. In addition, increasing the absorber saturation modulation depth at a constant gain coefficient has been shown to reduce the pulse peak power and energy. (control of laser radiation parameters)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Perpendicular crossing chains enable high mobility in a noncrystalline conjugated polymer

The nature of interchain π-system contacts, and their relationship to hole transport, are elucidated for the high-mobility, noncrystalline conjugated polymer C16-IDTBT by the application of scanning tunneling microscopy, molecular dynamics, and quantum chemical calculations. The microstructure is shown to favor an unusual packing motif in which paired chains cross-over one another at near-perpendicular angles. By linking to mesoscale microstructural features, revealed by coarse-grained molecular dynamics and previous studies, and performing simulations of charge transport, it is demonstrated that the high mobility of C16-IDTBT can be explained by the promotion of a highly interconnected transport network, stemming from the adoption of perpendicular contacts at the nanoscale, in combination with fast intrachain transport.

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

Statistical Modeling and Analysis of k -Layer Coverage of Two-Dimensional Materials in Inkjet Printing Processes

Two-dimensional layered materials/flakes, also known as crystalline atom-thick layer nanosheets, have recently been receiving great attention in electronics fabrication due to their unique and intriguing properties. The k-layer coverage area (i.e., the area covered by k number of overlapping layers) of the printed flake pattern significantly impacts on the properties of the printed electronics. In this work, we constructed a statistical model to describe the k-layer coverage of randomly distributed two-dimensional materials. A series of results are obtained to provide not only the expectation but also the variance of the coverage area. The boundary effects on the random flakes coverage are also studied. In addition, an approximated statistical testing approach is also developed in this work to detect abnormal coverage patterns. Furthermore, the case studies based on simulated data and real flakes images obtained from the inkjet printing process demonstrate the accuracy and effectiveness of the proposed model and analysis methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spectral transition of multiscale turbulence in the tokamak pedestal

The transition in the turbulence spectrum from ion-scale dominated regimes to multiscale transport regimes that couple ion and electron scales is studied with gyrokinetic simulations of turbulent transport. The simulations are based on DIII-D high-confinement mode (H-mode) plasma parameters in the tokamak pedestal. The transition is initiated by varying the ion temperature gradient. To our knowledge, no full multiscale simulations of pedestal-like transport have been done previously. The experimental parameters lie in a bifurcation region between the two regimes. At long wavelengths, a complex, ion-direction hybrid mode is the dominant linearly unstable drift wave, while an electron temperature gradient-driven mode is unstable at short wavelengths. In the transition from the multiscale branch to the ion-scale branch, the magnitude of the ion-scale poloidal wavenumber spectrum of the nonlinear turbulent energy flux increases and the magnitude of the high-wavenumber spectrum decreases. Here, the decrease in the electron-scale transport is due to nonlinear mixing with ion-scale fluctuations and the ion-scale-driven zonal flows. A shift in the total energy associated with the fluctuating electrostatic potential intensity from dominantly drift kinetic energy in the multiscale regime to dominantly potential intensity in the ion-scale regime is well-correlated with the trend in the total energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Planar laser-induced fluorescence system for spatiotemporal ion velocity distribution function measurements

In this work, we present a planar laser-induced fluorescence (PLIF) system for measuring two-dimensional (2D), spatiotemporally resolved ion velocity distribution functions (IVDFs). A continuous-wave tunable diode laser produces a laser sheet that irradiates the plasma, and the resulting fluorescence is captured by an intensified CCD (ICCD) camera. Fluorescence images recorded at varying laser wavelengths are converted into 2D IVDFs using the Doppler shift principle. The developed diagnostic is implemented in an electron beam generated E × B plasma with a bulk plasma density of $\sim\!\!{10^{10}}{\text{c}}{{\text{m}}^{ - 3}}$ . The developed diagnostic is validated against a conventional single-point LIF method using photomultiplier tube-based detection, while significantly reducing the total measurement time by the number of spatial positions measured. The time-resolving capability of this diagnostic is tested by oscillating the plasma between two nominal operating modes with different density profiles and triggering the ICCD camera by the externally driven plasma oscillation. The measured 2D IVDF maps reveal several signatures of ion dynamics in this plasma source, including radially outflowing ions across the electric field and anomalous ion heating at the periphery, consistent with recent kinetic simulations and theoretical studies. A possible correlation between these ion kinetic features and rotating spoke structures is discussed.

E×B plasma↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Impurity radiation seeding of neoclassical tearing mode growth

Abstract The physics of neoclassical tearing modes (NTMs) is of great concern to tokamak plasma stability and performance, especially in the burning plasma regime. Whereas in many situations the seed events can be clearly identified, such as sawteeth and edge localized modes, the potential seeding mechanism of NTMs due to the resistive tearing instability driven by the impurity radiation cooling still needs more study. Recent NIMROD simulations have demonstrated that local impurity radiation cooling can drive the seed island growth and trigger the subsequent onset of NTM instability. The seed island is mainly driven by the local helical perturbation of the diamagnetic current induced by the perturbed pressure gradient as a result of the impurity radiative cooling on the rational surface. A heuristic closure for the neoclassical viscosity is adopted, and the seed island is further driven by the perturbed bootstrap current induced from the neoclassical electron viscous stress in the extended Ohm’s law. The growth rate of the NTM in simulations is found proportional to the electron neoclassical viscosity, and a theoretical neoclassical driving term is adopted to account for the nonlinear neoclassical island growth in the simulations.

Zeng, Shiyong (ORCID:000000017098818X)↗
DOE OSTI2022

Simultaneous stoquasticity

Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Furthermore, solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Selective Enhancement of Spectroscopic Features by Quantum Optimal Control

Tailored light can be used to steer atomic motions into selected quantum pathways. In optimal control theory (OCT), the target is usually expressed in terms of the molecular wave function, a quantity that is not directly observable in experiment. In this study, we present simulations using OCT that optimize the spectroscopic signal itself. By shaping the optical pump, the x-ray stimulated Raman signal, which occurs solely during the passage through conical intersections, is temporally controlled and amplified by up to 2 orders of magnitude. This enhancement can be crucial in order to bring small coherence-based signatures above the detectable threshold. Our approach is applicable to any signal that depends on the expectation value of a positive definite operator.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Learning Grain-Boundary Segregation: From First Principles to Polycrystals

The segregation of solute atoms at grain boundaries (GBs) can strongly impact the structural and functional properties of polycrystals. Yet, due to the limited availability of simulation tools to study polycrystals at the atomistic scale (i.e., interatomic potentials), there is a minimal understanding of the variation of solute segregation tendencies across the very complex space of GB microenvironments and the large range of alloys in which it can occur. Here, we develop an algorithmic framework that can directly learn the full spectrum of segregation energies for a metal solute atom in a metal polycrystal from ab initio methods, bypassing the need for alloy interatomic potentials. This framework offers a pathway to a comprehensive catalog of GB solute segregation with quantum accuracy, for the entire alloy space. As an initial demonstration in this pursuit, we build an extensive GB segregation database for aluminum-based alloys across the periodic table, including dozens of alloys for which there are substantially no prior data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Essential role of magnetic frustration in the phase diagrams of doped cobaltites

Doped perovskite cobaltites (e.g., La 1–x Sr x CoO 3 ) have been extensively studied for their spin-state physics, electronic inhomogeneity, and insulator-metal transitions. Ferromagnetically interacting spin-state polarons emerge at low x in the phase diagram of these compounds, eventually yielding long-range ferromagnetism. The onset of long-range ferromagnetism (x≈0.18) is substantially delayed relative to polaron percolation (x≈0.05), however, generating a troubling inconsistency. Here, in this study, Monte Carlo simulations of a disordered classical spin model are used to establish that previously ignored magnetic frustration is responsible for this effect, enabling faithful reproduction of the magnetic phase diagram.

36 MATERIALS SCIENCE↗