SEARCH · Engineering Papers

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,639 records · Page 91

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Transient Cold Flow Simulation of Fast Fluidized Bed Fuel Reactors for Chemical-Looping Combustion

Circulating fluidized bed in chemical-looping combustion (CLC) is a recent technology that provides great advantage for gas–solid interaction and efficiency. In order to obtain a thorough understanding of this technology and to assess its effectiveness for industrial scale deployment, numerical simulations are conducted. Computational fluid dynamics (CFD) simulations are performed with dense discrete phase model (DDPM) to simulate the gas–solid interactions. CFD commercial software ansysfluent is used for the simulations. Two bed materials of different particle density and diameter, namely the molochite and Fe100, are used in studying the hydrodynamics and particle behavior in a fuel reactor corresponding to the experimental setup of Haider et al. (2016, “A Hydrodynamic Study of a Fast-Bed Dual Circulating Fluidized Bed for Chemical Looping Combustion,” Energy Technol., 4(10), pp. 1254–1262.) at Cranfield University in the UK. Both the simulations show satisfactory agreement with the experimental data for both the static pressure and volume fraction at various heights above the gas inlet in the reactor. It is found that an appropriate drag law should be used in the simulation depending on the particle size and flow conditions in order to obtain accurate results. Furthermore, the simulations demonstrate the ability of CFD/DDPM to accurately capture the physics of circulating fluidized bed-based CLC process at pilot scale which can be extended to industrial scale projects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2019

Predicting Cycle-to-Cycle Variation With Concurrent Cycles in a Gasoline Direct Injected Engine With Large Eddy Simulations

High cycle-to-cycle variation (CCV) is detrimental to engine performance, as it leads to poor combustion and high noise and vibration. In this work, CCV in a gasoline engine is studied using large eddy simulation (LES). The engine chosen as the basis of this work is a single-cylinder gasoline direct injection (GDI) research engine. Two stoichiometric part-load engine operating points (6 BMEP, 2000 RPM) were evaluated: a non-dilute (0% EGR) case and a dilute (18% EGR) case. The experimental data for both operating conditions had 500 cycles. The measured CCV in IMEP was 1.40% for the non-dilute case and 7.78% for the dilute case. To estimate CCV from simulation, perturbed concurrent cycles of engine simulations were compared to consecutively obtained engine cycles. The motivation behind this is that running consecutive cycles to estimate CCV is quite time-consuming. For example, running 100 consecutive cycles requires 2–3 months (on a typical cluster), however, by running concurrently one can potentially run all 100 cycles at the same time and reduce the overall turnaround time for 100 cycles to the time taken for a single cycle (2 days). The goal of this paper is to statistically determine if concurrent cycles, with a perturbation applied to each individual cycle at the start, can be representative of consecutively obtained cycles and accurately estimate CCV. 100 cycles were run for each case to obtain statistically valid results. The concurrent cycles began at different timings before the combustion event, with the motivation to identify the closest time before spark to minimize the run time. Only a single combustion cycle was run for each concurrent case. The calculated standard deviation of peak pressure and coefficient of variance (COV) of indicated mean effective pressure (IMEP) were compared between the consecutive and concurrent methods to quantify CCV. It was found that the concurrent method could be used to predict CCV with either a velocity or numerical perturbation. Here, a large and small velocity perturbation were compared and both produced correct predictions, implying that the type of perturbation is not important to yield a valid realization. Starting the simulation too close to the combustion event, at intake valve close (IVC) or at spark timing, under-predicted the CCV. When concurrent simulations were initiated during or before the intake even, at start of injection (SOI) or earlier, distinct and valid realizations were obtained to accurately predict CCV for both operating points. By simulating CCV with concurrent cycles, the required wall clock time can be reduced from 2–3 months to 1–2 days. Additionally, the required core-hours can be reduced up to 41%, since only a portion of each cycle needs to be simulated.

42 ENGINEERING↗

DOE OSTI2014

TH-E-BRE-09: TrueBeam Monte Carlo Absolute Dose Calculations Using Monitor Chamber Backscatter Simulations and Linac-Logged Target Current

Purpose: To simulate and measure radiation backscattered into the monitor chamber of a TrueBeam linac; establish a rigorous framework for absolute dose calculations for TrueBeam Monte Carlo (MC) simulations through a novel approach, taking into account the backscattered radiation and the actual machine output during beam delivery; improve agreement between measured and simulated relative output factors. Methods: The ‘monitor backscatter factor’ is an essential ingredient of a well-established MC absolute dose formalism (the MC equivalent of the TG-51 protocol). This quantity was determined for the 6 MV, 6X FFF, and 10X FFF beams by two independent Methods: (1) MC simulations in the monitor chamber of the TrueBeam linac; (2) linac-generated beam record data for target current, logged for each beam delivery. Upper head MC simulations used a freelyavailable manufacturer-provided interface to a cloud-based platform, allowing use of the same head model as that used to generate the publicly-available TrueBeam phase spaces, without revealing the upper head design. The MC absolute dose formalism was expanded to allow direct use of target current data. Results: The relation between backscatter, number of electrons incident on the target for one monitor unit, and MC absolute dose was analyzed for open fields, as well as a jaw-tracking VMAT plan. The agreement between the two methods was better than 0.15%. It was demonstrated that the agreement between measured and simulated relative output factors improves across all field sizes when backscatter is taken into account. Conclusion: For the first time, simulated monitor chamber dose and measured target current for an actual TrueBeam linac were incorporated in the MC absolute dose formalism. In conjunction with the use of MC inputs generated from post-delivery trajectory-log files, the present method allows accurate MC dose calculations, without resorting to any of the simplifying assumptions previously made in the TrueBeam MC literature. This work has been partially funded by Varian Medical Systems.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2025

Determining Levels of Detail for Simulators of Parallel and Distributed Computing Systems via Automated Calibration

There are two sources of inaccuracy when simulating parallel and distributed computing systems: (i) a simulator implemented at an insufficient level of detail; and (ii) incorrectly calibrated simulation parameter values. Increasing the simulator’s level of detail can improve accuracy, but at the cost of higher space, time, and/or software complexity. Furthermore, evaluating the intrinsic accuracy of a simulator requires that its parameters be well-calibrated. Making decisions regarding the level of detail is thus challenging. We propose a methodology for instantiating the simulation calibration process and a framework for automating this process, which makes it possible to pick appropriate levels of detail for any simulator. We demonstrate the usefulness of our approach via two case studies for two different domains.

McDonald, Jessie [University of Hawaii at Manoa, H↗

DOE OSTI2014

Framework for Network Co-Simulation

The Framework for Network Co-Simulation (FNCS) uses a federated approach to integrate simulations which may have differing time scales. Special consideration is given to integration with a communication network simulation such that inter-simulation messages may be optionally routed through and delayed by such a simulation. In addition, FNCS uses novel time synchronization algorithms to accelerate co-simulation including the application of speculative multithreading. FNCS accomplishes all of these improvements with minimal end user intervention. Simulations can be integrated using FNCS while maintaining their original model input files simply by linking with the FNCS library and making appropriate calls into the FNCS API.

Daily, Jeff↗

DOE OSTI2015

Parametric Sensitivity Analysis for the Asian Summer Monsoon Precipitation Simulation in the Beijing Climate Center AGCM Version 2.1

In this study, we apply an efficient sampling approach and conduct a large number of simulations to explore the sensitivity of the simulated Asian summer monsoon (ASM) precipitation, including the climatological state and interannual variability, to eight parameters related to the cloud and precipitation processes in the Beijing Climate Center AGCM version 2.1 (BCC_AGCM2.1). Our results show that BCC_AGCM2.1 has large biases in simulating the ASM precipitation. The precipitation efficiency and evaporation coefficient for deep convection are the most sensitive parameters in simulating the ASM precipitation. With optimal parameter values, the simulated precipitation climatology could be remarkably improved, e.g. increased precipitation over the equator Indian Ocean, suppressed precipitation over the Philippine Sea, and more realistic Meiyu distribution over Eastern China. The ASM precipitation interannual variability is further analyzed, with a focus on the ENSO impacts. It shows the simulations with better ASM precipitation climatology can also produce more realistic precipitation anomalies during El Niño decaying summer. In the low-skill experiments for precipitation climatology, the ENSO-induced precipitation anomalies are most significant over continents (vs. over ocean in observation) in the South Asian monsoon region. More realistic results are derived from the higher-skill experiments with stronger anomalies over the Indian Ocean and weaker anomalies over India and the western Pacific, favoring more evident easterly anomalies forced by the tropical Indian Ocean warming and stronger Indian Ocean-western Pacific tele-connection as observed. Our model results reveal a strong connection between the simulated ASM precipitation climatological state and interannual variability in BCC_AGCM2.1 when key parameters are perturbed.

Asian summer monsoon precipitation↗

DOE OSTI2020

Evaluation of Mesoscale Convective Systems in Climate Simulations: Methodological Development and Results from MPAS-CAM over the U.S.

In this study, a process-oriented approach is developed to evaluate warm-season mesoscale convective system (MCS) precipitation and their favorable large-scale meteorological patterns (FLSMPs) over the U.S. This approach features a novel observation-driven MCS-tracking algorithm using infrared brightness temperature and precipitation feature at 12, 25 and 50 km resolution and metrics to evaluate the model large-scale environment favorable for MCS initiation. The tracking algorithm successfully reproduces the observed MCS statistics from a reference 4-km radar MCS database. To demonstrate the utility of the new methodologies in evaluating MCS in climate simulations with mesoscale resolution, the process-oriented approach is applied to two climate simulations produced by the Variable-Resolution Model for Prediction Across Scales coupled to the Community Atmosphere Model physics, with refined horizontal grid spacing at 50 km and 25 km over North America. With the tracking algorithm applied to simulations and observations at equivalent resolutions, the simulated number of MCS and associated precipitation amount, frequency and intensity are found to be consistently underestimated in the Central U.S., particularly from May to August. The simulated MCS precipitation shows little diurnal variation and lasts too long, while MCS precipitation area is too large and intensity is too weak. The model is able to simulate four types of observed FLSMP associated with frontal systems and low-level jets (LLJ) in spring, but the frequencies are underestimated because of low-level dry bias and weaker LLJ. Precipitation simulated under different FLSMPs peak during daytime, in contrast to the observed nocturnal peak. Implications of these findings for future model development and diagnostics are discussed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Assessment of CMIP5 and CMIP6 AMIP Simulated Clouds and Surface Shortwave Radiation Using ARM Observations over Different Climate Regions

The simulations of clouds and surface radiation from 10 CMIP6 models and their CMIP5 predecessors are compared to the ARM ground-based observations over different climate regions. Compared to the ARM radar-lidar derived total cloud fractions (CFT) and cloud fraction vertical distributions over the six selected sites, both CMIP5 and CMIP6 significantly underestimated CFT and low-level CF over the Northern Hemispheric midlatitude sites (SGPC1 and ENAC1), although the biases are generally smaller in CMIP6. Over the tropical oceanic site (TWPC2), 5 out of 10 CMIP6 models better simulated low-level CF than their CMIP5 predecessors. CMIP6 simulations generally agreed well with the ARM observations in CFT and cloud fraction vertical distributions over the tropical continental (MAOM1) and coastal (TWPC3) sites but missed the transitions between dry and wet seasons, similar to CMIP5 simulations. The improvements in downwelling shortwave fluxes (SWdn) at the surface from the majority of CMIP6 compared to CMIP5 primarily resulted from the improved cloud fraction simulations, especially over the SGPC1, ENAC1, and TWPC3 sites. By contrast, both CMIP5 and CMIP6 models exhibited diverse performances of clouds and shortwave radiation over the Arctic site (NSAC1), where CMIP6 models produced more clouds than CMIP5 models, especially for the low-level clouds. The comparisons between observations and CMIP5 and CMIP6 simulations provide valuable quantitative assessments of the accuracy of mean states and variabilities in the model simulations and shed light on general directions to improve climate models in different regions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Understanding the Biases in Global Monsoon Simulations from the Perspective of Atmospheric Energy Transport

Understanding global monsoon (GM) variability and projecting its future changes rely heavily on climate models. However, climate models generally show pronounced biases in GM simulations, and the reasons for this remain unclear. Here, in this study, we evaluate the performance of 20 pairs of climate models that participated in both phase 5 of the Coupled Model Intercomparison Project (CMIP5) and phase 6 of CMIP (CMIP6) and identify the sources of their GM simulation biases from an energy transport perspective. The multimodel mean improvement in CMIP6 compared to CMIP5 is demonstrated by the increasing skill scores for various GM metrics from 0.20–0.79 to 0.48–0.83. More specifically, the dry biases in the Northern Hemisphere Summer Monsoon (NHSM) precipitation in CMIP5 [root-mean-square error (RMSE): 1.85 mm day −1 ] are reduced in CMIP6 (RMSE: 1.66 mm day −1 ). This higher simulation skill is associated with higher skill in simulating the precipitation-solstitial mode, monsoon intensity, and monsoon domains. The improvement in the NHSM precipitation simulation results from that in the meridional transport of atmospheric energy. Atmospheric energy budget analysis shows that the negative biases in downward surface longwave radiation and northward energy transport are smaller in CMIP6 than in CMIP5 in the boreal summer, resulting in a more realistic interhemispheric thermal contrast and meridional gradient of moist static energy. However, a major weakness of the CMIP6 models is found in the Southern Hemisphere Summer Monsoon precipitation simulation due to the positive bias in the top-of-the-atmosphere downward longwave radiation. This study shows that reasonably reproducing the meridional global atmospheric energy transportation is necessary for skillful GM simulation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

VERA Core Simulator methodology for pressurized water reactor cycle depletion

This paper describes the methodology developed and implemented in the Virtual Environment for Reactor Applications Core Simulator (VERA-CS) to perform high-fidelity, pressurized water reactor (PWR), multicycle, core physics calculations. Depletion of the core with pin-resolved power and nuclide detail is a significant advance in the state of the art for reactor analysis, providing the level of detail necessary to address the problems of the U.S. Department of Energy Nuclear Reactor Simulation Hub, the Consortium for Advanced Simulation of Light Water Reactors (CASL). VERA-CS has three main components: the neutronics solver MPACT, the thermal-hydraulic (T-H) solver COBRA-TF (CTF), and the nuclide transmutation solver ORIGEN. This paper focuses on MPACT and provides an overview of the resonance self-shielding methods, macroscopic-cross-section calculation, two-dimensional/one-dimensional (2-D/1-D) transport, nuclide depletion, T-H feedback, and other supporting methods representing a minimal set of the capabilities needed to simulate high-fidelity models of a commercial nuclear reactor. Results are presented from the simulation of a model of the first cycle of Watts Bar Unit 1. The simulation is within 16 parts per million boron (ppmB) reactivity for all state points compared to cycle measurements, with an average reactivity bias of <5 ppmB for the entire cycle. Comparisons to cycle 1 flux map data are also provided, and the average 2-D root-mean-square (rms) error during cycle 1 is 1.07%. To demonstrate the multicycle capability, a state point at beginning of cycle (BOC) 2 was also simulated and compared to plant data. The comparison of the cycle 2 BOC state has a reactivity difference of +3 ppmB from measurement, and the 2-D rms of the comparison in the flux maps is 1.77%. Lastly, these results provide confidence in VERA-CS’s capability to perform high-fidelity calculations for practical PWR reactor problems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Quantitative comparison between sub-millisecond time resolution single-molecule FRET measurements and 10-second molecular simulations of a biosensor protein

Molecular Dynamics (MD) simulations seek to provide atomic-level insights into conformationally dynamic biological systems at experimentally relevant time resolutions, such as those afforded by single-molecule fluorescence measurements. However, limitations in the time scales of MD simulations and the time resolution of single-molecule measurements have challenged efforts to obtain overlapping temporal regimes required for close quantitative comparisons. Achieving such overlap has the potential to provide novel theories, hypotheses, and interpretations that can inform idealized experimental designs that maximize the detection of the desired reaction coordinate. Here, we report MD simulations at time scales overlapping with in vitro single-molecule Förster (fluorescence) resonance energy transfer (smFRET) measurements of the amino acid binding protein LIV-BP SS at sub-millisecond resolution. Computationally efficient all-atom structure-based simulations, calibrated against explicit solvent simulations, were employed for sampling multiple cycles of LIV-BP SS clamshell-like conformational changes on the time scale of seconds, examining the relationship between these events and those observed by smFRET. The MD simulations agree with the smFRET measurements and provide valuable information on local dynamics of fluorophores at their sites of attachment on LIV-BP SS and the correlations between fluorophore motions and large-scale conformational changes between LIV-BP SS domains. We further utilize the MD simulations to inform the interpretation of smFRET data, including Förster radius (R 0 ) and fluorophore orientation factor (κ 2 ) determinations. The approach we describe can be readily extended to distinct biochemical systems, allowing for the interpretation of any FRET system conjugated to protein or ribonucleoprotein complexes, including those with more conformational processes, as well as those implementing multi-color smFRET.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2013

Assessment of Land Surface Model Performance in WRF for Simulating Wind at Heights Relevant to the Wind Energy Community

The Weather Research and Forecasting (WRF) model was used to investigate the impact of land surface model (LSM) physics on the simulated surface energy balance and the near surface wind profile, including heights reached by a modern wind turbine. Simulations were made for a domain centered over the Department of Energy (DOE) Southern Great Plains (SGP) Atmospheric Radiation Measurement Program (ARM) Central Facility in northern Oklahoma. High-resolution measurements of surface fluxes and profiles of wind speed, direction and turbulence were available for model validation. We tested 10 common LSM configurations in WRF including Noah, Noah-MP, SSiB, Pleim-Xiu, RUC, and others to evaluate (1) the accuracy of simulated surface energy fluxes, (2) the accuracy of forecasted rotor-disk height wind speeds, and (3) the sensitivity of rotor disk wind speed to the choice of LSM. WRF was run for four, two-week periods covering different canopy, surface meteorology and soil moisture conditions. LSM performances varied by simulation period and errors were greatest during periods with high soil moisture and live crop cover, indicating highest LSM sensitivity when PBL behavior is most coupled to land-atmosphere exchange. Higher wind speed errors also occurred at night during the spring, possibly from the difficulty with simulating nocturnal low-level jets. We found that choice of land surface model led to a ~10% improvement in simulating hub-height wind speed. Errors in wind shear were also sensitive to the choice of LSM. Overall the best performing models were Noah and Noah-MP. The variability of LSM performance across different soil moisture and vegetation canopy conditions suggests LSM representation of surface energy exchange processes in WRF remains a large source of model uncertainty. Furthermore, our results indicate the choice of LSM impacts simulated wind fields in WRF at heights relevant to the wind industry, at least for the Southern Great Plains region.

17 WIND ENERGY↗

DOE OSTI2013

RELAP-7 Simulation Resolving an SBO Scenario on a Simplified Geometry of a BWR

The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 will become the main reactor systems toolkit for the Risk-Informed Safety Margin Characterization Pathway of the Light Water Reactor Sustainability Program and the next generation tool in the RELAP reactor safety/systems analysis application series (i.e., the replacement for RELAP5). The code is being developed based on Idaho National Laboratory’s modern scientific software development framework – MOOSE (the Multi-Physics Object-Oriented Simulation Environment). During Fiscal Year 2013, a number of physical components with two-phase flow capability have been developed to support the simplified boiling water reactor (BWR) station blackout analyses. The case selected for demonstration calculation is built from the specifications documented in an Organization for Economic Cooperation and Development benchmark problem for BWR turbine trip analysis. The reference design for the benchmark problem was from the Peach Bottom-2 nuclear station, which is a General Electric BWR-4 design. The demonstration case includes the major components for the primary system of a BWR, as well as the safety system components for reactor core isolation cooling (RCIC) and the wet well of a BWR containment. The case was initially run to steady-state with RELAP-7. The station blackout transient simulations were subsequently initiated by using the INL-developed RAVEN code. Two scenarios for the station blackout simulations have been considered. Scenario I represents an extreme station blackout accident with no safety injection functioning. The reactor core could experience dry out fairly quickly with this scenario. Scenario II represents a more probable station blackout accident progression with the RCIC system functioning. In this scenario, the RCIC system is fully coupled with the reactor primary system and the safety injection to provide makeup cooling water to the reactor core from the suppression pool is dynamically simulated. With the RCIC system functioning, the core dry out is significantly postponed when compared to the results from Scenario I. This fully coupled RCIC system simulation capability represents the first-of-a-kind simulation capability. Sensitivity studies have also been carried out to study the effect of RCIC control strategy with varying core makeup cooling water mass flow rates. The next stage of development will be to demonstrate more refined BWR station blackout analyses with more realistic geometries, to be reported in the next demonstration simulation report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Simulations of Failure via Three-Dimensional Cracking in Fuel Cladding for Advanced Nuclear Fuels

Enhancing performance of fuel cladding and duct alloys is a key means of increasing fuel burnup. This project will address the failure of fuel cladding via three-dimensional cracking models. Researchers will develop a simulation code for the failure of the fuel cladding and validate the code through experiments. The objective is to develop an algorithm to determine the failure of fuel cladding in the form of three-dimensional cracking due to prolonged exposure under varying conditions of pressure, temperature, chemical environment, and irradiation. This project encompasses the following tasks: 1. Simulate 3D crack initiation and growth under instantaneous and/or fatigue loads using a new variant of the material point method (MPM); 2. Simulate debonding of the materials in the crack path using cohesive elements, considering normal and shear traction separation laws; 3. Determine the crack propagation path, considering damage of the materials incorporated in the cohesive elements to allow the energy release rate to be minimized; 4. Simulate the three-dimensional fatigue crack growth as a function of loading histories; 5. Verify the simulation code by comparing results to theoretical and numerical studies available in the literature; 6. Conduct experiments to observe the crack path and surface profile in unused fuel cladding and validate against simulation results; and 7. Expand the adaptive mesh refinement infrastructure parallel processing environment to allow adaptive mesh refinement at the 3D crack fronts and adaptive mesh merging in the wake of cracks. Fuel cladding is made of materials such as stainless steels and ferritic steels with added alloying elements, which increase stability and durability under irradiation. As fuel cladding is subjected to water, chemicals, fission gas, pressure, high temperatures, and irradiation while in service, understanding performance is essential. In the fast fuel used in advanced burner reactors, simulations of the nuclear fuels are critical to understand the burnup, and thus the fuel efficiency.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2013

Advanced Simulation and Computing Fiscal Year 14 Implementation Plan, Rev. 0.5

The Stockpile Stewardship Program (SSP) is a single, highly integrated technical program for maintaining the surety and reliability of the U.S. nuclear stockpile. The SSP uses nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program requires the continued use of experimental facilities and programs, and the computational enhancements to support these programs. The Advanced Simulation and Computing Program (ASC) is a cornerstone of the SSP, providing simulation capabilities and computational resources that support annual stockpile assessment and certification, study advanced nuclear weapons design and manufacturing processes, analyze accident scenarios and weapons aging, and provide the tools to enable stockpile Life Extension Programs (LEPs) and the resolution of Significant Finding Investigations (SFIs). This requires a balanced resource, including technical staff, hardware, simulation software, and computer science solutions. In its first decade, the ASC strategy focused on demonstrating simulation capabilities of unprecedented scale in three spatial dimensions. In its second decade, ASC is now focused on increasing predictive capabilities in a three-dimensional (3D) simulation environment while maintaining support to the SSP. The program continues to improve its unique tools for solving progressively more difficult stockpile problems (sufficient resolution, dimensionality, and scientific details), quantify critical margins and uncertainties, and resolve increasingly difficult analyses needed for the SSP. Moreover, ASC’s business model is integrated and focused on requirements-driven products that address long-standing technical questions related to enhanced predictive capability in the simulation tools.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Development of an Advanced Simulator to Model Mobility Control and Geomechanics during CO 2 Floods

The simulator is an isothermal, three-dimensional, four-phase, compositional, equation-of–state (EOS) simulator. We have named the simulator UTDOE-CO2 capable of simulating various recovery processes (i.e., primary, secondary waterflooding, and miscible and immiscible gas flooding). We include both the Peng-Robinson EOS and the Redlich-Kwong EOS models. A Gibbs stability test is also included in the model to perform a phase identification test to consistently label each phase for subsequent property calculations such as relative permeability, viscosity, density, interfacial tension, and capillary pressure. Our time step strategy is based on an IMPEC-type method (implicit pressure and explicit concentration). The gridblock pressure is solved first using the explicit dating of saturation-dependent terms. Subsequently, the material balance equations are solved explicitly for the total concentration of each component. The physical dispersion term is also included in the governing equations. The simulator includes (1) several foam model(s) for gas mobility control, (2) compositional relative permeability models with the hysteresis option, (3) corner point grid and several efficient solvers, (4) geomechanics module to compute stress field as the result of CO 2 injection/production, (5) the format of commercial visualization software, S3graf from Science-soft Ltd., was implemented for user friendly visualization of the simulation results. All tasks are completed and the simulator was fully tested and delivered to the DOE office including a user’s guide and several input files and the executable for Windows Pcs. We have published several SPE papers, presented several posters, and one MS thesis is completed (V. Pudugramam, 2013) resulting from this DOE funded project.

02 PETROLEUM↗

DOE OSTI2015

Fast Dynamic Simulation-Based Small Signal Stability Assessment and Control

Power grid planning and operation decisions are made based on simulation of the dynamic behavior of the system. Enabling substantial energy savings while increasing the reliability of the aging North American power grid through improved utilization of existing transmission assets hinges on the adoption of wide-area measurement systems (WAMS) for power system stabilization. However, adoption of WAMS alone will not suffice if the power system is to reach its full entitlement in stability and reliability. It is necessary to enhance predictability with "faster than real-time" dynamic simulations that will enable the dynamic stability margins, proactive real-time control, and improve grid resiliency to fast time-scale phenomena such as cascading network failures. Present-day dynamic simulations are performed only during offline planning studies, considering only worst case conditions such as summer peak, winter peak days, etc. With widespread deployment of renewable generation, controllable loads, energy storage devices and plug-in hybrid electric vehicles expected in the near future and greater integration of cyber infrastructure (communications, computation and control), monitoring and controlling the dynamic performance of the grid in real-time would become increasingly important. The state-of-the-art dynamic simulation tools have limited computational speed and are not suitable for real-time applications, given the large set of contingency conditions to be evaluated. These tools are optimized for best performance of single-processor computers, but the simulation is still several times slower than real-time due to its computational complexity. With recent significant advances in numerical methods and computational hardware, the expectations have been rising towards more efficient and faster techniques to be implemented in power system simulators. This is a natural expectation, given that the core solution algorithms of most commercial simulators were developed decades ago, when High Performance Computing (HPC) resources were not commonly available.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2015

High Fidelity Simulations of Large-Scale Wireless Networks (Plus-Up)

Sandia has built a strong reputation in scalable network simulation and emulation for cyber security studies to protect our nation’s critical information infrastructures. Georgia Tech has preeminent reputation in academia for excellence in scalable discrete event simulations, with strong emphasis on simulating cyber networks. Many of the experts in this field, such as Dr. Richard Fujimoto, Dr. George Riley, and Dr. Chris Carothers, have strong affiliations with Georgia Tech. The collaborative relationship that we intend to immediately pursue is in high fidelity simulations of practical large-scale wireless networks using ns-3 simulator via Dr. George Riley. This project will have mutual benefits in bolstering both institutions’ expertise and reputation in the field of scalable simulation for cyber-security studies. This project promises to address high fidelity simulations of large-scale wireless networks. This proposed collaboration is directly in line with Georgia Tech’s goals for developing and expanding the Communications Systems Center, the Georgia Tech Broadband Institute, and Georgia Tech Information Security Center along with its yearly Emerging Cyber Threats Report. At Sandia, this work benefits the defense systems and assessment area with promise for large-scale assessment of cyber security needs and vulnerabilities of our nation’s critical cyber infrastructures exposed to wireless communications.

Onunkwo, Uzoma↗