Engineering Papers⌕ Search

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.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2022

Spectroscopy of Twisted Bilayer Graphene Correlated Insulators

Here we analytically compute the scanning tunneling microscopy (STM) signatures of integer-filled correlated ground states of the magic angle twisted bilayer graphene (TBG) narrow bands. After experimentally validating the strong-coupling approach at ± 4 electrons/moiré unit cell, we consider the spatial features of the STM signal for 14 different many-body correlated states and assess the possibility of Kekulé distortion (KD) emerging at the graphene lattice scale. Remarkably, we find that coupling the two opposite graphene valleys in the intervalley-coherent (IVC) TBG insulators does not always result in KD. As an example, we show that the Kramers IVC state and its nonchiral U(4) rotations do not exhibit any KD, while the time-reversal-symmetric IVC state does. Our results, obtained over a large range of energies and model parameters, show that the STM signal and Chern number of a state can be used to uniquely determine the nature of the TBG ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Scale-dependent halo bias and the squeezed limit bispectrum in the presence of radiation

Here, we investigate the gravitational effect of large-scale radiation perturbations on small-scale structure formation. In addition to making the growth of matter perturbations scale dependent, the free streaming of radiation also affects the coupling between structure formation at small and large scales. We study this using separate universe N -body simulations to compute the (isotropized) squeezed limit matter bispectrum and the linear halo bias. Our results show that the scale dependence in the growth of long-wavelength matter perturbations, caused by radiation, translates into these quantities acquiring a nontrivial scale dependence at k ≲ 0.05 Mpc -1 . In a universe with radiation composed of cosmic microwave background photons and three species of massless neutrinos, the bias of halos with b = 2 at high k will decrease by 0.29%, 0.45%, and 0.8% between k = 0.05 Mpc -1 and k = 0.0005 Mpc -1 at redshifts z = 0, 1, and 3, respectively. For objects with b >> 1 , these differences approach 0.43%, 0.68%, and 1.2%, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Understanding Reliability Trade-Offs in 1T-nC and 2T-nC FeRAM Designs

Ferroelectric random access memory (FeRAM) is a promising candidate for energy-efficient nonvolatile memory, particularly for logic-in-memory and compute-in-memory (CIM) applications. Among the available cell architectures, One-Transistor–n-Capacitor (1T-nC) and two-transistor–n-capacitor (2T-nC) FeRAMs each offer distinct trade-offs in density, scalability, and reliability. In this work, we present a comparative study of these two architectures under both dimensional scaling ( XY/Z shrinkage) and vertical integration (increasing stacked capacitors per cell). Using technology computer-aided design (TCAD) and circuit-level simulations, we analyze how scaling impacts ferroelectric capacitance, parasitic coupling, and floating-node (FN) dynamics, which together dictate sense margin (SM) and read stability. A key mitigation strategy—floating unselected capacitors—is applied to both architectures, effectively decoupling the SM from the number of stacked capacitors and enabling tractable analysis across scaling regimes. Results show that 1T-nC suffers more from charge sharing with the bitline (BL), while 2T-nC benefits from transistor isolation and stronger low-voltage sensing at the cost of increased area. By systematically evaluating these behaviors across scaling directions, this work establishes the reliability trade-offs of 1T-nC and 2T-nC cells and provides design guidelines for high-density, vertically integrated FeRAM systems.

1T-nC↗
DOE OSTI2025

The Dynamical Role of Optical Phonons and Sublattice Screening in a Solid-State Ion Conductor

Solid-state electrolytes (SSEs) require ionic conductivities that are competitive with liquid electrolytes to realize applications in all-solid-state batteries. Although candidate SSEs have been discovered, the underlying mechanisms enabling superionic conduction (>1 mS cm –1 ) remain elusive. In particular, the role of ultrafast lattice dynamics in mediating ion migration, which involves couplings between ions, phonons, and electrons, is rarely explored experimentally at their corresponding time scales. To investigate the complex contributions of coupled lattice dynamics on ion migration, we modulate the charge density occupations within the crystal framework and then measure the time-resolved change in impedance on picosecond time scales for a candidate SSE, Li 0.5 La 0.5 TiO 3 (LLTO). Upon perturbation, we observe enhanced ion migration at ultrafast time scales. The respective transients match the time scales of optical and acoustic phonon vibrations, suggesting their involvement in ion migration. We further computationally evaluate the effect of a charge transfer from the O 2p to the Ti 3d band on the electronic and physical structure of LLTO. We hypothesize that the charge-transfer excitation distorts the TiO 6 polyhedra by altering the local charge density occupancy of the hopping site at the migration pathway saddle point, thereby causing a reduction in the migration barrier for the Li + hop. We rule out the contribution of photogenerated electron carriers and laser heating. Overall, our investigation introduces a new spectroscopic tool to probe fundamental ion hopping mechanisms transiently at ultrafast time scales, which has previously only been achieved in a time-averaged manner or solely via computational methods.

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

Investigation of Surface and Marine-Cloud Coupling and its Impact on Cloud Droplet Number Concentrations and Cloud Cover Over the Southern Ocean

The proposed work involves characterizing and quantifying both the thermodynamic and dynamical coupling of marine low clouds (MLC) with the sea surface using Atmospheric Radiation Measurement (ARM) observations during MARCUS field campaign and an LES model. ARM data from multiple sensors will be used (e.g., Doppler cloud radar, ceilometer, and microwave radiometer) to characterize the MLC-surface coupling by virtue of the vertical structure and integrated quantities of boundary-layer clouds, aerosols, as well as atmospheric profiles and surface meteorology. We have used a combination of case studies and statistics-based composite analyses to find any linkages between large-scale dynamics, MLC-surface coupling, and cloud and boundary layer properties, the sequence of which reflects the chain of causality. An LES model with explicit aerosol physics, which includes the cycle of aerosols by being consumed as CCN and be regenerated following cloud droplet evaporation, has been run to determine the sources and sinks of Nd and their dependence on the degree of MLC-surface coupling. We have examined the systematic differences in both Nd and cloud occurrence between the two clusters under different meteorological conditions and further examine their respective roles, as well as causal relationships by means of LES modeling. This study helped improve our understanding of the ACI by differentiating the dynamic role of the coupling and cloud physics denoted by Nd, bridging the linkage in the chain toward understanding mechanisms governing the persistence of MLC over the Southern Ocean, solving the long-lasting problem of the cloud cover underestimation over the SO by GCMs. Ample ARM data and LES model have been employed to achieve the objectives of the study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗
DOE OSTI2020

Stochastically Realized Observables for Excitonic Molecular Aggregates

In this work, we show that a stochastic approach enables calculations of the optical properties of large 2-dimensional and nanotubular excitonic molecular aggregates. Previous studies of such systems relied on numerically diagonalizing the dense and disordered Frenkel Hamiltonian, which scales approximately as $\mathcal{O}$($N$ 3 ) for $N$ dye molecules. Our approach scales much more efficiently as $\mathcal{O}$($N$ log($N$)), enabling quick study of systems with a million of coupled molecules on the micrometer size scale. We calculate several important experimental observables, including the optical absorption spectrum and density of states, and develop a stochastic formalism for the participation ratio. Quantitative agreement with traditional matrix diagonalization methods is demonstrated for both small- and intermediate-size systems. The stochastic methodology enables the study of the effects of spatial-correlation in site energies on the optical signatures of large 2D aggregates. Our results demonstrate that stochastic methods present a path forward for screening structural parameters and validating experiments and theoretical predictions in large excitonic aggregates.

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

Mechanical quantum sensing in the search for dark matter

Numerous astrophysical and cosmological observations are best explained by the existence of dark matter, a mass density which interacts only very weakly with visible, baryonic matter. Here, searching for the extremely weak signals produced by this dark matter strongly motivate the development of new, ultra-sensitive detector technologies. Paradigmatic advances in the control and readout of massive mechanical systems, in both the classical and quantum regimes, have enabled unprecedented levels of sensitivity. In this white paper, we outline recent ideas in the potential use of a range of solid-state mechanical sensing technologies to aid in the search for dark matter in a number of energy scales and with a variety of coupling mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Coupling physics in machine learning to predict properties of high-temperatures alloys

Abstract High-temperature alloy design requires a concurrent consideration of multiple mechanisms at different length scales. We propose a workflow that couples highly relevant physics into machine learning (ML) to predict properties of complex high-temperature alloys with an example of the 9–12 wt% Cr steels yield strength. We have incorporated synthetic alloy features that capture microstructure and phase transformations into the dataset. Identified high impact features that affect yield strength of 9Cr from correlation analysis agree well with the generally accepted strengthening mechanism. As a part of the verification process, the consistency of sub-datasets has been extensively evaluated with respect to temperature and then refined for the boundary conditions of trained ML models. The predicted yield strength of 9Cr steels using the ML models is in excellent agreement with experiments. The current approach introduces physically meaningful constraints in interrogating the trained ML models to predict properties of hypothetical alloys when applied to data-driven materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Coupling turbulent flow with blade aeroelastics and control modules in large-eddy simulation of utility-scale wind turbines

We present a large-eddy simulation framework capable of control co-design of large wind turbines, coupling the turbulent flow environment with blade aeroelastics and turbine controllers. The geometry and aerodynamics of the rotor blades and the turbine nacelle are parameterized using an actuator surface model. The baseline collective pitch control and individual pitch control (IPC) algorithms, consisting of a single-input, single-output proportional–integral controller and two integral controllers, respectively, are incorporated into the simulation framework. Furthermore, a second-order model based on the Euler–Bernoulli beam theory is implemented to describe the blade deformation. Simulations are carried out to investigate the impact of collective and individual pitch control strategies on the deflection of turbine blades. Our results show that the IPC reduces the blade tip deflection fluctuations in the out-of-plane direction, while the fluctuations of the blade tip deflection along the in-plane direction are barely affected by the IPC. Furthermore, the blade out-of-plane deformation fluctuation is underestimated by the one-way coupling approach compared to the two-way coupling approach. Furthermore, the findings of this study reveal the importance of advanced control systems in reducing the dynamic loads on wind turbine blades and underscore the potential of control co-design to reduce the levelized cost of wind energy.

42 ENGINEERING↗
DOE OSTI2021

Resurrecting low-mass axion dark matter via a dynamical QCD scale

In the framework where the strong coupling is dynamical, the QCD sector may confine at a much higher temperature than it would in the Standard Model, and the temperature-dependent mass of the QCD axion evolves in a non-trivial way. We find that, depending on the evolution of Λ QCD , the axion field may undergo multiple distinct phases of damping and oscillation leading generically to a suppression of its relic abundance. Such a suppression could therefore open up a wide range of parameter space, resurrecting in particular axion dark-matter models with a large Peccei-Quinn scale f a $\gg$ 10 12 GeV, i.e., with a lighter mass than the standard QCD axion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: LDRD Final Project Poster

Electric-field assisted sintering (EFAS) manufacturing techniques reduce energy usage requirements, fabricate near-net-shape parts, and produce bulk nanostructured materials through rapid heating and cooling rates achievable through Joule heating. Process parameters of pulsed DC voltage and constant pressure are applied through stainless steel rams and graphite tooling to the ceramic yttrium oxide powder to be sintered. Our high-fidelity modeling code application advances manufacturing for extreme environments by enabling process-informed design and improving microstructure consistency through multiphysics multiscale process-structure-property-performance correlation simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Subtropical Eastern North Pacific SST bias in earth system models

Karthik Balaguru (PNNL), Luke Van Roekel (LLNL), L. Ruby Leung (PNNL) and Milena Veneziani (LLNL) Abstract: This study systematically evaluates the warm SST bias in the Subtropical Eastern North Pacific (SENP), a problem plaguing most CMIP6 models, using the Energy Exascale Earth System Model version 1 (E3SM). Compared to the model at its standard resolution (1o atmosphere, 30-60 km ocean), the SST bias improves in several regions at high resolution (0.25o atmosphere, 18-6 km ocean) except in subtropical eastern boundary regions, such as the SENP. A mixed layer heat budget analysis reveals that errors in the surface radiative fluxes and horizontal heat advection are primarily responsible for the bias. Analysis of HighResMIP models indicates that the lack of improvement in the SENP SST bias with resolution is not native to E3SM alone. A larger CMIP6 inter-model correlation of SST bias with biases in surface winds than radiation highlights the important role of large-scale circulation biases in the SENP SST bias in coupled models.

Balaguru, Karthik↗
DOE OSTI2020

Explicitly correlated coupled cluster method for accurate treatment of open-shell molecules with hundreds of atoms

We present a near-linear scaling formulation of the explicitly correlated coupled-cluster singles and doubles with the perturbative triples method [CCSD(T)F12¯] for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g., the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the C160H322 n-alkane, with the measured size complexity of ≈1.2). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)F12¯ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580–6602 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)], the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)F12¯ yielded a RMSD of ∼0.4 kcal/mol with respect to the CBS CCSD(T) estimate.

Kumar, Ashutosh (ORCID:0000000175896030)↗