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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Eddy current flow meter model validation with a moving solid rod*

We report that the eddy current flow meter (ECFM) has been used to measure velocities and temperatures of conductive flows, such as liquid metal flows in a nuclear fission reactor. The goal of this paper is to develop a finite element electromagnetic model that can characterize the ECFM sensor performance and validate this finite element model with detailed velocity measurements of a controlled, well-characterized moving conductive solid rod. Both measurements and modeling were performed for various parameters that are important for ECFM performance such as rod velocity, rod material, ECFM sensor coil length, number of sensor coils, applied alternating current (AC) current amplitude, and applied ECFM AC frequency. For all parametric scans, the measurement and modeling agree well in both magnitude and trend. The normalized root-mean-square error between measurement and modeling is less than 10% for all cases. These results suggest that electromagnetic modeling could eventually be used to cost-effectively design future ECFM sensors in arbitrary geometry for more challenging applications such as liquid metal nuclear fission reactors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Surface-Functionalized Cellulose Nanocrystals as Nanofillers for Crosslinking Processes: Implications for Thermosetting Resins

Understanding the response of fillers in the epoxy resin crosslinking process and characterizing polymer–filler dynamics are the key features that guide the engineering of new thermosetting resin composites. Here, in this work, X-ray photon correlation spectroscopy (XPCS) is used as a thermal analysis tool to track the microscopic changes occurring during the cure of a functionalized cellulose nanocrystal (mCNC)–epoxy composite. In contrast, traditional methods such as differential scanning calorimetry (DSC) and curing rheology are used to understand the kinetics and properties on macroscopic length scales. Of interest is the influence of the mCNC on the curing kinetics and properties of the thermosetting resin. Two levels of modification (increasing hydrophobicity) were chosen to observe the effect of functionalization. Before cure, the highly functionalized CNC (mCNC3) shows a 44% increase in complex viscosity (η*), while the less functionalized CNC (mCNC2) shows a η* value similar to that of the neat resin. As the cure cycle progresses, results from DSC, rheology, and XPCS further show the enhancement in dispersion for mCNC3. The results show a clear difference in the maximum drift velocity, maximum heat flow, and complex viscosity during the ramp to the isothermal cure temperature (T cure ). Such results suggest that mCNC3 contains a well-dispersed network of particles due to the higher level of functionalization. During T cure , a transition in elastic modulus (G') occurs only for the highly functionalized CNC particle system. We believe that heat-induced aggregation occurs, and the crosslinked resin ultimately dominates the macroscopic properties of the final cured system for all samples. The results from the three techniques are in good agreement and showcase XPCS as a beneficial experimental tool for characterizing the microscopic dynamics of particulate-filled thermosetting resins. Hence, we envision this to be a fundamental curing study for the design of thermosetting resin composites.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Response Surface Methodology As a New Approach for Finding Optimal MALDI Matrix Spraying Parameters for Mass Spectrometry Imaging

Automated spraying devices have become ubiquitous in laboratories employing matrixassisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in part because they permit control of a number of matrix application parameters that can easily be reproduced for intra- and interlaboratory studies. Determining the optimal parameters for MALDI matrix application, such as temperature, flow rate, spraying velocity, number of spraying cycles, and solvent composition for matrix application, is critical for obtaining high-quality MALDI-MSI data. However, there are no established approaches for optimizing these multiple parameters simultaneously. Instead optimization is performed iteratively (i.e., one parameter at a time), which is time consuming and can lead to overall non-optimal settings. In this report, we demonstrate the use a novel experimental design and the response surface methodology to optimize five parameters of MALDI matrix application using a robotic sprayer. Thirtytwo combinations of MALDI matrix spraying conditions were tested, which allowed us to elucidate relationships between each of the application parameters as determined by MALDI-MS (specifically, using a 15 Tesla Fourier transform ion cyclotron resonance mass spectrometer). As such, we were able to determine the optimal automated spraying parameters that minimized signal de-localization and enabled high MALDI sensitivity. We envision this optimization strategy can be utilized for other MALDI-MSI applications, tissue types, and matrix application approaches.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Investigation of carbon incorporation in laser-assisted MOCVD of GaN

Background carbon (C) impurity incorporation in metalorganic chemical vapor deposition (MOCVD) grown gallium nitride (GaN) represents one of the major issues in further improving GaN vertical power device performance. This work presents a laser-assisted MOCVD (LA-MOCVD) technique to address the high-C issue in MOCVD homoepitaxial GaN under different growth rate (Rg) regimes and studies the correlations between [C] and Rg. [C] in LA-MOCVD GaN is reduced by 50%–90% as compared to the conventional MOCVD GaN for a wide growth rate range between 1 and 16 μm/h. A mass-transport based model is developed to understand the C incorporation at different Rg regimes. The results obtained from the developed model are in good agreement with experimental data. The model further reveals that LA-MOCVD effectively suppresses C incorporation by reducing the active C species in the gas phase. Moreover, high step velocity in step flow growth mode can facilitate C incorporation at fast Rg, exhibiting steeper C increase. The theoretical model indicates that [C] can be suppressed below 1016 cm−3 with a fast growth rate (Rg) of 10 μm/h by utilizing higher power LA-MOCVD and freestanding GaN substrates with larger off-cut angles.

Physics↗
DOE OSTI2023

Visualization of fast ion phase-space flow in plasmas well-below, near and well-above Alfvén eigenmode stability threshold in tokamak

Abstract Transport of fast ions along certain local phase space paths, referred to as fast ion phase-space flow, has been systematically measured by an imaging neutral particle analyzer (INPA) in three plasma regimes, which are well-below, near, and well-above the Alfvén eigenmode stability threshold (Du et al 2021 Phys. Rev. Lett. 127 235002). (1) In plasmas well-below the Alfvén eignenmodes (AE) stability threshold, fast ions are well-confined on passing particle orbits without noticeable transport over the phase space. The observed INPA images agree well with the synthetic INPA images, using the fast ion distribution predicted by neoclassical theory. (2) In plasmas near the AE stability threshold, INPA images in the presence of AE activity moderately deviate from those without AE activity. The image difference can be well interpreted by AE-driven, phase-space fast ion flow. Paths of this flow over the velocity space (or streamlines) are reconstructed by the intersection lines of curved E ʹ and µ surfaces, referred to as E ʹ and µ line (where E ′ ≡ E − ( ω / n ) P ζ ; E , P ζ and µ are the energy, canonical toroidal momentum and magnetic moment of ions; ω and n are the angular frequency and toroidal mode number of AEs, respectively). Resonant fast ions move radially inward by gaining energy and move radially outward by losing energy and the trajectory well aligns with the E ʹ and µ lines that pass through the mode resonances near the injection energy of neutral beams. (3) In plasmas well-above the AE stability threshold, fast ion phase-space dynamics shows additional features. Fast ions are transported out of the birth positions so promptly along the streamlines that the slowing-down process from the injection energy is not observable, exhibiting strong critical gradient behavior at local phase space. As a result, the increase of electron temperature is very small, in spite of an increase of beam power by ∼ 45 % . It should be emphasized that the directions of phase-space transport, induced by AEs with different frequencies, structures and mode numbers, do not largely differ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Linear Grains

LIGR (Linear Grains) implements a linear (Fourier mode) analysis of a shock propagating through density perturbations. It has two primary capabilities. First, it calculates the shock front distortion and post-shock flow variables (density, velocity, pressure) for single Fourier mode pairs in two dimensions (a pair of x and y modes). Second, it calculates these same quantities for a grains-like pre-shock density perturbation, which consists of a sum of many Fourier modes. This grains-like perturbation is a periodic grid of rectangles at a specifiable density and size (the grains), separated by regions of a different specifiable size and density (interstitial spaces between the grains). The single-mode analysis is based on the work of Velikovich et al., Phys. Plasmas 14 072706 (2007).

Li, GraceJ↗
DOE OSTI2023

Stable water isotopes and tritium data from porewater at Elkhorn Slough

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we collected geochemical tracers across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. The experimental transect is located at the Elkhorn Slough National Estuarine Research Reserve.Here we present a data table (.csv) with the results of electrical conductivity, stable water isotopes, and tritium in porewater across several elevations in the hillslope-to-salt marsh continuum (information in the data set). We analyzed the stable water isotope samples by cavity ring-down spectroscopy in a Picarro L2130-i at the University of California Santa Cruz (accuracy of 0.025‰ and 0.1‰ for δ18O and δ2H, respectively). We measured the electric conductivity of the samples with an Orion Star™ A329 multiparameter meter (0.5% accuracy, Thermo Fisher Scientific, Massachusetts, USA). We analyzed 3H samples at Lawrence Livermore National Laboratory by helium-3 accumulation (Clarke et al., 1976; Surano et al., 1992).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

NMDQi Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST)

The Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST) is developed to support the thermal-hydraulic aspects for the design of optimal irradiation vehicles for the Nuclear Materials Discovery and Qualification initiative (NMDQi). TOAST is currently a MATLAB-based tool that utilizes a simplified steady-state analytical model where a radial thermal resistance circuit is utilized to account for conduction through up to 2 capsules, a gas gap, a thermal bond, a sample cladding, and a sample, as well as the convection from a water coolant outside the capsules. The geometry is discretized axially along the user-defined height of a basket with a user-defined number of points/nodes, essentially turning TOAST into a semi-2D heat transfer analysis utility for cylindrical irradiation vehicles with practically any configuration, solving for temperatures radially at multiple axial locations. TOAST utilizes a Graphical User Interface (GUI) in-which a user can define the geometrical layout, the materials utilized for each component, the coolant characteristics, and the required solution. The user can define the geometry by inputting the diameters, thicknesses, and heights for each component, whereas the materials are defined via the user-provided constant or variable thermal conductivities. The user can select the coolant's inlet temperature, pressure, and the pressure drop across the height of the problem (which is used to calculate the velocity of the flow). Finally, the user can choose to solve for a sample heat generation rate limit by inputting temperature limits for the sample and capsules or can choose to purely solve for the axial temperature distributions of each component by inputting a pre-selected heat generation rate. Either way, TOAST provides the user with axial temperature distributions of all the components including the coolant, and the heat generation rate limit as well as the maximum outlet coolant temperature. The user can also choose to do one of 6 sensitivity analyses in TOAST. The sensitivity analyses yield plots of the sensitivity of the heat generation rate limit, the maximum coolant temperature, and the pressure limit for the annular components due to perturbations in 1-2 unknown variables based on the selected sensitivity analysis. Benchmarks between computations in TOAST and equivalent 2D axis symmetric ABAQUS models are presented, showing that TOAST results are within less than 3% of ABAQUS results in most cases, and a maximum of 8% difference in some cases. The benchmarks also revealed that this uncertainty is tied to the selection of an appropriate Nusselt number correlation and appropriate thermal conductivities, which the user can do from the GUI. Regardless, TOAST is demonstrated as a computationally efficient, highly accessible, and accurate utility for optimization and scoping s of studies different irradiation vehicles.

42 ENGINEERING↗
DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗
DOE OSTI2025

Influence of Reservoir Convection on Heat Extraction with Closed-Loop Geothermal Systems

Over the past few years, advancements in closed-loop geothermal systems (CLGS), also called advanced geothermal systems (AGS), have sparked a renewed interest in these types of designs. CLGS have certain advantages over traditional and enhanced geothermal systems (EGS), including not requiring in-situ reservoir permeability, conservation of the circulating fluid, and allowing for different fluids, including working fluids directly driving a turbine at the surface. CLGS may be attractive in environments where water resources are limited, rock contaminants must be avoided, and stimulation treatments are not available (e.g., due to regulatory or technical reasons). Despite these advantages, CLGS have some challenges, including limited surface area for heat transfer and requiring long wellbores and laterals to obtain multi-MW output in conduction-only reservoirs. CLGS have been investigated in conduction-only systems. In this paper, we explore the impact of both forced and natural convection on the levels of heat extraction with a CLGS deployed in a hot wet rock reservoir. We bound potential benefits of convection by investigating liquid reservoirs over a range of natural and forced convective coefficients. Additionally, we investigate the effects of permeability, porosity, and geothermal temperature gradient in the reservoir on CLGS outputs. Reservoir simulations indicate that reservoir permeabilities of at least ~100 mD are required for natural convection to increase the heat output with respect to a conduction-only scenario. The impact increases with increasing reservoir temperature. When subject to a forced convection flow field, Darcy velocities of at least 10-7 m/s are required to obtain an increase in heat output.

advanced geothermal systems↗
DOE OSTI2024

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Investigation of Anchor Link Magnetohydrodynamic Effects in the Toroidally Symmetric Lead-Lithium (TSLL) Blanket Concept

This paper presents a comprehensive study on the magnetohydrodynamic (MHD) flow in a slotted channel with a cylindrical anchor link, which is an imperative component of the toroidally symmetric lead-lithium (TSLL) blanket concept. Following the validation of the MHD solver implemented in COMSOL and a proposed subtraction approach to compute the anchor link pressure drop, the effects of computational domain size on the pressure drop and velocity distribution are examined. The results show that the pressure drop associated with the anchor link and the maximum flow velocity follow an asymptotic trend as the domain width increases, with wall-induced pressure drop being more dominant than that of the anchor link. The velocity distribution analysis revealed the formation of an internal boundary layer extending along the magnetic field direction, which is a unique feature of the investigated MHD flow. An estimation of the total MHD pressure drop associated with the array of anchor links under the TSLL blanket conditions suggests ~0.3 MPa, which is significantly lower than the recommended maximum allowable blanket pressure drop of 2 MPa. In conclusion, this work offers valuable insights into the anchor link–associated MHD phenomena and serves as a foundation for further development of the TSLL blanket concept for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2021

Analyzing Mechanical integrity of Microwave resonant Cavity flowmeter

The microwave resonant cavity flow sensor has been proposed to address the difficulties of measuring high-temperature molten salt flow in reactor vessels. The flow sensor consists of a metallic cylindrical cavity coupled to a microwave waveguide. The cylindrical cavity has one of the flat surfaces thin enough to undergo microscopic deflection due to dynamic pressure of the exerted by the flowing fluid. Deflection of the membrane causes a shift in the microwave resonant frequency. Fluid velocity can be estimated from the calibration curve. The objective of this report is to investigate mechanical integrity of the flowmeter’s membrane. The membrane is a thin stainless-steel circular plate with a clamped edge, subjected to uniform load. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. If the maximum stress is greater than the yield strength, the plate will undergo plastic deformation, and the sensor will be disabled. A stress greater than the ultimate tensile strength will cause the plate to fracture. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Stresses in both the analytic model and COMSOL numerical model were calculated for stainless steel 316 thin circular plates under room temperature conditions, where the fluid is water. The diameter of the circular plate is approximately one inch, and thickness is 10mil. Fluid flow velocity is in the range from 0.5m/s to 1m/s. The results from the two models are similar. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗
DOE OSTI2023

Characterizing the Impacts of Multi-Scale Heterogeneity on Solute Transport in Fracture Networks

We model flow and transport in fracture networks with varying degrees of fracture-to-fracture aperture heterogeneity and network intensity to show how changes in these properties can cause the emergence of anomalous flow and transport behavior. If fracture-to-fracture aperture heterogeneity is increased in sparse networks, velocity fluctuations inhibit high flow rates and solute transport can be delayed. Surprisingly, transport can be slowed even in cases where hydraulic aperture is monotonically increased. As the intensity of the networks is increased, more connected pathways allow for particles to bypass these effects. There exists transition behavior where with relatively few connected pathways in a network, first arrival times of particles are not heavily affected by fracture-to-fracture aperture heterogeneity, but the scaling behavior of the tails is strongly influenced. These results reinforce the importance of considering multi-scale effects in fractured systems and can inform flow and transport processes in both natural and engineered fractured systems.

58 GEOSCIENCES↗
DOE OSTI2023

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗