A HYBRID METHOD FOR STIFF REACTION-DIFFUSION EQUATIONS.

The second-order implicit integration factor method (IIF2) is effective at solving stiff reaction-diffusion equations owing to its nice stability condition. IIF has previously been applied primarily to systems in which the reaction contained no explicitly time-dependent terms and the boundary conditions were homogeneous. If applied to a system with explicitly time-dependent reaction terms, we find that IIF2 requires prohibitively small time-steps, that are relative to the square of spatial grid sizes, to attain its theoretical second-order temporal accuracy. Although the second-order implicit exponential time differencing (iETD2) method can accurately handle explicitly time-dependent reactions, it is more computationally expensive than IIF2. In this paper, we develop a hybrid approach that combines the advantages of both methods, applying IIF2 to reaction terms that are not explicitly time-dependent and applying iETD2 to those which are. The second-order hybrid IIF-ETD method (hIFE2) inherits the lower complexity of IIF2 and the ability to remain second-order accurate in time for large time-steps from iETD2. Also, it inherits the unconditional stability from IIF2 and iETD2 methods for dealing with the stiffness in reaction-diffusion systems. Through a transformation, hIFE2 can handle nonhomogeneous boundary conditions accurately and efficiently. In addition, this approach can be naturally combined with the compact and array representations of IIF and ETD for systems in higher spatial dimensions. Various numerical simulations containing linear and nonlinear reactions are presented to demonstrate the superior stability, accuracy, and efficiency of the new hIFE method

A survey on stationary problems, Green's functions and spectrum of Sturm–Liouville problem with nonlocal boundary conditions

In this paper, we present a survey of recent results on the Green's functions and on spectrum for stationary problems with nonlocal boundary conditions. Results of Lithuanian mathematicians in the field of differential and numerical problems with nonlocal boundary conditions are described. *The research was partially supported by the Research Council of Lithuania (grant No. MIP-047/2014)

Preventing blow up by convective terms in dissipative PDEs

We study the impact of the convective terms on the global solvability or finite time blow up of solutions of dissipative PDEs. We consider the model examples of 1D Burger's type equations, convective Cahn-Hilliard equation, generalized Kuramoto-Sivashinsky equation and KdV type equations, we establish the following common scenario: adding sufficiently strong (in comparison with the destabilizing nonlinearity) convective terms to equation prevents the solutions from blowing up in finite time and makes the considered system globally well-posed and dissipative and for weak enough convective terms the finite time blow up may occur similarly to the case when the equation does not involve convective term. This kind of result has been previously known for the case of Burger's type equations and has been strongly based on maximum principle. In contrast to this, our results are based on the weighted energy estimates which do not require the maximum principle for the considered problem

Well-posedness and exponential decay of solutions for the Blackstock-Crighton-Kuznetsov equation

The present work provides well-posedness and exponential decay results for the Blackstock-Crighton-Kuznetsov equation arising in the modeling of nonlinear acoustic wave propagation in thermally relaxing viscous fluids. First, we treat the associated linear equation by means of operator semigroups. Moreover, we derive energy estimates which we will use in a fixed-point argument in order to obtain well-posedness of the Blackstock-Crighton-Kuznetsov equation. Using a classical barrier argument we prove exponential decay of solutions.Comment: 18 page