Simulation of Flow Through Variable Permeability Darcy-Brinkman Layers

2021 ◽  
Author(s):  
SamerA Alokaily
Keyword(s):  
Pressure Gradient ◽  
Strain Rate ◽  
Layer Thickness ◽  
Transition Layer ◽  
Darcy Number ◽  
Fluid Viscosity ◽  
Channel Depth ◽  
Clear Fluid ◽  
Variable Permeability ◽  
Permeable Media

Abstract In this paper, coupled parallel flow in a triple layer channel is studied numerically. The channel consists of a clear fluid sandwiched between two Darcy-Brinkman permeable layers of variable porousness. A single binary equation is presented, in which the penetrability within transition porous layers, is portrayed by a nth degree objective capacity. However, because of the absence of explanatory arrangement of the issue, direct numerical simulations are performed in order to give a novel knowledge into the fluid dynamics inside permeable media of variable porousness. These simulations are carried out through utilizing a modified steady state finite volume solver from the open source programming bundle OpenFOAM. After check and approval of the solver and mathematical technique, parametric investigation is acted in which the Darcy number, intensity of the penetrability degree, transition layer thickness, channel depth, fluid viscosity, and pressure gradient vary. The findings of the current study show that velocity increases when: First, the Darcy number, the degree, or the channel depth increases. Second, when the transition layer thickness decreases. Also, strain rate is almost independent of both Darcy number and degree, and nearly doubles when either the thickness of transition layer halves or the channel depth doubles. In addition, velocity and strain rate are found to scale with viscosity and pressure gradient.

Study on the Flow Characteristics of Shengli Oilfield Super Heavy Oil

2017 ◽  
Vol 10 (1) ◽  
pp. 69-78 ◽  
Author(s):  
Wang Shou-long ◽  
Li Ai-fen ◽  
Peng Rui-gang ◽  
Yu Miao ◽  
Fu Shuai-shi
Keyword(s):  
Pressure Drop ◽  
Pressure Gradient ◽  
Heavy Oil ◽  
Flow Characteristics ◽  
Fluid Viscosity ◽  
Linear Flow ◽  
Start Up ◽  
Non Linear ◽  
Core Permeability ◽  
Velocity Pressure

Objective:The rheological properties of oil severely affect the determination of percolation theory, development program, production technology and oil-gathering and transferring process, especially for super heavy oil reservoirs. This paper illustrated the basic seepage morphology of super heavy oil in micro pores based on its rheological characteristics.Methods:The non-linear flow law and start-up pressure gradient of super heavy oil under irreducible water saturation at different temperatures were performed with different permeable sand packs. Meanwhile, the empirical formulas between start-up pressure gradient, the parameters describing the velocity-pressure drop curve and the ratio of gas permeability of a core to fluid viscosity were established.Results:The results demonstrate that temperature and core permeability have significant effect on the non-linear flow characteristics of super heavy oil. The relationship between start-up pressure gradient of oil, the parameters representing the velocity-pressure drop curve and the ratio of core permeability to fluid viscosity could be described as a power function.Conclusion:Above all, the quantitative description of the seepage law of super heavy oil reservoir was proposed in this paper, and finally the empirical diagram for determining the minimum and maximum start-up pressure of heavy oil with different viscosity in different permeable formations was obtained.

The structure of a highly decelerated axisymmetric turbulent boundary layer

2021 ◽  
Vol 929 ◽  
Author(s):  
N. Agastya Balantrapu ◽  
Christopher Hickling ◽  
W. Nathan Alexander ◽  
William Devenport
Keyword(s):  
Boundary Layer ◽  
Pressure Gradient ◽  
Shear Layer ◽  
Layer Thickness ◽  
Boundary Layer Thickness ◽  
Large Scale ◽  
Convection Velocity ◽  
Mean Flow ◽  
Mean Velocity ◽  
The Mean

Experiments were performed over a body of revolution at a length-based Reynolds number of 1.9 million. While the lateral curvature parameters are moderate ( $\delta /r_s < 2, r_s^+>500$ , where $\delta$ is the boundary layer thickness and r s is the radius of curvature), the pressure gradient is increasingly adverse ( $\beta _{C} \in [5 \text {--} 18]$ where $\beta_{C}$ is Clauser’s pressure gradient parameter), representative of vehicle-relevant conditions. The mean flow in the outer regions of this fully attached boundary layer displays some properties of a free-shear layer, with the mean-velocity and turbulence intensity profiles attaining self-similarity with the ‘embedded shear layer’ scaling (Schatzman & Thomas, J. Fluid Mech., vol. 815, 2017, pp. 592–642). Spectral analysis of the streamwise turbulence revealed that, as the mean flow decelerates, the large-scale motions energize across the boundary layer, growing proportionally with the boundary layer thickness. When scaled with the shear layer parameters, the distribution of the energy in the low-frequency region is approximately self-similar, emphasizing the role of the embedded shear layer in the large-scale motions. The correlation structure of the boundary layer is discussed at length to supply information towards the development of turbulence and aeroacoustic models. One major finding is that the estimation of integral turbulence length scales from single-point measurements, via Taylor's hypothesis, requires significant corrections to the convection velocity in the inner 50 % of the boundary layer. The apparent convection velocity (estimated from the ratio of integral length scale to the time scale), is approximately 40 % greater than the local mean velocity, suggesting the turbulence is convected much faster than previously thought. Closer to the wall even higher corrections are required.

scholarly journals Steady MHD Flow of Viscous Fluid between Two Parallel Porous Plates with Heat Transfer in an Inclined Magnetic Field

2015 ◽  
Vol 7 (3) ◽  
pp. 21-31 ◽  
Author(s):  
D. R. Kuiry ◽  
S. Bahadur
Keyword(s):  
Magnetic Field ◽  
Heat Transfer ◽  
Viscous Fluid ◽  
Pressure Gradient ◽  
Reverse Flow ◽  
Flow Behavior ◽  
Fluid Pressure ◽  
Mhd Flow ◽  
Fluid Viscosity ◽  
Temperature Dependent Viscosity

The steady flow behavior of a viscous, incompressible and electrically conducting fluid between two parallel infinite insulated horizontal porous plates with heat transfer is investigated along with the effect of an external uniform transverse magnetic field, the action of inflow normal to the plates, the pressure gradient on the flow and temperature. The fluid viscosity is supposed to vary exponentially with the temperature. A numerical solution for the governing equations for both the momentum transfer and energy transfer has been developed using the finite difference method. The velocity and temperature distribution graphs have been presented under the influence of different values of magnetic inclination, fluid pressure gradient, inflow acting perpendicularly on the plates, temperature dependent viscosity and the Hartmann number. In our study viscosity is shown to affect the velocity graph. The flow parameters such as viscosity, pressure and injection of fluid normal to the plate can cause reverse flow. For highly viscous fluid, reverse flow is observed. The effect of magnetic force helps to restrain this reverse flow.

scholarly journals Observations of the Transition Layer

2009 ◽  
Vol 39 (3) ◽  
pp. 780-797 ◽  
Author(s):  
T. M. Shaun Johnston ◽  
Daniel L. Rudnick
Keyword(s):  
Layer Thickness ◽  
Mixed Layer ◽  
Potential Vorticity ◽  
Transition Layer ◽  
Mixed Layer Depth ◽  
Density Ratio ◽  
Horizontal Resolution ◽  
Surface Mixed Layer ◽  
Layer Depth ◽  
Vertical Displacements

Abstract The transition layer is the poorly understood interface between the stratified, weakly turbulent interior and the strongly turbulent surface mixed layer. The transition layer displays elevated thermohaline variance compared to the interior and maxima in current shear, vertical stratification, and potential vorticity. A database of 91 916 km or 25 426 vertical profiles of temperature and salinity from SeaSoar, a towed vehicle, is used to define the transition layer thickness. Acoustic Doppler current measurements are also used, when available. Statistics of the transition layer thickness are compared for 232 straight SeaSoar sections, which range in length from 65 to 1129 km with typical horizontal resolution of ∼4 km and vertical resolution of 8 m. Transition layer thicknesses are calculated in three groups from 1) vertical displacements of the mixed layer base and of interior isopycnals into the mixed layer; 2) the depths below the mixed layer depth of peaks in shear, stratification, and potential vorticity and their widths; and 3) the depths below or above the mixed layer depth of extrema in thermohaline variance, density ratio, and isopycnal slope. From each SeaSoar section, the authors compile either a single value or a median value for each of the above measures. Each definition yields a median transition layer thickness from 8 to 24 m below the mixed layer depth. The only exception is the median depth of the maximum isopycnal slope, which is 37 m above the mixed layer base, but its mode is 15–25 m above the mixed layer base. Although the depths of the stratification, shear, and potential vorticity peaks below the mixed layer are not correlated with the mixed layer depth, the widths of the shear and potential vorticity peaks are. Transition layer thicknesses from displacements and the full width at half maximum of the shear and potential vorticity peak give transition layer thicknesses from 0.11× to 0.22× the mean depth of the mixed layer. From individual profiles, the depth of the shear peak below the stratification peak has a median value of 6 m, which shows that momentum fluxes penetrate farther than buoyancy fluxes. A typical horizontal scale of 5–10 km for the transition layer comes from the product of the isopycnal slope and a transition layer thickness suggesting the importance of submesoscale processes in forming the transition layer. Two possible parameterizations for transition layer thickness are 1) a constant of 11–24 m below the mixed layer depth as found for the shear, stratification, potential vorticity, and thermohaline variance maxima and the density ratio extrema; and 2) a linear function of mixed layer depth as found for isopycnal displacements and the widths of the shear and potential vorticity peaks.

Buoyancy-driven crack propagation: the limit of large fracture toughness

2007 ◽  
Vol 580 ◽  
pp. 359-380 ◽  
Author(s):  
S. M. ROPER ◽  
J. R. LISTER
Keyword(s):  
Fracture Toughness ◽  
Pressure Gradient ◽  
Viscous Flow ◽  
Numerical Solutions ◽  
Elastic Solid ◽  
Propagation Rate ◽  
Fluid Viscosity ◽  
Viscous Effects ◽  
Vertical Propagation ◽  
Viscous Pressure

We study steady vertical propagation of a crack filled with buoyant viscous fluid through an elastic solid with large effective fracture toughness. For a crack fed by a constant flux Q, a non-dimensional fracture toughness K=Kc/(3μQm3/2)1/4 describes the relative magnitudes of resistance to fracture and resistance to viscous flow, where Kc is the dimensional fracture toughness, μ the fluid viscosity and m the elastic modulus. Even in the limit K ≫ 1, the rate of propagation is determined by viscous effects. In this limit the large fracture toughness requires the fluid behind the crack tip to form a large teardrop-shaped head of length O(K2/3) and width O(K4/3), which is fed by a much narrower tail. In the head, buoyancy is balanced by a hydrostatic pressure gradient with the viscous pressure gradient negligible except at the tip; in the tail, buoyancy is balanced by viscosity with elasticity also playing a role in a region within O(K2/3) of the head. A narrow matching region of length O(K−2/5) and width O(K−4/15), termed the neck, connects the head and the tail. Scalings and asymptotic solutions for the three regions are derived and compared with full numerical solutions for K ≤ 3600 by analysing the integro-differential equation that couples lubrication flow in the crack to the elastic pressure gradient. Time-dependent numerical solutions for buoyancy-driven propagation of a constant-volume crack show a quasi-steady head and neck structure with a propagation rate that decreases like t−2/3 due to the dynamics of viscous flow in the draining tail.

Effect of Horizontal Alternating Current Electric Field on the Stability of Natural Convection in a Dielectric Fluid Saturated Vertical Porous Layer

2015 ◽  
Vol 137 (4) ◽  
Author(s):  
B. M. Shankar ◽  
Jai Kumar ◽  
I. S. Shivakumara
Keyword(s):  
Natural Convection ◽  
Electric Field ◽  
Prandtl Number ◽  
Traveling Wave ◽  
Effective Viscosity ◽  
Darcy Number ◽  
Wave Mode ◽  
Dielectric Fluid ◽  
Fluid Viscosity ◽  
The Stability

The stability of natural convection in a dielectric fluid-saturated vertical porous layer in the presence of a uniform horizontal AC electric field is investigated. The flow in the porous medium is governed by Brinkman–Wooding-extended-Darcy equation with fluid viscosity different from effective viscosity. The resulting generalized eigenvalue problem is solved numerically using the Chebyshev collocation method. The critical Grashof number Gc, the critical wave number ac, and the critical wave speed cc are computed for a wide range of Prandtl number Pr, Darcy number Da, the ratio of effective viscosity to the fluid viscosity Λ, and AC electric Rayleigh number Rea. Interestingly, the value of Prandtl number at which the transition from stationary to traveling-wave mode takes place is found to be independent of Rea. The interconnectedness of the Darcy number and the Prandtl number on the nature of modes of instability is clearly delineated and found that increasing in Da and Rea is to destabilize the system. The ratio of viscosities Λ shows stabilizing effect on the system at the stationary mode, but to the contrary, it exhibits a dual behavior once the instability is via traveling-wave mode. Besides, the value of Pr at which transition occurs from stationary to traveling-wave mode instability increases with decreasing Λ. The behavior of secondary flows is discussed in detail for values of physical parameters at which transition from stationary to traveling-wave mode takes place.

ANTI-BACTERIAL APPLICATIONS FOR NEW THERMAL CONDUCTIVITY MODEL IN ARTERIES WITH CNT SUSPENDED NANOFLUID

2016 ◽  
Vol 16 (05) ◽  
pp. 1650063 ◽  
Author(s):  
NOREEN SHER AKBAR ◽  
M. RAZA ◽  
R. ELLAHI
Keyword(s):  
Magnetic Field ◽  
Pressure Gradient ◽  
Mathematical Formulation ◽  
Pressure Rise ◽  
Darcy Number ◽  
Peristaltic Flow ◽  
Permeable Walls ◽  
Thermal Conductivity Model ◽  
Number Formula ◽  
Physical Quantities

The peristaltic flow of a carbon nanotubes (CNTs) water fluid investigate the effects of heat generation and magnetic field in permeable vertical diverging tube is studied. The mathematical formulation is presented, the resulting equations are solved exactly. The obtained expressions for pressure gradient, pressure rise, temperature, velocity profile are described through graphs for various pertinent parameters. The streamlines are drawn for some physical quantities to discuss the trapping phenomenon. It is observed that pressure gradient profile is decreasing by increase of Darcy number [Formula: see text] because Darcy number is due to porous permeable walls of the tube and when walls are porous fluid cannot easily flow in tube, so that will decrease the pressure gradient.

Non-Newtonian Fluid Flow and Heat Transfer in a Semicircular Microtube Induced by Electroosmosis and Pressure Gradient

2018 ◽  
Vol 140 (12) ◽  
Author(s):  
Mehdi Karabi ◽  
Ali Jabari Moghadam
Keyword(s):  
Nusselt Number ◽  
Pressure Gradient ◽  
Flow Behavior ◽  
Velocity Ratio ◽  
Shear Thinning ◽  
Fluid Viscosity ◽  
Flow Behavior Index ◽  
Shear Thinning Fluids ◽  
Power Law Fluids ◽  
Behavior Index

The hydrodynamic and thermal characteristics of electroosmotic and pressure-driven flows of power-law fluids are examined in a semicircular microchannel under the constant wall heat flux condition. For sufficiently large values of the electrokinetic radius, the Debye length is thin; the active flow within the electric double layer (EDL) drags the rest of the liquid due to frictional forces arising from the fluid viscosity, and consequently a plug-like velocity profile is attained. The velocity ratio can affect the pure electrokinetic flow as well as the flow rate depending on the applied pressure gradient direction. Since the effective viscosity of shear-thinning fluids near the wall is quite small compared to the shear-thickening fluids, the former exhibits higher dimensionless velocities than the later close to the wall; the reverse is true at the middle section. Poiseuille number increases with increasing the flow behavior index and/or the electrokinetic radius. Due to the comparatively stronger axial advection and radial diffusion in shear-thinning fluids, better temperature uniformity is achieved in the channel. Reduction of Nusselt number continues as far as the fully developed region where it remains unchanged; as the electrokinetic radius tends to infinity, Nusselt number approaches a particular value (not depending on the flow behavior index).

Transition layer thickness in microlaminar deposits

1993 ◽  
Vol 23 (6) ◽  
pp. 662-668 ◽  
Author(s):  
A. R. Despić ◽  
T. LJ. Trišović
Keyword(s):  
Layer Thickness ◽  
Transition Layer
Sign in / Sign up

Export Citation Format

Share Document