Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW - PLATES”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Catalytic converters having non-linear flow channels

Disclosed is a honeycomb catalyst substrate core having geometrically non-linear flow channels. In an embodiment, the honeycomb catalyst substrate core includes helical flow channels. In another embodiment, the honeycomb catalyst substrate core includes sinusoidal flow channels. In yet another embodiment, the honeycomb catalyst substrate core includes helical plus sinusoidal flow channels. The honeycomb catalyst substrate core comprises a plurality of parallel non-linear flow channels formed along a longitudinal axis of symmetry of the catalyst substrate core, each non-linear flow channel configured such that a turbulent vortical flow occurs during engine exhaust gas flow. Also disclosed is a method for manufacturing a ceramic honeycomb having non-linear flow channels, comprising the steps extrusion soft ceramic material through a die whilst the die moves through six degrees of freedom along its axis of symmetry. Disclosure includes a method for manufacturing a ceramic honeycomb having non-linear flow channels using three-dimensional printing.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

CTF Improved Drag Model and Flow Regime Transition Criteria

The demand for accurate prediction of two-phase flow behavior in a boiling water reactor (BWR) requires a comprehensive understanding of flow regime, void fraction, heat transfer, and pressure drop. The CTF subchannel code, which is used for the Thermal/Hydraulic (T/H) solution in the Consortium for Advanced Simulation of Light Water Reactors (CASL)-developed Virtual Environment for Reactor Application (VERA) core simulator, is being further developed for BWR applications. In support of this goal, the present work highlights some of the two-phase closure model developments towards improving the CTF void fraction prediction, especially for subcooled boiling. The drift-flux approach has been well-developed for upward dispersed two-phase flows and proven to be accurate in predicting void fraction in bubbly and slug flow regimes. In this work, these kinematic constitutive relations for the drift-flux velocity have been implemented into CTF to describe the interfacial drag of bubbly flow as an alternative to the existing model for better void fraction prediction. The success of these constitutive relations also relies on a good flow regime map that accounts for flow conditions and channel geometry. A more reliable flow regime transition criteria that account for the flow condition has also been implemented in this study for modeling the flow regime transition criteria. The newly implemented models are shown to give improved void fraction predictions in comparison to experimental data.

Hizoum, Belgacem↗
DOE OSTI2023

Flow control valve

Disclosed is a flow device including an inlet, an outlet, and a plurality of fluid flow paths hydraulically connected in parallel to the inlet and the outlet, wherein the plurality of fluid flow paths forms a first ring of fluid flow paths circumferentially arranged at a first radial distance from a centerline of the fluid flow device a second ring of fluid flow patch circumferentially arranged at a second radial distance from the centerline of the fluid flow device, each of the plurality of fluid flow paths has a first hydraulic resistance in a forward flow direction and a second hydraulic resistance in a reverse flow direction, and the second hydraulic resistance is greater than the first hydraulic resistance.

Siebert, Brett W.↗
DOE OSTI2022

Wedge-Shaped Hopper Design for Milled Woody Biomass Flow

The bioenergy industry has been challenged by unstable flow and transport of milled biomass in material handling operations. Handling issues such as hopper clogging and auger jamming are attributed to knowledge gaps between existing handling units designed for bulk solids and their suitability for milled biomass with high compressibility. This work investigates various flow behaviors of granular woody biomass in wedge-shaped hoppers. Herein, hopper flow physical experiments and numerical simulations are conducted to study the influence of the critical material attributes and critical processing parameters on the flow pattern, arching, and throughput. The results show that (1) the preferred flow pattern, mass flow, can be achieved by controlling the material’s internal friction angle, hopper inclination, and hopper wall friction; (2) hopper arching, governed by the competing gravity-driven force against flow resistance from material internal friction and material–wall friction, can be controlled by the hopper wall friction angle and the inclination angle; and (3) flow throughput can be accurately estimated from our empirical equation with inputs of hopper outlet geometry and particle-scale to bulk-scale material attributes. This study elucidates woody biomass flow physics and provides guidance for industrial equipment design.

09 BIOMASS FUELS↗
DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
DOE OSTI2021

Predicting Fluid Flow Regime, Permeability, and Diffusivity in Mudrocks from Multiscale Pore Characterisation

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H 2 , CH 4 ) or permanent (CO 2 , radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2 nm–5 μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size-dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size-dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of particle size and moisture on flow performance of loblolly pine anatomical fractions: Experimental findings and model predictions

The rising energy demand has highlighted biomass as a promising next-generation energy source. However, commercializing biomass-derived energy faces challenges, particularly in handling biomass feedstock. Factors like particle size, shape, moisture content, and surface roughness significantly impact biomass flowability. This study addresses a crucial knowledge gap by examining the effects of particle size and moisture content on the flow behavior and shear properties of different anatomical fractions of loblolly pine (Pinus taeda). The bulk shear behavior was examined using a Schulze ring shear tester, while flow performance was tested through gravity-driven flow experiments in a variable wedge-shape hopper. Results were incorporated into empirical and machine learning-based flow prediction models to evaluate their accuracy and limitations. The study found that samples with higher moisture content show higher unconfined yield strength. The critical arching distance increased with particle size, e.g., from approximately 13 and 33 mm for 2- and 6-mm whole chips, respectively at a 32-degree inclination angle. Conversely, the flow rate decreased for a given hopper opening as particle size increased. For instance, at a 60-mm hopper opening and a 32-degree inclination angle, the mass flow rates for 2- and 6-mm whole chips were 7.83 and 6.42 tonne/h, respectively. The empirical model consistently overpredicted the mass flow rate for all anatomical fractions, while the machine learning model more accurately predicted the central tendency of flow rate but was insensitive to varying tissue proportions. These novel findings provide comprehensive characterization of anatomical fractions, reveal significant combined effects of particle size and moisture content on biomass flow behavior, and demonstrate a better predictive accuracy of a machine learning model, all of which are useful for optimizing material handling strategies and biomass utilization technologies in the industry.

09 - BIOMASS FUELS↗
DOE OSTI2022

Three-dimensional realizations of flood flow in large-scale rivers using the neural fuzzy-based machine-learning algorithms

Machine learning methods have been extensively used to study the dynamics of complex fluid flows. One such algorithm, known as adaptive neural fuzzy inference system (ANFIS), can generate data-driven predictions for flow fields, but has not been applied to natural geophysical flows in large-scale rivers. Herein, we demonstrate the potential of ANFIS to produce three-dimensional (3D) realizations of the instantaneous flood flow field in several large-scale, virtual meandering rivers. The 3D dynamics of flood flow in large-scale rivers were obtained using large-eddy simulation (LES). The LES results, i.e., the 3D velocity components, were employed to train the learnable coefficients of an ANFIS. Further, the trained ANFIS, along with a few time-steps of LES results (precursor data) were then used to produce 3D realizations of flood flow fields in large-scale rivers with geometries other than the one the ANFIS was trained with. We also used the trained ANFIS to generate 3D realizations of river flow at a discharge other than that the ANFIS was trained with. The flow field results obtained from ANFIS were validated using separate LES runs to assess the accuracy of the 3D instantaneous realizations of the machine learning algorithm. An error analysis was conducted to quantify the discrepancies among the ANFIS and LES results for various flood flow predictions in large-scale rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Feasibility of using thermal response methods for nonintrusive compressed air flow measurement

Robust, low-cost nonintrusive flow meters are of interest in many industries. Particularly, a reliable nonintrusive flow measurement for the diagnosis of air leaks in compressed air systems is desirable. Measurement of the air flow due to leaks in the system ensures an accurate estimation of potential cost and energy savings. This study evaluates a novel method of using thermal responses to nonintrusively measure leakage rates in compressed air lines. The method uses heat and the resulting thermal response to calculate the flow rate inside the compressed air line. Compared to the current methods for flow measurement, this method can simplify flow measurement while decreasing the sensitivity to errors when measuring flow rates. In this study, the methodology of the proposed method is explained along with the potential advantages to the design. Two approaches are evaluated: a dynamic step response and sinusoidal frequency response. Simulated tests evaluate the feasibility of the proposed methods, followed by experiments that validate the simulation results. A clear correlation between the thermal step response and the flow rate indicate viability of the proposed method in simulation. Experimental results yielded similar results, confirming the validity of the proposed method. The results of a field test in an industrial environment demonstrate the capability of the approach to other flow rate measurement techniques.

42 ENGINEERING↗
DOE OSTI2020

Deep learning interfacial momentum closures in coarse-mesh CFD two-phase flow simulation using validation data

Multiphase flow phenomena have been widely observed in the industrial applications while it remains a challenging yet unsolved problems. Three-dimensional computational fluid dynamics (CFD) approaches resolve the flow fields on a finer special and temporal scales which can complement the dedicated experimental study. However, closures have to be introduced to reflect the underlying physics in multiphase flow. Among them, the interfacial forces, including drag, lift, turbulent dispersion and wall lubrication forces, play in important role on the bubble’s distribution and migration in liquid-vapor two-phase flow. Development of those closures traditionally rely on the experimental data and analytical derivation with simplified assumptions which usually cannot deliver a universal solution across wide range of flow conditions. In this paper, a data-driven approach, named as Feature Similarity Measurement (FSM), is developed and applied to improve the simulation capability of two-phase flow with coarse-mesh CFD approach. Interfacial momentum transfer in adiabatic bubbly flow serves as the focus of the present study. Both a mature and a simplified set of interfacial closures are taken as the low fidelity data. Experimental data and fine mesh CFD simulations results are adopted as high-fidelity data. Qualitative and quantitative analysis are performed in this paper which reveals that FSM can substantially improve the prediction of coarse mesh CFD model regardless of the choice of interfacial closures and it provides scalability and consistency across discontinuous flow regimes. Furthermore, it demonstrates that data-driven method can aid the multiphase flow modeling by exploring the connections between local physical features and simulation errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) for engineering flow measurements

Increasingly, fully 3D experimental measurements of flow in complex engineering geometries are required to validate computational fluid dynamics models that support and inform reactor design and licensing. One barrier to such measurements is the complexity of typical reactor components and subsequent lack of optical access in these systems. To overcome this, the deployment of coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) is explored for flow measurement in reactor thermal hydraulic components and model (scaled) systems. Through this methodology, fully 3D flow information (via PEPT) and detailed internal geometry (via CT) are captured in opaque systems such as pipes, rod bundles, packed beds, etc. The reconstructed flow field and geometry can then be overlain to reveal detailed flow features around internal structures within a given test section. This is enabled through the use of a combined preclinical PET/CT scanner with overlapping PET and CT fields of view. Such measurements are useful for characterizing flow inside such intricate nuclear thermal hydraulic components as core geometries and heat exchangers, among others, and providing valuable 3D validation data for CFD models. In this work, basic tests of this 3D flow/geometry mapping are presented, and the implications of such measurements are discussed. Further, preliminary measurements are made with both point sources and flow in a simple pipe flow geometry to evaluate the capabilities of this technique. PEPT and CT features are coregistered with up to 0.1 mm precision, and pipe flow mean velocity and Reynolds stresses are reconstructed with similar accuracy to previous PEPT demonstrations. The utility of PEPT/CT is shown herein, and suggestions for future measurements are made.

3D flow measurement↗
DOE OSTI2025

Effects of Particle Mixing and Gravitational Settling on Charge Transport in Carbon Flow-Electrode Cells

Flowable carbon slurries are actively studied and under development for charge transport in various electrochemical systems including flow capacitors, capacitive deionization cells, semi-solid flow batteries, and lithium extraction. However, much less is known about in operando slurry flow dynamics and their corresponding effect on charge transport. We performed an experimental study of mixing and settling dynamics of slurry electrodes within an electrochemical flow cell during continuous operations. The electrochemical cell consisted of two horizontal co-flowing channels, separated by a cation-exchange membrane (CEM). Here we used high-speed optical imaging of planes parallel to gravity and simultaneous electrochemical measurements. At low flow rates, dense yet dynamic particle beds formed on the bottom electrode in each channel, which unexpectedly yielded the highest currents. This approach enables the operation of the flow cell at low system-average particle concentrations while leveraging gravity-driven particle settling to locally enhance carbon concentrations precisely at the current collector sites. Conversely, high flow rates were characterized by thin particle beds and well-mixed particle flows. In the latter case, the electrodes in closest proximity (located on either side of the CEM) achieved a current higher than the other electrode pairs. The observations have implications for slurry control and electrode designs in electrochemical systems.

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

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗
DOE OSTI2024

Solvent-Dependent Dynamics of Cellulose Nanocrystals in Process-Relevant Flow Fields

Flow-assisted alignment of anisotropic nanoparticles is a promising route for the bottom-up assembly of advanced materials with tunable properties. While aligning processes could be optimized by controlling factors such as solvent viscosity, flow deformation, and the structure of the particles themselves, it is necessary to understand the relationship between these factors and their effect on the final orientation. In this study, we investigated the flow of surface-charged cellulose nanocrystals (CNCs) with the shape of a rigid rod dispersed in water and propylene glycol (PG) in an isotropic tactoid state. In situ scanning small-angle X-ray scattering (SAXS) and rheo-optical flow-stop experiments were used to quantify the dynamics, orientation, and structure of the assigned system at the nanometer scale. The effects of both shear and extensional flow fields were revealed in a single experiment by using a flow-focusing channel geometry, which was used as a model flow for nanomaterial assembly. Due to the higher solvent viscosity, CNCs in PG showed much slower Brownian dynamics than CNCs in water and thus could be aligned at lower deformation rates. Moreover, CNCs in PG also formed a characteristic tactoid structure but with less ordering than CNCs in water owing to weaker electrostatic interactions. The results indicate that CNCs in water stay assembled in the mesoscale structure at moderate deformation rates but are broken up at higher flow rates, enhancing rotary diffusion and leading to lower overall alignment. Albeit being a study of cellulose nanoparticles, the fundamental interplay between imposed flow fields, Brownian motion, and electrostatic interactions likely apply to many other anisotropic colloidal systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stress effects on flow and transport in three-dimensional fracture networks

We investigate the effects of various external stress regimes on fracture apertures, fluid flow, and solute transport in three-dimensional fracture networks. We use well-established geomechanics equations coupled with discrete fracture network modeling to characterize changes in primary flow paths within a complex network as a function of stress magnitude and orientation. These changes manifest in the alterations of the fluid flow field and are measured in terms of Eulerian and Lagrangian flow observables including solute transport, which is a key problem in many hydrologic applications. Changes in primary flow paths affect the solute transport in the network by promoting anomalously early arrival or long tailing behavior. However, early time arrival is not ubiquitous in anisotropically stressed networks and in most cases there is a delayed arrival of solute, which is attributed to (i) the presence of low-velocity zones normal to the flow direction, (ii) the angle between flow direction and major compressive principal stress directions, and (iii) changes in the primary flow paths (i.e., increases in tortuosity and active network structure). Overall, flows become more channelized in anisotropically stressed fracture networks than in unstressed and isotropically stressed networks.

58 GEOSCIENCES↗