Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW VELOCITY - DISTRIBUTION”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Noninvasive interrogation of local flow phenomena in twisted tape swirled flow via positron emission particle tracking ($\mathrm{PEPT}$)

Twisted tape inserts are often used to enhance the performance of heat exchangers. Previous simulations suggest the existence of secondary vortices and modified axial velocity profiles in these swirled flows, with only limited experimental evidence for these predictions. Here in this work, positron emission particle tracking (PEPT) is used to interrogate turbulent twisted tape swirl flow in a pipe. PEPT is a noninvasive, radiotracer based flow measurement technique that does not require optical access, making it suitable for measurements in a number of engineering flow systems. The existence of the modified velocity profile and secondary flows is confirmed using PEPT, and flow is seen to still be developing 20 diameters downstream of the twisted channel entrance. Secondary vortices are further confirmed by the presence of increased vorticity in these regions, and turbulent fluctuations are seen to increase in magnitude near these flow structures. These results are in line with predictions and point to possible regions of reduced heat transfer and hot spots in twisted tape based heat exchangers. Suggestions are made for future experiments to investigate heated flows, utilizing the ability of PEPT to image flows in opaque systems for validating simulation of flows in engineering equipment, such as nuclear heat transport equipment.

42 ENGINEERING↗
DOE OSTI2021

Pulse Detonation Engine for Advanced Oxy-Combustion of Coal-Based Fuel for Direct Power Extraction Applications

Pressure gain combustion and magnetohydrodynamic (MHD) systems have the potential to provide a step increase in the efficiency of combined-cycle power plants. Specific advantages include a net pressure increase to the system instead of a pressure drop, the high temperature of the detonation waves can increase the efficiency of power extraction compared to other processes, significant thermal energy can be released in a compact region, and the high velocities of the flow increase extraction of electrical power. In summary, a pressure gain combustor coupled with a MHD has the potential to be transformative. Despite the potential advantages, relatively little research has been conducted considering coupled pressure gain combustion systems with MHD systems. With this background and motivation, the overall goal of this effort was to advance the knowledge, technology, and computational tools associated with coupled detonation and MHD systems. A joint experimental and computational approach was used while seeking to accomplish the goals of this work. Specifically, two pulse detonation engines were developed and used for the experiments to produce detonations. Detonation speeds were measured for a variety of flow and fuel conditions (e.g., methane, propane, with coal particles). Preliminary electrical conductivity measurements were collected. An extensive amount of research was performed to identify sensitivities of detonation behaviors to the presence of combustion products. Computationally, a twofold approach was used in this work. First, a solver was developed for solving the governing equations for a reactive flow with coupled detonation and ionization chemistry. The solver was applied to study the impacts of seed material ionization on detonation. Second, a conservation element-solution element (CE-SE) based numerical solver for detonation studies with a reduced reaction mechanism for oxy-methane combustion was developed and verified on standard test cases. Key findings and contributions from this work are as follows. A system was developed for injecting powderized coal, or other seeding material, into a pulse-detonation system. The influence of a combustion product (i.e., CO 2 ) on detonation behavior was identified. Knowledge gained from this work is applicable to devices such as rotating detonation engines, where combustion products mix with fresh reactants. A system for measuring the electrical conductivity of the exhaust from a pulse-detonation engine was developed. The open-source solver, Clawpack, was extended to solve the reactive Euler equations for simulating detonations. A coupled combustion and ionization chemistry was developed in a single chemical kinetic model for methane oxidation. This model can be used to solve coupled MHD and detonation simulations. It was found that parasitic interactions from ionization chemistry with the magnetic field can reduce the detonation velocity by up to 8%, with a potential impact on power extraction of 15%. It is recommended that interactions between the detonation front and MHD field be considered. Higher gas temperatures and velocities were achieved owing to oxy-fuel detonations. Use of radical dissociation reactions in the reduced reaction mechanism, was found to be critical in predicting detonation temperature and velocity accurately.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Astronomy & Astrophysics↗
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 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 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 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 OSTI2021

Hydrokinetic Energy Conversion by Flow-Induced Oscillation of Two Tandem Cylinders of Different Stiffness

Abstract The vortex-induced vibration for aquatic clean energy (VIVACE) converter harnesses hydrokinetic energy by enhancing flow-induced oscillations (FIOs) of elastically supported rigid cylinders in a river, tide, or ocean current. The harnessing power depends on the intensity of the oscillation, which is a consequence of the flow–structure interaction. The inflow condition for the downstream (second) cylinder is slowed down and perturbed by the upstream (first) cylinder, due to the shielding effect. Therefore, the optimal structural parameters, i.e., stiffness and damping ratio, for the second cylinder may be different from the first cylinder, in terms of energy harnessing. To improve the performance of the VIVACE converter, a series of experiments are conducted in a recirculating water channel, with various stiffness combinations of two cylinders in tandem. Results show that the stiffness of the second cylinder, K2, does not affect the energy harnessing power in vortex-induced vibration (VIV) occurring at low speeds, because the oscillation of the downstream cylinder in this velocity range is completely dominated by the wake of the upstream cylinder. K2 has a great influence on the harnessing power at higher velocities in the transition region from VIV to galloping and in galloping. Changing K2 onsets and enhances galloping at lower flow velocity and harnesses up to 110% more energy than the case of K1 = K2.

Engineering↗