Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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

Enhancement of Distribution System State Estimation Using Pruned Physics-Aware Neural Networks: Preprint

Realizing complete observability in the three-phase distribution system remains a challenge that hinders the implementation of classical state estimation algorithms. In this paper, a new method so-called pruned physics-aware neural network (P2N2) is developed to improve the voltage estimation accuracy in the distribution system. The method relies on the physical grid topology, which is used to design the connections between different hidden layers of a neural network model. To verify the proposed method, a numerical simulation based on one-year smart meter data of load consumptions for threephase power flow is developed to generate the measurement and voltage state data. The IEEE 123 node system is selected as the test network to benchmark the proposed algorithm against the classical weighted least squares (WLS). Numerical results show that P2N2 outperforms WLS, in terms of data redundancy and estimation accuracy.

distribution systems state estimation↗
DOE OSTI2020

Reaching 1.5 and 2.0 °C global surface temperature targets using stratospheric aerosol geoengineering

A new set of stratospheric aerosol geoengineering (SAG) model experiments has been performed with Community Earth System Model version 2 (CESM2) with the Whole Atmosphere Community Climate Model (WACCM6) that are based on the Coupled Model Intercomparison Project phase 6 (CMIP6) overshoot scenario (SSP5-34-OS) as a baseline scenario to limit global warming to 1.5 or 2.0 °C above 1850–1900 conditions. The overshoot scenario allows us to applying a peak-shaving scenario that reduces the needed duration and amount of SAG application compared to a high forcing scenario. In addition, a feedback algorithm identifies the needed amount of sulfur dioxide injections in the stratosphere at four pre-defined latitudes, 30 °N, 15 °N, 15 °S, and 30 °S, to reach three surface temperature targets: global mean temperature, and interhemispheric and pole-to-Equator temperature gradients. These targets further help to reduce side effects, including overcooling in the tropics, warming of high latitudes, and large shifts in precipitation patterns. These experiments are therefore relevant for investigating the impacts on society and ecosystems. Comparisons to SAG simulations based on a high emission pathway baseline scenario (SSP5-85) are also performed to investigate the dependency of impacts using different injection amounts to offset surface warming by SAG. We find that changes from present-day conditions around 2020 in some variables depend strongly on the defined temperature target (1.5 °C vs. 2.0 °C). These include surface air temperature and related impacts, the Atlantic Meridional Overturning Circulation, which impacts ocean net primary productivity, and changes in ice sheet surface mass balance, which impacts sea level rise. Others, including global precipitation changes and the recovery of the Antarctic ozone hole, depend strongly on the amount of SAG application. Furthermore, land net primary productivity as well as ocean acidification depend mostly on the global atmospheric CO 2 concentration and therefore the baseline scenario. Multi-model comparisons of experiments that include strong mitigation and carbon dioxide removal with some SAG application are proposed to assess the robustness of impacts on societies and ecosystems.

42 ENGINEERING↗
DOE OSTI2022

Size-Controlled Nanoparticles Embedded in a Mesoporous Architecture Leading to Efficient and Selective Hydrogenolysis of Polyolefins

A catalytic architecture, comprising a mesoporous silica shell surrounding platinum nanoparticles (NPs) supported on a solid silica sphere (mSiO2/Pt-X/SiO2; X is the mean NP diameter), catalyzes hydrogenolysis of melt-phase polyethylene (PE) into a narrow C23-centered distribution of hydrocarbons in high yield using very low Pt loadings (~10 -5 g Pt/g PE). During catalysis, a polymer chain enters a pore and contacts a Pt NP where the C–C bond cleavage occurs and then the smaller fragment exits the pore. mSiO 2 /Pt/SiO 2 resists sintering or leaching of Pt and provides high yields of liquids; however, many structural and chemical effects on catalysis are not yet resolved. Here, we report the effects of Pt NP size on activity and selectivity in PE hydrogenolysis. Time-dependent conversion and yields and a lumped kinetics model based on the competitive adsorption of long vs short chains reveal that the activity of catalytic material is highest with the smallest NPs, consistent with a structure-sensitive reaction. Remarkably, the three mSiO 2 /Pt-X/SiO 2 catalysts give equivalent selectivity. We propose that mesoscale pores in the catalytic architecture template the C 23 -centered distribution, whereas the active Pt sites influence the carbon–carbon bond cleavage rate. This conclusion provides a framework for catalyst design by separating the C–C bond cleavage activity at catalytic sites from selectivity for chain lengths of the products influenced by the structure of the catalytic architecture. The increased activity, selectivity, efficiency, and lifetime obtained using this architecture highlight the benefits of localized and confined environments for isolated catalytic particles under condensed-phase reaction conditions.

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

Microstructure-based prediction of thermal aging strength reduction factors for grade 91 ferritic-martensitic steel

This study aimed to develop a microstructure-based mechanistic approach to address the long-term thermal aging effect on yield stress and ultimate tensile strength and to provide a physical basis for developing thermal aging factors for G91 for a design life of 60 years. Several heats of G91 steel were examined. Controlled aging experiments were conducted on two heats at temperatures of 550, 600, and 650 °C for times up to ~64,000 h. Specimens were also taken from archived creep-tested specimens of G91 and from a tube removed from the Kingston coal-fired power plant to obtain data that cover a wide range of temperatures and for aging times up to 155,000 h. Thermal aging caused significant subgrain recovery, coarsening of M 23 C 6 carbides at subgrain boundaries and MX carbonitrides within subgrains. The intermetallic Laves phase forms during aging and grows rapidly. The growth rate of the subgrain width and MX mean size during thermal aging were described by kinetic models. Thermal aging results in the reductions in the yield stress and the ultimate tensile strength of G91. The effects of thermal aging on the reductions in yield stress and ultimate tensile strength were well described by the microstructure-strength model that considers three superimposed strengthening mechanisms, namely, sub-boundary strengthening, MX precipitation hardening, and Mo solid solution strengthening. Furthermore, the model was independently validated by the ASME Code values and is being used by the ASME to develop the aging-induced strength reduction factors for G91 steel for a design life of 60 years.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modified Energy Span Analysis Reveals Heterogeneous Catalytic Kinetics

Mechanistic modeling is a cornerstone of catalyst development generally conducted with microkinetic models or density functional theory-based energy profiles. We extend the energy span model of homogeneous catalysis to heterogeneous systems by introducing the modified energy span analysis (MESA) model by implementing collision theory and gas-phase concentration effects. We determine analytically turnover frequencies, coverages, rate-determining steps, apparent activation energies, and reaction orders in agreement with microkinetic and kinetic Monte Carlo simulations. The model applies to discrete catalyst state systems, including single-atom catalysts, homogeneous systems, and microporous materials. A generalizable rate expression is derived for reactor modeling or mechanistic insights from solely experimentally measured reaction orders. Here, we illustrate MESA on three published mechanisms and reveal unexpected phenomena absent in traditional energy span profiles.

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

Evaluating Possible Formation Mechanisms of Criegee Intermediates during the Heterogeneous Autoxidation of Squalene

Organic molecules in the environment oxidatively degrade by a variety of free radical, microbial, and biogeochemical pathways. A significant pathway is heterogeneous autoxidation, in which degradation occurs via a network of carbon and oxygen centered free radicals. Recently, we found evidence for a new heterogeneous autoxidation mechanism of squalene that is initiated by hydroxyl (OH) radical addition to a carbon–carbon double bond and apparently propagated through pathways involving Criegee Intermediates (CI) produced from β-hydroxy peroxy radicals (β-OH-RO 2 •). It remains unclear, however, exactly how CI are formed from β-OH-RO 2 •, which could occur by a unimolecular or bimolecular pathway. Combining kinetic models and multiphase OH oxidation measurements of squalene, we evaluate the kinetic viability of three mechanistic scenarios. Scenario 1 assumes that CI are formed by the unimolecular bond scission of β-OH-RO 2 •, whereas Scenarios 2 and 3 test bimolecular pathways of β-OH-RO 2 • to yield CI. Scenario 1 best replicates the entire experimental data set, which includes effective uptake coefficients vs [OH] as well as the formation kinetics of the major products (i.e., aldehydes and secondary ozonides). Finally, although the unimolecular pathway appears to be kinetically viable, future high-level theory is needed to fully explain the mechanistic relationship between CI and β-OH-RO 2 • in the condensed phase.

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

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗
DOE OSTI2026

Revisiting the three-kaon interaction and its relation with 𝐾⁡(1460)

Here, we test the hypothesis of 𝐾⁡(1460) being a hadronic 𝐽 𝑃 = 0 − molecule through nonperturbative 𝑆-wave $KK\bar{⁢K}$, 𝐾⁢𝜋⁢𝜋, 𝐾⁢𝜋⁢𝜂 coupled-channel dynamics in a three-body unitary isobar approach. The scalar two-body coupled-channel resonances 𝑓 0 ⁡(500), 𝑓 0 ⁡(980), 𝑎 0 ⁡(980), and 𝐾$^*_0$⁡(700) are generated in the subsytems in amplitudes that match phase shifts from experiment. In a first step, previous results in the limit of mass-degenerate, stable 𝑓 0 and 𝑎 0 isobars are reproduced. Once the full seven-coupled channel model is switched on, other 𝑆-matrix effects obscure and modify the resonance signal, including complex thresholds, a two-body cusp, and a triangle singularity at threshold.

Döring, M. [George Washington Univ., Washington, D↗
DOE OSTI2021

Origin of the magnetic and orbital ordering in α-Sr 2 CrO 4

Motivated by recent experimental progress in transition metal oxides with the K 2 NiF 4 structure, in this study we investigate the magnetic and orbital ordering in α-Sr 2 CrO 4 . Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the ab initio band structure near the Fermi level. The unique reverse splitting of t 2g orbitals in α-Sr 2 CrO 4 , with the 3d 2 electronic configuration for the Cr 4+ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground state phase diagram for the two dimensional compound α-Sr 2 CrO 4 . We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of α-Sr 2 CrO 4 , we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound α-Sr 2 CrO 4 , in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with d xy driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the d xz and d yz orbitals driving the orbital order.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Genetic Basis of Chromate Adaptation and the Role of the Pre-existing Genetic Divergence during an Experimental Evolution Study with Desulfovibrio vulgaris Populations

Hexavalent chromium [Cr(VI)] is a common environmental pollutant. However, little is known about the genetic basis of microbial evolution under Cr(VI) stress and the influence of the prior evolution histories on the subsequent evolution under Cr(VI) stress. In this study, Desulfovibrio vulgaris Hildenborough (DvH), a model sulfate-reducing bacterium, was experimentally evolved for 600 generations. By evolving the replicate populations of three genetically diverse DvH clones, including ancestor (AN, without prior experimental evolution history), non-stress-evolved EC3-10, and salt stress-evolved ES9-11, the contributions of adaptation, chance, and pre-existing genetic divergence to the evolution under Cr(VI) stress were able to be dissected. Significantly decreased lag phases under Cr(VI) stress were observed in most evolved populations, while increased Cr(VI) reduction rates were primarily observed in populations evolved from EC3-10 and ES9-11. The pre-existing genetic divergence in the starting clones showed strong influences on the changes in lag phases, growth rates, and Cr(VI) reduction rates. Additionally, the genomic mutation spectra in populations evolved from different starting clones were significantly different. A total of 14 newly mutated genes obtained mutations in at least two evolved populations, suggesting their importance in Cr(VI) adaptation. An in-frame deletion mutation of one of these genes, the chromate transporter gene DVU0426, demonstrated that it played an important role in Cr(VI) tolerance. Overall, our study identified potential key functional genes for Cr(VI) tolerance and demonstrated the important role of pre-existing genetic divergence in evolution under Cr(VI) stress conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Light element ($$Z=1,2$$) production from spontaneous ternary fission of $$^{252}$$Cf

Abstract The yields of light elements ( $$Z=1,2$$ Z = 1 , 2 ) obtained from spontaneous ternary fission of $$^{252}$$ 252 Cf are treated within a nonequilibrium approach, and the contribution of unstable nuclei and excited bound states is taken into account. These light cluster yields may be used to probe dense matter, and to infer in-medium corrections such as Pauli blocking which is determined by the nucleon density. Continuum correlations are calculated from scattering phase shifts using the Beth-Uhlenbeck formula, and the effect of medium modification is estimated. The relevant distribution is reconstructed from the measured yields of isotopes. This describes the state of the nucleon system at scission and cluster formation, using only three Lagrange parameters which are the nonequilibrium counterparts of the temperature and chemical potentials, as defined in thermodynamic equilibrium. We concluded that a simple nuclear statistical equilibrium model neglecting continuum correlations and medium effects is not able to describe the measured distribution of H and He isotopes. Moreover, the freeze-out concept may serve as an important ingredient to the nonequilibrium approach using the relevant statistical operator concept.

Röpke, G.↗
DOE OSTI2023

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Abstract Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Environmental Sciences & Ecology↗
DOE OSTI2020

TRANSENSOR: Transformer Real-time Assessment INtelligent System with Embedded Network of Sensors and Optical Readout. Final Report

Utilities across the world are wrestling with evolving market dynamics, population growth and climate change. Distributed Energy Resources (DERs) such as solar photovoltaics, distributed generators and energy storage systems are becoming important parts of the U.S. energy mix. These factors are driving an increasing need for low-cost grid asset monitoring. Aside from being costly, traditional utility monitoring systems are not sufficiently robust and do not provide real-time visibility into the condition of grid assets such as transformers. Lack of accurate real-time measurements on performance has resulted in the use of lagging indicators, such as oil sample analysis. To address this need, an innovative embedded optical sensing technology, Transformer Real-time Assessment Intelligent System (TRANSENSOR) was developed, validated and demonstrated in this project. To date, TRANSENSOR has focused on transformers but is extendable to other grid assets. In addition to the technology being embedded into new transformers during manufacturing, a retrofit configuration that can be installed on existing transformers in the field was also developed. It is anticipated that the ability to retrofit existing transformers will help accelerate adoption of the underlying TRANSENSOR technology. Phase 1 of the project focused on laboratory development and qualification of the technology. Following iterations and exploration of relevant optical sensing modalities and multiplexed configurations of interest for the transformer environment, an effective candidate configuration was agreed upon, down-selected and custom-designed for embedding into General Electric (GE) network transformers. Two new GE 500 kilovolt ampere (kVA), 27 kilovolt (kV) distribution network transformers were built with embedded fiber-optic (FO) sensors and successfully qualified per industry standards at GE’s Shreveport, Louisiana facility during Phase 1 of the project. Following the successful completion of Phase 1, the team proceeded to Phase 2, which focused on a field demonstration of the technology. Over Phase 2, the first new GE transformer built with embedded fiber-optic sensors was installed in an above-ground cage and connected to the grid at Con Edison’s Astoria facility. A second transformer equipped with TRANSENSOR was installed in an underground vault. Additionally, an older (1982 year model) GE transformer in an above-ground cage was retrofitted with fiber- optic sensors and reconnected to the grid at ConEd’s facility. Analysis of data acquired from the sensors showed interesting correlations with transformer loading. Additionally, key events such as transformer low-voltage network connection, primary-side energizing, and a pressure loss event from an oil sampling were detected by the TRANSENSOR system. Online data processing/feature extraction algorithms for the second transformer in the underground vault detected key features and event alerts that were transmitted through a wireless 4G connection. The remote deployment concept showed promising results with data collection running for a total of 8 months across the three (3) GE transformers instrumented for the Phase 2 field demonstration. This sets the stage well for further development and commercialization.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Bidentate N‐based Ligands for Highly Reusable, Ligand‐coordinated, Supported Pt Hydrosilylation Catalysts

Abstract A significant challenge in designing supported metal‐ligand catalysts for solution‐phase reactions is the stabilization of the metal active sites against leaching into solution. Here, we examine alkene hydrosilylation reactions as model systems to improve the stability of highly dispersed Pt using a metal‐ligand coordination strategy on high surface area oxide supports. By evaluating a series of bidentate N‐based ligands, we demonstrate several design strategies to improve stability of the highly dispersed Pt 2+ centers against leaching, while maintaining a high level of catalytic activity, selectivity, and recyclability for alkene hydrosilylation batch reactions under mild conditions. These involve a bi‐functional approach to ligand design, which considers interaction to the support and a well‐defined coordination environment for the metal active site. Three strategies are reported: modifying ligands for stronger interaction with oxide surfaces, mixing ligands, and pre‐depositing an “anchoring ligand” to the support before loading the metal‐ligand catalyst. Each of these is successful in enhancing Pt recyclability. Particularly, two Pt‐phenanthroline catalysts exhibit excellent reusability for multiple batch reaction cycles, due to high stability of the active Pt species. Addressing the active site leaching problem significantly enhances the utility of ligand‐coordinated supported metal catalysts as highly stable and selective catalysts for solution‐phase reactions.

Chen, Linxiao↗
DOE OSTI2022

OC6 Phase Ia Definition Document: Validation of Nonlinear Hydrodynamic Loading on the DeepCwind Semisubmersible

The objective of OC6 Phase I is to investigate the persistent under-prediction (about 20% on average) of the structural loads in the OC5-DeepCwind floating semisubmersible wind system, as observed within the previous OC5 Phase II project. The OC5 results indicate that much of this under-prediction originates from the low-frequency response of the system at the surge/pitch natural frequencies, resulting from nonlinear hydrodynamic loading and/or overprediction of hydrodynamic damping. This work package therefore focuses on better understanding this low-frequency response behavior, the applicability of computational-fluid-dynamics (CFD) models and engineering-level hydrodynamic models for predicting the nonlinear hydrodynamic loading, and the level of uncertainty in the measured response characteristics. To address these objectives, two experimental campaigns were developed and performed in 2017 and 2018 by a sub-group of the OC5 project: Floating test campaign: A simplified OC5-DeepCwind semisubmersible model (no turbine, rigid tower) was moored to three taut-spring lines and subjected to motion decay and wave tests (regular and irregular waves), during which motion amplitudes of the model were measured (see Figure 1a). The purpose of this test campaign was to examine the repeatability of the low-frequency response of the system to wave excitation and determine the level of uncertainty in the response behavior.Constrained test campaign: The same model, but without the tower present, was fixed to a carriage and subjected to towing, forced oscillation and wave tests (regular and irregular waves), during which forces exerted upon the model were measured (see Figure 1b). The purpose of this test campaign was to examine the hydrodynamic loading contributions by separating the wave excitation and radiation forces. Both model test campaigns were performed in the concept basin of the Maritime Research Institute Netherlands (MARIN) in Wageningen, the Netherlands and the first test was within the framework of the MARINET2 project. Both experimental campaigns were examined within Phase Ia of the OC6 project. The purpose of this report is to supply the needed information to participants of the OC6 project to build and simulate the experimental campaigns.

17 WIND ENERGY↗
DOE OSTI2020

Deep Reinforcement Learning Based Volt-VAR Optimization in Smart Distribution Systems

This paper develops a model-free volt-VAR optimization (VVO) algorithm via multi-agent deep reinforcement learning (DRL) in unbalanced distribution systems. This method is novel since we cast the VVO problem in distribution networks to an intelligent deep Q-network (DQN) framework, which avoids solving a specific optimization model directly when facing time-varying operating conditions in the systems. We consider statuses/ratios of switchable capacitors, voltage regulators, and smart inverters installed at distributed generators as the action variables of the agents. A delicately designed reward function guides these agents to interact with the distribution system, in the direction of reinforcing voltage regulation and power loss reduction simultaneously. The forward-backward sweep method for radial three-phase distribution systems provides accurate power flow results within a few iterations to the DRL environment. The proposed method realizes the dual goals for VVO. We test this algorithm on the unbalanced IEEE 13-bus and 123-bus systems. Numerical simulations validate the excellent performance of this method in voltage regulation and power loss reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Towards a New μ→eγ Search with the MEG II Experiment: From Design to Commissioning

The MEG experiment represents the state of the art in the search for the Charged Lepton Flavour Violating μ + →e + γ decay. With its first phase of operations at the Paul Scherrer Institut (PSI), MEG set the most stringent upper limit on the BR (μ + →e + γ) ≤ 4.2×10 –13 at 90% confidence level, imposing one of the tightest constraints on models predicting LFV-enhancements through new physics beyond the Standard Model. An upgrade of the MEG experiment, MEG II, was designed and it is presently in the commissioning phase, aiming at a sensitivity level of 6×10 –14 . The MEG II experiment relies on a series of upgrades, which include an improvement of the photon detector resolutions, brand new detectors on the positron side with better acceptance, efficiency and performances and new and optimized trigger and DAQ electronics to exploit a muon beam intensity twice as high as that of MEG (7×10 7 μ + /s). This paper presents a complete overview of the MEG II experimental apparatus and the current status of the detector commissioning in view of the physics data taking in the upcoming three years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗