Search NASA⌕ Search

SEARCH · Search NASA

Results for “complex fluid”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Flow Boiling and Condensation Experiment (FBCE): Summary of Findings for Flow Boiling based on Completed ISS Experiments

Since 2011, researchers from Purdue University and NASA Glenn Research Center (GRC) have been collaborating to investigate the effects of gravity on several aspects of flow boiling and flow condensation. This massive research endeavor, termed the Flow Boiling and Condensation Experiment (FBCE), has culminated in development of NASA’s largest and most complex facility for investigation of two-phase fluid physics onboard the ISS. FBCE consists of two separate studies: flow boiling, using the Flow Boiling Module (FBM), and flow condensation, using the Condensation Module for Heat Transfer Measurements (CM-HT); this presentation concerns the FBM portion of FBCE. Fitted with FBM, FBCE was launched to the ISS in August 2021 and experiments were successfully performed starting in February 2022 to amass a large microgravity-flow-boiling database supported by extensive visualization of interfacial behavior using high-speed video. Based on FBCE experiments performed previously in Earth gravity and parabolic flight, along with the new ISS data, several empirical correlations, theoretical models, and CFD models have been developed, enabling the prediction of crucial flow boiling parameters, including heat transfer coefficient and critical heat flux (CHF). These predictive tools, which will be summarized in this presentation, are expected to serve as foundation for design of future space systems involving flow boiling.

Two-phase flow and phase transition↗
NASA NTRSDate not supplied

Flow Boiling and Condensation Experiment (FBCE): Latest Findings from the Summary ISS Experiments

Since 2011, researchers from Purdue University and NASA Glenn Research Center (GRC) have been collaborating to investigate the effects of gravity on several aspects of flow boiling and flow condensation. This massive research endeavor, termed the Flow Boiling and Condensation Experiment (FBCE), has culminated in development of NASA’s largest and most complex facility for investigation of two-phase fluid physics onboard the ISS. FBCE consists of two separate studies: flow boiling, using the Flow Boiling Module (FBM), and flow condensation, using the Condensation Module for Heat Transfer Measurements (CM-HT); the FBM experiments have already been completed while the condensation experiments began in 2024. This presentation will summarize mostly new results from the flow boiling experiments, with a focus on analysis of pressure drop and two-phase flow instabilities in microgravity using both experimental data and video records from the ISS experiments, as well as development of machine learning models. These new predictive tools are part of the arsenal of predictive methods developed by the Purdue-Glenn team for design of future space systems.

Microgravity↗
DOE OSTI2022

Simulation study of chemo-mechanical impacts of CO 2 injection in morrow b sandstone reservoir

Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.

58 GEOSCIENCES↗
DOE OSTI2022

The Geology of The Mt. Simon Sandstone Storage Complex at the Wabash #1 Well, Vigo Co., Indiana (Subtask 7.2, Technical Report)

The Wabash CarbonSAFE project drilled the Wabash #1 stratigraphic test well (ID# 168045) at the Wabash Valley Resources (WVR) IGCC facility in Vigo County, Indiana, to characterize and evaluate the basal Cambrian Mt. Simon Sandstone for commercial-scale CO 2 storage near the site. This report presents an extensive geologic characterization of the Mt. Simon storage complex and relevant data collected from the Wabash #1 well, such as lithologic data collected from cuttings and core, geophysical logging, geomechanical analysis of core samples, and well testing and fluid sampling within the Mt. Simon Sandstone. The Mt. Simon storage complex comprises two major sections: the Mt. Simon Sandstone as the potential reservoir and the overlying Eau Claire Formation as its primary seal. Within the report, an extensive depositional, sedimentological, and geochronologic characterization of the Mt. Simon is included with supportive chapters on the regional geology and the geophysical, petrophysical, and petrologic data collected during the project. An overview of 2D seismic reflection data collected from and around the test well is presented. Also presented are chapters on the characterization of the sealing Eau Claire Formation, including a chapter on the capacity of the primary and secondary seals to the Mt. Simon as well as a chapter on geomechanical testing results of the Eau Claire Formation and Mt. Simon Sandstone. Some of the information discussed in this report was used in the development of static and dynamic geologic models of the Mt. Simon Sandstone storage complex. The static and dynamic modeling of CO 2 injection in the Mt. Simon Sandstone are discussed in a separate report (Dessenberger et al., 2022) under the Wabash CarbonSAFE project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗
DOE OSTI2022

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗
DOE OSTI2019

Simulation of discrete cracks driven by nearly incompressible fluid via 2D combined finite-discrete element method

Here, we propose a novel hydraulic solver in order to simulate key mechanisms that control fluid-driven cracks in the framework of the combined finite-discrete element method (FDEM). The main innovative aspect of the present work is the independence of the fluid's critical time step size with the fracture opening. This advantage is extremely important because it means that very fine meshes can be used around areas of interest, such as boreholes, without penalizing the computational cost as fractures propagate (ie, open) and the fluid flows through them. This is a great advantage over other recently introduced approaches that exhibit a dependency of the time step in the form of Δt crit ∝ (l/a) 2 where l is the element size and a is the fracture aperture. This paper presents a series of benchmark cases for the proposed solver. The rationale adopted by the authors was to benchmark and validate the implementation of the hydraulic solver in an incremental fashion, starting from the simplest cases and building in complexity. The calculations shown in this work clearly demonstrate that the proposed approach is able to reproduce analytical results for fluid flow through a single crack. The results presented in this paper also demonstrate that the new approach is robust enough to deal with complex fracture patterns and complex geometries; the obtained fluid-driven fracture patterns in the vicinity of a borehole certainly stand to the scrutiny of human visual perception.

58 GEOSCIENCES↗
NASA NTRS1998

Overview of Sensitivity Analysis and Shape Optimization for Complex Aerodynamic Configurations

This paper presents a brief overview of some of the more recent advances in steady aerodynamic shape-design sensitivity analysis and optimization, based on advanced computational fluid dynamics. The focus here is on those methods particularly well- suited to the study of geometrically complex configurations and their potentially complex associated flow physics. When nonlinear state equations are considered in the optimization process, difficulties are found in the application of sensitivity analysis. Some techniques for circumventing such difficulties are currently being explored and are included here. Attention is directed to methods that utilize automatic differentiation to obtain aerodynamic sensitivity derivatives for both complex configurations and complex flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid computational fluid dynamics algorithms are demonstrated for different formulations of the sensitivity equations. Finally, the use of advanced, unstructured-grid computational fluid dynamics in multidisciplinary analyses and multidisciplinary sensitivity analyses within future optimization processes is recommended and encouraged.

Newman, Perry A.↗
DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗
NASA NTRS2015

Unsteady Three-Dimensional Simulation of a Shear Coaxial GO2/GH2 Rocket Injector with RANS and Hybrid-RAN-LES/DES Using Flamelet Models

Historically, the analysis and design of liquid rocket engines (LREs) has relied on full‐scale testing and one‐dimensional empirical tools. The testing is extremely expensive and the one‐dimensional tools are not designed to capture the highly complex, and multi‐dimensional features that are inherent to LREs. Recent advances in computational fluid dynamics (CFD) tools have made it possible to predict liquid rocket engine performance, stability, to assess the effect of complex flow features, and to evaluate injector‐driven thermal environments, to mitigate the cost of testing. Extensive efforts to verify and validate these CFD tools have been conducted, to provide confidence for using them during the design cycle. Previous validation efforts have documented comparisons of predicted heat flux thermal environments with test data for a single element gaseous oxygen (GO2) and gaseous hydrogen (GH2) injector. The most notable validation effort was a comprehensive validation effort conducted by Tucker et al. [1], in which a number of different groups modeled a GO2/GH2 single element configuration by Pal et al [2]. The tools used for this validation comparison employed a range of algorithms, from both steady and unsteady Reynolds Averaged Navier‐Stokes (U/RANS) calculations, large‐eddy simulations (LES), detached eddy simulations (DES), and various combinations. A more recent effort by Thakur et al. [3] focused on using a state‐of‐the‐art CFD simulation tool, Loci/STREAM, on a two‐dimensional grid. Loci/STREAM was chosen because it has a unique, very efficient flamelet parameterization of combustion reactions that are too computationally expensive to simulate with conventional finite‐rate chemistry calculations. The current effort focuses on further advancement of validation efforts, again using the Loci/STREAM tool with the flamelet parameterization, but this time with a three‐dimensional grid. Comparisons to the Pal et al. heat flux data will be made for both RANS and Hybrid RANSLES/ Detached Eddy simulations (DES). Computation costs will be reported, along with comparison of accuracy and cost to much less expensive two‐dimensional RANS simulations of the same geometry.

Westra, Doug G.↗
DOE OSTI2023

Intermediate Temperature Proton‐Conducting Solid Oxide Electrolysis Cells with Improved Performance and Durability

Technical understanding: A high-performance and stable intermediate temperature proton-conducting solid oxide electrolysis cell (PCEC) will add significantly to the hydrogen roadmap of EERE. However, the low electrocatalytic activity of the air electrode and the instability of the electrolyte are technical barriers. In this research, we utilized computational fluid dynamic modeling and materials/microstructure engineering to investigate the complex reaction chain and to improve the performance and reliability of the protonic ceramic electrochemical cells. An active and stable Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ (BZCYYb) scaffold conformal coating electrode was rationally designed where the robust Pr 1.8 Ba 0.2 NiO 4+δ functions both as a catalytic active and a protective layer. The PCECs with porosity-optimized scaffold conformal coating electrodes demonstrate 1.46 Acm -2 under 1.3 V at 600°C, which is among the highest in recent years. The final optimization by computational fluid dynamic modeling and conformal coating technique was applied to evolve the PCEC air electrode from Pr 2 NiO 4+δ to Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ scaffold conformal coating. The electrolysis current was boosted from 0.445A/cm 2 to 1.46A/cm 2 at 600°C and 1.3V. A superior low degradation rate of 1mV/kh was achieved in a 5000-hour of operation. Technical effectiveness and economic feasibility: The effectiveness of this TCO conformal coating on electrolyte backbone as air electrode has been proven by long-term demonstration. Based on the highly conductive BZCYYb electrolyte and the TCO steam electrode PNO, this project rationally took advantage of the existent materials to structure a PCEC device with high TRL and low risk. The solution infiltration method used to structure the conformal coating is a standard method used by the solid oxide cell community. Very minimal adoption is needed in the current SOC fabrication process. It is a cost-effective, scalable method for manufacturing, and brings up remarkable impact on the current technical status of PCEC. Benefit to the public: The findings in the research showed the industry that the severe challenge in PCEC can be overcome in a practical way. A various energy-hydrogen ecosystem is not beyond reach. The success of this technology provides sustainable and affordable H 2 to society. It will stimulate investment from the governmental and private sectors in response to the device manufacture and deployment, and the related infrastructure construction. It will generate high-pay domestic jobs and environmental benefits to the general public in the United States.

08 HYDROGEN↗
DOE OSTI2020

Wind Farm Simulation and Layout Optimization in Complex Terrain: Preprint

This work reports on incorporating complex terrain into wind farm simulations for the purpose of layout optimization. Adding complex terrain boundary conditions to NREL's medium fidelity computational fluid dynamics model, WindSE, produces significant separation, flow curvature, and speedup effects that would otherwise be difficult to capture with lower-fidelity models or a flat-terrain assumption. These flow features, in turn, can significantly impact the optimal turbine array layout. We demonstrate the impact of complex terrain on flow in both an idealized and real-world setting, and discuss modifications to the code that enable gradient-based optimization using terrain-aware adjoint gradients. Through several optimization case studies, we show that the layout optimization process takes advantage of speedup effects on terrain high points, and leverages flow curvature effects that modify wake trajectories. This yields substantial power improvements over gridded layouts, and hints at future research directions in simulation and optimization for wake trajectories in complex terrain.

17 WIND ENERGY↗
NASA NTRS1991

Analysis and representation of complex structures in separated flows

We discuss our recent work on extraction and visualization of topological information in separated fluid flow data sets. As with scene analysis, an abstract representation of a large data set can greatly facilitate the understanding of complex, high-level structures. When studying flow topology, such a representation can be produced by locating and characterizing critical points in the velocity field and generating the associated stream surfaces. In 3D flows, the surface topology serves as the starting point. The 2D tangential velocity field near the surface of the body is examined for critical points. The tangential velocity field is integrated out along the principal directions of certain classes of critical points to produce curves depicting the topology of the flow near the body. The points and curves are linked to form a skeleton representing the 2D vector field topology. This skeleton provides a basis for analyzing the 3D structures associated with the flow separation. The points along the separation curves in the skeleton are used to start tangent curve integrations. Integration origins are successively refined to produce stream surfaces. The map of the global topology is completed by generating those stream surfaces associated with 3D critical points.

Helman, James↗
DOE OSTI2018

Lagrangian particle method for compressible fluid dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Final Scientific/Technical Report: Correlations and Fluctuations in Weakly Collisional Plasma

Plasma is a state of matter that exhibits a very rich range of phenomena. To begin with, plasma is both electrical and mechanical - bringing together theories of particle motion and the electromagnetic field. Furthermore, and especially important for this project, a weakly-collisional plasma, such as is found in high-temperature (fusion energy) experiments on earth and the majority of contexts in space and astrophysics, has many moving parts. For example, sitting in earth’s atmosphere we are immersed in a mechanical wave field (sound), a possibly turbulent fluid motion (wind), and an electromagnetic vector wave field with two polarizations (light). This is already enough to produce a rich range of possibilities. In plasma, the electromagnetic field is coupled to the mechanical motion of the medium because it is ionized. Furthermore, a weakly-collisional plasma supports an infinite number of mechanically independent fluids. Thus, plasmas support an infinite number of independent electromechanical waves. Much has been done to describe plasmas with "reduced models" of various kinds. The goal of this project was to both explore the validity of reduced plasma models that are in use, and to propose and validate new models of plasma motion. The primary means to his end was laboratory experiments employing both electrical probes and laser spectroscopy. Laser spectroscopy enables many techniques which can separate the spectrum of independent fluid motions in the ion phase-space. The choice was to focus on low frequency electrostatic waves because the electron motion is relatively simple, the experiments can be on a spatial scale of a few meters, and all the relevant parameters can be measured with a few lasers systems. No study of this kind had previously been undertaken for the study of plasmas. The validation of theories required that the experimental descriptions be compared with theory and simulation in detail. It was found that even multi-fluid theories leave out a large part of the complexity of plasma motion. Reduced descriptions were found to fail under most circumstances. A new technique was developed that enabled a measurement of the phase-space resolved ion correlation function for the first time. The wide range of plasma dynamics possible became clear through this technique. It was found that collisionless (Vlasov) theory has a large field of application even when the plasma is weakly-collisional. A new approach, the kinetic wave expansion, was proposed, tested and found to be very useful for describing electrostatic ion waves. This project demonstrated a new way of looking at the "degrees-of-freedom" of plasmas and provided significant validation tests of fluid and kinetic plasma descriptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Evidence for Carbonate Surface Complexation during Forsterite Carbonation in Wet Supercritical Carbon Dioxide

Continental flood basalts are attractive formations for geologic sequestration of carbon dioxide because of their reactive divalent-cation containing silicates, such as forsterite (Mg2SiO4), suitable for long-term trapping of CO2 mineralized as metal carbonates. The goal of this study was to investigate at a molecular level the carbonation products formed during the reaction of forsterite with supercritical CO2 (scCO2) as a function of the concentration of H2O adsorbed to the forsterite surface. Experiments were performed at 50 °C and 90 bar using an in situ IR titration capability, and post-reaction samples were examined by ex situ techniques, including SEM, XPS, FIB-TEM, TGA-MS, and MAS-NMR. Carbonation products and reaction extents varied greatly with adsorbed H2O. We show for the first time evidence of Mg-carbonate surface complexation under wet scCO2 conditions. Carbonate is found to be coordinated to Mg at the forsterite surface in a predominately bidentate fashion at adsorbed H2O concentrations below 27 µmol/m2. Above this concentration and up to 76 µmol/m2, monodentate coordinated complexes become dominant. Beyond a threshold adsorbed H2O concentration of 76 µmol/m2, crystalline carbonates continuously precipitate as magnesite, and the particles that form are hundreds of times larger than the estimated thicknesses of the adsorbed water films of about 7 to 15 Å. At an applied level, these results suggest that mineral carbonation in scCO2 dominated fluids near the wellbore and adjacent to caprocks will be insignificant and limited to surface complexation, unless adsorbed H2O concentrations are high enough to promote crystalline carbonate formation. At a fundamental level, the surface complexes and their dependence on adsorbed H2O concentration give insights regarding forsterite dissolution processes and magnesite nucleation and growth.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2015

Modeling the Interaction Between Hydraulic and Natural Fractures Using Dual-Lattice Discrete Element Method

The interaction between hydraulic fractures (HF) and natural fractures (NF) will lead to complex fracture networks due to the branching and merging of natural and hydraulic fractures in unconventional reservoirs. In this paper, a newly developed hydraulic fracturing simulator based on discrete element method is used to predict the generation of complex fracture network in the presence of pre-existing natural fractures. By coupling geomechanics and reservoir flow within a dual lattice system, this simulator can effectively capture the poro-elastic effects and fluid leakoff into the formation. When HFs are intercepting single or multiple NFs, complex mechanisms such as direct crossing, arresting, dilating and branching can be simulated. Based on the model, the effects of injected fluid rate and viscosity, the orientation and permeability of NFs and stress anisotropy on the HF-NF interaction process are investigated. Combined impacts from multiple parameters are also examined in the paper. The numerical results show that large values of stress anisotropy, intercepting angle, injection rate and viscosity will impede the opening of NFs.

02 PETROLEUM↗
DOE OSTI2022

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗