Diffusion in poro‐plastic media

A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial‐boundary‐value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto‐visco‐plastic type. The variational form of this problem in Hilbert space is a non‐linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi‐static momentum equation. The essential sufficient conditions for the well‐posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd.     

Generalized Darcy–Oseen resolvent problem

In this paper, we study the well‐posedness of a coupled Darcy–Oseen resolvent problem, describing the fluid flow between free‐fluid domains and porous media separated by a semipermeable membrane. The influence of osmotic effects, induced by the presence of a semipermeable membrane, on the flow velocity is reflected in the transmission conditions on the surface between the free‐fluid domain and the porous medium. To prove the existence of a weak solution of the generalized Darcy–Oseen resolvent system, we consider two auxiliary problems: a mixed Navier–Dirichlet problem for the generalized Oseen resolvent system and Robin problem for an elliptic equation related to the general Darcy equations. © 2016 The Authors Mathematical Methods in the Applied Sciences Published by John Wiley & Sons Ltd.     

First‐Order system least squares for generalized‐Newtonian coupled Stokes‐Darcy flow

The coupled problem for a generalized Newtonian Stokes flow in one domain and a generalized Newtonian Darcy flow in a porous medium is studied in this work. Both flows are treated as a first‐order system in a stress‐velocity formulation for the Stokes problem and a volumetric flux‐hydraulic potential formulation for the Darcy problem. The coupling along an interface is done using the well‐known Beavers–Joseph–Saffman interface condition. A least squares finite element method is used for the numerical approximation of the solution. It is shown that under some assumptions on the viscosity the error is bounded from above and below by the least squares functional. An adaptive refinement strategy is examined in several numerical examples where boundary singularities are present. Due to the nonlinearity of the problem a Gauss–Newton method is used to iteratively solve the problem. It is shown that the linear variational problems arising in the Gauss–Newton method are well posed. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1150–1173, 2015     

Evolution of microstructure in unstable porous media flow: A relaxational approach

We study the flow of two immiscible fluids of different density and mobility in a porous medium. If the heavier phase lies above the lighter one, the interface is observed to be unstable. The two phases start to mix on a mesoscopic scale and the mixing zone grows in time—an example of evolution of microstructure. A simple set of assumptions on the physics of this two‐phase flow in a porous medium leads to a mathematically ill‐posed problem—when used to establish a continuum free boundary problem. We propose and motivate a relaxation of this “nonconvex” constraint of a phase distribution with a sharp interface on a macroscopic scale. We prove that this approach leads to a mathematically well‐posed problem that predicts shape and evolution of the mixing profile as a function of the density difference and mobility quotient. © 1999 John Wiley & Sons, Inc.     

On the well‐posedness of a mathematical model of quorum‐sensing in patchy biofilm communities

We analyze a system of reaction–diffusion equations that models quorum‐sensing in a growing biofilm. The model comprises two nonlinear diffusion effects: a porous medium‐type degeneracy and super diffusion. We prove the well‐posedness of the model. In particular, we present for the first time a uniqueness result for this type of problem. Moreover, we illustrate the behavior of model solutions in numerical simulations. Copyright © 2011 John Wiley & Sons, Ltd.     

About modelling chemical vapour infiltration of pyrocarbon

Chemical vapour infiltration (CVI) is an important method for producing carbon reinforced carbon fibres (CFC). Thereby, initially gaseous carbon is deposited on the surface of a porous substrate. Mathematically, one has to deal with a moving boundary problem formed by the interface between the gas phase and the substrate surface. Within the gas phase, a nonlinear convection‐diffusion‐reaction‐system (cdr‐system) with a reduced reaction scheme to model the chemical reactions has to be solved. One‐dimensional simulations of deposition profiles within cylindrical model pores including an explicit construction of the moving boundary are performed for different values of the process parameters. Based on these calculations geometries and conditions for complete infiltration of the pores can be identified.     

Hydromagnetic rotating flow of a fourth‐order fluid past a porous plate

The rotating flow in the presence of a magnetic field is a problem belonging to hydromagnetics and deserves to be more widely studied than it has been to date. In the non‐linear regime the literature is scarce. We develop the governing equations for the unsteady hydromagnetic rotating flow of a fourth‐order fluid past a porous plate. The steady flow is governed by a boundary value problem in which the order of differential equations is more than the number of available boundary conditions. It is shown that by augmenting the boundary conditions based on asymptotic structures at infinity it is possible to obtain numerical solutions of the nonlinear hydromagnetic equations. Effects of uniform suction or blowing past the porous plate, exerted magnetic field and rotation on the flow phenomena, especially on the boundary layer structure near the plate, are numerically analysed and discussed. The flow behaviours of the Newtonian fluid and second‐, third‐ and fourth‐order non‐Newtonian fluids are compared for the special flow problem, respectively. Copyright © 2004 John Wiley & Sons, Ltd.     

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 existence result for nonisothermal immiscible incompressible 2‐phase flow in heterogeneous porous media

In the present article, we study the temperature effects on two‐phase immiscible incompressible flow through a porous medium. The mathematical model is given by a coupled system of 2‐phase flow equations and an energy balance equation. The model consists of the usual equations derived from the mass conservation of both fluids along with the Darcy‐Muskat and the capillary pressure laws. The problem is written in terms of the phase formulation; ie, the saturation of one phase, the pressure of the second phase, and the temperature are primary unknowns. The major difficulties related to this model are in the nonlinear degenerate structure of the equations, as well as in the coupling in the system. Under some realistic assumptions on the data, we show the existence of weak solutions with the help of an appropriate regularization and a time discretization. We use suitable test functions to obtain a priori estimates. We prove a new compactness result to pass to the limit in nonlinear terms.     

Global existence,singular solutions,and ill‐posedness for the Muskat problem

The Muskat, or Muskat‐Leibenzon, problem describes the evolution of the interface between two immiscible fluids in a porous medium or Hele‐Shaw cell under applied pressure gradients or fluid injection/extraction. In contrast to the Hele‐Shaw problem (the one‐phase version of the Muskat problem), there are few nontrivial exact solutions or analytic results for the Muskat problem. For the stable, forward Muskat problem, in which the higher‐viscosity fluid expands into the lower‐viscosity fluid, we show global‐in‐time existence for initial data that is a small perturbation of a flat interface. The initial data in this result may contain weak (e.g., curvature) singularities. For the unstable, backward problem, in which the higher‐viscosity fluid contracts, we construct singular solutions that start off with smooth initial data but develop a point of infinite curvature at finite time. © 2004 Wiley Periodicals, Inc.     

Continuous Dependence on the Heat Source in Resonant Porous Penetrative Convection

I study the structural stability for a problem in a porous medium when the density of saturating liquid is a nonlinear function of temperature and an internal heat source is present. It has been shown that for this problem when one considers thermal convection in a plane infinite layer then resonance may occur between internal layers that arise. A key parameter is the internal heat source and its presence may lead to oscillatory instability inducing resonance. Therefore, in this paper, I analyze the general structural stability problem of continuous dependence on the heat source itself for a model of nonisothermal flow in a porous medium of Forchheimer type, in a general three‐dimensional domain.     

A second‐order backward difference time‐stepping scheme for penalized Navier–Stokes equations modeling filtration through porous media

We investigate the stability and convergence of a fully implicit, linearly extrapolated second‐order backward difference time‐stepping scheme for the penalized Navier–Stokes equations modeling filtration through porous media. In the penalization approach, an extended Navier–Stokes equation is used in the entire computational domain with suitable resistance terms to mimic the presence of porous medium. It is widely used as an alternative to the heterogeneous approach in which different types of partial differential equations (PDEs) are used in fluid and porous subregions along with suitable continuity conditions at the interface. However, the introduction of extra resistance terms makes the penalized Navier–Stokes equations more nonlinear. We prove that the linearly extrapolated scheme is unconditionally stable and derive optimal order error estimates without any stability condition. To show feasibility and applicability of the approach, it is used to numerically solve a passive control problem in which flow around a solid body is controlled by adding porous layers on the surface. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 681–705, 2016     

MODELING THE EFFECTS OF URBAN HEAT ISLANDS GENERATED MESOSCALE WIND ON AIR POLLUTION DISPERSION IN A PATCHY ATMOSPHERE

Abstract The paper presents a steady‐state two‐dimensional mathematical model to study the dispersion of air pollutants continuously emitting from a point source in a patchy atmosphere. The first patch, which extends to some downwind distance from the source (present in first patch), is the urban center, whereas the second patch is its adjacent countryside. For simulating the effects of urban heat islands on pollutant distribution, a local mesoscale wind is introduced in the first patch along with the usual large‐scale wind. Keeping in mind the real‐life situations, the meteorological, deposition, and removal parameters are taken to be different in different patches. The concentration of pollutants in both the patches has been computed under different stability conditions of the atmosphere through finite difference scheme. The results of the model show that the urban heat island (mesoscale wind) has a tremendous impact on air pollution distribution, both over the urban center as well as its adjacent countryside (rural). The model solutions also reveal about the concentration distribution of pollutants in different parts of the atmosphere in different atmospheric conditions.     

Flexural‐Gravity Wave Scattering by a Circular‐Arc‐Shaped Porous Plate

The interaction of flexural‐gravity waves with a thin circular‐arc‐shaped permeable plate submerged beneath the ice‐covered surface of water with uniform finite depth is considered under the assumption of linear theory. The problem is reduced to a second kind hypersingular integral equation for the potential difference across the plate which is solved approximately by an expansion–collocation method. Utilizing the solution, the reflection and the transmission coefficients and the hydrodynamic forces are evaluated numerically. The focus of the paper is to illustrate the effect of a porous curved plate submerged in finite depth water with an ice‐cover on the normally incident waves. Numerical results for a circular‐arc‐shaped plate for different configurations are derived and represented graphically. Also, by choosing an appropriate set of parameters, the known results for a circular‐arc‐shaped rigid plate submerged in deep water and a semicircular porous plate submerged in finite depth water with a free surface are recovered as special cases.     

On some local property of the shape of interfaces for the porous media equation

We consider the porous media equation with absorption for various conditions and prove that the shape of ist interface never becomes strongly upward convex. For this sake we derive an improperly posed estimate for solutions of the porous media equation for the non‐characteristic Cauchy problem (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

多孔介质中两相不可压缩可混溶驱动问题的动态混合元方法及其特征修正格式

许多依赖时间的问题涉及到局部化现象,如突出的前沿位置、激波、边界层等, 其位置随时间而变动.多孔介质中两相不可压缩可混溶驱动问题是一典型的、有代表性 的"局部化现象"问题,其数学模型为耦合非线性偏微分方程组的初边值问题.为减轻数 值解在局部前沿位置的数值振荡,提高解的精确性,本文给出了该问题的动态混合元格 式和沿特征线修正的动态混合元格式,证明了其收敛性,并给出了误差估计.     

Block‐centered upwind multistep difference method and convergence analysis for numerical simulation of oil reservoir

The three‐dimensional displacement of two‐phase flow in porous media is a preliminary problem of numerical simulation of energy science and mathematics. The mathematical model is formulated by a nonlinear system of partial differential equations to describe incompressible miscible case. The pressure is defined by an elliptic equation, and the concentration is defined by a convection‐dominated diffusion equation. The pressure generates Darcy velocity and controls the dynamic change of concentration. We adopt a conservative block‐centered scheme to approximate the pressure and Darcy velocity, and the accuracy of Darcy velocity is improved one order. We use a block‐centered upwind multistep method to solve the concentration, where the time derivative is approximated by multistep method, and the diffusion term and convection term are treated by a block‐centered scheme and an upwind scheme, respectively. The composite algorithm is effective to solve such a convection‐dominated problem, since numerical oscillation and dispersion are avoided and computational accuracy is improved. Block‐centered method is conservative, and the concentration and the adjoint function are computed simultaneously. This physical nature is important in numerical simulation of seepage fluid. Using the convergence theory and techniques of priori estimates, we derive optimal estimate error. Numerical experiments and data show the support and consistency of theoretical result. The argument in the present paper shows a powerful tool to solve the well‐known model problem.     

Numerical analysis of a finite volume scheme for a seawater intrusion model with cross‐diffusion in an unconfined aquifer

We consider a degenerate parabolic system modeling the flow of fresh and saltwater in a porous medium in the context of seawater intrusion. We propose and analyze a finite volume scheme based on two‐point flux approximation with upwind mobilities. The scheme preserves at the discrete level the main features of the continuous problem, namely the nonnegativity of the solutions, the decay of the energy and the control of the entropy and its dissipation. Based on these nonlinear stability results, we show that the scheme converges toward a weak solution to the problem. Numerical results are provided to illustrate the behavior of the model and of the scheme.     

Monotone iterative method of upper and lower solutions applied to a multilayer combustion model in porous media

This paper presents a new nonlinear reaction–diffusion–convection system coupled with a system of ordinary differential equations that models a combustion front in a multilayer porous medium. The model includes heat transfer between the layers and heat loss to the external environment. A few assumptions are made to simplify the model, such as incompressibility; then, the unknowns are determined to be the temperature and fuel concentration in each layer. When the fuel concentration in each layer is a known function, we prove the existence and uniqueness of a classical solution for the initial and boundary value problem for the corresponding system. The proof uses a new approach for combustion problems in porous media. We construct monotone iterations of upper and lower solutions and prove that these iterations converge to a unique solution for the problem, first locally and then, in time, globally.     

Global solutions and blow‐up problems to a porous medium system with nonlocal sources and nonlocal boundary conditions

This paper deals with a porous medium system with nonlocal sources and weighted nonlocal boundary conditions. The main aim of this paper is to study how the reaction terms, the diffusion terms, and the weight functions in the boundary conditions affect the global and blow‐up properties to a porous medium system. The conditions on the global existence and blow‐up in finite time for nonnegative solutions are given. Furthermore, the blow‐up rate estimates of the blow‐up solutions are also established. Copyright © 2011 John Wiley & Sons, Ltd.