Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions

When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A resistive MHD model and simulation on plasma flow evolution in the presence of resonant magnetic perturbation in a tokamak

Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Investigation of compression of indirect-drive targets under conditions of the NIF facility using one-dimensional modelling

The indirect compression dynamics of targets containing capsules with ablators of a plastic, high-density carbon and beryllium is simulated in the framework of a one-dimensional model based on the 1D RADIAN code. Experiments with such targets are performed on the NIF facility in the Livermore Laboratory (USA) in 2014 – 2018. The 1D simulation data are consistent with the results of experiments and calculations made at this Laboratory. The effect of the hard part of hohlraum radiation on capsule compression parameters is confirmed. We demonstrate the feasibility of eliminating this influence not only by selecting the hohlraum material but also by introducing admixtures into the capsule that absorb this radiation. It is shown how varying the amount of admixture in the capsule ablator varies the spectrum of the radiation that heats the DT fuel. (interaction of laser radiation with matter. laser plasma)

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of a Physics-Informed Convolutional Neural Network for Monitoring the Temperature Fields in High-Temperature Gas Reactors

Here, this work presents current advances in applying a physics-informed convolutional neural network (CNN) to evaluate temperature distributions in advanced reactors. Our goal is to demonstrate that the CNN can reconstruct temperature fields within the solid region of a prismatic fuel assembly in a high-temperature gas reactor (HTGR) with sensor data available in only a few cooling channels. Before that, we showcase the superior performance of the physics-informed CNN in comparison to a purely data-driven multilayer perceptron (MLP), considering a canonical heated channel setup. This analysis shows the advantages of our approach and justifies its choice. The datasets employed here are obtained upon numerical simulations performed with codes under the Nuclear Energy Advanced Modeling and Simulation program. This work is important, as industry experience indicates that the assembly material in HTGR concepts is prone to large thermal-mechanical loads nearing operational limits. This makes it crucial to characterize peak temperatures and their distributions near hot spots. Modern thermocouples are unreliable in these types of harsh environments because of the high neutron fluxes and elevated temperatures involved. The CNN-based field reconstruction represents an attractive solution, enabling sensor arrays in less aggressive locations and augmenting indirect predictions for less accessible regions. The results show that the CNN reduces prediction errors by orders of magnitude in comparison to the MLP, considering the simple yet well-representative heated channel case. In the case of the HTGR fuel assembly, the CNN can successfully reconstruct temperature fields over various cooling regimes. Furthermore, we also explore the algorithm’s ability to detect abnormalities. Interestingly, the CNN proves it has the capacity to detect blockage in one of the noninstrumented cooling channels.

Machine learning↗
DOE OSTI2026

Staged Z-pinch experiment at the double-EAGLE pulsed power facility

The staged Z-pinch is a potential high-energy gain fusion concept where a high atomic number liner implodes on a deuterium target using a pulsed multi-MA current source. Over the past several years this concept has been studied on 0.5–1.0 MA facilities using Ar and Kr gas puffs injected at R = 1.2 cm. Neutron yield up to (2.5 ± 0.34) x 10 10 was measured, Ruskov et al (2023 IEEE Trans. Plasma Sci. 51 3310–6). Here, in this study, we present experimental results from the 4 MA, 110 ns current rise time Double-EAGLE facility at L3Harris (currently, Fisica Inc.) where a larger radius nozzle created gas density profiles peaked at R = 2.5 cm. Modeling with the MACH2 and FLASH codes indicates that the larger radius allows stronger acceleration of the liner plasma and generation of shock waves that preheat a target plasma layer next to the liner to temperatures $T_i >$1 keV. The resulting counter thermal pressure on the liner plasma limits the growth of the Magneto-Rayleigh–Taylor (MRT) instability at the liner-vacuum boundary and the implosion proceeds in a relatively stable manner. Near bang time the enormous target plasma thermal pressure smoothens the MRT perturbations developed during the earlier implosion stages. Time integrated x-ray pinhole images with cutoff energy of 100 eV confirm that a long (∼3 cm), stable and uniform high energy density plasma column is formed in the final implosion stage. Consistent neutron yield in the 10 10 –10 11 range was measured for both Ar and Kr liners imploding on a deuterium target.

Physics - Plasma physics↗
DOE OSTI2020

Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks

Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling plasma toroidal flow profile control via NTV torque with n = 2 3D fields in MAST-U

Toroidal modeling utilizing the MARS-Q code shows that large (compared to electromagnetic torque) neoclassical toroidal viscous (NTV) torque can be achieved in the plasma core region in reference MAST(-U) L-mode plasmas, by applying the n=2 (n is the toroidal mode number) 3-D fields generated by the magnetic coils used for controlling the edge localized mode and/or by the error field correction coils. Large NTV torque, occurring at relatively slow plasma flow, in combination with strong variation of the torque amplitude versus the coil phasing of the 3-D coils, offers a tool to control the plasma toroidal rotation profile in spherical tokamaks. To effectively employ the toroidal NTV torque, the potentially fast initial flow needs to first be damped by 3-D fields. In conclusion, MARS-Q quasi-linear initial value simulations demonstrate that this is achievable in MAST-U, within the designed 3-D coil current capability for the device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

First evidence of dominant influence of E × B drifts on plasma cooling in an advanced slot divertor for tokamak power exhaust

Addressing power exhaust in tokamaks is presently recognized as one of the major remaining open issues for the development of fusion reactors. At the forefront of this endeavor is the effort to develop an advanced divertor by maximizing dissipation of plasma power and momentum inside the divertor. Here, we demonstrate, for the first time, that the electromagnetic (E × B) drifts exert a key influence on plasma and gas dynamics in a new advanced slot divertor in the DIII-D tokamak, named the small angle slot (SAS). SAS leverages the effect of drifts to achieve a highly dissipative divertor with electron temperature T e ≲ 10 eV over a wide range of plasma densities, for ion B × ∇B away from the divertor, as used for advanced tokamak operation on DIII-D. Modeling with the SOLPS-ITER code shows that for this drift direction, the E × B flow carries particles toward the outer common flux region (CFR) via the private flux region (PFR), reinforcing neutral recycling and enhancing divertor dissipation. In contrast, for the opposite field direction, the E × B flow carries particles away from the outer CFR into the PFR, offsetting the anticipated SAS geometric effects. Furthermore, this finding is an important step in the understanding of the behavior of advanced divertors for achieving a power exhaust solution for fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear two-fluid modeling of plasma response to RMPs for the ELM control in the ITER baseline

Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. Furthermore, GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as ${I}_{\text{RMP}}^{0.5}$ and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95 . The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. In this work, a set of nonlinear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The nonlinear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP β increases, the EP β threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the nonlinear destabilization of the 3/4 and 2/3 energetic particle modes, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗
DOE OSTI2022

Simulation of the TAEs’ saturation phase in the Large Helical Device: MHD burst

Here, the aim of the present study is to analyze the saturation regime of the toroidal Alfven eigenmodes (TAE) in the Large Helical Device plasma, particularly the MHD burst. The linear and nonlinear evolution of the TAEs are simulated by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The linear simulations indicate the overlapping of 1/2–1/1, 2/3–2/4 and 3/5–3/6 TAEs in the inner-middle plasma region and frequency range of 45–75 kHz, triggered by EPs with an energy of T f = 45 keV and EP β = 0.022. The nonlinear simulations show that 2/3–2/4 and 3/4–3/5 TAEs are further destabilized due to the energy transfer from the 1/1–1/2 TAE, leading to broad TAE radial overlapping and triggering of the MHD burst. The energy of the 1/1–1/2 TAE is also nonlinearly transferred to the thermal plasma destabilizing the 0/0 and 0/1 modes, inducing the generation of shear flows and zonal currents, as well as large deformations in the thermal pressure and EP density radial profiles. The nonlinear simulation reproduces the same succession of instabilities and the same frequency range with respect to the experiment. The instability propagates outward during the bursting phase, showing a large decrease of the EP density profile between the middle-outer plasma, indicating the loss of part of the EP population that explains the decrease in the plasma heating efficiency observed during the MHD burst.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SetBERT: the deep learning platform for contextualized embeddings and explainable predictions from high-throughput sequencing

MOTIVATION: High-throughput sequencing (HTS) is a modern sequencing technology used to profile microbiomes by sequencing thousands of short genomic fragments from the microorganisms within a given sample. This technology presents a unique opportunity for artificial intelligence to comprehend the underlying functional relationships of microbial communities. However, due to the unstructured nature of HTS data, nearly all computational models are limited to processing DNA sequences individually. This limitation causes them to miss out on key interactions between microorganisms, significantly hindering our understanding of how these interactions influence the microbial communities as a whole. Furthermore, most computational methods rely on post-processing of samples which could inadvertently introduce unintentional protocol-specific bias. RESULTS: Addressing these concerns, we present SetBERT, a robust pre-training methodology for creating generalized deep learning models for processing HTS data to produce contextualized embeddings and be fine-tuned for downstream tasks with explainable predictions. By leveraging sequence interactions, we show that SetBERT significantly outperforms other models in taxonomic classification with genus-level classification accuracy of 95%. Furthermore, we demonstrate that SetBERT is able to accurately explain its predictions autonomously by confirming the biological-relevance of taxa identified by the model. AVAILABILITY AND IMPLEMENTATION: All source code is available at https://github.com/DLii-Research/setbert. SetBERT may be used through the q2-deepdna QIIME 2 plugin whose source code is available at https://github.com/DLii-Research/q2-deepdna.

Ludwig, David W↗
DOE OSTI2020

A low-luminosity core-collapse supernova very similar to SN 2005cs

ABSTRACT We present observations and analysis of PSN J17292918+7542390, a low-luminosity Type II-P supernova (LL SN IIP). The observed sample of such events is still low, and their nature is still under debate. Such SNe are similar to SN 2005cs, a well-observed LL Type II-P event, having low expansion velocities, and small ejected 56Ni mass. We have developed a robust and relatively fast Monte Carlo code that fits semi-analytic models to light curves of core-collapse SNe. This allows the estimation of the most important physical parameters, like the radius of the progenitor star, the mass of the ejected envelope, the mass of the radioactive nickel synthesized during the explosion, among others. PSN J17292918+7542390 has $R_0 = 91_{-70}^{+119} \times 10^{11} \, \text{cm}$, $M_\text{ej} = 9.89_{-1.00}^{+2.10} \, \mathrm{ M}_{\odot }$, $E_{\mbox{kin}} = 0.65_{-0.18}^{+0.19} \, \text{foe}$, and $v_{\mbox{exp}} = 3332_{-347}^{+216}$ km s−1, for its progenitor radius, ejecta mass, kinetic energy, and expansion velocity, respectively. The initial nickel mass of the PSN J17292918+7542390 turned out to be $1.55_{-0.70}^{+0.75} \times 10^{-3} \, \mathrm{M}_{\odot }$. The measured photospheric velocity at the earliest observed phase is 7000 km s−1. As far as we can tell based on the small population of observed LL SNe IIP, the determined values are typical for these events.

Jäger, Jr, Zoltán↗
DOE OSTI2025

Enhancing quantum memory lifetime with measurement-free local error correction and reinforcement learning

Reliable quantum computation requires systematic identification and correction of errors that occur and accumulate in quantum hardware. To diagnose and correct such errors, standard quantum error-correcting protocols utilize global error information across the system obtained by mid-circuit readout of ancillary qubits. We investigate circuit-level error-correcting protocols that are measurement-free and based on local error information. Such a local error correction (LEC) circuit consists of faulty multi-qubit gates to perform both syndrome extraction and ancilla-controlled error removal. We develop and implement a reinforcement learning framework that takes a fixed set of faulty gates as inputs and outputs an optimized LEC circuit. To evaluate this approach, we quantitatively characterize an extension of logical qubit lifetime by a noisy LEC circuit. For the two-dimensional (2D) classical Ising model and four-dimensional toric code, our optimized LEC circuit performs better at extending a memory lifetime compared with a conventional LEC circuit based on Toom's rule in a subthreshold gate error regime. We further show that such circuits can be used to reduce the rate of mid-circuit readouts to preserve a 2D toric code memory. Lastly, we discuss the application of the LEC protocol on dissipative preparation of quantum states with topological phases.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Fission product yields from the 238 U ( n,f ) reaction at E n = 4.6 MeV

Here, we present cumulative fission product yields from the 238 U (n, f) reaction at an incident energy of 4.6 MeV for fission products with half-lives from less than one second to several hours. We employed direct γ-ray spectroscopy to measure unseparated fission products with a conventional neutron activation analysis and a cyclic activation technique using the RApid Belt-driven Irradiated Target Transfer System (RABITTS). Yields were determined for 78 isotopes and isomers, covering the mass range 84 ≤ A ≤ 148 and the charge range from 33(As) ≤ Z ≤ 58(Ce). Our experimental cumulative yields are compared to the literature and evaluated data at E n ≈0.5 MeV from ENDF/B-VIII.0 and JEFF-3.3 libraries and with theoretical models using the GEF and FREYA codes at E n =4.6 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗