Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic boundaries”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Lagrangian Study of Interfaces at the Edges of Cumulus Clouds

Abstract Interfaces at the edge of an idealized, nonprecipitating, warm cloud are studied using direct numerical simulation (DNS) complemented with a Lagrangian particle tracking routine. Once a shell has formed, four zones can be distinguished: the cloud core, visible shell, invisible shell, and the environment. The union of the visible and invisible regions is the shell commonly referred to in literature. The boundary between the invisible shell and the environment is the turbulent–nonturbulent interface (TNTI), which is typically not considered in cloud studies. Three million particles were seeded homogeneously across the domain and properties were recorded along individual trajectories. The results demonstrate that the traditional cloud boundary (separating cloudy and noncloudy regions using thresholds applied on liquid condensate or updraft velocity) are some distance away from the TNTI. Furthermore, there is no dynamic difference between the traditional liquid-condensate boundary and the region extending to the TNTI. However, particles crossing the TNTI exhibit a sharp jump in enstrophy and a smooth increase in buoyancy. The traditional cloud boundary coincides with the location of minimum buoyancy in the shell. The shell premixes the entraining and detraining air and analysis reveals a highly skewed picture of entrainment and detrainment at the traditional cloud boundary. A preferential entrainment of particles with velocity and specific humidity higher than the mean values in the shell is observed. Large-eddy simulation of a more realistic setup detects an interface with similar properties using the same thresholds as in the DNS, indicating that the DNS results extrapolate beyond their idealized conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗
DOE OSTI2024

Boundary correlators and the Schwarzian mode

The effective low temperature dynamics of near-extremal black holes is governed by the quantum fluctuations of the Schwarzian mode of JT gravity. Utilizing as a proxy a planar charged black hole in asymptotically Anti-de-Sitter spacetime, we investigate the effects of these fluctuations on a probe scalar field. The corresponding holographic real-time boundary correlators are computed following a holographic renormalization procedure, using the dubbed gravitational Schwinger-Keldysh geometry (grSK) and known exact results of boundary correlators from the near-horizon region. This analysis gives rise to a retarded Green’s function that decays as a power law for late Lorentzian times. Its analytic structure indicates the presence of a branch cut in the complex frequency domain at finite temperature. These features are a non-perturbative hallmark that prevails as long as the planar transverse space is kept compact.

2D Gravity↗
DOE OSTI2025

An open retail boundary dataset for South Korea using open data and computer vision technique

Although delineating retail boundaries is important to explore and comprehend the dynamics of the retail sector, it is hard to find studies specifically addressing it in the South Korean context. This study fills this gap by proposing new retail boundaries across South Korea. To achieve this goal, we employed a variety of retailers and building datasets and proposed a unique computer vision-based framework with a deep ensemble voting technique. As a result, we delineated 6,636 distinct retail boundaries that were validated against existing reference retail boundaries. These newly delineated retail boundaries provide valuable insights for researchers, governments, and other relevant stakeholders by enhancing their understanding of retail geography. This dataset can be used as a foundational resource for analyses on topics such as pandemic recovery, retail gentrification, and the resilience of retail spaces in response to e-commerce growth, ultimately contributing to more robust retail sector research in South Korea.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Real-Time Modal Analysis of Electric Power Grids - The Need for Dynamic State Estimation

We articulate the reason why dynamic state estimation is needed to push the boundaries of the modal analysis of electric power grids in real-time operation. Then, we demonstrate how to unravel linear and nonlinear modes by using the extended dynamic mode decomposition along with estimates of the synchronous generators’ rotor angles and rotor speed deviations from nominal speed. The estimated modes are associated with electromechanical oscillations that take place continuously in electric power grids because of imbalances between power generation and demand. The numerical simulations are performed on a synthetic, albeit realistic, 2,000-bus network that was designed to resemble the electric power grid of Texas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Real-Time Modal Analysis of Electric Power Grids - The Need for Dynamic State Estimation: Preprint

We articulate the reason why dynamic state estimation is needed to push the boundaries of the modal analysis of electric power grids in real-time operation. Then, we demonstrate how to unravel linear and nonlinear modes by using the extended dynamic mode decomposition along with estimates of the synchronous generators' rotor angles and rotor speed deviations from nominal speed. The estimated modes are associated with electromechanical oscillations that take place continuously in electric power grids because of imbalances between power generation and demand. The numerical simulations are performed on a synthetic, albeit realistic, 2,000-bus network that was designed to resemble the electric power grid of Texas.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2023

Friction-mediated phase transition in confined active nematics

Using a minimal continuum model, we investigate the interplay between circular confinement and substrate friction in active nematics. Upon increasing the friction from low to high, we observe a dynamical phase transition from a circulating flow phase to an anisotropic flow phase in which the flow tends to align perpendicular to the nematic director at the boundary. We demonstrate that both the flow structure and dynamic correlations in the latter phase differ from those of an unconfined, active turbulent system and may be controlled by the prescribed nematic boundary conditions. Here, our results show that substrate friction and geometric confinement act as valuable control parameters in active nematics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Excellent ballistic impact resistance of Al 0.3 CoCrFeNi multi-principal element alloy with unique bimodal microstructure

Multi-principal element alloys represent a new paradigm in structural alloy design with superior mechanical properties and promising ballistic performance. Here, the mechanical response of Al 0.3 CoCrFeNi alloy, with unique bimodal microstructure, was evaluated at quasistatic, dynamic, and ballistic strain rates. The microstructure after quasistatic deformation was dominated by highly deformed grains. High density of deformation bands was observed at dynamic strain rates but there was no indication of adiabatic shear bands, cracks, or twinning. The ballistic response was evaluated by impacting a 12 mm thick plate with 6.35 mm WC projectiles at velocities ranging from 1066 to 1465 m/s. The deformed microstructure after ballistic impact was dominated by adiabatic shear bands, shear band induced cracks, microbands, and dynamic recrystallization. The superior ballistic response of this alloy compared with similar Al x CoCrFeNi alloys was attributed to its bimodal microstructure, nano-scale L1 2 precipitation, and grain boundary B2 precipitates. Deformation mechanisms at quasistatic and dynamic strain rates were primarily characterized by extensive dislocation slip and low density of stacking faults. Deformation mechanisms at ballistic strain rates were characterized by grain rotation, disordering of the L1 2 phase, and high density of stacking faults.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Rapid Quantum Ground State Preparation via Dissipative Dynamics

Inspired by natural cooling processes, dissipation has become a promising approach for preparing low-energy states of quantum systems. However, the potential of dissipative protocols remains unclear beyond certain commuting Hamiltonians. This work provides significant analytical and numerical insights into the power of dissipation for preparing the ground state of noncommuting Hamiltonians. For quasi-free dissipative dynamics, including certain 1D spin systems with boundary dissipation, our results reveal a new connection between the mixing time in trace distance and the spectral properties of a non-Hermitian Hamiltonian, leading to an explicit and sharp bound on the mixing time that scales polynomially with system size. For more general spin systems, we develop a tensor network-based algorithm for constructing the Lindblad jump operator and for simulating the dynamics. Using this algorithm, we demonstrate numerically that dissipative ground state preparation protocols can achieve rapid mixing for certain 1D local Hamiltonians under bulk dissipation, with a mixing time that scales logarithmically with the system size. We then prove the rapid mixing result for certain weakly interacting spin and fermionic systems in arbitrary dimensions, extending recent results for high-temperature quantum Gibbs samplers to the zero-temperature regime. Together, these results show that dissipation can be a powerful tool for ground state preparation, with potential applications across condensed matter physics, quantum materials science, and beyond.

decoherence↗
DOE OSTI2020

Closed-form solutions towards real-time energy-optimal and collision-free speed planning with limited information

Under real-world driving conditions, available information is limited due to communication range, and connected and automated vehicles (CAVs) must plan and follow energy-optimal and collision-free speed trajectory with a high updating rate. To this end, this paper presents a holistic design of a speed planning algorithm using analytical closed-form optimal solutions. Using the simplest vehicle model, we derive closed-form solutions as a function of boundary conditions, and summarize them without and with state constraints imposed by speed limits and vehicle in front. Then, we introduce multiple modes (e.g., eco-approach to traffic signal) to properly respond dynamically changing situations, and present how to set boundary conditions for each mode, while guaranteeing nature of globally optimal solutions. Finally, we perform a large-scale simulation study to identify and quantify the energy impacts of CAVs for real-world driving routes. A simple but effective algorithm based on closed-form solutions shows a significant energy saving potential in comparison with human-driven vehicle and adaptive cruise controlled vehicle with connectivity.

Han, Jihun↗
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase segregation dynamics in mixed-halide perovskites revealed by plunge-freeze cryo-electron microscopy

Mixed-halide lead perovskites, with photoexcited charge-carrier properties suitable for high-efficiency photovoltaics, hold significant promise for high-efficiency tandem solar cells. However, phase segregation under illumination, where an iodide-rich phase forms carrier trap states, remains a barrier to application. This study employs plunge-freeze cryoelectron microscopy to visualize nanoscale phase segregation dynamics in CsPb(Br x I 1–x ) 3 films. By rapidly freezing the illuminated samples, we preserve transient photoexcited ion distributions for high-resolution structural and compositional analysis at the nanoscale. Cryogenic scanning transmission electron microscopy (STEM) techniques (electron energy loss spectroscopy [EELS] and 4D-STEM) captured the dynamics of photo-induced iodine migration from grain boundaries to centers, identified the buildup of anisotropic strain, and captured the heterogeneous evolution of this process within a single grain. These findings provide insight into microscopic phase segregation mechanisms and their dynamics, enhancing our understanding of mixed-halide perovskite photostability.

42 ENGINEERING↗
DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Microsecond Molecular Dynamics Simulations of Proteins Using a Quasi-Equilibrium Solvation Shell Model

In this work, we describe the development and implementation of a quasi-equilibrium hydration shell model of biomolecular solvation with adaptive boundaries. Applying the model to microsecond-long molecular dynamics simulations of several protein systems of varying complexity, we find that the model simulation results are of comparable quality to those obtained from simulations of fully solvated systems, but at a reduced computational cost. We discuss the dominant sources of error in the model and outline directions for future improvements.

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

Dislocation dynamics in polycrystalline materials

An extension of the dislocation dynamics method to polycrystalline materials is presented with emphasis on handling interactions between dislocations and grain boundaries. The key new features of the method include: (1) a grain boundary detection algorithm, (2) new dislocation emission criteria based on power dissipation and area growth, (3) mobility of dislocations at grain boundaries, and (4) extension of topological operations e.g. mesh adaption, collision, and dissociation to handle dislocations interacting with grain boundaries. Furthermore, verification and validation examples are executed to compare our calculations with existing analytic methods, molecular dynamics simulations, and experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Causality-respecting adaptive refinement for PINNs: enabling precise interface evolution in phase field modeling

Physics-informed neural networks (PINNs) have emerged as a powerful tool for solving physical systems described by partial differential equations (PDEs). However, their accuracy in dynamical systems, particularly those involving sharp moving boundaries with complex initial morphologies, remains a challenge. Here, this study introduces an approach combining residual-based adaptive refinement (RBAR) with causality-informed training to enhance the performance of PINNs in solving spatio-temporal PDEs. Our method employs a three-step iterative process: initial causality-based training, RBAR-guided domain refinement, and subsequent causality training on the refined mesh. Applied to the Allen-Cahn equation, a widely-used model in phase field simulations, our approach demonstrates significant improvements in solution accuracy and computational efficiency over traditional PINNs. Notably, we observe an ‘overshoot and relocate’ phenomenon in dynamic cases with complex morphologies, showcasing the method’s adaptive error correction capabilities. This synergistic interaction between RBAR and causality training enables accurate capture of interface evolution, even in challenging scenarios where traditional PINNs fail. Our framework not only resolves the limitations of uniform refinement strategies but also provides a generalizable methodology for solving a broad range of spatio-temporal PDEs. The enhanced performance of the RBAR–causality combined framework demonstrates its strong potential for advancing PINN-based modeling of physical systems characterized by complex, evolving interfaces.

Allen-Cahn equations↗
DOE OSTI2026

Experimental shock Hugoniot and melt boundary temperatures of aluminum

Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.

Hartsfield, Thomas Murray [Sandia National Laborat↗