Finite element modeling of convective pore-fluid flow in fluid-saturated porous rocks within upper crust:An overview

Convective pore-fluid flow (CPFF) plays a critical role in generating mineral deposits and oil reservoirs within the deep Earth. Therefore, theoretical understanding and numerical modeling of the thermodynamic process that triggers and controls the CPFF are extremely important for the exploration of new mineral deposits and underground oil resources. From the viewpoint of science, the CPFF within the upper crust can be treated as a kind of thermodynamic instability problem of pore-fluid in fluid-saturated porous media. The key issue of dealing with this kind of problem is to assess whether a nonlinear thermodynamic system under consideration is supercritical. To overcome limitations of using theoretical analysis and experimental methods in dealing with the CPFF problems within the upper crust, finite element modeling has been broadly employed for solving this kind of problem over the past two decades. The main purpose of this paper is to overview recent developments and applications of finite element modeling associated with solving the CPFF problems in large length-scale geological systems of complicated geometries and complex material distributions. In particular, two kinds of commonly-used finite element modeling approaches, namely the steady-state and transient-state approaches, and their advantages/disadvantages are thoroughly presented and discussed.

Finite element solution of nonlinear convective flow of Oldroyd-B fluid with Cattaneo-Christov heat flux model over nonlinear stretching sheet with heat generation or absorption

In this study,a two-dimensional boundary layer flow of steady incompressible nonlinear convective flow of Oldroyd-B fluid over a nonlinearly stretching sheet with Cattaneo-Christov heat flux model and heat generation or absorption is examined.The governing equations of the boundary layer flow which are highly nonlinear partial differential equations are converted to the ordinary differential equations using similarity transformations and then the Galerkin finite element method(GFEM)is used to solve the proposed problem.The effect of local Deborah numbers 0,and ft.local buoyancy parameter z,Prandtl number Pr,Deborah number y,and heat generation/absorption parameter<5 on the temperature and the velocity as well as heat transfer rate and shear stress are discussed both in graphical and tabular forms.The result shows the enlargement in the local buoyancy parameter A will improve the velocity field and the heat transfer rate of the boundary layer flow.Moreover,our present work evinced both local skin friction coefficient and heat transfer rate step up if we add the values of non-linear stretching sheet parameter and local heat generation/absorption parameter has quite the opposite effect.The numerically computed values of local skin friction coefficient and local Nusselt number are validated with available literature and evinced excellent agreement.

EFFECTS OF TRANSVERSE OSCILLATIONS ON NATURAL CONVECTIVE FLOW INDUCED BY COMBINED THERMAL AND MASS DIFFUSIONS IN POROUS MEDIA

This paper studies the unsteady heat and mass natural convection in a highly porous medium bounded by an infinite vertical porous wall. The unsteady source of the problem arises from the transverse oscillations in suction velocity of fluids, The analytical results for the problem are obtained based on the method of small parameter, and show that the natural circulation in the porous medium is affected by this kind of oscillation.

Convective heat and mass transfer in MHD mixed convection flow of Jeffrey nanofluid over a radially stretching surface with thermal radiation

Mixed convection flow of magnetohydrodynamic(MHD) Jeffrey nanofluid over a radially stretching surface with radiative surface is studied. Radial sheet is considered to be convectively heated. Convective boundary conditions through heat and mass are employed. The governing boundary layer equations are transformed into ordinary differential equations. Convergent series solutions of the resulting problems are derived. Emphasis has been focused on studying the effects of mixed convection, thermal radiation, magnetic field and nanoparticles on the velocity, temperature and concentration fields. Numerical values of the physical parameters involved in the problem are computed for the local Nusselt and Sherwood numbers are computed.

Peristaltic flow in an asymmetric channel with convective boundary conditions and Joule heating

The peristaltic transport of viscous fluid in an asymmetric channel is concentrated. The channel walls exhibit convective boundary conditions. Both cases of hydrodynamic and magnetohydrodynamic(MHD) fluids are considered. Mathematical analysis has been presented in a wave frame of reference. The resulting problems are non-dimensionalized. Long wavelength and low Reynolds number approximations are employed. Joule heating effect on the thermal equation is retained. Analytic solutions for stream function and temperature are constructed. Numerical integration is carried out for pressure rise per wavelength. Effects of influential flow parameters have been pointed out through graphs.

Impact of Artificial Compressibility on the Numerical Solution of Incompressible Nanofluid Flow

The numerical solution of compressible flows has become more prevalent than that of incompressible flows.With the help of the artificial compressibility approach,incompressible flows can be solved numerically using the same methods as compressible ones.The artificial compressibility scheme is thus widely used to numerically solve incompressible Navier-Stokes equations.Any numerical method highly depends on its accuracy and speed of convergence.Although the artificial compressibility approach is utilized in several numerical simulations,the effect of the compressibility factor on the accuracy of results and convergence speed has not been investigated for nanofluid flows in previous studies.Therefore,this paper assesses the effect of this factor on the convergence speed and accuracy of results for various types of thermo-flow.To improve the stability and convergence speed of time discretizations,the fifth-order Runge-Kutta method is applied.A computer program has been written in FORTRAN to solve the discretized equations in different Reynolds and Grashof numbers for various grids.The results demonstrate that the artificial compressibility factor has a noticeable effect on the accuracy and convergence rate of the simulation.The optimum artificial compressibility is found to be between 1 and 5.These findings can be utilized to enhance the performance of commercial numerical simulation tools,including ANSYS and COMSOL.

Mathematical Modeling of Biological Fluid Flow Through a Cylindrical Layer with Due Account for Barodiffusion

The work proposes a model of biological fluid flow in a steady mode through a cylindrical layer taking into account convection and diffusion.The model considers finite compressibility and concentration expansion connected with both barodiffusion and additional mechanism of pressure change in the pore volume due to the concentration gradient.Thus,the model is entirely coupled.The paper highlights the complexes composed of scales of physical quantities of different natures.The iteration algorithm for the numerical solution of the problem was developed for the coupled problem.The work involves numerical studies of the considered effects on the characteristics of the flow that can be convective or diffusive,depending on the relation between the dimensionless complexes.It is demonstrated that the distribution of velocity and concentration for different cylinder wall thicknesses is different.It is established that the barodiffusion has a considerable impact on the process in the convective mode or in the case of reduced cylinder wall thickness.

The effect of far field flow on a spherical crystal growth in the undercooled melt

The effect of convective flow on a spherical crystal growth in the undercooled melt with a moderate far field flow is studied. The asymptotic solution of the evolution of the interface of the spherical crystal growth is obtained by the matched asymptotic expansion method. The analytic result shows that the convective flow in the undercooled melt has a strong effect on the evolution of spherical crystal growth. The convective flow induced by the far field flow makes the interface of the growing spherical crystal enhance its growth velocity in the upstream direction of the far field flow and inhibit growth in the downstream direction, and the interface of the decaying spherical crystal further decay in the upstream direction and inhibit decay in the downstream direction. The maximum growth velocity of the interface of the spherical crystal influenced by the far field flow is obtained.

Flow Pattern and Temperature Distribution in HMCZ Silicon Melt

The axial and radial convective flow,temperature fluctuation and distribution in the HMCZ silicon melt are studied tentatively.The experimental results show that the axial and radial convective speeds,the tempera- ture variation and the radial temperature gradient,parallel to magnetic field and near melt surface,all decrease,but the axisymmetry of temperature distribution no longer exists when the magnetic field is applied.

Residence Time Distribution at Laminar Pulsatile Flow in a Straight Pipe

This paper deals with the problem of theoretical identification of the residence time distribution （RTD） characteristics of a straight pipe at laminar pulsatile flow, if tracer diffusion can be neglected. This situation is typical for micro-apparatuses （e.g. fluidic element） and also for flow in large arteries. Residence time distribution based on velocity profiles at pulsatile flow of a Newtonian liquid in a rigid pipe will be derived theoretically and compared with the well known results for a constant flow rate E（τ） = τ-^2/2τ^3 at τ 〉 τ^-/2, where E （τ） is differential distribution, x is residence time and τ^- is the mean residence time. The following part of the paper deals stimulus response experimental techniques using tracers. The principal problem related to laminar and convection dominated pulsatile flows is discussed： Can the impulse response also be identified with the actual residence time distribution in the case of variable flow？ The general answer is no, and differences between RTD and impulse responses are evaluated as a function of the frequency and amplitude of pulsatile flows.

EXPERIMENTAL INVESTIGATION ON THE CHAOTIC PHENOMENA IN THE WAKE OF A NATURAL THERMAL CONVECTION FLOW

Chaotic phenomena in the wake of thermal convection flow fields above a heating flat plate were investigated experimentally. A newly developed electron beam fluorescence technique (EBF) was used to simultaneously measure density fluctuation at 7 points in a cross section above the plate. Correlation dimensions, intermittence coefficients, Fourier spectrum have been obtained for different Grashof numbers. Spatial distribution of correlation dimensions are presented. The experimental result shows that there is a certain relationship between the density fluctuation and the Gr number. And time-spacial characteristic of chaos evolution is also given.

Unsteady natural convection Couette flow of heat generating/absorbing fluid between vertical parallel plates filled with porous material

The extended Brinkman Darcy model for momentum equations and an energy equation is used to calculate the unsteady natural convection Couette flow of a viscous incompressible heat generating/absorbing fluid in a vertical channel （formed by two infinite vertical and parallel plates） filled with the fluid-saturated porous medium. The flow is triggered by the asymmetric heating and the accelerated motion of one of the bounding plates. The governing equations are simplified by the reasonable dimensionless parameters and solved analytically by the Laplace transform techniques to obtain the closed form solutions of the velocity and temperature profiles. Then, the skin friction and the rate of heat transfer are consequently derived. It is noticed that, at different sections within the vertical channel, the fluid flow and the temperature profiles increase with time, which are both higher near the moving plate. In particular, increasing the gap between the plates increases the velocity and the temperature of the fluid, however, reduces the skin friction and the rate of heat transfer.

A simplified two-dimensional boundary element method with arbitrary uniform mean flow

To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.

A Study of Two Dimensional Unsteady MHD Free Convection Flow over a Vertical Plate in the Presence of Radiation

This paper compacts with an exact analysis of radiative effects on the magnetohydrodynamic (MHD) free convection flow of an electrically conducting incompressible viscous fluid over a vertical plate. The non-dimensional continuity, momentum, and energy equations are solved using appropriate transformation. The dimensionless momentum and energy equations are solved numerically through an explicit finite difference method. The stability and convergence analysis also discussed. Finally, outcomes of the parameters on velocity and temperature profiles are displayed graphically and qualitatively.

Three-dimensional mixed convection squeezing flow

The unsteady mixed convection squeezing flow of an incompressible Newtonian fluid between two vertical parallel planes is discussed. The fluid is electrically conducting. The governing equations are transformed into ordinary differential equations （ODEs） by appropriate transformations. The transformed equations are solved successfully by a modern and powerful technique. The effects of the emerging parameters on the flow and heat transfer characteristics are studied and examined. The values of the skin friction coefficient and the local Nusselt number are tabulated and analyzed.

Effects of transpiration on unsteady MHD flow of an upper convected Maxwell (UCM) fluid passing through a stretching surface in the presence of a first order chemical reaction

The aim of this article is to present the effects of transpiration on the unsteady two-dimensional boundary layer flow of non-Newtonian fluid passing through a stretching sheet in the presence of a first order constructive/destructive chemical reaction. The upper-convected Maxwell （UCM） model is used here to characterize the non-Newtonian behavior of the fluid. Using similarity solutions, the governing nonlinear partial differential equations are transformed into ordinary ones and are then solved numerically by the shooting method. The flow fields and mass transfer are significantly influenced by the governing parameters. The fluid velocity initially decreases as the unsteadiness parameter increases and the concentration decreases significantly due to the increase in the unsteadiness. The effect of increasing values of transpiration （suction） and the Maxwell parameter is to suppress the velocity field; however, the concentration is enhanced as transpiration （suction） and the Maxwell parameter increase. Also, it is found that the fluid velocity decreases as the magnetic parameter increases; however, the concentration increases in this case.

Analysis of Combined Radiation and Forced Convection Heat Transfer in 3D Laminar Flow over an Inclined Forward Facing Step

In the current study, a numerical investigation of three-dimensional combined convection-radiation heat transfer over an inclined forward facing step (FFS) in a horizontal rectangular duct is presented. The fluid is treated as a gray, absorbing, emitting and scattering medium. To simulate the incline surface of FFS, the blocked-off method is employed in this study. The set of governing equations for gas flow are solved numerically using the CFD technique to obtain the temperature and velocity fields. Since the gas is considered as a radiating medium, all of the convection, conduction and radiation heat transfer mechanisms are presented in the energy equation. For computation of radiative term in energy equation, the radiative transfer equation (RTE) is solved numerically by the discrete ordinates method (DOM) to find the divergence of radiative heat flux distribution inside the radiating medium. The effects of optical thickness, radiation-conduction parameter and albedo coefficient on heat transfer behavior of the system are carried out.

Flows’ Similarities between Tornadoes or Cyclones Kinematics and Motions Resulting from Weather Phenomena Coupling Geostrophic Wind with Passive Convection

Tornadoes and cyclones, as is stated in numerous literary and audiovisual works dedicated to these out of balance physical systems, are two powerful and spectacular atmospheric phenomena whose vertical and horizontal profiles of winds and temperatures are not yet well known. Indeed, data and routine observations accumulated in the World Meteorological Organization (WMO) databases, regardless of their diversity and perfection of the instruments used to achieve these data (e.g. satellites, onboard cameras, wind profilers, ultra modern calculators, etc.), offer mind-blowing performances on the extent of damage caused by these disturbances, but information provided by these ground and space based observations will never allow access to real profiles of winds associated with tornadoes and cyclones both at the ground’s surface and aloft. The works recently carried out by C. Mbane Biouele allow us to discover that winds associated with tornadoes and hurricanes result from vectors addition of troposphere’s horizontal geostrophic winds and vertical movements associated with passive convection. Unfortunately, geostrophic wind and passive convection are two familiar meteorological phenomena described with much awkwardness and monumental mistakes by all scientific books written by authors who have remained loyal to Hadley principle which states (for centuries) that hot air is lighter than cold air. It is very important to know that C. Mbane Biouele’s very recent publications demonstrate that Hadley principle is not valid in the troposphere’s regions occupied by Ferrell cells. Indeed, it is urgent for the development of meteorology to highlight with great insistence to everyone that there is a Physics principle diametrically opposed to popular Hadley one which provides thermodynamic reasons of the formation of Ferrell cells. This Principle will be named Mbane Biouele Principe and be clearly stated in this paper.

Soret and Dufour Effects on Unsteady Free Convection Fluid Flow in the Presence of Hall Current and Heat Flux

Unsteady MHD natural convective heat and mass transfer flow through a semi-infinite vertical porous plate in a rotating system have been investigated with the combined Soret and Dufour effects in the presence of Hall current and constant heat flux. It is considered that the porous plate is subjected to constant heat flux. The obtained non-dimensional, non-similar coupled non-linear and partial differential equations have been solved by explicit finite difference technique. Numerical solutions for velocities, temperature and concentration distributions are obtained for various parameters by the above mentioned technique. The local and average shear stresses, Nusselt number as well as Sherwood number are also investigated. The stability conditions and convergence criteria of the explicit finite difference scheme are established for finding the restriction of the values of various parameters to get more accuracy. The obtained results are illustrated with the help of graphs to observe the effects of various legitimate parameters.

Three-dimensional mixed convection stagnation-point fow past a vertical surface with second-order slip velocity

This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.