Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative systems”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Autonomous Emergency Landing for Fixed-Wing Aircraft with Energy-Constrained Closed-Loop Prediction

Here this paper presents a new approach for autonomous motion planning for aircraft suffering from a loss-of-thrust emergency. Specifically, we show how modifications to the Closed-Loop Rapidly exploring Random Trees (CL-RRT) framework combined with controlled energy dissipation can enable rapid and effective kinodynamic motion planning. This CL-RRT Glide algorithm uses closed-loop prediction not only for node connections but also to estimate the remaining energy and prune infeasible paths. This greatly speeds up the search process, which is essential for emergency situations. In addition, we improve the ability of the gliding aircraft to reach a goal position and energy state. We do so by creating a Dissipative Total Energy Control Scheme (TECS). Dissipative TECS enables the glider to lose excess altitude in order to reach a desired energy level. Simulation results illustrate how the proposed methods enable faster motion planning. We also integrate the system into a small unmanned aerial vehicle system and experimentally demonstrate autonomous glide planning and execution during a motor-failure event. This type of algorithm can primarily benefit unmanned aircraft but can also serve to assist pilots in stressful emergency situations.

42 ENGINEERING↗
DOE OSTI2023

MD Simulation of Water Using a Rigid Body Description Requires a Small Time Step to Ensure Equipartition

In simulations of aqueous systems, it is common to freeze the bond vibration and angle bending modes in water to allow for a longer time step δt for integrating the equations of motion. Thus, δt = 2 fs is often used in simulating rigid models of water. We simulate the SPC/E model of water using δt from 0.5 to 3.0 fs and up to 4 fs using hydrogen mass repartitioning. In these simulations, we find that for all but δt = 0.5 fs, equipartition is not obtained between translational and rotational modes, with the rotational modes exhibiting a lower temperature than the translation modes. To probe the reasons for the lack of equipartition, we study the autocorrelation of the translational velocity of the center of mass and the angular velocity of the rigid water molecule, respectively. Here, we find that the rotational relaxation occurs on a timescale comparable to vibrational periods, calling into question the original motivations for freezing the vibrations. Furthermore, a time step with δt ≥ 1 fs is not able to capture accurately the fast rotational relaxation, which reveals its impact as an effective slowing-down of rotational relaxation. The fluctuation–dissipation relation then leads to the conclusion that the rotational temperature should be cooler for δt greater than the reference value of 0.5 fs. Consideration of fluctuation–dissipation in equilibrium molecular dynamics simulations also emphasizes the need to capture the temporal evolution of fluctuations with fidelity and the role of δt in this regard. The time step also influences the solution thermodynamic properties: both the mean system potential energies and the excess entropy of hydration of a soft repulsive cavity are sensitive to δt.

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

A Twin Circuit Theory-Based Framework for Oscillation Event Analysis in Inverter-Dominated Power Systems With Case Study for Kaua‘i System

Here, this paper proposes a real-world oscillation event analysis framework for power systems that include inverter-based resources together with synchronous generators. Specifically, the proposed framework combines both measurement-and model-based techniques to readily identify potential oscillation sources, replay the oscillation event with numerical simulation, unveil the underlying oscillation mechanism, and suggest mitigation methods for a wide range of oscillation events. To strengthen the theoretical foundation of our analysis framework, this paper proposes a twin circuit theory that provides theoretical support for one key utilized but not well-proven measurement-based oscillation source identification method-Dissipating Energy Flow. Our twin circuit theory also shows that adopting well-tuned grid-forming inverters can be a potential mitigation method for oscillation events. Finally, the effectiveness of our proposed oscillation event analysis framework is demonstrated by addressing a real-world 18-20 Hz oscillation event in Kaua‘i's power system on November 21, 2021.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Open quantum systems and Schwinger-Keldysh holograms

We initiate the study of open quantum field theories using holographic methods. Specifically, we consider a quantum field theory (the system) coupled to a holographic field theory at finite temperature (the environment). We investigate the effects of integrating out the holographic environment with an aim of obtaining an effective dynamics for the resulting open quantum field theory. The influence functionals which enter this open effective action are determined by the real-time (Schwinger-Keldysh) correlation functions of the holographic thermal environment. To evaluate the latter, we exploit recent developments, wherein the semiclassical gravitational Schwinger-Keldysh saddle geometries were identified as complexified black hole spacetimes. We compute real-time correlation functions using holographic methods in these geometries, and argue that they lead to a sensible open effective quantum dynamics for the system in question, a question that hitherto had been left unanswered. In addition to shedding light on open quantum systems coupled to strongly correlated thermal environments, our results also provide a principled computation of Schwinger-Keldysh observables in gravity and holography. In particular, these influence functionals we compute capture both the dissipative physics of black hole quasi- normal modes, as well as that of the fluctuations encoded in outgoing Hawking quanta, and interactions between them. We obtain results for these observables at leading order in a low frequency and momentum expansion in general dimensions, in addition to determining explicit results for two dimensional holographic CFT environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Physics successfully implements Lagrange multiplier optimization

Optimization is a major part of human effort. While being mathematical, optimization is also built into physics. For example, physics has the Principle of Least Action; the Principle of Minimum Power Dissipation, also called Minimum Entropy Generation; and the Variational Principle. Physics also has Physical Annealing, which, of course, preceded computational Simulated Annealing. Physics has the Adiabatic Principle, which, in its quantum form, is called Quantum Annealing. Thus, physical machines can solve the mathematical problem of optimization, including constraints. Binary constraints can be built into the physical optimization. In that case, the machines are digital in the same sense that a flip–flop is digital. A wide variety of machines have had recent success at optimizing the Ising magnetic energy. We demonstrate in this paper that almost all those machines perform optimization according to the Principle of Minimum Power Dissipation as put forth by Onsager. Further, we show that this optimization is in fact equivalent to Lagrange multiplier optimization for constrained problems. We find that the physical gain coefficients that drive those systems actually play the role of the corresponding Lagrange multipliers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Universal time scalings of sensitivity in Markovian quantum metrology

Assuming Markovian time evolution of a quantum sensing system, we study the general characterization of the optimal sensitivity scalings with time, under most general quantum control protocols. We allow the estimated parameter to influence both the Hamiltonian as well as the dissipative part of the quantum master equation and focus on the asymptotic-time along with the short-time sensitivity scalings. We find that via simple algebraic conditions (in terms of the Hamiltonian, the jump operators as well as their parameter derivatives), one can characterize the four classes of metrological models that represent: quadratic-linear, quadratic-quadratic, linear-linear, and linear-quadratic time scalings. We also investigate the relevant time scales on which the transition between the two regimes appears. Additionally, we provide universal numerical methods to obtain quantitative bounds on sensitivity that are the tightest that exist in the literature. Simplicity and universality of our results make it suitable for diverse applications in quantum metrology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Shape optimization and fluid control of near isothermal compressor for transcritical carbon dioxide cycle

The compressor holds great promise for enhanced energy efficiency over a wide range of vapor compression systems (VCS). However, the majority of them focused on liquid refrigerant injection or inter-stage cooling technologies to improve heat transfer. Increasing the heat transfer area provides another perspective to promote heat dissipation, which commonly involves shape optimization whilst bringing some limitations, mainly from the manufacturing side. To fill the research gap, this paper compares the performance of isothermal compressors under different compression chamber shapes and chamber connection configurations. Additionally, the fluid control of the liquid piston was combined with shape optimization to explore the potential for further heat transfer enhancement. The simulation results have demonstrated that the COP improvements can reach up to 9% compared with the baseline case by optimizing the chamber shape and properly controlling the liquid piston speed, which provides a guideline for designing and customizing the isothermal compression chamber.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

LES simulations of a vacuum membrane distillation channel with geometric alterations

3D LES simulations were carried out to study the performance of a vacuum membrane distillation module. A wiggly wall profile with/without embedded stiffeners was considered to alleviate axial and radial temperature polarization, the cause of performance decrease in membrane distillation. Results of the flow field show the wiggles create unsteady vortex shedding inducing intense mixing in the channel. Vortex shedding intensity increases as the Reynolds number increases or stiffeners are added. 3D results show a 56% improvement in flux when moving from a flat sheet membrane to a wiggly membrane with stiffeners at a constant Reynolds number, corresponding to an improvement from 11.7 to 41.7 in the Nusselt number, showing that the alterations improved the flux performance by enhancing the heat transfer along the membrane surface and therefore alleviating temperature polarization, a critical bottleneck in membrane distillation systems. Further, a merit criterion was defined based on the Nusselt number and friction factor, and a 40% increase in merit was shown switching from a flat to a wiggly channel, while a 97% merit increase was seen going from flat to wiggly with stiffeners. A 42% enhancement in the flux was also seen moving from a straight channel to a wiggly channel at a higher Reynolds number which highlights the importance of strategically choosing the mass flow rate, as well as inducing flow separation and vortex shedding in the channel to promote mixing and dissipate the thermal boundary layer. The variation in results between the wiggly module with stiffeners and the wiggly module at high Reynolds numbers suggests enhancing mixing structures can be more impactful on flux for this geometry than increasing the flow rate to a turbulent/transitional regime. However, both are preferable for peak system flux performance. Furthermore, a 2D approximation was used to perform simulations on more extended channels to examine the length degradation. The modules with wiggly channels performed at the same flux level with a doubling of the length. On the other hand, the flat sheet modules experienced length degradation by temperature polarization and dropped in flux yield by around 13% of the short-channel value. This work illustrates that modeling the system and understanding how the performance decreases as the membrane surface area increases are critical for a larger module (scaling up from a lab to a prototype module) to maintain high flux performance.

42 ENGINEERING↗
DOE OSTI2026

Optical diagnostic design for measuring the radiation front with a mid-leg pumped divertor on DIII-D

The design of an optical diagnostic system to localize the radiation front in a mid-leg pumped divertor configuration on DIII-D is presented. Divertor detachment is a key mechanism for handling power exhaust in tokamaks, and mid-leg pumping offers a promising approach to radiatively dissipate power while maintaining acceptable core performance. To predict the location of the radiation front and design a spectroscopic diagnostic to measure its position, a database of SOLPS-ITER simulations across a range of input powers and gas puffing rates representative of DIII-D operation was employed. These simulations provide self-consistent plasma backgrounds for a Cherab-Raysect synthetic diagnostic framework, which incorporates detailed tokamak geometry and physically accurate ray-tracing. Synthetic simulations of optical sightlines and viewing cone geometry were used to calculate line-integrated emission from the plasma, particularly of the C III 465 nm line, to serve as a proxy for the T e ≈ 7 − 10 eV temperature region associated with the onset and evolution of divertor detachment. Neutral deuterium emission is also evaluated for comparison. The synthetic diagnostic is used to assess the performance of existing DIII-D optical systems, including filterscopes and the Multichord Divertor Spectrometer, and to optimize line-of-sight placement within mechanical and installation constraints. The results provide quantitative guidance for diagnostic implementation on DIII-D and demonstrate the advantage of integrated synthetic diagnostics for divertor design studies and future advanced divertor concepts.

Cherab↗
DOE OSTI2023

Electron-scale current sheets and energy dissipation in 3D kinetic-scale plasma turbulence with low electron beta

ABSTRACT Three-dimensional kinetic-scale turbulence is studied numerically in the regime where electrons are strongly magnetized (the ratio of plasma species pressure to magnetic pressure is βe = 0.1 for electrons and βi = 1 for ions). Such a regime is relevant in the vicinity of the solar corona, the Earth’s magnetosheath, and other astrophysical systems. The simulations, performed using the fluid-kinetic spectral plasma solver (sps) code, demonstrate that the turbulent cascade in such regimes can reach scales smaller than the electron inertial scale, and results in the formation of electron-scale current sheets (ESCS). Statistical analysis of the geometrical properties of the detected ESCS is performed using an algorithm based on the medial axis transform. A typical half-thickness of the current sheets is found to be on the order of electron inertial length or below, while their half-length falls between the electron and ion inertial length. The pressure–strain interaction, used as a measure of energy dissipation, exhibits high intermittency, with the majority of the total energy exchange occurring in current structures occupying approximately 20 per cent of the total volume. Some of the current sheets corresponding to the largest pressure–strain interaction are found to be associated with Alfvénic electron jets and magnetic configurations typical of reconnection. These reconnection candidates represent about 1 per cent of all the current sheets identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Nanoelectromechanical Control of Spin–Photon Interfaces in a Hybrid Quantum System on Chip

Color centers (CCs) in nanostructured diamond are promising for optically linked quantum technologies. Scaling to useful applications motivates architectures meeting the following criteria: C1 individual optical addressing of spin qubits; C2 frequency tuning of spin-dependent optical transitions; C3 coherent spin control; C4 active photon routing; C5 scalable manufacturability; and C6 low on-chip power dissipation for cryogenic operations. Here, we introduce an architecture that simultaneously achieves C1–C6. We realize piezoelectric strain control of diamond waveguide-coupled tin vacancy centers with ultralow power dissipation necessary. The DC response of our device allows emitter transition tuning by over 20 GHz, combined with low-power AC control. We show acoustic spin resonance of integrated tin vacancy spins and estimate single-phonon coupling rates over 1 kHz in the resolved sideband regime. Combined with high-speed optical routing, our work opens a path to scalable single-qubit control with optically mediated entangling gates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Computing equilibrium free energies through a nonequilibrium quench

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium “quench” process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications.

Chemistry↗
DOE OSTI2024

A Study of Intermittent Turbulence in Stable Arctic Boundary Layers

Abstract Polar boundary layers are difficult to model due to the existence of intermittent turbulence within stable layers. Here we present a case study evaluation of coherent structures in a stable boundary layer observed during a series of flights with an uncrewed aircraft system (DataHawk2) on 19 October 2016 at Oliktok Point, Alaska as part of the ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems) field campaign. During a sequence of five flights over a nine-hour period, 57 profiles of atmospheric properties (0–400 m a.g.l) were collected. Turbulence was identified using derived Richardson Number, temperature structure function parameter, and turbulence kinetic energy dissipation. Throughout all flights on this strongly stable day, intermittent turbulence was observed. These turbulent layers showed well-mixed potential temperature profiles embedded within otherwise stable potential temperature profiles; often resulting in a characteristic staircase pattern. Turbulent layers ranged from 1 to 30 m deep, with most individual layers being 1–2 m deep. Vertical propagation velocities of layers in the lower atmosphere were on the order of a few cm s −1 , typical of non-convective environments. In different regions of the profile, turbulence was driven by a different balance of buoyancy and shear forces, with turbulence in the near surface environment driven by strong shear forces overcoming strong resistance to buoyancy, while turbulence in elevated layers characterized by weaker shear forces overcoming weaker resistance to buoyancy. We discuss the potential of such datasets for improving subgrid parameterizations of small-scale turbulence embedded within stable boundary layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Substituent Effects on Transient, Carbodiimide-Induced Geometry Changes in Diphenic Acids

Nucleotide-induced conformational changes in motor proteins are key to many important cell functions. Inspired by this biological behavior, we report a simple chemically fueled system that exhibits carbodiimide-induced geometry changes. Bridging via transient anhydride formation leads to a significant reduction of the twist about the biaryl bond of substituted diphenic acids, giving a simple molecular clamp. The kinetics are well-described by a simple mechanism, allowing structure–property effects to be determined. The kinetic parameters can be used to derive important characteristics of the system such as the efficiencies (anhydride yields), maximum anhydride concentrations, and overall lifetimes. Transient diphenic anhydrides tolerate steric hindrance ortho to the biaryl bond but are significantly affected by electronic effects, with electron-deficient substituents giving lower yields, peak conversions, and lifetimes. Furthermore, the results provide useful guidelines for the design of functional systems incorporating diphenic acid units.

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

Effect of Salt Concentration on the Interfacial Solvation Structure and Early Stage of Solid–Electrolyte Interphase Formation in Ca(BH 4 ) 2 /THF for Ca Batteries

The Ca 2+ solvation structure at the electrolyte/electrode interface is of central importance to understand electroreduction stability and solid–electrolyte interphase (SEI) formation for the novel multivalent Ca battery systems. Here, using an exemplar electrolyte, the concentration-dependent solvation structure of Ca(BH 4 ) 2 -tetrahydrofuran on a gold model electrode has been investigated with various electrolyte concentrations via electrochemical quartz crystal microbalance with dissipation (EQCM-D) and X-ray photoelectron spectroscopy (XPS). For the first time, in situ EQCM-D results prove that the prevalent species adsorbed at the interface is CaBH 4 + across all concentrations. As the salt concentration increases, the number of BH 4 – anions associated with Ca 2+ increases, and much larger solvated complexes such as CaBH 4 + ·4THF or Ca(BH 4 ) 3 – ·4THF form at the interface at high concentrations prior to Ca plating. Different interfacial chemistries lead to the formation of SEIs with different components demonstrated by XPS. High electrolyte concentrations reduce the solvent decomposition and promote the formation of thick, uniform, and inorganic-rich (i.e., CaO) SEI layers, which contribute to improved Ca plating efficiency and current density in electrochemical measurements.

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

Data-driven learning of Mori–Zwanzig operators for isotropic turbulence

Developing reduced-order models for turbulent flows, which contain dynamics over a wide range of scales, is an extremely challenging problem. In statistical mechanics, the Mori–Zwanzig (MZ) formalism provides a mathematically exact procedure for constructing reduced-order representations of high-dimensional dynamical systems, where the effects due to the unresolved dynamics are captured in the memory kernel and orthogonal dynamics. Turbulence models based on MZ formalism have been scarce due to the limited knowledge of the MZ operators, which originates from the difficulty in deriving MZ kernels for complex nonlinear dynamical systems. In this work, we apply a recently developed data-driven learning algorithm, which is based on Koopman's description of dynamical systems and Mori's linear projection operator, on a set of fully resolved isotropic turbulence datasets to extract the Mori–Zwanzig operators. With data augmentation using known turbulence symmetries, the extracted Markov term, memory kernel, and orthogonal dynamics are statistically converged and the generalized fluctuation–dissipation relation can be verified. The properties of the memory kernel and orthogonal dynamics, and their dependence on the choices of observables are investigated to address the modeling assumptions that are commonly used in MZ-based models. A series of numerical experiments are then constructed using the extracted kernels to evaluate the memory effects on prediction. The results show that the prediction errors are strongly affected by the choice of observables and can be further reduced by including the past history of the observables in the memory kernel.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Heat generation due to the Anderson catastrophe in mesoscopic devices

Anderson's orthogonality catastrophe (AOC) theorem establishes that the ground state of the many-body fermion system is asymptotically orthogonal to the ground state of the same system perturbed by a scattering potential so that the overlap between the original and new ground states decays to zero with the system size. We adopt the AOC for a description of heat production in a complementary metal-oxide-semiconductor transistor. We find that the heat released in the transistor comprises two distinct components, contribution from the dissipation accompanying electron transmission under the applied voltage and purely quantum-mechanical AOC part due to the change in scattering matrix for electrons upon switching between high and low conductance regimes. We calculate the AOC-induced heat production, which we call switching heat.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗