Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic parameter”

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 631 records · Page 35

Dynamic Stabilization of Metastable States in Triple-Well Ferroelectric Sn 2 P 2 S 6

Polarization dynamics in ferroelectric materials is governed by the effective potential energy landscape of the order parameter. The unique aspect of ferroelectrics compared to many other transitions is the possibility of more than two potential wells, leading to complicated energy landscapes with new fundamental and functional properties. In this work, direct dynamic evidence is revealed of a triple-well potential in the metal thiophosphate Sn 2 P 2 S 6 compound using multivariate scanning probe microscopy combined with theoretical simulations. The key finding is that the metastable zero polarization state can be accessed through a gradual switching process and is stabilized over a broad range of electric fields. Simulations confirm that the observed zero polarization state originates from a kinetic stabilization of the nonpolar state of the triple-well, as opposed to domain walls. Dynamically, the triple-well of Sn 2 P 2 S 6 becomes equivalent to antiferroelectric hysteresis loops. Therefore, this material combines the robust and well-defined domain structure of a proper ferroelectric with dynamic hysteresis loops present in antiferroelectrics. Moreover, the triple-well enhances mem-capacitive effects in Sn 2 P 2 S 6 , which are forbidden for ideal double-well ferroelectrics. These findings provide a path to tunable electronic elements for beyond binary high-density computing devices and neuromorphic circuits based on dynamic properties of the triple-well.

36 MATERIALS SCIENCE↗

Ab initio generalized Langevin equation

We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

97 MATHEMATICS AND COMPUTING↗

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗

A physics-constrained neural network for multiphase flows

The present study develops a physics-constrained neural network (PCNN) to predict sequential patterns and motions of multiphase flows (MPFs), which includes strong interactions among various fluid phases. To predict the order parameters, which locate individual phases in the future time, a neural network (NN) is applied to quickly infer the dynamics of the phases by encoding observations. The multiphase consistent and conservative boundedness mapping algorithm (MCBOM) is next implemented to correct the predicted order parameters. This enforces the predicted order parameters to strictly satisfy the mass conservation, the summation of the volume fractions of the phases to be unity, the consistency of reduction, and the boundedness of the order parameters. Then, the density of the fluid mixture is updated from the corrected order parameters. Finally, the velocity in the future time is predicted by another NN with the same network structure, but the conservation of momentum is included in the loss function to shrink the parameter space. The proposed PCNN for MPFs sequentially performs (NN)-(MCBOM)-(NN), which avoids nonphysical behaviors of the order parameters, accelerates the convergence, and requires fewer data to make predictions. Numerical experiments demonstrate that the proposed PCNN is capable of predicting MPFs effectively.

Mechanics↗

Energy efficiency in industrial drying: A hybrid ultrasonic system with a novel dynamic optimization framework

Drying processes are among the most energy-consuming operations in industrial and manufacturing settings, demanding strategic selection, design, and control for enhanced efficiency. Advancing drying technologies is critical for improving sustainability, lowering energy use, reducing carbon emissions, and minimizing waste. This study explores two innovative strategies aimed at transforming drying processes into sustainable, low-carbon systems by reducing energy consumption, minimizing waste, and maintaining a strong emphasis on preserving product quality. The first strategy showcases a sub-pilot scale hybrid ultrasonic-convective dryer for agrifood products. This technology, powered by electricity (process electrification), integrates non-thermal ultrasonic dehydration with convective heating and is presented as a sustainable and energy-efficient solution that enhances eco-friendly practices. The second strategy involves introducing and implementing a novel, multiobjective, mixed integer dynamic optimization technique to determine the optimal time-dependent process parameter values for the drying operation. This optimization technique yields operating conditions that are piecewise constant in time aiming to maximize the energy efficiency of the hybrid ultrasonic-convective dryer while ensuring strict adherence to product quality constraints. By adopting the hybrid ultrasonic-convective dryer, a notable 35% improvement in energy efficiency was achieved compared to conventional hot-air drying systems for drying apple slices. The proposed optimization framework further enhanced energy efficiency by nearly 14% over the most efficient process on the identical testbed, under static operating conditions. The reported enhancements have been experimentally validated. Regarding drying time (thereby improving production yield), the developed hybrid ultrasonic-convective dryer demonstrates as much as a 41% reduction in total processing time, which is further optimized by an additional 10% using our proposed optimization framework. The research outcomes have profound implications for the design and operation of drying systems, encompassing crucial aspects such as process electrification, cost-effectiveness, energy savings, time efficiency, product yield, product quality, and process automation.

Dynamic optimization↗

Characterizing skyrmion flow phases with principal component analysis

Principal component analysis (PCA) is a powerful method that can identify patterns in large, complex data sets by constructing low-dimensional order parameters from higher-dimensional feature vectors. There are increasing efforts to use space-and-time-dependent PCA to detect transitions in nonequilibrium systems that are difficult to characterize with equilibrium methods. Here, we demonstrate that feature vectors incorporating the position and velocity information of driven skyrmions moving through random disorder permit PCA to resolve different types of disordered skyrmion motion as a function of driving force and the ratio of the Magnus force to the dissipation. Since the Magnus force creates gyroscopic motion and a finite Hall angle, skyrmions can exhibit a greater range of flow phases than what is observed in overdamped driven systems with quenched disorder. We show that in addition to identifying previously known skyrmion flow phases, PCA detects several additional phases, including different types of channel flow, moving fluids, and partially ordered states. Guided by the PCA analysis, we further characterize the disordered flow phases to elucidate the different microscopic dynamics and show that the changes in the PCA-derived order parameters can be connected to features in bulk transport measures, including the transverse and longitudinal velocity-force curves, differential conductivity, topological defect density, and changes in the skyrmion Hall angle as a function of drive. We discuss how asymmetric feature vectors can be used to improve the resolution of the PCA analysis, and how this technique can be extended to find disordered phases in other nonequilibrium systems with time-dependent dynamics.

36 MATERIALS SCIENCE↗

Dynamic phase transformations in additively manufactured Ti-6Al-4V during thermo-mechanical gyrations

A complex interaction of process parameters, geometry and scan strategies in Additive Manufacturing (AM), can bring about spatial and temporal transients, i.e., Σ T ( x, y, z, time ), within a part. Published literature focusses on fluctuating thermal cycles on the microstructure evolution. However, the microstructural variations have not been correlated to dynamic flow behavior due to the macro- and micro-scale phenomena, i.e., accumulated plastic strains brought about by large thermal gradients, transformational strains and crystallographic misfit strains. Therefore, here we studied the mechanical response of Ti6Al4V alloys produced by AM under externally imposed controlled thermo-mechanical reversals in a Gleeble® thermo-mechanical simulator. The stress-strain behaviors were correlated to phase fractions, lattice strains, and also limited information on crystallographic texture using neutron diffraction techniques at the VULCAN Beamline at SNS, ORNL and also metallographic studies. The results are discussed and rationalized based on theories of static and dynamic phase transformations.

36 MATERIALS SCIENCE↗

A coarse-grained simulation model for colloidal self-assembly via explicit mobile binders

Colloidal particles with mobile binding molecules constitute a powerful platform for probing the physics of self-assembly. Binding molecules are free to diffuse and rearrange on the surface, giving rise to spontaneous control over the number of droplet–droplet bonds, i.e., valence, as a function of the concentration of binders. This type of valence control has been realized experimentally by tuning the interaction strength between DNA-coated emulsion droplets. Optimizing for valence two yields droplet polymer chains, termed ‘colloidomers’, which have recently been used to probe the physics of folding. To understand the underlying self-assembly mechanisms, here we present a coarse-grained molecular dynamics (CGMD) model to study the self-assembly of this class of systems using explicit representations of mobile binding sites. Further, we explore how valence of assembled structures can be tuned through kinetic control in the strong binding limit. More specifically, we optimize experimental control parameters to obtain the highest yield of long linear colloidomer chains. Subsequently tuning the dynamics of binding and unbinding via a temperature-dependent model allows us to observe a heptamer chain collapse into all possible rigid structures, in good agreement with recent folding experiments. Our CGMD platform and dynamic bonding model (implemented as an open-source custom plugin to HOOMD-Blue) reveal the molecular features governing the binding patch size and valence control, and opens the study of pathways in colloidomer folding. This model can therefore guide programmable design in experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

DRIPS: A framework for dimension reduction and interpolation in parameter space

Reduced-order models are often used to describe the behavior of complex systems, whose simulation with a full model is too expensive, or to extract salient features from the full model’s output. We introduce a new model-reduction framework DRIPS (dimension reduction and interpolation in parameter space) that combines the offline local model reduction with the online parameter interpolation of reduced-order bases (ROBs). The offline step of this framework relies on dynamic mode decomposition (DMD) to build a low-rank linear surrogate model, equipped with a local ROB, for quantities of interest derived from the training data generated by repeatedly solving the (nonlinear) high-fidelity model for multiple parameter points. The online step consists of the construction of a parametric reduced-order model for each target/test point in the parameter space, with the interpolation of ROBs done on a Grassman manifold and the interpolation of reduced-order operators done on a matrix manifold. The DMD component enables DRIPS to model (typically low-dimensional) quantities of interest directly, without having to access the (typically high-dimensional and possibly nonlinear) operators in a high-fidelity model that governs the dynamics of the underlying high-dimensional state variables, as required in projection-based reduced-order modeling. A series of numerical experiments suggests that DRIPS yields a model reduction, which is computationally more efficient than the commonly used projection-based proper orthogonal decomposition; it does so without requiring a prior knowledge of the governing equation for quantities of interest. Furthermore, for the nonlinear systems considered, DRIPS is more accurate than Gaussian-process interpolation (Kriging).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular dynamics simulations of fission gas xenon (Xe) diffusion at UO 2 grain-boundaries (Rev.1)

The diffusivity of fission gas xenon (Xe) at UO 2 grain-boundaries is one of the most important parameters in mechanistic modeling of fission gas diffusion in UO 2 based nuclear fuels. In this report, we use molecular dynamics simulations to investigate the Xe diffusivity in UO 2 grain-boundaries, employing the many-body potential developed by Cooper, Rushton and Grimes for UO 2 . Three different types of grain-boundaries are investigated, twist Σ5, tilt Σ5, and a random grain-boundary. Diffusion activation energies in the range of 0.39 – 1.46 eV are obtained for the Xe diffusivity. Comparison to results for the uranium vacancy diffusivity from MD simulations employing the same methodology suggests a weak to moderate attractive Xe-uranium vacancy binding energy depending on the grain-boundary type.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Operational resilience of additively manufactured parts to stealthy cyberphysical attacks using geometric and process digital twins

Cyberphysical attacks on the digital backbone of Additive Manufacturing (AM) can compromise the printed part’s functionality. They can alter features in the digital geometry to introduce geometric defects (e.g., missing fillets) or alter process parameters to create local defects (e.g., voids). Addressing the downtime, waste, and quality deterioration associated with existing solutions requires operational resilience, i.e., rapid elimination or disruption of defect formation (to retain part function) without production stoppage or part disposal (to retain yield). This need is unmet due to the inherently unpredictable nature of attack-induced alterations, lack of access to the original geometric model for identification of altered geometric features, and in-process imposition of unknown process dynamics via attack-driven alteration of real-time-uncontrolled (or exogenous) parameters. This work establishes the above-mentioned operational resilience for the first time by creating two Digital Twins (DT). The Geometric DT (Geo-DT) is based on a unique physical-field-driven soft sensor and topology optimization method. The Process Digital Twin (Pro-DT) combines local defect quantification with a novel Reinforcement Learning formulation and training method. The importance of these methodological advances and the scalability of our approach are examined on a real AM testbed. It is shown that Geo-DT can correct geometric defects without access to the original digital geometry or explicit knowledge of attack-altered geometric features. Further, Pro-DT can accelerate real-time disruption of local defects despite attack-driven imposition of unknown process dynamics. We discuss how our framework goes beyond the contemporary focus on pre-attack security and in-attack detection towards resilience for AM and beyond.

Additive Manufacturing↗

Dynamic deuterium retention properties of tungsten measured using laser-induced breakdown spectroscopy

Dynamic deuterium (D) retention properties of tungsten (W) are investigated under steady-state D plasma exposure in the PISCES-A linear plasma device. In contrast to static retention, dynamic retention is quickly released from the material at the termination of the incident plasma flux, and thus in-operando laser-induced breakdown spectroscopy (LIBS) measurements have been conducted during steady-state D plasma exposure. A procedure is, first, established to extract the dynamic retention component from the in-operando LIBS D I 656.1 nm line intensity, which can contain contributions from dynamic and static retention as well as background D/D 2 gas, excited by both steady-state and laser-induced plasmas. Using the developed procedure, the dynamic D retention in W is systematically examined while scanning the following plasma exposure parameters: incident ion energy, E i , sample temperature, T s , incident ion flux, Γ i . No clear E i dependence is seen in the range of E i ~ 45–175 eV, as expected from a small variation in the average implantation depth (~3-6 nm) of D in W. The dynamic retention is found to monotonously decrease with increasing T s from 348 to 573 K, while the total static retention is reported to peak at T s ~ 500–600 K. It is revealed that the dynamic retention linearly increases with increasing Γ i , and then saturates at Γ i ≥ 0.75 × 10 21 m –2 s –1 . Possible physical mechanisms for the observed dependence of the dynamic retention on T s and Γ i are discussed. Furthermore, sequential pure He, followed by pure D, plasma exposures show that the dynamic D retention is not strongly affected by He bubbles in the near-surface region.

super-saturated surface layer↗

GridLAB-D Technical Support Document: Residential Equivalent Thermal Parameter Model

GridLAB-D is a power distribution systems simulation and analysis tool developed by Pacific Northwest National Laboratory. Although GridLAB-D and its underlying residential load model has been used in a variety of power systems analyses, comprehensive documentation of the derivation of the residential house model has not been formally published. The purpose of this technical support document is to serve as a comprehensive resource for the thermal dynamics modeling and house definition implemented in GridLAB-D. It covers model derivation, parameter definitions, and implementation steps in the codebase.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mathematical analysis of laboratory microbial experiments demonstrating deterministic chaotic dynamics

Presented is a system of four ordinary differential equations and a mathematical analysis of microbiological experiments in a four-component chemostat—nutrient n, rods r, cocci c, and predators p. The analysis is consistent with the conclusion that previous experiments produced features of deterministic chaotic and classical dynamics depending on dilution rate. The surrogate model incorporates as much experimental detail as possible, but necessarily contains unmeasured parameters. The objective is to understand better the differences between model simulations and experimental results in complex microbial populations. The key methodology for simulation of chaotic dynamics, consistent with the measured dilution rate and microbial volume averages, was to cause the preference of p for r vs. c to vary with the r and c concentrations, to make r more competitive for nutrient than c, and to recycle some dying p biomass, leading to a modified version of the Monod kinetics model. Our mathematical model demonstrated that the occurrence of chaotic dynamics requires a predator, p, preference for r versus c to increase significantly with increases in r and c populations. Also included is a discussion of several generalizations of the existing model and a possible involvement of the minimum energy dissipation principle. This principle appears fundamental to thermodynamic systems including living systems. Several new experiments are suggested.

59 BASIC BIOLOGICAL SCIENCES↗

On the use of vibronic coherence to identify reaction coordinates for ultrafast excited-state dynamics of transition metal-based chromophores

In this study, the question of whether one can use information from quantum coherence as a means of identifying vibrational degrees of freedom that are active along an excited-state reaction coordinate is discussed. Specifically, we are exploring the notion of whether quantum oscillations observed in single-wavelength kinetics data exhibiting coherence dephasing times that are intermediate between that expected for either pure electronic or pure vibrational dephasing are vibronic in nature and therefore may be coupled to electronic state-to-state evolution. In the case of a previously published Fe(ιι) polypyridyl complex, coherences observed subsequent to 1 A 1 → 1 MLCT excitation were linked to large-amplitude motion of a portion of the ligand framework; dephasing times on the order of 200–300 fs suggested that these degrees of freedom could be associated with ultrafast (~100 fs) conversion from the initially formed MLCT excited state to lower-energy, metal-centered ligand-field excited state(s) of the compound. Incorporation of an electronically benign but sterically restrictive Cu(ι) ion into the superstructure designed to interfere with this motion yielded a compound exhibiting a ~25-fold increase in the compound’s MLCT lifetime, a result that was interpreted as confirmation of the initial hypothesis. However, new data acquired on a different chemical system – Cr(acac') 3 (where acac' represents various derivatives of acetylacetonate) – yielded results that call into question this same hypothesis. Coherences observed subsequent to 4 A 2 → 4 T 2 ligand-field excitation on a series of molecules implicated similar vibrational degrees of freedom across the series, but exhibited dephasing times ranging from 340 fs to 2.5 ps without any clear correlation to the dynamics of excited-state evolution in the system. Taken together, the results obtained on both of these chemical platforms suggest that while identification of coherences can indeed point to degrees of freedom that should be considered as candidate modes for defining reaction trajectories, our understanding of the factors that determine the interplay across coherences, dephasing times, and electronic and geometric structure is insufficient at the present time to view this parameter as a robust metric for differentiating active versus spectator modes for ultrafast dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of McMillan mappings III. Symmetric map with mixed nonlinearity

This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.

43 PARTICLE ACCELERATORS↗

In-situ sensor monitoring of multi-class gas porosity formation in laser powder bed fusion using convolutional neural network

In-situ monitoring of defect formation remains a significant challenge in the laser powder bed fusion (LPBF) process. Recent advances have enabled real-time defect detection with machine learning and in-situ sensing technologies; however, most studies focus on binary classification of keyhole pores, limiting nuanced multi-class pore differentiation and formation mechanisms. This work introduces a multi-class pore detection framework (no pore, small pores < 15 µm, and large pores > 15 µm) by leveraging photodiode sensor data alongside high-fidelity synchrotron X-ray imaging. The 15 µm threshold is selected to distinguish between two fundamentally different defect mechanisms, following the physical size-mechanism boundary established by prior high-resolution synchrotron X-ray characterization of Al6061 LPBF. Distinguishing these classes is critical because large keyhole pores are structurally detrimental, whereas small gas pores are often benign, requiring different process control strategies. Thermal emission monitoring data collected simultaneously with high-speed X-ray imaging at the Stanford Synchrotron Radiation Lightsource (SSRL), are correlated with subsurface melt pool dynamics to establish ground truth. Continuous Wavelet Transform (CWT) with optimized parameters converts the photodiode time-series signals into time–frequency images, facilitating feature extraction. Convolutional Neural Networks (CNN) are then applied for real-time multi-class pore classification in an average inference time of 1 ms per signal window. It achieves 79% accuracy and an Area Under the Receiver Operating Characteristic curve (AUC ROC) score of 0.89 with five-fold cross-validation. The results demonstrate that coupling CWT-based feature engineering with CNN architecture enables reliable multi-class pore detection in Al6061 builds using affordable in-situ sensors. This approach advances scalable and affordable quality assurance in additive manufacturing by moving beyond binary defect detection toward more nuanced classification of porosity mechanisms with in-situ sensors and machine learning.

Laser powder bed fusion, Multi-class pores, In-sit↗