Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow velocity”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A study of methods to investigate nozzle boundary layer transition

To further investigate nozzle flow, numerical computations are employed. The computations produce complete flow velocity and temperature fields within the nozzle. As a check, these results can be compared with experimental data at the wall. Once an accurate numerical scheme has been validated, it can be used as a design tool to predict the performance of other nozzle designs without the cost of experimental testing. Typically, the numerical analysis assumes either a laminar boundary layer or a fully turbulent boundary layer which is steady and two-dimensional. Boundary layer transition is not considered. Computing both the completely laminar boundary layer and the completely turbulent boundary layer conditions gives the minimum and maximum wall heat flux possible for a specified geometry. When the experimental heat flux measurements lie between these two values, the nature of the boundary layer is unknown. The boundary layer may have transitioned from laminar to turbulent; three-dimensional structures may be present in the boundary layer, or the inlet flow conditions may not be correctly specified in the computation.

Pauley, Laura L.↗
DOE OSTI2021

Radiative Heat Transport and Optical Characterization of High Temperature Molten Salts

The goal of this project was to advance the field of molten fluoride radiative heat transfer (RHT), through optical property measurements, heat transfer loop construction and operation, and computational fluid dynamics simulations. The following outcomes have been accomplished. A comprehensive review has been conducted to compile all available RHT-relevant optical absorption data for molten fluorides. The various measurement methods have also been compiled and critiqued. Thermal emission measurement methods have been refined and documented. A unique thermal emissivity measurement facility has been constructed and used to measure the emissivity of flowing FLiBe-exposed SS316 samples. A slight increase has been observed in the emissivity compared to unexposed SS316. A dual reflectance-emission setup was partially demonstrated for the purpose of molten fluoride absorption measurement. Preliminary infrared edge data for FLiNaK was obtained. A visible near-infrared (NIR) measurement setup has been constructed and will be used to measure absorption of corrosion products in FLiBe. A new multi-band thermal radiation method using blackbody-weighted spectral bands for the absorption-coefficient was deployed to investigate the thermal performance of the UW Natural Circulation FLiBe Loop under various power levels and operating conditions. The results showed an increasing and decreasing trend in the overall heat transfer with optical thickness with a maximum at optical thickness τ D ≈Ο(1). A forced convection FLiBe loop was constructed for the purpose of studying corrosion and thermal hydraulic phenomena. A vertical cantilever pump is used to achieve up to 5 gpm and 5 m/s flow velocity in hot and cold test sections. A methodology has been established for modeling radiative heat transfer (RHT) in molten salt coolants for the purpose of quantifying the contribution of thermal radiation in conditions relevant to FHR or MSR systems. The comparison between the results obtained for forced and mixed convection showed that in the studied ranges of optical thickness, the maximum effect of RHT on overall heat transfer was 60% improvement observed at z/D = 250 for forced convection; whereas for mixed convection, that value was higher at 170% observed at z/D = 30. The highest effect of RHT on the overall energy distribution was observed at τ D ≈3.4. Negligible effect (<20%) was observed above τ D =18 and below τ D ≈0.002. A CFD code-to-code comparison was conducted using COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT, for a forced and mixed convection laminar flow for a grey medium in a vertical heated using both COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT. The relative error was the largest at τ D =40 for the radiative heat flux, bulk temperature, and Nusselt number. On the other hand, the relative error was the least at τ D =0.002 for the conductive heat flux, bulk temperature, and Nusselt number.

Anderson, Mark↗
NASA NTRS1999

NASA GRC's High Pressure Burner Rig Facility and Materials Test Capabilities

The High Pressure Burner Rig (HPBR) at NASA Glenn Research Center is a high-velocity. pressurized combustion test rig used for high-temperature environmental durability studies of advanced materials and components. The facility burns jet fuel and air in controlled ratios, simulating combustion gas chemistries and temperatures that are realistic to those in gas turbine engines. In addition, the test section is capable of simulating the pressures and gas velocities representative of today's aircraft. The HPBR provides a relatively inexpensive. yet sophisticated means for researchers to study the high-temperature oxidation of advanced materials. The facility has the unique capability of operating under both fuel-lean and fuel-rich gas mixtures. using a fume incinerator to eliminate any harmful byproduct emissions (CO, H2S) of rich-burn operation. Test samples are easily accessible for ongoing inspection and documentation of weight change, thickness, cracking, and other metrics. Temperature measurement is available in the form of both thermocouples and optical pyrometery. and the facility is equipped with quartz windows for observation and video taping. Operating conditions include: (1) 1.0 kg/sec (2.0 lbm/sec) combustion and secondary cooling airflow capability: (2) Equivalence ratios of 0.5- 1.0 (lean) to 1.5-2.0 (rich), with typically 10% H2O vapor pressure: (3) Gas temperatures ranging 700-1650 C (1300-3000 F): (4) Test pressures ranging 4-12 atmospheres: (5) Gas flow velocities ranging 10-30 m/s (50-100) ft/sec.: and (6) Cyclic and steady-state exposure capabilities. The facility has historically been used to test coupon-size materials. including metals and ceramics. However complex-shaped components have also been tested including cylinders, airfoils, and film-cooled end walls. The facility has also been used to develop thin-film temperature measurement sensors.

Robinson, R. Craig↗
NASA NTRSDate not supplied

Determining the Cause of Reduced Concurrent Flame Spread over Thin Solid Fuels in Low Pressure and Low Gravity

The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.

Fire Safety↗
NASA NTRSDate not supplied

Determining the Cause of Reduced Concurrent Flame Spread over Thin Solid Fuels in Low Pressure and Low Gravity

The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.

microgravity↗
DOE OSTI2025

Machine Learning Prediction of Tritium‐Helium Groundwater Ages in the Central Valley, California, USA

Abstract Groundwater ages provides insight into recharge rates, flow velocities, and vulnerability to contaminants. The ability to predict groundwater ages based on more accessible parameters via Machine Learning (ML) would advance our ability to guide sustainable management of groundwater resources. In this study, ML models were trained and tested on a large data set of tritium concentrations and tritium‐helium groundwater ages from the California Central Valley, a large groundwater basin with complex land use, irrigation, and water management practices. The ML models were trained on 63 features, including location, well construction information, landscape characteristics, and climate variables, water chemistry, and stable isotopes. The Bagging regressor method can accurately classify (F1‐score = 0.91) groundwater samples as either modern or pre‐modern whereas the accuracy of the ML prediction of continuous tritium‐helium groundwater ages is limited and explains only of the variability in this data set. In general, ML groundwater age prediction relies mostly on features related to (a) the source of groundwater recharge, (b) contaminant history, (c) aquifer materials, (d) well construction, and (e) geochemical reactions along flow paths.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

The cerebral hemodynamics of normotensive hypovolemia during lower-body negative pressure

Although severe hypovolemia can lead to hypotension and neurological decline, many patients with neurosurgical disorders experience a significant hypovolemia while autonomic compensatory mechanisms maintain a normal blood pressure. To assess the effects of normotensive hypovolemia upon cerebral hemodynamics, transcranial Doppler ultrasound monitoring of 13 healthy volunteers was performed during graded lower-body negative pressure of up to -50 mm Hg, an accepted laboratory model for reproducing the physiological effects of hypovolemia. Middle cerebral artery flow velocity declined by 16% +/- 4% (mean +/- standard error of the mean) and the ratio between transcranial Doppler ultrasound pulsatility and systemic pulsatility rose 22% +/- 8%, suggesting cerebral small-vessel vasoconstriction in response to the sympathetic activation unmasked by lower-body negative pressure. This vasoconstriction may interfere with the autoregulatory response to a sudden fall in blood pressure, and may explain the common observation of neurological deficit during hypovolemia even with a normal blood pressure.

Non-NASA Center↗
NASA NTRS1998

Autoregulation of cerebral blood flow in orthostatic hypotension

BACKGROUND AND PURPOSE: We sought to evaluate cerebral autoregulation in patients with orthostatic hypotension (OH). METHODS: We studied 21 patients (aged 52 to 78 years) with neurogenic OH during 80 degrees head-up tilt. Blood flow velocities (BFV) from the middle cerebral artery were continuously monitored with transcranial Doppler sonography, as were heart rate, blood pressure (BP), cardiac output, stroke volume, CO2, total peripheral resistance, and cerebrovascular resistance. RESULTS: All OH patients had lower BP (P<.0001), BFV_diastolic (P<.05), CVR (P<.007), and TPR (P<.02) during head-up tilt than control subjects. In control subjects, no correlations between BFV and BP were found during head-up tilt, suggesting normal autoregulation. OH patients could be separated into those with normal or expanded autoregulation (OH_NA; n=16) and those with autoregulatory failure (OH_AF; n=5). The OH_NA group showed either no correlation between BFV and BP (n=8) or had a positive BFV/BP correlation (R2>.75) but with a flat slope. An expansion of the "autoregulated" range was seen in some patients. The OH_AF group was characterized by a profound fall in BFV in response to a small reduction in BP (mean deltaBP <40 mm Hg; R2>.75). CONCLUSIONS: The most common patterns of cerebral response to OH are autoregulatory failure with a flat flow-pressure relationship or intact autoregulation with an expanded autoregulated range. The least common pattern is autoregulatory failure with a steep flow-pressure relationship. Patients with patterns 1 and 2 have an enhanced capacity to cope with OH, while those with pattern 3 have reduced capacity.

Non-NASA Center↗
DOE OSTI2024

Radiatively cooled magnetic reconnection experiments driven by pulsed power

We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly driven plasma flows (MA≈7) with anti-parallel magnetic fields, which form a reconnection layer (SL≈120) at the mid-plane. The net cooling rate far exceeds the Alfvénic transit rate (τcool−1/τA−1≫1), leading to strong cooling of the reconnection layer. We determine the advected magnetic field and flow velocity using inductive probes positioned in the inflow to the layer, and inflow ion density and temperature from analysis of visible emission spectroscopy. A sharp decrease in x-ray emission from the reconnection layer, measured using filtered diodes and time-gated x-ray imaging, provides evidence for strong cooling of the reconnection layer after its initial formation. X-ray images also show localized hotspots, regions of strong x-ray emission, with velocities comparable to the expected outflow velocity from the reconnection layer. These hotspots are consistent with plasmoids observed in 3D radiative resistive magnetohydrodynamic simulations of the experiment. X-ray spectroscopy further indicates that the hotspots have a temperature (170 eV) much higher than the bulk layer (≤75 eV) and inflow temperatures (about 2 eV) and that these hotspots generate the majority of the high-energy (>1 keV) emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Radiative cooling effects on reverse shocks formed by magnetized supersonic plasma flows

We study the structure of reverse shocks formed by the collision of supersonic, magnetized plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminum wire arrays, the shock is sharp and well-defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, and instead we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer that is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities, which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminum plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Cerebral versus systemic hemodynamics during graded orthostatic stress in humans

BACKGROUND: Orthostatic syncope is usually attributed to cerebral hypoperfusion secondary to systemic hemodynamic collapse. Recent research in patients with neurocardiogenic syncope has suggested that cerebral vasoconstriction may occur during orthostatic hypotension, compromising cerebral autoregulation and possibly contributing to the loss of consciousness. However, the regulation of cerebral blood flow (CBF) in such patients may be quite different from that of healthy individuals, particularly when assessed during the rapidly changing hemodynamic conditions associated with neurocardiogenic syncope. To be able to interpret the pathophysiological significance of these observations, a clear understanding of the normal responses of the cerebral circulation to orthostatic stress must be obtained, particularly in the context of the known changes in systemic and regional distributions of blood flow and vascular resistance during orthostasis. Therefore, the specific aim of this study was to examine the changes that occur in the cerebral circulation during graded reductions in central blood volume in the absence of systemic hypotension in healthy humans. We hypothesized that cerebral vasoconstriction would occur and CBF would decrease due to activation of the sympathetic nervous system. We further hypothesized, however, that the magnitude of this change would be small compared with changes in systemic or skeletal muscle vascular resistance in healthy subjects with intact autoregulation and would be unlikely to cause syncope without concomitant hypotension. METHODS AND RESULTS: To test this hypothesis, we studied 13 healthy men (age, 27 +/- 7 years) during progressive lower body negative pressure (LBNP). We measured systemic flow (Qc is cardiac output; C2H2 rebreathing), regional forearm flow (FBF; venous occlusion plethysmography), and blood pressure (BP; Finapres) and calculated systemic (SVR) and forearm (FVR) vascular resistances. Changes in brain blood flow were estimated from changes in the blood flow velocity in the middle cerebral artery (VMCA) using transcranial Doppler. Pulsatility (systolic minus diastolic/mean velocity) normalized for systemic arterial pressure pulsatility was used as an index of distal cerebral vascular resistance. End-tidal PACO2 was closely monitored during LBNP. From rest to maximal LBNP before the onset of symptoms or systemic hypotension, Qc and FBF decreased by 29.9% and 34.4%, respectively. VMCA decreased less, by 15.5% consistent with a smaller decrease in CBF. Similarly, SVR and FVR increased by 62.8% and 69.8%, respectively, whereas pulsatility increased by 17.2%, suggestive of a mild degree of small-vessel cerebral vasoconstriction. Seven of 13 subjects had presyncope during LBNP, all associated with a sudden drop in BP (29 +/- 9%). By comparison, hyperventilation alone caused greater changes in VMCA (42 +/- 2%) and pulsatility but never caused presyncope. In a separate group of 3 subjects, superimposition of hyperventilation during highlevel LBNP caused a further decrease in VMCA (31 +/- 7%) but no change in BP or level of consciousness. CONCLUSIONS: We conclude that cerebral vasoconstriction occurs in healthy humans during graded reductions in central blood volume caused by LBNP. However, the magnitude of this response is small compared with changes in SVR or FVR during LBNP or other stimuli known to induce cerebral vasoconstriction (hypocapnia). We speculate that this degree of cerebral vasoconstriction is not by itself sufficient to cause syncope during orthostatic stress. However, it may exacerbate the decrease in CBF associated with hypotension if hemodynamic instability develops.

NASA Discipline Cardiopulmonary↗
DOE OSTI2022

Experimental and computational investigation of the influence of iso -butanol on autoignition of n -decane and n -heptane in non-premixed flows

Experimental and computational investigation is carried out to elucidate the influence of iso-butanol on critical conditions of autoignition of n-decane and n-heptane. The counterflow configuration is employed. In this configuration an axisymmetric stream of air is directed over the surface of an evaporating pool of a liquid fuel. A stagnation plane is established. The temperature of the air is increased until autoignition is achieved in the mixing layer in the vicinity of the stagnation plane. The temperature of the air stream at autoignition, T ig , is measured at various values of the strain rate, which is defined as the axial gradient of the axial component of the flow velocity at the stagnation plane. Fuels tested are n-decane, n-heptane, iso-butanol and various mixtures of n-decane/iso-butanol and n-heptane/iso-butanol. Kinetic modeling is carried out using the recently updated version of the comprehensive CRECK chemical-kinetic mechanism. Critical conditions of autoignition are predicted for all fuels and fuel mixtures and the results are compared with the measurements. Computations show that low-temperature chemistry plays a significant role in promoting autoignition of n-decane and n-heptane, and the influence of low-temperature chemistry decreases with increasing strain rates because there is insufficient residence for the low-temperature reactions to take place. Experimental data and numerical simulations show that addition of even small amounts of iso-butanol to n-decane or n-heptane increases the value of T ig at low strain rates indicating that iso-butanol strongly inhibits the low-temperature chemistry of n-decane and n-heptane. Furthermore, the simulations show that when iso-butanol is present in the liquid fuel, only a small amount of n-decane is available in the gas phase and its low concentration cannot effectively sustain the low temperature degenerate branching oxidation process. Here, predicted flame structures show that the peak values of mole fraction of ketohydroperoxide are significantly reduced when iso-butanol is added to n-decane indicating that the kinetic pathway to low temperature ignition is blocked. This observation is confirmed from sensitivity analysis.

42 ENGINEERING↗
DOE OSTI2022

Ion temperature anisotropy model with cross‐field drifts in the scrape‐off layer

Abstract An ion temperature anisotropy model including cross‐field drifts has been implemented in the 2D transport code UEDGE. Two ion temperature equations are solved in the model to give separate ion temperatures perpendicular and parallel to the magnetic field, which results in an additional mirror force for ion parallel transport and also leads to an additional curvature drift in the perpendicular direction. In our previous work, an ion temperature anisotropy model was implemented and applied to the full 2D version of UEDGE including recycled neutrals. The simulation results showed that the mirror force due to anisotropy reduces ion parallel flow velocity, however, the simulations were presented without cross‐field electric and magnetic drifts. Nevertheless, in present medium‐size tokamaks or future compact fusion devices, cross‐fields are important. In this work, the anisotropy model is extended to be compatible with the inclusion of cross‐field drifts. Preliminary results on comparison between the model solutions with and without cross‐field drifts are presented.

Zhao, Menglong↗
DOE OSTI2025

Sediment transport and morphodynamic changes due to upstream variation and sea level rise on the Western Mekong Delta, Vietnam

The hydrological and morphological dynamics across the Vietnam Mekong Delta (VMD) coastline are increasingly shaped by human interventions and climate change, affecting the sustainability of socio-economic environments. This study employs the TELEMAC-MASCARET modeling system (TELEMAC-2D, TOMAWAC, and SISYPHE modules) to simulate hydrodynamics, sediment transport, and morphological changes for a 10-year period under baseline (SC1) and future climate (SC2) scenarios, considering wave climatology of Northeast (NEM) and Southwest monsoons (SWM). The results of sedimentation and erosion were analyzed across five key coastal zones from Ca Mau Cape (KV1) to Ha Tien (KV5). Model results highlight distinct seasonal patterns of sediment redistribution shaped by monsoonal flows, wave action, and seabed morphology. During the NEM, sediment deposition is concentrated near estuarine zones, while the SWM favors offshore sediment accumulation due to perpendicular wave actions, and shallow, gently sloping bathymetry that dissipates wave energy. Coastal erosion predominantly occurs at convex shoreline features subjected to stronger coastal currents, particularly in fine-grained sediment areas. Under future sea level rise and land subsidence scenarios, reduced flow velocities lead to lower bottom shear stress, promoting deposition and mitigating erosion problems despite increased wave energy. Besides that, sediment transport analyses across 14 cross-sections (C-S) reveal a southward trend during the NEM and a northward reversal during the SWM, with upstream sediment reductions primarily influencing near-estuarine zones and exerting minimal impact beyond the mid-coastal regions. Furthermore, these findings advance the understanding of sediment redistribution processes under climate change and seasonal change, and provide a scientific basis for adaptive coastal management and sustainable development.

Climate change↗
DOE OSTI2022

Oxidation Lifetime Modeling of 625 and 120 Foils After Long-term Exposure in Flowing Air + 10% ${\hbox {H}_2\hbox {O}}$ at 700 and 800 $^\circ$C

A model describing oxidation kinetics considering simultaneous oxidation and volatilization of external oxide scales was integrated with a CALPHAD-based coupled thermodynamic and kinetic approach to predict the oxidation-induced lifetime of foils of alloys 625 and 120. Long-term exposures in flowing air +10% H 2 O at 700 and 800 °C under two flow velocities were used to validate the modeling results. The time to a critical Cr concentration of 10 wt% at the oxide/alloy interface, t 10 , was shown to correlate with the breakdown of the Cr 2 O 3 scale for the studied experimental conditions. Finally, the modeling approach was extended to consider time-dependent oxidation and volatilization rate constants to develop a more appropriate lifetime criterion for microturbine recuperator foils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Transport Upscaling under Flow Heterogeneity and Matrix-Diffusion in Three-Dimensional Discrete Fracture Networks

For this work, we investigate the combined effects of network scale flow variability and retention due to matrix-diffusion on the scaling behavior of transport through fractured media. Two of the principal mechanisms controlling the transport of solutes through fractured low-permeability media are broad distributions of flow velocities and retention times in the solid matrix. We study the relative impact of these two processes under different initial conditions using a set of three-dimensional discrete fracture network simulations. We use these simulations to develop and calibrate an upscaled continuous time random walk (CTRW) approach for advective transport based on an Ornstein-Uhlenbeck model for the particle velocities that accounts for the fracture-matrix coupling using a compound Poisson process. This CTRW model can be conditioned on the initial solute distribution and allows to observe late-time scaling behavior at distances beyond what is feasible to observe using high-fidelity direct numerical simulations. We determine that the initial distribution of particles leads to marked differences in the persistent long-term scale behavior in the solute travel time distributions, even those undergoing retention due to matrix diffusion through implementation and analysis of the model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗