Surface Forces Driven Wrinkling of an Elastic Half-Space Coated with a Thin Stiff Surface Layer

This paper studies surface instability of a coated semi-infinite linear elastic body interacting with another flat rigid body through surface van der Waals (vdW) forces under plane strain conditions. The emphasis is on the effect of the surface coating layer on the wavelength of surface wrinkling. It is shown that the surface of the coated elastic half-plane is always unstable even in the presence of a very stiff coating layer. However, the numerical results show that the stiff coating layer has a significant effect on the wavelength of the surface instability mode and can effectively prevent the surface from short-wavelength wrinkling. In particular, the surface tangential displacement associated with the surface instability vanishes when the elastic half-plane is incompressible. In this case, the in-plane rigidity of the coating layer has no effect on surface instability while the bending stiffness of the coating layer has an effect on the wave-length of the surface instability mode. Furthermore, the Poisson’s ratio of elastic half-plane has a significant role in the surface instability and the associated wave-length. C. Q. Ru is on leave from the University of Alberta, Edmonton, T6G. 2G8, Canada.     

Exact Solutions for a Thick Elastic Plate with a Thin Elastic Surface Layer

A procedure has been developed in previous papers for constructing exact solutions of the equations of linear elasticity in a plate (not necessarily thin) of inhomogeneous isotropic linearly elastic material in which the elastic moduli depend in any specified manner on a coordinate normal to the plane of the plate. The essential idea is that any solution of the classical equations for a hypothetical thin plate or laminate (which are two-dimensional theories) generates, by straightforward substitutions, a solution of the three-dimensional elasticity equations for the inhomogeneous material. In this paper we consider a thick plate of isotropic elastic material with a thin surface layer of different isotropic elastic material. It is shown that the interface tractions and in-plane stress discontinuities are determined only by the initial two-dimensional solution, without recourse to the three-dimensional elasticity theory. Two illustrative examples are described.     

Homogenization of Wall-Slip Gas Flow Through Porous Media

The permeability of reservoir rocks is most commonly measured with an atmospheric gas. Permeability is greater for a gas than for a liquid. The Klinkenberg equation gives a semi-empirical relation between the liquid and gas permeabilities. In this paper, the wall-slip gas flow problem is homogenized. This problem is described by the steady state, low velocity Navier–Stokes equations for a compressible gas with a small Knudsen number. Darcy's law with a permeability tensor equal to that of liquid flow is shown to be valid to the lowest order. The lowest order wall-slip correction is a local tensorial form of the Klinkenberg equation. The Klinkenberg permeability is a positive tensor. It is in general not symmetric, but may under some conditions, which we specify, be symmetric. Our result reduces to the Klinkenberg equation for constant viscosity gas flow in isotropic media.     

T-RANS Simulations of Subcritical Flow with Heat Transfer Past a Circular Cylinder Surrounded by a Thin Porous Layer

We study flow and heat transfer to a cylinder in cross flow at Re = 3,900–80,000 by means of three-dimensional transient RANS (T-RANS) simulations, employing an RNG k − ε turbulence model. Both the case of a bare solid cylinder and that of a solid cylinder surrounded at some fixed distance by a thin porous layer have been studied. The latter configuration is a standard test geometry for measuring the insulating and protective performance of garments. In this geometry, the flow in the space between the solid cylinder and the porous layer is laminar but periodic, whereas the outer flow is transitional and characterized by vortex shedding in the wake of the cylinder. The results from the T-RANS simulations are validated against data from Direct Numerical Simulations and experiments. It is found that T-RANS is very well suited for simulating this type of flow. The transient nature of the flow underneath the porous layer is well reproduced, as well as the influence of vortex shedding on the heat transfer in the downstream stagnation zone. T-RANS results are found to be in much better agreement with DNS and experimental data than results from steady-state RANS.     

Free Convection Boundary Layer Flow Past a Vertical Surface in a Porous Medium with Temperature-Dependent Properties

Numerical solutions for the free convection heat transfer in a viscous fluid at a permeable surface embedded in a saturated porous medium, in the presence of viscous dissipation with temperature-dependent variable fluid properties, are obtained. The governing equations for the problem are derived using the Darcy model and the Boussinesq approximation (with nonlinear density temperature variation in the buoyancy force term). The coupled non-linearities arising from the temperature-dependent density, viscosity, thermal conductivity, and viscous dissipation are included. The partial differential equations of the model are reduced to ordinary differential equations by a similarity transformation and the resulting coupled, nonlinear ordinary differential equations are solved numerically by a second order finite difference scheme for several sets of values of the parameters. Also, asymptotic results are obtained for large values of | f _w|. Moreover, the numerical results for the velocity, the temperature, and the wall-temperature gradient are presented through graphs and tables, and are discussed. It is observed that by increasing the fluid variable viscosity parameter, one could reduce the velocity and thermal boundary layer thickness. However, quite the opposite is true with the non-linear density temperature variation parameter.     

Deformation of the Free Surface in a Moving Locally-Heated Thin Liquid Layer

Liquid film flow on a vertical surface is studied experimentally and theoretically under the determining influence of the thermocapillary forces. In the two-dimensional steady-state case the shape of the film surface is calculated numerically within the thin layer approximation with allowance for the temperature dependence of the viscosity of the liquid and redistribution of the heat flux in the heating element. A local heat source was used in the experiments to produce temperature gradients up to 10 K/mm on the liquid surface. The film thickness was determined by means of the schlieren method with reflection. The relative thickness of the roller in the upper heater edge zone, characteristic of the formation of regular structures, is measured. The thickness is h/h ₀=1.32 ±0.07, which agrees satisfactorily with the results of numerical calculations.     

Brinkmann Model and Double Penalization Method for the Flow Around a Porous Thin Layer

In this paper we study a penalization method used to compute the flow of a viscous fluid around a thin layer of porous material. Using a BKW method, we perform an asymptotic expansion of the solution when a little parameter, measuring the thickness of the thin layer and the inverse of the penalization coefficient, tends to zero. We compare then this numerical method with a Brinkman model for the flow around a porous thin layer.      

Stability of the Boundary Layer with Internal Heat Release and Gas Injection through a Porous Wall

Fluid Dynamics - The problem of stability of a subsonic boundary layer is solved under the conditions of heat supply inside the boundary layer with injection of a homogeneous gas through a porous...     

Onset of Buoyancy-Driven Convection in a Liquid-Saturated Cylindrical Anisotropic Porous Layer Supported by a Gas Phase

A theoretical analysis of convective instability driven by buoyancy forces under the transient concentration fields is conducted in an initially quiescent, liquid-saturated, and anisotropic cylindrical porous layer supported by a gas phase. Darcy’s law and Boussinesq approximation are used to explain the characteristics of fluid motion, and linear stability theory is employed to predict the onset of buoyancy-driven motion. Under the quasi-steady-state approximation, the stability equations are derived in a similar boundary layer coordinate and solved by the numerical shooting method. The critical $Ra_D$ is determined as a function of the anisotropy ratio. Also, the onset time and corresponding wavelength are obtained for the various anisotropic ratios. The onset time becomes smaller with increasing $Ra_D$ and follows the asymptotic relation derived in the infinite horizontal porous layer. Anisotropy effect makes the system more stable by suppressing the vertical velocity.     

Gas Capture in a Cavity with Porous Walls

A configuration like an upside-down bell made of porous material is considered which is initially dry but then subjected to a rising pool of liquid. As liquid touches the rim of the bell, capillary transport is initiated. Starting with a vertical wicking phase, the imbibing liquid will eventually reach the ceiling of the bell and switch over to horizonal wicking. At the end of the horizontal wicking, the cavity inside the porous bell is enclosed by liquid and the gas inside it is captured. We present a model to describe the capillary transport in the bell for both Cartesian and cylindrical geometry. As far as possible, we derive analytical solutions to the normalized differential equations that describe the problem. Beyond analytical solutions, we use Runge–Kutta shooting method to obtain numerical results. We calculate the normalized closure time to capture the gas, the amount of captured gas, and reflect on the pressure development in the gas chamber.     

Peristaltic Flow of a Maxwell Model Through Porous Boundaries in a Porous Medium

In the present investigation we have presented the peristaltic flow of a linear Maxwell model through porous boundaries in a porous medium. The governing non-dimensional partial differential are solved in wave frame by using regular perturbation method and assumed form of solution. We have discussed the problem only for free pumping case. The effects of various physical parameters involved in the problem have been investigated and shown graphically.     

Nonlinear Stability of Convection in a Porous Layer with Solid Partitions

We show that for many classes of convection problem involving a porous layer, or layers, interleaved with finite but non-deformable solid layers, the global nonlinear stability threshold is exactly the same as the linear instability one. The layer(s) of porous material may be of Darcy type, Brinkman type, possess an anisotropic permeability, or even be such that they are of local thermal non-equilibrium type where the fluid and solid matrix constituting the porous material may have different temperatures. The key to the global stability result lies in proving the linear operator attached to the convection problem is a symmetric operator while the nonlinear terms must satisfy appropriate conditions.     

Free Fluid Flow Through a Visco-Deformable Porous Medium with a Moving Boundary

One-dimensional Darcy-law flow through a porous matrix representing a high-viscosity liquid is investigated. The flow develops in a region which depends on time due to sedimentation. The problem considered simulates the geological process of sedimentation in a basin. In accordance with geological data, the permeability and viscosity coefficients of the matrix are assumed to depend nonlinearly on the porosity. The asymptotic properties of the flow are described for large times. The agreement between the results of asymptotic and numerical solutions is satisfactory at intermediate times and good at large times under the realistic sedimentary basin conditions. The simplicity of the asymptotic solution obtained makes it possible to vary the problem parameters and determine the porosity, pressure, and velocities for particular geological conditions by means of simple calculations.     

Thermal Instability in a Porous Medium Layer Saturated with a Viscoelastic Nanofluid

The onset of convection in a horizontal layer of a porous medium saturated with a viscoelastic nanofluid was studied in this article. The modified Darcy model was applied to simulate the momentum equation in porous media. An Oldroyd-B type constitutive equation was used to describe the rheological behavior of viscoelastic nanofluids. The model used for the viscoelastic nanofluid incorporates the effects of Brownian motion and thermophoresis. The onset criterion for stationary and oscillatory convection was analytically derived. The effects of the concentration Rayleigh number, Prandtl number, Lewis number, capacity ratio, relaxation, and retardation parameters on the stability of the system were investigated. Oscillatory instability is possible in both bottom- and top-heavy nanoparticle distributions. Results indicated that there is competition among the processes of thermophoresis, Brownian diffusion, and viscoelasticity that causes the convection to set in through oscillatory rather than stationary modes. Regimes of stationary and oscillatory convection for various parameters were derived and are discussed in detail.     

Mixed Convection in a Ventilated Cavity Filled with a Triangular Porous Layer

Accurate prediction of the macroscopic flow parameters needed to describe flow in porous media relies on a good knowledge of flow field distribution at a much smaller scale—in the pore spaces. The extent of the inertial effect in the pore spaces cannot be underestimated yet is often ignored in large-scale simulations of fluid flow. We present a multiscale method for solving Oseen’s approximation of incompressible flow in the pore spaces amid non-periodic grain patterns. The method is based on the multiscale finite element method [MsFEM Hou and Wu in J Comput Phys 134:169–189, 1997)] and is built in the vein of Crouzeix and Raviart elements (Crouzeix and Raviart in Math Model Numer Anal 7:33–75, 1973). Simulations of inertial flow in highly non-periodic settings are conducted and presented. Convergence studies in terms of numerical errors relative to the reference solution are given to demonstrate the accuracy of our method. The weakly enforced continuity across coarse element edges is shown to maintain accurate solutions in the vicinity of the grains without the need for any oversampling methods. The penalisation method is employed to allow a complicated grain pattern to be modelled using a simple Cartesian mesh. This work is a stepping stone towards solving the more complicated Navier–Stokes equations with a nonlinear inertial term.     

Mixed Convection Boundary Layer Flow Near the Stagnation Point on a Vertical Surface Embedded in a Porous Medium with Anisotropy Effect

The steady boundary-layer flow near the stagnation point on a vertical flat plate embedded in a fluid-saturated porous medium characterized by an anisotropic permeability is investigated. Using appropriate similarity transformation, the governing system of partial differential equations is transformed into a system of ordinary differential equations. This system is then solved numerically. The features of the flow and the heat transfer characteristics for different values of the governing parameters, namely, the modified mixed convection parameter Λ, and the anisotropy parameter A are analyzed and discussed. It is found that dual solutions exist for both assisting and opposing flows. Moreover, the range of Λ for which the solution exists increases with A.     

Weakly Nonlinear Convection in a Porous Layer with Multiple Horizontal Partitions

We consider convection in a horizontally uniform fluid-saturated porous layer which is heated from below and which is split into a number of identical sublayers by impermeable and infinitesimally thin horizontal partitions. Rees and Genç (Int J Heat Mass Transfer 54:3081–3089, 2010) determined the onset criterion by means of a detailed analytical and numerical study of the corresponding dispersion relation and showed that this layered system behaves like the single-sublayer constant-heat-flux Darcy–Bénard problem when the number of sublayers becomes large. The aim of the present work is to use a weakly nonlinear analysis to determine whether the layered system also shares the property of the single-sublayer constant-heat-flux Darcy–Bénard problem by having square cells, as opposed to rolls, as the preferred planform for convection.     

Mixed Convection Boundary Layer Flow on a Vertical Surface in a Saturated Porous Medium: New Perturbation Solutions

Transport in Porous Media -     

Flow and Heat Transfer Through a Polygonal Duct filled with a Porous Medium

The fully developed flow and H1 heat transfer in a polygonal duct filled with a Darcy–Brinkman medium is studied. The efficient method of boundary collocation is used. The problem is governed by the duct shape and a non-dimensional parameter s which characterizes the inverse square root of permeability. Asymptotic formulas for small and large s are derived.