Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time Evolution”

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

Analysis of silica fouling on nonwetting surfaces

Ground water sources used as coolant fluids in a variety of thermal systems such as heat exchangers and power plant condensers contain silica particles that accrete on heat transfer surfaces over time leading to reduction in thermal performance, a problem that is particularly exacerbated with temperature. Nonwetting superhydrophobic, lubricant-infused, and a new class of solid-infused surfaces introduced in this work are candidates for fouling mitigation, by virtue of their water repellency, but little is known about fouling of silica on the surfaces, especially under dynamic flow conditions and as a function of temperature. This article presents, for the first time, a systematic study of dynamic flow fouling of silica on nonwetting surfaces vis-à-vis conventional copper surface over a temperature range 20–50°C. The mechanism of silica aggregate formation and its adherence to the different surfaces is elucidated by scanning electron microscope (SEM) imaging. Sigmoidal growth model is used to describe the time evolution of fouling thermal resistance and an Arrhenius model is presented for the temperature-dependent increase in the asymptotic fouling resistance on nonwetting and conventional surfaces alike. Lubricant-infused and solid-infused surfaces are shown to reduce fouling thermal resistance by up to 25% and 13%, respectively, compared to a conventional surface, whereas superhydrophobic surfaces lose their non-wettability under flow conditions, leading to an adverse increase in the fouling resistance by up to 13%. Considering the possible lubricant depletion in lubricant-infused surfaces over prolonged exposure to a flowing fluid, solid-infused surfaces present a robust alternative.

42 ENGINEERING↗
DOE OSTI2025

Quantum dynamics of the temporary capture of light atoms by superfluid helium nanodroplets at very low collision energies (≈1–13 meV): the case of the hydrogen atom and its isotopes

The capture dynamics of a H atom and isotopic variants [D, T and Q (hypothetical isotope of mass equal to four times the mass of H)] by a superfluid helium nanodroplet (HeND) has been investigated theoretically. The HeND (T = 0.37 K) is ( 4 He) N=400 and a mean field quantum hybrid approach [TDDFT (helium) + quantum wave packet (H, D, T or Q)] at zero angular momentum, is used to explore a rather wide range of very low initial kinetic energies (E k,0 ≈ 10–150 K). The analysis of the capture mechanism shows the existence of a dynamical barrier and a dynamical minimum that play key roles to understand the time evolution of the capture, especially the former property. In general, the H atom shows a different behavior from the other isotopes, with the behavior of T and Q being very similar to each other and the D atom behaving inbetween H and T. Besides, it is worth noting that, in principle, at the very low initial kinetic energies considered only “short” and “long” lived atom⋯HeND collision complexes are formed, i.e., in the atom-helium nanodroplet collision only the temporary capture of the atom takes place. The different behaviors observed have been interpreted considering the faster motion of the H atom when colliding with ( 4 He) N=400 and the more quantum character of the H behavior both due to its significantly lower mass. As far as we know, this is the first quantum dynamics study carried out on the collision of light atoms with HeNDs at very low energies.

Sternberg, Michael [Argonne National Laboratory (A↗
DOE OSTI2020

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

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

Determining population densities in bimodal micellar solutions using contrast-variation small angle neutron scattering

Self-assembly of amphiphilic polymers in water is of fundamental and practical importance. Significant amounts of free unimers and associated micellar aggregates often coexist over a wide range of phase regions. The thermodynamic and kinetic properties of the microphase separation are closely related to the relative population density of unimers and micelles. Although the scattering technique has been employed to identify the structure of micellar aggregates as well as their time-evolution, the determination of the population ratio of micelles to unimers remains a challenging problem due to their difference in scattering power. In this work, using small-angle neutron scattering (SANS), we present a comprehensive structural study of amphiphilic n-dodecyl-PNIPAm polymers, which shows a bimodal size distribution in water. By adjusting the deuterium/hydrogen ratio of water, the intra-micellar polymer and water distributions are obtained from the SANS spectra. The micellar size and number density are further determined, and the population densities of micelles and unimers are calculated to quantitatively address the degree of micellization at different temperatures. Our method can be used to provide an in-depth insight into the solution properties of microphase separation, which are present in many amphiphilic systems.

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

Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas

A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.

Solmaz, Evrim (ORCID:0000000211830881)↗
DOE OSTI2021

Constructing tensor network influence functionals for general quantum dynamics

Here, we describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of its bond dimension and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation.

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

A multi-chord, two-color interferometer using Hilbert transform phase detection for measuring electron density in spheromak plasmas

A new, four-chord, CO 2 /He–Ne heterodyne interferometer has been designed and built for measuring line-averaged plasma density in the HIT-SI3 and subsequent HIT-SIU sustained spheromak devices. The two-color system successfully eliminates vibration-induced errors caused by mirrors that are secured to the vacuum chamber and is able to resolve electron densities ne in the full operating range of 10 18 –10 20 m –3 in both experiments with an integrated error of 4.7 × 10 17 m –2 . Here, data are presented from high toroidal current plasma discharges, showing the time evolution of electron densities n e and j Φ /n e along multiple chords.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Resonant three-wave interactions with linear source/sink terms: Exploring the parameter space

Here, it is shown that the resonance three wave interactions in the presence of linear source/sink terms exhibit a wide variety of time-evolution dynamics, such as exponential growth, limit cycles, and strange attractors. The particular features of the dynamics are determined primarily by the values of linear growth rates, forming a 2D parameter space. An exploration of this parameter space is the main focus of this present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quasilinear theory of general electromagnetic fluctuations including discrete particle effects for magnetized plasmas: General analysis

The general quasilinear Fokker–Planck kinetic equation for the gyrophase-averaged plasma particle distribution functions in magnetized plasmas is derived, making no restrictions on the energy of the particles and on the frequency of the electromagnetic fluctuations and avoiding the often made Coulomb approximation of the electromagnetic interactions. Here, the inclusion of discrete particle effects breaks the dichotomy of nonlinear kinetic plasma theory divided into the test particle and the test fluctuation approximation because it provides expression of both the non-collective and collective electromagnetic fluctuation spectra in terms of the plasma particle distribution functions. Within the validity of the quasilinear approach, the resulting full quasilinear transport equation can be regarded as a determining nonlinear equation for the time evolution of the plasma particle distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Flexocaloric effect in superelastic materials

We present a combined theoretical–experimental study of flexocaloric effects in superelastic materials exhibiting structural transitions. We study a Ginzburg–Landau model combined with a vibrational model for a beam near a ferroelastic transition loaded with a three-point bending setup. We also perform experiments on a Cu–Al–Ni single crystal undergoing a martensitic transition. We measure bent beam profiles, vertical force vs vertical deflection during a slow isothermal process, time evolution of the bending and unbending amplitudes, and the evolution of temperature profiles. We also compute the evolution of heat source and heat sink profiles. Finally, we study the location of acoustic emission events during the bending/unbending experiment. Our observations are consistent with the model predictions and allow us to identify the main physical parameters relevant for flexocaloric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗
DOE OSTI2023

Disruption halo current rotation scaling on Alcator C-Mod and HBT-EP

Asymmetric halo currents (HCs) can exert large net forces on the vacuum vessel and other components during disruptions on tokamaks. The displacements caused by these forces can then be amplified if these asymmetric forces rotate at frequencies resonant with the vessel. This paper reports on the investigation of a recently proposed scaling law for the disruption HC rotation frequency [Saperstein et al., “Halo current rotation scaling in post-disruption plasmas,” Nucl. Fusion 62, 026044 (2022)] that combines measurements on Alcator C-Mod with those on HBT-EP. We find that a new non-circular version of the scaling law [ ⟨ f rot ⟩ m / ⟨ m ⟩ ∝ 1 B T ( S / π )] takes into consideration the dependence of frot on the poloidal structure of the MHD instability (m) driving the asymmetry and describes the disruption-averaged rotation frequency on C-Mod. Disruption rotation is also found to be insensitive to the vertical position and impurity content of the plasma at the onset of the disruption. However, a stagnation in the time evolution of frot is occasionally observed. In conclusion, observations are consistent with the dominance of poloidal rotation during the disruption, which is motivated by the poloidal drift nature of the scaling law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

NiGa 2 O 4 interfacial layers in NiO/Ga 2 O 3 heterojunction diodes at high temperature

NiO/Ga 2 O 3 heterojunction diodes have attracted attention for high-power applications, but their high temperature performance and reliability remain underexplored. Here, we report the time evolution of the electrical properties in the widely studied p-NiO/n-Ga 2 O 3 heterojunction diodes and formation of NiGa 2 O 4 interfacial layers at high temperatures. Results of our thermal cycling experiment show an initial leakage current increase which stabilizes after sustained thermal load, due to reactions at the NiO–Ga 2 O 3 interface. High-resolution TEM microstructure analysis of the devices after thermal cycling indicates that the NiO–Ga 2 O 3 interface forms a ternary compound at high temperatures, and thermodynamic calculations suggest the formation of the spinel NiGa 2 O 4 layer between NiO and Ga 2 O 3 . First-principles defect calculations find that NiGa 2 O 4 shows low p-type intrinsic doping and hence can serve to limit electric field crowding at the interface. Vertical NiO/ Ga 2 O 3 diodes with intentionally grown 5 nm thin spinel-type NiGa 2 O 4 interfacial layers show an excellent device ON/OFF ratio of >10 10 (± 3 V), V ON of ~1.9 V, and increased breakdown voltage of ~1.2 kV for an initial unoptimized 300 lm diameter device. These p–n heterojunction diodes are promising for high-voltage, high temperature applications.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum–classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

Chemistry↗
DOE OSTI2024

A new method for solving the linearized 1D Vlasov–Poisson system yielding a new class of solutions

We describe a new method for solving the linearized 1D Vlasov–Poisson system by using properties of Cauchy-type integrals. Our method remedies critical flaws of the two standard methods, reveals a previously unrecognized Gaussian-in-time-like decay, and can also account for an externally applied electric field. The Landau approximation involves deforming the Bromwich contour around the poles closest to the real axis due to the analytically continued dielectric function, finding the long-time behavior for a stable system: Landau damping. Jackson's generalization encircles all poles while sending the contour to infinity, assuming its contribution vanishes, which is not true in general. This gives incorrect solutions for physically reasonable configurations and can exhibit pathological behavior, of which we show examples. The van Kampen method expresses the solution for a stable equilibrium as a continuous superposition of waves, resulting in an opaque integral. Case's generalization includes unstable systems and predicts a decaying discrete mode for each growing discrete mode, an apparent contradiction to both the Jackson solution and ours. We show, without imposing additional constraints, that the decaying modes are never present in the time evolution due to an exact cancellation with part of the continuum. Our solution is free of integral expressions, is obtained using algebra and Laurent series expansions, does not rely on analytic continuations, and results in a correct asymptotically convergent form in the case of infinite sums. The analysis used can be readily applied in higher-dimensional, electromagnetic systems and also provides a new technique for evaluating certain inverse Laplace transforms.

Physics↗
DOE OSTI2025

Planar shock-induced bubble collapse and jetting in water captured via X-ray phase contrast imaging

Shock wave–bubble interactions in water manifest rich dynamics driven by a combination of strong pressure and density mismatches. They have a wide variety of applications, including the injection of pharmaceuticals, and through scaling, enable the exploration of various aspects of high-energy-density systems such as inertial confinement fusion. In this work, the interaction between a micrometric nitrogen bubble and a planar shock wave, characterized by a Mach number of M = 1.24 and a peak pressure of p max = 0.57 , is experimentally recorded using ultra-high-speed x-ray phase contrast imaging. Highly resolved radiographs provide access to all phase discontinuities along the beam path, offering quantities such as the time-varying bubble size, the speed of a jet produced during the bubble collapse, and the time evolution of the shock wave front, which are critical benchmark data for numerical scheme validation. This study addresses the lack of well-characterized, repeatable, and high spatiotemporal resolution experiments at negative Atwood numbers by providing shock–bubble visualization and corresponding numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗