The equation of pattern formation induced by buoyancy or by surface-tension gradient in finite systems confined between horizontal poor heat conductors is introduced by Knobloch[1990] where u is the planform function,...The equation of pattern formation induced by buoyancy or by surface-tension gradient in finite systems confined between horizontal poor heat conductors is introduced by Knobloch[1990] where u is the planform function, μ is the scaled Rayleigh number, K = 1 and α represents the effects of a heat transfer finite Blot number. The cofficients β, δ and γ do not vanish when the boundary, conditions at top and bottom are not identical (β / 0, δ / 0) or nonBoussinesq effects are taked into account (γ / 0). In this paper, the Knobloch equation with α > 0 is considered, the global existence in L2-space and the finite existence time of solution in V2-space have been obtained respectively.展开更多
In the present paper,the local existence of classical solutions to the periodic boundary problem and the Cauchy problem of a quasilinear evolution equation are studied under the assumptions that do not require the mon...In the present paper,the local existence of classical solutions to the periodic boundary problem and the Cauchy problem of a quasilinear evolution equation are studied under the assumptions that do not require the monotonicity of σi(s) (i= 1,…, n). The nonexistence of global solutions to the initial-boundary value problem of the equation is also discussed, a blowup theorem is proved and a concrete example is given.展开更多
In this paper, we extend the reliable modification of the Adomian Decom-position Method coupled to the Lesnic’s approach to solve boundary value problems and initial boundary value problems with mixed boundary condit...In this paper, we extend the reliable modification of the Adomian Decom-position Method coupled to the Lesnic’s approach to solve boundary value problems and initial boundary value problems with mixed boundary conditions for linear and nonlinear partial differential equations. The method is applied to different forms of heat and wave equations as illustrative examples to exhibit the effectiveness of the method. The method provides the solution in a rapidly convergent series with components that can be computed iteratively. The numerical results for the illustrative examples obtained show remarkable agreement with the exact solutions. We also provide some graphical representations for clear-cut comparisons between the solutions using Maple software.展开更多
This paper is concerned with the initial-boundary value problem of scalar conservation laws with weak discontinuous flux, whose initial data are a function with two pieces of constant and whose boundary data are a con...This paper is concerned with the initial-boundary value problem of scalar conservation laws with weak discontinuous flux, whose initial data are a function with two pieces of constant and whose boundary data are a constant function. Under the condition that the flux function has a finite number of weak discontinuous points, by using the structure of weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux-Nedelec, we give a construction method to the global weak entropy solution for this initial-boundary value problem, and by investigating the interaction of elementary waves and the boundary, we clarify the geometric structure and the behavior of boundary for the weak entropy solution.展开更多
An initial boundary-value problem for the Hirota equation on the half-line,0<x<∞, t>0, is analysed by expressing the solution q(x, t) in terms of the solution of a matrix Riemann-Hilbert(RH) problem in the c...An initial boundary-value problem for the Hirota equation on the half-line,0<x<∞, t>0, is analysed by expressing the solution q(x, t) in terms of the solution of a matrix Riemann-Hilbert(RH) problem in the complex k-plane. This RH problem has explicit(x, t) dependence and it involves certain functions of k referred to as the spectral functions. Some of these functions are defined in terms of the initial condition q(x,0) = q0(x), while the remaining spectral functions are defined in terms of the boundary values q(0, t) = g0(t), qx(0, t) = g1(t) and qxx(0, t) = g2(t). The spectral functions satisfy an algebraic global relation which characterizes, say, g2(t) in terms of {q0(x), g0(t), g1(t)}.The spectral functions are not independent, but related by a compatibility condition, the so-called global relation.展开更多
The object of this work is to investigate the initial-boundary value problem for coupled Hirota equation on the half-line.We show that the solution of the coupled Hirota equation can be expressed in terms of the solut...The object of this work is to investigate the initial-boundary value problem for coupled Hirota equation on the half-line.We show that the solution of the coupled Hirota equation can be expressed in terms of the solution of a 3×3 matrix Riemann-Hilbert problem formulated in the complex k-plane.The relevant jump matrices are explicitly given in terms of the matrix-valued spectral functions s(k)and S(k)that depend on the initial data and boundary values,respectively.Then,applying nonlinear steepest descent techniques to the associated 3×3 matrix-valued Riemann-Hilbert problem,we can give the precise leading-order asymptotic formulas and uniform error estimates for the solution of the coupled Hirota equation.展开更多
In the hyperbolic research community,there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability.In the first par...In the hyperbolic research community,there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability.In the first part of the series[6],the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly.By applying this technique,the authors demonstrate that a pure continu-ous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way.In this work,we extend this investigation to the nonlinear case and focus on entropy conservation.By switching to entropy variables,we provide an estimation of the boundary operators also for nonlinear problems,that guarantee conservation.In numerical simulations,we verify our theoretical analysis.展开更多
A Legendre-Legendre spectral collocation scheme is constructed for Korteweg-de Vries(KdV)equation on bounded domain by using the Legendre collocation method in both time and space,which is a nonlinear matrix equation ...A Legendre-Legendre spectral collocation scheme is constructed for Korteweg-de Vries(KdV)equation on bounded domain by using the Legendre collocation method in both time and space,which is a nonlinear matrix equation that is changed to a nonlinear systems and can be solved by the usual fixed point iteration.Numerical results demonstrate the efficiency of the method and spectral accuracy.展开更多
Using the Fokas unified method, we consider the initial boundary value problem for the Fokas-Lenells equation on the finite interval. We present that the Neumann boundary data can be explicitly expressed by Dirichlet ...Using the Fokas unified method, we consider the initial boundary value problem for the Fokas-Lenells equation on the finite interval. We present that the Neumann boundary data can be explicitly expressed by Dirichlet boundary conditions prescribed, and extend the idea of the linearizable boundary conditions for equations on the half line to Fokas-Lenells equation on the finite interval.展开更多
The Gelfand-Levitan-Marchenko representation is used to analyze the initialboundary value problem of two-component nonlinear Schr¨odinger equation on the half-line.It has shown that the global relation can be eff...The Gelfand-Levitan-Marchenko representation is used to analyze the initialboundary value problem of two-component nonlinear Schr¨odinger equation on the half-line.It has shown that the global relation can be effectively analyzed by the Gelfand-LevitanMarchenko representation. we also derive expressions for the Dirichlet-to-Neumann map to characterize the unknown boundary values.展开更多
The Method Of Lines (MOL) and Scheifele’s G-functions in the design of algorithms adapted for the numeric integration of parabolic Partial Differential Equations (PDE) in one space dimension are applied. The semi-dis...The Method Of Lines (MOL) and Scheifele’s G-functions in the design of algorithms adapted for the numeric integration of parabolic Partial Differential Equations (PDE) in one space dimension are applied. The semi-discrete system of ordinary differential equations in the time direction, obtained by applying the MOL to PDE, is solved with the use of a method of Adapted Series, based on Scheifele’s G-functions. This method integrates exactly unperturbed linear systems of ordinary differential equations, with only one G-function. An implementation of this algorithm is used to approximate the solution of two test problems proposed by various authors. The results obtained by the Dufort-Frankel, Crank-Nicholson and the methods of Adapted Series versus the analytical solution, show the results of mistakes made.展开更多
This paper is concerned with the initial-boundary value problem of a nonlinear conservation law in the half space R+= {x |x > 0} where a>0 , u(x,t) is an unknown function of x ∈ R+ and t>0 , u ± , um ar...This paper is concerned with the initial-boundary value problem of a nonlinear conservation law in the half space R+= {x |x > 0} where a>0 , u(x,t) is an unknown function of x ∈ R+ and t>0 , u ± , um are three given constants satisfying um=u+≠u- or um=u-≠u+ , and the flux function f is a given continuous function with a weak discontinuous point ud. The main purpose of our present manuscript is devoted to studying the structure of the global weak entropy solution for the above initial-boundary value problem under the condition of f '-(ud) > f '+(ud). By the characteristic method and the truncation method, we construct the global weak entropy solution of this initial-boundary value problem, and investigate the interaction of elementary waves with the boundary and the boundary behavior of the weak entropy solution.展开更多
By using fluid dynamics theory with the effects of adsorption and reaction, the chromatography model with a reaction A →B was established as a system of two hyperbolic partial differential equations (PDE’s)....By using fluid dynamics theory with the effects of adsorption and reaction, the chromatography model with a reaction A →B was established as a system of two hyperbolic partial differential equations (PDE’s). In some practical situations, the reaction chromatography model was simplified a semi-coupled system of two linear hyperbolic PDE’s. In which, the reactant concentration wave model was the initial-boundary value problem of a self-closed hyperbolic PDE, while the resultant concentration wave model was the initial-boundary value problem of hyperbolic PDE coupling reactant concentration. The general explicit expressions for the concentration wave of the reactants and resultants were derived by Laplace transform. The δ-pulse and wide pulse injections were taken as the examples to discuss detailedly, and then the stability analysis between the resultant solutions of the two modes of pulse injection was further discussed. It was significant for further analysis of chromatography, optimizing chromatographic separation, determining the physical and chemical characters.展开更多
The ideal reaction chromatography model can be regarded as a semi-coupled system of two hyperbolic partial differential equations, in which, one is a self-closed nonlinear equation for the reactant concentration and a...The ideal reaction chromatography model can be regarded as a semi-coupled system of two hyperbolic partial differential equations, in which, one is a self-closed nonlinear equation for the reactant concentration and another is a linear equation coupling the reactant concentration for the resultant concentration. This paper is concerned with the initial-boundary value problem for the above model. By the characteristic method and the truncation method, we construct the global weak entropy solution of this initial initial-boundary value problem for Riemann type of initial-boundary data. Moreover, as examples, we apply the obtained results to the cases of head-on and wide pulse injections and give the expression of the global weak entropy solution.展开更多
We study an finite-difference time-domain (FDTD) system of uniaxial perfectly matched layer (UPML) method for electromagnetic scattering problems. Particularly we analyze the discrete initial-boundary value problems o...We study an finite-difference time-domain (FDTD) system of uniaxial perfectly matched layer (UPML) method for electromagnetic scattering problems. Particularly we analyze the discrete initial-boundary value problems of the transverse magnetic mode (TM) to Maxwell's equations with Yee's algorithm. An exterior domain in two spacial dimension is truncated by a square with a perfectly matched layer filled by a certain artificial medium. Besides, an artificial boundary condition is imposed on the outer boundary of the UPML. Using energy method, we obtain the stability of this FDTD system on the truncated domain. Numerical experiments are designed to approve the theoretical analysis.展开更多
In this paper, the Fokas unified method is used to analyze the initial-boundary value problem of a complex Sharma–Tasso–Olver(c STO) equation on the half line. We show that the solution can be expressed in terms of ...In this paper, the Fokas unified method is used to analyze the initial-boundary value problem of a complex Sharma–Tasso–Olver(c STO) equation on the half line. We show that the solution can be expressed in terms of the solution of a Riemann–Hilbert problem. The relevant jump matrices are explicitly given in terms of the matrix-value spectral functions spectral functions {a(λ), b(λ)} and {A(λ), B(λ)}, which depending on initial data u_0(x) = u(x, 0) and boundary data g_0(y) = u(0, y), g_1(y) = ux(0, y), g_2(y) = u_(xx)(0, y). These spectral functions are not independent, they satisfy a global relation.展开更多
In this paper we construct the suitable weak solution for the initial-boundary valueproblem of the Boussinesq equations and obtain some properties for these solutions.Also in the caseof a two dimensional space the uni...In this paper we construct the suitable weak solution for the initial-boundary valueproblem of the Boussinesq equations and obtain some properties for these solutions.Also in the caseof a two dimensional space the uniqueness of weak solution is proved.展开更多
Initial-boundary value problems for a class of systems of semilinear parabolicequations with singular terms are investigated, conditions for the existence and nonexistence of the global solutions to the above problems...Initial-boundary value problems for a class of systems of semilinear parabolicequations with singular terms are investigated, conditions for the existence and nonexistence of the global solutions to the above problems are given, and upper bounds of thequenching time and the critical length are obtained.展开更多
This paper is concerned with the asymptotic behavior of the solution to the Euler equations with time-depending damping on quadrant(x,t)∈R^+×R^+,with the null-Dirichlet boundary condition or the null-Neumann bou...This paper is concerned with the asymptotic behavior of the solution to the Euler equations with time-depending damping on quadrant(x,t)∈R^+×R^+,with the null-Dirichlet boundary condition or the null-Neumann boundary condition on u. We show that the corresponding initial-boundary value problem admits a unique global smooth solution which tends timeasymptotically to the nonlinear diffusion wave. Compared with the previous work about Euler equations with constant coefficient damping, studied by Nishihara and Yang(1999), and Jiang and Zhu(2009, Discrete Contin Dyn Syst), we obtain a general result when the initial perturbation belongs to the same space. In addition,our main novelty lies in the fact that the cut-off points of the convergence rates are different from our previous result about the Cauchy problem. Our proof is based on the classical energy method and the analyses of the nonlinear diffusion wave.展开更多
基金Project supported by the National Natural Science Foundation of China!(No:19861004)
文摘The equation of pattern formation induced by buoyancy or by surface-tension gradient in finite systems confined between horizontal poor heat conductors is introduced by Knobloch[1990] where u is the planform function, μ is the scaled Rayleigh number, K = 1 and α represents the effects of a heat transfer finite Blot number. The cofficients β, δ and γ do not vanish when the boundary, conditions at top and bottom are not identical (β / 0, δ / 0) or nonBoussinesq effects are taked into account (γ / 0). In this paper, the Knobloch equation with α > 0 is considered, the global existence in L2-space and the finite existence time of solution in V2-space have been obtained respectively.
基金Natural Science Foundation of Henan Province!(Grant No.98405070) National Natural Science Foundation of China (Grant No.19
文摘In the present paper,the local existence of classical solutions to the periodic boundary problem and the Cauchy problem of a quasilinear evolution equation are studied under the assumptions that do not require the monotonicity of σi(s) (i= 1,…, n). The nonexistence of global solutions to the initial-boundary value problem of the equation is also discussed, a blowup theorem is proved and a concrete example is given.
文摘In this paper, we extend the reliable modification of the Adomian Decom-position Method coupled to the Lesnic’s approach to solve boundary value problems and initial boundary value problems with mixed boundary conditions for linear and nonlinear partial differential equations. The method is applied to different forms of heat and wave equations as illustrative examples to exhibit the effectiveness of the method. The method provides the solution in a rapidly convergent series with components that can be computed iteratively. The numerical results for the illustrative examples obtained show remarkable agreement with the exact solutions. We also provide some graphical representations for clear-cut comparisons between the solutions using Maple software.
文摘This paper is concerned with the initial-boundary value problem of scalar conservation laws with weak discontinuous flux, whose initial data are a function with two pieces of constant and whose boundary data are a constant function. Under the condition that the flux function has a finite number of weak discontinuous points, by using the structure of weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux-Nedelec, we give a construction method to the global weak entropy solution for this initial-boundary value problem, and by investigating the interaction of elementary waves and the boundary, we clarify the geometric structure and the behavior of boundary for the weak entropy solution.
基金supported by the China Postdoctoral Science Foundation(No.2015M580285).
文摘An initial boundary-value problem for the Hirota equation on the half-line,0<x<∞, t>0, is analysed by expressing the solution q(x, t) in terms of the solution of a matrix Riemann-Hilbert(RH) problem in the complex k-plane. This RH problem has explicit(x, t) dependence and it involves certain functions of k referred to as the spectral functions. Some of these functions are defined in terms of the initial condition q(x,0) = q0(x), while the remaining spectral functions are defined in terms of the boundary values q(0, t) = g0(t), qx(0, t) = g1(t) and qxx(0, t) = g2(t). The spectral functions satisfy an algebraic global relation which characterizes, say, g2(t) in terms of {q0(x), g0(t), g1(t)}.The spectral functions are not independent, but related by a compatibility condition, the so-called global relation.
基金supported by the China Postdoctoral Science Foundation(Grant No.2019TQ0041)。
文摘The object of this work is to investigate the initial-boundary value problem for coupled Hirota equation on the half-line.We show that the solution of the coupled Hirota equation can be expressed in terms of the solution of a 3×3 matrix Riemann-Hilbert problem formulated in the complex k-plane.The relevant jump matrices are explicitly given in terms of the matrix-valued spectral functions s(k)and S(k)that depend on the initial data and boundary values,respectively.Then,applying nonlinear steepest descent techniques to the associated 3×3 matrix-valued Riemann-Hilbert problem,we can give the precise leading-order asymptotic formulas and uniform error estimates for the solution of the coupled Hirota equation.
基金funded by the SNF Grant(Number 200021175784)the UZH Postdoc grant+1 种基金funded by an SNF Grant 200021_153604The Los Alamos unlimited release number is LA-UR-19-32411.
文摘In the hyperbolic research community,there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability.In the first part of the series[6],the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly.By applying this technique,the authors demonstrate that a pure continu-ous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way.In this work,we extend this investigation to the nonlinear case and focus on entropy conservation.By switching to entropy variables,we provide an estimation of the boundary operators also for nonlinear problems,that guarantee conservation.In numerical simulations,we verify our theoretical analysis.
基金Supported by National Natural Science Foundation of China(Grant Nos.11771299,11371123)Natural Science Foundation of Henan Province(Grant No.202300410156).
文摘A Legendre-Legendre spectral collocation scheme is constructed for Korteweg-de Vries(KdV)equation on bounded domain by using the Legendre collocation method in both time and space,which is a nonlinear matrix equation that is changed to a nonlinear systems and can be solved by the usual fixed point iteration.Numerical results demonstrate the efficiency of the method and spectral accuracy.
基金supported by grants from the National Natural Science Foundation of China(11271079,11626090)
文摘Using the Fokas unified method, we consider the initial boundary value problem for the Fokas-Lenells equation on the finite interval. We present that the Neumann boundary data can be explicitly expressed by Dirichlet boundary conditions prescribed, and extend the idea of the linearizable boundary conditions for equations on the half line to Fokas-Lenells equation on the finite interval.
基金support by NSFC(11671095,11501365)Shanghai Sailing Program supported by Science and Technology Commission of Shanghai Municipality under Grant NO.15YF1408100+1 种基金Shanghai Youth Teacher Assistance program NO.ZZslg15056the Hujiang Foundation of China(B14005)
文摘The Gelfand-Levitan-Marchenko representation is used to analyze the initialboundary value problem of two-component nonlinear Schr¨odinger equation on the half-line.It has shown that the global relation can be effectively analyzed by the Gelfand-LevitanMarchenko representation. we also derive expressions for the Dirichlet-to-Neumann map to characterize the unknown boundary values.
文摘The Method Of Lines (MOL) and Scheifele’s G-functions in the design of algorithms adapted for the numeric integration of parabolic Partial Differential Equations (PDE) in one space dimension are applied. The semi-discrete system of ordinary differential equations in the time direction, obtained by applying the MOL to PDE, is solved with the use of a method of Adapted Series, based on Scheifele’s G-functions. This method integrates exactly unperturbed linear systems of ordinary differential equations, with only one G-function. An implementation of this algorithm is used to approximate the solution of two test problems proposed by various authors. The results obtained by the Dufort-Frankel, Crank-Nicholson and the methods of Adapted Series versus the analytical solution, show the results of mistakes made.
文摘This paper is concerned with the initial-boundary value problem of a nonlinear conservation law in the half space R+= {x |x > 0} where a>0 , u(x,t) is an unknown function of x ∈ R+ and t>0 , u ± , um are three given constants satisfying um=u+≠u- or um=u-≠u+ , and the flux function f is a given continuous function with a weak discontinuous point ud. The main purpose of our present manuscript is devoted to studying the structure of the global weak entropy solution for the above initial-boundary value problem under the condition of f '-(ud) > f '+(ud). By the characteristic method and the truncation method, we construct the global weak entropy solution of this initial-boundary value problem, and investigate the interaction of elementary waves with the boundary and the boundary behavior of the weak entropy solution.
文摘By using fluid dynamics theory with the effects of adsorption and reaction, the chromatography model with a reaction A →B was established as a system of two hyperbolic partial differential equations (PDE’s). In some practical situations, the reaction chromatography model was simplified a semi-coupled system of two linear hyperbolic PDE’s. In which, the reactant concentration wave model was the initial-boundary value problem of a self-closed hyperbolic PDE, while the resultant concentration wave model was the initial-boundary value problem of hyperbolic PDE coupling reactant concentration. The general explicit expressions for the concentration wave of the reactants and resultants were derived by Laplace transform. The δ-pulse and wide pulse injections were taken as the examples to discuss detailedly, and then the stability analysis between the resultant solutions of the two modes of pulse injection was further discussed. It was significant for further analysis of chromatography, optimizing chromatographic separation, determining the physical and chemical characters.
基金supported by the State Key Program of National Natural Science Foundation of China(Grants No.11731008)the National Natural Science Foundation of China(Grants No.10771087)。
文摘The ideal reaction chromatography model can be regarded as a semi-coupled system of two hyperbolic partial differential equations, in which, one is a self-closed nonlinear equation for the reactant concentration and another is a linear equation coupling the reactant concentration for the resultant concentration. This paper is concerned with the initial-boundary value problem for the above model. By the characteristic method and the truncation method, we construct the global weak entropy solution of this initial initial-boundary value problem for Riemann type of initial-boundary data. Moreover, as examples, we apply the obtained results to the cases of head-on and wide pulse injections and give the expression of the global weak entropy solution.
文摘We study an finite-difference time-domain (FDTD) system of uniaxial perfectly matched layer (UPML) method for electromagnetic scattering problems. Particularly we analyze the discrete initial-boundary value problems of the transverse magnetic mode (TM) to Maxwell's equations with Yee's algorithm. An exterior domain in two spacial dimension is truncated by a square with a perfectly matched layer filled by a certain artificial medium. Besides, an artificial boundary condition is imposed on the outer boundary of the UPML. Using energy method, we obtain the stability of this FDTD system on the truncated domain. Numerical experiments are designed to approve the theoretical analysis.
基金Supported by National Natural Science Foundation of China under Grant Nos.11271008 and 61072147
文摘In this paper, the Fokas unified method is used to analyze the initial-boundary value problem of a complex Sharma–Tasso–Olver(c STO) equation on the half line. We show that the solution can be expressed in terms of the solution of a Riemann–Hilbert problem. The relevant jump matrices are explicitly given in terms of the matrix-value spectral functions spectral functions {a(λ), b(λ)} and {A(λ), B(λ)}, which depending on initial data u_0(x) = u(x, 0) and boundary data g_0(y) = u(0, y), g_1(y) = ux(0, y), g_2(y) = u_(xx)(0, y). These spectral functions are not independent, they satisfy a global relation.
文摘In this paper we construct the suitable weak solution for the initial-boundary valueproblem of the Boussinesq equations and obtain some properties for these solutions.Also in the caseof a two dimensional space the uniqueness of weak solution is proved.
文摘Initial-boundary value problems for a class of systems of semilinear parabolicequations with singular terms are investigated, conditions for the existence and nonexistence of the global solutions to the above problems are given, and upper bounds of thequenching time and the critical length are obtained.
基金supported by National Natural Science Foundation of China (Grant Nos. 11331005,11771150,11601164 and 11601165)
文摘This paper is concerned with the asymptotic behavior of the solution to the Euler equations with time-depending damping on quadrant(x,t)∈R^+×R^+,with the null-Dirichlet boundary condition or the null-Neumann boundary condition on u. We show that the corresponding initial-boundary value problem admits a unique global smooth solution which tends timeasymptotically to the nonlinear diffusion wave. Compared with the previous work about Euler equations with constant coefficient damping, studied by Nishihara and Yang(1999), and Jiang and Zhu(2009, Discrete Contin Dyn Syst), we obtain a general result when the initial perturbation belongs to the same space. In addition,our main novelty lies in the fact that the cut-off points of the convergence rates are different from our previous result about the Cauchy problem. Our proof is based on the classical energy method and the analyses of the nonlinear diffusion wave.