Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Damping”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Here we developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

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

Dynamics of electrified liquid metal surface using shallow water model

A shallow water model that incorporates surface tension and electric field effects is developed to investigate the dynamics of an electrified liquid surface. The computational model is verified against the Zakharov–Kuznetsov equation and is applied to study the growth and damping of the electrified liquid surface. A linear wave analysis is performed under a shallow water theory assuming an analytic solution of the electric field, similar to the Tonks–Frenkel instability. The electrified liquid surface grows or dampens based on the balance of the electric field, surface tension, and gravitational forces. As a result, the numerical results obtained from the electrified shallow water solver are in good agreement with the theoretical analysis.

Computational models↗
DOE OSTI2023

Variational deep learning of equilibrium transition path ensembles

Here, we present a time-dependent variational method to learn the mechanisms of equilibrium reactive processes and efficiently evaluate their rates within a transition path ensemble. This approach builds off of the variational path sampling methodology by approximating the time-dependent commitment probability within a neural network ansatz. The reaction mechanisms inferred through this approach are elucidated by a novel decomposition of the rate in terms of the components of a stochastic path action conditioned on a transition. This decomposition affords an ability to resolve the typical contribution of each reactive mode and their couplings to the rare event. The associated rate evaluation is variational and systematically improvable through the development of a cumulant expansion. We demonstrate this method in both over- and under-damped stochastic equations of motion, in low-dimensional model systems, and in the isomerization of a solvated alanine dipeptide. In all examples, we find that we can obtain quantitatively accurate estimates of the rates of the reactive events with minimal trajectory statistics and gain unique insights into transitions through the analysis of their commitment probability.

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

Superconducting microwave cavities and qubits for quantum information systems

Superconducting microwave cavities featuring ultrahigh Q-factors, which measure the efficiency of energy storage in relation to energy loss in a system, are revolutionizing quantum computing by providing long coherence times exceeding 1 ms, crucial for the development of scalable multi-qubit quantum systems with low error rates. Here, in this work, we provide an in-depth analysis of recent advances in ultrahigh Q-factor cavities, integration of Josephson junction-based qubits, and bosonic-encoded qubits in 3D cavities. We examine the sources of quantum state dephasing caused by damping and noise mechanisms in cavities and qubits, highlighting the critical challenges that need to be addressed to achieve even higher coherence times. We critically survey the latest progress made in implementing single 3D qubits using superconducting materials, normal metals, and multi-qubit and multi-state quantum systems. Our work sheds light on the promising future of this research area, including novel materials for cavities and qubits, modes with nontrivial topological properties, error correction techniques for bosonic qubits, and new light–matter interaction effects.

3D printing↗
DOE OSTI2023

Error-Field Penetration Thresholds in Ohmically Heated ITER and SPARC Plasmas

The critical n=1, n=2, and n=3 error-field amplitudes needed to trigger error-field penetration in ITER and SPARC, steady-state, ohmically heated plasmas are calculated using a standard asymptotic matching approach. The calculation incorporates plasma impurities, trapped particles, the bootstrap current, and neoclassical poloidal flow-damping. The energy confinement time is specified by the neo-Alcator scaling law in the low-density linear ohmic confinement (LOC) regime, and by the ITER-89P L-mode scaling law in the high-density saturated ohmic confinement regime (SOC) regime. The response of the plasma in the inner region is calculated using a linearized version of the four-field model. At the normal operating electron number density, diamagnetic levels of rotation are found to be sufficient to protect ITER and SPARC ohmically heated plasmas from m=2/n=1 error-field penetration. On other hand, SPARC, and especially ITER, ohmically heated plasmas may be vulnerable to n>1 error-field penetration. ITER and SPARC ohmically heated plasmas are also slightly more susceptible to error-field penetration when the electron fluid at the rational surface rotates in the ion diamagnetic direction, rather than the electron diamagnetic direction. At electron number densities that are sufficiently low that the plasma lies in the LOC confinement regime, the error-field penetration threshold increases with increasing density. However, as soon as the electron number density becomes large enough that the plasma enters the SOC regime, the increase in the error-field penetration threshold with increasing density levels off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry

Extensive experimental evidence has shown that the presence of poloidal flow in tokamaks can dramatically improve transport properties. However, theory indicates that poloidal flows are damped by poloidal viscosity, thus necessitating external drivers, such as neutral beam injection or radio frequency heating. In this work, ideal magnetohydrodynamic equilibria are calculated via the FORTRAN code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)] and a postprocessor is used to estimate the neoclassical poloidal viscosity. The equilibrium inputs, which correspond to intuitive physical quantities, are then numerically optimized to reduce a viscosity figure of merit. We present supersonic equilibria in tokamak geometry with minimized neoclassical poloidal viscosities for various velocity free function inputs, plasma aspect ratios, and collisionality regimes. Benchmarks are made against an analytic theory as well as a classical expression of poloidal viscosity. Numerical confirmation of the analytic theory is obtained in the high aspect ratio and high collisionality limit. Good agreement is also seen near the plasma core and edge, with discrepancies arising in the intermediate region. Outside of these limits, rotation input function profiles are found that provide ∼order of magnitude improvements over the analytic theory, with additional progress being made toward predictions for tokamak-relevant equilibria.

Physics↗
DOE OSTI2024

Modeling stimulated Brillouin backscatter from outer-cone quads across multiple inertial confinement fusion hohlraum designs

Stimulated Brillouin scattering (SBS) is a potential risk for laser damage in the experiments carried out at the National Ignition Facility (NIF), and by altering the energy deposition pattern in hohlraums, it affects the symmetry of indirect-drive inertial confinement fusion implosions. We have surveyed backward SBS on outer-cone quads across NIF integrated hohlraums of various platforms numerically, using three-dimensional (3D) simulations with the backscatter code pF3D [Berger et al. Phys. Plasmas 5, 4337 (1998)] and ray-based gain calculations. Simulated reflected powers and energies, as well as the spectrum of reflected light all compare favorably with measurements. Ray-based calculations of exponential SBS amplification (“gain”), which assume a strongly damped plasma wave and steady-state response, are performed using a novel method that includes the 3D speckled field of the laser that drives SBS. This approach is useful for understanding qualitative differences between hohlraum designs and identifying regions susceptible to SBS within hohlraums. Quantitatively, gains are not found to correlate with SBS reflectivities in 3D, necessitating fully wave-based calculations that naturally include diffraction and various temporal dependencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The theory of kinetic effects on resistive wall mode stability in tokamaks

Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear dynamics of multiple Alfvén modes driven by trapped energetic ions in tokamaks: A triplet paradigm

Non-linear dynamics of multiple infernal Alfvén eigenmodes—a subset of global Alfvén eigenmodes in tokamak plasmas with extended low-shear central core [Marchenko et al., Phys. Plasmas 16, 092502 (2009)]—is studied. The analysis is carried out for a mode triplet with toroidal mode-numbers n=1, 2, 3. It was assumed that the n = 1 mode was linearly unstable due to precession resonance with trapped fast ions, whereas the other modes were linearly damped. The modes were coupled due to a non-linearity in a bounce-averaged drift kinetic equation for the distribution function of fast ions. Nonlinear equations for the mode amplitudes and phases are derived and solved numerically. It is found that the temporal evolution of the amplitudes and the phase (responsible for the frequency chirping) of the modes exhibit Hopf bifurcations to stable limit cycles. This can explain a synchronous cyclic destabilization of multiple modes in Alfvén avalanches (sudden growth of amplitudes of the mode cluster with different n and approximately equal frequency spacing) in NSTX and bursting modes in MAST—events, which resulted in enhanced loss of fast ions.

Marchenko, V. S. (ORCID:0000000261480533)↗
DOE OSTI2024

In situ electric-field control of ferromagnetic resonance in the low-loss organic-based ferrimagnet V[TCNE] x ∼2

We demonstrate indirect electric-field control of ferromagnetic resonance (FMR) in devices that integrate the low-loss, molecule-based, room-temperature ferrimagnet vanadium tetracyanoethylene (V[TCNE]x∼2) mechanically coupled to PMN-PT piezoelectric transducers. Upon straining the V[TCNE]x films, the FMR frequency is tuned by more than 6 times the resonant linewidth with no change in Gilbert damping for samples with α = 6.5 × 10−5. We show this tuning effect is due to a strain-dependent magnetic anisotropy in the films and find the magnetoelastic coefficient |λs| ∼ (1–4.4) ppm, backed by theoretical predictions from density-functional theory calculations and magnetoelastic theory. Noting the rapidly expanding application space for strain-tuned FMR, we define a new metric for magnetostrictive materials, magnetostrictive agility, given by the ratio of the magnetoelastic coefficient to the FMR linewidth. This agility allows for a direct comparison between magnetostrictive materials in terms of their comparative efficacy for magnetoelectric applications requiring ultra-low loss magnetic resonance modulated by strain. With this metric, we show V[TCNE]x is competitive with other magnetostrictive materials, including YIG and Terfenol-D. This combination of ultra-narrow linewidth and magnetostriction, in a system that can be directly integrated into functional devices without requiring heterogeneous integration in a thin film geometry, promises unprecedented functionality for electric-field tuned microwave devices ranging from low-power, compact filters and circulators to emerging applications in quantum information science and technology.

Materials Science↗
DOE OSTI2024

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

We present the development and verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations. We apply the corrected binary collision model originally developed for charged-particles collisions to all considered scattering channels, including Coulomb collisions, elastic neutral–neutral and charged–neutral collisions, ionization, excitation, and fusion. The model's implementation is described and verified through a series of simulations, including charged- and neutral-particle thermal equilibration, slowing of electrons in warm solid-density aluminum, collisional damping of a Langmuir wave, helium gas breakdown in an applied electric field, and thermonuclear and beam–target fusion. Then, we demonstrate the model within simulations of hydrogen plasma formation in the Princeton Field-Reversed Configuration as well as of the burning of aneutronic fusion fuel p-11B. The latter includes measurement of the fusion power density in a low-density plasma and fusion production due to the stopping of a proton ignitor beam in a compressed boron target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2024

Vortex splitting in two-dimensional fluids and non-neutral electron plasmas with smooth vorticity profiles

Initially, elliptical, quasi-two-dimensional (2D) fluid vortices can split into multiple pieces if the aspect ratio is sufficiently large due to the growth and saturation of perturbations known as Love modes on the vortex edge. Presented here are experiments and numerical simulations, showing that the aspect ratio threshold for vortex splitting is significantly higher for vortices with realistic, smooth edges than that predicted by a simple “vortex patch” model, where the vorticity is treated as piecewise constant inside a deformable boundary. The experiments are conducted by exploiting the E × B drift dynamics of collisionless, pure electron plasmas in a Penning–Malmberg trap, which closely model 2D vortex dynamics due to an isomorphism between the Drift–Poisson equations describing the plasmas and the Euler equations describing ideal fluids. The simulations use a particle-in-cell method to model the evolution of a set of point vortices. The aspect ratio splitting threshold ranges up to about twice as large as the vortex patch prediction and depends on the edge vorticity gradient. This is thought to be due to spatial Landau damping, which decreases the vortex aspect ratio over time and, thus, stabilizes the Love modes. Near the threshold, asymmetric splitting events are observed in which one of the split products contains much less circulation than the other. These results are relevant to a wide range of quasi-2D fluid systems, including geophysical fluids, astrophysical disks, and drift-wave eddies in tokamak plasmas.

Physics↗
DOE OSTI2024

A new method for solving the linearized 1D Vlasov–Poisson system yielding a new class of solutions

We describe a new method for solving the linearized 1D Vlasov–Poisson system by using properties of Cauchy-type integrals. Our method remedies critical flaws of the two standard methods, reveals a previously unrecognized Gaussian-in-time-like decay, and can also account for an externally applied electric field. The Landau approximation involves deforming the Bromwich contour around the poles closest to the real axis due to the analytically continued dielectric function, finding the long-time behavior for a stable system: Landau damping. Jackson's generalization encircles all poles while sending the contour to infinity, assuming its contribution vanishes, which is not true in general. This gives incorrect solutions for physically reasonable configurations and can exhibit pathological behavior, of which we show examples. The van Kampen method expresses the solution for a stable equilibrium as a continuous superposition of waves, resulting in an opaque integral. Case's generalization includes unstable systems and predicts a decaying discrete mode for each growing discrete mode, an apparent contradiction to both the Jackson solution and ours. We show, without imposing additional constraints, that the decaying modes are never present in the time evolution due to an exact cancellation with part of the continuum. Our solution is free of integral expressions, is obtained using algebra and Laurent series expansions, does not rely on analytic continuations, and results in a correct asymptotically convergent form in the case of infinite sums. The analysis used can be readily applied in higher-dimensional, electromagnetic systems and also provides a new technique for evaluating certain inverse Laplace transforms.

Physics↗
DOE OSTI2024

Enhanced spin pumping in heterostructures of coupled ferrimagnetic garnets

Spin pumping has significant implications for spintronics, providing a mechanism to manipulate and transport spins for information processing. Understanding and harnessing spin currents through spin pumping is critical for the development of efficient spintronic devices. The use of a magnetic insulator with low damping enhances the signal-to-noise ratio in crucial experiments such as spin-torque ferromagnetic resonance (FMR) and spin pumping. A magnetic insulator coupled with a heavy metal or quantum material offers a more straightforward model system, especially when investigating spin-charge interconversion processes to greater accuracy. This simplicity arises from the absence of unwanted effects caused by conduction electrons unlike in ferromagnetic metals. Here, we investigate the spin pumping in coupled ferrimagnetic (FiM) Y 3 Fe 5 O 12 (YIG)/Tm 3 Fe 5 O 12 (TmIG) bilayers combined with heavy-metal (Pt) using the inverse spin Hall effect. It is observed that magnon transmission occurs at both of the FiMs FMR positions. The enhancement of spin pumping voltage (V sp ) in the FiM garnet heterostructures is observed. The plausible reason might be the interfacial exchange coupling between FiMs. The modulation of V sp is achieved by tuning the bilayer structure. Further, the spin mixing conductance for these coupled systems is found to be ≈ 10 18 m –2 . Furthermore, our findings describe a coupled FiM system for the investigation of magnon coupling providing opportunities for magnonic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The nonlinear behavior of generic tail fins for small wind turbines

This paper describes analysis and measurements of the yaw response of tail fins for small wind turbines. It is based on an extension of unsteady slender body theory (USBT) to cover non-slender fins and high angles of incidence, both of which make the theory nonlinear. We provide three main additions to the substantial literature on linearized USBT for tail fins. First, USBT is extended to high angles by modeling the nonlinear vortex dynamics. Second, the restriction to slender bodies is removed by modeling the chordwise load variation. Third, we consider the effect of time-varying wind speed. Further, the extended theory is compared to wind tunnel measurements of the yaw behavior of delta, elliptical, and rectangular tail fins without a rotor and nacelle. The fins were released from initial yaw angles of -40° and -80°; the latter is of sufficient magnitude to show the importance of the nonlinear yaw dynamics. Generally good agreement was found between the theory and measurements, and the theory was shown to be more accurate than a “polar” or quasi-steady model which uses only the lift and drag of a delta planform. Of the three planforms, the rectangular one showed the lowest accuracy in terms of frequency but the damping was accurately predicted. Overall, the results demonstrate the importance of nonlinearity in the response of a yawing tail fin, particularly for the higher aspect ratio fins at large yaw angles.

17 WIND ENERGY↗
DOE OSTI2024

Toroidal Alfvén mode instability driven by plasma current in low-density Ohmic plasmas of the spherical tori

Due to intrinsically low magnetic fields, at low density the drift speed of the current-carrying electrons in spherical tokamaks can exceed fraction of the Alfvén speed sufficient for the excitation of the Alfvén gap modes. A particular case of the toroidal mode, observed during minor disruptions in Ohmic shots on SUNIST [Liu et al., Phys. Plasmas 23, 120706 (2016)], is considered in the present communication. Due to the negligible effect of the electron pressure gradient, the growth rate scales linearly with the drift speed, with slope inversely proportional to electron thermal velocity. In the absence of continuum damping, the threshold value of the drift speed for TAE excitation is independent of electron temperature.

Marchenko, V. S. (ORCID:0000000261480533)↗
DOE OSTI2024

Spectral broadening and vibronic dynamics of the S2 state of canthaxanthin in the orange carotenoid protein

We have performed a series of broadband multidimensional electronic spectroscopy experiments to probe the electronic and vibrational dynamics of the canthaxanthin chromophore of the Orange Carotenoid Protein (OCP) from Synechocystis sp. PCC 6803 in its photoactivated red state, OCPR. Cross-peaks observed below the diagonal of the two-dimensional electronic spectrum indicate that absorption transitions prepare the bright S2 state of the ketocarotenoid canthaxanthin near to a sequence of conical intersections, allowing passage to the dark S1 state via the Sx intermediate in <50 fs. Rapid damping of excited-state coherent wavepacket motions suggests that the branching coordinates of the conical intersections include out-of-plane deformation and C=C stretching coordinates of the π-conjugated isoprenoid backbone. The unusual proximity of the Franck–Condon S2 state structure to the conical intersections with Sx and S1 suggests that the protein surroundings of canthaxanthin prepare it to function as an excitation energy trap in the OCPR–phycobilisome complex. Numerical simulations using the multimode Brownian oscillator model demonstrate that the ground-state absorption spectrum of OCPR overlaps with the fluorescence emission spectrum of allophycocyanin due to spectral broadening derived especially from the intramolecular motions of the canthaxanthin chromophore in its binding site.

Yang, Kunwei↗