Global solutions for quasilinear parabolic problems

Results on the global existence of classical solutions for quasilinear parabolic equations in bounded domains with homogeneous Dirichlet or Neumann boundary conditions are presented. Besides quasilinear parabolic equations, the method is also applicable to some weakly-coupled reaction-diffusion systems and to elliptic equations with nonlinear dynamic boundary conditions. Received December 21, 2000; accepted August 30, 2001.     

Degenerate Second Order Differential Operators Generating Analytic Semigroups in Lp and W1,p

We deal with the problem of analyticity for the semigroup generated by the second order differential operator Au ≔ αu″ + βu′ (or by some restrictions of it) in the spaces L^p(0, 1), with or without weight, and in W^1,p(0, 1), 1 < p < ∞. Here α and β are assumed real‐valued and continuous in [0, 1], with α(x) > 0 in (0, 1), and the domain of A is determined by the generalized Neumann boundary conditions and by Wentzell boundary conditions.     

Finite energy solutions of self‐adjoint elliptic mixed boundary value problems

This paper describes existence, uniqueness and special eigenfunction representations of H¹‐solutions of second order, self‐adjoint, elliptic equations with both interior and boundary source terms. The equations are posed on bounded regions with Dirichlet conditions on part of the boundary and Neumann conditions on the complement. The system is decomposed into separate problems defined on orthogonal subspaces of H¹(Ω). One problem involves the equation with the interior source term and the Neumann data. The other problem just involves the homogeneous equation with Dirichlet data. Spectral representations of the solution operators for each of these problems are found. The solutions are described using bases that are, respectively, eigenfunctions of the differential operator with mixed null boundary conditions, and certain mixed Steklov eigenfunctions. These series converge strongly in H¹(Ω). Necessary and sufficient conditions for the Dirichlet part of the boundary data to have a finite energy extension are described. The solutions for a problem that models a cylindrical capacitor is found explicitly. Copyright © 2009 John Wiley & Sons, Ltd.     

Conditional symmetries and exact solutions of nonlinear reaction-diffusion systems with non-constant diffusivities

Q-conditional symmetries (nonclassical symmetries) for the general class of two-component reaction-diffusion systems with non-constant diffusivities are studied. Using the recently introduced notion of Q-conditional symmetries of the first type, an exhausted list of reaction-diffusion systems admitting such symmetry is derived. The results obtained for the reaction-diffusion systems are compared with those for the scalar reaction-diffusion equations. The symmetries found for reducing reaction-diffusion systems to two-dimensional dynamical systems, i.e., ODE systems, and finding exact solutions are applied. As result, multiparameter families of exact solutions in the explicit form for a nonlinear reaction-diffusion system with an arbitrary diffusivity are constructed. Finally, the application of the exact solutions for solving a biologically and physically motivated system is presented.     

The lower bounds of life span of classical solutions to one‐dimensional initial‐Neumann boundary value problems for general quasilinear wave equations

In this paper, we will study the lower bounds of the life span (the maximal existence time) of solutions to the initial‐boundary value problems with small initial data and zero Neumann boundary data on exterior domain for one‐dimensional general quasilinear wave equations u_tt?u_xx=b(u,Du)u_xx+F(u,Du). Our lower bounds of the life span of solutions in the general case and special case are shorter than that of the initial‐Dirichlet boundary value problem for one‐dimensional general quasilinear wave equations. We clarify that although the lower bounds in this paper are same as that in the case of Robin boundary conditions obtained in the earlier paper, however, the results in this paper are not the trivial generalization of that in the case of Robin boundary conditions because the fundamental Lemmas 2.4, 2.5, 2.6, and 2.7, that is, the priori estimates of solutions to initial‐boundary value problems with Neumann boundary conditions, are established differently, and then the specific estimates in this paper are different from that in the case of Robin boundary conditions. Another motivation for the author to write this paper is to show that the well‐posedness of problem 1.1 is the essential precondition of studying the lower bounds of life span of classical solutions to initial‐boundary value problems for general quasilinear wave equations. The lower bound estimates of life span of classical solutions to initial‐boundary value problems is consistent with the actual physical meaning. Finally, we obtain the sharpness on the lower bound of the life span 1.8 in the general case and 1.10 in the special case. Copyright © 2015 John Wiley & Sons, Ltd.     

Fast direct solvers for Poisson equation on 2D polar and spherical geometries

A simple and efficient class of FFT‐based fast direct solvers for Poisson equation on 2D polar and spherical geometries is presented. These solvers rely on the truncated Fourier series expansion, where the differential equations of the Fourier coefficients are solved by the second‐ and fourth‐order finite difference discretizations. Using a grid by shifting half mesh away from the origin/poles, and incorporating with the symmetry constraint of Fourier coefficients, the coordinate singularities can be easily handled without pole condition. By manipulating the radial mesh width, three different boundary conditions for polar geometry including Dirichlet, Neumann, and Robin conditions can be treated equally well. The new method only needs O(MN log₂ N) arithmetic operations for M × N grid points. © 2002 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 18: 56–68, 2002     

Finite difference schemes for the fourth-order parabolic equations with different boundary value conditions

In this paper, the fourth-order parabolic equations with different boundary value conditions are studied. Six kinds of boundary value conditions are proposed. Several numerical differential formulae for the fourth-order derivative are established by the quartic interpolation polynomials and their truncation errors are given with the aid of the Taylor expansion with the integral remainders. Effective difference schemes are presented for the third Dirichlet boundary value problem, the first Neumann boundary value problem and the third Neumann boundary value problem, respectively. Some new embedding inequalities on the discrete function spaces are presented and proved. With the method of energy analysis, the unique solvability, unconditional stability and unconditional convergence of the difference schemes are proved. The convergence orders of derived difference schemes are all O(τ² + h²) in appropriate norms. Finally, some numerical examples are provided to confirm the theoretical results.     

Bifurcation analysis of amathematical model for the growth of solid tumors in the presence of external inhibitors

We study bifurcations from radial solution of a free boundary problem modeling the dormant state of nonnecrotic solid tumors in the presence of external inhibitors. This problem consists in three linear elliptic equations with two Dirichlet and one Neumann boundary conditions and a fourth boundary condition coupling surface tension effects on free boundary. In this paper, surface tension coefficient γ plays the role of bifurcation parameter. We prove that in certain situations there exists a positive null point sequence for γ where bifurcation occurs from radial solution, while in the other situations, either bifurcation occurs at only finite many points of γ or even it does not occur for any γ > 0. Copyright © 2014 John Wiley & Sons, Ltd.     

Cosine families generated by second-order differential operators on W1,1(0, 1) with generalized Wentzell boundary conditions

Using the abstract framework [Bátkai, A. and Engel, K.-J., 2004, Abstract wave equations with generalized Wentzell boundary conditions. Journal of Differential Equations, 207, 1–20.] we show that certain second-order differential operators with generalized Wentzell boundary conditions generate cosine families and hence also analytic semigroups on W^1,1(0,1). This complements the main result [Favini, A., Ruiz Goldstein, G., Goldstein, J.A., Obrecht, E. and Romanelli, S., 2003, General Wentzell boundary conditions and analytic semigroups on W^{1, p} (0,1). Applicable Analysis, 82, 927–935.] on the generation of an analytic semigroup by the second derivative with generalized Wentzell boundary conditions on W^{1, p} (0,?1) for 1<p<∞.     

Regularity for Fully Nonlinear Elliptic Equations with Neumann Boundary Data

We obtain local C ^α, C ^1,α, and C ^2,α regularity results up to the boundary for viscosity solutions of fully nonlinear uniformly elliptic second order equations with Neumann boundary conditions.     

Scaling Methods and Approximative Equations for Homogeneous Reaction—Diffusion Systems and Applications to Epidemics

Summary. We consider a reaction-diffusion equation that is homogeneous of degree one. This homogeneity is a symmetry. The dynamics is factorized into trivial evolution due to symmetry and nontrivial behavior by a projection to an appropriate hypermanifold. The resulting evolution equations are rather complex. We examine the bifurcation behavior of a stationary point of the projected system. For these purposes we develop techniques for dimension reduction similar to the Ginzburg-Landau (GL) approximation, the modulation equations. Since we are not in the classical GL situation, the remaining approximative equations have a quadratic nonlinearity and the amplitude does not scale with ε but with ε ² where ε ² denotes the bifurcation parameter. Moreover, the symmetry requires that not only one but two equations are necessary to describe the behavior of the system. We investigate traveling fronts for the modulation equations. This result is used to analyze an epidemic model. Received April 9, 1996; second revision received January 3, 1997; final revision received October 7, 1997; accepted January 19, 1998     

Spatial inhomogenity due to turing bifurcation in a system of Gierer-Meinhardt type

This paper treats the conditions for the existence and stability properties of stationary solutions of reaction-diffusion equations of Gierer-Meinhardt type, subject to Neumann boundary data. The domains in which diffusion takes place are of three types: a regular hexagon, a rectangle and an isosceles rectangular triangle. Considering one of the relevant features of the domains as a bifurcation parameter it will be shown that at a certain critical value a diffusion driven instability occurs and Turing bifurcation takes place: a pattern emerges.     

Fast Fourier Transform Solvers and Preconditioners for Quadratic Spline Collocation

Quadratic Spline Collocation (QSC) methods of optimal order of convergence have been recently developed for the solution of elliptic Partial Differential Equations (PDEs). In this paper, linear solvers based on Fast Fourier Transforms (FFT)are developed for the solution of the QSC equations. The complexity of the FFT solvers is O(N ² log N), where N is the gridsize in one dimension. These direct solvers can handle PDEs with coefficients in one variable or constant, and Dirichlet, Neumann, alternating Dirichlet-Neumann or periodic boundary conditions, along at least one direction of a rectangular domain. General variable coefficient PDEs are handled by preconditioned iterative solvers. The preconditioner is the QSC matrix arising from a constant coefficient PDE. The convergence analysis of the preconditioner is presented. It is shown that, under certain conditions, the convergence rate is independent of the gridsize. The preconditioner is solved by FFT techniques, and integrated with one-step or acceleration methods, giving rise to asymptotically almost optimal linear solvers, with complexity O(N ² log N). Numerical experiments verify the effectiveness of the solvers and preconditioners, even on problems more general than the analysis assumes. The development and analysis of FFT solvers and preconditioners is extended to QSC equations corresponding to systems of elliptic PDEs.     

Complete symmetry group and nonlocal symmetries for some two-dimensional evolution equations

The complete symmetry group of an 1+1 evolution equation of maximal symmetry has been demonstrated to be represented by the six-dimensional Lie algebra of point symmetries sl(2,R)s_⊕W, where W is the three-dimensional Heisenberg-Weyl algebra. We construct a complete symmetry group of a 1+2 evolution equation ut=(Fy(u)ux) for some functions F using the point symmetries admitted by the equation. The 1+2 equation is not completely specifiable by point symmetries alone for some specific functions F. We make use of Ansätze already reported by Myeni and Leach [S.M. Myeni, P.G.L. Leach, Nonlocal symmetries and complete symmetry groups of evolution equations, J. Nonlinear Math. Phys. 13 (2006) 377-392] which provide a route to the determination of the required generic nonlocal symmetries necessary to supplement the point symmetries for the complete specification of these 1+2 evolution equations. Further we find that taking some suitable linear combination of Lie point symmetries helps to optimise the procedure of specifying the equation. A general result concerning the number of symmetries required to form a complete symmetry group of evolution is presented in the Conclusion.     

Planforms in three dimensions

The spatially periodic, steady-state solutions to systems of partial differential equations (PDE) are calledplanforms. There already exists a partial classification of the planforms for Euclidean equivariant systems of PDE inR ² (see [6, 7]), In this article we attempt to give such a classification for Euclidean equivariant systems of PDE inR ³. Based on the symmetry and spatial periodicity of each planform, 59 different planforms are found.We attempt to find the planforms on all lattices inR ³ that are forced to exist near a steady-state bifurcation from a trivial solution. The proof of our classification uses Liapunov-Schmidt reduction with symmetry (which can be used if we assume spatial periodicity of the solutions) and the Equivariant Branching Lemma. The analytical problem of finding planforms for systems of PDE is reduced to the algebraic problem of computing isotropy subgroups with one dimensional fixed point subspaces.The Navier-Stokes equations and reaction-diffusion equations (with constant diffusion coefficients) are examples of systems of PDE that satisfy the conditions of our classifications. In this article, we show that our classification applies to the Kuramoto-Sivashinsky equation.     

Finite difference schemes for the singularly perturbed reaction-diffusion equation in the case of spherical symmetry

The boundary value problem for the singularly perturbed reaction-diffusion parabolic equation in a ball in the case of spherical symmetry is considered. The derivatives with respect to the radial variable appearing in the equation are written in divergent form. The third kind boundary condition, which admits the Dirichlet and Neumann conditions, is specified on the boundary of the domain. The Laplace operator in the differential equation involves a perturbation parameter ?², where ? takes arbitrary values in the half-open interval (0, 1]. When ? → 0, the solution of such a problem has a parabolic boundary layer in a neighborhood of the boundary. Using the integro-interpolational method and the condensing grid technique, conservative finite difference schemes on flux grids are constructed that converge ?-uniformly at a rate of O(N ^?2ln² N + N ₀ ^?1 ), where N + 1 and N ₀ + 1 are the numbers of the mesh points in the radial and time variables, respectively.     

The Kramers–;Wannier Symmetry and S-Duality in the Two-Dimensional g03A6;4 Theory

We show that the exact beta function of the two-dimensional g⁴ theory possesses two dual symmetries. These are the Kramers–Wannier symmetry d(g) and the strong–weak-coupling symmetry, or the S-duality f(g), connecting the strong- and weak-coupling domains lying above and below the fixed point g ^*. We obtain explicit representations for the functions d(g) and f(g). The S-duality transformation f(g) allows using the high-temperature expansions to approximate the contributions of the higher-order Feynman diagrams. From the mathematical standpoint, the proposed scheme is highly unstable. Nevertheless, the approximate values of the renormalized coupling constant g ^* obtained from the duality equations agree well with the available numerical results.     

The blow-up rate for a system of heat equations with neumann boundary conditions

This paper deals with the blow-up properties of solutions to a system of heat equations u _t=Δu, v _t=Δv in B _R×(0, T) with the Neumann boundary conditions εu/εη=e ^v, εv/εη=e ^u on S _R×[0, T). The exact blow-up rates are established. It is also proved that the blow-up will occur only on the boundary. This work is supported by the National Natural Science Foundation of China     

Gidas-Ni-Nirenberg results for finite difference equations: Estimates of approximate symmetry

Are positive solutions of finite difference boundary value problems Δhu=f(u) in Ωh, u=0 on ∂Ωh as symmetric as the domain? To answer this question we first show by examples that almost arbitrary non-symmetric solutions can be constructed. This is in striking difference to the continuous case, where by the famous Gidas-Ni-Nirenberg theorem [B. Gidas, Wei-Ming Ni, L. Nirenberg, Symmetry and related problems via the maximum principle, Comm. Math. Phys. 68 (1979) 209-243] positive solutions inherit the symmetry of the underlying domain. Then we prove approximate symmetry theorems for solutions on equidistantly meshed n-dimensional cubes: explicit estimates depending on the data are given which show that the solutions become more symmetric as the discretization gets finer. The quality of the estimates depends on whether or not f(0)<0. The one-dimensional case stands out in two ways: the proofs are elementary and the estimates for the defect of symmetry are O(h) compared to O(1/|log(h)|) in the higher-dimensional case.     

Reduktionsverfahren für Differenzengleichungen bei Randwertaufgaben II

Summary In a recent paper [11], two of the authors investigated a fast reduction method for solving difference equations which approximate certain boundary value problems for Poisson's equation. In this second paper, we prove the numerical stability of the reduction method, and also report on further developments of the method. For the general case, the provided bounds for the numerical errors behave roughly like the condition numberO(n ²) of the linear system; for more realistic model problems estimates of order less thanO(n) are obtained (n ^–1=h=mesh width). The number of operations required for the reduction method isO(n ² ), for the usual five-point difference formula, as well as for the common nine-point formula with discretization error of orderh ⁴.