A parameter‐uniform scheme for the parabolic singularly perturbed problem with a delay in time

In this paper, a parameter‐uniform numerical scheme for the solution of singularly perturbed parabolic convection–diffusion problems with a delay in time defined on a rectangular domain is suggested. The presence of the small diffusion parameter ? leads to a parabolic right boundary layer. A collocation method consisting of cubic B ‐spline basis functions on an appropriate piecewise‐uniform mesh is used to discretize the system of ordinary differential equations obtained by using Rothe's method on an equidistant mesh in the temporal direction. The parameter‐uniform convergence of the method is shown by establishing the theoretical error bounds. The numerical results of the test problems validate the theoretical error bounds.     

Regular Degenerate Separable Differential Operators and Applications

The boundary value problems for differential-operator equations with variable coefficients, degenerated on all boundary are studied. Several conditions for the separability, fredholmness and resolvent estimates in L _p -spaces are given. In applications degenerate Cauchy problem for parabolic equation, boundary value problems for degenerate partial differential equations and systems of degenerate elliptic equations on cylindrical domain are studied.     

Stochastic partial differential equations in bounded domains with Dirichlet boundary conditions

In this paper we prove the existence of a unique solution for a class of stochastic parabolic partial differential equations in bounded domains, with Dirichlet boundary conditions. The main tool is an equivalence result, provided by the stochastic characteristics method, between the stochastic equations under investigation and a class of deterministic parabolic equations with moving boundaries, depending on random coefficients. We show the existence of the solution to this last problem, thus providing a solution to the former.     

A layer adaptive B‐spline collocation method for singularly perturbed one‐dimensional parabolic problem with a boundary turning point

In this article, we develop a parameter uniform numerical method for a class of singularly perturbed parabolic equations with a multiple boundary turning point on a rectangular domain. The coefficient of the first derivative with respect to x is given by the formula a₀(x, t)x^p, where a₀(x, t) ≥ α > 0 and the parameter p ∈ [1,∞) takes the arbitrary value. For small values of the parameter ε, the solution of this particular class of problem exhibits the parabolic boundary layer in a neighborhood of the boundary x = 0 of the domain. We use the implicit Euler method to discretize the temporal variable on uniform mesh and a B‐spline collocation method defined on piecewise uniform Shishkin mesh to discretize the spatial variable. Asymptotic bounds for the derivatives of the solution are established by decomposing the solution into smooth and singular component. These bounds are applied in the convergence analysis of the proposed scheme on Shishkin mesh. The resulting method is boundary layer resolving and has been shown almost second‐order accurate in space and first‐order accurate in time. It is also shown that the proposed method is uniformly convergent with respect to the singular perturbation parameter ε. Some numerical results are given to confirm the predicted theory and comparison of numerical results made with a scheme consisting of a standard upwind finite difference operator on a piecewise uniform Shishkin mesh. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1143–1164, 2011     

Homogenization of the Signorini boundary‐value problem in a thick junction and boundary integral equations for the homogenized problem

We consider a mixed boundary‐value problem for the Poisson equation in a thick junction Ω_ε which is the union of a domain Ω₀ and a large number of ε—periodically situated thin cylinders. The non‐uniform Signorini conditions are given on the lateral surfaces of the cylinders. The asymptotic analysis of this problem is done as ε→0, i.e. when the number of the thin cylinders infinitely increases and their thickness tends to zero. We prove a convergence theorem and show that the non‐uniform Signorini boundary conditions are transformed in the limiting variational inequalities in the region that is filled up by the thin cylinders as ε→0. The convergence of the energy integrals is proved as well. The existence and uniqueness of the solution to this non‐standard limit problem is established. This solution can be constructed by using a penalty formulation and successive iteration. For some subclass, these problems can be reduced to an obstacle problem in Ω₀ and an appropriate postprocessing. The equations in Ω₀ finally are also treated with boundary integral equations. Copyright © 2010 John Wiley & Sons, Ltd.     

Partial compactness for Landau-Lifshitz maxwell equation in two-dimension

We study the partial regularity of weak solutions to the 2-dimensional LandauLifshitz equations coupled with time dependent Maxwell equations by Ginzburg-Landau type approximation. Outside an energy concentration set of locally finite 2-dimensional parabolic Hausdorff measure, we prove the uniform local C ∞ bounds for the approaching solutions and then extract a subsequence converging to a global weak solution of the Landau-Lifshitz-Maxwell equations which are smooth away from finitely many points.     

An inverse problem of identifying the coefficient of first-order in a degenerate parabolic equation

This work studies an inverse problem of determining the first-order coefficient of degenerate parabolic equations using the measurement data specified at a fixed internal point. Being different from other ordinary parameter identification problems in parabolic equations, in our mathematical model there exists degeneracy on the lateral boundaries of the domain, which may cause the corresponding boundary conditions to go missing. By the contraction mapping principle, the uniqueness of the solution for the inverse problem is proved. A numerical algorithm on the basis of the predictor-corrector method is designed to obtain the numerical solution and some typical numerical experiments are also performed in the paper. The numerical results show that the proposed method is stable and the unknown function is recovered very well. The results obtained in the paper are interesting and useful, and can be extended to other more general inverse coefficient problems of degenerate PDEs.     

Complete flux scheme for parabolic singularly perturbed differential‐difference equations

In this study, we investigate the concept of the complete flux (CF) obtained as a solution to a local boundary value problem (BVP) for a given parabolic singularly perturbed differential‐difference equation (SPDDE) with modified source term to propose an efficient complete flux‐finite volume method (CF‐FVM) for parabolic SPDDE which is μ‐ and ?‐uniform method where μ, ? are shift and perturbation parameters, respectively. The proposed numerical method is shown to be consistent, stable, and convergent and has been successfully implemented on three test problems.     

On global regular solutions of third order partial differential equations

Diffusion in the presence of high-diffusivity paths is an important issue of current technology. In metals high-diffusivity paths are identified with dislocations, grain boundaries, free surfaces and internal microcracks. Diffusion in a media with two distinct families of diffusion paths is modelled by two coupled linear partial differential equations of parabolic type with diffusivities D₁ and D₂. Physically the situation D₂ ? D₁ is of some considerable interest and previously established results, for D₂ non-zero, for the solution of boundary value problems, are not applicable to the idealized theory characterized by D₂ vanishing. An integral equation, which arises in the solution of boundary value problems for this idealized theory, is formally solved.     

一类特殊空间上映射的梯度几何流柯西问题解的存在唯一性

通过几何分析方法与抛物型方程组解的逼近理论,研究特殊空间(一维球面S~1到二维球面S~2)上映射的梯度几何流柯西问题解的存在唯一性.利用能量法和空间本身特有的性质来解决能量守恒的问题,并利用适当的抛物型方程组逼近该梯度几何流,在适当的Sobolev空间中建立先验估计,找到其时间的一致正下界和抛物型方程组一列解的Sobo1ev范数的一致边界,借助于抛物型偏微分方程的理论,以此决定该柯西问题解的存在唯一性.     

Determination of an unknown function in a parabolic equation with an overspecified condition

In this paper we consider the inverse problems of identifying some space-dependent unknown coefficients in parabolic equations subject to initial boundary value conditions along with an overspecified condition at the final time t = T. We use the overspecified information to transform the problems into non-linear parabolic equations involving a functional of the solution with respect to the time variable. This transformation allows us to establish existence theorems for these inverse problems by employing the Schauder fixed-point theorem.     

On the combination of some semi‐discretization methods and boundary integral equations for the numerical solution of initial boundary value problems

We consider initial boundary value problems for the homogeneous differential equation of hyperbolic or parabolic type in the unbounded two‐ or three‐dimensional spatial domain with the homogeneous initial conditions and with Dirichlet or Neumann boundary condition. The numerical solution is realized in two steps. At first using the Laguerre transformation or Rothe's method with respect to the time variable the non‐stationary problem is reduced to the sequence of boundary value problems for the non‐homogeneous Helmholtz equation. Further we construct the special integral representation for solutions and obtain the sequence of boundary integral equations (without volume integrals). For the full‐discretization of integral equations we propose some projection methods.     

DISCRETE APPROXIMATIONS FOR SINGULARLY PERTURBED BOUNDARY VALUE PROBLEMS WITH PARABOLIC LAYERS,Ⅱ

l)ThisworkwassupportedbyNWOthroughgrantIBo7-3Go12.BOUNDAarv^LUEPRoBLEMFORELLIPTICEQUMIONwiTHMIXEDBOUNDAavCONDITION1.IntroductionInthispedwesketchavarietyofspecialmethodswhichareusedforconstructinge-unifornilyconvergelltschemes-WeshaJldemonstrateamethodwhichachieveshaprovedaccuracyforsolvingsingularlyperturbedb0undaryvalueproblemforeiliPicequatiouswithparabolicboundarylayers-InSecti0n4weshallintroduceanaturalclass,B,oftritefferenceschemes,inwhich(bytheabovementi0nedaP…     

Discrete Approximations for Singularly Perturbed Boundary Value Problems with Parabolic Layers,III

1.IntroductionThesolution0fpartialdifferentiaJequationsthataresingularlyperturbedand/orhavediscontinu0usboundaryconditionsgenerallyhave0nlylimitedsmoothness.DuetothisfaCtdndcultiesaPpearwhenwesolvethesepr0blemsbynumericalmethods.Forexampleforregularparab0licequationswithdiscontinuousboundaryconditions,classicalmethods(FDMorFEM)onregularrectangulargridsd0n0tconvergeintheIoo-normonadomainthatincludesaneighbourhood0fthediscontinulty[8,9,4].Iftheparametermultiplyingthehighest-orderderivativeva…     

The use of cubic B‐spline scaling functions for solving the one‐dimensional hyperbolic equation with a nonlocal conservation condition

Problems for parabolic partial differential equations with nonlocal boundary conditions have been studied in many articles, but boundary value problems for hyperbolic partial differential equations have so far remained nearly uninvestigated. In this article a numerical technique is presented for the solution of a nonclassical problem for the one‐dimensional wave equation. This method uses the cubic B‐spline scaling functions. Some numerical results are reported to support our study. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Discrete Approximations for Singularly Perturbed Boundary Value Problems with Parabolic Layers,I

In this series of three papers we study singularly perturbed (SP) boundary value problems for equations of elliptic and parabolic type. For small values of the perturbation parameter parabolic boundary and interior layers appear in these problems. If classical discretisation methods are used, the solution of the finite difference scheme and the approximation of the diffusive flux do not converge uniformly with respect to this parameter. Using the method of special, adapted grids, we can construct difference schemes that allow approximation of the solution and the normalised diffusive flux uniformly with respect to the small parameter. We also consider singularly perturbed boundary value problems for convection-diffusion equations. Also for these problems we construct special finite difference schemes, the solution of which converges $ε$-uniformly. We study what problems appear, when classical schemes are used for the approximation of the spatial derivatives. We compare the results with those obtained by the adapted approach. Results of numerical experiments are discussed. In the three papers we first give an introduction on the general problem, and then we consider respectively (i) Problems for SP parabolic equations, for which the solution and the normalised diffusive fluxes are required; (ii) Problems for SP elliptic equations with boundary conditions of Dirichlet, Neumann and Robin type; (iii) Problems for SP parabolic equation with discontinuous boundary conditions.     

Some remarks on a class of inverse problems related to the parabolic approximation to the Maxwell equations: a controllability approach

We consider a class of inverse source problems for the parabolic approximation to the Maxwell equations.We relate this to an exact controllability problem; the regularisation of the considered source problems is studied with an optimal control method. Copyright © 2014 John Wiley & Sons, Ltd.     

Natural frequencies of thin rectangular plates using homotopy-perturbation method

The homotopy perturbation method (HPM) was developed to search for asymptotic solutions of nonlinear problems involving parabolic partial differential equations with variable coefficients. This paper illustrates that HPM be easily adapted to solve parabolic partial differential equations with constant coefficients. Natural frequencies of a rectangular plate of uniform thickness, simply-supported on all sides, are obtained with minimum amount of computation. The solution is shown to converge rapidly to a combination of sine and cosine functions. Truncating the series solution by using only the first three terms of the sine and cosine functions as compared to the exact solution results in an absolute error not exceeding 2 × 10⁻⁴ and 9×10⁻⁴ for the trigonometric functions respectively. HPM is then applied to solve the nonlinear problem of a rectangular plate of variable thickness. A direct expression for the eigenvalues (natural frequencies) of the rectangular plate is obtained as compared to determining its eigenvalues by solving the characteristic equation using the conventional method. Comparison of results for the frequency parameter with existing literature show that HPM is highly efficient and accurate. Natural frequencies of a simply-supported guitar soundboard were obtained using an equivalent rectangular plate with the same boundary condition.     

Vanishing viscosity limit of the Navier‐Stokes equations to the euler equations for compressible fluid flow

We establish the vanishing viscosity limit of the Navier‐Stokes equations to the isentropic Euler equations for one‐dimensional compressible fluid flow. For the Navier‐Stokes equations, there exist no natural invariant regions for the equations with the real physical viscosity term so that the uniform sup‐norm of solutions with respect to the physical viscosity coefficient may not be directly controllable. Furthermore, convex entropy‐entropy flux pairs may not produce signed entropy dissipation measures. To overcome these difficulties, we first develop uniform energy‐type estimates with respect to the viscosity coefficient for solutions of the Navier‐Stokes equations and establish the existence of measure‐valued solutions of the isentropic Euler equations generated by the Navier‐Stokes equations. Based on the uniform energy‐type estimates and the features of the isentropic Euler equations, we establish that the entropy dissipation measures of the solutions of the Navier‐Stokes equations for weak entropy‐entropy flux pairs, generated by compactly supported C² test functions, are confined in a compact set in H^?1, which leads to the existence of measure‐valued solutions that are confined by the Tartar‐Murat commutator relation. A careful characterization of the unbounded support of the measure‐valued solution confined by the commutator relation yields the reduction of the measurevalued solution to a Dirac mass, which leads to the convergence of solutions of the Navier‐Stokes equations to a finite‐energy entropy solution of the isentropic Euler equations with finite‐energy initial data, relative to the different end‐states at infinity. © 2010 Wiley Periodicals, Inc.     

On stable implicit difference scheme for hyperbolic–parabolic equations in a Hilbert space

The first‐order of accuracy difference scheme for approximately solving the multipoint nonlocal boundary value problem for the differential equation in a Hilbert space H, with self‐adjoint positive definite operator A is presented. The stability estimates for the solution of this difference scheme are established. In applications, the stability estimates for the solution of difference schemes of the mixed type boundary value problems for hyperbolic–parabolic equations are obtained. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009