Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “lack”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

n-Butanol production by Rhodopseudomonas palustris TIE-1

Anthropogenic carbon dioxide (CO 2 ) release in the atmosphere from fossil fuel combustion has inspired scientists to study CO 2 to biofuel conversion. Oxygenic phototrophs such as cyanobacteria have been used to produce biofuels using CO 2 . However, oxygen generation during oxygenic photosynthesis adversely affects biofuel production efficiency. To produce n-butanol (biofuel) from CO 2 , here we introduce an n-butanol biosynthesis pathway into an anoxygenic (non-oxygen evolving) photoautotroph, Rhodopseudomonas palustris TIE-1 (TIE-1). Using different carbon, nitrogen, and electron sources, we achieve n-butanol production in wild-type TIE-1 and mutants lacking electron-consuming (nitrogen-fixing) or acetyl-CoA-consuming (polyhydroxybutyrate and glycogen synthesis) pathways. The mutant lacking the nitrogen-fixing pathway produce the highest n-butanol. Coupled with novel hybrid bioelectrochemical platforms, this mutant produces n-butanol using CO 2 , solar panel-generated electricity, and light with high electrical energy conversion efficiency. Overall, this approach showcases TIE-1 as an attractive microbial chassis for carbon-neutral n-butanol bioproduction using sustainable, renewable, and abundant resources.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

CryoTRANS: predicting high-resolution maps of rare conformations from self-supervised trajectories in cryo-EM

Cryogenic electron microscopy (cryo-EM) has revolutionized structural biology, enabling efficient determination of structures at near-atomic resolutions. However, a common challenge arises from the severe imbalance among various conformations of vitrified particles, leading to low-resolution reconstructions in rare conformations due to a lack of particle images in these quasi-stable states. We introduce CryoTRANS, a method that predicts high-resolution maps of rare conformations by constructing a self-supervised pseudo-trajectory between density maps of varying resolutions. This trajectory is represented by an ordinary differential equation parameterized by a deep neural network, ensuring retention of detailed structures from high-resolution density maps. By leveraging a single high-resolution density map, CryoTRANS significantly improves the reconstruction of rare conformations and has been validated on four real-world datasets: alpha-2-macroglobulin, actin-binding protein complexes, SARS-CoV-2 spike glycoprotein, and the 70S ribosome. CryoTRANS can also predict high-resolution structures in cryogenic electron tomography maps using a high-resolution cryo-EM map.Cryogenic electron microscopy (cryo-EM) has revolutionized structural biology, enabling efficient determination of structures at near-atomic resolutions. However, a common challenge arises from the severe imbalance among various conformations of vitrified particles, leading to low-resolution reconstructions in rare conformations due to a lack of particle images in these quasi-stable states. We introduce CryoTRANS, a method that predicts high-resolution maps of rare conformations by constructing a self-supervised pseudo-trajectory between density maps of varying resolutions. This trajectory is represented by an ordinary differential equation parameterized by a deep neural network, ensuring retention of detailed structures from high-resolution density maps. By leveraging a single high-resolution density map, CryoTRANS significantly improves the reconstruction of rare conformations and has been validated on four real-world datasets: alpha-2-macroglobulin, actin-binding protein complexes, SARS-CoV-2 spike glycoprotein, and the 70S ribosome. CryoTRANS can also predict high-resolution structures in cryogenic electron tomography maps using a high-resolution cryo-EM map.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Mycobacterium tuberculosis CitA activity is modulated by cysteine oxidation and pyruvate binding

As an adaptation for survival during infection, Mycobacterium tuberculosis becomes dormant, reducing its metabolism and growth. Two types of citrate synthases have been identified in Mycobacterium tuberculosis, GltA2 and CitA. Previous work shows that overexpression of CitA, the secondary citrate synthase, stimulates the growth of Mycobacterium tuberculosis under hypoxic conditions without showing accumulation of triacylglycerols and makes mycobacteria more sensitive to antibiotics, suggesting that CitA may play a role as a metabolic switch during infection and may be an interesting TB drug target. To assess the druggability and possible mechanisms of targeting CitA with small-molecule compounds, the CitA crystal structure was solved to 2.1 Å by X-ray crystallography. The solved structure shows that CitA lacks an NADH binding site that would afford allosteric regulation, which is atypical of most citrate synthases. However, a pyruvate molecule is observed within the analogous domain, suggesting pyruvate may instead be the allosteric regulator for CitA. Further, the R149 and R153 residues forming the charged portion of the pyruvate binding pocket were mutated to glutamate and methionine, respectively, to assess the effect of mutations on activity. Protein thermal shift assay shows thermal stabilization of CitA in the presence of pyruvate compared to the two CitA variants designed to decrease pyruvate affinity. Solved crystal structures of both variants show no significant structural changes. However, the catalytic efficiency of the R153M variant increases by 2.6-fold. Additionally, we show that covalent modification of C143 of CitA by Ebselen completely arrests enzyme activity. Similar inhibition is observed using two spirocyclic Michael acceptor containing compounds, which inhibit CitA with IC app 50 values of 6.6 and 10.9 μM. A crystal structure of CitA modified by Ebselen was solved, but significant structural changes were lacking. Considering that covalent modification of C143 inactivates CitA and the proximity of C143 to the pyruvate binding site, this suggests that structural and/or chemical changes in this sub-domain are responsible for regulating CitA enzymatic activity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Open-source generation of sigma profiles: impact of quantum chemistry and solvation treatment on machine learning performance

The combination of machine learning (ML) models with chemistry-related tasks requires the description of molecular structures in a machine-readable way. The nature of these so-called molecular descriptors has a direct and major impact on the performance of ML models and remains an open problem in the field. Structural descriptors like SMILES strings or molecular graphs lack size-independence and can be memory intensive. Machine-learned descriptors can be of low dimensionality and constant size but lack physical significance and human interpretability. Sigma profiles, which are unnormalized histograms of the surface charge distributions of solvated molecules, combine physical significance with low dimensionality and size-independence, making them a suitable candidate for a universal molecular descriptor. However, their widespread adoption in ML applications requires open access to sigma profile generation, which is currently not available. This work details the development of OpenSPGen – an open-source tool for generating sigma profiles. Also presented are studies on the effect of different settings on the efficacy of the generated sigma profiles at predicting thermophysical material properties when used as inputs to a Gaussian process as a simple surrogate ML model. We find that a higher level of theory does not translate to more accurate results. We also provide further recommendations for sigma profile calculation and use in ML models.

Salih, Fathya Y. M. [University of Notre Dame, IN ↗
DOE OSTI2021

Conservation of endo -glucanase 16 (EG16) activity across highly divergent plant lineages

Plant cell walls are highly dynamic structures that are composed predominately of polysaccharides. As such, endogenous carbohydrate active enzymes (CAZymes) are central to the synthesis and subsequent modification of plant cells during morphogenesis. The endo-glucanase 16 (EG16) members constitute a distinct group of plant CAZymes, angiosperm orthologs of which were recently shown to have dual β-glucan/xyloglucan hydrolase activity. Molecular phylogeny indicates that EG16 members comprise a sister clade with a deep evolutionary relationship to the widely studied apoplastic xyloglucan endo-transglycosylases/hydrolases (XTH). A cross-genome survey indicated that EG16 members occur as a single ortholog across species and are widespread in early diverging plants, including the non-vascular bryophytes, for which functional data were previously lacking. Remarkably, enzymological characterization of an EG16 ortholog from the model moss Physcomitrella patens (PpEG16) revealed that EG16 activity and sequence/structure are highly conserved across 500 million years of plant evolution, vis-à-vis orthologs from grapevine and poplar. Ex vivo biomechanical assays demonstrated that the application of EG16 gene products caused abrupt breakage of etiolated hypocotyls rather than slow extension, thereby indicating a mode-of-action distinct from endogenous expansins and microbial endo-glucanases. As a result, the biochemical data presented here will inform future genomic, genetic, and physiological studies of EG16 enzymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Issues in Bi-Directional Interoperability between BIM and BEM

The state-of-the-art interoperability between Building Information Modeling (BIM) and Building Energy Modeling (BEM) is mainly one-way. BEM software like EnergyPlus takes design input from BIM and simulates the energy consumption of the building. BIM, however, is unable to leverage energy simulation results from BEM smoothly. This lack of bi-directional interoperability between BIM and BEM prevents buildings from seeking more energy-efficient design options and limits BIM applications in optimal generative design. In this research, the authors have conducted the following activities related to two-way interoperability between BIM and BEM: (1) surveyed its current status through literature review, (2) investigated barriers to the bi-directional interoperability, and (3) suggested possible approaches to enable bi-directional interoperability and research outlook in the field. As a result, the authors have identified three main barriers to BIM BEM interoperability, namely, inefficient collaboration, limited scope in information transformation, and lack of BEM to BIM transformation. Future directions are suggested to address these barriers.

building design↗
DOE OSTI2025

Life Cycle Cost Analysis of Prestressed Concrete Poles Subjected to Wind, Surges, and Waves

Prestressed concrete (PC) poles are becoming popular choices to support coastal power transmission systems. However, the existing literature does not offer a detailed analysis of the effectiveness of PC poles in terms of long-term vulnerabilities and the direct and indirect costs. This is due to (1) lack of fragility models for PC poles and (2) lack of probabilistic wind, storm surge, and wave models in coastal settings. In this study, we address these gaps through a series of Monte Carlo simulations to estimate fragility of PC poles as a function of age and hazard (wind, surges, and waves) intensity, and the development of a probabilistic hazard model based on 10,000 years of synthetic tropical cyclone data. The probabilistic hazard model is used in conjunction with high-resolution hydrodynamic models to generate realizations of coastal wind, storm surges, and waves for the Louisiana and Mississippi coasts. A comprehensive life cycle cost analysis for a service life of 70 years considering direct and indirect losses is conducted to compare the performance of a transmission line located in Pascagoula, Mississippi, when wood poles are replaced by PC poles. Results showed that aging has a minor effect on the reliability of PC poles, highlighting the advantages of replacing wood poles with PC poles, especially in coastal areas. In addition, PC poles are significantly more cost-effective compared with wood poles over their life cycle, leading to an estimated saving of $11.55 million (68.17% reduction). The results of this study provide key insight to inform decision-making processes to keep the coastal power grids resilient and cost-effective against future storm hazards.

natural disasters↗
DOE OSTI2020

Modeling q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

A toroidal asymptotic matching model of the response of a tokamak plasma to a static resonant magnetic perturbation (RMP) is used to simulate the n = 3 RMP-induced edge-localized-mode-suppression windows in q 95 that are evident when the plasma current is slowly ramped in DIII-D discharge #145380. All quantities employed in the simulation are derived from experimental measurements, apart from the neutral particle data. Three cases are considered. In the first case, the natural frequencies of tearing modes resonant in the plasma are determined by the ion flows at the corresponding resonant surfaces, which is the prediction of nonlinear tearing mode theory. In the second case, the natural frequencies are determined by the local E x B velocities at the resonant surfaces. In the third case, the natural frequencies are determined by the electron flows at the resonant surfaces, which is the prediction of linear tearing mode theory. The second case gives the best agreement between the simulations and the experimental observations. The first and third cases only lead to partial agreement between the simulations and the observations. In the first case, the lack of complete agreement may be a consequence of using an inaccurate assumption for the neutral particle distribution in the pedestal. In the third case, the lack of complete agreement is probably due to the fact that the response of a tokamak plasma to an RMP is not accurately described by linear tearing mode theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Capsular polysaccharide correlates with immune response to the human gut microbe Ruminococcus gnavus

Active inflammatory bowel disease (IBD) often coincides with increases of Ruminococcus gnavus, a gut microbe found in nearly everyone. It was not known how, or if, this correlation contributed to disease. We investigated clinical isolates of R. gnavus to identify molecular mechanisms that would link R. gnavus to inflammation. Here, we show that only some isolates of R. gnavus produce a capsular polysaccharide that promotes a tolerogenic immune response, whereas isolates lacking functional capsule biosynthetic genes elicit robust proinflammatory responses in vitro. Germ-free mice colonized with an isolate of R. gnavus lacking a capsule show increased measures of gut inflammation compared to those colonized with an encapsulated isolate in vivo. These observations in the context of our earlier identification of an inflammatory cell-wall polysaccharide reveal how some strains of R. gnavus could drive the inflammatory responses that characterize IBD.

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

Green diatom mutants reveal an intricate biosynthetic pathway of fucoxanthin

Fucoxanthin is a major light-harvesting pigment in ecologically important algae such as diatoms, haptophytes, and brown algae (Phaeophyceae). Therefore, it is a major driver of global primary productivity. Species of these algal groups are brown colored because the high amounts of fucoxanthin bound to the proteins of their photosynthetic machineries enable efficient absorption of green light. While the structure of these fucoxanthin-chlorophyll proteins has recently been resolved, the biosynthetic pathway of fucoxanthin is still unknown. Here, we identified two enzymes central to this pathway by generating corresponding knockout mutants of the diatom Phaeodactylum tricornutum that are green due to the lack of fucoxanthin. Complementation of the mutants with the native genes or orthologs from haptophytes restored fucoxanthin biosynthesis. We propose a complete biosynthetic path to fucoxanthin in diatoms and haptophytes based on the carotenoid intermediates identified in the mutants and in vitro biochemical assays. It is substantially more complex than anticipated and reveals diadinoxanthin metabolism as the central regulatory hub connecting the photoprotective xanthophyll cycle and the formation of fucoxanthin. Moreover, our data show that the pathway evolved by repeated duplication and neofunctionalization of genes for the xanthophyll cycle enzymes violaxanthin de-epoxidase and zeaxanthin epoxidase. Brown algae lack diadinoxanthin and the genes described here and instead use an alternative pathway predicted to involve fewer enzymes. Our work represents a major step forward in elucidating the biosynthesis of fucoxanthin and understanding the evolution, biogenesis, and regulation of the photosynthetic machinery in algae.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Trajectories for the evolution of bacterial CO 2 -concentrating mechanisms

Cyanobacteria rely on CO 2 -concentrating mechanisms (CCMs) to grow in today’s atmosphere (0.04% CO 2 ). These complex physiological adaptations require ≈15 genes to produce two types of protein complexes: inorganic carbon (Ci) transporters and 100+ nm carboxysome compartments that encapsulate rubisco with a carbonic anhydrase (CA) enzyme. Mutations disrupting any of these genes prohibit growth in ambient air. If any plausible ancestral form—i.e., lacking a single gene—cannot grow, how did the CCM evolve? Here, we test the hypothesis that evolution of the bacterial CCM was “catalyzed” by historically high CO 2 levels that decreased over geologic time. Using an E. coli reconstitution of a bacterial CCM, we constructed strains lacking one or more CCM components and evaluated their growth across CO 2 concentrations. We expected these experiments to demonstrate the importance of the carboxysome. Instead, we found that partial CCMs expressing CA or Ci uptake genes grew better than controls in intermediate CO 2 levels (≈1%) and observed similar phenotypes in two autotrophic bacteria, Halothiobacillus neapolitanus and Cupriavidus necator . To understand how CA and Ci uptake improve growth, we model autotrophy as colimited by CO 2 and HCO 3 − , as both are required to produce biomass. Our experiments and model delineated a viable trajectory for CCM evolution where decreasing atmospheric CO 2 induces an HCO 3 − deficiency that is alleviated by acquisition of CA or Ci uptake, thereby enabling the emergence of a modern CCM. This work underscores the importance of considering physiology and environmental context when studying the evolution of biological complexity.

Flamholz, Avi I.↗
DOE OSTI2024

A minimal complex of KHNYN and zinc-finger antiviral protein binds and degrades single-stranded RNA

Detecting viral infection is a key role of the innate immune system. The genomes of some RNA viruses have a high CpG dinucleotide content relative to most vertebrate cell RNAs, making CpGs a molecular marker of infection. The human zinc-finger antiviral protein (ZAP) recognizes CpG, mediates clearance of the foreign CpG-rich RNA, and causes attenuation of CpG-rich RNA viruses. While ZAP binds RNA, it lacks enzymatic activity that might be responsible for RNA degradation and thus requires interacting cofactors for its function. One of these cofactors, KHNYN, has a predicted nuclease domain. Using biochemical approaches, we found that the KHNYN NYN domain is a single-stranded RNA ribonuclease that does not have sequence specificity and digests RNA with or without CpG dinucleotides equivalently in vitro. We show that unlike most KH domains, the KHNYN KH domain does not bind RNA. Indeed, a crystal structure of the KH region revealed a double-KH domain with a negatively charged surface that accounts for the lack of RNA binding. Rather, the KHNYN C-terminal domain (CTD) interacts with the ZAP RNA-binding domain (RBD) to provide target RNA specificity. We define a minimal complex composed of the ZAP RBD and the KHNYN NYN-CTD and use a fluorescence polarization assay to propose a model for how this complex interacts with a CpG dinucleotide-containing RNA. In the context of the cell, this module would represent the minimum ZAP and KHNYN domains required for CpG-recognition and ribonuclease activity essential for attenuation of viruses with clusters of CpG dinucleotides.

Yeoh, Zoe C. (ORCID:0000000226949068)↗
DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗
DOE OSTI2025

Baryon number violation: from nuclear matrix elements to BSM physics

Processes that violate baryon number, most notably proton decay and $n\bar{n}$ transitions, are promising probes of physics beyond the Standard Model (BSM) needed to understand the lack of antimatter in the Universe. To interpret current and forthcoming experimental limits, theory input from nuclear matrix elements to UV complete models enters. Thus, an interplay of experiment, effective field theory, lattice QCD, and BSM model building is required to develop strategies to accurately extract information from current and future data and maximize the impact and sensitivity of next-generation experiments. Here, we briefly summarize the main results and discussions from the workshop ‘INT-25-91W: Baryon Number Violation: From Nuclear Matrix Elements to BSM Physics,’ held at the Institute for Nuclear Theory, University of Washington, Seattle, WA, 13–17 January 2025.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Variation in supercurrent connectivity and vortex pinning of state-of-the-art pulsed laser deposited RE Ba 2 Cu 3 O 7- x coated conductors

REBa 2 Cu 3 O 7-x (REBCO, RE: rare earth) coated conductors (CCs) suffer from great critical current I c differences between different manufacturers, I c variations within individual manufacturers, and often significant lengthwise fluctuations. The understanding of such variations is complicated by the lack of a direct correlation between I c and the critical current density J c . In fact, although J c is the fundamental property determined by the local vortex pinning landscape, I c is often limited by variable current blocking mechanisms. An important practical complexity is commercial practices of performing J c and I c evaluations based on I c at 77 K and self-field (sf), where connectivity variation dominates over vortex pinning variations. However, at higher fields and lower temperatures, vortex pinning becomes more complex and highly variable, making predictions of I c and J c at arbitrary temperature T, magnetic field H, and field orientation θ, quite uncertain. To address some aspects of this problem, we conducted detailed spool-to-spool performance characterization on recently manufactured REBCO CCs. Despite J c (77 K, sf) varies by only ∼11%, J c (77 K, 1 T) of its minimum and maximum (for H//ab-plane) show variations of ∼21% and ∼32%, respectively. This emphasizes the shortcoming in using J c (77 K, sf) as parameter for evaluating even the low field performance. An even more remarkable spool-to-spool J c variation of ∼68% was observed at 20 K and 15 T for H//ab-plane. To identify the origin of such lack of reproducibility we performed microstructural characterizations, which revealed, within the REBCO layer, large variation in the density of copper oxide (CuO x ) particles ranging from 0.1 to 2 μm in size. We believe that they play a decisive role in reducing the effective cross-section of the REBCO layer by not simply blocking current themselves, but also by nucleating off-axis REBCO grains, whose misoriented grain boundaries adversely impact REBCO grain-to-grain connectivity. The REBCO growth associated with high density CuO x particles also leads to the more disordered spatial arrangement of BaHfO 3 precipitate arrays, which, when self-aligning along the ab-planes, generate stronger pinning enhancing J c (H//ab) at all temperatures. In this way, we established that variations of both connectivity and vortex pinning are thus directly coupled. Our results also explain why the so-called ‘lift-factor’, typically defined by the ratio between I c (T, H) and I c (77 K, sf), frequently turns out to be unreliable.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Connecting primordial gravitational waves and dark energy

Cosmic acceleration manifested in the early universe as inflation, generating primordial gravitational waves detectable in the cosmic microwave background (CMB) radiation. Cosmic acceleration is occurring again at present as dark energy, detectable in cosmic distance and structure surveys. We explore the intriguing idea of connecting the two occurrences through quintessential inflation by an α-attractor potential without a cosmological constant. For this model we demonstrate robustness of the connection 1 + w 0 ≈ 4/(3N 2 r) between the present day dark energy equation of state parameter w 0 and the primordial tensor to scalar ratio r for a wide range of initial conditions. Analytic and numerical solutions produce current thawing behavior, resulting in a tight relation w a ≈ -1.53(1 + w 0 )≈ -0.2 (4 × 10 -3 /r). Upcoming CMB and galaxy redshift surveys can test this consistency condition. In conclusion, within this model, lack of detection of a dark energy deviation from Λ predicts a higher r, and lack of detection of r predicts greater dark energy dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Toroidal modeling of runaway avalanche in DIII-D discharges

Abstract A toroidal modeling tool is developed to study the runaway electron (RE) avalanche production process in tokamak plasmas, by coupling the Rosenbluth–Putvinski avalanche model (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355) with an n = 0 magneto-hydrodynamic (MHD) solver. Initial value numerical simulations are carried out for two DIII-D discharges with different plasma shapes (one near circular, and the other with high elongation). It is found that, assuming the same level of about 1% seed current level, the Rosenbluth–Putvinski model somewhat underestimates the RE plateau current for the circular-shaped plasma, as compared with that measured in DIII-D experiments. For an elongated, higher current plasma, simulations find strong runaway current avalanche production despite the lack of measured plateau RE current in experiments. A possible reason for this discrepancy is a lack of additional RE dissipation physics in the present two-dimensional model. Systematic scans of the plasma boundary shape, at fixed pre-disruption plasma current, find that the plasma elongation helps to reduce the RE avalanche production, confirming recent results obtained with an analytic model (Fülöp et al 2020 J. Plasma Phys. 86 474860101). The effect of the plasma triangularity (either positive or negative), on the other hand, has a minor effect. On the physics side, the avalanche process involves two competing mechanisms associated with the electric field. On the one hand, a stronger electric field produces a higher instantaneous avalanche growth rate. On the other hand, a fast growing RE current quickly reduces the fraction of the conduction current together with the electric field, and hence a faster dissipation of the poloidal flux. As a final result of these two dynamic processes, the runaway plateau current is not always the largest with the strongest initial electric field. These results lay the foundation for future self-consistent inclusion of the MHD dynamics and the RE amplification process.

Physics↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗