Cooling of a cylindrical container blasted into the rock

An upright cylindrical container blasted into the rock is filled instantaneously with a warm liquid. Heat is transferred from the liquid into the surrounding rock and the open air. The temperature of the liquid and the surrounding rock is determined as a function of time.The differential equation and the auxiliary conditions of the transient heat conduction problem are Laplace transformed, the subsidiary equations are solved by two-dimensional relaxation, and the resulting temperature is obtained by means of numerical inversion of the Laplace transform.The results are presented numerically and graphically.     

Propagation of shear waves in a poroelastic layer constrained between two elastic layers

In this paper, we investigated the propagation of shear waves in a transversely isotropic poroelastic layer constrained between two elastic layers. Following Biot’s theory, the dispersion equation for shear waves in this structure was derived. The numerical values on the dimensionless phase velocities are calculated and presented graphically to illustrate the dependences upon geometry, anisotropy and porosity comparatively. It is observed that the phase velocities increase with the increase of the porosity and the decrease of the anisotropy. In addition, the geometry in this structure has a significant effect on the phase velocity of the shear waves.     

Error analysis of a Collocation method for numerically inverting a Laplace transform in case of real samples

In [S. Cuomo, L. D’Amore, A. Murli, M.R. Rizzardi, Computation of the inverse Laplace transform based on a collocation method which uses only real values, J. Comput. Appl. Math., 198 (1) (2007) 98–115] the authors proposed a Collocation method (C-method) for real inversion of Laplace transforms (Lt), based on the truncated Laguerre expansion of the inverse function:

where σ, b are parameters and c_k, kN, are the MacLaurin coefficients of a function depending on the Lt. The computational kernel of a C-method is the solution of a Vandermonde linear system, where the right hand side is obtained evaluating the Lt on the real axis. The Bjorck Pereira algorithm has been used for solving the Vandermonde linear system, providing a computable componentwise error bound on the solution.

For an inversion problem on discrete data F is known on a pre-assigned set of points (we refer to these points as samples of F) only and the major challenge is to deal with a significative loss of information. A natural approach to overcome this intrinsic difficulty is to construct a suitable fitting model that approximates the given data. In this case, we show that such approach leads to a C-method with perturbed right hand side, and then we use again the Bjorck Pereira algorithm.

Starting from the error introduced by the fitting model, we study its propagation in order to determine the maximum attainable accuracy on f_N. Moreover we derive a computable error bound that allows to get the suitable value of the parameter N that gives the maximum attainable accuracy.     

Love waves in a fluid-saturated porous layer under a rigid boundary and lying over an elastic half-space under gravity

In this paper, mathematical modeling of the propagation of Love waves in a fluid-saturated porous layer under a rigid boundary and lying over an elastic half-space under gravity has been considered. The equations of motion have been formulated separately for different media under suitable boundary conditions at the interface of porous layer, elastic half-space under gravity and rigid layer. Following Biot, the frequency equation has been derived which contain Whittaker’s function and its derivative that have been expanded asymptotically up to second term (for approximate result) for large argument due to small values of Biot’s gravity parameter (varying from 0 to 1). The effect of porosity and gravity of the layers in the propagation of Love waves has been studied. The effect of hydrostatic initial stress generated due to gravity in the half-space has also been shown in the phase velocity of Love waves. The phase velocity of Love waves for first two modes has been presented graphically. Frequency equations have also been derived for some particular cases, which are in perfect agreement with standard results. Subsequently the lower and upper bounds of Love wave speed have also been discussed.     

A common error made in investigation of boundary layer flows

The concept of boundary layer flow, introduced in 1904 by Prandtl, is a popular field in Fluid Mechanics for engineers, physicists and mathematicians. The present work is a critique to many papers published in the last 15 years in the field of boundary layer flow. The critique concerns the shape of velocity, temperature and concentration profiles which are truncated due to small calculation domain used during the numerical solution procedure. These truncated profiles are not compatible with the boundary layer theory and introduce errors in wall shear stress and wall heat transfer values.     

Steady flow of a viscous liquid in a porous elliptic tube

Flow of a viscous-liquid in a porous tube of elliptic cross-section is studied using the generalized momentum equation. As a particular case, flow of the liquid in a tube of a circular cross-section is obtained. It is observed that the classical Darcian effect is realized only in a core very near to the axis of the tube while the non-Darcian phenomenon is felt predominantly near the boundary of the tube.     

Three-dimensional flow over a stretching surface in a viscoelastic fluid

This article looks at the hydrodynamic elastico-viscous fluid over a stretching surface. The equations governing the flow are reduced to ordinary differential equations, which are analytically solved by applying an efficient technique namely the homotopy analysis method (HAM). The solutions for the velocity components are computed. The numerical values of wall skin friction coefficients are also tabulated. The present HAM solution is compared with the known exact solution for the two-dimensional flow and an excellent agreement is found.     

On boundary layer flow of a Sisko fluid over a stretching sheet

Abstract

In this paper, the steady boundary layer flow of a non-Newtonian fluid over a nonlinear stretching sheet is investigated. The Sisko fluid model, which is combination of power-law and Newtonian fluids in which the fluid may exhibit shear thinning/thickening behaviors, is considered. The boundary layer equations are derived for the two-dimensional flow of an incompressible Sisko fluid. Similarity transformations are used to reduce the governing nonlinear equations and then solved analytically using the homotopy analysis method. In addition, closed form exact analytical solutions are provided for n = 0 and n = 1. Effects of the pertinent parameters on the boundary layer flow are shown and solutions are contrasted with the power-law fluid solutions.     

Time-harmonic acoustic propagation in the presence of a shear flow

This work deals with the numerical simulation, by means of a finite element method, of the time-harmonic propagation of acoustic waves in a moving fluid, using the Galbrun equation instead of the classical linearized Euler equations. This work extends a previous study in the case of a uniform flow to the case of a shear flow. The additional difficulty comes from the interaction between the propagation of acoustic waves and the convection of vortices by the fluid. We have developed a numerical method based on the regularization of the equation which takes these two phenomena into account. Since it leads to a partially full matrix, we use an iterative algorithm to solve the linear system.     

The velocity and shear stress functions of the Falkner-Skan equation arising in boundary layer theory

New compactness results on the velocity functions and shear stress functions of the well-known Falkner-Skan equation are obtained. The methodology is to utilize the equivalence between the Falkner-Skan equation and a singular integral equation established recently by Lan and Yang.     

Derivation of a contact law between a free fluid and thin porous layers via asymptotic analysis methods

We describe the asymptotic behaviour of an incompressible viscous free fluid in contact with a porous layer flow through the porous layer surface. This porous layer has a small thickness and consists of thin channels periodically distributed. Two scales are present in this porous medium, one associated to the periodicity of the distribution of the channels and the other to the size of these channels. Proving estimates on the solution of this Stokes problem, we establish a critical link between these two scales. We prove that the limit problem is a Stokes flow in the free domain with further boundary conditions on the basis of the domain which involve an extra velocity, an extra pressure and two second-order tensors. This limit problem is obtained using Γ-convergence methods. We finally consider the case of a Navier–Stokes flow within this context.     

A nonlinear two phase fluid flow through a porous medium in presence of a well

We study a flow of fresh and salt water in a two dimensional axially symmetric coastal aquifer with a well on the central axis. The flow is governed by a nonlinear Darcy's law. We also show the behaviour of the solution when the out flow of salt water at well goes to 0. Received May 1999     

N-dimensional fractional Fokker-Planck equation and its solutions for anomalous radial two-phase flow in porous media

In this paper we present a time fractional Fokker-Planck equation (fFPE) for radial two-phase flow of liquid and gas in porous media. The fFPE of order α is solved for both two- and three-dimensional flow patterns using the Laplace transform method. The general solutions of the fFPE for both two- and three- dimensional flows are given as a convolution integral of the input and a kernel in the Laplace domain. Special solutions for a large value and a periodic boundary condition are also given in the time domain when the inverse Laplace transform can be found analytically. The fFPE for two-phase flow in porous media presented in this paper is the first report of its kind.     

Approximate analytical solutions of a class of boundary layer equations over nonlinear stretching surface

Third order nonlinear ordinary differential equations, subject to appropriate boundary conditions arising in fluid dynamics, are solved using three different methods viz., the Dirichlet series, method of stretching of variables, and asymptotic function method. Similarity transformations are used to convert the governing partial differential equations into nonlinear ordinary differential equations. The numerical results obtained from the above methods for various problems are given in terms of skin friction. Our study revealed that the results obtained from these methods agree well with those of direct numerical simulation of ordinary differential equations. Also, these methods have advantages over pure numerical methods in obtaining derived quantities such as velocity profile accurately for various values of the parameters at a stretch.     

Thermo-elastodynamic response of a spherical cavity in saturated poroelastic medium

In this paper, an attempt has been made to investigate the thermo-hydro-elastodynamic response of a spherical cavity in isotropic saturated poroelastic medium when subjected to a time dependent thermal/mechanical source. The fully coupling thermo-hydro-elastodynamic model is presented on the basis of equations of motion, fluid flow, feat flow and constitutive equation with effective stress and temperature change. Solutions of displacement, temperature and stresses are obtained by using a semi-analytical approach in the domain of Laplace transform. Numerical results are also performed for portraying the nature of variations of the field variables, i.e. the coefficient of thermo-osmosis, the permeability. In addition, comparisons are presented with the corresponding partially thermo-hydro-elastodynamic model and thermo-elastodynamic model to ascertain the validity and the difference between these models.     

Convergence analysis of a DDFV scheme for a system describing miscible fluid flows in porous media

In this article, we prove the convergence of a discrete duality finite volume scheme for a system of partial differential equations describing miscible displacement in porous media. This system is made of two coupled equations: an anisotropic diffusion equation on the pressure and a convection‐diffusion‐dispersion equation on the concentration. We first establish some a priori estimates satisfied by the sequences of approximate solutions. Then, it yields the compactness of these sequences. Passing to the limit in the numerical scheme, we finally obtain that the limit of the sequence of approximate solutions is a weak solution to the problem under study. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 723–760, 2015     

An analytical study of linear and nonlinear double diffusive convection in a fluid saturated anisotropic porous layer with Soret effect

The double diffusive convection in a horizontal anisotropic porous layer saturated with a Boussinesq fluid, which is heated and salted from below in the presence of Soret coefficient is studied analytically using both linear and nonlinear stability analyses. The normal mode technique is used in the linear stability analysis while a weak nonlinear analysis based on a minimal representation of double Fourier series method is used in the nonlinear analysis. The generalized Darcy model including the time derivative term is employed for the momentum equation. The critical Rayleigh number, wavenumber for stationary and oscillatory modes and frequency of oscillations are obtained analytically using linear theory. The effect of anisotropy parameters, solute Rayleigh number, Soret parameter and Lewis number on the stationary, oscillatory, finite amplitude convection and heat and mass transfer are shown graphically.     

On over-reflection of acoustic-gravity waves incident upon a magnetic shear layer in a compressible fluid

A study is made of over-reflection of acoustic-gravity waves incident upon a magnetic shear layer in an isothermal compressible electrically conducting fluid in the presence of an external magnetic field. The reflection and transmission coefficients of hydromagnetic acoustic-gravity waves incident upon magnetic shear layer are calculated. The invariance of wave-action flux is used to investigate the properties of reflection, transmission and absorption of the waves incident upon the shear layer, and then to discuss how these properties depend on the wavelength, length scale of the shear layers, and the ratio of the flow speed and phase speed of the waves. Special attention is given to the relationship between the wave-amplification and critical-level behaviour. It is shown that there exists a critical level within the shear layer and the wave incident upon the shear layer is over-reflected, that is, more energy is reflected back towards the source than was originally emitted. The mechanism of the over-reflection (or wave amplification) is due to the fact that the excess reflected energy is extracted by the wave from the external magnetic field. It is also found that the absence of critical level within the shear layer leads to non-amplification of waves. For the case of very large vertical wavelength of waves, the coefficients of incident, reflected and transmitted energy are calculated. In this limiting situation, the wave is neither amplified nor absorbed by the shear layer. Finally, it is shown that resonance occurs at a particular value of the phase velocity of the wave.     

Propagation of torsional surface waves in a homogeneous layer of finite thickness over an initially stressed heterogeneous half-space

In the present paper, the dispersion equation which determines the velocity of torsional surface waves in a homogeneous layer of finite thickness over an initially stressed heterogeneous half-space has been obtained. The dispersion equation obtained is in agreement with the classical result of Love wave when the initial stresses and inhomogeneity parameters are neglected. Numerical results analyzing the dispersion equation are discussed and presented graphically. The result shows that the initial stresses have a pronounced influence on the propagation of torsional surface waves. It has also been shown that the effect of density, directional rigidities and non-homogeneity parameter on the propagation of torsional surface waves is prominent.