Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid simulations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermal ion kinetic effects and Landau damping in fishbone modes

We report the kinetic–magnetohydrodynamic (MHD) hybrid simulation approach for macroscopic instabilities in plasmas can be extended to include the kinetic effects of both thermal ions and energetic ions. The new coupling scheme includes synchronization of the density and parallel velocity between thermal ions and MHD, in addition to pressure coupling, to ensure the quasineutrality condition and avoid numerical errors. The new approach has been implemented in the kinetic-MHD code M3D-C1-K, and was used to study the thermal ion kinetic effects and Landau damping in fishbone modes in both DIII-D and NSTX. It is found that the thermal ion kinetic effects can cause an increase of the frequencies of the non-resonant n = 1 fishbone modes driven by energetic particles for q min >1, and Landau damping can provide additional stabilization effects. A nonlinear simulation for n = 1 fishbone mode in NSTX is also performed, and the perturbation on magnetic flux surfaces and the transport of energetic particles are calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Magnetopause Reconnection and Indents Induced by Foreshock Turbulence

Abstract Based on global hybrid simulation results, we predict that foreshock turbulence can reach the magnetopause and lead to reconnection as well as Earth‐sized indents. Both the interplanetary magnetic field (IMF) and solar wind are constant in our simulation, and hence, all dynamics are generated by foreshock instabilities. The IMF in the simulation is mostly Sun‐Earth aligned with a weak northward and zero dawn‐dusk component, such that subsolar magnetopause reconnection is not expected without foreshock turbulence modifying the magnetosheath fields. We show a reconnection example to illustrate that the turbulence can create large magnetic shear angles across the magnetopause to induce local bursty reconnection. Magnetopause reconnection and indents developed from the impact of foreshock turbulence can potentially contribute to dayside loss of planetary plasmas.

58 GEOSCIENCES↗
DOE OSTI2021

Magnetic Reconnection Inside Solar Wind Rotational Discontinuity During Its Interaction With the Quasi-Perpendicular Bow Shock and Magnetosheath

Here, using a three-dimensional global hybrid simulation, we investigate the formation and evolution of ion-scale magnetic reconnection inside an interplanetary rotational discontinuity (RD) owing to its interaction with the quasi-perpendicular (Q-⟂) bow shock and the magnetosphere. The interplanetary magnetic field (IMF) is initially predominantly northward, while it changes to purely southward across the RD. A significantly thinned RD current layer with a width ~ ion skin depth and normal magnetic field B n ≈0 is formed by a shock compression process as the RD interacts with the Q-⟂ shock. Magnetic reconnection thereupon takes place inside the thin RD current layer, where Hall magnetic and electric fields, reconnection electric field, and high-speed ion outflow jets are identified. Simultaneously, flux ropes form with an extension of a few ion inertial lengths. As the RD is transmitted into the magnetosheath, multiple reconnection sites lead to the formation of longer flux ropes. Moreover, magnetosheath reconnection is also found at new reconnection sites inside the RD. No reconnection is found in the Q-⟂ shock alone outside the RD. The flux ropes propagate poleward and tailward in the magnetosheath. In addition, magnetopause reconnection takes place under the southward IMF on the sunward side of the RD after the RD passes through the magnetopause.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Laser-produced plasmas as drivers of laboratory collisionless quasi-parallel shocks

The creation of a repeatable collisionless quasi-parallel shock in the laboratory would provide a valuable platform for experimental studies of space and astrophysical shocks. However, conducting such an experiment presents substantial challenges. Scaling the results of hybrid simulations of quasi-parallel shock formation to the laboratory highlights the experimentally demanding combination of dense, fast, and magnetized background and driver plasmas required. One possible driver for such experiments is high-energy laser-produced plasmas (LPPs). Preliminary experiments at the University of California, Los Angeles, have explored LPPs as drivers of quasi-parallel shocks by combining the Phoenix Laser Laboratory [Niemann et al., J. Instrum. 7, P03010 (2012)] with a large plasma device [Gekelman et al., Rev. Sci. Instrum. 87, 025105 (2016)]. Beam instabilities and waves characteristic of the early stages of shock formation are observed, but spatial dispersion of the laser-produced plasma prematurely terminates the process. This result is illustrated by experimental measurements and Monte Carlo calculations of LPP density dispersion. The experimentally validated Monte Carlo model is then applied to evaluate several possible approaches to mitigating LPP dispersion in future experiments.

Heuer, P. V. (ORCID:0000000150506606)↗
DOE OSTI2021

Bursty magnetic reconnection at the Earth's magnetopause triggered by high-speed jets

The impact of high-speed jets—dynamic pressure enhancements in the magnetosheath—on the Earth's magnetopause has been observed to trigger local magnetic reconnection. We perform a three-dimensional hybrid simulation to study the magnetosheath and magnetopause under turbulent conditions using a quasi-radial southward interplanetary magnetic field (IMF). In contrast to quasi-steady reconnection with a strong southward IMF, we show that after the impact of a jet on the magnetopause, the magnetopause moves inward, the current sheet is compressed and intensified and signatures of local magnetic reconnection are observed, showing similarities to spacecraft measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Symmetry-specific orientational order parameters for complex structures

A comprehensive framework of characterizing complex self-assembled structures with a set of orientational order parameters is presented. It is especially relevant in the context of using anisotropic building blocks with various symmetries. Two classes of tensor order parameters are associated with polyhedral nematic and bond orientational order. For the latter, a variation of classical bond order parameters (BOPs) is introduced, which takes advantage of the symmetry of constituent particles and/or expected crystalline phases. These symmetrized BOPs can be averaged over an entire system or assigned locally to an individual bond. By combining that with the bond percolation procedure, one is able to identify coherent domains within a self-assembled structure. As a demonstration of the proposed framework, we apply it to a simulated hybrid system that combines isotropic and patchy particles with octahedral symmetry. In conclusion, not only does the methodology allow one to identify individual crystalline domains but also it detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.

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

Velocity-space signatures of shock-drift acceleration at quasi-perpendicular collisionless shocks

The shock-drift acceleration of ions at quasi-perpendicular shocks is a well-known kinetic mechanism for the acceleration of a small fraction of incoming ions to high energy. Here, we use a suite of sixteen hybrid simulations of quasi-perpendicular collisionless shocks over the range of Alfvén Mach number 4.3 ≤ M A ≤ 15.8 (corresponding to a range of fast magnetosonic Mach numbers 2.6 ≤ M f ≤ 9.4) and shock-normal angle 45°≤ θ Bn ≤ 90° to identify the velocity-space signature of shock-drift acceleration using the field-particle correlation technique. We show that the features of the ion velocity distribution in the shock foot and ramp regions can be clearly interpreted by analysis of the single-particle trajectory of a reflected ion through the full 3D-3V phase space. The characteristic features of the velocity-space signature of shock-drift acceleration remain qualitatively robust over the full parameter range of our simulations, providing a potential means for its identification using single-point spacecraft measurements. At higher Alfvén Mach numbers M A ≳8 (M f ≳5), kinetic instabilities generate fluctuations of the electromagnetic fields within the shock transition region, leading to fluctuations in and smearing out of the resulting velocity-space signatures, but the signature remains generally robust and identifiable. The results on the shock-drift acceleration of ions presented here represent a novel means to determine more completely the partitioning of upstream bulk flow kinetic energy into plasma heating, particle acceleration, and electromagnetic fields in collisionless shocks.

Howes, Gregory G. [Univ. of Iowa, Iowa City, IA (U↗
DOE OSTI2021

Studies of a modulated Hall thruster

A typical Hall thruster is powered from a DC power supply and operates with a constant discharge voltage. In operation, the discharge current is subject to strong low frequency oscillations (so-called breathing oscillations). Recent studies have shown that not only can these breathing oscillations be correlated with improved performance, but these oscillations can be induced and controlled by modulating the anode voltage. In this work, a systematic experimental study of the plasma flow in a modulated cylindrical Hall thruster was performed to characterize the effect of natural and modulated breathing oscillations on thruster performance. Measurements suggest that modulating the anode voltage in resonance with the natural breathing frequency does increase the thrust, but a corresponding phase alignment of discharge current and discharge voltage causes the efficiency gains to be insignificant. In addition, the outward shift of the acceleration region causes the plasma plume divergence to increase at the resonance condition and thereby, limit the thrust increase. Further, mechanisms underlying the relative phase between discharge current, ion current, and discharge voltage are investigated experimentally and corroborated with one-dimensional hybrid simulations of the thruster discharge.

modelling↗
DOE OSTI2025

Electron density measurements and calculations in a helium capacitively-coupled radio-frequency plasma

We report a comparison of inferred electron density (n e ) in a He capacitively-coupled plasma, deduced from laser-collision induced fluorescence measurements, with values computed using a hybrid simulation framework based on particle-in-cell/Monte Carlo collisions simulations and a fluid model for excited He atoms. The studies were carried out for gas pressures between 50 mTorr and 1000 mTorr and peak-to-peak radio-frequency (13.56 MHz) voltages between 150 V and 350 V, in a highly symmetric source equipped with plane-parallel electrodes. A good agreement is found between the experimental and modeling results for n e except at the lowest operating voltages and gas pressures. The (effective) electron temperature (T e ) values derived by the two methods agree as well reasonably within the plasma bulk. The simulation results are used to compare the density distributions of He + and various He excited levels and their major populating and de-populating channels at 100 mTorr and 1000 mTorr.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On the Effect of Beating during Nonlinear Frequency Chirping

Spectroscopic analyses of energetic particle (EP) driven bursts of MHD fluctuations in magnetically confined plasmas often exhibit chirps that occur simultaneously in groups of two or more. While the superposition of oscillations at multiple frequencies necessarily causes beating in the signal acquired by a localized external probe, self-consistent hybrid simulations of chirping EP modes in a JT-60U tokamak plasma have demonstrated the possibility of global beating, where the mode's electromagnetic field vanishes globally between beats and reappears with opposite phase [Bierwage et al., Nucl. Fusion 57, 016036 (2017)]. This implies that there can be a single coherent field mode that oscillates at multiple frequencies simultaneously when it is resonantly driven by multiple density waves in EP phase space. Conversely, this means that the EP density waves are mutually coupled and interfere with each other via the jointly driven field, a mechanism ignored in some theories of chirping. In this thesis-style treatise, we study the role of field pulsations in general and beating in particular using the Hamiltonian guiding center orbit-following code ORBIT with a reduced wave-particle interaction model in realistic geometry. Beating is found to drive the evolution of EP phase space structures. Additionally, a key mechanism is the pulsation of effective phase space islands combined with the alternation of their effective O- and X-points due to phase jumps between each beat. Observations: (1) Beating causes density wave fronts to advance radially in a pulsed manner and the resulting chirps become staircase-like. (2) The pulsations facilitate convective transfer of material between neighboring layers of phase space density waves. On the one hand, this may inhibit the early detachment of solitary phase space vortices. On the other hand, it facilitates the accumulation of hole and clump fragments into larger structures. (3) Long-range chirping is observed when massive holes or clumps detach and drift away from the turbulent belt around the seed resonance. It is remarkable that the detached vortices remain robust and,on average, maintain their concentric nested layers while being visibly perturbed by the field's continued beating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Particle Acceleration during Magnetic Reconnection in Solar Flares and Its Parameter Dependence [Slide]

We performed 3D hybrid simulations (particle ions and fluid electrons) of magnetic reconnection and we have achieved, for the first time, nonthermal acceleration of all ion species (up to Iron) into power-law energy spectra. The spectra have a scaling of maximum energy per nucleon as a function of species consistent with satellite observations by NASA’s Parker Solar Probe mission.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Development of a Subchannel Capability for Liquid-Metal Fast Reactors in Pronghorn

This report details the development and demonstration of an entirely new capability in Pronghorn, namely the ability to model liquid-metal fast reactor (LMFR) flow conditions on the engineering scale. We developed two modeling approaches for LMFR that can be used separately or be combined into a hybrid simulation: (1) a modern subchannel capability called Pronghorn-Subchannel for square and hexagonal lattices, and (2) a porous flow capability for LMFR geometries. The report emphasizes the novel aspects of the developed subchannel capability and the interoperability of the subchannel capability, porous flow capability, and multiphysics tools within the multiphysics object oriented simulation environment (MOOSE). Here we demonstrate the ability to: (1) Accurately model subchannel flow in hexagonal lattices; (2) Couple the subchannel flow model to multidimensional finite-element method (FEM) or finite-volume method (FVM) heat conduction models; (3) Model LMFRs using Pronghorn’s porous media FVM approach; (4) Couple porous flow FVM and subchannel models in a single simulation; (5) Explicitly model inter-wrapper flows along with conjugate heat transfer from the intra-element flow; and (6) Demonstrate the numerical robustness of the subchannel algorithm by simulating intra-element flow recirculation in a high-buoyancy, low-flow fuel element.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗
DOE OSTI2022

Parametric Decay of Alfvénic Wave Packets in Nonperiodic Low-beta Plasmas

The parametric decay of finite-size Alfvén waves in nonperiodic low-beta plasmas is investigated using one-dimensional (1D) hybrid simulations. Compared with the usual small periodic system, a wave packet in a large system under the absorption boundary condition shows different decay dynamics, including reduced energy transfer, localized density cavitation, and ion heating. The resulting Alfvén wave dynamics are influenced by several factors relating to this instability, including the growth rate, central wave frequency, and unstable bandwidth. A final steady state of the wave packet may be achieved when the instability does not have enough time to develop within the residual packet, and the packet size shows well-defined scaling dependencies on the growth rate, wave amplitude, and plasma beta. Under the proper conditions, enhanced secondary decay can also be excited in the form of a narrow, amplified wave packet. These results may help to interpret laboratory and spacecraft observations of Alfvén waves, and to refine our understanding of the associated energy transport and ion heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

The Role of Three-dimensional Effects on Ion Injection and Acceleration in Perpendicular Shocks

Understanding the conditions that enable particle acceleration at nonrelativistic collisionless shocks is essential to unveil the origin of cosmic rays. We employ 2D and 3D hybrid simulations (with kinetic ions and fluid electrons) to explore particle acceleration and magnetic field amplification in nonrelativistic perpendicular shocks, focusing on the role of shock drift acceleration and its dependence on the shock Mach number. We perform an analysis of the ion injection process and demonstrate why efficient acceleration is only observed in 3D. In particular, we show that ion injection critically depends on the “porosity” of the magnetic turbulence in the downstream region near the shock, a property describing how easily the postshock region allows particles to traverse it and return upstream without being trapped. This effect can only be properly captured in 3D. Additionally, we explore the impact of numerical resolution on ion energization, highlighting how resolving small-scale turbulence—on scales below the thermal ion gyroradius—is essential for accurately modeling particle injection. Overall, our results emphasize the necessity of high-resolution 3D simulations to capture the fundamental microphysics driving particle acceleration at perpendicular shocks.

Orusa, Luca [Princeton Univ., NJ (United States); ↗
DOE OSTI2024

MFANS 2024 - Formally Proving Characteristics of Cyber-Physical Systems

Cyber-physical systems (CPS) are engineered systems that rely on the smooth integration of computational algorithms and physical elements. This integration presents new challenges for verifying that systems will behave as expected. The goal of this presentation is to present current challenges and potential solutions for the formal verification of cyber-physical systems. For cyber systems, formal methods refer to systematically rigorous mathematical techniques employed in the specification, development, analysis, and verification of both software and hardware systems. Recent advancements in computer science have yielded sophisticated tools specifically designed to address challenges associated with formal methods in complex systems. These tools leverage various foundational concepts such as logic, formal languages, program semantics, type systems, type theory, and automata theory. A notable achievement in the application of formal methods is the seL4 microkernel, claimed to be the first general-purpose operating-system kernel to be verified. Its proof implies the absence of bugs and guarantees that the kernel meets specifications. For physical systems, dynamic and control theory has a history of using rigorous analytic techniques to prove functional correctness. Lyapunov, optimal, classical, modern, and robust control theories all provide rigorous mathematical methods both to analyze system performance and to design controller that can be guaranteed to meet certain objectives. Recent computational techniques like level set theory and reachability analysis provide assertions that a system's state will avoid unsafe regions. Even though success has been independently achieved for cyber systems and physical systems, the integration of such systems creates new challenges. In particular, there is an obvious discrepancy between finite-state machines and infinite-state systems, resulting in different approaches for modeling and analyzing these system. While it is possible to simulate hybrid systems, this provides only a demonstration of a performance and not proof. For hybrid systems, current formal methods and system analysis approaches typically require a workarounds to work on hybrid systems like CPS. This paper will outline the state of the art and limits of current practice for formally verifying CPS and will identify possible research directions that require attention.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Simulation budgeting for hybrid effective field theories

In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model. Our accuracy goals, determined with careful consideration of statistical and systematic uncertainties, are 1% accurate in the high-likelihood range of cosmological parameters, and 2% accurate over a broader parameter space volume for k < 1 h Mpc -1 and z < 3. Focusing in part on the 8-parameter w 0 w a CDM+ m ν cosmological model, we find that < 225 simulations are required to meet our error goals over our wide parameter space, including models with rapidly evolving dark energy, given our simulation and emulator recommendations. For a more restricted parameter space volume, as few as 80 simulations are sufficient. We additionally present simulation forecasts for example use cases, and make the code used in our analyses publicly available. These results offer practical guidance for efficient emulator design and simulation budgeting in future cosmological analyses.

cosmological parameters from LSS↗
DOE OSTI2020

Hybrid EMT-TS Simulation Strategies to Study High Bandwidth MMC-Based HVdc Systems

Modular multilevel converters (MMCs) are widely used in the design of modern high-voltage direct current (HVdc) transmission system. High-fidelity dynamic models of MMCs-based HVdc system require small simulation time step and can be accurately modeled in electro-magnetic transient (EMT) simulation programs. The EMT program exhibits slow simulation speed and limitation on the size of the model and brings certain challenges to test the high-fidelity HVdc model in system-level simulations. This paper presents the design and implementation of a hybrid simulation framework, which enables the co-simulation of the EMT model of Atlanta-Orlando HVdc line and the transient stability (TS) model of the entire Eastern Interconnection system. This paper also introduces the implementation of two high-fidelity HVdc line models simulated at different time steps and discusses a dedicated method for sizing the buffer areas on both sides of the HVdc line. The simulation results of the two HVdc models with different sizes of buffer areas are presented and compared.

Liu, Yuan↗