Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equations of motion”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ion Transport in Concentrated Crosslinked Solid Polymer Electrolytes

Crosslinking polymers is a common approach to create mechanically stable solid materials such as polymer electrolytes for lithium batteries. In conventional liquid electrolytes, the solvent molecules move freely to accommodate the field-induced motion of ions. However, in crosslinked polymer electrolytes, the rearrangement of polymer segments is constrained by the deformation limits of the network. Herein, we develop a new transport model that accounts for both the formation of concentration gradients and the elasticity of the electrolyte. The elasticity is incorporated by adding an additional term related to the entropy of crosslinked strands to the electrochemical potential of the salt. The resulting Crosslink Model contains two adjustable parameters: $\mathcal{N}$, the average number of monomers in a strand, and λ crit , the maximum strain the network can sustain. These solid-like constraints produce singularities in the governing transport equations, fundamentally altering the concentration profiles. Plateaus in salt concentrations emerge near the electrodes, and network elasticity introduces a strain overpotential. When compared to a Baseline Model ($\mathcal{N}$ → ∞, equivalent to concentrated solution theory), which predicts steepest gradients near the electrodes, both models yield similar current–voltage relationships. Model predictions are compared to electrochemical data for a poly(ethylene oxide)-based crosslinked polymer electrolyte.

Patel, Vivaan [University of California, Berkeley,↗
DOE OSTI2023

An efficient protocol for excited states of large biochromophores

Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Q y transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory’s accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.

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

Magnetic reconnection and thermal equilibration

When a magnetic field is forced to evolve on a time scale sev, as by footpoint motions driving the solar corona or non-axisymmetric instabilities in tokamaks, the magnetic field lines undergo large-scale changes in topology on a time scale approximately an order of magnitude longer than $\tau$ ev . But the physics that allows such changes operates on a time scale eight or more orders of magnitude slower. An analogous phenomenon occurs in air. Temperature equilibration occurs on a time scale approximately an order of magnitude longer than it takes air to cross a room, $\tau$ ev , although the physical mechanism that allows temperature equilibration is approximately four orders of magnitude slower than $\tau$ ev . The use of Lagrangian coordinates allows the fundamental equations to be solved and both phenomena explained. The theories of thermal equilibration and magnetic reconnection are developed in parallel to help readers obtain an understanding of the importance and implications of analyses using Lagrangian coordinates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Noise-specific beating in the higher-level Ramsey curves of a transmon qubit

We report: in the higher levels of superconducting transmon devices and more generally charge sensitive devices, T 2 * measurements made in the presence of low-frequency time-correlated 1/f charge noise and quasiparticle-induced parity flips can give an underestimation of the total dephasing time. The charge variations manifest as beating patterns observed in the overlay of several Ramsey fringe curves and are reproduced with a phenomenological Ramsey curve model, which accounts for the charge variations. T 2 * dephasing times, which more accurately represent the total dephasing time, are obtained. The phenomenological model is compared with a Lindblad master equation model. Both models are found to be in agreement with one another and the experimental data. Finally, the phenomenological formulation enables a simple method in which the power spectral density for the low-frequency noise can be inferred from the overlay of several Ramsey curves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Revisiting Theoretical Limits for One Degree-of-Freedom Wave Energy Converters

This work revisits the theoretical limits of one degree-of-freedom wave energy converters (WECs). This work considers the floating sphere used in the OES Task 10 WEC modeling and verification effort for analysis. Analytical equations are derived to determine bounds on displacement amplitude, time-averaged power (TAP), and power-take-off (PTO) force. A unique result found shows that the TAP absorbed by a WEC can be defined solely by the inertial properties and radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force was derived that provides lower and upper bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression only provides a lower bound, as there is theoretically no upper bound. These bounds assist in comparing the performance of the floating sphere if it were to extract energy using surge or heave motion. The analysis shows because of differences in hydrodynamic coefficients for each oscillating mode, there are different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on TAP while utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate TAP with a nonideal PTO is modified by the electrical conversion efficiency and the ratio of the PTO spring and damping coefficients. Furthermore, the PTO spring and damping coefficients were separated in the expression, allowing for limits to be set on the PTO coefficients to ensure net power generation.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

The Spiral Galaxies Flat Rotational Velocity Curve Explained by the Constant Group Velocity of a Nonlinear Density Wave

The rotation velocity curves of stars in galaxies, the motions of pairs of galaxies, and the behavior of galaxies in clusters and super-clusters all indicate that there is a lack of mass on different scales in the universe. In this paper, we derive the expression for rotational velocity using the nonlinear density wave theory considering only stellar components and we show that such theory can support the observed flat rotational velocity curve due to the main property of the soliton wave, which is a constant group velocity of the wave. The surface mass density (SMD) function, used to derive gravitational potential gradient and rotational velocity, is not assumed but rather derived as a solution of the nonlinear Srödinger equation, on the contrary to the widely used, in the literature, exponential disk approximation. Three parameters relevant to the curve shape are the intensities of equilibrium SMD, the amplitude of the wave, and total angular velocity or differential rotation, equivalently. Since the shape of the rotational velocity is highly sensitive to the mentioned parameters, this theory eventually provides a method for a very accurate estimation of galaxy mass and angular velocity as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum Gravity and Laser Interferometry: Towards Observable Predictions

Understanding quantum gravity remains one of the deepest challenges in modern physics, as direct experimental access to Planck-scale effects is beyond current technological reach. However, recent theoretical advances indicate that quantum fluctuations of spacetime may produce measurable effects in precision experiments, particularly near causal horizons. This opens new avenues for testing quantum gravity phenomena through high-precision measurement techniques. This dissertation develops multiple theoretical models to characterize these effects and examines their potential observational signatures in future gravitational wave interferometers. We begin by investigating the role of quantum fluctuations in near-horizon geometries through the lens of the AdS/CFT correspondence, which provides a powerful framework for understanding the interplay between quantum field theory and general relativity via holographic principles. By modeling stochastic energy-momentum sources in Rindler-AdS spacetime, we demonstrate that vacuum fluctuations transform the Einstein equations into a Langevin-type stochastic differential equation, leading to potentially observable fluctuations in photon traversal times. Extending this approach to Minkowski spacetime, we establish a correspondence between gravitational shockwaves and fluid dynamics, showing that near-horizon perturbations satisfy an equation analogous to that governing incompressible fluids, thereby reinforcing the membrane paradigm and hydrodynamic analogies in the context of the fluid/gravity duality. Furthermore, we construct the covariant phase space of a spherically symmetric causal diamond in Minkowski spacetime, identifying two fundamental charges that govern its evolution. These results provide a foundation for quantizing causal horizons and understanding their microscopic degrees of freedom. Building upon these theoretical developments, we further examine a related stochastic phenomenon: the gravitational wave memory background arising from the cumulative memory steps produced by supermassive black hole mergers. After reviewing the standard stochastic gravitational wave background, gravitational memory effects, and BMS symmetries, we model the stochastic memory background using a Brownian motion framework. We show that while the cumulative memory background initially appears above the sensitivity curve of space-based interferometers like LISA, the realistic subtraction of individually resolvable merger events substantially suppresses the residual signal, making its detection more challenging. This highlights the critical importance of source subtraction when evaluating the detectability of gravitational memory effects. By bridging fundamental theory with experimental prospects, this dissertation contributes to the ongoing effort to uncover the quantum nature of spacetime through precision measurement techniques. Whether through detecting quantum spacetime fluctuations, gravitational memory backgrounds, or probing the symmetries of causal horizons, the pursuit of observable quantum gravity phenomena continues to expand the frontiers of both theory and experiment.

Zhang, Yiwen [Caltech] (ORCID:0000000323559416)↗
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae

Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T /∣ W ∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

Astronomy & Astrophysics↗
DOE OSTI2024

Nature of polariton transport in a Fabry-Perot cavity

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (nonzero) in-plane wave vectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present the results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. Here, we address this problem with two approaches: the standard transfer-matrix method (TMM) for planar structures and a second simulation based on the numerical solution of the Maxwell-Bloch equations that can be used for general configurations (complex metasurfaces, for instance) and excitation modes. The agreement between the TMM and the full numerical calculation with the Maxwell-Bloch equations when applied to multilayer planar configurations provides proof of the validity of both approaches for the present analyses. Overall, we find that the transport of the molecular excitons inside the cavity synchronize with the evolution of the enhanced electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. Finally, we relate the group velocity to the Hopfield coefficient obtained from quantum modeling and suggest that the dependence of light-matter coupling on the in-plane wave vector can be an important but overlooked factor for understanding the transport behavior of polaritons. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps us to better understand how experimental observations of polaritons should be interpreted.

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

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2022

The Pantheon+ analysis: Improving the redshifts and peculiar velocities of Type Ia supernovae used in cosmological analyses

We examine the redshifts of a comprehensive set of published Type Ia supernovae, and provide a combined, improved catalogue with updated redshifts. We improve on the original catalogues by using the most up-to-date heliocentric redshift data available; ensuring all red shifts have uncertainty estimates; using the exact formulae to convert heliocentric redshifts into the Cosmic Microwave Background (CMB) frame; and utilising an improved peculiar velocity model that calculates local motions in redshift-space and more realistically accounts for the external bulk flow at high-redshifts. We review 2607 supernova redshifts; 2285 are from unique supernovae and 322 are from repeat observations of the same supernova. In total, we updated 990 unique heliocentric redshifts, and found 5 cases of missing or incorrect heliocentric corrections, 44 incorrect or missing supernova coordinates, 230 missing heliocentric or CMB frame redshifts, and 1200 missing redshift uncertainties. The absolute corrections range between 10 -8 ≤ Δz ≤ 0.038, and RMS(Δz) ~ 3 × 10 -3 . The sign of the correction was essentially random, so the mean and median corrections are small: 4 × 10 -4 and 4 × 10 -6 respectively. We examine the impact of these improvements for H 0 and the dark energy equation of state w and find that the cosmological results change by ΔH 0 = -0.12 km s -1 Mpc -1 and Δw = 0.003, both significantly smaller than previously reported uncertainties for H 0 of 1.0 km s -1 Mpc -1 and w of 0.04 respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

X‑ray Coherent Diffractive Imaging of Large Helium Nanodroplets Doped with Small Molecules

We report the first X-ray coherent diffractive imaging experiment on molecule-doped helium nanodroplets. It complements previous work, where we reported single-shot X-ray coherent diffractive imaging studies of Xe dopant clusters formed in 4He and 3He droplets. These noble gas clusters were used to visualize the impact of rotational excitation of the droplets on the spatial distribution of atomic dopants within the droplets, and to study the differences and connections between quantum and classical droplet rotational motion. Here, we expand our studies to the molecular dopants CF4, CHF3, CH3CN, and SF6, imaged with 1.5 keV photons. We find multiple Bragg spots in the diffraction patterns of molecule-doped droplets with radii of approximately 600 nm, which provide evidence that molecules form elongated clusters with preferential alignment along the angular momentum axis of the 4He droplets, in agreement with our previous results on the aggregation of Xe clusters on quantum vortices. Real-space reconstructions of molecular dopant cluster density profiles are obtained for droplets with smaller radii of approximately 300 nm. The diffuse images suggest the formation of low-density, potentially porous, molecular clusters upon aggregation at T = 0.4 K in 4He droplets. In the normal fluid 3He droplets, molecules aggregate into loose clusters on the droplets' equator, similar to previous observations for Xe atoms. Time-of-flight mass spectra reveal that the doped helium nanodroplet moieties fragment extensively into constituent atomic ions, producing only a small fraction of molecular fragment ions. The findings are discussed in the context of previously proposed schemes to use He droplets as potential tamper materials for ultrafast X-ray imaging experiments.

Feinberg, AlexandraJ↗
DOE OSTI2024

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The weakly nonlinear response and nonaffine interpretation of the Johnson–Segalman/Gordon–Schowalter model

We derive new analytical solutions for the non-affine Johnson-Segalman/Gordon-Schowalter (JS/GS) constitutive equation with a general relaxation kernel in medium-amplitude oscillatory shear (MAOS) deformation. The results show time-strain separable (TSS) nonlinearity, therefore providing new physically-meaningful interpretation to the heuristic TSS nonlinear parameter in MAOS. The upper-convected, lower-convected, and corotational Maxwell models are all subsets of the results presented here. The model assumes that the microscale elements causing stress in the material slip compared to the continuum deformation. We introduce a visualization of the non-affine deformation field that acts on stress-generating elements to reinforce the hysical interpretation of the JS/GS class of models. Lastly, a case study is presented where previously published results, from fitting TSS models to MAOS data, can be re-interpreted based on the concept of non-affine motion of the JS/GS framework.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Thermalization of the Kerr index of refraction in acetone and methanol using femtosecond pump-probe scattering spectroscopy

Femtosecond time-resolved pump-probe scattering is used to temporally investigate the index thermalization Δn(t) of the molecular motions associated with the Kerr nonlinear index n 2 . This is done using an intense 300 fs 1034 nm pump pulse and temporally probed by a 517 nm and supercontinuum pulse in acetone and methanol. The optically pumped molecular states change the index of refraction and through the processes of Raman and Rayleigh scattering are shown to engage in a 4-step energy relaxation process. Here, the pump first initiates the Kerr effect as well as populate resonant vibrational molecular motions, which act as a “mother” energy state. The electronic component of the Kerr effect is shown to couple with these vibrational modes through the process of Born-Oppenheimer coupling, which can be measured through the Raman scattering of the probe. The “mother” coupled state then decays into “daughter” anharmonic non-resonant states, which then further decay into “granddaughter” and thermal bath states. The population of these states alter the index of refraction in time, Δn(t), and can be measured through the changing in the probe beam’s Rayleigh scattering. This thermalization process takes ~5.5 ps in acetone and ~3.7 ps in methanol. The scattering signals from each Δn(t) thermalization stage in the localized region, (“mother”, to “daughter”, to “granddaughter”, to the bath states, to the ground state), are fitted to theoretical decay equations that show the rise and fall times of the energy decay routes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Observation and Numerical Simulation of Cold Ions Energized by EMIC Waves

This is the first report of significant energization (up to 7,000 eV) of low-energy He + ions, which occurred simultaneously with H-band electromagnetic ion cyclotron (EMIC) wave activity, in a direction mostly perpendicular to the ambient magnetic field. The event was detected by the Arase satellite in the dayside plasmatrough region off the magnetic equator on 15 May 2019. The peak energy of the He + flux enhancements is mostly above 1,000 eV. At some interval, the He + ions are energized up to ~7,000 eV. The H-band waves are excited in a frequency band between the local crossover and helium gyrofrequencies and are close to a linear polarization state with weakly left-handed or right-handed polarization. The normal angle of the waves exhibits significant variation between 0° and 80°, indicating a non-parallel propagation. Here, we run a hybrid code with parameters estimated from the Arase observations to examine the He + energization. The simulations show that cold He + ions are energized up to more than 1,000 eV, similar to the spacecraft observations. From the analysis of the simulated wave fields and cold plasma motions, we found that the ratio of the wave frequency to He + gyrofrequency is a primary factor for transverse energization of cold He + ions. As a consequence of the numerical analysis, we suggest that the significant transverse energization of He + ions observed by Arase is attributed to H-band EMIC waves excited near the local helium gyrofrequency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Enriched immersed boundary method (EIBM) for interface-coupled multi-physics and applications to convective conjugate heat transfer

It is challenging to develop numerical methods that simultaneously maintain accuracy of resolving boundary conditions and mesh flexibility to handle the interface in interface-coupled multi-physics problems involving large property discontinuity. On the one hand, boundary-fitted methods possess high accuracy in capturing the interfacial phenomenon but involve complicated volumetric mesh generation, mesh-motion, and even re-meshing procedures. On the other hand, immersed boundary methods (IBM) provide mesh flexibility but sometimes suffer from inferior interface representations, leading to poor enforcement of boundary conditions. This paper presents an enriched immersed boundary method (EIBM) to overcome this challenge and demonstrates its efficacy in conjugate heat transfer, a representative example in interface-coupled multi-physics systems that have implications for many industrial processes. Further, the core technique of the method is to enhance IBM’s accuracy of the fluid–solid interface by enriching the degrees of freedom of the cut elements to enforce temperature and flux compatibilities and resolve all the physical unknowns on the background mesh to simplify the volumetric mesh generation. We implement the EIBM under the framework of a variational multiscale formulation for coupled Navier–Stokes and thermodynamics equations. The enriched DoFs enable better enforcement of temperature and flux compatibilities with large conductivity ratios across the fluid–solid interface. At the same time, the immersed nature of the proposed method still attains mesh flexibility. The EIBM’s accuracy is thoroughly evaluated through a set of examples, ranging from benchmark problems with analytical solutions to real-world cooling processes of a moving metallic structure with complex geometry.

42 ENGINEERING↗