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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Simulations of the fault current limiting operation of a long-length REBCO CORC ® superconducting cable cooled by helium gas

Conductor-on-round-core (CORC®) cables composed of rare-earth-barium-copper-oxide high-temperature superconducting (HTS) tapes are of great interest for power transmission applications due to their many advantages such as high power density, light weight, and low loss. Closed circulation loops of cryogenic helium gas can be used to cool HTS cables down to low temperatures to significantly improve their current-carrying capacity. Coupled circuit-electromagnetic-thermal finite element simulations implemented in the COMSOL Multiphysics package were developed, validated, and then used for simulating the fault current limiting (FCL) performance and the cooling processes of an 8-layer CORC® cable cooled with a flow of cryogenic helium gas. In the simulations, the temperature dependence of the electrical and thermal properties of all component materials is implemented for improved accuracy. To overcome computational challenges caused by the considerable difference in geometrical scales (i.e. few-µm-thick HTS layers versus 10 m-long HTS cable), the model is divided into two separate simulations. The first simulation is performed on the transverse cross-section of the cable to calculate the electric field, heating power and temperature rise in each component of a CORC® cable during FCL operation. The heating power calculated in the first simulation is transferred to the second model to simulate the cooling of a 10 m-long cable after the fault is cleared. The effect of the helium gas flow rate on the cooling process is also investigated to develop strategic approaches for optimizing cooling systems for HTS cables with FCL capability. The simulations indicated that a 40 ms fault with a voltage drop of 20 V m −1 along the cable can result in a temperature increase from 60 K to about 165 K inside the cable, and it takes about 500 s to cool the cable back to nearly 60 K with a flow of cold helium gas at a rate of 5 g s −1 .

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2019

Verification and validation of integrated simulation of energetic particles in fusion plasmas

Here, we report on verification and validation of linear simulations of Alf´ven eigenmodes in the current ramp phase of DIII-D L-mode discharge #159243 using gyrokinetic, gyrokinetic-MHD hybrid, and eigenvalue codes. Using a classical fast ion profile, all simulation codes find that reversed shear Alf´ven eigenmodes (RSAE) are the dominant instability. The real frequencies from all codes have a coefficient of variation of less than 5% for the most unstable modes with toroidal mode number n = 4 and 5. Due to uncertainty in the experimentally measured minimum safety factor (q min ) and the sensitivity of the RSAE’s frequency to q min , the average difference in frequency between simulated and measured is df/f = 14.7%, 22.6%, for n = 4 and 5, respectively. If qmin is adjusted slightly, within experimental errors, an agreement of within 4.5% and 3.1%, for n = 4 and 5 respectively, is found for all codes. The simulated growth rates exhibit greater variation, and simulations find that pressure gradients of thermal plasmas make a significant contribution to the growth rates. Mode structures of the dominant modes agree well among all codes. Using a calculated fast ion profile that takes into account the action of the spectrum of multiple unstable modes, a toroidal Alf´ven eigenmode (TAE) with n = 6 is found to be unstable in the outer edge, consistent with the experimental data. Variations of the real frequencies and growth rates of the TAE are slightly larger than those of the RSAE. Lastly, electron temperature fluctuations and radial phase shifts from simulations show no significant differences with the experimental data for the strong n = 4 RSAE, but significant differences for the weak n = 6 TAE. The verification and validation for the linear simulations of Alf´ven eigenmodes is the first step to determine an integrated simulation of energetic particles confinement in burning plasmas incorporating multiple physical processes with kinetic effects of both energetic and thermal particles on an equal footing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Selecting durable building envelope systems with machine learning assisted hygrothermal simulations database

Hygrothermal simulations provide insight into the energy performance and moisture durability of building envelope components under dynamic conditions. The inputs required for hygrothermal simulations are extensive, and carrying out simulations and analyses requires expert knowledge. An expert system, the Building Science Advisor (BSA), has been developed to predict the performance and select the energy-efficient and durable building envelope systems for different climates. The BSA consists of decision rules based on expert opinions and thousands of parametric simulation results for selected wall systems. The number of potential wall systems results in millions, too many to simulate all of them. We present how machine learning can help predict durability data, such as mold growth, while minimizing the number of simulations needed to run. The simulation results are used for training and validation of machine learning tools for predicting wall durability. We tested Artificial Neural Network (ANN) and Gradient Boosted Decision Trees (GBDT) for their applicability and model accuracy. Models developed with both methods showed adequate prediction performance (root mean square error of 0.195 and 0.209, respectively). Finally, we introduce how the information supports guidance for envelope design via an easy-to-use web-based tool that does not require the end-user to run hygrothermal simulations.

Salonvaara, Mikael↗

DOE OSTI2024

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗

DOE OSTI2022

Hydrodynamic simulations of electron-capture supernovae: progenitor and dimension dependence

ABSTRACT We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in flash with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with six progenitor models that, at the time of collapse, span a range of propagating flame front radii. For comparison, we also perform a simulation with the same set-up using the canonical, spherically symmetrical progenitor model n8.8. We found that the variations in the progenitor models inherited from simulations of stellar evolution and flame propagation do not significantly alter the global properties of the neutrino-driven ECSN explosion, such as the explosion energy (∼1.36–1.48 × 1050 erg) and the mass (∼0.017–0.018 M⊙) and composition of the ejecta. Due to aspherical perturbations induced by the 2D flame, the ejecta contains a small amount (≲1.8 × 10−3 M⊙) of low-Ye (0.35 < Ye < 0.4) component. The baryonic mass of the protoneutron star is ∼1.34 M⊙ (∼1.357 M⊙) with the new (n8.8) progenitor models when simulations end at ∼400 ms and the discrepancy is due to updated weak-interaction rates in the progenitor evolutionary simulations. Our results reflect the nature of ECSN progenitors containing a strongly degenerate oxygen–neon–magnesium (ONeMg) core and suggest a standardized ECSN explosion initialized by ONeMg core collapse. Moreover, we carry out a rudimentary three-dimensional simulation and find that the explosion properties are fairly compatible with the 2D counterpart. Our paper facilitates a more thorough understanding of ECSN explosions following the ONeMg core collapse, though more three-dimensional simulations are still needed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Cosmological Evidence Modelling: a new simulation-based approach to constrain cosmology on non-linear scales

Extracting accurate cosmological information from galaxy–galaxy and galaxy–matter correlation functions on non-linear scales (${\lesssim } 10 \, h^{-1}{\rm {Mpc}}$) requires cosmological simulations. Additionally, one has to marginalize over several nuisance parameters of the galaxy–halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal {O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realizations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. In this work, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalized over the galaxy–halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalize over highly complex models of the galaxy–halo connection and greater accuracy, thereby reducing the number of simulations required.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Direct numerical simulation and Reynolds-averaged Navier-Stokes modeling of the sudden viscous dissipation for multicomponent turbulence

Simulations of a turbulent multicomponent fluid mixture undergoing isotropic deformations are carried out to investigate the sudden viscous dissipation. This dissipative mechanism was originally demonstrated using simulations of an incompressible single-component fluid [S. Davidovits and N. J. Fisch, Phys. Rev. Lett. 116, 105004 (2016)]. By accounting for the convective and diffusive transfer of various species, the current work aims to increase the physical fidelity of previous simulations and their relevance to inertial confinement fusion applications. Direct numerical simulations of the compressed fluid show that the sudden viscous dissipation of turbulent kinetic energy is unchanged from the single-component scenario. More importantly, the simulations demonstrate that the mass fraction variance and covariance for the various species also exhibit a sudden viscous decay. Reynolds-averaged Navier-Stokes simulations were carried out using the $k - l$ model to assess its ability to reproduce the sudden viscous dissipation. Results show that the standard $k - l$ formulation does not capture the sudden decay of turbulent kinetic energy, mass-fraction variance, and mass-fraction covariance for simulations with various compression and expansion rates, or different exponents for the power-law model of viscosity. A new formulation of the $k - l$ model that is based on previous improvements to the $k - ε$ family of models is proposed, which leads to consistently good agreement with the direct numerical simulations for all the isotropic deformations under consideration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

A Versatile Simulated Data Transport Layer for in Situ Workflows Performance Evaluation

In situ processing does not only allow scientific applications to face the explosion in data volume and velocity but also to address the time constraints of many simulation-analysis workflows by providing scientists with early insights about their applications at runtime. Multiple frameworks implement the concept of a data transport layer (DTL) to enable such in situ workflows. These tools are very versatile, directly or indirectly access the data generated on the same node, another node of the same compute cluster, or a completely distinct node, and allow data publishers and subscribers to run on the same computing resources or not. This versatility puts on researchers the onus of taking key decisions related to resource allocation and how to transport data to ensure the most efficient execution of their in situ workflows. However, domain scientists and workflow practitioners lack the appropriate tools to assess the respective performance of particular design and deployment options. In this paper we introduce a versatile simulated DTL designed to provide researchers with insights on the respective performance of different execution scenarios of in situ workflows. This open-source, standalone library builds on the SimGrid toolkit and can be linked to any SimGrid-based simulator. It facilitates the evaluation of the performance behavior, at scale, of different data transport configurations and the study of the effects of resource allocation strategies. We demonstrate the scalability, versatility, and accuracy of this simulated DTL by reproducing the execution of two synthetic benchmarks and of a real-world in situ workflow composed of an MPI application and a parallel data analysis. Results of simulations run on a single core show that the proposed library can simulate the interactions of tens of thousands of simulated processes deployed on two interconnected commodity clusters in a few seconds, and the execution by a thousand simulated processes of an in situ workflow in less than three minutes.

Suter, Fred [ORNL] (ORCID:0000000319021955)↗

DOE OSTI2019

Evaluating Carbon Extremes in a Coupled Climate-Carbon Cycle Simulation

Gross primary production (GPP) measures the photosynthetic update of carbon by terrestrial ecosystems. Accurately quantifying and simulating GPP and its extremes remains a challenge in global carbon cycle sciences. Here, we evaluate GPP extremes in a coupled biogeochemistry (BGC) simulation by the Department of Energy's Energy Exascale Earth System Model (E3SMv1.1) using the Generalized Extreme Value (GEV) distribution statistical model. The simulation is evaluated against the Global Bio-Atmosphere Flux (GBAF) data. Temporal trends and ENSO dependence are also investigated by using GEV models where time and the Niño3.4 index are introduced as linear covariates. The E3SMv1.1 model simulation generally predicts stronger negative and positive GPP extremes as compared to GBAF data. It also tends to simulate stronger temporal trends of GPP extremes than GBAF data. While negative GPP extreme trends are not significant in either E3SM or GBAF, positive GPP trends are statistically significant over several regions only for the E3SMv1.1 model simulation. ENSO dependence is generally stronger in the E3SMv1.1 model simulation, but ENSO dependence is found not to be significant for the time period analyzed (1980-2006) to match GBAF data. For the longer simulation period of 1900-2006, ENSO dependence is found to be statistically significant over Amazon, the maritime continent and Northern Australia for both negative and positive extremes.

Xu, Min↗

DOE OSTI2023

High-Performance Computing Based EMT Simulation of Large PV or Hybrid PV Plants

Faults in the transmission grid have led to reduced power generation from power electronics resources that are typically not connected to the faulted transmission line. In many of the cases, partial loss of power is observed within the power electronics resources like large photovoltaic (PV) power plants. This phenomena is not captured in existing simulation models and/or simulators. High-fidelity switched system electromagnetic transient (EMT) dynamic models of PV power plants can improve the fidelity of models available for accurate analysis of the impact on PV plants during simulation of faults. However, these models are extremely computationally expensive and take a long time to simulate. Long simulation times limit the ability to use these models as larger regions are studied in EMT simulations with more power electronics resources. In this paper, numerical simulation algorithms are combined with high-performance computing techniques and applied to the high-fidelity switched system EMT model of PV plants. Using these techniques, a speed-up of up to 58x is obtained, while preserving the accuracy of the simulation at greater than 98%.

Debnath, Suman↗

DOE OSTI2022

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Monte Carlo simulations of secondary electron emission due to ion beam milling

We present a Monte Carlo simulation study of secondary electron emission resulting from focused ion beam milling of a copper target. The basis of this study is a simulation code which simulates ion induced excitation and emission of secondary electrons, in addition to simulating focused ion beam sputtering and milling. This combination of features permits the simulation of the interaction between secondary electron emission, and the evolving target geometry as the ion beam sputters material. Previous ion induced SE Monte Carlo simulation methods have been restricted to predefined target geometries, while the dynamic target in the presented simulations makes this study relevant to image formation in ion microscopy, and chemically assisted ion beam etching, where the relationship between sputtering, and its effects on secondary electron emission, is important. We focus on a copper target, and validate our simulation against experimental data for a range of: noble gas ions, ion energies, ion/substrate angles and the energy distribution of the secondary electrons. We then provide a detailed account of the emission of secondary electrons resulting from ion beam milling; we quantify both the evolution of the yield as high aspect ratio valleys are milled, as well as the emission of electrons within these valleys that do not escape the target, but which are important to the secondary electron contribution to chemically assisted ion induced etching.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Tensor Network Quantum Virtual Machine for Simulating Quantum Circuits at Exascale

The numerical simulation of quantum circuits is an indispensable tool for development, verification, and validation of hybrid quantum-classical algorithms intended for near-term quantum co-processors. The emergence of exascale high-performance computing (HPC) platforms presents new opportunities for pushing the boundaries of quantum circuit simulation. Here, we present a modernized version of the Tensor Network Quantum Virtual Machine (TNQVM) that serves as the quantum circuit simulation backend in the eXtreme-scale ACCelerator (XACC) framework. The new version is based on the scalable tensor network processing library ExaTN (Exascale Tensor Networks). It provides multiple configurable quantum circuit simulators that perform either an exact quantum circuit simulation via the full tensor network contraction or an approximate simulation via a suitably chosen tensor factorization scheme. Upon necessity, stochastic noise modeling from real quantum processors is incorporated into the simulations by modeling quantum channels with Kraus tensors. By combining the portable XACC quantum programming frontend and the scalable ExaTN numerical processing backend, we introduce an end-to-end virtual quantum development environment that can scale from laptops to future exascale platforms. We report initial benchmarks of our framework, which include a demonstration of the distributed execution, incorporation of quantum decoherence models, and simulation of the random quantum circuits used for the certification of quantum supremacy on Google’s Sycamore superconducting architecture.

Nguyen, Thien↗

DOE OSTI2019

The Application of an Evolutionary Programming Process to a Simulation of the ETEX Large-Scale Airborne Dispersion Experiment

Abstract Airborne tracer simulations are typically performed using a dispersion model driven by a high-resolution meteorological model. Besides solving the dynamic equations of momentum, heat, and moisture on the resolved model grid, mesoscale models must account for subgrid-scale fluxes and other unresolved processes. These are estimated through parameterization schemes of eddy diffusion, convection, and surface interactions, and they make use of prescribed parameters set by the user. Such “free” model parameters are often poorly constrained, and a range of plausible values exists for each. Evolutionary programming (EP) is a process to improve the selection of the parameters. A population of simulations is first run with a different set of parameter values for each member, and the member judged most accurate is selected as the “parent” of a new “generation.” After a number of iterations, the simulations should approach a configuration that is best adapted to the atmospheric conditions. We apply the EP process to simulate the first release of the 1994 European Tracer Experiment (ETEX) project, which comprised two experiments in which a tracer was released in western France and sampled by an observing network. The EP process is used to improve a simulation of the RAMS mesoscale weather model, with weather data collected during ETEX being used to “score” the individual members according to how well each simulation matches the observations. The meteorological simulations from before and after application of the EP process are each used to force a dispersion model to create a simulation of the ETEX release, and substantial improvement is observed when these are validated against sampled tracer concentrations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

High fidelity blade-resolved and actuator line data from a 16 turbine wind farm simulation using ExaWind

This data was generated with the ExaWind code suite (https://github.com/Exawind) as a demonstration of a large, 16 turbine wind farm simulation, calculated using two different levels of fidelity. The lower level of fidelity approach uses an actuator line approach to represent the turbines, and was simulated with AMR-Wind (https://github.com/Exawind/amr-wind/) as the background flow solver, coupled to OpenFAST (https://github.com/OpenFAST/openfast). The higher level of fidelity simulation uses a blade-resolved approach, and is done using AMR-Wind, Nalu-Wind (https://github.com/Exawind/nalu-wind), OpenFAST, and TIOGA (https://github.com/Exawind/tioga). In the blade-resolved simulation, ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. OpenFAST handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. In the actuator line simulation, a mesh of 295M elements was used for a 5km x 5km domain, and it was simulated using 256 nodes (2048 GPU's) on the Oak Ridge Leadership Computing Facility Frontier supercomputer. For the blade-resolved simulation, 1.5B element mesh was used in the AMR-Wind background 5km x 5km domain, and 16M elements were used for each turbine in the Nalu-Wind domains, for a total of 1.7B elements. This was simulated using 384 nodes on Frontier, with each node using 56 cores for Nalu-Wind and 8 GPU cores. The data in this archive includes the turbine outputs from OpenFAST, 2D sampling planes from AMR-Wind, and full-field solution files from AMR-Wind and Nalu-Wind.

17 WIND ENERGY↗

DOE OSTI2023

Earth System Model Simulations of Ameriflux sites using dynamic roots and dynamic leaf to root allocation

This data package contains point-scale Earth System Model simulations run with DOE's E3SM Model. The package includes four sets of simulations, one set of simulations run in the default configuration, a second set of simulations run with the dynamic root module enabled, the third set with dynamic leaf to fine root allocation, and the fourth with the combination of dynamic roots and dynamic allocation. These simulations were run in point mode a Ameriflux sites, US-MMS, US-Moz, US-Syv, and US-NR1. These datasets can be used to assess the impact that dynamic roots have on Land Surface fluxes such as carbon and water. The simulations were run with atmospheric forcing from the Global Soil Wetness Project. We provide specifically here Gross Primary Production, Transpiration, relative root fraction per soil layer as well as Air Temperature, Precipitation and all the parameterizations used for the simulations. Our purpose for generating and analyzing these simulations was to assess the role that root dynamics and foraging for water has in the recovery timescale of ecosystems following climate perturbation. All data provided are in the "netcdf" format using standard CF-1 (Climate and Forecast) Convention. The data are all machine readable using various software platforms including Matlab, NCO, Panoply and R.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Dataset for scientific paper "Simulated plant‑mediated oxygen input has strong impacts on fine‑scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands", a modeling study based on field observation at the tidal salt marshes of the Parker River Estuary, Massachusetts, United States

This dataset is the raw and processed data for the paper "Simulated plant ‑ mediated oxygen input has strong impacts on fine ‑ scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands". This study investigated how plant-mediated oxygen input affects subsurface biogeochemical reactions of organic carbon degradation and the resulting methane emissions of coastal wetlands by model simulation. We used the subsurface geochemical simulator PFLOTRAN for the modeling, which produced the simulated changes in porewater chemical substances and methane emissions over 10 days under different scenarios of plant-mediated oxygen input.Specifically, this dataset contains: 1) the input files for PFLOTRAN of all simulation runs conducted in this study. Those files are with an extension of ".in", containing information of the biogeochemical reaction network (stoichiometry, reaction rate, Monod constants, etc), fluid flow rate and oxygen concentration in the fluid which together simulated the plant-mediated oxygen input, the configuration of artificial reactions that simulated the methane fluxes, etc. The PFLOTRAN input files are text files, which can be opened by NotePad, but running these input files will require proper installation of PFLOTRAN (instruction: https://documentation.pflotran.org/user_guide/how_to/installation/installation.html). 2) the raw and processed model output from PFLOTRAN of all simulation runs, and 3) the python scripts used to process the raw model output, including random allocation of root cells, converting raw data into organized formats, calculating the methane fluxes based on the model output, data visualization, etc. The raw and processed model output from PFLOTRAN are in .spydata format, which can be viewed with Python. and 3) the python scripts for data processing and analysis are programming scripts, which can be opened with Python.This modeling work, in particular the model parameterization of root density and initial conditions of porewater concentrations of biogeochemical substances, was based on field measurements at the salt marsh of the Upper Parker River Estuary, Massachusetts, United States.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Oil Recovery Prediction for Polymer Flood Field Test of Heavy Oil on Alaska North Slope Via Machine Assisted Reservoir Simulation

The first ever polymer flood field pilot to enhance the recovery of heavy oils on the Alaska North Slope is ongoing. This study constructs and calibrates a reservoir simulation model to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding, as indicated by a significant decrease in water cut. Later, transmissibility contrasts are reinstated to replicate a water surge event observed in one of the producing wells during polymer flooding. This event may represent decreased injection conformance from fracture overextension; its anticipated occurrence in the other production well is included in the final forecast. The definition of polymer retention in the simulator incorporates the tailing effect reported in laboratory studies; this tailing effect is useful to the simultaneous history match of producing water cut and produced polymer concentration. The top 24 best-matched simulation models produced at each stage of the history matching process are used to forecast oil recovery. The final forecast clearly demonstrates that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone. This exercise displays that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the benefits of an enhanced oil recovery process such as polymer flooding. Therefore, caution is recommended in using basic waterflood simulation models to scope potential enhanced oil recovery projects.

Keith, Cody Douglas↗