Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 1,045 records · Page 58

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗

HERO WEC V1.0 2024 - WEC-Sim Detailed Simulation Runs and Summary Data

This dataset includes results from simulations of NREL's hydraulic and electric reverse osmosis wave energy converter (HEREO WEC). Simulation runs include 135 wave cases that were based on the updated WEC-Sim model, which is linked below. The data represented in this repository is based on an updated WEC-Sim model using laboratory data to tune and refine the original WEC-Sim model for the V1.0 HERO WEC. The 135 wave cases represent waves with the following wave height and wave period ranges: - Significant Wave Height: 0.25 - 3.75m in 0.25m increments - Wave Period: 5 - 13 sec in 1 sec increments Each run was simulated using a Pierson-Moskowitz irregular wave spectrum with a 100 second ramp time, a total simulation time of 3,100 seconds, and a simulation time-step of 0.005s. A reference table has been included to map each multi condition run (MCR) case with each wave condition. Summary data set includes a spreadsheet and image files with matrices that are associated with data from simulation runs. All matrices cover the same significant wave height and wave periods from the simulation runs, in the same increments. The following matrices are included: - Power Abs: The average absorbed power from the WEC (calculated from anchor reaction force and heave velocity) - Power Hyd: The average hydraulic power output at pump (calculated from pump output flow and pressure) - Power - Hyd ROi: The average hydraulic power measured at the RO system inlet (calculated from RO system pressure and flow (pre-accumulator)) - Flow - Pump out: The average flowrate measured at the pump outlet - Flow - Perm: The average permeate (clean water) production - Flow - RO (pre): The average flowrate measured at the inlet of the RO system before the accumulators - Flow - RO (post): The average flowrate measured after the accumulator bank in the RO system - Pressure - RO: The average pressure measured at the inlet of the RO system This data set has been developed by the National Renewable Energy Laboratory, operated by the Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. Funding provided by the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Water Power Technologies Office.

16 TIDAL AND WAVE POWER↗

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hardware-in-the-Loop Simulation of a Distribution System with Air Conditioners under Model Predictive Control: Preprint

Many have proposed that responsive load provided by distributed energy resources (DERs) and demand response (DR) are an option to provide flexibility to the grid and especially to distribution feeders. However, because responsive load involves a complex interplay between tariffs and DER and DR technologies, it is challenging to test and evaluate options without negatively impacting customers. This paper describes a hardware-in-the-loop (HIL) simulation system that has been developed to reduce the cost of evaluating the impact of advanced controllers (e.g., model predictive controllers) and technologies (e.g., responsive appliances). The HIL simulation system combines large-scale software simulation with a small set of representative building equipment hardware. It is used to perform HIL simulation of a distribution feeder and the loads on it under various tariff structures. In the reported HIL simulation, loads include many simulated air conditioners and one physical air conditioner. Independent model predictive controllers manage operations of all air conditioners under a time-of-use tariff. Results from this HIL simulation and a discussion of future development work of the system are presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Early implementation of digital upgrades on simulators - 347

An important aspect of successful Nuclear Power Plant (NPP) digital control system upgrade projects is to upgrade the Operator Training Simulator (OTS) or a related Engineering Simulator (ES), prior to installation of the Distributed Control System (DCS) on the plant allowing for early testing of the digital controls/HMI and for early training of the plant operators while maintaining OTS configuration management (i.e. OTS in sync with the plant's current design). Experience has shown that early implementation of digital upgrades on the OTS and/or an ES can reduce plant startup duration (related to digital systems commissioning) and ensure operators are proficient with the new digital systems to ensure safe and reliable plant operation. Including simulators in the design process offers many benefits including early identification of design deficiencies, testing deficiency fixes, developing and validating plant procedures and supporting human factors engineering, to name a few. Implementing digital upgrades on the OTS is also a regulatory requirement driven by the ANSI/ANS-3.5-2009 standard (Nuclear Power Plant Simulators for Use in Operator Training and Examination) if the new DCS affects actions to be taken by the plant operator in the main control room. This paper presents the benefits of implementing digital upgrades on the OTS and/or ES prior to plant implementation. The various simulator implementation strategies (stimulation, emulation, simulation or a hybrid solution) and selection factors are discussed. The paper also covers four (4) case studies demonstrating how digital upgrades can be implemented on simulators. The first case study covers the implementation of a Siemens offgas Programmable Logic Controller (PLC) system on the Enrico Fermi OTS (Michigan). The second case study covers the implementation of a GE main generator excitation system on the Cernavoda OTS (Romania). The third case study covers the implementation of a Siemens main feedwater control system on the Callaway OTS and ES (Missouri). The final case study covers the implementation of the non-safety related and safety related DCS on the Nuclear New Build (NNB) Hongyanhe Phase I dual-purpose OTS/ES (China). (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extension of the ADAPT Framework for Multiple Simulators

Probabilistic Risk Assessment (PRA) uses fault tree/event tree analysis to evaluate the environmental impacts of nuclear power plants due to internal or external initiating events. In a fault tree, basic failure events are assembled using primarily And/Or logic to determine the combinations of failures that lead to failure of the system as a whole. An event tree is forward-facing, and begins with a single initiating event. From there the analysis considers what event may occur next, and branches out among the possible configurations of occurrence and non-occurrence of that event. This branching continues until end states are reached. End states may include a safe and stable configuration of the plant, or any number of differing failure states. Both fault and event trees require basic event probabilities as inputs to provide insight on the likelihood of different outcomes. PRA has been applied to nuclear power plants as an analysis tool since the 1975 Reactor Safety Study. In a traditional PRA, the order of events is prescribed by the analyst and each event is typically a binary aleatory uncertainty: occurrence or non-occurrence. Discrete dynamic event tree (DDET) analysis eliminates this subjectivity in ordering of events by using the output of a dynamic system model (simulator) to inform the branching. Branching conditions are triggered by the existence of a relevant plant state in the code, and therefore only occur as physically appropriate. The DDET approach also allows consideration of both epistemic and aleatory uncertainties on a phenomenologically and stochastically consistent platform. One limitation of the DDET driver codes developed to date (namely ADAPT, ADS, and MCDET ) is that the tree is driven by a single simulator. The goal of this work is a general multi-simulator driver. This limits the parameter space that may be explored in a single DDET to that which may be covered by a single simulator. In this work, the Analysis of Dynamic Accident Progression Trees (ADAPT) code is extended to support multiple arbitrary simulators. This capability will increase the depth and breadth of phenomena that may be analyzed in a single DDET. For example, a large-scale system code may pass the analysis to a more detailed code if certain conditions are met, in order to capture relevant phenomena that would be lost if only the system code were to be used. In addition, the time required to cover a given uncertainty space may be reduced by dynamically identifying and running only the situations that are appropriate given the plant conditions in each simulator. The use of multiple simulators is expected to further reduce subjectivity in PRA by enabling high fidelity modeling of reality. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sources and pathways of the upscale effects on the Southern Hemisphere jet in MPAS-CAM4 variable-resolution simulations

Impacts of regional grid refinement on large-scale circulations (“upscale effects”) were detected in a previous study that used the Model for Prediction Across Scales-Atmosphere coupled to the physics parameterizations of the Community Atmosphere Model version 4. The strongest upscale effect was identified in the Southern Hemisphere jet during austral winter. This study examines the detailed underlying processes by comparing two simulations at quasi-uniform resolutions of 30 and 120 km to three variable-resolution simulations in which the horizontal grids are regionally refined to 30 km in North America, South America, or Asia from 120 km elsewhere. In all the variable-resolution simulations, precipitation increases in convective areas inside the high-resolution domains, as in the reference quasi-uniform high-resolution simulation. With grid refinement encompassing the tropical Americas, the increased condensational heating expands the local divergent circulations (Hadley cell) meridionally such that their descending branch is shifted poleward, which also pushes the baroclinically unstable regions, momentum flux convergence, and the eddy-driven jet poleward. This teleconnection pathway is not found in the reference high-resolution simulation due to a strong resolution sensitivity of cloud radiative forcing that dominates the aforementioned teleconnection signals. The regional refinement over Asia enhances Rossby wave sources and strengthens the upper level southerly flow, both facilitating the cross-equatorial propagation of stationary waves. Evidence indicates that this teleconnection pathway is also found in the reference high-resolution simulation. Lastly, the result underlines the intricate diagnoses needed to understand the upscale effects in global variable-resolution simulations, with implications for science investigations using the computationally efficient modeling framework.

54 ENVIRONMENTAL SCIENCES↗

Portable interactive visualization of large-scale simulations in geotechnical engineering using Unity3D

Development in large-scale geotechnical engineering simulation places tremendous demand for efficient visualization of such simulation data. This study presents a lightweight software tool, i.e. Geotechnical Interactive Visualization (GIV), as a solution to this challenge, which achieves efficient interactive visualization of large-scale simulations data in geotechnical engineering. Visualization data flow and algorithms specifically optimized for common geotechnical engineering applications are implemented in GIV. GIV can visualize geotechnical structure models with time-varying attributes attached to mesh with fixed topology, and also models with time-varying mesh topologies but no attributes attached, the two most common visualization tasks in geotechnical engineering. Furthermore, challenges for large-scale simulation data visualization, including parallel simulation data redundancy, massive data size, dynamic user interaction, and portability are overcome via specifically designed algorithms for simulation data preprocessing and optimized visualization modules using the powerful 3D rendering and interactive game engine Unity3D. Comparison of GIV with several widely used visualization tools for the visualization of large-scale idealized datasets and realistic geotechnical simulations highlights the visualization efficiency, smooth interactivity, and lightweight features of GIV.

42 ENGINEERING↗

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗

Closed-loop simulation with Grad-Shafranov equilibrium evolution for plasma control system development

Plasma control design increasingly depends on fast simulations able to connect to operational plasma control system (PCS) software for iterative development. GSevolve is a nonlinear free-boundary simulation that evolves the Grad-Shafranov equilibrium including current and pressure profiles, and can connect to all versions of the DIII-D PCS operational on devices around the world. Its ability to simulate the evolution of current continuously from open to closed field lines enables GSevolve to simulate full-discharge control from breakdown through plasma termination. GSevolve is compliant with the Plasma Control System Simulation Platform (PCSSP), a Matlab/Simulink®-based environment that serves as the primary framework for development of the ITER Plasma Control System (PCS). The code employs a novel response function-updating scheme to achieve fast simulations. The plasma evolves linearly until a significant shape change triggers a nonlinear update of the response while enforcing a prescribed accuracy. Old responses are recalled if the plasma changes back to an old shape. The equilibrium is a function of a state vector containing axisymmetric conductor currents and parameters related to parallel and perpendicular currents, which evolves according to an efficient state space circuit formulation. Control design and validation results are presented for simulations of DIII-D and EAST, connected to their corresponding versions of the DIII-D PCS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D xRAGE simulation of inertial confinement fusion implosion with imposed mode 2 laser drive asymmetry

Low-mode asymmetries represent an important obstacle to achieving high-gain inertial confinement fusion implosions. As a step in learning how to control such effects, an OMEGA experiment with imposed mode 2 laser drive asymmetries was done to study the expected signatures of this type of asymmetry [M. Gatu Johnson et al., PRE 2018]. In the present work, a 3D xRAGE simulation including the stalk mount has been brought to bear on the data from that experiment. Furthermore, comprehensive comparisons between simulated and measured observables are made. Good agreement between simulated and measured x-ray image-inferred shell trajectories, bang times and neutron emission widths are seen, showing that the hydrodynamics are well captured in the simulation. Asymmetries seen in simulated and measured time-resolved and time-integrated x-ray images and areal densities also compare well, showing impact of both stalk and mode 2. On the other hand, important differences in measured and simulated neutron emission histories, yield, and ion temperature (T ion ) asymmetries are seen, suggesting that the simulation is overestimating shock yield. Overall, the results clearly demonstrate the importance of considering all asymmetry sources when interpreting measured signatures of asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of Geant4 Monte Carlo simulation in dose calculations for small radiosurgical fields

The Geant4 toolkit was used to develop a Monte Carlo (MC)-based engine for accurate dose calculations in small radiation field sizes. The Geant4 toolkit (version 10.1.p02) was used to simulate 6-MV photon beam of a Varian2100C linear accelerator that is being used for stereotactic radiosurgery (SRS) treatment with small radiation fields. Geometric models of 3 in-house designed radiosurgical divergent cones, with the diameters of their projections at the isocenter being 10, 20, and 30 mm, were simulated. The accuracy of the MC simulation technique was examined by reproducing several different simulated dosimetric parameters of the primary beams with the experimental data. The dose distributions are first checked for single beams for each cone, then standard multiple field (SMF) techniques are applied. A sample set of DICOM files from computed tomography (CT) scan imaging of a patient's head was converted to the Geant4 geometry format to implement MC-based engine for a clinical test. To validate the accuracy of the MC-based calculations for SMF arrangements, the isodose lines from MC simulation in water phantom were compared with the measured isodose lines using EBT3 Gafchromic film in Solid Water phantoms. Agreements between measured and simulated depth dose values and beam profiles for SRS cones were generally within 2%/2 mm. For output factors, the largest discrepancy was observed for 10 mm SRS cone, which was 1.7%. For SMF techniques, in SRS cones, the MC simulation and EBT3 Gafchromic film dosimetry were in acceptable agreement (5%/5 mm). Excellent agreement between the results of the MC-based and measured dose values for both single and SMF techniques in SRS cones indicates the ability of the Geant4 toolkit to be applied as the platform for treatment planning of advanced radiotherapy techniques.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

POLARRIS: A POLArimetric Radar Retrieval and Instrument Simulator

This paper introduces a synthetic polarimetric radar simulator and retrieval package, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), for evaluating cloud-resolving models (CRMs). POLARRIS is composed of forward (POLARRIS-f) and inverse (retrieval and diagnostic) components (iPOLARRIS) to generate not only polarimetric radar observables ($Z_h, Z_{dr}, K_{dp}, ρ_{hv})$ but also radar-consistent geophysical parameters such as hydrometeor identification, vertical velocity, and rainfall rates retrieved from CRM data. To demonstrate its application and uncertainties, POLARRIS is applied to simulations of a mesoscale convective system over the Southern Great Plains on 23 May 2011, using the Weather Research and Forecasting model with both spectral bin microphysics (SBM) and the Goddard single-moment bulk 4ICE microphysics. Statistical composites reveal a significant dependence of simulated polarimetric observables $(Z_{dr}, K_{dp})$ on the assumptions of the particle axis ratio (oblateness) and orientation angle distributions. We note the simulated polarimetric variables differ considerably between the SBM and 4ICE microphysics in part due to the differences in their ice particle size distributions as revealed by comparisons with aircraft measurements. Regardless of these uncertainties, simulated hydrometeor identification distributions overestimate graupel and hail fractions, especially from the simulation with SBM. To minimize uncertainties in forward model, the particle shape and orientation angle distributions of frozen particles should be predicted in a microphysics scheme in addition to the size distributions and particle densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Long-time simulations for fixed input states on quantum hardware

Publicly accessible quantum computers open the exciting possibility of experimental dynamical quantum simulations. While rapidly improving, current devices have short coherence times, restricting the viable circuit depth. Despite these limitations, we demonstrate long-time, high fidelity simulations on current hardware. Specifically, we simulate an XY-model spin chain on Rigetti and IBM quantum computers, maintaining a fidelity over 0.9 for 150 times longer than is possible using the iterated Trotter method. Our simulations use an algorithm we call fixed state Variational Fast Forwarding (fsVFF). Recent work has shown an approximate diagonalization of a short time evolution unitary allows a fixed-depth simulation. fsVFF substantially reduces the required resources by only diagonalizing the energy subspace spanned by the initial state, rather than over the total Hilbert space. We further demonstrate the viability of fsVFF through large numerical simulations, and provide an analysis of the noise resilience and scaling of simulation errors.

97 MATHEMATICS AND COMPUTING↗

Adaptive time steps in trajectory surface hopping simulations

Trajectory surface hopping (TSH) simulations are often performed in combination with active-space multi-reference configuration interaction (MRCI) treatments. Technical problems may arise in such simulations if active and inactive orbitals strongly mix and switch in some particular regions. We propose to use adaptive time steps when such regions are encountered in TSH simulations. For this purpose, we present a computational protocol that is easy to implement and increases the computational effort only in the critical regions. We test this procedure through TSH simulations of a GFP chromophore model (OHBI) and a light-driven rotary molecular motor (F-NAIBP) on semiempirical MRCI potential energy surfaces, by comparing the results from simulations with adaptive time steps to analogous ones with constant time steps. For both test molecules, the number of successful trajectories without technical failures rises significantly, from 53% to 95% for OHBI and from 25% to 96% for F-NAIBP. The computed excited-state lifetime remains essentially the same for OHBI and increases somewhat for F-NAIBP, and there is almost no change in the computed quantum efficiency for internal rotation in F-NAIBP. We recommend the general use of adaptive time steps in TSH simulations with active-space CI methods because this will help to avoid technical problems, increase the overall efficiency and robustness of the simulations, and allow for a more complete sampling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Speeding up N -body simulations of modified gravity: chameleon screening models

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied f ( R ) gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of f ( R ) simulations. For example, a test simulation with 512{sup 3} particles in a box of size 512 Mpc/ h is now 5 times faster than before, while a Millennium-resolution simulation for f ( R ) gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

79 ASTRONOMY AND ASTROPHYSICS↗

Divertor heat flux simulations in ELMy H-mode discharges of EAST

This paper presents heat flux simulations for the ELMy H-mode on the Experimental Advanced Superconducting Tokamak (EAST) using a six-field two-fluid model in BOUT++. Three EAST ELMy H-mode discharges with different plasma currents I p and geometries are studied. The trend of the scrape-off layer width $\lambda_{q}$ with I p is reproduced by the simulation. The simulated width is only half of that derived from the EAST scaling law, but agrees well with the international multi-machine scaling law. Note that there is no radio-frequency (RF) heating scheme in the simulations, and RF heating can change the boundary topology and increase the flux expansion. Anomalous electron transport is found to contribute to the divertor heat fluxes. A coherent mode is found in the edge region in simulations. The frequency and poloidal wave number $k_{\theta}$ are in the range of the edge coherent mode in EAST. The magnetic fluctuations of the mode are smaller than the electric field fluctuations. Statistical analysis of the type of turbulence shows that the turbulence transport type (blobby or turbulent) does not influence the heat flux width scaling. Here, the two-point model differs from the simulation results but the drift-based model shows good agreement with simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗