高阶Schroedinger型方程的两层高精度恒稳差分格式

众所周知，高阶Schroedinger方程在量子力学、非线性光学及流体力学中都有广泛的应用。本文对高阶Schroedinger型方程δu/δt=i(-1)^mδ2m/δx^2m(其中i=√-1，m为正整数)，利用待定系数法，构造出一个两层高精度的隐式差分格式。其截断误差阶为O((△t)^2 (Δx)^6)，比同类格式精度高2～4阶，并用Fourier分析法证明了它是绝对稳定的。最后，数值例子表明本文格式比著名的Crank-Nicolson格式精度高10^-2～10^-7，这说明我们的格式是有效的，理论分析与实际计算相吻合。     

高阶Schr(o)dinger型方程的两层高精度恒稳差分格式

众所周知, 高阶Schrodinger方程在量子力学、非线性光学及流体力学中都有广泛的应用.本文对高阶Schrodinger型方程(eu/et)＝I(-1)m(e2mu)/(ex2m)(其中I＝-1,m为正整数),利用待定系数法,构造出一个两层高精度的隐式差分格式.其截断误差阶为O((Δt)2+(Δx)6),比同类格式精度高2～4阶,并用Fourier分析法证明了它是绝对稳定的.最后,数值例子表明本文格式比著名的Crank-Nicolson格式精度高10-2～10-7,这说明我们的格式是有效的,理论分析与实际计算相吻合.     

污染扩散方程高精度紧致格式及其稳定性分析

针对污染扩散方程提出了时间任意阶精度的显式格式,并对该格式的稳定性和精度进行了分析,理论结果表明:一阶精度的计算格式是传统的显格式,其稳定条件为:s≤1/2(s=D.Δt/Δx2,D为扩散系数,Δt为时间步长,Δx为空间步长),随着保留精度阶数的增加,稳定性范围也会随之增大;当保留无穷阶精度时,格式是无条件稳定的。这也就从一个侧面揭示了稳定性与时间精度之间的关系,为高性能数值计算格式的构思提供了可以借鉴的原则。数值算例的结果表明,本文格式具有一定的实用性。     

CONCENTRATION OF COUPLED CUBIC NONLINEAR SCHRDINGER EQUATIONS

IntroductionThe following coupled cubic nonlinear Schr dinger equations inR2are considered:it-Δ=(||2-|ψ|2),t>0,x∈R2,iψt-Δψ=(-||2 |ψ|2)ψ,t>0,x∈R2,(1)where=(t,x):[0,T)×R2→C,ψ=ψ(t,x):[0,T)×R2→C,Δis the Laplaceoperator onR2.The e     

Asymptotic theory of initial value problems for nonlinear perturbed Klein-Gordon equations

Introduction InthispaperasymptotictheoryofthefollowinginitialvalueproblemforanonlinearKlein Gordonequationisconsidered.tt-Δ =εF(t,x,,ε),t>0,x∈R2,(0,x,ε)=0(x,ε),t(0,x,ε)=1(x,ε),x∈R2,(1)where(t,x)isarealvaluedunknownfunction,Δ=2i     

求解浅水波方程的半离散中心迎风方法

将Jin's的界面方法应用到求解双曲守恒型方程的半离散中心迎风方法中,给出了一种新的求解浅水波方程的半离散中心迎风差分方法。对于源项,不是采用传统的单元均值而是采用单元界面处的值来近似,使所得格式对稳定态的求解是均衡的。且已证明所给的二阶精度的求解格式保持水深的非负性,这一特性使其能够较好的处理干河床问题。使用该方法产生的数值粘性(与O(Δ2r-1)同阶)要比交错的中心格式小(与O(Δx2r/Δt)同阶),而且由于数值粘性与时间步长无关,从而时间步长可根据稳定性需要尽可能的小,因此适用于稳定态的求解。     

NEWTON METHOD FOR SOLVING A CLASS OF SMOOTH CONVEX PROGRAMMING

1ConvexProgrammingandExactPenaltyFunction Weconsiderthefollowingconvexprogramming:(P)minf(x)s.t.x∈S={x∈Rn:gi(x)≤0,i=1,…,m}.SupposethatSisacompactset.ThusthereexistsalargeboundedboxX,suchthatS={x∈Rn:gi(x)≤0,i=1,…,m}intX.Assumethatf(x),gi(x),i=1,…,m     

On the existence of periodic solutions for nonlinear system with multiple delays

IntroductionInthispaper,westudyT_periodicsolutionsofthefollowingnonlinearsystemwithmultipledelays x(t) =f(t,x(t) ,x(t-τ1(t) ) ,… ,x(t -τm(t) ) ) ,(1 )wherex(t) ∈C(R ,R) ,fiscontinuous,f(t+T ,·) =f(t,·) ,τi(t) (i=1 ,2 ,… ,m)arecontinuousperiodicfunctionsofperiodT .AlemmaisintroducedfordiscussingtheexistenceofT_periodicsolutionofsystem (1 ) .LetXbeaBanachSpace ,considerthefollowingoperatorequation :Lx =λNx　　 (λ∈ [0 ,1 ] ) ,whereL :DomL∩X→Xisalinearoperator,λ∈ [0 ,1 ]isapa…     

SINGULARLY PERTURBED BOUNDARY VALUE PROBLEMS FOR SEMI-LINEAR RETARDED DIFFERENTIAL EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS

IntroductionInthispaper,westudiedakindofboundaryvalueproblems (BVPs)forsemi_linearretardeddifferentialequationwithnonlinearboundarycondition :　　　　εx″(t) =f(t,x(t) ,x(t-ε) ,ε) ,　　t∈(0 ,1 ) ,(1 )　　　　x(t) =φ(t,ε) ,　t∈[-ε0 ,0 ] ,h(x(1 ) ,x′(1 ) ,ε) =A(ε) ,(2 )whereε>0isasmallparameterandε0 isasufficientlysmallpositiveconstant.ThereweremanyresultsofstudyingonsingularlyperturbedboundaryvalueproblemforretardeddifferentialequationinRefs.[1～5] .Butthosestudiespossessedanesse…     

Difference scheme and numerical simulation based on mixed finite element method for natural convection problem

ntroductionLetΩ R2 beaboundeddomain .Weconsiderthefollowingnon_stationarynaturalconvectionproblem :Problem (Ⅰ )　Findu =(u1,u2 ) ,p ,andTsuchthat,foranyt1>0 ,ut- μΔu +(u· )u + p=λjT　　 ((x ,y ,t) ∈Ω× (0 ,t1) ) ,divu =0　　　　　　　　　 ((x ,y,t) ∈Ω× (0 ,t1) ) ,Tt-ΔT +λu· T =0　　 ((x,y,t) ∈Ω× (0 ,t1) ) ,u =0 ,T =0　　　　　　 ((x,y,t)∈ Ω× (0 ,t1) ) ,u(x ,y ,0 ) =0 ,　T(x,y,0 ) =f(x,y)　　 ((x,y) ∈Ω) ,whereuisthefluidvelocityvectorfield ,pthepressurefield ,Tthet…     

DYNAMICAL CHARACTER FOR A PERTURBED COUPLED NONLINEAR SCHRDINGER SYSTEM

Introduction WeconsidertheperturbedcouplednonlinearSchr dingerequations(CNLS)iq1t=q1xx 2[|q1|2 |q2|2-ω2]q1 iε(^D1q1-r1),iq2t=q2xx 2[|q1|2 |q2|2-ω2]q2 iε(^D2q2-r2),(1)whereqjis2πperiodicandeveninx.rjisconstant,^Djisboundeddissipativeoperatorand assumedtotaketheform^Djqj=-αjqj βj^Bqj,(2)αjandβjarepositiveconstants,j=1,2.^BisaFouriertruncationofLaplaceoperator xx,i.e.,^Bcos(kx)=-k2cos(kx),k0isasmallperturbation parameter.Th…     

CONTROLLABILITY OF A CLASS OF HYBRID DYNAMIC SYSTEMS (Ⅱ)—SINGLE TIME_DELAY CASE

IntroductionConsideraswitchedlinearsystemwithsingletime_delayinthecontrolfunctiongivenby x(t) =Ar(t) x(t) +Br(t) u(t) +Dr(t) u(t -τ) ,(1)wherex(t) ∈Rnisthestate,u(t) ∈R mtheinput,thepiecewiseconstantscalarfunctionr(t) :R+ → 1,2 ,… ,Nistheswitchingpath ,(Ai,Bi,Di) ｜i =1,… ,Nisafinitefamilyofsystemrealizations.Moreover,r(t) =iimpliesthat(Ai,Bi,Di)ischosenasthesystemrealizationattimet.τ >0isthefixedtimedelayincontrol.Inthispaper ,wewillstudythecontrollabilityofsystem (1) .Theone_per…     

CONTROLLABILITY OF A CLASS OF HYBRID DYNAMIC SYSTEMS (Ⅲ)—MULTIPLE TIME_DELAY CASE

IntroductionConsidertheswitchedlinearsystemwithmultipletime_delaysincontrolfunctiongivenby x(t) =Ar(t) x(t) +Br(t) u(t) + Kk=1Dr(t) ,ku(t -τk) ,(1)wherex(t) ,u(t)andr(t)isdefinedasbefore ,(Ai,Bi,Di,1,… ,Di,K)｜i=1,… ,Nisafinitefamilyofsystemrealizations.Moreover,r(t) =iimplies (Ai,Bi,Di,1,… ,Di,K)ischosenasthesystemrealizationattimet.K <∞isthenumberoftime_delaysinthesystem .0<τ1<… <τKareKfixedtime_delays .Remark 1　Similarly ,wecandescribeswitchingpathbyintroducingswitchingsequ…     

On the problem of preventing blowing-up and quenching for semilinear heat equation

In this paper,the global existence of solutions to the IVP=Δu+g(t)f(u) (t>0),u|_(t=0)=u_0(x)and the (?)PVPu_t=Δu-g(t,x)f(u)(t>0,x∈Ω),u|_(t-0)=u|_(?)(?)is investigated. As has been done in [6]the (?)duction of factor g(t) or g(t.x) innonlinear term is to prevent(?) occurrance of blowing-up or quenching of solutions.It isshown in this paper that most of the restrictions on f,g and u_0 in the theorems of[6] maybe cancelled or relaxed,that the smallness of g is required only for t large,and thatunder certain conditions controlling initial state can avoid blowing-up.     

非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

On the boundedness and the stability properties of solution of certain fourth order differential equations

This paper investigates equation(1)in two cases:(i)P≡0,(ii)P(≠O)satisfies|P(t,x,y,z,ω)|≤(A |y| |z| |ω|)q(t),where q(t)is a nonnegative function of t.For case(i)the asymptotic stability in the large of the trivial solution x=0 is investigatedand for case(ii)the boundedness result is obtained for solutions of equation(1).Theseresults improve and include several well-known results.     

非稳态Hamilton-Jacobi方程的7阶加权紧致非线性格式

Hamilton-Jacobi (HJ) 方程是一类重要的非线性偏微分方程, 在物理学、流体力学、图像处理、微分几何、金融数学、最优化控制理论等方面有着广泛的应用. 由于HJ方程的弱解存在但不唯一, 且解的导数可能出现间断, 导致其数值求解具有一定的难度. 本文提出了非稳态HJ方程的7阶精度加权紧致非线性格式 (WCNS). 该格式结合了Hamilton函数的Lax-Friedrichs型通量分裂方法和一阶空间导数左、右极限值的高阶精度混合节点和半节点型中心差分格式. 基于7点全局模板和4个4点子模板推导了半节点函数值的高阶线性逼近和4个低阶线性逼近, 以及全局模板和子模板的光滑度量指标. 为避免间断附近数值解产生非物理振荡以及提高格式稳定性, 采用WENO型非线性插值方法计算半节点函数值. 时间离散采用3阶TVD型Runge-Kutta方法. 通过理论分析验证了WCNS格式对于光滑解具有最佳的7阶精度. 为方便比较, 经典的7阶WENO格式也被推广用于求解HJ方程. 数值结果表明, 本文提出的WCNS格式能够很好地模拟HJ方程的精确解, 且在光滑区域能够达到7阶精度; 与经典的同阶WENO格式相比, WCNS格式在精度、收敛性和分辨率方面更优, 计算效率略高.      

Analysis of chebyshev pseudospectral method for multi-dimensional generalized srlw equations

IntroductionAsymmetricregularizedlongwaveequation (SRLWE) 2 x2 -1 u t = x ρ+ 12 u2 ,　 ρ t+ u x=0 (1 )hasbeeninvestigatedinRef.[1 ] .Thesystem (1 )ofequationsisshowntodescribeweaklynonlinearion_acousticwaveandspace_chargewaves.Thehuperbolicsecantsquaredsolitarywaves ,thefourconservationlaws,andsomenumericalresultshavebeenobtainedinRef.[1 ] .Obviously ,eliminatingρin (1 ) ,weobtainaclassofregularizedlongwaveequationutt-uxx+ 12 u2xt-uxxtt =0 . (2 )TheSRLWequationisexplicitlysymmetric…     

Global linear and quadratic one-step smoothing newton method for vertical linear complementarity problems

IntroductionLetN ∈Rm0 ×n(m0 ≥n)beaverticalblockmatrixoftype(m1,… ,mn) ,q∈Rm0 beaconstantvectorpartitionedconformablywithN ,thatareN =N1Nn,　q =q1qn,　m0 = ni=1mi,whereNi ∈Rmi×n,qi ∈Rmi,i ∈I =1 ,… ,n .WeconsidertheverticallinearomplementarityproblemVLCP(N ,q)associatedwith (N ,q)offindingavectorx∈Rnsuchthatx≥ 0 ,　si(x) :=Nix+qi ≥ 0 ,　xi∏mij=1sij(x) =0　　 (i∈I) ,wherexi,sij(x)denotetheithcomponentofxandthejthcomponentofsi(x) ,respectively ;∏ mj=1ajdenotesa1…am.Iti…     

非结构网格上p型多重网格法流场数值模拟

在二维、三维非结构网榕上,针对间断Galerkin方法计算量大、收敛慢的缺点将p型多重网格方法应用于该方法求解跨音速Euler方程,提高计算效率。p型多重网格方法是通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛。文中对高阶近似(p＞0)使用显式格式,最低阶近似(p=0)采用隐式格式。NACA0012翼型和O...