非线性双曲型方程的混合元方法

1　引　　言考虑下述非线性双曲型方程的混合问题:c(x,u)utt-.(a(x,u)u)=f(x,u,t),　　x∈Ω,t∈J,(1.1)u(x,0)=u0(x),　　x∈Ω,(1.2)ut(x,0)=u1(x),　　x∈Ω,(1.3)u(x,t)=-g(x,t),　　(x,t)∈Ω×J,(1.4)其中ΩR2是一具有Lipschitz边界Ω的有界区域,J=[0,T],0相似文献   

积分微分方程有限元逼近的强超收敛性

考虑下面的抛物型积分微分方程初边值问题: 　（a） ut+A（t）u+∫0tB（t,s）u（s）ds=f, （x,t）∈Q=Ω×J,J=（0,T] （b） u=0,（x,t）∈　Ω×J,（1） （c）　u（x,0）=u0,x∈Ω,其中Ω为Rd（d≤4）中具有分片光滑边界　Ω的有界域,A（t）是一致正定的二阶椭圆微分算子     

AN EFFECT ITERATION ALGORITHM FOR NUMERICAL SOLUTION OF DISCRETE HAMILTON-JACOBIBELLMAN EQUATIONS

§1Introduction ConsidertheHamilton-Jacobi-Bellmanequation max1≤v≤m[A(v)u(x)-f(v)(x)]=0,x∈Ω(1.1)withtheboundarycondition u(x)=0,x∈Ω(1.2)whereΩisabounded,smoothdomaininEuclideanspaceRd,d∈N;f(v)(x)aregiven functionsfromC2(Ω);A(v)aresecond-orderuniformlyellipticoperatorsoftheform A(v)=-d i,j=1a(v)ij2xixj+di=1b(v)ixi+c(v).(1.3)Intheaboveexpression(1.3)therearecoefficientsa(v)ij,b(v)i,c(v)∈C2(Ω)satisfying,forall1≤v≤m,a(v)ij(x)=a(v)ji(x),1≤i,j≤d,c(v)≥c0≥0,x∈Ω,a…     

一类带临界指数增长的椭圆型方程组两个正解的存在性

该文主要讨论带临界指数的椭圆型方程组{-Δu + a(x)u =2α/α+βuα-1vβ + f(x),x ∈Ω,-Δv+b(x)v=2β/α+βuαvβ-1+ g(x),x ∈ Ω,(*)u ＞ 0,v ＞ 0,x ∈Ω,u=v=0,x ∈(a)Ω解的存在性,其中Ω是RN中一个光滑有界区域,N=3,4,a≥2,β≥2...     

地下水污染中一离子反应数学模型的有限元方法

1　引　　言在地下水含水层中 ,污染物随地下水运移并常常发生各种化学反应 [1 ] .描述地下水含水层中一类阳离子交换反应 m M1 +r M2 k2k1 r M2 +m M1 的数学模型[2 ] 为 : s1 t- dΔs1 =f1 ,　　 x∈Ω ,t∈ J (1.1a) s2 t- dΔs2 =f2 ,　　 x∈Ω ,t∈ J (1.1b) c1 t+ρ s1 t- DΔc1 =0 ,　　 x∈Ω ,t∈ J (1.2 a) c2 t+ρ s2 t- DΔc2 =0 ,　　 x∈Ω ,t∈ J (1.2 b)其中 Ω R2为具有光滑边界的有界区域 ,J=(0 ,T].这里 D>d>0为扩散系数 ,ρ>0为固体颗粒密度 ,均为常数 .根据 [1,2 ]应有 :f1 =m[k1 Rm1 Rr2 cm1 sr2 -…     

一类带弱奇异核非线性偏积分微分方程的全离散有限元

1引言我们将研究下面一类带弱奇异核非线性偏积分微分方程的数值解:u_t-▽·(a(u)▽u)-integral from n=0 to tβ(t-s)△u(s)ds=f(u),x∈Ω,t∈(?),(1.1) u(·,t)=0,x∈(?)Ω,t∈J,(1.2) u(·,0)=v(x),x∈Ω,(1.3)其中Ω为平面上的凸角域,J=(0,T],α和f为R上的光滑函数,满足0相似文献   

三维可压核废料污染问题的交替方向特征有限元格式的敛速分析

1　引　　言三维可压核废料污染问题的数学模型为[1 ] :(a)　φ1 p t+ .u =-q +R′s(c)(b)　φ c t+u . c- .(Ec c) =g(c)(c)　φKi ci t+u . ci - .(Ec ci) +d3(ci) p t=fi(c,c1 ,c2 ,… ,c Nc) ,(i =1 ,… ,Nc)(d)　 d2 T t+cpu . T - .(EH T) +d1 (p) p t=Q(u,T,c,p) (1 .1 )其中 :u=-a(c) p=-k(x)μ(c) p.(x,t)∈Ω×J,Ω=I×I×I,I=(0 ,1 ) ,J=(0 ,T] .假设问题 (1 .1 )满足周期边界条件 ,p(x,t) .c(x,t) .ci(x,t) .T(x,t)的初始条件分别取为 p0 (x) ,c0 (x) ,c0i(x) ,T0 (x) ,(i=1 ,… ,Nc) .假设 (1 .1 )的系数可关…     

具边界反馈时滞的粘弹方程的能量衰减(英文)

考虑如下具边界反馈时滞的粘弹方程ut(x,t)-Δu(x,t)+∫0tg(t-s)Δu(x,s)ds=0,x∈Ω,t0,u(x,t)=0,x∈Γ0,t0,?u /?v=∫0tg(t-s)/vu(s)ds-μ1ut(x,t)-μ2ut(x,t-τ),x∈Γ1,t0,u(x,0)=u0(x),ut(x,0)=u1(x),x∈Ω,ut(x,t-τ)=f0(x,t-τ),x∈Ω,0tτ,其中Ω∈Rn(n≥1)是具C2类边界Ω的有界域.此外,g是所谓的"记忆核",μ1,μ2是两个实数,τ为时滞.在假设|μ2|μ1下,通过构造合适的Lyapunov函数,证明上述问题能量的一般衰减性,使得指数型衰减和多项式衰减仅仅是其特殊情况.     

认真当好一个医生

考虑如下具边界反馈时滞的粘弹方程ut(x,t)-Δu(x,t)+∫0tg(t-s)Δu(x,s)ds=0,x∈Ω,t>0,u(x,t)=0,x∈Γ0,t>0,?u /?v=∫0tg(t-s)/vu(s)ds-μ1ut(x,t)-μ2ut(x,t-τ),x∈Γ1,t>0,u(x,0)=u0(x),ut(x,0)=u1(x),x∈Ω,ut(x,t-τ)=f0(x,t-τ),x∈Ω,0相似文献   

Sobolev方程的基于H~1-Galerkin混合方法的新分裂格式

正1引言考虑如下Sobolev方程u_t-▽·(a(x)▽u_t+a(x)▽u)+u=f(x,t),(x,t)∈Ω×J,u(x,t)=0,(x,t)∈аΩ×J,(1)u(x,0)=u_0(x),x∈Ω.其中Ω是R~d(d=1,2,3)中具有边界     

ERROR ESTIMATES FOR SPARSE OPTIMAL CONTROL PROBLEMS BY PIECEWISE LINEAR FINITE ELEMENT APPROXIMATION

Optimization problems with L¹-control cost functional subject to an elliptic partial differential equation(PDE)are considered.However,different from the finite dimensiona l¹-regularization optimization,the resulting discretized L¹norm does not have a decoupled form when the standard piecewise linear finite element is employed to discretize the continuous problem.A common approach to overcome this difficulty is employing a nodal quadrature formula to approximately discretize the L¹-norm.In this paper,a new discretized scheme for the L¹-norm is presented.Compared to the new discretized scheme for L¹-norm with the nodal quadrature formula,the advantages of our new discretized scheme can be demonstrated in terms of the order of approximation.Moreover,finite element error estimates results for the primal problem with the new discretized scheme for the L¹-norm are provided,which confirms that this approximation scheme will not change the order of error estimates.To solve the new discretized problem,a symmetric Gauss-Seidel based majorized accelerated block coordinate descent(sGS-mABCD)method is introduced to solve it via its dual.The proposed sGS-mABCD algorithm is illustrated at two numerical examples.Numerical results not only confirm the finite element error estimates,but also show that our proposed algorithm is efficient.     

热传导方程一个新的H~1-Galerkin非协调混合有限元格式

对热传导方程提出了一个新的H~1-Galerkin非协调混合有限元格式,其逼近空间不需满足LBB相容性条件,且在不引进传统的Rutz投影的情况下,得到了与以往协调有限元方法相同的L~2-模和H~1-模的误差估计.     

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper,a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure.The convergence analysis is presented and optimal error estimates of both broken H1-norm and L2-norm for velocity as well as the L2-norm for the pressure are derived.     

定常Navier-Stokes方程的一个二阶非协调三角型混合元格式（英文）

In this paper, a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure. The convergence analysis is presented and optimal error estimates of both broken H~1-norm and L~2-norm for velocity as well as the L~2-norm for the pressure are derived.     

非稳态奇异系数微分方程的时间间断时空有限元方法

针对一类二维单轴奇异系数非稳态问题构造了一种时间间断时空有限元格式,利用以Radau点为节点的Lagrange插值多项式的特性,结合有限差分法和有限元法的技巧证明了格式的稳定性和有限元解的时间最大模、空间加权L²（?）-模误差估计.最后列举了一些数值试验结果,验证了理论结果和格式的可行性.     

A symmetric finite volume scheme for selfadjoint elliptic problems

Based on a linear finite element space, a symmetric finite volume scheme for a self-adjoint elliptic boundary-value problem is proposed. Error estimates in L²-norm, H¹-norm, and L^∞-norm are derived. Some post-processing techniques are also provided.     

非定常Navier-Stokes方程的质量集中各向异性非协调有限元逼近

该文将一个低阶Crouzeix-Raviart型非协调三角形元应用到非定常Navier-Stokes方程,给出了其质量集中有限元逼近格式.在不需要传统Ritz-Volterra投影下,通过引入两个辅助有限元空间对边界进行估计的技巧,在各向异性网格下导出了速度的L~2模和能量模及压力的L~2模的误差估计.     

ALTERNATING DIRECTION FINITE ELEMENT METHOD FOR SOME REACTION DIFFUSION MODELS

This paper is concerned with some nonlinear reaction - diffusion models. To solve this kind of models, the modified Laplace finite element scheme and the alternating direction finite element scheme are established for the system of patrical differential equations. Besides, the finite difference method is utilized for the ordinary differential equation in the models. Moreover, by the theory and technique of prior estimates for the differential equations, the convergence analyses and the optimal L^2- norm error estimates are demonstrated.     

Finite volume element approximation and analysis for a kind of semiconductor device simulation

In this paper, we present a finite volume element scheme for a kind of two dimensional semiconductor device simulation. A general framework is developed for finite volume element approximation of the semiconductor problems. We construct a fully discrete finite volume element scheme based on triangulations with a piecewise linear finite element space and a general type of control volume. Optimal-order convergence in H ¹-norm is derived.     

A NEW NONCONFORMING MIXED FINITE ELEMENT SCHEME FOR THE STATIONARY NAVIER-STOKES EQUATIONS

In this article, a new stable nonconforming mixed finite element scheme is proposed for the stationary Navier-Stokes equations, in which a new low order CrouzeixRaviart type nonconforming rectangular element is taken for approximating space for the velocity and the piecewise constant element for the pressure. The optimal order error estimates for the approximation of both the velocity and the pressure in L2-norm are established, as well as one in broken H1-norm for the velocity. Numerical experiments are given which are consistent with our theoretical analysis.