Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational physics”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hybrid gyrokinetic ion/fluid electron simulation of toroidal tearing modes

The effects of toroidicity and kinetic ions on the resistive tearing mode are systematically studied with the gyrokinetic particle-in-cell simulation code GEM and compared with analytic theory. A new field solver in toroidal geometry has been developed for the simulation of low-n (n = 1, 2) modes in tokamaks. It is found that the toroidal effect significantly reduces the growth rate of the tearing mode. The toroidal effect can also increase the radial width of the tearing mode and change the scaling between the radial mode width and resistivity due to the toroidal pressure term in the electron continuity equation. The kinetic effects of ions can decrease the growth rate of the tearing mode. The plasma flux-surface shaping is found to have significant effect on the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fast estimation of ion-pairing for screening electrolytes: A cluster can approximate a bulk liquid

The propensity for ion-pairing can often dictate the thermodynamic and kinetic properties of electrolyte solutions. Fast and accurate estimates of ion-pairing can thus be extremely valuable for supplementing design and screening efforts for novel electrolytes. We introduce an efficient cluster model to estimate the local ion-pair potential-of-mean-force between ionic solutes in electrolytes. The model incorporates an ion-pair and a few layers of explicit solvent in a gas-phase cluster and leverages an enhanced sampling approach to achieve high efficiency and accuracy. We employ harmonic restraints to prevent solvent escape from the cluster and restrict sampling of large inter-ion distances. We develop a cluster ion-pair sampling tool that implements our cluster model and demonstrate its potential utility for screening simple and poly-electrolyte systems.

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

Fluid structure interaction of a subaqueous pendulum: Analyzing the effect of wake correction via large eddy simulations

The dynamic behavior of a subaqueous cylindrical pendulum and corresponding flow dynamics are investigated. The objectives were twofold: (i) to examine whether the two-dimensional model equations sufficiently capture the three-dimensional dynamics and (ii) to investigate the emerging three-dimensional vortical flow structures. Large eddy simulations with two-way coupling fluid structure interaction were carried out using the immersed boundary method to simulate the motion of the pendulum and its interactions with the initially stagnant water. The resulting pendulum motion is compared against measured data obtained in a series of experimental tests to validate the simulation results and the model equations with and without wake corrections. An analysis of the flow vorticity revealed the development of a vortex ring during the first swing and the formation of tip vortices. The evolution of the vortex rings emerging from the motion of the subaqueous cylindrical pendulum was visualized using Q-criteria showing a reasonable agreement with vortical structures observed in the experiment using particle imaging velocimetry. The hydrodynamic moments acting on the simulated pendulum and the moments calculated from the model equations are analyzed. Using the insights from these numerical simulations, a modification of the wake correction is proposed to enhance the accuracy of the rate of decay and period. The transient effect of coherent flow on pendulum dynamics, especially the added mass effect, is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

An optoelectronic heterostructure for neuromorphic computing: CdS/V 3 O 5

Nonvolatile resistive switching is one of the key phenomena for emerging applications in optoelectronics and neuromorphic computing. In most of the cases, an electric field is applied to a two terminal dielectric material device and leads to the formation of a low resistance filament due to ion migration. However, the stochastic nature of the ion migration can be an impediment for the device robustness and controllability, with uncontrolled variations of high and low resistance states or threshold voltages. Here, we report an optically induced resistive switching based on a CdS/V 3 O 5 heterostructure which can overcome this issue. V 3 O 5 is known to have a second order insulator to metal transition around T c ≈ 415 K, with an electrically induced threshold switching at room temperature. Upon illumination, the direct transfer of the photoinduced carriers from the CdS into V 3 O 5 produces a nonvolatile resistive switching at room temperature. The initial high resistance can be recovered by reaching the high temperature metallic phase, i.e., temperatures above T c . Interestingly, this resistive switching becomes volatile around the T c . In conclusion, by locally manipulating the volatile and nonvolatile resistive switching using electric field and light, this system is a promising platform for hardware based neuromorphic computing implementations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Microgap breakdown with floating metal rod perturbations

Here, we report the characterization of microgap breakdown with perturbations from a metal rod floating between anode and cathode electrodes. The effects of the metal rod on the electric field distribution and the field enhancement factor are evaluated by numerical simulation and the conformal mapping method, and they indicate that the field emission regime is not reached. The breakdown voltages in the Townsend discharge regime are determined based on the voltage–current characteristics, which are obtained from two-dimensional fluid simulations. It is found that the breakdown characteristics can be significantly modulated by the floating metal rod, and the breakdown curve (breakdown voltage vs the net gap distance) is no longer U-shaped, which deviates from the conventional Paschen's law. The underlying physical mechanisms are related to the electric field enhancement, curved breakdown path, and nonuniform ion flux caused by the electric shielding effect. The results provide insights into breakdown characteristics in microscale discharges, which may promote conventional investigation of simplified clean gaps toward more complex conditions (e.g., with floating microparticles) in miniaturized plasma devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Investigation of boosted proton energies through proton radiography of target normal sheath acceleration fields in the multi-ps regime

Multi-kilojoule, multi-picosecond short-pulse lasers, such as the National Ignition Facility-Advanced Radiographic Capability laser and the OMEGA-Extended Performance laser, which have been constructed over the last two decades, enable exciting opportunities to produce high-brightness, high-energy laser-driven proton sources for applications in high-energy-density science like proton fast ignition for inertial fusion energy, particle radiography, and materials science studies. Results on these platforms have demonstrated enhanced accelerated proton energies and electron temperatures when compared to established scaling laws. Recent work has developed a new scaling for proton TNSA in the multi-ps regime. However, this new physics in the multi-ps regime motivates the need to understand the origin of the enhancement in proton energies. Toward this goal, here this work presents the first measurements of the TNSA accelerating sheath field in the multi-ps regime for pulse durations of 0.6, 5, and 10 ps. This measurement was achieved by using a separate TNSA proton source to radiograph the spatiotemporal profile of the accelerating sheath that is responsible for proton acceleration. The use of stacked radiochromic film detectors allows for a discrete time profile of the radiographs, thus enabling the measurement of the temporal and spatial evolution of the accelerating field. In performing this measurement, we extract quantities such as the sheath strength as a function of time and pulse duration, which shows that longer pulse durations sustain a stronger electric field for a longer duration when compared to sub-ps laser pulses, which may enable the observed boosted proton energies and proton conversion efficiencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

A nonlinear journey from structural phase transitions to quantum annealing

Motivated by an exact mapping between equilibrium properties of a one-dimensional chain of quantum Ising spins in a transverse field (the transverse field Ising (TFI) model) and a two-dimensional classical array of particles in double-well potentials (the “$\phi$ 4 model”) with weak inter-chain coupling, we explore connections between the driven variants of the two systems. Here, we argue that coupling between the fundamental topological solitary waves in the form of kinks between neighboring chains in the classical $\phi$ 4 system is the analog of the competing effect of the transverse field on spin flips in the quantum TFI model. As an example application, we mimic simplified measurement protocols in a closed quantum model system by studying the classical $\phi$ 4 model subjected to periodic perturbations. This reveals memory/loss of memory and coherence/decoherence regimes, whose quantum analogs are essential in annealing phenomena. In particular, we examine regimes where the topological excitations control the thermal equilibration following perturbations. This paves the way for further explorations of the analogy between lower-dimensional linear quantum and higher-dimensional classical nonlinear systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Results from a synthetic model of the ITER XRCS-Core diagnostic based on high-fidelity x-ray ray tracing

A high-fidelity synthetic diagnostic has been developed for the ITER core x-ray crystal spectrometer diagnostic based on x-ray ray tracing. This synthetic diagnostic has been used to model expected performance of the diagnostic, to aid in diagnostic design, and to develop engineering tolerances. The synthetic model is based on x-ray ray tracing using the recently developed xicsrt ray tracing code and includes a fully three-dimensional representation of the diagnostic based on the computer aided design. The modeled components are: plasma geometry and emission profiles, highly oriented pyrolytic graphite pre-reflectors, spherically bent crystals, and pixelated x-ray detectors. Plasma emission profiles have been calculated for Xe 44+ , Xe 47+ , and Xe 51+ , based on an ITER operational scenario available through the Integrated Modelling & Analysis Suite database, and modeled within the ray tracing code as a volumetric x-ray source; the shape of the plasma source is determined by equilibrium geometry and an appropriate wavelength distribution to match the expected ion temperature profile. All individual components of the x-ray optical system have been modeled with high-fidelity producing a synthetic detector image that is expected to closely match what will be seen in the final as-built system. Particular care is taken to maintain preservation of photon statistics throughout the ray tracing allowing for quantitative estimates of diagnostic performance.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Posterior comparison of model dynamics in several hybrid turbulence model forms

Hybrid turbulence models that can accurately reproduce unsteady three-dimensional flow physics across the entire range of grid scales and turbulence dynamics from Reynolds-averaged Navier–Stokes (RANS), through large-eddy simulation (LES), down to direct numerical simulations (DNS) are of increasing interest to the turbulence modeling community. However, despite decades of research and development, the basic tasks of eliminating poor-performing hybrid RANS-LES models and accelerating adoption of superior models through well-designed validation and verification have yet to occur. As a step in this direction, in this work we evaluate thirteen different hybrid RANS-LES models via systematic grid refinement of decaying homogeneous isotropic turbulence. We further derive a novel mathematical framework for assessing the energy partitioning dynamics of each Hybrid RANS-LES model, wherein model-to-model variations in energy partitioning can be interpreted as different feedback mechanisms operating on a low-dimensional nonlinear dynamical system. We found that model forms similar to the flow simulation methodology—also often termed very-large eddy simulation—are dynamically inconsistent with DNS at all resolutions. Additionally, we found a strong dynamical similarity in the feedback mechanisms of all models related to detached eddy simulation and partially averaged Navier–Stokes that is inherent to their general model forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Evaluating nonlocal heat transport in directly driven chromium spheres using x-ray spectroscopy

We report on experiments investigating heat transport in laser-generated plasmas using directly driven chromium spheres. The spheres are fielded at the OMEGA laser facility and are driven with laser intensities of 5×10 14 Wcm −2 . Plasma conditions in the corona and scattered light are measured experimentally and compared against predictions from two-dimensional (2D) radiation-hydrodynamic simulations using different heat transport models. Spectroscopic analysis of x-ray self-emission is used as an additional diagnostic. X-ray emission is integrated over a large region of the plasma, probing regions that are not observed by localized optical Thomson scattering. In particular, x-ray emission peaks near the plasma critical density, so emission from optically thin lines provides information on plasma conditions where nonlocal transport is most likely to be significant. Three common heat transport models are considered: local transport with flux limiters f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. Consistent with previous work, both the high-flux (f = 0.15) and SNB models show good agreement with experimentally measured plasma conditions in the corona despite overpredicting laser absorption, whereas the low-flux (f = 0.03) model fails to match any experimental data. Conditions inferred from x-ray self-emission line ratios support this conclusion during the period of laser peak power, although synthetic spectra for all models fail to match the experiment during the transient portions of the pulse. For these reasons, the low-flux model is again rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Carbon in GaN as a nonradiative recombination center

Trap-assisted nonradiative recombination has been shown to limit the efficiency of optoelectronic devices. While substitutional carbon (⁠C N ) has been suggested to be a nonradiative recombination center in GaN devices, a complete recombination cycle including the two charge-state transition levels has not been previously described. In this work, we investigate the trap-assisted recombination process due to C N in GaN, including multiphonon emission, radiative recombination, trap-assisted Auger–Meitner (TAAM) recombination as well as thermal emission of holes. Our study shows the key role of TAAM processes at the high carrier densities relevant for devices. We also reveal the carrier-density regimes where thermal emission and radiative recombination are expected to play an observable role. Our results highlight that carbon concentrations exceeding ~10 17 cm −3 can have a noticeable impact on device efficiency, not just in GaN active layers but also in InGaN and AlGaN. Furthermore, our comprehensive formalism not only offers detailed results for carbon but also provides a general framework for assessing the multiple processes that participate in trap-assisted recombination in semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

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

SHarmonic: A fast and accurate implementation of spherical harmonics for electronic-structure calculations

The authors present SHarmonic, a new implementation of the spherical harmonics targeted for electronic-structure calculations. Their approach is to use explicit formulas for the harmonics written in terms of normalized Cartesian coordinates. This approach results in a code that is as precise as other implementations while being at least one order of magnitude more computationally efficient. The library can run on graphics processing units as well, achieving an additional order of magnitude in execution speed. This new implementation is simple to use and is provided under an open-source license; it can be readily used by other codes to avoid the error-prone and cumbersome implementation of the spherical harmonics.

Mathematics and Computing↗
DOE OSTI2015

Thermodynamic approach to the stability of multi-phase systems. Application to the Y 2 O 3 –Fe system

Oxide-metal systems (OMSs) are important in many practical applications, and therefore, are under extensive studies using a wide range of techniques. The most accurate theoretical approaches are based on density functional theory (DFT), which are limited to ~10 2 atoms. Multi-scale approaches, e.g., DFT+Monte Carlo, are often used to model OMSs at the atomic level. These approaches can describe qualitatively the kinetics of some processes but not the overall stability of OMSs. In this paper, we propose a thermodynamic approach to study equilibrium in multiphase systems, which can be sequentially enhanced by considering different defects and microstructures. We estimate the thermodynamic equilibrium by minimization the free energy of the whole multiphase system using a limited set of defects and microstructural objects for which the properties are calculated by DFT. As an example, we consider Y 2 O 3 +bcc Fe with vacancies in both the Y 2 O 3 and bcc Fe phases, Y substitutions and O interstitials in Fe, Fe impurities and antisite defects in Y 2 O 3 . The output of these calculations is the thermal equilibrium concentration of all the defects for a particular temperature and composition. The results obtained confirmed the high temperature stability of yttria in iron. As a result, model development towards more accurate calculations is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Quantum eigenvalue estimation via time series analysis

We present an efficient method for estimating the eigenvalues of a Hamiltonian H from the expectation values of the evolution operator for various times. For a given quantum state ρ, our method outputs a list of eigenvalue estimates and approximate probabilities. Each probability depends on the support of ρ in those eigenstates of H associated with eigenvalues within an arbitrarily small range. The complexity of our method is polynomial in the inverse of a given precision parameter ε, which is the gap between eigenvalue estimates. Unlike the well-known quantum phase estimation algorithm that uses the quantum Fourier transform, our method does not require large ancillary systems, large sequences of controlled operations, or preserving coherence between experiments, and is therefore more attractive for near-term applications. The output of our method can be used to estimate spectral properties of H and other expectation values efficiently, within additive error proportional to ε.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗