Boundary equations in the finite transfer method for solving differential equation systems

The finite transfer method is going to be used to solve a p system of linear ordinary differential equations. The complete problem is extended by adding the p boundary equations involved. It is chosen a fourth order scheme to obtain finite transfer expressions. A recurrence strategy is used in these equations and permits one to relate different points in the domain where boundary equations are defined. Finally a 2p algebraic system of equations is noted and solved. To show the efficiency and accuracy, the method is applied to determine the structural behavior of a bending beam with different supports and to solve a differential equation of second degree with different boundary conditions.     

A two-grid finite element method for the Allen-Cahn equation with the logarithmic potential

In this paper, we present a two-grid finite element method for the Allen-Cahn equation with the logarithmic potential. This method consists of two steps. In the first step, based on a fully implicit finite element method, the Allen-Cahn equation is solved on a coarse grid with mesh size H. In the second step, a linearized system whose nonlinear term is replaced by the value of the first step is solved on a fine grid with mesh size h. We give the energy stabilities of the traditional finite element method and the two-grid finite element method. The optimal convergence order of the two-grid finite element method in H¹ norm is achieved when the mesh sizes satisfy h = O(H²). Numerical examples are given to demonstrate the validity of the proposed scheme. The results show that the two-grid method can save the CPU time while keeping the same convergence rate.     

On positive solutions for a fourth order asymptotically linear elliptic equation under Navier boundary conditions

A result on existence of positive solution for a fourth order nonlinear elliptic equation under Navier boundary conditions is established. The nonlinear term involved is asymptotically linear both at the origin and at infinity. We exploit topological degree theory and global bifurcation.     

On the existence and Morse index of solutions to the Allen-Cahn equation in two dimensions

In this paper we study transition layers in the solutions to the Allen-Cahn equation in two dimensions. We show that for any straight line segment intersecting the boundary of the domain orthogonally there exists a solution to the Allen-Cahn equation, whose transition layer is located near this segment. In addition we analyze stability of such solutions and show that it is completely determined by a geometric eigenvalue problem associated to the transition layer. We prove the existence of both stable and unstable solutions. In the case of the stable solutions we recover a result of Kohn and Sternberg [13]. As for the unstable solutions we show that their Morse index is either 1 or 2. Mathematics Subject Classification (2000) 35J60, 35Q72, 35J20, 35P15, 35P20, 35B25, 35B35, 35B40, 35B41     

Vanishing viscosity limit for a coupled Navier-Stokes/Allen-Cahn system

In this paper, we study the vanishing viscosity limit for a coupled Navier-Stokes/Allen-Cahn system in a bounded domain. We first show the local existence of smooth solutions of the Euler/Allen-Cahn equations by modified Galerkin method. Then using the boundary layer function to deal with the mismatch of the boundary conditions between Navier-Stokes and Euler equations, and assuming that the energy dissipation for Navier-Stokes equation in the boundary layer goes to zero as the viscosity tends to zero, we prove that the solutions of the Navier-Stokes/Allen-Cahn system converge to that of the Euler/Allen-Cahn system in a proper small time interval. In addition, for strong solutions of the Navier-Stokes/Allen-Cahn system in 2D, the convergence rate is cν1/2.     

Modified combined grid method for solving the Dirichlet problem for the Laplace equation on a rectangular parallelepiped

A modified combined grid method is proposed for solving the Dirichlet problem for the Laplace equation on a rectangular parallelepiped. The six-point averaging operator is applied at next-to-the-boundary grid points, while the 18-point averaging operator is used instead of the 26-point one at the remaining grid points. Assuming that the boundary values given on the faces have fourth derivatives satisfying the Hölder condition, the boundary values on the edges are continuous, and their second derivatives obey a matching condition implied by the Laplace equation, the grid solution is proved to converge uniformly with the fourth order with respect to the mesh size.     

Multiplicity results on a nonlinear biharmonic equation

Let ω be a bounded open set in Rn with smooth boundary ω We are concerned with a fourth order semilinear elliptic boundary value problem Δ²u + cΔu = bu⁺ + s inω under Dirichlet boundary condition. We investigate the existence of solutions of the fourth order nonlinear equation (0.1) when the nonlinearity bu⁺ crosses eigenvalues of Δ² + cΔ under Dirichlet boundary condition.     

One-Dimensional Symmetry and Liouville Type Results for the Fourth Order Allen-Cahn Equation in \RN

In this paper, the authors prove an analogue of Gibbons' conjecture for the extended fourth order Allen-Cahn equation in R~N, as well as Liouville type results for some solutions converging to the same value at infinity in a given direction. The authors also prove a priori bounds and further one-dimensional symmetry and rigidity results for semilinear fourth order elliptic equations with more general nonlinearities.     

An iterative method based on equation decomposition for the fourth‐order singular perturbation problem

In this article, we propose an iterative method based on the equation decomposition technique ( 1 ) for the numerical solution of a singular perturbation problem of fourth‐order elliptic equation. At each step of the given method, we only need to solve a boundary value problem of second‐order elliptic equation and a second‐order singular perturbation problem. We prove that our approximate solution converges to the exact solution when the domain is a disc. Our numerical examples show the efficiency and accuracy of our method. Our iterative method works very well for singular perturbation problems, that is, the case of 0 < ε ? 1, and the convergence rate is very fast. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A high‐order finite difference method for 1D nonhomogeneous heat equations

In this article a sixth‐order approximation method (in both temporal and spatial variables) for solving nonhomogeneous heat equations is proposed. We first develop a sixth‐order finite difference approximation scheme for a two‐point boundary value problem, and then heat equation is approximated by a system of ODEs defined on spatial grid points. The ODE system is discretized to a Sylvester matrix equation via boundary value method. The obtained algebraic system is solved by a modified Bartels‐Stewart method. The proposed approach is unconditionally stable. Numerical results are provided to illustrate the accuracy and efficiency of our approximation method along with comparisons with those generated by the standard second‐order Crank‐Nicolson scheme as well as Sun‐Zhang's recent fourth‐order method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

一维Allen-Cahn方程紧差分格式的离散最大化原则和能量稳定性研究

本文主要研究相场模拟中的Allen-Cahn模型,考虑一维Allen-Cahn方程紧差分方法的数值逼近.建立具有O(∫τ2+h4)精度的全离散紧差分格式,证明在合理的步长比和时间步长的约束下,其数值解满足离散最大化原则,在此基础上,研究了全离散格式的能量稳定性.最后给出数值算例.     

On a combined grid method for solving the Dirichlet problem for the Laplace equation in a rectangular parallelepiped

A combined grid method for solving the Dirichlet problem for the Laplace equation in a rectangular parallelepiped is proposed. At the grid points that are at the distance equal to the grid size from the boundary, the 6-point averaging operator is used. At the other grid points, the 26-point averaging operator is used. It is assumed that the boundary values have the third derivatives satisfying the Lipschitz condition on the faces; on the edges, they are continuous and their second derivatives satisfy the compatibility condition implied by the Laplace equation. The uniform convergence of the grid solution with the fourth order with respect to the grid size is proved     

Morse index of layered solutions to the heterogeneous Allen-Cahn equation

Let u? be a single layered radially symmetric unstable solution of the Allen-Cahn equation −?²Δu=u(u−a(|x|))(1−u) over the unit ball with Neumann boundary conditions. We estimate the small eigenvalues of the linearized eigenvalue problem at u? when ? is small. As a consequence, we prove that the Morse index of u? is asymptotically given by [μ^∗+o(1)]?−(N−1)/2 with μ^∗ a certain positive constant expressed in terms of parameters determined by the Allen-Cahn equation. Our estimates on the small eigenvalues have many other applications. For example, they may be used in the search of other non-radially symmetric solutions, which will be considered in forthcoming papers.     

On an initial‐boundary value problem for a wide‐angle parabolic equation in a waveguide with a variable bottom

We consider the third‐order Claerbout‐type wide‐angle parabolic equation (PE) of underwater acoustics in a cylindrically symmetric medium consisting of water over a soft bottom B of range‐dependent topography. There is strong indication that the initial‐boundary value problem for this equation with just a homogeneous Dirichlet boundary condition posed on B may not be well‐posed, for example when B is downsloping. We impose, in addition to the above, another homogeneous, second‐order boundary condition, derived by assuming that the standard (narrow‐angle) PE holds on B, and establish a priori H² estimates for the solution of the resulting initial‐boundary value problem for any bottom topography. After a change of the depth variable that makes B horizontal, we discretize the transformed problem by a second‐order accurate finite difference scheme and show, in the case of upsloping and downsloping wedge‐type domains, that the new model gives stable and accurate results. We also present an alternative set of boundary conditions that make the problem exactly energy conserving; one of these conditions may be viewed as a generalization of the Abrahamsson–Kreiss boundary condition in the wide‐angle case. Copyright © 2008 John Wiley & Sons, Ltd.     

Inverse viscosity problem for the Navier-Stokes equation

We consider an inverse problem of determining a viscosity coefficient in the Navier-Stokes equation by observation data in a neighborhood of the boundary. We prove the Lipschitz stability by the Carleman estimates in Sobolev spaces of negative order.     

Strong solutions to the Richards equation in the unsaturated zone

We study the unsaturated case of the Richards equation in three space dimensions with Dirichlet boundary data. We first establish an a priori L^∞-estimate. With its help, by means of a fixed point argument we prove global in time existence of a unique weak solution in Sobolev spaces. Finally, we are able to improve the regularity of this weak solution in order to gain a strong one.     

一类非线性偏微分方程的四阶格子Boltzmann模型

通过Chapman-Enskog展开技术和多尺度分析,建立了一种新的D1Q4带修正项的四阶格子Boltzmann模型,一类非线性偏微分方程从连续的Boltzmann方程得到正确恢复．统一了KdV和Burgers等已知方程类型的格子BGK模型,还首次给出了组合KdV-Burgers,广义Burgers—Huxley等方程...     

Fast direct solver for Poisson equation in a 2D elliptical domain

In this article, we extend our previous work 3 for developing some fast Poisson solvers on 2D polar and spherical geometries to an elliptical domain. Instead of solving the equation in an irregular Cartesian geometry, we formulate the equation in elliptical coordinates. The solver relies on representing the solution as a truncated Fourier series, then solving the differential equations of Fourier coefficients by finite difference discretizations. Using a grid by shifting half mesh away from the pole and incorporating the derived numerical boundary value, the difficulty of coordinate singularity can be elevated easily. Unlike the case of 2D disk domain, the present difference equation for each Fourier mode is coupled with its conjugate mode through the numerical boundary value near the pole; thus, those two modes are solved simultaneously. Both second‐ and fourth‐order accurate schemes for Dirichlet and Neumann problems are presented. In particular, the fourth‐order accuracy can be achieved by a three‐point compact stencil which is in contrast to a five‐point long stencil for the disk case. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 72–81, 2004     

Conditions at infinity for boltzmann's equation.

We consider here realistic conditions at infinity for solutions of the Boltzmann's equation, such as a pure Maxwellian equilibrium at infinity possibly with suitable boundary conditions on an exterior domain, different Maxwellian equilibria at +∞ and -;∞ in a tube-like situation and more generally conditions at infinity obtained from a fixed solution. In order to adapt the recent global existence and compactness results due to R.J. DiPerna and the author, we have to obtain some local a priori estimates on the mass, kinetic energy and entropy. And this is precisely what we achieve here by two different and new methods. The first one consists in using the relative entropy of solutions with respect to a fixed, possibly local, Maxwellian. This method allows to treat general collision kernels with angular cut-off and some of the conditions at infinity mentioned above. The second method is based upon a L¹ estimate and an extension of the entropy identity which uses a truncated H-functional. This method requires a “uniform integrability” condition on the collision kernel but allows to consider the most general conditions at infinity.