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.

SourceResearch connectionFrom yearThrough year

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

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Enhanced Severe Transient Analysis for Prevention Technical Program Plan

This document outlines the development of a high fidelity, best estimate nuclear power plant severe transient simulation capability that will complement or enhance the integral system codes historically used for licensing and analysis of severe accidents. As with other tools in the Risk Informed Safety Margin Characterization (RISMC) Toolkit, the ultimate user of Enhanced Severe Transient Analysis and Prevention (ESTAP) capability is the plant decision-maker; the deliverable to that customer is a modern, simulation-based safety analysis capability, applicable to a much broader class of safety issues than is traditional Light Water Reactor (LWR) licensing analysis. Currently, the RISMC pathway’s major emphasis is placed on developing RELAP-7, a next-generation safety analysis code, and on showing how to use RELAP-7 to analyze margin from a modern point of view: that is, by characterizing margin in terms of the probabilistic spectra of the “loads” applied to systems, structures, and components (SSCs), and the “capacity” of those SSCs to resist those loads without failing. The first objective of the ESTAP task, and the focus of one task of this effort, is to augment RELAP-7 analyses with user-selected multi-dimensional, multi-phase models of specific plant components to simulate complex phenomena that may lead to, or exacerbate, severe transients and core damage. Such phenomena include: coolant crossflow between PWR assemblies during a severe reactivity transient, stratified single or two-phase coolant flow in primary coolant piping, inhomogeneous mixing of emergency coolant water or boric acid with hot primary coolant, and water hammer. These are well-documented phenomena associated with plant transients but that are generally not captured in system codes. They are, however, generally limited to specific components, structures, and operating conditions. The second ESTAP task is to similarly augment a severe (post-core damage) accident integral analyses code with high fidelity simulations that would allow investigation of multi-dimensional, multi-phase containment phenomena that are only treated approximately in established codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2001

Sedimentation along the Eastern Chenier Plain Coast: Down Drift Impact of a Delta Complex Shift

The Mississippi River Chenier Plain is a shore parallel landform (down-drift from the Atchafalaya distributary of the Mississippi River) consisting of an alternating series of transgressive sand-shell ridges and regressive, progradational mudflats. The late 1940s shift of 1/3 of the flow of the Mississippi to the newly developing Atchafalaya delta complex to the west has resulted in injection of the river waters and suspended sediment into the westward flowing currents of the coastal current system. This has reactivated the dormant processes of mud accumulation along this coast. These environmental circumstances have provided the opportunity to: (1) investigate the depositional processes of the prograding, fine grained, mud flat facies of the open Chenier main coast and (2) to test the hypothesis that the impacts of the frequent cold front passages of fall, winter and spring exceed those of the occasional and more localized hurricane in shaping the coast and powering the dominant sedimentary processes. We conducted field investigations with the benefit of multi - scale, time series environmental surveillance by remote sensing systems, including airborne and satellite sensors. These systems provided invaluable new information on areal geomorphic patterns and the behavior of the coastal waters. This is a classic case of weather impacting inner shelf waters and sediments and causing the development of a new landform. It is clear that mud flats of the eastern chenier plain are prograding seaward, as well as progressively growing in a westerly direction.

Huh, Oscar K.↗

NASA NTRSDate not supplied

Effect of Injection Dynamics on Wave Propagation in a Linear Detonation Combustor

Injector response to self-excited and sustained detonation wave propagation in an optically-accessible nonpremixed natural gas-oxygen linear detonation combustor is investigated. 100 kHz measurements of planar laser induced fluorescence (PLIF) of trace quantities of acetone injected into the fuel supply and simultaneous OH∗ chemiluminescence measurements were used to correlate the effect on injection dynamics on detonation wave structure and propagation. Measurements over a range of equivalence ratios (0.61-1.48) revealed wave propagation through regions of high fuel stratification in the axial and transverse direction of reactant injection. Consequently, the propagation of the detonation wave along the channel is through a flow-field with significant variation in fuel mass fraction, reactivity, and sound speeds. At lean conditions, wave propagation frequency is generally lower (∼6 kHz) resulting in longer time between wave passages for mixing and pre-heat resulting in steeper detonation wave-fronts anchored near the injector exit plane. Conversely, at rich equivalence ratios the wave frequency is higher (∼10 kHz) resulting in shorter intra-cycle period for reactant fill and mixing resulting waves lifted from the injection plane that are more compact. Higher chemiluminescence intensity in the post-wave region with deflagrative combustion is observed in these cases. Phase averaged representation of the acetone-PLIF measurements indicate that the fraction of time between wave passages available for reactant fill, mixing and pre-heat is consistent across all equivalence ratios.

Propulsion↗

NASA NTRSDate not supplied

Effect of Injection Dynamics on Detonation Wave Propagation in a Linear Detonation Combustor

Propulsion↗

DOE OSTI2020

A Combined Day-ahead and Real-time approach for Scheduling Real and Reactive power dispatch of Battery Energy Storage System

A combined day-ahead (look-ahead) and real-time optimization approach is proposed to schedule BESS real and reactive power dispatch. The look ahead scheduling agent maximizes profit and essentially performs the function of energy arbitrage by looking at forecasted prices, network constraints and BESS operating conditions to generate real power dispatch schedule. While this schedule is followed, the real time optimization agent computes an optimum reactive power set-point (for a particular objective function) using optimal power flow framework. The concepts are implemented by co-simulating legacy tools and demonstrated on a large real world distribution system. In addition to the profits from energy arbitrage, optimum reactive power set point computation results in achieving desirable network attributes like flat voltage profiles and minimizing voltage deviations from desired value.

Bhatti, Bilal A.↗

NASA NTRS1996

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Dugan, E. T.↗

DOE OSTI2022

Transmission-distribution long-term volt-var planning considering reactive power support capability of distributed PV

High penetration of grid-edge, inverter-based photovoltaic (PV) can cause significant voltage fluctuations not only at the distribution but also at the sub-transmission levels due to PV output intermittency. Traditional reactive power planning approaches do not consider intermittency, nor the possibility of coordinating the control of existing and future volt-ampere reactive resources. This paper proposes a reactive power planning tool for sub-transmission systems to mitigate voltage violations and fluctuations caused by high PV penetration and intermittency with a minimum investment cost. The planning tool coordinates with an optimization-based volt-var operational tool for: a) modeling the coordination of all existing var assets in both sub-transmission and distribution systems to reduce the need of new equipment, and b)selecting a set of scenarios with voltage violations, derived from PV intermittency c) testing the final investment decision. The tool obtains an investment need for each intermittency scenario with a proposed optimal power-flow framework with efficient techniques to handle a high number of discrete variables. Two options are provided for final planning decision: i) a conservative direct combination of investment need solutions and ii) a machine learning-based selection of representative investment needs at most time steps. The final investment decision options are verified using a realistic large-scale sub-transmission system and 5-minute PV and load data. The results show a significant voltage performance improvement with a lower investment cost for additional var equipment compared to conventional approaches.

14 SOLAR ENERGY↗

DOE OSTI2018

Microbial Sulfate Reduction and Perchlorate Inhibition in a Novel Mesoscale Tank Experiment

Typically, microbial sulfate reduction occurs ubiquitously in natural environments. In oil and gas reservoirs, the generation of sulfide (also known as souring) can result in the corrosion of steel infrastructure and downgrading of oil quality, among other environmental and health-related concerns. The complex interplay between hydrological, geochemical, and biological processes during souring is poorly understood, preventing effective treatment and mitigation especially in naturally heterogeneous subsurfaces. In this work, three-dimensional flow tank experiments are utilized as a mesoscale experiment that links well-constrained batch and column experiments to field measurements. The mesoscale tank experiment investigating perchlorate treatment of souring is coupled with reactive transport modeling to understand the effects of heterogeneity on souring and effectiveness of perchlorate treatment. Tracer experiments were performed at the start and end of the experiment to constrain flow pathways and heterogeneities. Isotopic, geochemical, and microbial data revealed that perchlorate effectively inhibited sulfidogenesis and the growth of dominant sulfate reducing Desulfobacteraceae . Perchlorate treatment enriched Desulfobulbaceae , a sulfur-oxidizing group of bacteria, and Sulfurimonas , a potential perchlorate reducer. More organisms, including sulfate reducing bacteria, were observed closer to the influent. Results from the three-dimensional reactive transport model indicate horizontal preferential flows, as a result of the permeability contrast, led to faster bacteria growth (sulfate reducing bacteria) and sulfate reduction in fast flow regions. Our work highlights the control that spatial distributions of hydrologic characteristics exert over reservoir souring and treatment.

58 GEOSCIENCES↗

DOE OSTI2017

Effects of lignin structure on hydrodeoxygenation reactivity of pine wood lignin to valuable chemicals

Hydrodeoxygenation (HDO) of two dilute acid flow through pretreated softwood lignin samples, including residual lignin in pretreated solid residues (ReL) and recovered insoluble lignin in pretreated liquid (RISL), with apparent different physical and chemical structures, was comprehensively studied. A combination of catalysts (HY zeolite and Ru/Al 2 O 3 ) was employed to investigate the effects of lignin structures, especially condensed structures, on the HDO upgrading process. Results indicated that the condensed structure and short side chains in lignin hindered its HDO conversion under different reaction conditions, including catalyst loading and composition, hydrogen pressure, and reaction time. In addition to lignin structure, HY zeolite was found crucial for lignin depolymerization, while Ru/Al 2 O 3 and relatively high hydrogen pressure (4 MPa) were necessary for upgrading unstable oxy-compounds to cyclohexanes at high selectivity (>95 wt %). Since the lignin structure essentially affects its reactivity during HDO conversion, the yield and selectivity of HDO products can be predicted by detailed characterization of the lignin structure. Furthermore, the insights gained from this study in the fundamental reaction mechanisms based on the lignin structure will facilitate upgrading of lignin to high-value products for applications in the production of both fuels and chemicals.

09 BIOMASS FUELS↗

NASA NTRS1999

Hydrodynamics of Spherical Flows and Geometry of Premixed Flames near the Stagnation Point of Axisymmetric Viscous Counterflows

Counterflow premixed flames play a significant role in the modeling of laminar flames. This is in part motivated by the fact that stretched premixed flames simulate local flamelet dynamics within turbulent premixed flames. In the present study, the modified form of the Navier-Stokes equation for reactive fields introduced earlier is employed to investigate the hydrodynamics of spherical flows embedded within counterflows. The geometry of premixed flames near the stagnation point is also determined. The predictions are in favorable agreement with the experimental observations and prior numerical studies.

Sohrab, Siavash H.↗

DOE OSTI2015

Catalyst synthesis and evaluation using an integrated atomic layer deposition synthesis–catalysis testing tool

An integrated atomic layer deposition-catalysis (I-ALD-CAT) tool was developed, combining an ALD manifold with a plug-flow reactor system for the synthesis of supported catalytic materials by ALD and immediate evaluation of catalyst reactivity using gas-phase probe reactions. The I-ALD-CAT system can deliver gaseous reagents comprised of 12 different metal ALD precursors, 4 oxidizing or reducing agents, and 4 catalytic reaction feeds to either of the two plug-flow reactors. The system can employ reactor pressures and temperatures in the range of 10 -3 to 1 bar and 300–1000 K, respectively. Furthermore, the instrument is also equipped with a gas chromatograph and a mass spectrometer unit for the detection and quantification of volatile species from ALD and catalytic reactions. In this report, we demonstrate the use of the I-ALD-CAT tool for the ALD of platinum active sites and Al 2 O 3 overcoats, and evaluation of catalyst propylene hydrogenation activity.

47 OTHER INSTRUMENTATION↗

DOE OSTI2016

Investigation of Natural Circulation Instability and Transients in Passively Safe Small Modular Reactors

The NEUP funded project, NEUP-3496, aims to experimentally investigate two-phase natural circulation flow instability that could occur in Small Modular Reactors (SMRs), especially for natural circulation SMRs. The objective has been achieved by systematically performing tests to study the general natural circulation instability characteristics and the natural circulation behavior under start-up or design basis accident conditions. Experimental data sets highlighting the effect of void reactivity feedback as well as the effect of power ramp-up rate and system pressure have been used to develop a comprehensive stability map. The safety analysis code, RELAP5, has been used to evaluate experimental results and models. Improvements to the constitutive relations for flashing have been made in order to develop a reliable analysis tool. This research has been focusing on two generic SMR designs, i.e. a small modular Simplified Boiling Water Reactor (SBWR) like design and a small integral Pressurized Water Reactor (PWR) like design. A BWR-type natural circulation test facility was firstly built based on the three-level scaling analysis of the Purdue Novel Modular Reactor (NMR) with an electric output of 50 MWe, namely NMR-50, which represents a BWR-type SMR with a significantly reduced reactor pressure vessel (RPV) height. The experimental facility was installed with various equipment to measure thermalhydraulic parameters such as pressure, temperature, mass flow rate and void fraction. Characterization tests were performed before the startup transient tests and quasi-steady tests to determine the loop flow resistance. The control system and data acquisition system were programmed with LabVIEW to realize the realtime control and data storage. The thermal-hydraulic and nuclear coupled startup transients were performed to investigate the flow instabilities at low pressure and low power conditions for NMR-50. Two different power ramps were chosen to study the effect of startup power density on the flow instability. The experimental startup transient results showed the existence of three different flow instability mechanisms, i.e., flashing instability, condensation induced flow instability, and density wave oscillations. In addition, the void-reactivity feedback did not have significant effects on the flow instability during the startup transients for NMR-50. ii Several initial startup procedures with different power ramp rates were experimentally investigated to eliminate the flow instabilities observed from the startup transients. Particularly, the very slow startup transient and pressurized startup transient tests were performed and compared. It was found that the very slow startup transients by applying very small power density can eliminate the flashing oscillations in the single-phase natural circulation and stabilize the flow oscillations in the phase of net vapor generation. The initially pressurized startup procedure was tested to eliminate the flashing instability during the startup transients as well. The pressurized startup procedure included the initial pressurization, heat-up, and venting process. The startup transient tests showed that the pressurized startup procedure could eliminate the flow instability during the transition from single-phase flow to two-phase flow at low pressure conditions. The experimental results indicated that both startup procedures were applicable to the initial startup of NMR. However, the pressurized startup procedures might be preferred due to short operating hours required. In order to have a deeper understanding of natural circulation flow instability, the quasi-steady tests were performed using the test facility installed with preheater and subcooler. The effect of system pressure, core inlet subcooling, core power density, inlet flow resistance coefficient, and void reactivity feedback were investigated in the quasi-steady state tests. The experimental stability boundaries were determined between unstable and stable flow conditions in the dimensionless stability plane of inlet subcooling number and Zuber number. To predict the stability boundary theoretically, linear stability analysis in the frequency domain was performed at four sections of the natural circulation test loop. The flashing phenomena in the chimney section was considered as an axially uniform heat source. And the dimensionless characteristic equation of the pressure drop perturbation was obtained by considering the void fraction effect and outlet flow resistance in the core section. The theoretical flashing boundary showed some discrepancies with previous experimental data from the quasi-steady state tests. In the future, thermal non-equilibrium was recommended to improve the accuracy of flashing instability boundary. As another part of the funded research, flow instabilities of a PWR-type SMR under low pressure and low power conditions were investigated experimentally as well. The NuScale reactor design was selected as the prototype for the PWR-type SMR. In order to experimentally study the natural circulation behavior of NuScale iii reactor during accidental scenarios, detailed scaling analyses are necessary to ensure that the scaled phenomena could be obtained in a laboratory test facility. The three-level scaling method is used as well to obtain the scaling ratios derived from various non-dimensional numbers. The design of the ideally scaled facility (ISF) was initially accomplished based on these scaling ratios. Then the engineering scaled facility (ESF) was designed and constructed based on the ISF by considering engineering limitations including laboratory space, pipe size, and pipe connections etc. PWR-type SMR experiments were performed in this well-scaled test facility to investigate the potential thermal hydraulic flow instability during the blowdown events, which might occur during the loss of coolant accident (LOCA) and loss of heat sink accident (LOHS) of the prototype PWR-type SMR. Two kinds of experiments, normal blowdown event and cold blowdown event, were experimentally investigated and compared with code predictions. The normal blowdown event was experimentally simulated since an initial condition where the pressure was lower than the designed pressure of the experiment facility, while the code prediction of blowdown started from the normal operation condition. Important thermal hydraulic parameters including reactor pressure vessel (RPV) pressure, containment pressure, local void fraction and temperature, pressure drop and natural circulation flow rate were measured and analyzed during the blowdown event. The pressure and water level transients are similar to the experimental results published by NuScale [51], which proves the capability of current loop in simulating the thermal hydraulic transient of real PWR-type SMR. During the 20000s blowdown experiment, water level in the core was always above the active fuel assemble during the experiment and proved the safety of natural circulation cooling and water recycling design of PWR-type SMR. Besides, pressure, temperature, and water level transient can be accurately predicted by RELAP5 code. However, the oscillations of natural circulation flow rate, water level and pressure drops were observed during the blowdown transients. This kind of flow oscillations are related to the water level and the location upper plenum, which is a path for coolant flow from chimney to steam generator and down comer. In order to investigate the transients start from the opening of ADS valve in both experimental and numerical way, the cold blow-down experiment is conducted. For the cold blowdown event, different from setting both reactor iv pressure vessel (RPV) and containment at high temperature and pressure, only RPV was heated close to the highest designed pressure and then open the ADS valve, same process was predicted using RELAP5 code. By doing cold blowdown experiment, the entire transients from the opening of ADS can be investigated by code and benchmarked with experimental data. Similar flow instability observed in the cold blowdown experiment. The comparison between code prediction and experiment data showed that the RELAP5 code can successfully predict the pressure void fraction and temperature transient during the cold blowdown event with limited error, but numerical instability exists in predicting natural circulation flow rate. Besides, the code is lack of capability in predicting the water level related flow instability observed in experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2017

Methodology of Experiment Design for Improving Transient Critical Heat Flux Model Development

One phenomenon that contributes to nuclear fuel and clad failure mechanisms is the departure from nucleate boiling (DNB). Accident scenarios such as reactivity-initiated accidents (RIAs) provide conditions where the frequency of boiling crisis is much more likely to occur and on much shorter time scales than loss of coolant (LOCA) or loss of flow (LOFA) accidents. Here, boiling crisis in the form of DNB is attributed to transient-induced heat flux values that exceed the critical heat flux (CHF). The resulting effects, such as bubble crowding and eventual dryout, can affect flow boiling heat transfer and compromise the structural integrity of nuclear reactor materials. The focus of the work proposed here lies in developing a rigorous process using out-of-pile testing, under simulated RIA conditions, to determine the transient critical heat flux of cladding materials, without contributions from aforementioned radiation-induced effects that have been observed during in-pile experiments. This study provides a road-map for developing and performing experiments that simulate the power profile and energy deposition aspects of transient events, specific to RIAs, on the laboratory scale. The goal of our project is to observe and measure DNB in an annular flow channel around a single heater rod that is operating and nominal thermal hydraulic conditions of a PWR. It is also an aim of these experiments to determine how to quantify the gap between SS and transient CHF cases, based on state-dependent properties of the cladding material type that is being interrogated. Traditionally, CHF is determined as a combination of four factors: hydraulic diameter, system pressure, mass flow, fluid quality. The authors propose that an additional measureable quantity, which is dependent on the structural properties of the cladding material itself, is partly reasonable for the difference between measured and model-developed CHF values. Using a carefully metered approach will allow for several key concepts to be realized: compensation for transient-based effects to current models for CHF; adaptation of current models for different cladding material types, including Gen-IV concept materials; and implementation of a correction factor to account for both power transient-based and radiation-induced effects on the CHF that have not yet been accounted for.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗

DOE OSTI2025

Permeate fluxes from desalination of brines and produced waters: A reactive transport modeling study

The increasing interest in the use of membrane systems to desalinate inland brackish water, agricultural drainage, and industrially produced wastewater demands improved means of predicting desalination system performance under variable feedwater compositions. The interaction among water flow, solute transport, and chemical composition in these systems impacts permeate flux evolution. Here, an established multicomponent reactive transport simulator that accounts for these coupled processes is applied to compute osmotic pressure and permeate fluxes in reverse osmosis (RO) systems. The model is first validated by predicting permeate fluxes for a set of benchtop crossflow experiments subject to a range of feed flow rates and compositions, under fouling and non-fouling conditions. Results compare favorably with measured data that show that solutions with similar total dissolved solids concentrations but different compositions result in different permeate fluxes. The model is then applied to predict permeate fluxes from the desalination of produced waters using a commercial spiral wound RO module. For NaCl-dominant brines, at total dissolved salt concentrations (TDS) below about 70 g/L, permeate fluxes are inversely proportional to water mole fraction as the latter is a reasonable approximation of water activity (i.e. ideal mixing). In the case of Ca–Cl-, Na–CO3- and Na–SO4-dominant brines below about 70 g/L TDS, this relationship does not hold as well and tends to overpredict osmotic pressure and thus underpredict permeate fluxes. However, the opposite becomes true at higher TDS values for typical produced waters. The scaling potential of these waters is also computed by allowing the precipitation of minerals above their saturation limit on the RO membrane. This work demonstrates how reactive transport models developed for the analysis of waters from geological systems can be extended to improve process design, optimization, and control in desalination systems from produced waters and beyond.

Molins, Sergi↗

DOE OSTI2026

Fault Slip and Fluid Flow: Seismic Source Analysis to Assess Role of Multiple Slip Patches in Fault Permeability

The relationship between fault reactivation, microearthquakes (MEQs), and permeability evolution during fluid injection plays a critical role in energy harvesting and waste disposal. Recent studies have demonstrated the possibility of predicting fault permeability using cumulative seismic moments of MEQs quantitatively. To understand the underlying physical processes, we conduct fault reactivation experiments using Utah FORGE granitoid and analyze acoustic emission (AE) signals generated during stepwise increases in fluid injection pressure. Frequency analysis of thousands of calibrated AE signals reveals that fault reactivation produces multiple AE source patches with millimeter-scale radii—smaller than the sample fault radius. The cumulative area of the reactivated patches covers the fault multiple times over (∼10x–50x area) for each pressure step. These findings provide mechanistic insight that measured permeability enhancement is not driven by a single large slip event, but by the sequential and interacting activation of multiple slip patches that create a continuous flow pathway.

Nurshal, M. E. M. [Pennsylvania State University, ↗

DOE OSTI2017

Linear permeability evolution of expanding conduits due to feedback between flow and fast phase change: Linear Permeability Evolution of Conduit

Conduits are ubiquitous and critical pathways for many fluids relevant for geophysical processes such as magma, water, and gases. Predicting flow through conduits is challenging when the conduit geometry coevolves with the flow. We theoretically show that the permeability (k) of a conduit whose walls are eroding due to fast phase change increases linearly with time because of a self-reinforcing mechanism. This simple result is surprising given complex feedbacks between flow, transport, and phase change. The theory is congruent with previous experimental observations of fracture dissolution in calcite. Supporting computational fracture dissolution experiments showed that k only slightly increases until the dissolution front reaches the narrowest conduit constriction, after which the linear evolution of k manifests. Finally, the theory holds across multiple scales and a broad range of Peclet and Damkohler numbers and thus advances the prediction of dynamic mass fluxes through expanding conduits in various geologic and environmental settings

58 GEOSCIENCES↗