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

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相似文献   

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)中具有边界     

INVARIANT SETS AND THE HUKUHARA-KNESER PROPERTY FOR PARABOLIC SYSTEMS

In this paper we are concerned with the nonlinear boundary value problem forparabolic system(Lu=f(x,t,u,▽u),x∈Ω,0相似文献   

带有流量限制的拟线性趋化模型解的整体有界性

本文考虑如下一类带流量限制的趋化模型:{ut=▽·(up▽u/√u2+|▽u|2)-x▽(uq▽v/√1+|▽v|2),x∈Ω,t＞0,0=△v-μ+u,x∈Ω,t＞0的初边值问题,其中Ω=BR(0)(∈)Rn(n≥2,R＞0),x＞0,μ:=1/|Ω|∫Ωu0(x)dx.本文主要通过区域划分,建立新的先验估计,在p＞...     

线性抛物型积分微分方程的扩展混合体积元方法

1 引言 考虑线性抛物型积分微分方程初边值问题: {pt(x,t)-▽.{A(x,t)▽p(x,t) +∫t0 B(x,t,τ)▽p(x,τ)dτ}=f(x,t),(x,t)∈Ω×(0,T],(1.1) p(x,0):p0(x), x∈Ω, p(x,t)=0, (x,t)∈(a)Ω×(0,T]. 这里x=(x,y),Ω=(a,b)×(c,d),(e)Ω是区域Ω的边界,p为未知函数,A=(aij)2×2为已知的对称正定矩阵,B=(bij)2×2为已知矩阵,而且aij,bij,(aij)t(i,j=1,2)光滑有界,f∈L2(Ω).     

GLOBAL EXISTENCE AND BLOW-UP OF SOLUTIONS FOR GENERALIZED POCHHAMMER-CHREE EQUATIONS

In this article,we study the initial boundary value problem of generalized Pochhammer-Chree equation u_(tt)-u_(xx)-u_(xxt)-u_(xxtt)=f(u) xx,x ∈Ω,t 0,u(x,0) = u0(x),u t(x,0)=u1(x),x ∈Ω,u(0,t) = u(1,t) = 0,t≥0,where Ω=(0,1).First,we obtain the existence of local W k,p solutions.Then,we prove that,if f(s) ∈ΩC k+1(R) is nondecreasing,f(0) = 0 and |f(u)|≤C1|u| u 0 f(s)ds+C2,u 0(x),u 1(x) ∈ΩW k,p(Ω) ∩ W 1,p 0(Ω),k ≥ 1,1 p ≤∞,then for any T 0 the problem admits a unique solution u(x,t) ∈ W 2,∞(0,T;W k,p(Ω) ∩ W 1,p 0(Ω)).Finally,the finite time blow-up of solutions and global W k,p solution of generalized IMBq equations are discussed.     

一类奇异p-Laplace方程无穷多解的存在性

讨论了一类具有奇异系数的p-Laplace问题-Δpu-μ|u|u|x|p=u|x|tu+λuq-2u,x∈Ω,u=0,x∈Ω无穷多解的存在性,其中N≥3,Ω是RN中一有界光滑区域,0∈Ω,Δpu=-div(|▽u|p-2▽u),0≤μ<μ=(N-p)ppp,10,1相似文献   

凸区域上拟线性椭圆型方程的尖峰解

本文讨论奇异扰动的拟线性椭圆型方程-ε△pu(x)=f(u(x)),u(x)≥0,x∈Ω;u=0,x∈Ω在Dirichlet边值条件下极小能量解的存在性和结构.其中ε＞0是小参数,p＞2,△pu=div(|Du|p-2Du),f(s)=sq-sp-1,p-1＜q＜Np/N-p-1.Ω RN(N≥2)是有界光滑区域.当ε→0时,方程存在一个极小能量解,应用移动平面方法可以证明此解在凸区域上会变成一个尖峰解.     

带非奇扰动项的(2, p)-Laplace方程无穷多解的存在性

本文研究带非奇扰动项的(2,p)-Laplace方程{-△u-△pu=a(x)|u|q-2u+f(x,u),u=0, x∈Ω,x∈(e)Ω,其中Ω (∈) RN是有界光滑区域,1＜q＜2＜p＜N,a∈C((Ω))可变号,f关于u不必是奇函数.利用变分方法,本文获得该方程无穷多解的存在性.     

抛物型方程非齐次边值问题的推广型LOD有限差分及有限元格式

1引言本文考虑区域Ω=[0,1]~d(d=2,3)上的非齐次抛物型方程第一边值问题(?)-C_1△u C_2u=f(x,t),x∈Ω,t∈(0,T],(1．1) u(x,0)=u_0(x),x∈Ω,(1．2) u(x,t)=(?)(x,t),x∈(?)Ω,t∈(0,T],(1．3)其中C_1,C_2为常数且C_1＞0,C_2≥0．对于以上问题,可以使用有限差分方法及有限元方法进行离散,并采用交替方向方法求解．交替方向方法能够将高维问题转化为一系列的一维问题进行计算,具有计算量少,计算稳定且易于并行实现等优点,在大规模科学计算中起着非常重要的作用,一直是计算数     

Upwind compact finite difference scheme for time-accurate solution of the incompressible Navier-Stokes equations

This article presents a time-accurate numerical method using high-order accurate compact finite difference scheme for the incompressible Navier-Stokes equations. The method relies on the artificial compressibility formulation, which endows the governing equations a hyperbolic-parabolic nature. The convective terms are discretized with a third-order upwind compact scheme based on flux-difference splitting, and the viscous terms are approximated with a fourth-order central compact scheme. Dual-time stepping is implemented for time-accurate calculation in conjunction with Beam-Warming approximate factorization scheme. The present compact scheme is compared with an established non-compact scheme via analysis in a model equation and numerical tests in four benchmark flow problems. Comparisons demonstrate that the present third-order upwind compact scheme is more accurate than the non-compact scheme while having the same computational cost as the latter.     

Stability and convergence of staggered Runge-Kutta schemes for semilinear wave equations

A staggered Runge-Kutta (staggered RK) scheme is a Runge-Kutta type scheme using a time staggered grid, as proposed by Ghrist et al. in 2000 [6]. Afterwards, Verwer in two papers investigated the efficiency of a scheme proposed by Ghrist et al. [6] for linear wave equations. We study stability and convergence properties of this scheme for semilinear wave equations. In particular, we prove convergence of a fully discrete scheme obtained by applying the staggered RK scheme to the MOL approximation of the equation.     

一种登陆作战兵力上船装载方案优化计算方法

建立登陆作战兵力上船装载方案优化模型,针对大规模登陆作战兵力上船装载方案确定问题特点,提出一种用需装载兵力、可用舰船和可用泊位类型排列表示方案的方案编码方法,称为方案的类型排列基因码.给出排列基因码生成及进化方法,实现装载方案的遗传优化.计算结果表明,使用该方法可以在较短时间内得到较好的上船装载方案.     

An efficient TVD‐WENO method for conservation laws

In this article, we present a high‐resolution hybrid scheme for solving hyperbolic conservation laws in one and two dimensions. In this scheme, we use a cheap fourth order total variation diminishing (TVD) scheme for smooth region and expensive seventh order weighted nonoscillatory (WENO) scheme near discontinuities. To distinguish between the smooth parts and discontinuities, we use an efficient adaptive multiresolution technique. For time integration, we use the third order TVD Runge‐Kutta scheme. The accuracy of the resulting hybrid high order scheme is comparable with these of WENO, but with significant decrease of the CPU cost. Numerical demonstrates that the proposed scheme is comparable to the high order WENO scheme and superior to the fourth order TVD scheme. Our scheme has the added advantage of simplicity and computational efficiency. Numerical tests are presented which show the robustness and effectiveness of the proposed scheme.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A higher-order hybrid spline difference method on adaptive mesh for solving singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions

We propose a hybrid numerical scheme to discretize a class of singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions on an equidistributed grid. The hybrid difference scheme is developed by using a modified backward difference scheme in time, a combination of the cubic spline and exponential spline difference scheme in space. The proposed scheme uses a cubic spline difference scheme for the discretization of robin-boundary conditions. For the time discretization of the problem, we use the standard uniform mesh while a layer adapted equidistributed grid is generated for the spatial discretization. By equidistributing a curvature-based monitor function, the spatial adaptive grid is able to capture the presence of parabolic boundary layers without using any prior information about the solution. Parameter uniform error estimates are derived to illustrate an optimal convergence of first-order in time and second-order in space for the proposed discretization. The accuracy of the proposed scheme is confirmed by the numerical experiments that underpin the theoretical analysis.     

Solution of double nonlinear problems in porous media by a combined finite volume–finite element algorithm

The combined finite volume–finite element scheme for a double nonlinear parabolic convection-dominated diffusion equation which models the variably saturated flow and contaminant transport problems in porous media is extended. Whereas the convection is approximated by a finite volume method (Multi-Point Flux Approximation), the diffusion is approximated by a finite element method. The scheme is fully implicit and involves a relaxation-regularized algorithm. Due to monotonicity and conservation properties of the approximated scheme and in view of the compactness theorem we show the convergence of the numerical scheme to the weak solution. Our scheme is applied for computing two dimensional examples with different degrees of complexity. The numerical results demonstrate that the proposed scheme gives good performance in convergence and accuracy.     

A fourth-order method of the convection-diffusion equations with Neumann boundary conditions

In this paper, we have developed a fourth-order compact finite difference scheme for solving the convection-diffusion equation with Neumann boundary conditions. Firstly, we apply the compact finite difference scheme of fourth-order to discrete spatial derivatives at the interior points. Then, we present a new compact finite difference scheme for the boundary points, which is also fourth-order accurate. Finally, we use a Padé approximation method for the resulting linear system of ordinary differential equations. The presented scheme has fifth-order accuracy in the time direction and fourth-order accuracy in the space direction. It is shown through analysis that the scheme is unconditionally stable. Numerical results show that the compact finite difference scheme gives an efficient method for solving the convection-diffusion equations with Neumann boundary conditions.     

ON THE CENTRAL RELAXING SCHEMES I:SINGLE CONSERVATION LAWS

1. IntroductionIn [6], Jin and adn constructed a class of uPWind relaxing schemes for nonlinearconservation lawswith initial data u(0, x) ~ "o(x), x ~ (xl, ...t -cd), by using the idea of the local relaxation approximation [2,3,6,10].The relaxing scheme is obtained in the following way: A linear hyperbolic systemwith a stiff source term is first constructed to approximate the original equation (1.1)with a small dissipative correction. Then this linear hyperbolic system is solved easilyby und…     

二维半线性反应扩散方程的交替方向隐格式

本文研究一类二维半线性反应扩散方程的差分方法.构造了一个二层线性化交替方向隐格式.利用离散能量估计方法证明了差分格式解的存在唯一性、差分格式在离散H~1模下的二阶收敛性和稳定性.最后给出两个数值例子验证了理论分析结果.     

A well-balanced van Leer-type numerical scheme for shallow water equations with variable topography

A well-balanced van Leer-type numerical scheme for the shallow water equations with variable topography is presented. The model involves a nonconservative term, which often makes standard schemes difficult to approximate solutions in certain regions. The construction of our scheme is based on exact solutions in computational form of local Riemann problems. Numerical tests are conducted, where comparisons between this van Leer-type scheme and a Godunov-type scheme are provided. Data for the tests are taken in both the subcritical region as well as supercritical region. Especially, tests for resonant cases where the exact solutions contain coinciding waves are also investigated. All numerical tests show that each of these two methods can give a good accuracy, while the van Leer -type scheme gives a better accuracy than the Godunov-type scheme. Furthermore, it is shown that the van Leer-type scheme is also well-balanced in the sense that it can capture exactly stationary contact discontinuity waves.