This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called...This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.展开更多
This paper deals with Mckean-Vlasov backward stochastic differential equations with weak monotonicity coefficients.We first establish the existence and uniqueness of solutions to Mckean-Vlasov backward stochastic diff...This paper deals with Mckean-Vlasov backward stochastic differential equations with weak monotonicity coefficients.We first establish the existence and uniqueness of solutions to Mckean-Vlasov backward stochastic differential equations.Then we obtain a comparison theorem in one-dimensional situation.展开更多
In this paper,we prove the existence of martingale solutions of a class of stochastic equations with a monotone drift of polynomial growth of arbitrary order and a continuous diffusion term with superlinear growth.Bot...In this paper,we prove the existence of martingale solutions of a class of stochastic equations with a monotone drift of polynomial growth of arbitrary order and a continuous diffusion term with superlinear growth.Both the nonlinear drift and diffusion terms are not required to be locally Lipschitz continuous.We then apply the abstract result to establish the existence of martingale solutions of the fractional stochastic reaction-diffusion equation with polynomial drift driven by a superlinear noise.The pseudo-monotonicity techniques and the Skorokhod-Jakubowski representation theorem in a topological space are used to pass to the limit of a sequence of approximate solutions defined by the Galerkin method.展开更多
A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue prob...A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).展开更多
The neutron diffusion equation plays a pivotal role in nuclear reactor analysis.Nevertheless,employing the physics-informed neural network(PINN)method for its solution entails certain limitations.Conventional PINN app...The neutron diffusion equation plays a pivotal role in nuclear reactor analysis.Nevertheless,employing the physics-informed neural network(PINN)method for its solution entails certain limitations.Conventional PINN approaches generally utilize a fully connected network(FCN)architecture that is susceptible to overfitting,training instability,and gradient vanishing as the network depth increases.These challenges result in accuracy bottlenecks in the solution.In response to these issues,the residual-based resample physics-informed neural network(R2-PINN)is proposed.It is an improved PINN architecture that replaces the FCN with a convolutional neural network with a shortcut(S-CNN).It incorporates skip connections to facilitate gradient propagation between network layers.Additionally,the incorporation of the residual adaptive resampling(RAR)mechanism dynamically increases the number of sampling points.This,in turn,enhances the spatial representation capabilities and overall predictive accuracy of the model.The experimental results illustrate that our approach significantly improves the convergence capability of the model and achieves high-precision predictions of the physical fields.Compared with conventional FCN-based PINN methods,R 2-PINN effectively overcomes the limitations inherent in current methods.Thus,it provides more accurate and robust solutions for neutron diffusion equations.展开更多
This paper is concerned with the stability of traveling wavefronts for a population dynamics model with time delay. Combining the weighted energy method and the comparison principle, the global exponential stability o...This paper is concerned with the stability of traveling wavefronts for a population dynamics model with time delay. Combining the weighted energy method and the comparison principle, the global exponential stability of noncritical traveling wavefronts （waves with speeds c 〉 c＊, where c=c＊ is the minimal speed） is established, when the initial perturbations around the wavefront decays to zero exponentially in space as x → -∞, but it can be allowed arbitrary large in other locations, which improves the results in[9, 18, 21].展开更多
A new approach, is established to show that the semigroup {S（t）≥0 generated by a reaction-diffusion equation with supercritical exponent is uniformly quasi-differentiable in L^q（Ω） （2 ≤ q 〈 ∞） with respect ...A new approach, is established to show that the semigroup {S（t）≥0 generated by a reaction-diffusion equation with supercritical exponent is uniformly quasi-differentiable in L^q（Ω） （2 ≤ q 〈 ∞） with respect to the initial value. As an application, this proves the upper-bound of fractal dimension for its global attractor in the corresponding space.展开更多
In this letter, a class of reaction-diffusion equations, which arise in chemical reaction or ecology and other fields of physics, are investigated. A more general analytical solution of the equation is obtained by usi...In this letter, a class of reaction-diffusion equations, which arise in chemical reaction or ecology and other fields of physics, are investigated. A more general analytical solution of the equation is obtained by using the first integral method.展开更多
In this paper, we prove the existence of random attractors for a stochastic reaction-diffusion equation with distribution derivatives on unbounded domains. The nonlinearity is dissipative for large values of the state...In this paper, we prove the existence of random attractors for a stochastic reaction-diffusion equation with distribution derivatives on unbounded domains. The nonlinearity is dissipative for large values of the state and the stochastic nature of the equation appears spatially distributed temporal white noise. The stochastic reaction-diffusion equation is recast as a continuous random dynamical system and asymptotic compactness for this demonstrated by using uniform estimates far-field values of solutions. The results are new and appear to be optimal.展开更多
Deep neural networks(DNNs)are effective in solving both forward and inverse problems for nonlinear partial differential equations(PDEs).However,conventional DNNs are not effective in handling problems such as delay di...Deep neural networks(DNNs)are effective in solving both forward and inverse problems for nonlinear partial differential equations(PDEs).However,conventional DNNs are not effective in handling problems such as delay differential equations(DDEs)and delay integrodifferential equations(DIDEs)with constant delays,primarily due to their low regularity at delayinduced breaking points.In this paper,a DNN method that combines multi-task learning(MTL)which is proposed to solve both the forward and inverse problems of DIDEs.The core idea of this approach is to divide the original equation into multiple tasks based on the delay,using auxiliary outputs to represent the integral terms,followed by the use of MTL to seamlessly incorporate the properties at the breaking points into the loss function.Furthermore,given the increased training dificulty associated with multiple tasks and outputs,we employ a sequential training scheme to reduce training complexity and provide reference solutions for subsequent tasks.This approach significantly enhances the approximation accuracy of solving DIDEs with DNNs,as demonstrated by comparisons with traditional DNN methods.We validate the effectiveness of this method through several numerical experiments,test various parameter sharing structures in MTL and compare the testing results of these structures.Finally,this method is implemented to solve the inverse problem of nonlinear DIDE and the results show that the unknown parameters of DIDE can be discovered with sparse or noisy data.展开更多
In this paper, we consider the existence of pullback random exponential attractor for non-autonomous random reaction-diffusion equation driven by nonlinear colored noise defined onR^(N) . The key steps of the proof ar...In this paper, we consider the existence of pullback random exponential attractor for non-autonomous random reaction-diffusion equation driven by nonlinear colored noise defined onR^(N) . The key steps of the proof are the tails estimate and to demonstrate the Lipschitz continuity and random squeezing property of the solution for the equation defined on R^(N) .展开更多
This paper investigates global solutions and long-time dynamics for the stochastic reaction-diffusion equation du=(Δu+f(u)+g(x,t))dt+σ(u)dW on a bounded domain,where the drift term f(u),with polynomial growth rate ...This paper investigates global solutions and long-time dynamics for the stochastic reaction-diffusion equation du=(Δu+f(u)+g(x,t))dt+σ(u)dW on a bounded domain,where the drift term f(u),with polynomial growth rate β,is strongly dissipative and the diffusion term σ(u)has growth rate γ,satisfying β+1>2γ.Under this condition,we establish the existence,uniqueness,and regularity of solutions in Bochner spaces.Our analysis relies only on weak monotonicity conditions and requires no further growth restrictions on f andσ.Moreover,we prove the existence of a weak mean random attractor for the system.These results offer new insights into the balance mechanism between stochastic perturbations and dissipative effects in superlinear regimes.展开更多
This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their...This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.展开更多
The methods of L-p estimation are used to discuss the extinction phenomena of the solutions to the following reaction-diffusion equations with initial-boundary values partial derivativeu/partial derivativet = Deltau -...The methods of L-p estimation are used to discuss the extinction phenomena of the solutions to the following reaction-diffusion equations with initial-boundary values partial derivativeu/partial derivativet = Deltau - lambda \u\(gamma-1) u - betau ((x, t) is an element of Omega x (0, + infinity)), u(x, t) \(partial derivativeOmegax (0, +infinity)) = 0, u(x, 0) = u(0) (x) is an element of H-0(1) (Omega) boolean AND L1+gamma(Omega) (x is an element of Omega). Sufficient and necessary conditions about the extinction of the solutions is given. Here lambda > 0, gamma > 0, beta > 0 are constants, Omega is an element of R-N is bounded with smooth boundary partial derivativeOmega. At last, it is simulated with a higher order equation by using the present methods.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
This paper deals with quasilinear elliptic equations of singular growth like-Δu-uΔ(u^(2))=a(x)u^(-1).We establish the existence of positive solutions for general a(x)∈L^(p)(Ω),p>2,whereΩis a bounded domain inℝ...This paper deals with quasilinear elliptic equations of singular growth like-Δu-uΔ(u^(2))=a(x)u^(-1).We establish the existence of positive solutions for general a(x)∈L^(p)(Ω),p>2,whereΩis a bounded domain inℝ^(N)with N≥1.展开更多
The initial boundary value problem of wave equations and reaction-diffusion equations with several nonlinear source terms in a bounded domain is studied by potential well method. The invarianee of some sets under the ...The initial boundary value problem of wave equations and reaction-diffusion equations with several nonlinear source terms in a bounded domain is studied by potential well method. The invarianee of some sets under the flow of these problems and the vacuum isolation of solutions are obtained by introducing a family of potential wells. Then the threshold result of global existence and nonexistence of solutions are given. Finally, the problem with critical initial conditions are discussed.展开更多
In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with ...In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with prescribed 2-norm has some normalized solutions by introducing variational methods.展开更多
In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,t...In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.展开更多
We study the Cauchy problem of the Kolmogorov-Fokker-Planck equations and show that the solution enjoys an analytic smoothing effect with L?initial datum for positive time.
文摘This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.
基金Supported by the National Natural Science Foundation of China(12001074)the Research Innovation Program of Graduate Students in Hunan Province(CX20220258)+1 种基金the Research Innovation Program of Graduate Students of Central South University(1053320214147)the Key Scientific Research Project of Higher Education Institutions in Henan Province(25B110025)。
文摘This paper deals with Mckean-Vlasov backward stochastic differential equations with weak monotonicity coefficients.We first establish the existence and uniqueness of solutions to Mckean-Vlasov backward stochastic differential equations.Then we obtain a comparison theorem in one-dimensional situation.
文摘In this paper,we prove the existence of martingale solutions of a class of stochastic equations with a monotone drift of polynomial growth of arbitrary order and a continuous diffusion term with superlinear growth.Both the nonlinear drift and diffusion terms are not required to be locally Lipschitz continuous.We then apply the abstract result to establish the existence of martingale solutions of the fractional stochastic reaction-diffusion equation with polynomial drift driven by a superlinear noise.The pseudo-monotonicity techniques and the Skorokhod-Jakubowski representation theorem in a topological space are used to pass to the limit of a sequence of approximate solutions defined by the Galerkin method.
基金supported by the National Natural Science Foundation of China(12401482)the second author was supported by the National Natural Science Foundation of China(12371371,12261160361,11971366)supported by the Open Research Fund of Hubei Key Laboratory of Computational Science,Wuhan University.
文摘A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).
基金supported by the Science and Technology on Reactor System Design Technology Laboratory(No.LRSDT12023108)supported in part by the Chongqing Postdoctoral Science Foundation(No.cstc2021jcyj-bsh0252)+2 种基金the National Natural Science Foundation of China(No.12005030)Sichuan Province to unveil the list of marshal industry common technology research projects(No.23jBGOV0001)Special Program for Stabilizing Support to Basic Research of National Basic Research Institutes(No.WDZC-2023-05-03-05).
文摘The neutron diffusion equation plays a pivotal role in nuclear reactor analysis.Nevertheless,employing the physics-informed neural network(PINN)method for its solution entails certain limitations.Conventional PINN approaches generally utilize a fully connected network(FCN)architecture that is susceptible to overfitting,training instability,and gradient vanishing as the network depth increases.These challenges result in accuracy bottlenecks in the solution.In response to these issues,the residual-based resample physics-informed neural network(R2-PINN)is proposed.It is an improved PINN architecture that replaces the FCN with a convolutional neural network with a shortcut(S-CNN).It incorporates skip connections to facilitate gradient propagation between network layers.Additionally,the incorporation of the residual adaptive resampling(RAR)mechanism dynamically increases the number of sampling points.This,in turn,enhances the spatial representation capabilities and overall predictive accuracy of the model.The experimental results illustrate that our approach significantly improves the convergence capability of the model and achieves high-precision predictions of the physical fields.Compared with conventional FCN-based PINN methods,R 2-PINN effectively overcomes the limitations inherent in current methods.Thus,it provides more accurate and robust solutions for neutron diffusion equations.
基金supported by NSF of China(11401478)Gansu Provincial Natural Science Foundation(145RJZA220)
文摘This paper is concerned with the stability of traveling wavefronts for a population dynamics model with time delay. Combining the weighted energy method and the comparison principle, the global exponential stability of noncritical traveling wavefronts （waves with speeds c 〉 c＊, where c=c＊ is the minimal speed） is established, when the initial perturbations around the wavefront decays to zero exponentially in space as x → -∞, but it can be allowed arbitrary large in other locations, which improves the results in[9, 18, 21].
基金Supported by NSFC Grant(11401100,10601021)the foundation of Fujian Education Department(JB14021)the innovation foundation of Fujian Normal University(IRTL1206)
文摘A new approach, is established to show that the semigroup {S（t）≥0 generated by a reaction-diffusion equation with supercritical exponent is uniformly quasi-differentiable in L^q（Ω） （2 ≤ q 〈 ∞） with respect to the initial value. As an application, this proves the upper-bound of fractal dimension for its global attractor in the corresponding space.
文摘In this letter, a class of reaction-diffusion equations, which arise in chemical reaction or ecology and other fields of physics, are investigated. A more general analytical solution of the equation is obtained by using the first integral method.
文摘In this paper, we prove the existence of random attractors for a stochastic reaction-diffusion equation with distribution derivatives on unbounded domains. The nonlinearity is dissipative for large values of the state and the stochastic nature of the equation appears spatially distributed temporal white noise. The stochastic reaction-diffusion equation is recast as a continuous random dynamical system and asymptotic compactness for this demonstrated by using uniform estimates far-field values of solutions. The results are new and appear to be optimal.
文摘Deep neural networks(DNNs)are effective in solving both forward and inverse problems for nonlinear partial differential equations(PDEs).However,conventional DNNs are not effective in handling problems such as delay differential equations(DDEs)and delay integrodifferential equations(DIDEs)with constant delays,primarily due to their low regularity at delayinduced breaking points.In this paper,a DNN method that combines multi-task learning(MTL)which is proposed to solve both the forward and inverse problems of DIDEs.The core idea of this approach is to divide the original equation into multiple tasks based on the delay,using auxiliary outputs to represent the integral terms,followed by the use of MTL to seamlessly incorporate the properties at the breaking points into the loss function.Furthermore,given the increased training dificulty associated with multiple tasks and outputs,we employ a sequential training scheme to reduce training complexity and provide reference solutions for subsequent tasks.This approach significantly enhances the approximation accuracy of solving DIDEs with DNNs,as demonstrated by comparisons with traditional DNN methods.We validate the effectiveness of this method through several numerical experiments,test various parameter sharing structures in MTL and compare the testing results of these structures.Finally,this method is implemented to solve the inverse problem of nonlinear DIDE and the results show that the unknown parameters of DIDE can be discovered with sparse or noisy data.
基金supported by the NSFC(12271141)supported by the Fundamental Research Funds for the Central Universities(B240205026)the Postgraduate Research and Practice Innovation Program of Jiangsu Province(KYCX24_0821).
文摘In this paper, we consider the existence of pullback random exponential attractor for non-autonomous random reaction-diffusion equation driven by nonlinear colored noise defined onR^(N) . The key steps of the proof are the tails estimate and to demonstrate the Lipschitz continuity and random squeezing property of the solution for the equation defined on R^(N) .
基金supported by the National Natural Science Foundation of China(No.12271399)the Fundamental Research Funds for the Central Universities(No.3122025090)。
文摘This paper investigates global solutions and long-time dynamics for the stochastic reaction-diffusion equation du=(Δu+f(u)+g(x,t))dt+σ(u)dW on a bounded domain,where the drift term f(u),with polynomial growth rate β,is strongly dissipative and the diffusion term σ(u)has growth rate γ,satisfying β+1>2γ.Under this condition,we establish the existence,uniqueness,and regularity of solutions in Bochner spaces.Our analysis relies only on weak monotonicity conditions and requires no further growth restrictions on f andσ.Moreover,we prove the existence of a weak mean random attractor for the system.These results offer new insights into the balance mechanism between stochastic perturbations and dissipative effects in superlinear regimes.
文摘This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.
文摘The methods of L-p estimation are used to discuss the extinction phenomena of the solutions to the following reaction-diffusion equations with initial-boundary values partial derivativeu/partial derivativet = Deltau - lambda \u\(gamma-1) u - betau ((x, t) is an element of Omega x (0, + infinity)), u(x, t) \(partial derivativeOmegax (0, +infinity)) = 0, u(x, 0) = u(0) (x) is an element of H-0(1) (Omega) boolean AND L1+gamma(Omega) (x is an element of Omega). Sufficient and necessary conditions about the extinction of the solutions is given. Here lambda > 0, gamma > 0, beta > 0 are constants, Omega is an element of R-N is bounded with smooth boundary partial derivativeOmega. At last, it is simulated with a higher order equation by using the present methods.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金Supported by National Science Foundation of China(11971027,12171497)。
文摘This paper deals with quasilinear elliptic equations of singular growth like-Δu-uΔ(u^(2))=a(x)u^(-1).We establish the existence of positive solutions for general a(x)∈L^(p)(Ω),p>2,whereΩis a bounded domain inℝ^(N)with N≥1.
基金the National Natural Science Foundation of China(No.10271034)the Basic Research Foundation of Harbin Engineering University(No.HEUF04012)
文摘The initial boundary value problem of wave equations and reaction-diffusion equations with several nonlinear source terms in a bounded domain is studied by potential well method. The invarianee of some sets under the flow of these problems and the vacuum isolation of solutions are obtained by introducing a family of potential wells. Then the threshold result of global existence and nonexistence of solutions are given. Finally, the problem with critical initial conditions are discussed.
基金Supported by the National Natural Science Foundation of China(11671403,11671236,12101192)Henan Provincial General Natural Science Foundation Project(232300420113)。
文摘In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with prescribed 2-norm has some normalized solutions by introducing variational methods.
基金Supported by Research Project Supported by Shanxi Scholarship Council of China(2021-029)International Cooperation Base and Platform Project of Shanxi Province(202104041101019)+2 种基金Basic Research Plan of Shanxi Province(202203021211129)Shanxi Province Natural Science Research(202203021212249)Special/Youth Foundation of Taiyuan University of Technology(2022QN101)。
文摘In this paper,we construct two fully decoupled,second-order semi-discrete numerical schemes for the Boussinesq equations based on the scalar auxiliary variable(SAV)approach.By introducing a scalar auxiliary variable,the original Boussinesq system is transformed into an equivalent one.Then we discretize it using the second-order backward di erentiation formula(BDF2)and Crank-Nicolson(CN)to obtain two second-order time-advanced schemes.In both numerical schemes,a pressure-correction method is employed to decouple the velocity and pressure.These two schemes possess the desired property that they can be fully decoupled with satisfying unconditional stability.We rigorously prove both the unconditional stability and unique solvability of the discrete schemes.Furthermore,we provide detailed implementations of the decoupling procedures.Finally,various 2D numerical simulations are performed to verify the accuracy and energy stability of the proposed schemes.
基金Supported by NSFC (No.12031006)Fundamental Research Funds for the Central Universities of China。
文摘We study the Cauchy problem of the Kolmogorov-Fokker-Planck equations and show that the solution enjoys an analytic smoothing effect with L?initial datum for positive time.
