Engineering Papers⌕ Search

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.

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

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Ground-state properties of the Hubbard model in one and two dimensions from the Gutzwiller conjugate gradient minimization theory

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Benchmarks and results of the two-band Hubbard model from the Gutzwiller conjugate gradient minimization theory

Ground-state properties, such as energies and double occupancies, of a one-dimensional two-band Hubbard model are calculated using a first-principles Gutzwiller conjugate gradient minimization theory. The favorable agreement with the results from the density matrix renormalization group theory demonstrates the accuracy of our method. A rotationally invariant approach is further incorporated into the method to greatly reduce the computational complexity with a speedup of approximately 50 times. Moreover, we investigate the Mott transition between a metal and a Mott insulator by evaluating the charge gap. In conclusion, with greatly reduced computational effort, our method reproduces the phase diagram in reasonable agreement with the density matrix renormalization group theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

𝐹 𝐾 ⁡/𝐹 𝜋 from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of 𝐹 𝐾 ⁡/𝐹 𝜋 using Möbius domain-wall fermions computed on gradient-flowed 𝑁 𝑓 =2 +1 +1 highly improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from 130 ≲ 𝑚 𝜋 ≲ 400 MeV , four lattice spacings of 𝑎 ∼ 0.15, 0.12, 0.09 and 0.06 fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the 𝑎 ∼ 0.06 fm ensemble is helpful, but not necessary to achieve a subpercent determination of 𝐹 𝐾 ⁡/𝐹 𝜋 . We also include an estimate of the strong isospin breaking corrections and arrive at a final result of 𝐹 $\hat{K}$ + ⁡ /𝐹 $\hat{𝜋}$ + = 1.1942⁢(45) with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of |𝑉 𝑢⁢𝑠 |/|𝑉 𝑢⁢𝑑 | = 0.2311⁢(10).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Laser-direct-drive fusion target design with a high- Z gradient-density pusher shell

Laser-direct-drive fusion target designs with solid deuterium-tritium (DT) fuel, a high-Z gradient-density pusher shell (GDPS), and a Au-coated foam layer have been investigated through both 1D and 2D radiation-hydrodynamic simulations. Compared with conventional low-Z ablators and DT-push-on-DT targets, these GDPS targets possess certain advantages of being instability-resistant implosions that can be high adiabat (α ≥ 8) and low hot-spot and pusher-shell convergence (CR hs ≈ 22 and CR PS ≈ 17), and have a low implosion velocity (v imp < 3 × 10 7 cm/s). Using symmetric drive with laser energies of 1.9 to 2.5 MJ, 1D LILAC simulations of these GDPS implosions can result in neutron yields corresponding to >~50–MJ energy, even with reduced laser absorption due to the cross-beam energy transfer (CBET) effect. Two-dimensional DRACO simulations show that these GDPS targets can still ignite and deliver neutron yields from 4 to ~10 MJ even if CBET is present, while traditional DT-push-on-DT targets normally fail due to the CBET-induced reduction of ablation pressure. If CBET is mitigated, these GDPS targets are expected to produce neutron yields of >20 MJ at a driven laser energy of ~2 MJ. Here are the key factors behind the robust ignition and moderate energy gain of such GDPS implosions are as follows: (1) The high initial density of the high-Z pusher shell can be placed at a very high adiabat while the DT fuel is maintained at a relatively low-entropy state; therefore, such implosions can still provide enough compression ρR >1 g/cm 2 for sufficient confinement; (2) the high-Z layer significantly reduces heat-conduction loss from the hot spot since thermal conductivity scales as ~1/Z; and (3) possible radiation trapping may offer an additional advantage for reducing energy loss from such high-Z targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Measurement of pressure gradients near the interface in the viscous fingering instability

The viscous fingering instability, which forms when a less-viscous fluid invades a more-viscous one within a confined geometry, is an iconic system for studying pattern formation. For both miscible and immiscible fluid pairs the growth dynamics change after the initial instability onset and the global structures, typical of late-time growth, are governed by the viscosity ratio. Here we introduce an experimental technique to measure flow throughout the inner and outer fluids. This probes the existence of a new length scale associated with the local pressure gradients around the interface and allows us to compare our results to the predictions of a previously proposed model for late-time finger growth. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

Robust Avoidance of Edge-Localized Modes alongside Gradient Formation in the Negative Triangularity Tokamak Edge

In a series of high performance diverted discharges on DIII-D, we demonstrate that strong negative triangularity (NT) shaping robustly suppresses all edge-localized mode (ELM) activity over a wide range of plasma conditions: ⟨n⟩ = 0.1 − 1.5 × 1020 m−3, Paux = 0 − 15 MW and |Bt| = 1 − 2.2 T, corresponding to Ploss/PLH08 ∼ 8. The full dataset is consistent with the theoretical prediction that magnetic shear in the NT edge inhibits access to ELMing H-mode regimes; all experimental pressure profiles are found to be at or below the infinite-n ballooning stability limit. Our present dataset also features edge pressure gradients in strong NT that are closer to an H–mode than a typical L-mode plasma, supporting the consideration of NT for reactor design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient quantum circuits based on the quantum natural gradient

Efficient preparation of arbitrary entangled quantum states is crucial for quantum computation. This is particularly important for noisy intermediate-scale quantum simulators relying on variational hybrid quantum-classical algorithms. To that end, we propose symmetry-conserving modified quantum approximate optimization algorithm (SCom-QAOA) circuits. The depths of these circuits depend not only on the desired fidelity to the target state but also on the amount of entanglement the state contains. The parameters of the SCom-QAOA circuits are optimized using the quantum natural gradient method based on the Fubini-Study metric. The SCom-QAOA circuit transforms an unentangled state into a ground state of a gapped one-dimensional Hamiltonian with a circuit depth that depends not on the system size but rather on the finite correlation length. In contrast, the circuit depth grows proportionally to the system size for preparing low-lying states of critical one-dimensional systems. Even in the latter case, SCom-QAOA circuits with depth less than the system size were sufficient to generate states with fidelity in excess of 99%, which is relevant for near-term applications. The proposed scheme enlarges the set of the initial states accessible for variational quantum algorithms and widens the scope of investigation of nonequilibrium phenomena in quantum simulators. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Energy Level Gradients from Surface to Bulk in Hybrid Metal-Halide Perovskite Thin Films

Variations in local strain, defect densities, and composition of hybrid metal-halide perovskites have been reported to create heterogeneous energy landscapes in thin films, which impact charge-carrier diffusion and recombination dynamics. Here, we employ one- and two-photon transient absorption spectroscopy to selectively probe the dynamics of charge carriers from surface and bulk regions of methylammonium lead bromide thin films. Differences in the transient absorption spectra indicate that an energy gradient of approximately 100 meV is formed between the higher band-gap surface and lower band-gap bulk regions. Thus, during their lifetime, photoexcited carriers move away from the surface to recombine in the bulk, where our experiments detect long-lived charge populations despite the significant band splitting that has conventionally been assumed to inhibit efficient radiative recombination. Supported by first-principles calculations, we demonstrate that bright emission can still arise from the bulk with states that occupy a wide range of momenta in the vicinity of the band extrema, which show strong dipole transitions. Our results report that photoexcitations in the hybrid perovskites avoid defect-rich surface regions, and that particularly strong emission is generated from accumulated excitation populations in the bulk. Published by the American Physical Society 2024

Bourelle, Sean A.↗
DOE OSTI2026

High-precision scale setting with the Ω-baryon mass and gradient flow

The gradient-flow scale 𝑤 0 in lattice QCD is determined using the mass of the Ω − baryon to set the physical scale. Nine ensembles using the highly improved staggered quark (HISQ) action with lattice spacings of 0.15 fm down to 0.04 fm are used, seven of which have nearly physical light-quark masses. Electromagnetic corrections to the Ω − mass are defined in order to compute a pure-QCD Ω mass. The final result is 𝑤 0 = 0.17187⁢(68) fm, corresponding to a relative uncertainty of 0.40% and a central value in good agreement with previous calculations in the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Gradient Flow for Parton Distribution Functions: First Application to the Pion

Parton distribution functions (PDFs) are central to precision QCD phenomenology. Their Mellin moments can be computed on the lattice, but direct determinations using local operators, besides ⟨𝑥⟩, face severe challenges from reduced hypercubic symmetry, limiting results to the lowest moments. A recently proposed method resolves these issues using gradient flow. We demonstrate the efficacy of this method by computing ratios of flavor nonsinglet pion PDF moments up to ⟨𝑥 5 ⟩ on four lattice spacings at 𝑚 𝜋 ≃ 411 MeV. The moments and reconstructed PDF agree quantitatively with recent phenomenological extractions.

Francis, Anthony [National Yang Ming Chiao Tung Un↗
DOE OSTI2022

Optimal Coordination of Distributed Energy Resources Using Deep Deterministic Policy Gradient

Recent studies showed that reinforcement learning (RL) is a promising approach for coordination and control of distributed energy resources (DER) under uncertainties. Many existing RL approaches, including Q-learning and approximate dynamic programming, are based on lookup table methods, which become inefficient when the problem size is large and infeasible when continuous states and actions are involved. In addition, when modeling battery energy storage system (BESS), the loss of life is not reasonably considered into the decision-making process. This paper proposes an innovative deep RL method for DER coordination considering BESS degradation. The proposed deep RL is designed based on an adaptive actor-critic architecture and employs an off-policy deterministic policy gradient method for determining the dispatch operation that minimizes the operation cost and BESS life loss. Case studies were performed to validate the proposed method and demonstrate the effects of incorporating degradation models into control design.

Das, Avijit↗
DOE OSTI2022

Data-Driven PMU Noise Emulation Framework using Gradient-Penalty-Based Wasserstein GAN

Availability of phasor measurement unit (PMUs) data has led to research on data-driven algorithms for event monitoring, control and ensuring stability of the grid. Unavailability of infrequent critical event field PMU data with component failures is driving the need to generate realistic synthetic PMU data for research. The synthetic data from power system simulation software often neglect noise profiles of received phasors, thus creating some discrepancies between real PMU data and synthetic ones. To address this issue, this work presents an initial study on the noise characteristics of PMUs, as well as presenting models for recreating their unique noise signatures. The proposed method, utilizing the Wasserstein generative adversarial network with gradient penalty (WGAN-GP) architecture, provides an excellent benchmark for matching the noise distribution. One can use a well-learned GAN model to draw noise signatures from a distribution that seemingly mirrors the real PMU noise distribution, while also being able to be detached from the PMU data once the training is done. Based on the observed results and employed data-driven methodology, it is expected that the proposed methods can be adapted to replicate the behavior of other sensors, providing research and other applications with a tool for data synthesis and sensor characterization.

PMU↗
DOE OSTI2023

An Incremental Gradient Method for Optimization Problems With Variational Inequality Constraints

We consider minimizing a sum of agent-specific nondifferentiable merely convex functions over the solution set of a variational inequality (VI) problem in that each agent is associated with a local monotone mapping. This problem finds an application in computation of the best equilibrium in nonlinear complementarity problems arising in transportation networks. We develop an iteratively regularized incremental gradient method where at each iteration, agents communicate over a directed cycle graph to update their solution iterates using their local information about the objective and the mapping. The proposed method is single-timescale in the sense that it does not involve any excessive hard-to-project computation per iteration. We derive nonasymptotic agent-wise convergence rates for the suboptimality of the global objective function and infeasibility of the VI constraints measured by a suitably defined dual gap function. Finally, the proposed method appears to be the first fully iterative scheme equipped with iteration complexity that can address distributed optimization problems with VI constraints over cycle graphs.

convergence↗
DOE OSTI2024

Field Harmonic Measurement of High-Gradient Permanent Magnet Quadrupoles for NSLS- II Upgrade

The National Synchrotron Light Source II (NSLS-II) is preparing for a major upgrade based on a new lattice concept, known as “Complex Bend,” to replace its existing double bend achromat (DBA) lattice. Here, this novel lattice will use high-gradient permanent magnet quadrupoles (PMQs) with small apertures in the range of 16–22 mm. A prototype complex bend branch with PMQs has been installed in the NSLS-II linac to evaluate the lattice's performance. An accurate and precise magnetic measurement system is required for field harmonic characterization and magnetic tuning of these PMQs. As part of this project a Rotating-Coil bench, based on a Printed Circuit Board (PCB) coil with a 12 mm diameter, has been configured to measure nine PMQs with bore diameter of 12.7 mm. This PCB coil has an active length of 270 mm, and is able to measure fields up to the 15th harmonic of the main field at a reference radius of 5 mm, while maintaining field quality levels within 10 ppm of the main field. This article presents the PCB coil setup, the harmonic measurement results and their temperature dependencies and repeatability. Simulation results on compensation of undesired multipoles using “IDBuilder”, a genetic algorithm-based optimizer code for magnetic tuning is included for field harmonic corrections.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Magnetic Measurements of Fermilab Rapid-Cycling Booster Gradient Magnets

Fermilab is upgrading its Booster synchrotron to increase ramp rate and intensity. This is part of the Proton Improvement Plan (PIP-II) that will allow the Main Injector to achieve proton beam power of 1.2 MW within the next few years. This upgrade includes running the 55-year-old Booster magnets at 20 Hz instead of the usual 15 Hz, and construction of some shorter and wider aperture versions of these combined-function gradient magnets. Magnetic measurements were performed to characterize the present 15 Hz AC performance, and then again with 20 Hz ramp cycle to ensure performance and compatibility in this new operating regime. A 3 m-long curved flat-coil was developed for these measurements using Printed Circuit Board (PCB) technology. The probe also has a separate 0.5 m-long body-field probe, allowing integral, body, and end fields to be measured across 100 mm of the magnet aperture. The sampling rate for these measurements during the AC cycle was 200 kHz, and field resolution was better than 0.01%. Details of the probe, measurements, and results are presented.

DiMarco, J. [Fermilab] (ORCID:0000000228400036)↗
DOE OSTI2025

Thermal Performance of a Conduction-Cooled CCT Dipole ReBCO Magnet: Several Cycles of Cool-Down and Thermal Gradient Measurements

Here, this paper presents experimental results from conduction-cooled thermal testing of a ReBCO canted cosine theta (CCT) magnet (C2), originally designed and fabricated at LBNL using CORC cables. While the performance of the coil under liquid helium and nitrogen environments has been previously established, this study explores its behavior under conduction cooling using a large test cryostat at The Ohio State University. The magnet, measuring 613 mm in length and weighing ∼75 kg, consists of four helical layers wound with ReBCO-based CORC wire and was thermally anchored to a copper cold ring supported by a G-10 strongback. Cooling was provided by two Sumitomo RDK-415D cryocoolers, offering a combined 3 W at 4.2 K and 150 W at 77 K. Multiple thermal cycles were performed, with cooldown durations of up to 45 hours. Final base temperatures of approximately 10.8 K (at the coil edge) and 12.0 K (at the coil center) were achieved, with an axial temperature difference of approximately 1.2 K. The warm-up period extended over approximately 24.6 hours. Voltage measurements from all four layers were recorded during cooldown and warmup. The system demonstrated stable cooldown performance, repeatable gradients, and good thermal anchoring. These results support the feasibility of conduction cooling in large-scale HTS magnets, aligning with broader goals for “green” cryogen-free accelerator technologies and paving the way for more sustainable, scalable, and energy-efficient high-field magnet systems in next-generation particle accelerators.

Canted cosine theta magnets↗