SEARCH · Engineering Papers
Results for “Reactive Flow”
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.
PyFLOTRAN Documentation Release 1.0.0
This is the documentation for the PyFLOTRAN library. PyFLOTRAN is a set of classes and methods that the enable use of the massively parallel subsurface flow and reactive transport code PFLOTRAN (http://www.pflotran.org) within the Python scripting environment. This allows the use of a wide range of libraries available in Python for pre- and post-processing.
2016 KIVA-hpFE Development. A Robust and Accurate Engine Modeling Software
Everything we are doing in R&D is to develop methods and a code for robust, accurate and efficient algorithms in a Parallel (MPI) Modular Object-Oriented code for Industry and Researchers.
Physical, Hydraulic, and Transport Properties of Sediments and Engineered Materials Associated with Hanford Immobilized Low-Activity Waste
Current plans for treatment and disposal of immobilized low-activity waste (ILAW) from Hanford’s underground waste storage tanks include vitrification and storage of the glass waste form in a nearsurface disposal facility. This Integrated Disposal Facility (IDF) is located in the 200 East Area of the Hanford Central Plateau. Performance assessment (PA) of the IDF requires numerical modeling of subsurface flow and reactive transport processes over very long periods (thousands of years). The models used to predict facility performance require parameters describing various physical, hydraulic, and transport properties. This report provides updated estimates of physical, hydraulic, and transport properties and parameters for both near- and far-field materials, intended for use in future IDF PA modeling efforts. Previous work on physical and hydraulic property characterization for earlier IDF PA analyses is reviewed and summarized. For near-field materials, portions of this document and parameter estimates are taken from an earlier data package. For far-field materials, a critical review is provided of methodologies used in previous data packages. Alternative methods are described and associated parameters are provided.
FEARCE Development - Institutional Computing Report for 2019
This report details the development, uses ,and methods of FEARCE (Fast, Easy, Accurate, and Robust Continuum Engineering) software.
Remote Off-grid Microgrid Design Support Tool (ROMDST) - An Optimal Design Support Tool for Remote, Resilient, and Reliable Microgrids Phase II - Final Report
The objective of the project is to develop an advanced optimization-based design support tool for AC or DC microgrids in remote locations, where utility grids may not be accessible. The mathematical model and the interface are being developed such that multiple design objectives and criteria/constraints can be easily enabled or disabled, to deliver a flexible tool relevant to a large, diverse user base, and to facilitate future feature developments. The new tool delivered by this project leverages the team’s extensive expertise in the development, testing, deployment, and commercialization of the state-of-the-art Distributed Energy Resources Customer Adoption Model (DER-CAM) – which is the foundation for the new tool. Phase I of the project leverages the microgrid design and decision support tool DER-CAM. Key Phase I developments include the following adaption to the DER-CAM model: a) AC or DC microgrid architectures, b) multiple economic objectives (e.g. single customer vs community vs utility perspective), c) constraints on fuel availability (e.g. total fuel available during outage periods), d) active and reactive power flow constraints for normal and N-1 contingency cases, e) component part load efficiencies (e.g. power electronic devices), and f) interactive data visualization capabilities (e.g. design topology elements). A component library including components’ default technical and economical characteristics was also developed. In Phase II, the project is focused on tool testing and validation, use for real-world microgrid designs, and outreach. Fully developing this tool will have the following impacts: 1) optimal off-grid microgrid designs that replace existing back of the envelope or non-optimal calculations, reducing capital costs and risk of microgrid deployment; 2) removing barriers to microgrid assessments by lowering microgrid soft costs, as the tool is freely usable; 3) reliable and resilient designs that reduce the cost of critical load shedding due to component outages.
Remote Off-grid Microgrid Design Support Tool (ROMDST) - An Optimal Design Support Tool for Remote, Resilient, and Reliable Microgrids (Phase II, Final Report)
The objective of the project is to develop an advanced optimization-based design support tool for AC or DC microgrids in remote locations, where utility grids may not be accessible. The mathematical model and the interface are being developed such that multiple design objectives and criteria/constraints can be easily enabled or disabled, to deliver a flexible tool relevant to a large, diverse user base, and to facilitate future feature developments. The new tool delivered by this project leverages the team’s extensive expertise in the development, testing, deployment, and commercialization of the state-of-the-art Distributed Energy Resources Customer Adoption Model (DER-CAM) – which is the foundation for the new tool. Phase I of the project leverages the microgrid design and decision support tool DER-CAM. Key Phase I developments include the following adaption to the DER-CAM model: a) AC or DC microgrid architectures, b) multiple economic objectives (e.g. single customer vs community vs utility perspective), c) constraints on fuel availability (e.g. total fuel available during outage periods), d) active and reactive power flow constraints for normal and N-1 contingency cases, e) component part load efficiencies (e.g. power electronic devices), and f) interactive data visualization capabilities (e.g. design topology elements). A component library including components’ default technical and economical characteristics was also developed. In Phase II, the project is focused on tool testing and validation, use for real-world microgrid designs, and outreach. Fully developing this tool will have the following impacts: 1) optimal off-grid microgrid designs that replace existing back of the envelope or non-optimal calculations, reducing capital costs and risk of microgrid deployment; 2) removing barriers to microgrid assessments by lowering microgrid soft costs, as the tool is freely usable; 3) reliable and resilient designs that reduce the cost of critical load shedding due to component outages.
High Voltage Modeling and Testing of Transformer, Line Interface Devices, and Bulk System Components Under Electromagnetic Pulse, Geomagnetic Disturbance, and other Abnormal Transients
This project aims at quantitatively analyzing possible knowledge gaps related to the risk of Electromagnetic Pulse (EMP) impact on high-voltage, power transformers (transmission-class), that represent the most critical power grid assets. In relation to the electric infrastructure, some past studies indicate that the grid is overall fairly robust and resilient against most EMP threats. This assessment stems, in part, from estimates of the peak voltage that can be induced on a transmission power line for a nominal EMP fast-transient, compared with the typically very high basic insulation level that is featured by the high-voltage infrastructure. The project provides additional insight toward a more comprehensive risk evaluation of the grid EMP and Geomagnetic Disturbances (GMD) vulnerabilities. This work is aligned with the U.S. Dept. of Energy current research thrust in support of the development of a highly secure and resilient power infrastructure. Project milestones have been devised to address the fundamental technical issues in support of protection of large power transformer against EMP threats. More specifically, this project included parallel research tasks focused on the following topics: analysis of EMP-induced, fast-transient, current and voltages on a variety of power line configurations and EMP impact scenarios, detailed characterization of EMP damage potential due to the voltage surge, from the impact on the power line, to the transformer windings, GMD and EMP slow transients and impact from harmonics generation and reactive power flow on a complex grid layout, and identification of risk mitigation paths, with emphasis on GIC damage prevention for large power transformers, and on impact protection options from a grid stability perspective.
EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN
This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.
FY20 SAM Code Developments and Validations for Transient Safety Analysis of Advanced non-LWRs
The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report summarizes major progress in SAM code development, capability enhancements, demonstration, and validation to support transient safety analysis of advanced non-LWRs. Rapid developments continued in fiscal year 2020 (FY20) to support various needs of the advanced reactor community, especially the NRC and industry on the licensing safety analysis of advanced reactor designs. Significant code changes were made to provide various capability enhancements, bug fixes, and user friendliness improvements. Major code updates are summarized in Section 1, while four important enhancements are detailed in Sections 2-5, including a multi-dimension flow model; reactivity feedback and decay heat models; control and trip system modeling, and additional fluid and solid thermophysical property models. Code validation activities in FY20 include using test data from the Fast Flux Test Facility (FFTF), the High Temperature Test Facility (HTTF), and several separate effects test facilities for pebble-bed modeling.
Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems
An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.
Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code
CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.
Synopsis of SURF and SURFplus
Developed by Menikoff and Shaw, the SURF reactive burn model builds on the Ignition and Growth concept by incorporating the lead shock pressure directly into the volumetric hot-spot burn rate. The plus-extension augments the model with a late time surface burn rate due to carbon clustering. Together, SURF and SURFplus allow for robust modeling of both conventional and insensitive high explosives. Practically, the SURF burn models require shock detection and the advection of the lead shock pressure as an additional material field. Calibrated parameters for given equations of state and thermodynamic closure are dependent on initial temperature and density. An implementation in Python is given.
Modeling commercial-scale CO 2 storage in the gas hydrate stability zone with PFLOTRAN v6.0
Abstract. Safe and secure carbon dioxide (CO2) storage is likely to be critical for mitigating some of the most dangerous effects of climate change. In the last decade, there has been a significant increase in activity associated with reservoir characterization and site selection for large-scale CO2 storage projects across the globe. These prospective storage sites tend to be selected for their optimal structural, petrophysical, and geochemical trapping potential. However, it has also been suggested that storing CO2 in reservoirs within the CO2 hydrate stability zone (GHSZ), characterized by high pressures and low temperatures (e.g., Arctic or marine environments), could provide a natural thermodynamic barrier to gas leakage. Evaluating the prospect of commercial-scale, long-term storage of CO2 in the GHSZ requires reservoir-scale modeling capabilities designed to account for the unique physics and thermodynamics associated with these systems. We have developed the HYDRATE flow mode and the accompanying fully implicit parallel well model in the massively parallel subsurface flow and reactive transport simulator PFLOTRAN to model CO2 injection into the marine GHSZ. We have applied these capabilities to a set of CO2 injection scenarios designed to reveal the challenges and opportunities for commercial-scale CO2 storage in the GHSZ.
Three-phase numerical model for subsurface hydrology in permafrost-affected regions (PFLOTRAN-ICE v1.0)
Degradation of near-surface permafrost due to changes in the climate is expected to impact the hydrological, ecological and biogeochemical responses of the Arctic tundra. From a hydrological perspective, it is important to understand the movement of the various phases of water (gas, liquid and ice) during the freezing and thawing of near-surface soils. We present a new non-isothermal, single-component (water), three-phase formulation that treats air as an inactive component. This single component model works well and produces similar results to a more complete and computationally demanding two-component (air, water) formulation, and is able to reproduce results of previously published laboratory experiments. A proof-of-concept implementation in the massively parallel subsurface flow and reactive transport code PFLOTRAN is summarized, and parallel performance of that implementation is demonstrated. When water vapor diffusion is considered, a large effect on soil moisture dynamics is seen, which is due to dependence of thermal conductivity on ice content. A large three-dimensional simulation (with around 6 million degrees of freedom) of seasonal freezing and thawing is also presented.
Photovoltaic power conditioning subsystem: State of the art and development opportunities
Photovoltaic systems, the state of the art of power conditioning subsystem components, and the design and operational interaction between photovoltaic systems and host utilities are detailed in this document. Major technical issues relating to the design and development of power conditioning systems for photovoltaic application are considered; these include: (1) standards, guidelines, and specifications; (2) cost effective hardware design; (3) impact of advanced components on power conditioning development; (4) protection and safety; (5) quality of power; (6) system efficiency; and (7) system integration with the host utility. Theories of harmonic distortion and reactive power flow are discussed, and information about power conditioner hardware and manufacturers is provided.
Comparison of numerical oblique detonation solutions with an asymptotic benchmark
In order to have confidence in a numerical method, the verification of its reproduction of known benchmark analytic solutions for simple model problems is of great importance. Attention is presently given to a novel benchmarking procedure for numerical models of high speed, reactive 2D flows. The procedure is illustrated by comparing asymptotic and numerical solutions for oblique detonations in which an attached oblique shock is followed by an exothermic reaction with a thick reaction zone.
Validation of the RPLUS3D Code for Supersonic Inlet Applications Involving Three-Dimensional Shock Wave-Boundary Layer Interactions
A three-dimensional computational fluid dynamics code, RPLUS3D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for glancing shock wave-boundary layer interactions. Both laminar and turbulent flows were studied. A supersonic flow over a wedge mounted on a flat plate was numerically simulated. For the laminar case, the static pressure distribution, velocity vectors, and particle traces on the flat plate were obtained. For turbulent flow, both the Baldwin-Lomax and Chien two-equation turbulent models were used. The static pressure distributions, pitot pressure, and yaw angle profiles were computed. In addition, the velocity vectors and particle traces on the flat plate were also obtained from the computed solution. Overall, the computed results for both laminar and turbulent cases compared very well with the experimentally obtained data.