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.     

一类带有铁磁材料参数的非线性涡流问题的A-φ有限元法

针对带有铁磁材料的非线性涡流问题,其非线性性通常体现在磁场强度和磁感应强度的关系上.本文提出了一种全离散的有限元A-φ格式,分别在时间和空间上采用向后欧拉公式以及节点有限元进行离散.首先,在合适的函数空间里给出时间上的半离散格式,通过考察其弱形式建立相应的适定性理论,并证明近似解收敛于弱解.其次,给出全离散格式并讨论其误差估计.最后,给出两个数值算例以验证理论结果.     

守恒型扩散方程非线性离散格式的性质分析和快速求解

对于守恒型扩散方程,研究其二阶时间精度非线性全隐有限差分离散格式的性质,证明了其解的存在唯一性.研究了二阶时间精度的Picard-Newton迭代格式,证明了迭代解对原问题真解的二阶时间和空间收敛性,以及对非线性离散解的二次收敛速度,实现了非线性问题的快速求解.本文中方法也适用于一阶时间精度格式的分析,并可推广至对流扩散问题.数值实验验证了二阶时间精度Picard-Newton迭代格式的高精度和高效率.     

Existence and asymptotic stability of relaxation discrete shock profiles

In this paper we study the asymptotic nonlinear stability of discrete shocks of the relaxing scheme for approximating the general system of nonlinear hyperbolic conservation laws. The existence of discrete shocks is established by suitable manifold construction, and it is shown that weak single discrete shocks for such a scheme are nonlinearly stable in , provided that the sums of the initial perturbations equal zero. These results should shed light on the convergence of the numerical solution constructed by the relaxing scheme for the single shock solution of the system of hyperbolic conservation laws. These results are proved by using both a weighted norm estimate and a characteristic energy method based on the internal structures of the discrete shocks.

     

Modeling and analysis of the unilateral contact of a piezoelectric body with a conductive support

We consider a mathematical model which describes the quasistatic process of contact between a piezoelectric body and an electrically conductive support, the so-called foundation. We model the material's behavior with a nonlinear electro-viscoelastic constitutive law; the contact is frictionless and is described with the Signorini condition and a regularized electrical conductivity condition. We derive a variational formulation for the problem and then we prove the existence of a unique weak solution to the model. The proof is based on arguments of nonlinear equations with multivalued maximal monotone operators and fixed point. Then we introduce a fully discrete scheme, based on the finite element method to approximate the spatial variable and the backward Euler scheme to discretize the time derivatives. We treat the unilateral contact conditions by using an augmented Lagrangian approach. We implement this scheme in a numerical code then we present numerical simulations in the study of two-dimensional test problems, together with various comments and interpretations.     

An efficient nonlinear iteration scheme for nonlinear parabolic-hyperbolic system

A nonlinear iteration method named the Picard-Newton iteration is studied for a two-dimensional nonlinear coupled parabolic-hyperbolic system. It serves as an efficient method to solve a nonlinear discrete scheme with second spatial and temporal accuracy. The nonlinear iteration scheme is constructed with a linearization-discretization approach through discretizing the linearized systems of the original nonlinear partial differential equations. It can be viewed as an improved Picard iteration, and can accelerate convergence over the standard Picard iteration. Moreover, the discretization with second-order accuracy in both spatial and temporal variants is introduced to get the Picard-Newton iteration scheme. By using the energy estimate and inductive hypothesis reasoning, the difficulties arising from the nonlinearity and the coupling of different equation types are overcome. It follows that the rigorous theoretical analysis on the approximation of the solution of the Picard-Newton iteration scheme to the solution of the original continuous problem is obtained, which is different from the traditional error estimate that usually estimates the error between the solution of the nonlinear discrete scheme and the solution of the original problem. Moreover, such approximation is independent of the iteration number. Numerical experiments verify the theoretical result, and show that the Picard-Newton iteration scheme with second-order spatial and temporal accuracy is more accurate and efficient than that of first-order temporal accuracy.     

非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅱ):向后一步的Euler全离散化格式

1．引言本文的工作主要是讨论非定常的热传导一对流问题的向后一步的Euler全离散化的非线性Galerkin混合元解的存在性及其误差估计．该工作是对山中的同一问题研究的第二部分．在第一部分［1］，我们已经讨论了此问题的半离散化的情形．由于所研究的目标都是非定常的热传导一对流问题，其背景是相同的，在此将不重复了，请参考［1］．本文的安排如下，52先回顾非定常的热传导一对流问题的混合元解的经典性质．53回顾半离散化的非线性Galerkin混合元解的性质，并导出后续讨论需要的一些关于时间导数的估计．54讨论向后一步的Euler全离散化…     

Fully discrete T‐ψ finite element method to solve a nonlinear induction hardening problem

We study an induction hardening model described by Maxwell's equations coupled with a heat equation. The magnetic induction field is assumed a nonlinear constitutional relation and the electric conductivity is temperature‐dependent. The T ‐ψ method is to transform Maxwell's equations to the vector–scalar potential formulations and to solve the potentials by means of the finite element method. In this article, we present a fully discrete T ‐ψ finite element scheme for this nonlinear coupled problem and discuss its solvability. We prove that the discrete solution converges to a weak solution of the continuous problem. Finally, we conclude with two numerical experiments for the coupled system.     

非定常对流扩散方程保正格式解的存在性

本文发展了非定常对流扩散方程的非线性保正格式.该格式为单元中心型有限体积格式,保持局部通量的守恒性,适用于任意星形多边形网格,本文证明了该离散格式解的存在性,并给出数值结果,表明该格式具有二阶精度.     

Numerical analysis for generalized Forchheimer flows of slightly compressible fluids

In this article, we will consider the generalized Forchheimer flows for slightly compressible fluids. Using Muskat's and Ward's general form of Forchheimer equations, we describe the fluid dynamics by a nonlinear degenerate parabolic equation for density. The long‐time numerical approximation of the nonlinear degenerate parabolic equation with time dependent boundary conditions is studied. The stability for all time is established in a continuous time scheme and a discrete backward Euler scheme. A Gronwall's inequality‐type is used to study the asymptotic behavior of the solution. Error estimates for the solution are derived for both continuous and discrete time procedures. Numerical experiments confirm the theoretical analysis regarding convergence rates.     

An efficient nonlinear iteration scheme for nonlinear parabolic–hyperbolic system

A nonlinear iteration method named the Picard–Newton iteration is studied for a two-dimensional nonlinear coupled parabolic–hyperbolic system. It serves as an efficient method to solve a nonlinear discrete scheme with second spatial and temporal accuracy. The nonlinear iteration scheme is constructed with a linearization–discretization approach through discretizing the linearized systems of the original nonlinear partial differential equations. It can be viewed as an improved Picard iteration, and can accelerate convergence over the standard Picard iteration. Moreover, the discretization with second-order accuracy in both spatial and temporal variants is introduced to get the Picard–Newton iteration scheme. By using the energy estimate and inductive hypothesis reasoning, the difficulties arising from the nonlinearity and the coupling of different equation types are overcome. It follows that the rigorous theoretical analysis on the approximation of the solution of the Picard–Newton iteration scheme to the solution of the original continuous problem is obtained, which is different from the traditional error estimate that usually estimates the error between the solution of the nonlinear discrete scheme and the solution of the original problem. Moreover, such approximation is independent of the iteration number. Numerical experiments verify the theoretical result, and show that the Picard–Newton iteration scheme with second-order spatial and temporal accuracy is more accurate and efficient than that of first-order temporal accuracy.     

非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅱ):向后一步的Euler全离散化格式

In this paper, a fully discrete format of nonlinear Galerkin mixed element method with backward one-step Euler discretization of time for the non stationary conduction-convection problems is presented. The scheme is based on two finite element spaces XH and Xh for the approximation of the velocity, defined respectively on a coarse grid with grids size H and another fine grid with grid size h<< H, a finite element space Mh for the approximation of the pressure and two finite element spaces AH and Wh, for the approximation of the temperature,also defined respectivply on the coarse grid with grid size H and another fine grid with grid size h. The existence and the convergence of the fully discrete mixed element solution are shown. The scheme consists in using standard backward one step Euler-Galerkin fully discrete format at first L0 steps (L0 2) on fine grid with grid size h, but using nonlinear Galerkin mixed element method of backward one step Euler-Galerkin fully discrete format through L0 + 1 step to end step. We have proved that the fully discrete nonlinear Galerkin mixed element procedure with respect to the coarse grid spaces with grid size H holds superconvergence.     

Long Time Asymptotic Behavior of Solution of Difference Scheme for a Semilinear Parabolic Equation

In this paper we prove that the solution of implicit difference scheme for a semilinear parabolic equation converges to the solution of difference scheme for the corresponding nonlinear stationary problem as $t\rightarrow\infty$. For the discrete solution of nonlinear parabolic problem, we get its long time asymptotic behavior which is similar to that of the continuous solution. For simplicity, we consider one-dimensional problem.     

浅水长波近似方程组的拟小波解

以浅水长波近似方程组为例,提出了拟小波方法求解(1 1)维非线性偏微分方程组数值解,该方程用拟小波离散格式离散空间导数,得到关于时间的常微分方程组,用四阶Runge-K utta方法离散时间导数,并将其拟小波解与解析解进行比较和验证.     

Integrable Discrete Model for One‐Dimensional Soil Water Infiltration

We propose an integrable discrete model of one‐dimensional soil water infiltration. This model is based on the continuum model by Broadbridge and White, which takes the form of nonlinear convection–diffusion equation with a nonlinear flux boundary condition at the surface. It is transformed to the Burgers equation with a time‐dependent flux term by the hodograph transformation. We construct a discrete model preserving the underlying integrability, which is formulated as the self‐adaptive moving mesh scheme. The discretization is based on linearizability of the Burgers equation to the linear diffusion equation, but the naïve discretization based on the Euler scheme which is often used in the theory of discrete integrable systems does not necessarily give a good numerical scheme. Taking desirable properties of a numerical scheme into account, we propose an alternative discrete model that produces solutions with similar accuracy to direct computation on the original nonlinear equation, but with clear benefits regarding computational cost.     

A Finite Volume Method Preserving Maximum Principle for the Conjugate Heat Transfer Problems with General Interface Conditions

In this paper, we present a unified finite volume method preserving discrete maximum principle (DMP) for the conjugate heat transfer problems with general interface conditions. We prove the existence of the numerical solution and the DMP-preserving property. Numerical experiments show that the nonlinear iteration numbers of the scheme in [24] increase rapidly when the interfacial coefficients decrease to zero. In contrast, the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero, which reveals that the unified scheme is more robust than the scheme in [24]. The accuracy and DMP-preserving property of the scheme are also veried in the numerical experiments.     

A finite difference approach to the equations of age-dependent population dynamics

We establish the convergence of the finite difference scheme for the nonlinear equations of population dynamics proposed by Guertin and MacCamy. The applicability of the discrete equations to establish qualitative properties of the solution to the continuous problem is also illustrated.     

Finite Element Solution of the Fundamental Equations of Semiconductor Devices II

In part I of the paper (see Zlamal [13]) finite element solutions of the nonstationary semiconductor equations were constructed. Two fully discrete schemes were proposed. One was nonlinear, the other partly linear. In this part of the paper we justify the nonlinear scheme. We consider the case of basic boundary conditions and of constant mobilities and prove that the scheme is unconditionally stable. Further, we show that the approximate solution, extended to the whole time interval as a piecewise linear function, converges in a strong norm to the weak solution of the semiconductor equations. These results represent an extended and corrected version of results announced without proof in Zlamal [14].     

Full discretization of a nonlinear parabolic problem containing volterra operators and an unknown dirichlet boundary condition

The reconstruction of an unknown solely time‐dependent Dirichlet boundary condition in a nonlinear parabolic problem containing a linear and a nonlinear Volterra operator is considered. The inverse problem is converted into a variational problem in which the unknown Dirichlet condition is eliminated using a given integral overdetermination. A time‐discrete recurrent approximation scheme is designed, using Backward Euler's method. The convergence of the approximations towards a solution of the variational problem is proved under appropriate assumptions on the data and on the Volterra operators. The uniqueness of this solution is shown in the case that the nonlinear Volterra operator satisfies a particular inequality. Moreover, the Finite Element Method is used to discretize the time‐discrete approximation scheme in space. Finally, full‐discrete error estimates are derived for a particular choice of the finite elements. The corresponding convergence rates are supported by a numerical experiment. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1444–1460, 2015     

一类非线性反应扩散方程组的小波数值方法

本文研究了生态学中一类非线性反应扩散方程组的小波Galerkin方法,利用多尺度分析的尺度空间作为试探函数空间,建立显式离散模型,证明了小波逼近解的存在唯一性,并进行了误差分析,最后给出数值模拟的例子.