Search NASA⌕ Search

SEARCH · Search NASA

Results for “heat simulation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling of thermal pressurization in tight claystone using sequential THM coupling: Benchmarking and validation against in-situ heating experiments in COx claystone

We apply thermoporoelasticity and a sequentially coupling technique for modeling thermally-driven coupled Thermo-Hydro-Mechanical (THM) processes in tight claystone. A THM benchmark case with a corresponding analytic solution for thermoporoelasticity under a constant heat loading verifies the model. Thereafter, two in situ heating experiments are simulated for model validation: a smaller-scale heating experiment (TED experiment) and a larger-scale experiment (ALC experiment) in Callovo-Oxfordian (COx) claystone at the Meuse/Haute-Marne underground research laboratory in France. The model exhibits good performance to match the observed temperature and pore pressure evolution for the smaller-scale TED experiment. For the larger-scale ALC experiment, general trends of thermal-pressurization are captured in the modeling, but pressure is underestimated at some monitoring points during cool-down. This indicates that the THM response in the field may be affected by the variability of rock's properties or irreversible or time-dependent mechanical processes that are not included in the current thermoporoelastic model. The main contributions of this work are as follows: (1) we verify and validate the numerical simulator, TOUGH-FLAC, to be a valuable coupled THM modeling tool; (2) prove that the laboratory determined material parameters can be used as reference values for upscaling experiments. However, to better identify and quantify THM processes with modeling of in situ tests, more emphasize should be dedicated to obtaining high-quality mechanical deformation data.

58 GEOSCIENCES↗
DOE OSTI2020

Benchmark Simulation of the Natural Convection Shutdown Heat Removal Test Facility Using SAM

Natural convection systems are a promising method to passively remove heat from reactor cavities during loss of forced flow accident scenarios. In this work at Argonne National Laboratory (ANL), a highly instrumented Natural Convection Shutdown Heat Removal Test Facility (NSTF) was used to demonstrate the effectiveness of air-cooled natural convection systems. In previous work, RELAP5-3D simulations were performed on this facility with favorable comparisons to experiment for mass flow rate, pressure drop, air temperature increase, and air velocity. Both experimental and simulation efforts with this facility present a useful opportunity to perform a benchmark study with the System Analysis Module (SAM). SAM is an advanced thermal-hydraulic system code currently in development at ANL for advanced non-light water reactor safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Potts Model parameter study of particle size, Monte Carlo temperature, and “Particle-Assisted Abnormal Grain Growth”

A Potts Model was proposed to account for temperature and particle heterogeneity-dependent Zener Pinning. This was accomplished by relating the random fluctuations of grain boundary position allowed by the Potts Model switching probability at finite simulation temperature to experimentally observed grain growth stagnation behavior. We assume that these fluctuations arise from random fluctuations in the thermodynamic and kinetic properties of the grain boundaries and/or interfaces as they change with temperature. As an application of this model, the grain growth kinetics of U-10 wt. % Mo nuclear fuels were simulated with the input of microstructural images during different heat treatment processes. Simulated average grain growth behavior is in good agreement with experiments.

Frazier, William E.↗
DOE OSTI2022

The Hanford Lead Canister Collaboration

The Hanford Site, Management of the Cesium and Strontium Capsules (MCSC) Project is preparing to place cesium and strontium material into dry storage using stainless steel canisters that are housed inside vertical concrete casks. Part of the site’s aging management program is to use a spare canister system as a leading indicator canister to provide advance warning for signs of chloride induced stress corrosion cracking (CISCC) or any other kind of canister degradation that might take place over the system’s 300-year design life. The Hanford Lead Canister (HLC) will use electric heaters to simulate the decay heat of cesium and strontium to ensure the HLC simulates the environment and conditions the capsule-loaded systems will experience. In addition to being part of the site’s aging management program, the HLC is being made available for research and development that advances the industry regarding canister mitigation and repair (M&R). The HLC system closely resembles current spent nuclear fuel (SNF) dry storage systems, making the HLC useful for demonstrating technologies or performing research that is relevant to SNF canisters . Several organizations are collaborating to prepare the HLC, along with technical strategy, to perform the needed research and development activities, to include the Pacific Northwest National Laboratory (PNNL), Central Plateau Cleanup Company (CPCCo), and the Electric Power Research Institute (EPRI). In one example of this collaborative research, EPRI is planning to use the HLC to demonstrate in-situ M&R technologies under controlled conditions before the HLC is deployed at Hanford. There is a period of several years when the HLC will be available for canister M&R activity before it is deployed. This paper describes the development of the HLC, its role in the Hanford site’s MCSC Project aging management program, and some of the research activities that are planned for the HLC.

Klymyshyn, Nicholas A.↗
DOE OSTI2021

A fast reduced model for a shell-and-tube based latent heat thermal energy storage heat exchanger and its application for cost optimal design by nonlinear programming

Numerical simulation of latent heat thermal energy storage (LHTES) systems plays a fundamental role in studying the physical process and guiding the engineering design. Discretization of the PDEs describing the nonlinear solidification/melting process of phase change materials (PCMs) leads to a large-scale complex dynamics system, where the system behavior depends on a set of parameters. In a design setting, repeated model evaluations are required over the set of parameters results in significant computational burden. In this paper, an explicit analytic solution was built for the propagation of the solidification front in a cylindrical coordinate. The analytic solution approach is further employed to develop a low computational reduced model (RM) as a module for a shell-and-tube based LHTES heat exchanger. The levelized Cost of Energy (LCOE) is used as a design metric and the RM model is used to apply system-level constraints in the nonlinear programming formulation that facilitates efficient global optimal design of the PCM properties, flow conditions and tube geometries. The use of LCOE as the design metric prevents over design of the heat transfer rate and also establishes a fair ground for evaluation of different thermal storage technologies and their integrated applications with other systems. Optimal results showed that a higher effectiveness results in a higher LCOE; the velocity of the HTF and the length of the channel are highly correlated with each other; both larger PCM latent energy and conductivity result in lower LCOE.

25 ENERGY STORAGE↗
DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗
DOE OSTI2021

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tusqh

SAND2025-00675O Tusqh is a software tool that generates cubical meshes in 2D and 3D and computes the homology of these meshes using persistent homology. It includes a grid cell in the output if its volume-fraction is above a selectable threshold, estimated by sampling points within the cell. Tusqh incorporates anti-aliasing algorithms to mitigate grid orientation and scale effects. It is designed for creating finite element meshes for simulations and can be used in various applications such as heat diffusion, mechanical simulations, and computer graphics rendering. The software outputs meshes in an open format compatible with downstream software. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2025

Anisotropic Heating and Parallel Heat Flux in Electron-only Magnetic Reconnection with Intense Guide Fields

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Analysis of HolosGen Sub-Scale Simulator with Plant Dynamics Code

The Subcritical Power Module Sub-scale Simulator (SPM-SS) has been designed and constructed by HolosGen LLC under the ARPA-E MEITNER program to simulate the thermal-hydraulic and heat transfer behavior of the full-scale Holos-Quad Subcritical Power Modules (SPMs). Four coupled SPMs, each rated at 5.5MW, form the Holos-Quad gas-cooled microreactor design. The SPM-SS represents a substantially scaled-down system with a power rating less than 40 kW, equipped with an electrically heated fuel cartridge heat exchanger, an electrically heated compressor heat exchanger, and a valve actuated turbine heat exchanger, in addition to a recuperator and a cooler heat exchanger. The fuel cartridge represents a portion of the full-scale SPM core, the compressor heat exchanger mimics the temperature changes resulting from the compressor’s turbomachinery inefficiencies, the turbine heat exchanger mimics the expansion process normally occurring through the turbine, while the recuperator and cooler heat exchangers complete the subscale simulator loop. The heaters equipping the fuel cartridge and the compressor heat exchangers are electronically controlled to simulate normal and off-normal SPM operating conditions. The full-scale Holos-Quad SPM design eliminates the traditional balance of plant and executes thermal-to-electric energy conversion by means of an intercooled Brayton cycle with decoupled compressor-turbine turbomachinery. The Holos-Quad full-scale design is equipped with a multi-stage axial Low- and High-Pressure compressor, and a multistage axial turbine. The SPM-SS is designed for testing and validation of selected components which are instead coupled by a traditional balance of plant. The SPMSS is not equipped with turbo-machinery (compressor and turbine) as the development of these components were excluded from the scope of work under the ARPA-E MEITNER funding program. The SPM-SS balance of plant enables modifications, replacement and testing of individual components with different working fluids and is designed to include the turbo-machinery components that will be developed in future research . The SPM-SS can be operated with different gases, variable mass-flow-rates, pressures, and temperatures to obtain test data for selected components, whose performance can be scaled to validate the computer model of the full-scale SPM at various conditions (e.g., start-up, transients conditions). The SPM-SS can operate at the maximum Holos-Quad design pressure of 7 MPa, and a maximum temperature limited to 650 °C by the electrical heaters. Several SPM-SS tests have been conducted and analyzed with the Plant Dynamics Code (PDC) developed at the Argonne National Laboratory (ANL). These tests aimed at validating the PDC modeled predictions of the full-scale Holos-Quad design with data from selected SPM-SS components. In order to address SPM-SS specific characteristics, such as components heat losses and absence of turbomachinery components, some modifications to the PDC have been implemented to factor the design differences from the full-scale Holos-Quad SPM to the SPM-SS. As the PDC offers capabilities to analyze systems with different working fluids, air, nitrogen, and helium were utilized as the SPM-SS working fluids. Air was utilized to fine-tune the SPM-SS Systems Structures and Components (SSCs), nitrogen was utilized to pressure test the SPM-SS loop at the SPMs design maximum pressure of 7MPa. Helium was utilized as the working fluid circulating through the SPM-SS SSCs for specific tests to validate the PDC predictions of the fuel cartridge heat exchanger. SPM-SS tests data were also analyzed with both the steady-state and transient analysis capabilities offered by the PDC. This report describes the PDC analysis of the SPM-SS tests data, including the necessary code modifications and comparison of the code results with the experimental data. Based on the results, a discussion is presented on how the analysis supports design and transient calculations of the full-scale Holos-Quad microreactor. Also based on the results of this work, recommendations are made for future optimizations of the SPM-SS components and PDC model development needs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Energy and Exergy Analysis of Low-Global Warming Potential Refrigerants as Replacement for R410A in Two-Speed Heat Pumps for Cold Climates

Heat pumps (HPs) are being developed with a new emphasis on cold climates. To lower the environmental impact of greenhouse gas (GHG) emissions, alternate low global warming potential (GWP) refrigerants must also replace the exclusive use of the refrigerant R410A, preferably without re-engineering the mechanical hardware. In this paper, we analyze the performance of four low-GWP alternative refrigerants (R32, R452B, R454B, and R466A) relative to the conventional R410A and draw conclusions on the relative performances for providing heating in cold climates based on the Air-Conditioning, Heating, and Refrigeration Institute (AHRI) 210/240 standard for two-speed heat pumps. The simulations are carried using the Department of Energy, Oak Ridge National Laboratory (DOE/ORNL) Heat Pump Design Model (HPDM), a well-known heating, ventilation, and air conditioning (HVAC) modeling and design tool in the public domain and the HVAC research and development community. The results of the simulation are further scrutinized using exergy analysis to identify sources of systemic inefficiency, the root cause of lost work. This rigorous approach provides an exhaustive analysis of alternate low-GWP refrigerants to replace R410A using available compressors and system components, without compromising performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Parameterizing Subgrid Variations of Land Surface Heat Fluxes to the Atmosphere Improves Boreal Summer Land Precipitation Simulation With the NCAR CESM1.2

Subgrid horizontal variations of land surface heat fluxes to the atmosphere resulting from subgrid land cover heterogeneity are important in land-atmosphere interaction in global climate models (GCMs). To incorporate it, a parameterization using stochastic sampling based on truncated normal distributions diagnosed from the land model and internal ensemble mean of multiple calls to the planetary boundary layer and deep convection schemes is developed. After its implementation in the National Center for Atmospheric Research (NCAR) Community Earth System Model version 1.2 (CESM1.2), large changes of simulated land precipitation are found to occur in regions where large subgrid variations of surface heat fluxes exist. The simulated precipitation is improved in boreal summer, especially over eastern China and the coastal areas of the Bay of Bengal. The improved precipitation is mainly a result of improved large-scale moisture convergence and advection by altered vertical diffusion and convection.

58 GEOSCIENCES↗
DOE OSTI2022

Simulating a solar parabolic trough collector plant used for industrial process heat using an optimized operating scheme that utilizes flexible heat integration

This study focuses on the modeling and simulation of a novel design and operation of a solar industrial process heat (SIPH) plant that uses a parabolic trough collector system for generating process heat. The SIPH plant incorporates flexible heat integration (FHI) by having two options for heat sinks as well as a flexible collection temperature for the heat transfer fluid. Leveraging the degrees of freedom created by FHI allows an optimized operating scheme to be created which maximizes the solar heat delivered to the industrial processes. Yearly results for a 27 MW-t plant located in Salt Lake City, UT are simulated for a base case and a case utilizing optimized FHI. For the main process temperature case presented, FHI increases the total solar share by 13.0% and reduces the levelized cost of heat of the SIPH plant by 10.4% relative to the base case. Furthermore, three case study days are presented which highlight how FHI maximizes solar share, specifically during mornings, afternoons, and days with intermittent solar conditions. An environmental analysis is also performed showing that the SIPH plant can reduce harmful emissions by 15.4% for the base case and 17.4% for the FHI case compared to a plant that only utilizes a natural gas furnace for process heat. Overall, the work thoroughly explores an effective operating scheme for a SIPH plant, helping to make it more economical and environmentally beneficial.

14 SOLAR ENERGY↗
DOE OSTI2025

Turbulence simulations of drift-Alfvén wave effects on edge-localized modes and divertor target heat flux width

A series of BOUT++ turbulence simulations are performed for two Pre-Fusion Power Operation (PFPO-1,2) phases of the ITER Research plan proclaimed in 2019. Linear simulation results show that PFPO-1 is unstable to Peeling-ballooning modes (PB), while the PFPO-2 is unstable to the coupling of PB and Drift-Alfvén wave (DAW) instabilities. The linear results are qualitatively consistent with the dispersion relation of DAW. To study the influence of the DAW on the edge-localized mode (ELM) crash, the PFPO-2 is adopted in this section. Different from the grassy ELM in [Li et al 2022 Nucl. Fusion 62, 096030], nonlinear simulations show that the ELM size of PFPO-2 is almost one-third of the grassy ELM, representing a distinct small ELM. However, simulations then show that if the PB instability is removed, the fluctuation amplitude drops by an order of magnitude and the ELM crash disappears, which is in accord with the theory in [Xu et al 2010 Physical Review Letters, Vol. 105, 175005] and the results in [Li et al 2022 Nucl. Fusion 62, 096030], confirming that the PB instability is a necessary condition for an ELM crash. Furthermore, removing the DAW drive also suppresses ELM crashes, implying that PB instability is necessary but insufficient for PFPO-2 ELM and that DAW could amplify PB-driven turbulence. In addition, by integrating results of both PFPO phases by BOUT++ turbulence analysis with heat flux width (λ q ) and electron thermal diffusivity (χ e ), these simulations establish a predictive framework for heat flux width and its dominant mechanisms, thereby providing a comprehensive scheme for the prediction of the heat flux width. Moreover, simulations indicate that DAW driving can increase the transport coefficient by enhancing the turbulent transport, leading to a broadened heat flux width once the transport coefficient exceeds

Drift-Alfvén wave↗