SOME TYPICAL BEHAVIORS OF WEAK SOLUTIONS OF LAYERED POROUS MEDIA EQUATION

1.IntroductionThefiltrationproblemoflayeredporousmediacanbedescribedasfollows(see[21):whereHisHeavisidefunction.u(x,t)isthecapillarypiezometrichead,oi(u)(i~l,2)thevolumetricmoisturecontentandK4(u)thehydraulicconductivityinupperlayer(x6(1,21).Theassumptionsarethesameasthosein[2]:(0KI)oi(r),Ki(r)EC' '(R).(0K2)oi(r)>0,Ki(r)>0,6:(r)>0andK:(r)20forrE(--co,u'.),whereu"thesaturatedcapillarypiezometricheadoftheithlayer.Whenr2u",oi(r)=osiandKi(r)~K'i,whereof'andK8iaresaturatedvaluesofthemoisturec…     

Some typical behaviors of weak solutions of layered proous media equation

This paper is the continuation of [2]. Some typical behaviors of weak solutions of layered porous media equations with boundary conditions will be discussed in this paper. For example, asymptotically, the saturated regions can appear only either near the layered interface, or near the boundaries. The necessary and sufficient conditions for the occurrence of such phenomena will be given.     

Fast integral equation solver for Maxwell’s equations in layered media with FMM for Bessel functions

The paper presents a new fast integral equation solver for Maxwell’s equations in 3-D layered media. First, the spectral domain dyadic Green’s function is derived, and the 0-th and the 1-st order Hankel transforms or Sommerfeld-type integrals are used to recover all components of the dyadic Green’s function in real space. The Hankel transforms are performed with the adaptive generalized Gaussian quadrature points and window functions to minimize the computational cost. Subsequently, a fast integral equation solver with O(N _z ² N _x N _y log(N _x N _y )) in layered media is developed by rewriting the layered media integral operator in terms of Hankel transforms and using the new fast multipole method for the n-th order Bessel function in 2-D. Computational cost and parallel efficiency of the new algorithm are presented.     

层状多孔介质中渗流湿峰的不稳定分析

在部分饱和与部分干的层状多孔介质中，建立了一维渗流问题的数学模型．在自然的条件下，证明了一类退化抛物方程弱解的存在唯一性．揭示了一个重要性质，即当下层介质的孔隙比上层更粗时，层状交界面具有隔水作用．这一作用将会在展状交界面不均匀时导致指状湿峰的出现．同时也说明，较粗孔隙的上层具有防止水分蒸发的作用．     

Preconditioning iterative algorithm for the electromagnetic scattering from a large cavity

A preconditioning iterative algorithm is proposed for solving electromagnetic scattering from an open cavity embedded in an infinite ground plane. In this iterative algorithm, a physical model with a vertically layered medium is employed as a preconditioner of the model of general media. A fast algorithm developed in (SIAM J. Sci. Comput. 2005; 27 :553–574) is applied for solving the model of layered media and classical Krylov subspace methods, restarted GMRES, COCG, and BiCGstab are employed for solving the preconditioned system. Our numerical experiments on cavity models with large numbers of mesh points and large wave numbers show that the algorithm is efficient and the number of iterations is independent of the number of mesh points and dependent upon the wave number. Copyright © 2008 John Wiley & Sons, Ltd.     

A Second-Order Method for the Electromagnetic Scattering from a Large Cavity

In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.     

Numerical solution of electromagnetic scattering from a large partly covered cavity

The paper focuses on the numerical study of electromagnetic scattering from two-dimensional (2D) large partly covered cavities, which is described by the Helmholtz equation with a nonlocal boundary condition on the aperture. The classical five-point finite difference method is applied for the discretization of the Helmholtz equation and a linear approximation is used for the nonlocal boundary condition. We prove the existence and uniqueness of the numerical solution when the medium in the cavity is y-direction layered or the number of the mesh points on the aperture is large enough. The fast algorithm proposed in Bao and Sun (2005) [2] for open cavity models is extended to solving the partly covered cavity problem with (vertically) layered media. A preconditioned Krylov subspace method is proposed to solve the partly covered cavity problem with a general medium, in which a layered medium model is used as a preconditioner of the general model. Numerical results for several types of partly covered cavities with different wave numbers are reported and compared with those by ILU-type preconditioning algorithms. Our numerical experiments show that the proposed preconditioning algorithm is more efficient for partly covered cavity problems, particularly with large wave numbers.     

On 3D orthogonal prolate spheroidal monogenics

S. G. Georgiev, Complete orthogonal systems of monogenic polynomials over 3D prolate spheroids have recently experienced an upsurge of interest because of their many remarkable properties. These generalized polynomials and their applications to the theory of quasi‐conformal mappings and approximation theory have played a major role in this development. In particular, the underlying functions of three real variables take on values in the reduced quaternions (identified with ) and are generally assumed to be null‐solutions of the well‐known Riesz system in . The present paper introduces and explores a new complete orthogonal system of monogenic functions as solutions to this system for the space exterior of a 3D prolate spheroid. This will be made in the linear spaces of square integrable functions over . The representations of these functions are explicitly given. Some important properties of the system are briefly discussed, from which several recurrence formulae for fast computer implementations can be derived. Copyright © 2015 John Wiley & Sons, Ltd.     

Diffraction of a Plane Discontinuos Wave in Layered Anisotropic Elastic Media

The problem of transformation of the front of a plane longitudinal shock wave by curvilinear interfaces of layered anisotropic elastic media with different mechanical properties is considered. To solve the nonlinear Snell equations, the approach based on a synthesis of the Newton method and the algorithm of parametric continuation of a solution is used. The cases of focusing and scattering of rays by convex and concave interfaces of layered composite media and by convex lenses are investigated. Numerical examples are given for various ratios between the parameters of the elastic media.     

Dynamic response of a stratified transversely isotropic half-space with a poroelastic interlayer due to a buried moving source

Response of a ground subjected to dynamic sources is a classical problem in many fields. This paper develops a transversely isotropic (TI) half-space ground model which considers the alternating distribution of different media: a TI elastic surface layer, a TI poroelastic interlayer and a TI elastic half-space from top to bottom. The model can consider soils with different properties and the existence of groundwater. Dynamic responses of the model subjected to a buried moving source are investigated. General solutions to the TI poroelastic interlayer in the wavenumber domain are derived by Fourier transform and a potential function method in a Cartesian coordinate system. The layered system with different media is solved adopting an ‘adapted stiffness matrix method’ for a buried source scenario. Dynamic responses in the time-spatial domain are obtained by inverse fast Fourier transform (iFFT) and numerical integration. The performance of the method in this study is first tested by comparison with existing research, and the influences of the existence of pore water in the interlayer, transverse isotropy, and source features are comprehensively investigated. The effect of pore water is more significant for a TI ground compared with isotropic ground, and the model in this study cannot be simplified by an ‘equivalent’ layered TI elastic model. Peak values of displacement and pore pressure exist with increasing source speed. Source speeds corresponding to them are related to the wave speeds of the medium and affected by transverse isotropy.     

The Schur algorithm and its applications

The Schur algorithm and its time-domain counterpart, the fast Cholseky recursions, are some efficient signal processing algorithms which are well adapted to the study of inverse scattering problems. These algorithms use a layer stripping approach to reconstruct a lossless scattering medium described by symmetric two-component wave equations which model the interaction of right and left propagating waves. In this paper, the Schur and fast Chokesky recursions are presented and are used to study several inverse problems such as the reconstruction of nonuniform lossless transmission lines, the inverse problem for a layered acoustic medium, and the linear least-squares estimation of stationary stochastic processes. The inverse scattering problem for asymmetric two-component wave equations corresponding to lossy media is also examined and solved by using two coupled sets of Schur recursions. This procedure is then applied to the inverse problem for lossy transmission lines.The work of this author was supported by the Exxon Education FoundationThe work of this author was supported by the Air Force Office of Scientific Research under Grant AFOSR-82-0135A.     

Solution of and bounding in a linearly constrained optimization problem with convex,polyhedral objective function

A dual method is presented to solve a linearly constrained optimization problem with convex, polyhedral objective function, along with a fast bounding technique, for the optimum value. The method can be used to solve problems, obtained from LPs, where some of the constraints are not required to be exactly satisfied but are penalized by piecewise linear functions, which are added to the objective function of the original problem. The method generalizes an earlier solution technique developed by Prékopa (1990). Applications to stochastic programming are also presented.This research was supported by the National Science Foundation, Grant No. DMS-9005159.Corresponding author.     

On methods of deriving equations describing effective models of layered media

Methods of deriving equations describing effective models of layered periodic media are presented. Elastic and fluid media, as well as porous Biot media, may be among these media. First, effective models are derived by a rigorous method, and then some operations in the derivation are replaced by simpler ones providing correct results. As a consequence, a comparatively simple and justified method of deriving equations of an effective model is established. In particular, this method allows us to simplify to a degree and justify the derivation of an effective model for media containing Biot layers; this method also produces equations of an effective model of a porous layered medium intersected by fractures with slipping contacts. Bibliography: 15 titles. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 250, 1998 pp. 219–243. Translated by L. A. Molotkov.     

Plateaued functions,partial geometric difference sets,and partial geometric designs

Plateaued functions on finite fields have been studied in many papers in recent years. As a generalization of plateaued functions on finite fields, we introduce the notion of a plateaued function on a finite abelian group. We will give a characterization of a plateaued function in terms of an equation of the matrix associated to the function. Then we establish a one‐to‐one correspondence between the ${\mathbb{Z}}_2$‐valued plateaued functions and partial geometric difference sets (with specific parameters) in finite abelian groups. We will also discuss two general methods (extension and lifting) for the construction of new partial geometric difference sets from old ones in (abelian or nonabelian) finite groups, and construct many partial geometric difference sets and plateaued functions. A one‐to‐one correspondence between partial geometric difference sets (in arbitrary finite groups) and partial geometric designs will be proved.     

Surface acoustic waves in the testing of layered media. The waves’ sensitivity to variations in the properties of the individual layers

Research on the use of surface acoustic waves for the nondestructive testing of layered media is reviewed. A model to describe horizontally polarized surface acoustic waves in layered anisotropic (monoclinic) media is constructed. A modified transfer-matrix method is developed to obtain a solution. Non-canonical type waves with horizontal transverse polarization are investigated. Dispersion curves are constructed for a multilayer composite in contact with an anisotropic half-space. It is shown that the variation of the physical characteristics and the geometry of any of the internal layers leads to a variation in the dispersion curves. This opens up the possibility of using dispersion analysis for the nondestructive testing of the properties of the individual layers.     

Some problems of the correlation theory of structurally inhomogeneous elastic media

The correlation theory of random functions is used to describe structurally inhomogeneous elastic media with randomly distributed inhomogeneities. Expressions are obtained for the elastic moduli of averaged layered, fibrous and granular media.Mekhanika Polimerov, Vol. 2, No. 3, pp. 365–371, 1966     

On some length problems for univalent functions

Let be the class of functions which are analytic in the unit disk . Let C(r) be the closed curve that is the image of the circle |z|=r < 1 under the mapping w = f(z), L(r) the length of C(r), and let A(r) be the area enclosed by the curve C(r). In 1968 D. K. Thomas shown that if , f is starlike with respect to the origin, and for 0≤r < 1, A(r) < A, an absolute constant, then Later, in 1969 Nunokawa has shown that if f is convex univalent, then This paper is devoted to obtaining a related correspondence between f(z) and L(r) for the case when f is univalent. Copyright © 2016 John Wiley & Sons, Ltd.     

A note on anisotropic,inhomogeneous, poro‐elastic media

A modification of the material law associated with the well‐known Biot system as suggested by Murad and Cushman (Int. J. Eng. Sci. 1996; 34(3):313–338) and first investigated by Showalter (J. Math. Anal. Appl. 2000; 251(1):310–340) is reconsidered, generalized and analysed in the light of a new approach to a comprehensive class of evolutionary problems. The framework allows a uniform approach to problems involving general anisotropic, inhomogeneous, non‐smooth media thus covering, for example, transmission problems in layered materials. Copyright © 2009 John Wiley & Sons, Ltd.     

Bicompact Finite-Difference Scheme for Maxwell’s Equations in Layered Media

Doklady Mathematics - In layered media, the solution of Maxwell’s equations suffers a strong or weak discontinuity at the layer boundaries. Finite-difference schemes providing convergence on...     

Tensor Green’s functions for axisymmetrical models of magnetotelluric marine sounding of locally nonhomogeneous layered media

A method for computing the electromagnetic fields of point (dipole) sources is proposed for magnetotelluric sounding (MTS) problems in axisymmetrical conducting layered media. The method expands the tensor Green’s functions of the layered medium in Fourier series in the azimuthal coordinate. For an arbitrary system of point sources we construct algorithms to compute the electromagnetic fields propagating across the plane interface of two conducting half-spaces with different constant conductivities. Translated from Prikladnaya Matematika i Informatika, No. 30, pp. 5–17, 2008.