SEARCH · Engineering Papers

Results for “scale coupling”

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 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗

DOE OSTI2021

On steady alternate bars forced by a localized asymmetric drag distribution in erodible channels

Studying the effect of different in-stream fluvial turbines siting on river morphodynamics allowed us to witness the onset of a time-averaged, large-scale, alternate distortion of bed elevations, which could not be exclusively related to the turbine rotor blockage. The longitudinal profiles of this two-dimensional bathymetric perturbation resemble those of steady fluvial bars. In this contribution we generalize the problem addressing a spatially impulsive, asymmetric distribution of drag force in the channel cross-section. Furthermore, this is experimentally investigated through the deployment of differently sized grids perpendicular to the flow, and analytically explored as a finite perturbation of an open channel flow over an erodible sediment layer, as described by a coupled flow–sediment shallow water equation. The steady solutions of this fluvial morphodynamic problem, physically represented by alternate bars scaling with the channel width, highlight the importance of the resonant conditions in defining the spatial extent of the bed deformation. The equations further suggest that in very shallow flows any asymmetric obstruction may lead to an upstream propagation of the steady bars, consistent with previous studies on the effects of channel curvature. In broad terms, this study provides the preliminary framework to control the onset of river meandering through imposed finite perturbations of the cross-section. In a more applied sense, it provides a tool to predict non-local scour–deposition patterns associated with the deployment of energy converters or other flow obstructions.

42 ENGINEERING↗

DOE OSTI2026

Neural-Network Inverse Design of SRF Cavities and Transmons for Bosonic Quantum Computation

Three-dimensional superconducting radio-frequency (SRF) cavities provide exceptionally long-lived electromagnetic modes and, when coupled to nonlinear elements such as transmon qubits, become promising architectures for bosonic quantum information processing. The inverse design of such systems, i.e., recovering device geometries that produce specified electromagnetic and coupling targets, is generally a one-to-many problem. The qubit-cavity coupling strength depends sensitively on both the transmon geometry and its position within the cavity's electromagnetic field. As these systems scale up and their design parameter spaces grow, the cost of conventional iterative simulation becomes prohibitive. We present two deep neural network (DNN) approaches that address this inverse-design problem at complementary levels of the design stack. The first proposes SRF cavity geometries that produce target cavity observables. The second proposes transmon qubit designs that produce target qubit-cavity parameters - the coupling rate, qubit frequency, and anharmonicity $(g, ν_q, α)$. The recovered candidate designs match the targets to within ~5% (cavity) and ~2% (transmon), confirmed by end-to-end re-simulation. Both approaches map desired device behavior directly to candidate designs, a fast alternative to the iterative simulation studies usually required.

Yaker, Joseph [Fermilab; Northwestern U.]↗

DOE OSTI2021

National direct-drive program on OMEGA and the National Ignition Facility

A major advantage of the laser direct-drive (DD) approach to ignition is the increased fraction of laser drive energy coupled to the hot spot and relaxed hot-spot requirements for the peak pressure and convergence ratios relative to the indirect-drive approach at equivalent laser energy. With the goal of a successful ignition demonstration using DD, the recently established national strategy has several elements and involves multiple national and international institutions. These elements include the experimental demonstration on OMEGA cryogenic implosions of hot-spot conditions relevant for ignition at MJ-scale energies available at the National Ignition Facility (NIF) and developing an understanding of laser-plasma interactions and laser coupling using DD experiments on the NIF. DD designs require reaching central stagnation pressures in excess of 100 Gbar. The current experiments on OMEGA have achieved inferred peak pressures of 56 Gbar (Regan et al 2016 Phys. Rev. Lett. 117 025001). Extensive analysis of the cryogenic target experiments and two- and three-dimensional simulations suggest that power balance, target offset, and target quality are the main limiting factors in target performance. In addition, cross-beam energy transfer (CBET) has been identified as the main mechanism reducing laser coupling. Reaching the goal of demonstrating hydrodynamic equivalence on OMEGA includes improving laser power balance, target position, and target quality at shot time. CBET must also be significantly reduced and several strategies have been identified to address this issue. (2016-01-01)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Dynamical attractors in contracting spacetimes dominated by kinetically coupled scalar fields

We present non-perturbative numerical relativity simulations of slowly contracting spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear σ model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble density fluctuations fully consistent with cosmic microwave background observations. We show that the non-linear evolution rapidly approaches a homogeneous, isotropic and flat Friedmann-Robertson-Walker (FRW) geometry for a wide range of inhomogeneous and anisotropic initial conditions. Ultimately, we find, the kinetic coupling causes the evolution to deflect away from flat FRW and towards a novel Kasner-like stationary point, but in general this occurs on time scales that are too long to be observationally relevant.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗

DOE OSTI2023

Physics Opportunities at a Beam Dump Facility at PIP-II and Beyond

The accelerator complex at Fermilab is currently undergoing improvements which will increase the available beam power to the complex and is known as Proton Improvement Plan-II (PIP-II). The PIP-II Linac is slated for operation near the end of this decade and will be the main proton driver for Fermilab experiments moving forward and provide the beam to LBNF/DUNE. However, the DUNE physics program requires only ~1% of the available protons provided by PIP-II. The Accelerator Complex Enhancement, or ACE, will provide further up grades in the 2030s in the form of increased power to LBNF along with a replacement for the Fermilab Booster which could also include an accumulator ring. A workshop was recently held at Fermilab in order to explore the physics possible when PIP-II is coupled with a fixed target or beam dump. The beam dump scenario pro vides an exciting opportunity to search for dark sector physics across detector threshold energy scales with examples being accelerator-produced dark matter, active-to-sterile neutrino oscillations, millicharged parti cles, and axion-like particles, which can be produced in the proton collisions with a fixed target. Additionally, experiments located at a PIP-II facility that is coupled to an accumulator ring would also be able to study coherent elastic neutrino-nucleus scattering, or CEvNS. In this talk, I will summarize the physics possible at a PIP-II beam dump facility and sensitivities to different dark sector physics and other beyond the Standard Model searches with detector concepts spanning from an eV-scale to MeV-scale detection threshold.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Ponderomotive force driven density modifications parallel to B on the LAPD

The ponderomotive force has previously been identified as a possible driver of observed density modifications close to radio frequency actuators during operation. This nonlinear force redistributes density in regions of gradients in the magnitude of an oscillating electric field and describes the influence of the fast time scale RF wave dynamics on slow time scale plasma transport. Depletion of the saturation current (a proxy for the density) measured at the Large Plasma Device (LAPD) was 30–35% during ion cyclotron range of frequencies operation. A coupled 1D plasma transport and cold plasma frequency domain wave solver was developed to self-consistently describe ponderomotive effects and was used to compare with results obtained from the LAPD experiment. The scaled current density driver for the wave model yielded an RF B field in close agreement with two components of the experimental data. However, the 1D parallel model did not accurately reproduce the amplitude or spatial distribution observed in experimental measurements of By. Within the limitations of the 1D model, initial simulation results showed that the ponderomotive force depleted up to 8% for high power (1 MW) and around 1% for the experimental power of 120 kW. This could suggest that the ponderomotive force is not the main driver of density modification for the LAPD experiments presented in this paper. Higher fidelity tools of at least 2D will be required to give a more realistic description of the RF E fields and the effect of the ponderomotive force on the LAPD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Powering ≈50 µm Motion by a Molecular Event in DNA Crystals

A major challenge in material design is to couple nanoscale molecular and supramolecular events into desired chemical, physical, and mechanical properties at the macroscopic scale. Here, a novel self-assembled DNA crystal actuator is reported, which has reversible, directional expansion and contraction for over 50 μm in response to versatile stimuli, including temperature, ionic strength, pH, and redox potential. The macroscopic actuation is powered by cooperative dissociation or cohesion of thousands of DNA sticky ends at the designed crystal contacts. The increase in crystal porosity and cavity in the expanded state dramatically enhances the crystal capability to accommodate/encapsulate nanoparticles/proteins, while the contraction enables a “sponge squeezing” motion for releasing nanoparticles. This crystal actuator is envisioned to be useful for a wide range of applications, including powering self-propelled robotics, sensing subtle environmental changes, constructing functional hybrid materials, and working in drug controlled-release systems.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Higgs Inflation and the Electroweak Gauge Sector

Abstract We introduce a method that allows the Higgs to be the inflaton. The Higgs is considered as a pseudo‐Nambu‐Goldstone (pNG) boson of a global coset symmetry , which is spontaneously breaks at an energy scale . A suitable Chern−Simons (CS) interaction is given to it, with representing the dimensionless CS coupling strength and an decay constant. As a result, slow‐roll inflation occurs via ‐induced friction down a steep sinusoidal potential. To obey electroweak symmetry, the lowest‐order CS interaction is required to be quadratic in the Higgs, with the coupling strength . Higher‐order interaction terms keep the full Lagrangian nearly invariant under the approximate pNG shift symmetry. Employing the simplest symmetry coset , ‐folds of inflation occur when . Successfully explaining inflation necessitates small values of the decay constant, ; this in turn requires large , which is ruled out by electric dipole measurements. Although the electroweak hierarchy problem while achieving successful inflation, the real benefit is found in providing a different path to identifying the Higgs as the inflaton, outside the standard modified‐gravity framework.

Alexander, Stephon [Brown Center for Theoretical P↗

DOE OSTI2022

Inverted c-functions in thermal states

We first compute the effect of a chiral anomaly, charge, and a magnetic field on the entanglement entropy in $\mathcal{N}$ = 4 Super-Yang-Mills theory at strong coupling via holography. Depending on the width of the entanglement strip the entanglement entropy probes energy scales from the ultraviolet to the infrared energy regime of this quantum field theory (QFT) prepared in a given state. From the entanglement entropy, we compute holographic c-functions and demonstrate an inverted c-theorem for them. That is, these c-functions in generic thermal states monotonically increase towards the infrared (IR) energy regime. This is in contrast to the c-functions in vacuum states which decrease along the renormalization group flow towards the IR regime of a renormalizable QFT. Furthermore, in thermal states and in the IR limit, the c-functions behave thermally, growing proportionally to the value of the thermal entropy. The chiral anomaly affects the c-functions more in the IR regime, and its effect is peaked at an intermediate value of the magnetic field at a fixed chemical potential and temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Mass production of high-performance single atomic FeNC electrocatalysts via sequenced ultrasonic atomization and pyrolysis process

Mass production of highly efficient, durable, and inexpensive single atomic catalysts is currently the major challenge associated with the oxygen reduction reaction (ORR) for fuel cells. In this study, we develop a general strategy that uses a simple ultrasonic atomization coupling with pyrolysis and calcination process to synthesize single atomic FeNC catalysts (FeNC SACs) at large scale. The microstructure characterizations confirm that the active centers root in the single atomic Fe sites chelating to the four-fold pyridinic N atoms. The identified specific Fe active sites with the variable valence states facilitate the transfer of electrons, endowing the FeNC SACs with excellent electrochemical ORR activity. The FeNC SACs were used as cathode catalysts in a homemade Zn-air battery, giving an open-circuit voltage (OCV) of 1.43 V, which is substantially higher than that of commercial Pt/C catalysts. This study provides a simple approach to the synthesis of single atomic catalysts at large scale.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effect of chemical disorder on the electronic stopping of solid solution alloys

The electronic stopping power of nickel-based equiatomic solid solutions alloys NiCr, NiFe and NiCo for protons and alpha projectiles is investigated in detail using real-time time-dependent density functional theory over a wide range of velocities. Recently developed numerical electronic structure methods are used to probe fundamental aspects of electron-ion coupling non-perturbatively and in a fully atomistic context, capturing the effect of the atomic scale disorder. The effects of particular electronic band structures and density of states reflect in the low velocity limit behavior. We compare our results for the alloys with those of a pure nickel target to understand how alloying affects the electronic stopping. We discover that NiCo and NiFe have similar stopping behavior as Ni while NiCr has an asymptotic stopping power that is more than a factor of two larger than its counterparts for velocities below 0.1 a.u.. Overall, we show that the low-velocity limit of electronic stopping power can be manipulated by controlling the broadening of the d-band through the chemical disorder. In this regime, the Bragg’s additive rule for the stopping of composite materials also fails for NiCr.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Steric effects of central dogma processes on the compaction and segregation of bacterial nucleoids

The bacterial cytoplasm is characterized by a distinctive membrane-less organelle, the nucleoid, which harbors the chromosomal DNA. Here, we investigate the steric effects of dynamic processes associated with transcription and translation on the structure of this organelle using coarse-grained molecular dynamics simulations that incorporate out-of-equilibrium reactions. Our model captures the scale of the entire cell and incorporates a reaction-diffusion system for ribosomes and polyribosomes, coupling their nonequilibrium kinetics to DNA excluded volume interactions. Our findings demonstrate that out-of-equilibrium reactions increase the size of the nucleoid and the number of ribosomes in this subcellular region. In addition, we show that nucleoid size is proportional to transcriptional activity. Our model reproduces the time-dependent change in nucleoid size observed in rifampicin treatment experiments, where the pool of polyribosomes is depleted. Furthermore, we find that these active processes are essential for complete sister chromosome separation and correct nucleoid positioning within the cell. Overall, our study reveals the effects of the central dogma processes on the internal organization and localization of bacterial nucleoids.

Chang, Mu-Hung [Univ. of Tennessee, Knoxville, TN ↗

DOE OSTI2025

Multiscale electro-chemo-mechanical model of high-capacity amorphous silicon anode material in sodium-ion batteries

The growing interest in sodium-ion batteries (SIBs) is fueled by their abundant resources and environmentally friendly nature, with amorphous silicon (a-Si) emerging as a promising anode material for enhancing capacity. However, the key challenge lies in sustaining reversible capacity during cycling. Here, in this work, we developed a multiscale electrochemical model incorporating an a-Si anode to elucidate the performance parameters of SIBs. Additionally, we integrated an electro-chemo-mechanical model at the particle level to capture stress generation, an essential factor in the degradation of high-capacity electrodes. Unlike existing models, our approach accounts for large-deformation chemo-mechanics at the particle scale and includes simulations under varying charge rates to explore multiscale behavior. The results reveal that coupled sodiation significantly prolongs complete cycling times and influences discharge dynamics, indicating that neglecting this coupling leads to an underestimation of actual capacity. Furthermore, we observed pronounced polarization effects at higher charge rates, resulting in heterogeneous stress distributions across the anode. With the identification of critical failure parameters for both active particles and binder materials, offering novel insights for mitigating degradation in high-capacity electrode systems.

Amorphous silicon anode↗

DOE OSTI2023

Multi–length scale characterization of point defects in thermally oxidized, proton irradiated iron oxides

A key for the success of safe nuclear power generation system is to consider structural materials that are economical, meet mechanical property needs, possess good corrosion resistance, and are radiation tolerant. Nevertheless, fundamental insights that elucidate the details of radiation damage on materials corrosion performance are lacking. This includes the behavior of surface oxides which often regulate corrosion. For example, it is unclear how non-equilibrium point defects, oxide structure, mass transport in oxides, and subsequent oxidation behavior are altered by the radiation. Here, in this work, some of the effects of proton irradiation on the attributes of point defects, iron oxide microstructures, and the physical nature of the oxidation product were correlated with corrosion behavior. Iron oxides, fabricated by thermal oxidation in air at 400°C and 800°C for 1 h, were subjected to 200 keV, 0.03 dpa (displacements per atom) of proton irradiation, and subjected to corrosion reactivity assessment using AC and DC electrochemical methods. Experimental methods that target materials information at different length scales, such as positron annihilation spectroscopy (atomistic), transmission electron microscopy (mesoscopic), and electrochemical methods (macroscopic), were coupled to shed light on the impact of radiation-induced defect modifications and structural changes in oxides on corrosion reactivity which ultimately affects durability in harsh environments.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Explicit modeling of pebble temperature in the porous-media model for pebble-bed reactors

In this study, we developed a multiscale model to include an explicit pebble-temperature model nested in the porous-media model for pebble-bed reactor applications. The multiscale solid-phase energy balance model, including the pebble surface energy balance equation and an explicit modeling of pebble temperature, can predict the macroscopic (pebble bed) and microscopic (pebble) temperature distributions under both steady-state and transient conditions. The proposed multiscale model is solved in a fully coupled manner using the Newton- Krylov method, and therefore iterations between the macroscopic (pebble-bed-scale) and microscopic (pebble- scale) model are avoided. Extensive code verifications, validation, and demonstrations have been performed for this newly developed model. By explicitly modeling pebble temperatures, this new model addresses a major deficiency of the basic porous-media model, which assumes homogeneous solid-phase temperature and is not appropriate for pebble-bed reactor design and safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗