SEARCH · Engineering Papers

Results for “gradients”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Phase-space dynamics of Alfvén mode chirping

In tokamak discharges, toroidal Alfvén eigenmodes often experience complex semi-periodic frequency modulation known as chirping. These events modify the local high energy particle distribution and are expected to occur in many future fusion devices, which include energetic beams or fusion products. In this paper, we present a study of simulations of mode chirping made in order to better understand its phase-space properties in a realistic tokamak configuration. We find a mechanism that permits rapid repeated chirping with strong amplitude variation in each chirp. Each chirp is then associated with an amplitude crash to low magnitude and local manipulation of the density gradients through a shift of mode phase through π. The chirping produces high density clumps, which propagate down the fast ion density gradient and low density holes that propagate up the density gradient away from the resonance. This flow of particles across the resonance provides an energy source and local gradients for repeated chirping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Closed-form theoretical model of the secondary drop size in partial coalescence—Capturing pertinent timescales and viscous forces

One of the most important outcomes of partial drop coalescence is the ratio of the secondary drop radius to the primary drop radius, known as the drop ratio, ri. This ratio is thought to be approximately constant and independent of physical parameters of the fluids involved. However, here we reveal that ambient fluid viscosity can alter the size of the secondary drop and the drop ratio consequently. Using scaling analysis, we derive a model that predicts the behavior of the drop ratio as a function of the Ohnesorge number, a dimensionless ratio of viscous to inertial forces. In addition, we present our experimental results of coalescing drops on a planar interface under the influence of surface tension gradients. A high-speed digital camera is used to observe the evolution of drops as they coalesce with a bulk liquid. We show that this process is influenced by the surface tension gradient between the drop and the bulk liquid. The ratio of the secondary drop to the primary drop in partial coalescence is smaller than the reported values for coalescence without a surface tension gradient. The analytical model derived through this study is based on a new modified Ohnesorge number that includes surface tension gradients. Our analytical model is compared against other models and the results illustrate good agreement with our experimental findings and experimental data in the literature.

42 ENGINEERING↗

DOE OSTI2020

Pair plasma instability in homogeneous magnetic guide fields

Pair plasmas, collections of both matter and antimatter particles of equal mass, represent a paradigm for the study of basic plasma science, and many open questions exist regarding these unique systems. They are found in many astrophysical settings, such as gamma-ray bursts, and have recently also been produced in carefully designed laboratory experiments. A central research topic in plasma physics is instability; however, unlike their more common ion–electron siblings, pair plasmas are generally thought to be stable to cross field pressure gradients in homogeneous magnetic fields. It is shown here by means of kinetic full-f simulations that, when a pressure gradient is first established, the Gradient-driven Drift Coupling mode is destabilized and becomes turbulent. Force balance is eventually achieved by a combination of flattened pressure profiles due to turbulent transport and establishment of a magnetic field gradient, saturating the growth. During the unstable phase, key physics can be captured by a δf gyrokinetic description, where it is shown analytically and numerically that parallel particle motion results in a coupling of all electromagnetic field components. A fluid model derived therefrom accurately predicts linear eigenmodes and is used to resolve global profile effects. For laser-based electron–positron plasma experiments, prompt instability is predicted with growth times much shorter than plasma lifetimes. Similarly, growth rates are calculated for the planned APEX experiment as well as gamma-ray burst scenarios, suggesting that the instability may contribute to the early evolution of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Weibel-like instability in magnetohydrodynamics

In magnetohydrodynamics (MHD), a density perturbation perpendicular to an electron temperature gradient generates a magnetic field around itself that acts to increase the perturbation, which can lead to instability. An MHD dispersion relation is obtained for perturbations perpendicular to a fixed electron temperature gradient with an initial in-plane magnetic field, including resistivity, viscosity, and the electrothermal coefficient. Instability occurs for sufficiently small electron temperature-gradient scale lengths determined by the ion collisionless skin depth. Both viscosity and resistivity are required to prevent growth at arbitrarily small spatial scales and to give a physical result for the fastest growing mode. The perpendicular electrothermal coefficient is only significant for a narrow range of low electron Hall parameters, causing a modest reduction in magnetic field growth and modifying the criteria for instability in the presence of viscosity. If the definition of the Weibel instabilityis extended to include all instabilities due to anisotropy in the electron velocity distribution, then this is a Weibel-like instability because an electron temperature gradient implies an anisotropic electron velocity distribution. Here, the implications for the formation of filaments in laser-produced plasmas and for the verification of MHD codes are considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Using the Metropolis algorithm to explore the loss surface of a recurrent neural network

In the limit of small trial moves the Metropolis Monte Carlo algorithm is equivalent to gradient descent on the energy function in the presence of Gaussian white noise. This observation was originally used to demonstrate a correspondence between Metropolis Monte Carlo moves of model molecules and overdamped Langevin dynamics, but it also applies in the context of training a neural network: making small random changes to the weights of a neural network, accepted with the Metropolis probability, with the loss function playing the role of energy, has the same effect as training by explicit gradient descent in the presence of Gaussian white noise. We explore this correspondence in the context of a simple recurrent neural network. We also explore regimes in which this correspondence breaks down, where the gradient of the loss function becomes very large or small. In these regimes the Metropolis algorithm can still effect training, and so can be used as a probe of the loss function of a neural network in regimes in which gradient descent struggles. We also show that training can be accelerated by making purposely-designed Monte Carlo trial moves of neural-network weights.

Casert, Corneel↗

DOE OSTI2026

Non-classical ion-indiffusion processes in quasi-one-dimensional ionic conductor KTiOPO4

Developments in coercive field engineering and waveguide writing via ion-exchange in the ferroelectric ionic-conductor KTiOPO4 (KTP) have led to breakthroughs within the field of nonlinear optics and quantum photonics. However, the microscopic diffusion dynamics of ion-exchange in KTP are poorly understood and show signs of non-classical diffusion dynamics. The precise control of ion-exchange gradients is critical for developing KTP into a mature, integrated nonlinear optical material but requires a level of control that, in turn, demands a deeper understanding of the ion-exchange process. We employ direct ToF-SIMS mapping to untangle the complex diffusion dynamics of the commonly used ion-exchange dopant ions: Rb, Ba, and K. We map periodic diffusion gradients of two exchanges, one dominated by Rb-ions and the other by Ba-ions. For the first time, we show the interdependent nature of the Ba- and Rb ion in the ion-exchange diffusion gradient. The low incorporation of the Ba ion is mapped and quantified, and dopant–crystal interactions are identified. This work highlights the complex diffusion dynamics of dopants in the KTP structure and lays the foundation for deterministic control of diffusion gradients for the next generation of integrated nonlinear and quantum photonic devices.

Barrett, Laura [ORNL]↗

DOE OSTI2026

Neoclassical transport and profile prediction in transport barriers

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Trinczek, Silvia [Princeton Plasma Physics Laborat↗

DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗

DOE OSTI2023

Spontaneous skyrmion conformal lattice and transverse motion during dc and ac compression

Abstract We use atomistic-based simulations to investigate the behavior of ferromagnetic skyrmions being continuously compressed against a rigid wall under dc and ac drives. The compressed skyrmions can be annihilated close to the wall and form a conformal crystal with both a size and a density gradient, making it distinct from conformal crystals observed previously for superconducting vortices and colloidal particles. For both dc and ac driving, the skyrmions can move transverse to the compression direction due to a combination of density and size gradients. Forces in the compression direction are converted by the Magnus force into transverse motion. Under ac driving, the amount of skyrmion annihilation is reduced and we find a skyrmion Magnus ratchet pump. We also observe shear banding in which skyrmions near the wall move up to twice as fast as skyrmions further from the wall. When we vary the magnitude of the applied drive, we find a critical current above which the skyrmions are completely annihilated during a time scale that depends on the magnitude of the drive. By varying the magnetic parameters, we find that the transverse motion is strongly dependent on the skyrmion size. Smaller skyrmions are more rigid, which interferes with the size gradient and destroys the transverse motion. We also confirm the role of the size gradient by comparing our atomistic simulations with a particle-based model, where we find that the transverse motion is only transient. Our results are relevant for applications where skyrmions encounter repulsive magnetic walls, domain walls, or interfaces.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Heat pulse propagation and anomalous electron heat transport measurements on the optimized stellarator W7-X

The optimized stellarator Wendelstein 7-X (W7-X) is designed to have an approximately quasi-isodynamic magnetic configuration with reduced neoclassical transport in comparison to a classical stellarator, and turbulent transport is expected to be a significant source of anomalous heat transport across the plasma minor radius. The ion temperature gradient driven mode and the trapped electron mode (TEM) are thought to be responsible for the ion-scale turbulence in W7-X plasmas with volume averaged pressure below 1%. In this work, the electron temperature gradient driven turbulence is shown to be a good candidate for the explanation of the observed electron heat flux, in the inner plasma region where the density gradient is weak (in the outer region, a relatively stronger density gradient would drive additional TEM turbulence). Further, the experimental electron heat transport measured during electron cyclotron resonant heating power and plasma density scans is compared to neoclassical predictions, and the stiffness in the electron heat transport measured during transient transport experiments is presented in three common magnetic configurations of W7-X. In low-(β) plasma discharges, the stiffness in the electron heat flux, quantified by the ratio of the heat pulse to power balance diffusivity, ${\chi }_{\mathrm{e}}^{\mathrm{H}\mathrm{P}}/{\chi }_{\mathrm{e}}^{\mathrm{P}\mathrm{B}}$, is measured to be less than 2, and trend downwards with increasing collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

On the role of mode resonances in regulating zonal-flow-moderated plasma microturbulence

Abstract The onset of turbulent heat transport at a higher temperature gradient than the critical gradient of linear instability, known as the Dimits shift, is a recurring feature of nonlinear simulations for magnetically confined fusion plasmas. Resonance in the nonlinear coupling between the modes that dominate energy transfer can lead to suppression of turbulence and transport above the linear critical gradient. As an expression of this resonance, gyrokinetic simulations show a quasi-coherent interaction between streamers and sidebands coupled through the zonal flow within the Dimits regime. This mechanism is further confirmed by use of artificial complex frequencies which break the resonance. By incorporating corresponding saturation physics, the standard quasilinear model for rapid head flux prediction is improved, which can now predict reduced heat flux in the Dimits regime. In particular, the triplet correlation time, the lifetime of the nonlinear interaction, is shown to be well-approximated by combinations of linear eigenvalues, and yields good representations of the heat flux variation both in and above the Dimits regime. Thus, a reduced but predictive model for transport near the critical gradient of zonal-flow saturated turbulence now exists.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Importance of $\delta B_{\|}$ on ETG stability, turbulence, and transport in NSTX

This study employs electron-scale gyrokinetic simulations to investigate the electron temperature gradient (ETG) driven instabilities, turbulence, and transport in the pedestal region of the National Spherical Torus Experiment, comparing non-lithiated (narrow pedestal) and lithiated (wide pedestal) scenarios. Our findings reveal that, in the non-lithiated case, a branch of strongly unstable ETG modes exhibiting finite parallel magnetic field fluctuations ($\delta B_{\parallel} \neq 0$) emerges at the pedestal top and upper density pedestal region. This branch is uncovered only when $\delta B_{\parallel}$ is retained in the simulations and is associated with substantial electrostatic electron heat flux. This region of strong ETG transport corresponds to the only region in the plasma where the pressure gradient is far below the critical gradient for kinetic ballooning modes. We investigated the origin of this finite $\delta B_{\parallel}$ ETG branch by analyzing the gyrokinetic field equations. Nonlinear saturation is also analyzed and contrasted for simulations with and without $\delta B_{\parallel}$. In contrast with the nonlithiated case, ETG modes in the lithiated case produce substantial transport in the steep gradient region, but are negligible at the pedestal top.

ETG↗

DOE OSTI2024

Robust wind farm layout optimization

Wake interactions in wind farms cause losses in annual energy production (AEP) on the order of 10%. Wind farm designers optimize the layout of the farm to mitigate wake losses, especially in the dominant site-specific wind directions. As wind turbines and wind farms grow in scale, optimization becomes more complex. Offshore wind farms regularly comprise more than 100 wind turbines and are characterized by complex boundaries due to shipping lanes, neighboring wind farms, and other constraints. Layout optimization methods are broadly split between gradient-based and gradient-free approaches. Gradient-based approaches can converge quickly and perform well for smaller, academic problems but are often sensitive to initial conditions and tuning parameters and require expert knowledge to use. On the other hand, gradient-free approaches can be more robust to problem complexities. We present a robust layout optimization approach based on a random search algorithm. The algorithm is intended for those who are not optimization experts and has few tuning parameters that need specification to achieve satisfactory results. Unlike off-the-shelf methods, which use generally available, non-domain-specific optimization routines that accept as inputs an optimization function and constraint definitions, this approach takes advantage of the relative computational costs of the different evaluations by evaluating cheaper computations first (boundary and minimum distance constraints) and running expensive AEP evaluations only if all other checks pass. Moreover, an outer genetic algorithm allows multiple solutions to evolve in parallel, enabling rapid solution development on high-performance computers. We discuss the relative ease of selecting necessary tuning parameters and demonstrate the efficacy of the genetic random search on a complex layout problem consisting of placing 70 turbines in a nonconvex and unconnected boundary region.

17 WIND ENERGY↗

DOE OSTI2024

The science case for an intermediate energy advanced and novel accelerator linear collider facility

It is widely accepted that the next lepton collider beyond a Higgs factory would require center-of-mass energy of the order of up to 15 TeV. Since, given reasonable space and cost restrictions, conventional accelerator technology reaches its limits near this energy, high-gradient advanced acceleration concepts are attractive. Advanced and novel accelerators (ANAs) are leading candidates due to their ability to produce acceleration gradients on the order of 1–100 GV/m, leading to compact acceleration facilities. However, intermediate energy facilities (IEF) are required to test the critical technology elements on the way towards multi-TeV-class colliders. Here a science case for a 20–100 GeV center-of-mass energy ANA-based lepton collider that can be a candidate for an intermediate energy facility is presented. The IEF can provide numerous opportunities for high energy physics studies including precision Quantum Chromodynamics and Beyond the Standard Model physics measurements, investigation of charged particle interactions with extreme electromagnetic fields, and exploring muon and proton beam acceleration. Possible applications of this collider include the studies of γγ and electron beam-fixed target/beamdump collider designs. Thus, the goal of the proposed IEF is to both carry out particle physics measurements in the 20-100 GeV ranges as well as to serve as an ANA demonstrator facility.

43 PARTICLE ACCELERATORS↗