Search NASA⌕ Search

SEARCH · Search NASA

Results for “step 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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling Geologic Waste Repository Systems Below Residual Saturation

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Hydraulic Response to Thermal Stimulation Efforts at Raft River Based on Stepped Rate Injection Testing

The injection well stimulation project at the Raft River geothermal field tests the effect of long-term cold water injection and high pressure injection on well injectivity, improvements to which could reduce operating costs. The primary data for analysis and interpretation of the injection test are step-rate flow tests run before each new phase of the injection. These tests were analyzed using a combination of standard pump-test analytical solution methods and methods developed expressly for the observed conditions. The stepped rate injection tests, combined with long-term flow and pressure response data suggest that the well is located within a fractured formation of low transmissivity but high storativity. These calculated parameters appeared to increase with pressure during the first injection test and the higher values were reproduced during the second stepped rate test. Calculated transmissivity and storativity are on the order of 4E-5 m cm 2 and 1E-4 m Pa -1 , respectively. The apparent pressure dependence of fitted hydraulic parameters may reflect near-well fracture compliance that increased the effective radius of the wellbore during the first test. While the type curve fit analysis also suggests that the reservoir behaves as a uniformly fractured reservoir with a radial flow regime, the hydraulic parameters indicate that condition may exist only a very limited distance (<10 m) from the well. Longer-term pressure response suggests that flow in the system effectively reaches steady state in a period of less than a day, which may reflect pressure stabilization resulting from pressure-dependent permeability or a region of much higher permeability located with a few meters of the well. The transmissivity estimates obtained from this analysis, converted to approximate fracture density and aperture, provide useful constraints on the distance to which the thermal front may migrate from the well during the cold water injection phase of the stimulation project. We estimate that the cooling front will migrate less than a tenth of a kilometer over an approximately one-year injection period. Here, the effects of that cooling, however, may be substantial, because increases in permeability have maximum effect nearest the well.

cold water injection↗
DOE OSTI2024

Three-step LES-C models for flows at high Reynolds numbers

Here, we investigate the need for the second correction step in the recently proposed LES-C models for fluid flows at high Reynolds numbers. These models use a predictor-corrector idea to enhance the efficiency of the existing Large Eddy Simulation models. Different three-step(one defect step, two corrections) LES-C models, based on the Leray-α, ADM and NS-ω LES models, are tested in three different situations. The new Leray-α-C2 model (C2 stands for two correction steps) is applied to the Navier–Stokes equations; the ADC2 is applied to the MagnetoHydroDynamic flow; and the NS-ω-C2 is used in the fluid-fluid interaction problem. We evaluate the effectiveness of the second correction step in all these settings, using qualitative and quantitative numerical tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

LEAKAGE CHARACTERISTICS OF A STEPPED BALANCE PISTON LABYRINTH SEAL IN SUPERCRITICAL CO2

In an integrally geared compressor, there is a net thrust load on the pinion shaft from gear loads and aerodynamic forces. The net thrust on the pinion is resolved with either a thrust bearing or a thrust collar (rider ring). As the pressure increases, this increases the aerodynamic force proportionally. In highpressure systems this can overwhelm the thrust bearing or thrust collar if thrust is mismanaged. In a recent test program for an sCO2 integrally geared compander (compressor/turbine), a balance piston was created on the opposite end of the turbine to provide an equalizing thrust; while this is not novel, the balance piston provided a test apparatus for measuring the leakage on a 20 tooth stepped labyrinth seal operating at pressures in excess of 250 bar (3600 psi). The stepped balance piston seal is a 20 tooth tooth-onrotor balance piston seal. This type of seal was selected since it was predicted to have substantially less leakage than a comparably sized straight through labyrinth seal. Comparisons between the two designs will be made, showing that the leakage of the stepped balance piston seal is considerably lower. The tooth clearance on the seal was selected to be 30 mils radial. While this seems large, a stack-up tolerance assessment showed that this magnitude of clearance was appropriate to prevent rotor stator contact. During the testing campaign, data was recorded showing the leakage across this stepped balance piston seal at pressures in excess of 250 bar (3600 psi) with a pressure ratio of 2. This leakage data was measured, and will be compared to leakage data predicted with CFD. Two measurement methods are presented for the stepped balance piston seal: an orifice plate and valve. The predicted seal leakage from a straight through seal was four times that of the predicted stepped labyrinth seal. The measured flow from the stepped seal was on the same order as the CFD prediction and offered adequate rotordynamic performance

wilkes, jason↗
DOE OSTI2020

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

The multi-region relaxed magnetohydrodynamics (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The stepped pressure equilibrium code (SPEC) was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad–Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Lastly, we provide an example equilibrium with a 3D helical field structure as the favored lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

3D‐Mapping and Manipulation of Photocurrent in an Optoelectronic Diamond Device

Abstract Establishing connections between material impurities and charge transport properties in emerging electronic and quantum materials, such as wide‐bandgap semiconductors, demands new diagnostic methods tailored to these unique systems. Many such materials host optically‐active defect centers which offer a powerful in situ characterization system, but one that typically relies on the weak spin‐electric field coupling to measure electronic phenomena. In this work, charge‐state sensitive optical microscopy is combined with photoelectric detection of an array of nitrogen‐vacancy (NV) centers to directly image the flow of charge carriers inside a diamond optoelectronic device, in 3D and with temporal resolution. Optical control is used to change the charge state of background impurities inside the diamond on‐demand, resulting in drastically different current flow such as filamentary channels nucleating from specific, defective regions of the device. Conducting channels that control carrier flow, key steps toward optically reconfigurable, wide‐bandgap optoelectronics are then engineered using light. This work might be extended to probe other wide‐bandgap semiconductors (SiC, GaN) relevant to present and emerging electronic and quantum technologies.

Wood, Alexander A.↗
DOE OSTI2023

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Unlocking the potential of biogas systems for energy production and climate solutions in rural communities

On-site conversion of organic waste into biogas to satisfy consumer energy demand has the potential to realize energy equality and mitigate climate change reliably. However, existing methods ignore either real-time full supply or methane escape when supply and demand are mismatched. Here, we show an improved design of community biogas production and distribution system to overcome these and achieve full co-benefits in developing economies. We take five existing systems as empirical examples. Mechanisms of synergistic adjusting out-of-step biogas flow rates on both the plant-side and user-side are defined to obtain consumption-to-production ratios of close to 1, such that biogas demand of rural inhabitants can be met. Furthermore, carbon mitigation and its viability under universal prevailing climates are illustrated. Coupled with manure management optimization, Chinese national deployment of the proposed system would contribute a 3.77% reduction towards meeting its global 1.5 °C target. Additionally, fulfilling others’ energy demands has considerable decarbonization potential.

09 BIOMASS FUELS↗
DOE OSTI2022

ALDWP Operational Highlights 2021

The Weapons Pit production manufacturing team completed pit build 21-1 assembly operations on December 23, 2020. This build utilized all essential flow sheet steps and included the second -01 long tube cut off extension and welding activity. This is the best build produced in the MC4597-01 pit program to date. Final out-of-line operations, non-destructive evaluation, and radiography was completed during the week of January 4. Depending on the results from high energy radiography, this pit will be considered to be a candidate for CERT-02 and shipment to Lawrence Livermore National Laboratory in support of pit certification. It is notable that this build was completed while concurrently making Plutonium (Pu) metal, producing additional castings, machining Pu parts for Build 22, completing disassembly of two pits, supporting Nightshade B Pu target characterization, and precision machining of Jasper certification base plates.

42 ENGINEERING↗
DOE OSTI2026

Arsenic Removal by Continuous Flow-Through Iron Electrocoagulation and Downstream Steps under Environmentally Relevant Conditions

Arsenic (As) contamination in drinking water is a major health concern. This study investigated As(III) and As(V) removal by continuous flow-through iron electrocoagulation (EC) under varying iron doses and water chemistry conditions. Effective arsenic removal (>90% of 100 µg/L) was achieved with a low iron dose of 5 mg/L and a short EC reactor residence time of 11 s that was followed by flocculation and settling. Removal was primarily through adsorption of As(III) and As(V) to the EC-generated Fe(III) oxyhydroxide solids ferrihydrite and lepidocrocite. For influents containing As(III), oxidation of As(III) to As(V) also occurred. The removal processes occurred in both the EC reactor and downstream treatment. EC-generated reactive solids with high surface area and adsorption capacity provided better arsenic removal than with preformed adsorbents. The removed arsenic remained largely immobile in the EC-generated solids. The extent of arsenic removal was generally similar across pH 6 to 8. Phosphate (1 mg/L as P) inhibited the removal of both arsenic species. Dissolved silica (20 mg/L as SiO2) suppressed As(III) removal but had little effect on the amount of As(V) removed. In addition to thermodynamic effects on adsorption equilibria, pH and co-occurring phosphate and silica influenced removal rates by affecting iron oxidation, nucleation, and aggregation kinetics. The results integrate mechanistic understanding with treatment performance, providing insight relevant to both fundamental research and the application of EC in water treatment.

Arsenic, electrocoagulation, iron oxide, phosphate↗
DOE OSTI2020

Toward particle-resolved accuracy in Euler–Lagrange simulations of multiphase flow using machine learning and pairwise interaction extended point-particle (PIEP) approximation

This study presents two different machine learning approaches for the modeling of hydrodynamic force on particles in a particle-laden multiphase flow. Results from particle-resolved direct numerical simulations (PR-DNS) of flow over a random array of stationary particles for eight combinations of particle Reynolds number (Re) and volume fraction (φ) are used in the development of the models. The first approach follows a two-step process. In the first flow prediction step, the perturbation flow due to a particle is obtained as an axisymmetric superposable wake using linear regression. In the second force prediction step, the force on a particle is evaluated in terms of the perturbation flow induced by all its neighbors using the generalized Faxén form of the force expression. In the second approach, the force data on all the particles from the PR-DNS simulations are used to develop an artificial neural network (ANN) model for direct prediction of force on a particle. Due to the unavoidable limitation on the number of fully resolved particles in the PR-DNS simulations, direct force prediction with the ANN model tends to over-fit the data and performs poorly in the prediction of test data. In contrast, due to the millions of grid points used in the PR-DNS simulations, accurate flow prediction is possible, which then allows accurate prediction of particle force. This hybridization of multiphase physics and machine learning is particularly important, since it blends the strength of each, and the resulting pairwise interaction extended point-particle model cannot be developed by either physics or machine learning alone.

42 ENGINEERING↗
DOE OSTI2025

Reaction Pathways over ZnZrO 2 -Based Catalysts and Catalytic Sorbents

Reactive capture and conversion (RCC) is a process intensification approach that integrates CO 2 capture and hydrogenation within a single unit, removing the CO 2 purification and storage steps of traditional process flow schemes. This alters the catalytic step from a traditional steady-state (SS) flow process to a transient capture and conversion cycle, which could lead to product distributions distinct from those observed in conventional SS experiments. Such differences are investigated in the combined capture and hydrogenation of carbon dioxide to methanol over a ZnZrO 2 catalyst and a ZnZrO 2 + NaNO 3 /Mg 3 AlO x catalytic sorbent (CS) using fixed-bed kinetic measurements, in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and steady-state isotopic transient kinetic analysis-DRIFTS (SSITKA-DRIFTS). Under SS conditions, ZnZrO 2 produced methanol through sequential hydrogenation of HCOO* and CH 3 O* intermediates. On the contrary, CO was attributed primarily to CO 2 dissociation at oxygen vacancies, as supported by isotopic shifts and measured reaction orders. For the CS, isotopic switching experiments suggested that monodentate carbonate species (CO 3 2− , abbreviated as m-CO 3 2− ) act as active intermediates that can be hydrogenated to HCOO* and subsequently to CH 3 O. Under RCC conditions, in situ DRIFTS and isotopic experiments reveal that m-CO 3 2− species formed during the CO 2 capture step follow two competing routes upon H 2 exposure: (i) direct hydrogenation to methane on the sorbent domain or (ii) migration of m-CO 3 2− to the ZnZrO 2 domain, where they are hydrogenated to methanol through the HCOO pathway. Overall, RCC enables carbonate hydrogenation routes not observed under SS cofeed conditions. Thus, the reaction pathways and rates during RCC can be different from operation under conventional SS conditions, and the product distribution is determined here by competition between carbonate hydrogenation on sorbent sites and migration to ZnZrO 2 for methanol synthesis.

CCUS↗
DOE OSTI2020

Efficient and Flexible Sensitivity Matrix Computation for Adaptive Electrical Capacitance Volume Tomography

Electrical capacitance tomography is a widely used sensor modality for flow imaging in many industrial settings. Adaptive Electrical Capacitance Volume Tomography (AECVT) extends the capabilities of traditional ECT by enabling direct volumetric imaging and an improved resolution. Construction of the sensitivity matrix is a necessary step to obtain flow images. This step requires computation of the electric field inside the sensing domain, which is done via a typical field solver such as the finite element method. In this work, we present an efficient and flexible method to construct the sensitivity matrix for Adaptive Electrical Capacitance Volume Tomography (AECVT) based on individual electrode segment excitations and their judicious combination to form desired matrix elements. We illustrate how the proposed method yields the same sensitivity matrix as the traditional method but at a much lower computational cost. Once all segment contributions are obtained, we also indicate how the proposed method, unlike the traditional approach, can generate the sensitivity matrix on demand for an arbitrary combination of synthetic electrodes and obviating the need for any additional field computations. Finally, we present image reconstruction results for two different experimental scenarios where the mutual capacitance data and the corresponding sensitivity vectors are obtained through the proposed measurement combination scheme.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

UAE6 - Wind Tunnel Tests Data - UAE6 - Sequence K - Raw Data

Sequence K: Step AOA, Probes (P) This sequence was designed to quantify the 3-D blade static angle-of-attack response in the presence of rotational influences by varying the blade pitch angle. Sequence K used an upwind, rigid turbine with a 0° cone angle. The wind speeds ranged from 6 m/s to 20 m/s, and data were collected at yaw angles of 0° and 30°. The rotor rotated at 72 RPM. Blade and probe pressure measurements were collected. The teeter dampers were replaced with rigid links, and these two channels were flagged as not applicable by setting the measured values in the data file to -99999.99 Nm. The teeter link load cell was pre-tensioned to 40,000 N. The blade pitch angle ramped continuously at 0.18°/s over a wide range of increasing and decreasing pitch angles. A step sequence was also performed. The blade pitch was stepped 5°, the flow was allowed to stabilize, and the pitch angle was held for 8 seconds. Then the pitch angle step was repeated. Again, a wide range of pitch angles was obtained, both increasing and decreasing. The file lengths for this sequence varied from 96 seconds to 6 minutes, depending on the pitch angle range. Some short points were collected at 0° yaw and 3° pitch to verify the functionality of the instrumentation. The file name convention used the initial letter K, followed by two digits specifying wind speed, followed by two digits for yaw angle, followed by RU, RD, or ST, followed by the repetition digit. The angle-of-attack motion was differentiated by RU (ramp up), RD (ramp down), and ST (step down then step up). This sequence is related to sequences L and R.

17 WIND ENERGY↗