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,765 records · Page 98

Interdecadal Connection Between Artic Temperature and Summer Precipitation Over the Yangtze River Valley in the CMIP5 Historical Simulations

This study assesses the ability of the Phase 5 Coupled Model Intercomparison Project (CMIP5) simulations in capturing the interdecadal precipitation enhancement over the Yangtze River valley (YRV) and investigates the contributions of Arctic warming to the interdecadal variability of the East Asian summer monsoon rainfall. Six CMIP5 historical simulations including models from Canada (CCCma), China (BCC), Germany (MPI-M), Japan (MRI), United Kingdom (MOHC), and United States (NCAR) are used. The NCEP/NCAR reanalysis and observed precipitation are also used for comparison. Among the six CMIP5 simulations, only CCCma can approximately simulate the enhancement of interdecadal summer precipitation over the YRV in 1990-2005 relative to 1960-1975, and the relationships between the summer precipitation with surface temperature (Ts), the 850hPa winds, and 500hPa height field (H500), and between Ts and H500 using regression, correlation, and SVD analyses. It is found that CCCma can reasonably simulate the interdecadal surface warming over the boreal mid-to high latitudes and the Arctic in winter, spring and summer. The summer Baikal blocking appears to be the bridge that links the winter and spring surface warming over the mid-to high latitude and Arctic with the enhancement of summer precipitation over the YRV. Models that missed some or all of these relationships found in CCCma and the reanalysis failed to simulate the interdecadal enhancement of precipitation over the YRV. This points to the importance of high latitude and Arctic processes on interdecadal variability of the East Asian summer monsoon and the challenge for global climate models to correctly simulate the linkages.

arctic↗

Race conditions and data partitioning: risks posed by common errors to reproducible parallel simulations

When parallel algorithms for simulation were introduced in the 1970s, their development and use interested only experts in parallel computation. This circumstance changed as multi-core processors became commonplace, putting a parallel computer into the hands of every modeler. A natural outcome is growing interest in parallel simulation among persons not intimately familiar with parallel computing. At the same time, parallel simulation tools continue to be developed with the implicit assumption that the modeler is knowledgeable about parallel programming. The unintended consequence is a rapidly growing number of users of parallel simulation tools that are unlikely to recognize when the interaction of race conditions, partitioning strategies, and simultaneous action in their simulation models make results non-reproducible, thereby calling into question the validity of conclusions drawn from the simulation data. Here, we illustrate the potential dangers of exposing parallel algorithms to users who are not experts in parallel computation with example models constructed using existing parallel simulation tools. By doing so, we hope to refocus tool developers on usability, even if this new focus incurs loss of some performance.

97 MATHEMATICS AND COMPUTING↗

A scalable framework for the global offline community land model ensemble simulation

Current earth system models have a large range of uncertainty, owing to differences in the simulation of feedbacks and insufficient information to constrain model parameters. Parameter disturbance experiment provides a straightforward method to quantify the variation (uncertainty) outputs caused by model inputs. Owing to the software complexity and computational intensity of earth system models, a large-scale simulation framework is needed to support ensemble simulation required by parameter disturbance experiment. This paper presents a parallel framework for the community land model ensemble simulation. After a software structure review of the community land model simulation, a single factor parameter disturbance experiment of a reference computational experiment design is used to demonstrate the software design principles, computational characteristics of individual application, parallel ensemble simulation implementation, as well as the weak scalability of this simulation framework on a high-end computer. Lastly, the paper discusses some preliminary diagnostic analysis results of the single factor parameter disturbance experiments. The framework design considerations and implementation details described in this paper can be beneficial to many other research programmes involving large scale, legacy modelling system.

58 GEOSCIENCES↗

Evaluation of AP-Farm Simulant Composition for Rotary Microfilter Testing

This document identifies the feed composition of a Hanford AP tank farm simulant for rotary microfiltration testing. The composition is based on an Hanford Tank Waste Operations Simulator (HTWOS) model run in combination with Tank Waste Information Network (TWINS) data and mineralogical studies of actual waste solids. The feed simulant is intended to be used in test runs at SRNL. The simulant will be prepared in two parts: (1) A supernate, composed of water-soluble salts and (2) The undissolved (actually, undissolvable) solids. Test slurries with distinct solids concentrations (e.g., 0.5, 5 and 10 wt%) are then prepared as needed. The base for the composition of supernate and solids is the modeled feed sequence for a deployment scenario of the Supplemental Pretreatment units within AP-farm. These units comprise a filtration part, the RMF, and a Cesium-removal part, a Small Column Ion Exchange. The primary use of this simulant is for filtration testing - however, in case that it is also used for ion-exchange tests, the amount of Cs-137 that would need to be added is available in Table 1 and Attachment 3. A modified model run (MMR-049) of the Hanford Tank Waste Operations Simulator (HTWOS) system plan 6 case 3 was performed to identify the feed sequence. Case 3 assumed supplemental treatment besides the low activity waste (LAW) melter with supplemental pretreatment supporting the pretreatment facility. The MMR did not cap the duration of supplemental pretreatment to 15 months, but rather used it throughout the entire treatment mission as an add-on option to the pretreatment facility at the Waste Treatment and Immobilization Plant (WTP). Tank 241-AP-105 (AP-105) was chosen as the feed tank to the filtration unit. Other parameters included a fixed minimum of 0.5 wt% solids in the feed and a maximum Na-concentration of 5M in the supernate. The solids rejection from the filtration unit was set to 99.99% and the maximum allowed amount of solids within tank AP-105 was set to 10 wt%. A comprehensive description of the run and the full suite of results were issued as SVF-2364-00. The list of individual feed events including the amounts of liquid and solids transferred for the first five years is added as Attachment 2; the chemical composition of the supernate feed comprises Attachment 3. For the simulant composition, only the first five years of proposed feed delivery were taken into account. The main outcome of MMR-049 was that for the first five years, the feed would come mostly from AP-farms. Multiple delivery campaigns to AP-105 are included in this average feed, while minimizing the amount of contributing tanks to the solids in the feed mix.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparison of Waste Feed Delivery Small Scale Mixing Demonstration Simulant to Hanford Waste

The Hanford double-shell tank (DST) system provides the staging location for waste feed delivery to the Hanford Tank Waste Treatment and Immobilization Plant (WTP). Hall (2008) includes WTP acceptance criteria that describe physical and chemical characteristics of the waste that must be certified as acceptable before the waste is transferred from the DSTs to the WTP. One of the more challenging requirements relates to the sampling and characterization of the undissolved solids (UDS) in a waste feed DST. The objectives of Washington River Protection Solutions' (WRPS) Small Scale Mixing Demonstration (SSMD) project are to understand and demonstrate the DST sampling and batch transfer performance at multiple scales using slurry simulants comprised of UDS particles and liquid (Townson 2009). The SSMD project utilizes geometrically scaled DST feed tanks to generate mixing, sampling, and transfer test data. In Phase 2 of the testing, RPP-49740, the 5-part simulant defined in RPP-48358 was used as the waste slurry simulant. The Phase 2 test data are being used to estimate the expected performance of the prototypic systems in the full-scale DSTs. As such, understanding of the how the small-scale systems as well as the simulant relate to the full-scale DSTs and actual waste is required. The focus of this report is comparison of the size and density of the 5-part SSMD simulant to that of the Hanford waste. This is accomplished by computing metrics for particle mobilization, suspension, settling, transfer line intake, and pipeline transfer from the characterization of the 5-part SSMD simulant and characterizations of the Hanford waste. In addition, the effects of the suspending fluid characteristics on the test results are considered, and a computational fluid dynamics tool useful to quantify uncertainties from simulant selections is discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

User's Guide of TOUGH2-EGS. A Coupled Geomechanical and Reactive Geochemical Simulator for Fluid and Heat Flow in Enhanced Geothermal Systems Version 1.0

TOUGH2-EGS is a numerical simulation program coupling geomechanics and chemical reactions for fluid and heat flows in porous media and fractured reservoirs of enhanced geothermal systems. The simulator includes the fully-coupled geomechanical (THM) module, the fully-coupled geochemical (THC) module, and the sequentially coupled reactive geochemistry (THMC) module. The fully-coupled flow-geomechanics model is developed from the linear elastic theory for the thermo-poro-elastic system and is formulated with the mean normal stress as well as pore pressure and temperature. The chemical reaction is sequentially coupled after solution of flow equations, which provides the flow velocity and phase saturation for the solute transport calculation at each time step. In addition, reservoir rock properties, such as porosity and permeability, are subjected to change due to rock deformation and chemical reactions. The relationships between rock properties and geomechanical and chemical effects from poro-elasticity theories and empirical correlations are incorporated into the simulator. This report provides the user with detailed information on both mathematical models and instructions for using TOUGH2-EGS for THM, THC or THMC simulations. The mathematical models include the fluid and heat flow equations, geomechanical equation, reactive geochemistry equations, and discretization methods. Although TOUGH2-EGS has the capability for simulating fluid and heat flows coupled with both geomechanical and chemical effects, it is up to the users to select the specific coupling process, such as THM, THC, or THMC in a simulation. There are several example problems illustrating the applications of this program. These example problems are described in details and their input data are presented. The results demonstrate that this program can be used for field-scale geothermal reservoir simulation with fluid and heat flow, geomechanical effect, and chemical reaction in porous and fractured media.

15 GEOTHERMAL ENERGY↗

User's guide of TOUGH2-EGS-MP: A Massively Parallel Simulator with Coupled Geomechanics for Fluid and Heat Flow in Enhanced Geothermal Systems VERSION 1.0

TOUGH2-EGS-MP is a parallel numerical simulation program coupling geomechanics with fluid and heat flow in fractured and porous media, and is applicable for simulation of enhanced geothermal systems (EGS). TOUGH2-EGS-MP is based on the TOUGH2-MP code, the massively parallel version of TOUGH2. In TOUGH2-EGS-MP, the fully-coupled flow-geomechanics model is developed from linear elastic theory for thermo-poro-elastic systems and is formulated in terms of mean normal stress as well as pore pressure and temperature. Reservoir rock properties such as porosity and permeability depend on rock deformation, and the relationships between these two, obtained from poro-elasticity theories and empirical correlations, are incorporated into the simulation. This report provides the user with detailed information on the TOUGH2-EGS-MP mathematical model and instructions for using it for Thermal-Hydrological-Mechanical (THM) simulations. The mathematical model includes the fluid and heat flow equations, geomechanical equation, and discretization of those equations. In addition, the parallel aspects of the code, such as domain partitioning and communication between processors, are also included. Although TOUGH2-EGS-MP has the capability for simulating fluid and heat flows coupled with geomechanical effects, it is up to the user to select the specific coupling process, such as THM or only TH, in a simulation. There are several example problems illustrating applications of this program. These example problems are described in detail and their input data are presented. Their results demonstrate that this program can be used for field-scale geothermal reservoir simulation in porous and fractured media with fluid and heat flow coupled with geomechanical effects.

15 GEOTHERMAL ENERGY↗

Extended Leach Testing of Simulated LAW Cast Stone Monoliths

This report describes the results from long-term laboratory leach tests performed at Pacific Northwest National Laboratory (PNNL) for Washington River Protection Solutions (WRPS) to evaluate the release of key constituents from monoliths of Cast Stone prepared with four simulated low-activity waste (LAW) liquid waste streams. Specific objectives of the Cast Stone long-term leach tests described in this report focused on four activities: 1. Extending the leaching times for selected ongoing EPA-1315 tests on monoliths made with LAW simulants beyond the conventional 63-day time period up to 609 days reported herein (with some tests continuing that will be documented later) in an effort to evaluate long-term leaching properties of Cast Stone to support future performance assessment activities. 2. Starting new EPA-1315 leach tests on archived Cast Stone monoliths made with four LAW simulants using two leachants (deionized water [DIW] and simulated Hanford Integrated Disposal Facility (IDF) Site vadose zone pore water [VZP]). 3. Evaluating the impacts of varying the iodide loading (starting iodide concentrations) in one LAW simulant (7.8 M Na Hanford Tank Waste Operations Simulator (HTWOS) Average) by manufacturing new Cast Stone monoliths and repeating the EPA-1315 leach tests using DIW and the VZP leachants. 4. Evaluating the impacts of using a non-pertechnetate form of Tc that is present in some Hanford tanks. In this activity one LAW simulant (7.8 M Na HTWOS Average) was spiked with a Tc(I)-tricarbonyl gluconate species and then solidified into Cast Stone monoliths. Cured monoliths were leached using the EPA-1315 leach protocol with DIW and VZP. The leach results for the Tc-Gluconate Cast Stone monoliths were compared to Cast Stone monoliths pertechnetate.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluation of CFD Methods for Simulation of Two-Phase Boiling Flow Phenomena in a Helical Coil Steam Generator

The U.S. Department of Energy, Office of Nuclear Energy charges participants in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program with the development of advanced modeling and simulation capabilities that can be used to address design, performance and safety challenges in the development and deployment of advanced reactor technology. The NEAMS has established a high impact problem (HIP) team to demonstrate the applicability of these tools to identification and mitigation of sources of steam generator flow induced vibration (SGFIV). The SGFIV HIP team is working to evaluate vibration sources in an advanced helical coil steam generator using computational fluid dynamics (CFD) simulations of the turbulent primary coolant flow over the outside of the tubes and CFD simulations of the turbulent multiphase boiling secondary coolant flow inside the tubes integrated with high resolution finite element method assessments of the tubes and their associated structural supports. This report summarizes the demonstration of a methodology for the multiphase boiling flow analysis inside the helical coil steam generator tube. A helical coil steam generator configuration has been defined based on the experiments completed by Polytecnico di Milano in the SIET helical coil steam generator tube facility. Simulations of the defined problem have been completed using the Eulerian-Eulerian multi-fluid modeling capabilities of the commercial CFD code STAR-CCM+. Simulations suggest that the two phases will quickly stratify in the slightly inclined pipe of the helical coil steam generator. These results have been successfully benchmarked against both empirical correlations for pressure drop and simulations using an alternate CFD methodology, the dispersed phase mixture modeling capabilities of the open source CFD code Nek5000.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Relativistic Modeling Capabilities in Perseus Extended MHD Simulation Code for HED Plasmas

This project is submitted in response to DE-FOA-0001153, and will be carried out by Dr. Charles Seyler and Dr. Nathaniel Hamlin as post doc under the direction of Professor Seyler. The objectives are to enable PERSEUS to self-consistently model relativistic high-energy-density (HED) laboratory and astrophysical phenomena. Electrons and ions in a plasma are generally modeled either as fluids using a two-fluid code, like PERSEUS, or as particles using a Particle-In- Cell (PIC) code. Nonrelativistic PERSEUS solves the two-fluid equations, formulated in terms of a generalized Ohm’s law (GOL), so as to model about nine orders of magnitude in density variation using a local semi-implicit method. Relativistic PERSEUS preserves this structure and therefore retains these advantageous properties, enabling it to model a broader range of relativistic HED phenomena than previous codes. We have also overcome a major technical challenge associated with solving a relativistic system of equations relating conserved quantities (momentum, energy, etc.) to primitive variables (density, velocity, and pressure). We have thus far developed a relativistic version of PERSEUS and demonstrated that it recovers expected nonrelativistic results. We are currently using relativistic PERSEUS to model two HED phenomena, namely laser-plasma interactions and X-pinches, both of which show spectroscopic evidence of relativistic electrons. In an X-pinch, a section of wire carrying sufficiently large current forms a plasma column that pinches down to a point-like source, generating X-rays for use in radiographic imaging. Laser-plasma interactions have been simulated using PIC codes, and we use their results as a benchmark for comparison. However, neither PIC codes nor existing two-fluid schemes can simulate X-pinches, due to the broad range of densities which have thus far made such a simulation computationally prohibitive. In both laser-plasma interaction and Xpinch simulations, we have observed relativistic phenomena that PERSEUS could not previously model. This includes relativistic channeling of a laser into an over-dense deuterium gas (i.e. induced transparency), consistent with PIC simulation results, along with energetic electrons and ions following an X-pinch and likely accompanied by strong X-ray emission, consistent with experimental observations of energetic X-rays. This proposal includes our latest simulation results. We will improve computational efficiency with techniques that allow the use of a less refined numerical grid. This will enable extensions to three-dimensional geometries. We will improve modeling of laser-plasma interactions and X-pinches, including radiative transport processes, in order to improve comparisons with PIC simulations and experimental diagnostics. We will also model power loss through electrode surface plasmas in magnetically insulated transmission lines (MITLs). We will model relativistic astrophysical phenomena in both astrophysical and HED regimes, which could shed light on the validity of laboratory astrophysics for certain problems. These phenomena are currently modeled using PIC codes and existing relativistic two-fluid codes, neither of which can model the range of length scales of which relativistic PERSEUS is capable. There is strong potential for this project to have a broad impact on both the astrophysical and laboratory plasma communities. In the laboratory community, our work is applicable to those concerned with laser-plasma interactions such as laser-driven fusion, and those in the dense Z pinch community concerned with X-pinches and electrode surface plasmas. An improved understanding of the mechanisms of hard (i.e. high-energy) X-ray generation will likely make important contributions toward the use of X-pinches for radiography through the harnessing or mitigation of hard X-rays, depending on the application. The modeling of electrode surface plasmas, including vacuum electron flow, will improve our understanding of power loss to various loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Roadmap for IES Modeling and Simulation Activities: Update 2023

Since inception of the DOE-NE Integrated Energy Systems (IES) program, several program and modeling visions and roadmaps have been published. In 2020, a last comprehensive roadmap for the DOE-NE Integrated Energy Systems program was published. In 2017, a separate modeling and simulation capability development plan was published that was extensively updated in 2022. The current report incrementally updates the 2022 modeling and simulation capability roadmap. This is not a new report, but individual sections have been updated to reflect recent accomplishments and new directions of the program. The role of modeling and simulation within IES is to support the demonstration of new coupled integrated technologies along every step of the technology maturation from the strategic analysis of preferred system architecture to preliminary design, to laboratory testing, up to full commercial testing and integration. To achieve this role, modeling and simulation must be able to assess the technical performance and the economic viability of potential IES. Also, modeling and simulation must provide support for experimental evaluations (i.e., support component design and real-time operations of experimental demonstration systems). Finally, modeling and simulation can help scale-up experimental systems to the final commercial systems. This report details the current state of the modeling and simulation efforts within IES for the above-mentioned areas, provides a gap analysis, and proposes next steps for a time horizon over the next years.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Collaborative: in situ visual analytics technologies for extreme scale combustion simulations

This project aims to drastically enhance the usability of in situ analysis and visualization for extreme-scale scientific simulations. Current exascale computing capabilities promise to offer greater predictive ability of simulations and to further push the frontiers of science and technology. However, to validate the simulation output at extreme scale, examine the modeled phenomena, and discover previously unknowns from the output data, the output must be reduced or transformed in situ as it is being generated during the simulation such that the amount of data to examine and store is kept to a minimum. Such in situ approaches allow us to process and analyze the data and any embedded geometry to an extent that would be prohibitively expensive, if not impossible, to perform as a post hoc task. While in situ processing has been demonstrated to be a feasible and promising approach, its full potential has not yet been leveraged. In this project, we have developed comprehensive enhancements to in situ technology based on probability distributions in data. Our research focuses on jointly developing new ways of interacting with massive statistical samples while creatively utilizing new state-of-the-art computational resources to push the boundaries of in situ exploration. Moreover, we have developed new time-dependent techniques to enable previously unattainable capabilities in areas such as intelligent simulation steering and precise feature identification. We have experimentally studied our design and implementation at NERSC and OLCF, and are able to leverage existing in situ infrastructures whenever possible. While the exemplar in this project is combustion, many other fields for which turbulent transport is important, e.g., fusion, climate, astrophysics among others, encounter similar issues as simulations scale up to the exascale. This project shows its potential to generate high impact on DOE missions since the resulting technology promises to improve scientists’ ability to rapidly and correctly interpret and tune extreme-scale simulations, leading to new scientific understanding and advancements.

97 MATHEMATICS AND COMPUTING↗

T&D LIBRA: Simulation Framework for Transmission and Distribution Load and IBR Assessment

This paper introduces T&D LIBRA, a HELICS-based dynamic co-simulation framework designed for transmission and distribution (T&D) load and inverter-based resource assessment. This framework is designed to model interactions between single-phase induction motors (i.e., Motor-D), distributed energy resources (DERs) and the transmission network, emphasizing system stability during fault conditions. The primary objective of this research is to study the fault-induced delayed voltage recovery (FIDVR) behavior of the T&D networks using the co-simulation framework and to compare the results with conventional PSS/E only simulations. Simulations are validated using T&D co-simulations with a single distribution node and the IEEE 123-bus distribution test case. The co-simulation results demonstrate the value of using T&D co-simulations to capture motor stalling, DER tripping and other sources of voltage diversity on distribution networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An Evaluation of Co-Simulation for Modeling Coupled Natural Gas and Electricity Networks

Reliance on natural gas for power generation has increased the coupling between gas and power networks. While this coupling can bring operational and economic benefits, it can also yield challenges, as the constraints in one system can impact the other. Co-simulation can capture the constraints and interactions between these systems, but so far, there has been limited comparison of co-simulation results to those of an integrated model. In this work, we develop a new co-simulation framework using the HELICS platform and the SAInt tool for modeling transient gas and AC optimal power flow. We evaluate this co-simulation framework against a fully integrated version of the SAInt power and gas simulators, thus providing a benchmarking of the co-simulation approach. We compare results across the two approaches for two test networks and a network representing the Belgian power and gas networks, testing both normal operating conditions and cases with compressor disruptions. In each of the cases tested, we find nearly identical results from the two approaches across various metrics of interest, such as nodal pressure, gas flow rates, and active power generation. This alignment suggests that co-simulation can yield comparable results to fully integrated models for modeling coupled gas and electricity networks.

03 NATURAL GAS↗

GLOBAL HIGH-RESOLUTION N-BODY SIMULATION OF PLANET FORMATION. I. PLANETESIMAL-DRIVEN MIGRATION

We investigated whether outward planetesimal-driven migration (PDM) takes place or not in simulations when the self-gravity of planetesimals is included. We performed N-body simulations of planetesimal disks with a large width (0.7–4 au) that ranges over the ice line. The simulations consisted of two stages. The first-stage simulations were carried out to see the runaway growth phase using the planetesimals of initially the same mass. The runaway growth took place both at the inner edge of the disk and at the region just outside the ice line. This result was utilized for the initial setup of the second-stage simulations, in which the runaway bodies just outside the ice line were replaced by the protoplanets with about the isolation mass. In the second-stage simulations, the outward migration of the protoplanet was followed by the stopping of the migration due to the increase of the random velocity of the planetesimals. Owing to this increase of random velocities, one of the PDM criteria derived in Minton and Levison was broken. In the current simulations, the effect of the gas disk is not considered. It is likely that the gas disk plays an important role in PDM, and we plan to study its effect in future papers.

79 ASTRONOMY AND ASTROPHYSICS↗

Determining pitch-angle diffusion coefficients from test particle simulations

The transport and acceleration of charged particles in turbulent media are topics of great interest in space physics and interstellar astrophysics. These processes are dominated by the scattering of particles off magnetic irregularities. The scattering process itself is usually described by small-angle scattering, with the pitch-angle coefficient D{sub μμ} playing a major role. Since the diffusion coefficient D{sub μμ} can be determined analytically only for the approximation of quasilinear theory, the determination of this coefficient from numerical simulations has become more important. So far these simulations have yielded particle tracks for small-scale scattering, which can then be interpreted using the running diffusion coefficients. This method has a limited range of validity. This paper presents two new methods that allow for the calculation of the pitch-angle diffusion coefficient from numerical simulations. These methods no longer analyze particle trajectories and instead examine the change of particle distribution functions. It is shown that these methods provide better resolved results and allow for the analysis of strong turbulence. The application of these methods to Monte Carlo simulations of particle scattering and hybrid MHD-particle simulations is presented. Both analysis methods are able to recover the diffusion coefficients used as input for the Monte Carlo simulations and provide better results in MHD simulations, especially for stronger turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

The Effect of Corner Modes in the Initial Conditions of Cosmological Simulations

In view of future high-precision large-scale structure surveys, it is important to quantify the percent and subpercent level effects in cosmological N-body simulations from which theoretical predictions are drawn. One such effect involves deciding whether to zero all modes above the one-dimensional Nyquist frequency, the so-called “corner” modes, in the initial conditions. We investigate this effect by comparing power spectra, density distribution functions, halo mass functions, and halo profiles in simulations with and without these modes. For a simulation with a mass resolution of m{sub p}∼10{sup 11} h{sup −1} M{sub ⊙}, we find that at z>6, the difference in the matter power spectrum is large at wavenumbers above ∼80% of k{sub Ny}, reducing to below 2% at all scales by z∼3. Including corner modes results in a better match between low- and high-resolution simulations at wavenumbers around the Nyquist frequency of the low-resolution simulation, but the effect of the corner modes is smaller than the effect of particle discreteness. The differences in mass functions are 3% for the smallest halos at z = 6 for the m{sub p}∼10{sup 11} h{sup −1} M{sub ⊙} simulation, but we find no significant difference in the stacked profiles of well-resolved halos at z⩽6. Thus removing power at |k|>k{sub Ny} in the initial conditions of cosmological simulations has a small effect on small scales and high redshifts, typically below a few percent.

79 ASTRONOMY AND ASTROPHYSICS↗