Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transitional flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

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

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗
DOE OSTI2025

Systems Analysis and Optimization of Circular PET Packaging Supply Chains in the United States: Environmental and Socioeconomic Impacts

Many actions are underway at global, national, and local levels to address the plastic waste problem and transition toward a circular economy of plastics. Studies evaluating environmental and socioeconomic impacts of such a transition are lacking. Here, the purpose of this study is to conduct a national systems analysis of polyethylene terephthalate (PET) packaging supply chains in the United States. Material flow data was combined with environmental and socioeconomic indicators to evaluate and compare the sustainability of the linear PET packaging supply chain, current (2019) supply chain, and possible future circular supply chain options in the United States. Environmentally optimal circular US PET packaging material flows showed 31% and 38% savings of GHG emissions and energy demand, respectively, with a circularity of 77% when compared with a linear supply chain. Additionally, the environmentally optimal system showed higher employment (29%) and wages (31%) than a linear system, but with a 5% decrease in revenue generation. A socioeconomically optimal circular PET supply chain showed increased employment (by 52%), wages (by 67%), and revenues (by 1%), with a circularity of 59% when compared with the linear system. However, it showed 14% higher GHG emissions than a linear system, indicating a trade-off between environmentally and socioeconomically optimal circular PET packaging systems. Overall, linear-to-circular material flow transition may not necessarily lead to increased revenues and decreased environmental impacts of the entire system, but it does benefit society due to increased employment and wages. Future systems analysis work should focus on improving data quality for environmental and socioeconomic dimensions.

PET↗
DOE OSTI2022

Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?

Creep due to ice flow is generally thought to be the main cause for the formation of crystallographic preferred orientations (CPOs) in polycrystalline anisotropic ice. However, linking the development of CPOs to the ice flow history requires a proper understanding of the ice aggregate's microstructural response to flow transitions. In this contribution the influence of ice deformation history on the CPO development is investigated by means of full-field numerical simulations at the microscale. We simulate the CPO evolution of polycrystalline ice under combinations of two consecutive deformation events up to high strain, using the code VPFFT (visco-plastic fast Fourier transform algorithm) within ELLE. A volume of ice is first deformed under coaxial boundary conditions, which results in a CPO. The sample is then subjected to different boundary conditions (coaxial or non-coaxial) in order to observe how the deformation regime switch impacts the CPO. The model results indicate that the second flow event tends to destroy the first, inherited fabric with a range of transitional fabrics. However, the transition is slow when crystallographic axes are critically oriented with respect to the second imposed regime. Therefore, interpretations of past deformation events from observed CPOs must be carried out with caution, particularly in areas with complex deformation histories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Challenge Problem 1: Preliminary Model Development and Assessment of Flexible Heat Transfer Modeling Approaches

This report presents the modeling progress of a university consortium studying flexible modeling for turbulent heat transfer. In the context of advanced nuclear reactor design, where challenges include non-unity Prandtl fluids, free and mixed convection, and transitional flow, computational fluid dynamics methods are underutilized due to low confidence in modeling approaches and high computational cost. This report evaluates the performance of engineering CFD models in diabatic, buoyant turbulent flow. It finds that all turbulence models including several variants of the k-ε and k-ω models struggle to predict accurate turbulent momentum and heat transfer in such flows. The Nusselt numbers have been compared between the models and the DNS data, where calculations have been performed for each case and trends have shown a good agreement between DNS estimated Nusselt numbers and available correlations as well as experimental data. A novel DNS correlation for high Pr cases as they are transitioning from mixed convection to natural convection has been proposed. While Nusselt number errors relative to DNS range from 20% to 50%, the models capture similar trends to DNS with respect to Nusselt suppression and amplification under varying levels of buoyancy effect. This report also highlights model form error as a significant contributor to CFD predictions and proposes a framework for quantifying model error and improving confidence in CFD calculations. To improve the predictive capability of engineering CFD models, data-driven approaches for turbulence models are investigated. Theoretical frameworks based on the invariant tensor / vector basis neural networks for prediction of Reynolds stresses and turbulent heat fluxes are employed. The models are developed using direct numerical simulations data for forced convection flows of different fluids in vertical planar channel domain. The framework is implemented in spectral element solvers Nek5000 / nekRS and has shown a potential for future development and consideration of mixed convection flows.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗
DOE OSTI2025

Uniting Surface Properties With Hydrodynamic Roughness in Shallow Overland Flow Models

Describing flow resistance from the properties of an underlying surface is a challenge in surface hydrology. Runoff models must specify a resistance formulation or “roughness scheme”—describing the functional relationship between flow resistance and flow depth/velocity—and its parameters. Uncertainty in runoff predictions derives from both the selected roughness scheme (e.g., Darcy Weisbach, Manning's, or laminar flow equations), and its parameterization with a roughness coefficient (e.g., Manning's n). Both choices are informed by model calibration to data, usually discharge, and, if available, velocity. In this study, a Saint Venant Equation-based runoff model is calibrated to discharge and velocity data from 112 rainfall simulator experiments. The results are used to identify the optimal roughness scheme among four widely-used options for each experiment, and to explore whether surface properties can be used to select the optimal roughness scheme and its coefficient. Among the tested roughness schemes, a transitional flow equation provided the best fit to the plurality of experiments. The most suitable roughness scheme for a given experiment was not related to measured surface properties. Regression models predicted the calibrated roughness coefficients with adjusted values between 0.48 and 0.54, depending on the roughness scheme used. Litter cover was the best predictor of the roughness coefficient, followed by soil cover and average canopy gap size. The results suggest that selection of an optimal roughness scheme based on surface properties alone remains difficult, but that once a scheme is selected, roughness coefficients can be estimated from surface properties.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2023

Adaptations for gas exchange enabled the elongation of lepidopteran proboscises

The extensive biodiversification of butterflies and moths (Lepidoptera) is partly attributed to their unique mouthparts (proboscis [Pr]) that can span in length from less than 1 mm to over 280mm in Darwin's sphinx moths. Lepi-doptera, similar to other insects, are believed to inhale and exhale respiratory gases only through valve-like spi-racles on their thorax and abdomen, making gas exchange through the narrow tracheae (Tr) challenging for the elongated Pr. How Lepidoptera overcome distance effects for gas transport to the Pr is an open question that is important to understanding how the Pr elongated over evolutionary time. Here, we show with scanning electron microscopy and X-ray imaging that distance effects on gas exchange are overcome by previously unreported micropores on the Pr surface and by superhydrophobic Tr that prevent water loss and entry. We find that the density of micropores decreases monotonically along the Pr length with the maxima proportional to the Pr length and that micropore diameters produce a Knudsen number at the boundary between the slip and transition flow regimes. By numerical estimation, we further show that the respiratory gas exchange for the Pr predominantly occurs via diffusion through the micropores. In conclusion, these adaptations are key innovations vital to Pr elongation, which likely facilitated lepidopteran biodiversification and the radiation of angiosperms by coevolutionary processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A Review of Infrasound and Seismic Observations of Sample Return Capsules since the End of the Apollo Era in Anticipation of the OSIRIS-REx Arrival

Advancements in space exploration and sample return technology present a unique opportunity to leverage sample return capsules (SRCs) towards studying atmospheric entry of meteoroids and asteroids. Specifically engineered for the secure transport of valuable extraterrestrial samples from interplanetary space to Earth, SRCs offer unexpected benefits that reach beyond their intended purpose. As SRCs enter the Earth’s atmosphere at hypervelocity, they are analogous to naturally occurring meteoroids and thus, for all intents and purposes, can be considered artificial meteors. Furthermore, SRCs are capable of generating shockwaves upon reaching the lower transitional flow regime, and thus can be detected by strategically positioned geophysical instrumentation. NASA’s OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC is one of only a handful of artificial objects to re-enter the Earth’s atmosphere from interplanetary space since the end of the Apollo era and it will provide an unprecedented observational opportunity. This review summarizes past infrasound and seismic observational studies of SRC re-entries since the end of the Apollo era and presents their utility towards the better characterization of meteoroid flight through the atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Scalable Simulation of Pressure Gradient-Driven Transport of Rarefied Gases in Complex Permeable Media Using Lattice Boltzmann Method

Accurate representations of slip and transitional flow regimes present a challenge in the simulation of rarefied gas flow in confined systems with complex geometries. In these regimes, continuum-based formulations may not capture the physics correctly. This work considers a regularized multi-relaxation time lattice Boltzmann (LB) method with mixed Maxwellian diffusive and halfway bounce-back wall boundary treatments to capture flow at high Kn. The simulation results are validated against atomistic simulation results from the literature. We examine the convergence behavior of LB for confined systems as a function of inlet and outlet treatments, complexity of the geometry, and magnitude of pressure gradient and show that convergence is sensitive to all three. The inlet and outlet boundary treatments considered in this work include periodic, pressure, and a generalized periodic boundary condition. Compared to periodic and pressure treatments, simulations of complex domains using a generalized boundary treatment conserve mass but require more iterations to converge. Convergence behavior in complex domains improves at higher magnitudes of pressure gradient across the computational domain, and lowering the porosity deteriorates the convergence behavior for complex domains.

42 ENGINEERING↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part II: The Effect of Roughness on the Performance of a High-Pressure Turbine Blade

Abstract Turbine blades generally present surface roughness introduced in the manufacturing process or caused by in-service degradation, which can have a significant impact on aero-thermal performance. A better understanding of the fundamental physical mechanisms arising from the interaction between the roughness and the turbine flow at engine-relevant conditions can provide insights for the design of blades with improved efficiency and longer operational life. To this end, a high-fidelity numerical framework combining a well-validated solver for direct numerical simulation and a second-order accurate immersed boundary method is employed to predict roughness-induced aero-thermal effects on an LS89 high-pressure turbine (HPT) blade at engine-relevant conditions. Different amplitudes and distributions of surface roughness are investigated and a reference smooth-blade simulation under the same flow conditions is conducted for comparison. Roughness of increasing amplitude progressively shifts the blade suction side boundary layer transition upstream, producing larger values of the turbulent kinetic energy and higher total wake losses. The on-surface data-capturing capabilities of the numerical framework provide direct measurements of the heat flux and the skin friction coefficient, hence offering quantitative information between the surface topology and engineering-relevant performance parameters. This work may provide a benchmark for future numerical studies of turbomachinery flows with roughness.

Engineering↗
DOE OSTI2023

Transition to turbulence in viscoelastic channel flow of dilute polymer solutions

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

Mechanics↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Non-universal current flow near the metal-insulator transition in an oxide interface

In systems near phase transitions, macroscopic properties often follow algebraic scaling laws, determined by the dimensionality and the underlying symmetries of the system. The emergence of such universal scaling implies that microscopic details are irrelevant. Here, we locally investigate the scaling properties of the metal-insulator transition at the LaAlO 3 /SrTiO 3 interface. We show that, by changing the dimensionality and the symmetries of the electronic system, coupling between structural and electronic properties prevents the universal behavior near the transition. By imaging the current flow in the system, we reveal that structural domain boundaries modify the filamentary flow close to the transition point, preventing a fractal with the expected universal dimension from forming.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Macrovoid resolved simulations of transport through HPRO relevant membrane geometries

Modeling the transport properties such as diffusivity and permeability of high pressure reverse osmosis (HPRO) membranes is critical for the selection and manufacture of membranes suitable for operation under high pressure. These properties can be significantly affected by the changes in heterogeneous pore structures due to compaction. Here, the modeling platform presented in this work resolves the two scale porosity in HPRO relevant membranes. A synthetic membrane geometry is constructed based on available experimental visualizations of pore structures. The simulations directly capture the flow channeling that results from a combination of material properties and geometric features of the macrovoids. A parametric study is presented to account for the transition of flow characteristics from a material governed regime to a macrovoid governed regime. Permeability of the membrane is evaluated using the simulation data and compared with an existing model that scales with the square of tortuosity over a range of material properties. The model is found to perform well under a narrow range of tortuosity while deviating from the calculated permeabilities at extreme conditions. The effective permeability of the membrane is found to vary by at least two orders of magnitude between the two flow regimes. It is also observed that tortuosity is a bounded property with its upper limit determined by the macrovoid geometry. Consequently, the tortuosity based correlations fail near a flow regime that is mainly governed by the macrovoids. The modeled permeability can be more than an order of magnitude smaller than the simulation result. A new model based on flux-weighted porosity of a membrane is introduced and its correlation with tortuosity is studied. The model agrees with the simulated data as, in addition to tortuosity, it also accounts for the flux partition within and outside the flow channels. Such correlations enable extending the existing understanding of flow characteristics to enhance predictability of porous media models.

42 ENGINEERING↗
DOE OSTI2021

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗