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.
New insights on carbon black suspension rheology—Anisotropic thixotropy and antithixotropy
We report a detailed experimental study of peculiar thixotropic dynamics of carbon black (CB, Vulcan XC-72) suspensions in mineral oil, specifically the observation of sequential stress increase then decrease at a fixed shear rate in a step-down test. We verify that such dynamics, though peculiar, come from a true material response rather than experimental artifacts. We also reveal how this long-time stress decay is associated with antithixotropic structural change rather than viscoelastic stress relaxation by using orthogonal superposition (OSP) rheometry to probe viscoelastic moduli during the step-down tests. The orthogonal storage and loss moduli are present, showing this two-timescale recovery then decay response, which demonstrates that this response is antithixotropic, and it involves shear-induced structuring. We further show a mechanical anisotropy in the CB suspension under shear using OSP. Based on the rheological results, a microstructural schematic is proposed, considering qualitatively thixotropic structure build-up, antithixotropic densification, and anisotropic structure evolution. Our observation for these CB suspensions is outside the standard paradigm of thixotropic structure-parameter models, and the elastic response provides us with new insight into the transient dynamics of CB suspensions.
Pressure-tuned quantum criticality in the large-D antiferromagnet DTN
Strongly correlated spin systems can be driven to quantum critical points via various routes. In particular, gapped quantum antiferromagnets can undergo phase transitions into a magnetically ordered state with applied pressure or magnetic field, acting as tuning parameters. These transitions are characterized by z = 1 or z = 2 dynamical critical exponents, determined by the linear and quadratic low-energy dispersion of spin excitations, respectively. Employing high-frequency susceptibility and ultrasound techniques, we demonstrate that the tetragonal easy-plane quantum antiferromagnet NiCl 2 · 4SC(NH 2 ) 2 (aka DTN) undergoes a spin-gap closure transition at about 4.2 kbar, resulting in a pressure-induced magnetic ordering. The studies are complemented by high-pressure-electron spin-resonance measurements confirming the proposed scenario. Powder neutron diffraction measurements revealed that no lattice distortion occurs at this pressure and the high spin symmetry is preserved, establishing DTN as a perfect platform to investigate z = 1 quantum critical phenomena. The experimental observations are supported by DMRG calculations, allowing us to quantitatively describe the pressure-driven evolution of critical fields and spin-Hamiltonian parameters in DTN.
Molecular dynamics studies of grain boundary mobility and anisotropy in BCC γ-uranium
Grain morphologies such as grain size and aspect ratio in uranium-based metallic fuels are important microstructural features that can impact various fuel performance properties such as fission-gas-induced swelling, thermal transport, high burnup structure formation, and radiation resistance. Accurate prediction of the fuel grain morphologies requires knowledge of critical grain growth parameters such as grain boundary (GB) mobility and anisotropy. In this work, molecular dynamics simulations were performed to study the GB mobility and its anisotropy in pure body-centered-cubic (BCC) γ uranium. Nine GBs with different combinations of misorientation angles (20°, 30°, 45°) and rotation axes (<100>, <110>, <111>), as well as an additional <111> 38.2° GB were studied using three interatomic potentials. Here, it is found that the GB mobility anisotropy has complex trends, depending on both rotation axis and misorientation. However, in general the <110> rotation axis has the fastest GB mobility at the same misorientation. The results of this work can be used as not only baselines for future studies of GB mobility in uranium-based alloys such as uranium-molybdenum (U-Mo) fuels, but also input for mesoscale modeling of grain growth in uranium-based alloys.
Exponential Decay in the Sensitivity Analysis of Nonlinear Dynamic Programming
In this report, we study the sensitivity of discrete-time dynamic programs with nonlinear dynamics and objective to perturbations in the initial conditions and reference parameters. Under uniform controllability and boundedness assumptions for the problem data, we prove that the directional derivative of the optimal state and control at time $k$, $x^*_k$, and $u^*_k$, with respect to the reference signal at time $i$, $d_i$, will have exponential decay in terms of $|k-i|$ with a decay rate $\rho$ independent of the temporal horizon length. The key technical step is to prove that a version of the convexification approach proposed by Verschueren et al. can be applied to the KKT conditions and results in a convex quadratic program with uniformly bounded data. In turn, Riccati techniques can be further employed to obtain the sensitivity result, borne from the observation that the directional derivatives are solutions of quadratic programs with structure similar to the KKT conditions themselves. We validate our findings with numerical experiments on a small nonlinear, nonconvex, dynamic program.
Characterization of Impedance and Stability for Doubly-Fed Induction Generator Based on Voltage-Modulated Direct Power Control
Voltage-modulated direct power control (VM-DPC) applied to the doubly-fed induction generator (DFIG) offers superior steady-state and transient performance but remains underexplored for suppressing wideband oscillations. Here, this article proposes a comprehensive impedance for the VM-DPC-based DFIG, analyzing its impedance characteristics and stability mechanisms compared with the DFIG based on vector-oriented control (VOC). The unified power transfer function is defined for DFIGs employing VM-DPC or VOC to ensure consistent comparison benchmarks. The comprehensive impedance of VM-DPC-based DFIG, incorporating mechanical and grid-side converter (GSC) dynamics, is derived using complex vector modeling in the αβ -frame. Furthermore, the influence of VM-DPC parameters and grid strength on the stability of grid-connected DFIG systems is assessed through eigenvalue trajectory analysis. Impedance analysis reveals the significant contributions of mechanical and GSC dynamics to DFIG impedance, as well as the narrower frequency range of negative resistance in the VM-DPC-based DFIG compared to the VOC-based DFIG. Stability analysis identifies the VM-DPC parameters of the rotor-side converter as dominant factors affecting system stability and confirms that the VM-DPC-based DFIG achieves better stability under weak grid conditions than its VOC-based counterpart. These findings are validated through simulations and experiments.
Response of the mode Grüneisen parameters with anisotropic compression: A pressure and temperature dependent Raman study of β-Sn
The lattice dynamic response of body-centered tetragonal beta-Sn (I4 1 /amd) under high-pressure and temperature conditions is determined using experimental optical vibration modes. Here, Raman scattering is used to map the phase stability region of beta-Sn to perform mode Gruneisen analysis, and we demonstrate the necessity of an optical intensity calibration for Raman thermometry. The Gruneisen tensor is evaluated along a set of isotherms to address shortcomings of single-mode Gruneisen parameters with respect to anisotropic deformations of this tetragonal structured soft metal. The changes observed here in the Gruneisen tensor as a function of temperature are related to anharmonicity and denote potential criteria for the onset of premelting.
Activated relaxation in supercooled monodisperse atomic and polymeric WCA fluids: Simulation and ECNLE theory
Here, we combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks–Chandler–Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, a c , that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, N c , that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of N c ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.
Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions
In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.
Evaluation of Converter Performance Considering Static and Dynamic Device Part-to-Part Variability
This paper presents a methodology to incorporate and analyze the impact of semiconductor device part-to-part variation on power converter performance. By integrating extensive static and dynamic device characterization data with an automated compact model generation process that reflects manufacturing variability, device models with inherent variability features are utilized in converter simulations for a comprehensive assessment of performance impacts. The traditional converter performance evaluation process typically yields fixed efficiency values, often dismissing the inherent part-to-part variability caused by the manufacturing process of semiconductor devices. To address this limitation, a large population of devices was characterized to capture variations in static parameters-such as transfer, output, and capacitance characteristics-as well as dynamic behaviors, including switching losses. This data-driven approach enables the development of individual compact models, which were then integrated into converter simulations to evaluate efficiency ranges rather than single point estimated values. The converter simulation results show that part-to-part component variation can lead to significant efficiency deviations, exceeding several percentage points in high-power conversion applications. By offering a more accurate representation of converter behavior under real-world manufacturing conditions, this methodology enables designers to anticipate performance variability, improving the robustness of power converter designs.
Computational Tools for Additive Manufacture of Tailored Microstructure and Properties
Additive manufacturing has the potential to revolutionize industrial hardware and unlock efficiency gains through the fabrication of geometries and architectures not possible by conventional processing. Currently most additive builds use a single set of process parameters (e.g. laser power and scan speed) which results in a part with a homogenous microstructure that provides a singular performance level. To move beyond this state, Raytheon Technologies Research Center worked to create a set of computational tools to track material evolution through each step of the additive process. Computational fluid dynamics and phase field models for microstructure evolution as a function of processing parameters, and crystal plasticity models fully coupling microstructure and mechanical properties for performance predictions were leveraged to establish a connection between additive parameters and the final microstructure. This framework was utilized to tailor spatially-varying mechanical properties in a part by appropriately controlling the microstructure evolution during the additive process. Specifically, a turbine blade was 3D printed from nickel superalloy IN718 using laser powder bed fusion with coarse grains in the airfoil section which experiences the highest temperatures and is creep limited while finer grains were printed in the root of the blade which experience higher stresses but at lower temperatures and is therefore fatigue limited. The benefit of being able to intentionally insert coarse grains in the high temperature region of the blade was showcased with a microstructure sensitive creep model that indicates longer creep life for coarser grains.
Coupled Induction Machine and HVAC Models for Simulating HVAC Performance Considering Grid Dynamics in Buildings
This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.
A Flux Detection Probe to Quantify Dynamic Groundwater-Surface Water Exchange in the Hyporheic Zone
A new probe was designed to quantify groundwater-surface water exchange in the hyporheic zone under dynamic stage condition. Current methods focus on either vertical pore water velocity or Darcy flux measurements. Both parameters must be understood to evaluate residence time and mass flux of constituents. Furthermore, most instruments are not well suited for monitoring instantaneous velocity or flux under dynamic exchange conditions. For this reason, the flux detection probe (FDP) was designed that employs electrogeophysical measurements to estimate in situ sediment porosity, which can be used to convert pore water velocity to Darcy flux. Dynamic pore water velocity is obtained by monitoring fluid conductivity and temperature along the FDP probe. Pressure sensors deployed at the top and bottom of the probe provide the additional information necessary to estimate vertical permeability. This study focuses on the use of a geophysical method to estimate pore water velocity, porosity, and permeability within a controlled soil column where simulated river water displaces simulated groundwater. The difference between probe derived and theoretical pore water velocity using natural tracers such as electrical conductivity and temperature was -4.9 and 3.9% for downward flow and 1.1 and 12.8% for upward flow, respectively. The difference in porosity calculated from mass and volume packed in the soil column and probe measure porosity ranged between -3.2% and 1.5%. Also, the calculated hydraulic conductivity differed from probe derived values by -8.9%.
Changes of deformed shell gaps at N ~ 100 in light rare-earth, neutron-rich nuclei
There have been compelling evidences indicating that deformed single-particle states in neutron- rich regions are different from those in the stable region. It was pointed out that for the light rare- earth (60Nd, 62Sm, and 64Gd), neutron-rich (N = 98-102) nuclei, the Woods-Saxon potential, the Nilsson modified oscillator potential with "universal" parameters, and the folded Yukawa potential all failed to describe the correct ordering of neutron single-particle states. The location and size of deformed shell gaps in this mass region are under current debate. We propose a modification for the "standard" Nilsson parameters of Bengtsson and Ragnarsson introduced in 1985. The proposed N-dependent spin-orbit interaction brings an effect of the tensor force for the deformed potential, which causes dynamic changes in deformed shell gaps with neutron number and deformation. By applying the modified Nilsson parameters to generate deformed bases for the Projected Shell Model, we demonstrate that the calculation can explain consistently the current experimental data, including the ground state configuration in odd-neutron nuclei, upbending of the yrast moment of inertia at higher spins and the energies of 2-quasineutron 6- and 4- isomers in even-even nuclei.
Dynamic Carrier Modulation via Nonlinear Acoustoelectric Transport in van der Waals Heterostructures
Dynamically manipulating carriers in van der Waals heterostructures could enable solid-state quantum simulators with tunable lattice parameters. A key requirement is the formation of deep potential wells to reliably trap excitations. Here, we report the observation of nonlinear acoustoelectric transport and dynamic carrier modulation in boron nitride-encapsulated graphene devices coupled to intense surface acoustic waves (SAWs) on LiNbO 3 substrates. SAWs generate strong acoustoelectric current densities ( J AE ), transitioning from linear to nonlinear regimes with increasing SAW intensity. In the nonlinear regime, periodic carrier (electrons, holes, or their mixtures) stripes emerge. Using counter-propagating SAWs, we create standing SAWs (SSAWs) to dynamically manipulate charge distributions without static gates. The saturation of J AE , attenuation transitions, and tunable resistance peaks confirms strong carrier localization. Finally, these results establish SAWs as a powerful tool for controlling carrier dynamics in two-dimensional (2D) materials, paving the way for the development of time-dependent quantum systems and acoustic lattices for quantum simulation.
A sui generis whipping-instability-based self-sequencing multi-monodisperse 2D spray from an anisotropic microfluidic liquid jet device
Well-defined aerosols pave the way for versatile basic and applied research. Here, we demonstrate a unique whipping instability that generates from a high-aspect-ratio microfluidic device resulting in a unique steady-state gas-focused whipping jet (WJ) without any need for electrification. This WJ device emanates a multi-monodisperse whipping spray jet with a two-dimensional (2D) profile. We demonstrate this phenomenon based on various fluidic parameters theoretically and experimentally. The 2D WJ’s unique behavior is derived using analytical fluid dynamics to explain jet diameter, whipping regime, and spreading angle. The phenomenon is further characterized experimentally by measuring the angle with respect to the flow rate, the distances between droplets, the droplet shapes, and the reproducibility of these parameters. We also explain the precise fabrication of such inexpensive devices. Lastly, we highlight these devices’ potential use as sample environments in versatile applications ranging from cryoelectron microscopy over mass spectrometry to drug formulation and structural studies at X-ray free-electron lasers.
Neural Lyapunov Differentiable Predictive Control
We present a learning-based predictive control methodology using the differentiable programming framework with sufficient Lyapunov-based guarantees. The neural Lyapunov differentiable predictive control (NLDPC) learns the policy by constructing a computational graph encompassing the system dynamics, state, and input constraints, and the necessary Lyapunov certification constraints, and thereafter using the automatic differentiation to update the neural policy parameters. In conjunction, our approach jointly learns a Lyapunov function that certifies the regions of state-space with stable dynamics. We also provide an initial condition sampling-based statistical guarantee for the training of NLDPC. Our offline training approach provides a computationally efficient and scalable alternative to classical model predictive control solutions, which now directly comes with Lyapunov based guarantees as shown in this paper. We substantiate the advantages of the proposed approach with simulations to stabilize the double integrator model and on an example of controlling an aircraft model.
Ion density waves driving the formation of filamentary dust structures
The PlasmaKristall-4 (PK-4) experiment on the International Space Station allows for the study of the 3-dimensional interaction between plasma and dust particles. Previous simulations of the PK-4 environment have discovered fast moving ionization waves in the dc discharge. These ionization waves vary the plasma parameters by up to an order of magnitude, which may affect the mechanisms responsible for the self-organization of chains seen in the PK-4 experiment. Here, we adapt a molecular dynamics simulation to employ temporally varying plasma conditions in order to investigate the effect on the dust charging and electrostatic potential. In order to describe the differences between the average of the plasma conditions and the time-varying plasma condition, we present a model to reproduce the interaction that takes into account the negative potential from the dust grain and the positive potential from the ion wake.