Anomalous scaling for three‐dimensional Cahn‐Hilliard fronts

We prove the stability of the one‐dimensional kink solution of the Cahn‐Hilliard equation under d‐dimensional perturbations for d ≥ 3. We also establish a novel scaling behavior of the large‐time asymptotics of the solution. The leading asymptotics of the solution is characterized by a length scale proportional to t^1/3 instead of the usual t^1/2 scaling typical to parabolic problems. © 2004 Wiley Periodicals, Inc.     

Well‐posedness for the incompressible magneto‐hydrodynamic system

This paper is concerned with well‐posedness of the incompressible magneto‐hydrodynamics (MHD) system. In particular, we prove the existence of a global mild solution in BMO^?1 for small data which is also unique in the space C([0, ∞); BMO^?1). We also establish the existence of a local mild solution in bmo^?1 for small data and its uniqueness in C([0, T); bmo^?1). In establishing our results an important role is played by the continuity of the bilinear form which was proved previously by Kock and Tataru. In this paper, we give a new proof of this result by using the weighted L^p‐boundedness of the maximal function. Copyright © 2006 John Wiley & Sons, Ltd.     

The Vlasov‐Poisson‐Boltzmann system near Maxwellians

The dynamics of dilute electrons can be modeled by the Vlasov‐Poisson‐Boltz‐mann system, where electrons interact with themselves through collisions and with their self‐consistent electric field. It is shown that any smooth, periodic initial perturbation of a given global Maxwellian that preserves the same mass, momentum, and total energy (including both kinetic and electric energy), leads to a unique global‐in‐time classical solution. The construction of global solutions is based on an energy method with a new estimate of dissipation from the collision: ∫₀^t〈Lf(s), f(s)〉ds is positive definite for solution f(t,x,v) with small amplitude to the Vlasov‐Poisson‐Boltzmann system (1.4). © 2002 Wiley Periodicals, Inc.     

Superconvergence of the direct discontinuous Galerkin method for a time‐fractional initial‐boundary value problem

A time‐fractional reaction–diffusion initial‐boundary value problem with periodic boundary condition is considered on Q ? Ω × [0, T] , where Ω is the interval [0, l] . Typical solutions of such problem have a weak singularity at the initial time t = 0. The numerical method of the paper uses a direct discontinuous Galerkin (DDG) finite element method in space on a uniform mesh, with piecewise polynomials of degree k ≥ 2 . In the temporal direction we use the L1 approximation of the Caputo derivative on a suitably graded mesh. We prove that at each time level of the mesh, our L1‐DDG solution is superconvergent of order k + 2 in L²(Ω) to a particular projection of the exact solution. Moreover, the L1‐DDG solution achieves superconvergence of order (k + 2) in a discrete L²(Q) norm computed at the Lobatto points, and order (k + 1) superconvergence in a discrete H¹(Q) seminorm at the Gauss points; numerical results show that these estimates are sharp.     

Asymptotic analysis and estimates of blow‐up time for the radial symmetric semilinear heat equation in the open‐spectrum case

We estimate the blow‐up time for the reaction diffusion equation u_t=Δu+ λf(u), for the radial symmetric case, where f is a positive, increasing and convex function growing fast enough at infinity. Here λ>λ^*, where λ^* is the ‘extremal’ (critical) value for λ, such that there exists an ‘extremal’ weak but not a classical steady‐state solution at λ=λ^* with ∥w(?, λ)∥_∞→∞ as 0<λ→λ^*?. Estimates of the blow‐up time are obtained by using comparison methods. Also an asymptotic analysis is applied when f(s)=e^s, for λ?λ^*?1, regarding the form of the solution during blow‐up and an asymptotic estimate of blow‐up time is obtained. Finally, some numerical results are also presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Computable analysis of a boundary‐value problem for the generalized KdV–Burgers equation

In this paper, we investigate the computability of the solution operator of the generalized KdV‐Burgers equation with initial‐boundary value problem. Here, the solution operator is a nonlinear map H^{3m ? 1}(R⁺) × H^m(0,T)→C([0,T];H^{3m ? 1}(R⁺)) from the initial‐boundary value data to the solution of the equation. By a technique that is widely used for the study of nonlinear dispersive equation, and using the type 2 theory of effectivity as computable model, we prove that the solution map is Turing computable, for any integer m ≥ 2, and computable real number T > 0. Copyright © 2014 John Wiley & Sons, Ltd.     

On the decay and stability of global solutions to the 3‐D inhomogeneous Navier‐Stokes equations

In this paper, we investigate the large‐time decay and stability to any given global smooth solutions of the 3‐D incompressible inhomogeneous Navier‐Stokes equations. In particular, we prove that given any global smooth solution (a,u) of (1.2), the velocity field u decays to 0 with an explicit rate, which coincides with the L² norm decay for the weak solutions of the 3‐D classical Navier‐Stokes system [26,29] as t goes to ∞. Moreover, a small perturbation to the initial data of (a,u) still generates a unique global smooth solution to (1.2), and this solution keeps close to the reference solution (a,u) for t > 0. We should point out that the main results in this paper work for large solutions of (1.2). © 2010 Wiley Periodicals, Inc.     

Combination of standard Galerkin and subspace methods for the time‐dependent Navier‐Stokes equations with nonsmooth initial data

We consider a combination of the standard Galerkin method and the subspace decomposition methods for the numerical solution of the two‐dimensional time‐dependent incompressible Navier‐Stokes equations with nonsmooth initial data. Because of the poor smoothness of the solution near t = 0, we use the standard Galerkin method for time interval [0, 1] and the subspace decomposition method time interval [1, ∞). The subspace decomposition method is based on the solution into the sum of a low frequency component integrated using a small time step Δt and a high frequency integrated using a larger time step pΔt with p > 1. From the H¹‐stability and L²‐error analysis, we show that the subspace decomposition method can yield a significant gain in computing time. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2009     

Global existence to a reaction–diffusion equation through a self‐similar‐like upper solution

A demonstration method is presented, which will ensure the existence of positive global solutions in time to the reaction–diffusion equation ?u_t+Δu+u^p=0 in ?ⁿ×[0, ∞), for exponents p?3 and space dimensions n?3. This method does not require the initial value to have a specific uniform smallness condition, but rather to satisfy a bell‐like form. The method is based on a specific upper solution, which models the diffusion process of the heat equation. The upper solution is not self‐similar, but does have a self‐similar‐like form. After transforming the reaction–diffusion problem into an equivalent one, whose initial value is uniformly very small, a local solution is obtained in the time interval [0, 1] by the use of this upper solution. This local solution is then extended to [0, ∞) through an infinite sequence of extensions. At each step, an appropriate change of variables will transform the extension into a problem nearly identical to the local problem in [0, 1]. These transformations exploit the diffusive and self‐similar‐like nature of the upper solution. Copyright © 2009 John Wiley & Sons, Ltd.     

Large time behavior for the incompressible magnetohydrodynamic equations in half‐spaces

The weighted L^r‐asymptotic behavior of the strong solution and its first‐order spacial derivatives to the incompressible magnetohydrodynamic (MHD) equations is established in a half‐space. Further, the L^∞‐decay rates of the second‐order spatial derivatives of the strong solution are derived by using the Stokes solution formula and employing a decomposition for the nonlinear terms in MHD equations. Copyright © 2014 John Wiley & Sons, Ltd.     

A two‐grid characteristic finite volume element method for semilinear advection‐dominated diffusion equations

A two‐grid finite volume element method, combined with the modified method of characteristics, is presented and analyzed for semilinear time‐dependent advection‐dominated diffusion equations in two space dimensions. The solution of a nonlinear system on the fine‐grid space (with grid size h) is reduced to the solution of two small (one linear and one nonlinear) systems on the coarse‐grid space (with grid size H) and a linear system on the fine‐grid space. An optimal error estimate in H¹ ‐norm is obtained for the two‐grid method. It shows that the two‐grid method achieves asymptotically optimal approximation, as long as the mesh sizes satisfy h = O(H²). Numerical example is presented to validate the usefulness and efficiency of the method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Convergence analysis for second‐order accurate schemes for the periodic nonlocal Allen‐Cahn and Cahn‐Hilliard equations

In this paper, we provide a detailed convergence analysis for fully discrete second‐order (in both time and space) numerical schemes for nonlocal Allen‐Cahn and nonlocal Cahn‐Hilliard equations. The unconditional unique solvability and energy stability ensures ? ⁴ stability. The convergence analysis for the nonlocal Allen‐Cahn equation follows the standard procedure of consistency and stability estimate for the numerical error function. For the nonlocal Cahn‐Hilliard equation, because of the complicated form of the nonlinear term, a careful expansion of its discrete gradient is undertaken, and an H ^?1 inner‐product estimate of this nonlinear numerical error is derived to establish convergence. In addition, an a priori bound of the numerical solution at the discrete level is needed in the error estimate. Such a bound can be obtained by performing a higher order consistency analysis by using asymptotic expansions for the numerical solution. Following the technique originally proposed by Strang (eg, 1964), instead of the standard comparison between the exact and numerical solutions, an error estimate between the numerical solution and the constructed approximate solution yields an O (s ³+h ⁴) convergence in norm, in which s and h denote the time step and spatial mesh sizes, respectively. This in turn leads to the necessary bound under a standard constraint s ≤C h . Here, we also prove convergence of the scheme in the maximum norm under the same constraint.     

Asymptotic growth bounds for the 3‐D Vlasov‐Poisson system with point charges

In this article, we consider the asymptotic behavior of the classical solution to the 3‐dimensional Vlasov‐Poisson plasma interacting repulsively with N point charges. The large time behavior in terms of diameters of its velocity‐spatial supports is improved to O(t^2/3+ϵ) for any ϵ>0.     

Stepanov‐like doubly weighted pseudo almost automorphic processes and its application to Sobolev‐type stochastic differential equations driven by G‐Brownian motion

This paper investigates the properties of the p‐mean Stepanov‐like doubly weighted pseudo almost automorphic (S^pDWPAA) processes and its application to Sobolev‐type stochastic differential equations driven by G‐Brownian motion. We firstly prove the equivalent relation between the S^pDWPAA and Stepanov‐like asymptotically almost automorphic stochastic processes based on ergodic zero set. We further establish the completeness of the space and the composition theorem for S^pDWPAA processes. These results obtained improve and extend previous related conclusions. As an application, we show the existence and uniqueness of the S^p DWPAA solution for a class of nonlinear Sobolev‐type stochastic differential equations driven by G‐Brownian motion and present a decomposition of this unique solution. Moreover, an example is given to illustrate the effectiveness of our results.     

A regularity criterion for two‐and‐half‐dimensional magnetohydrodynamic equations with horizontal dissipation and horizontal magnetic diffusion

We study the global regularity of classical solution to two‐and‐half‐dimensional magnetohydrodynamic equations with horizontal dissipation and horizontal magnetic diffusion. We prove that any possible finite time blow‐up can be controlled by the L^∞‐norm of the vertical components. Copyright © 2016 John Wiley & Sons, Ltd.     

A high‐order B‐spline collocation scheme for solving a nonhomogeneous time‐fractional diffusion equation

A high‐accuracy numerical approach for a nonhomogeneous time‐fractional diffusion equation with Neumann and Dirichlet boundary conditions is described in this paper. The time‐fractional derivative is described in the sense of Riemann‐Liouville and discretized by the backward Euler scheme. A fourth‐order optimal cubic B‐spline collocation (OCBSC) method is used to discretize the space variable. The stability analysis with respect to time discretization is carried out, and it is shown that the method is unconditionally stable. Convergence analysis of the method is performed. Two numerical examples are considered to demonstrate the performance of the method and validate the theoretical results. It is shown that the proposed method is of order O(Δx⁴ + Δt^{2 ? α}) convergence, where α ∈ (0,1) . Moreover, the impact of fractional‐order derivative on the solution profile is investigated. Numerical results obtained by the present method are compared with those obtained by the method based on standard cubic B‐spline collocation method. The CPU time for present numerical method and the method based on cubic B‐spline collocation method are provided.     

Global Well‐Posedness of the Three‐Dimensional Primitive Equations with Only Horizontal Viscosity and Diffusion

In this paper, we consider the initial boundary value problem of the three‐dimensional primitive equations for planetary oceanic and atmospheric dynamics with only horizontal eddy viscosity in the horizontal momentum equations and only horizontal diffusion in the temperature equation. Global well‐posedness of the strong solution is established for any H² initial data. An N‐dimensional logarithmic Sobolev embedding inequality, which bounds the L^∞‐norm in terms of the L^q‐norms up to a logarithm of the L^p‐norm for p > N of the first‐order derivatives, and a system version of the classic Grönwall inequality are exploited to establish the required a~priori H² estimates for global regularity.© 2016 Wiley Periodicals, Inc.     

Least‐squares solutions and least‐rank solutions of the matrix equation AXA* = B and their relations

A Hermitian matrix X is called a least‐squares solution of the inconsistent matrix equation AXA^* = B, where B is Hermitian. A^* denotes the conjugate transpose of A if it minimizes the F‐norm of B ? AXA^*; it is called a least‐rank solution of AXA^* = B if it minimizes the rank of B ? AXA^*. In this paper, we study these two types of solutions by using generalized inverses of matrices and some matrix decompositions. In particular, we derive necessary and sufficient conditions for the two types of solutions to coincide. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite difference discretization of the Benjamin‐Bona‐Mahony‐Burgers equation

Numerical solutions of the Benjamin‐Bona‐Mahony‐Burgers equation in one space dimension are considered using Crank‐Nicolson‐type finite difference method. Existence of solutions is shown by using the Brower's fixed point theorem. The stability and uniqueness of the corresponding methods are proved by the means of the discrete energy method. The convergence in L^∞‐norm of the difference solution is obtained. A conservative difference scheme is presented for the Benjamin‐Bona‐Mahony equation. Some numerical experiments have been conducted in order to validate the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

The generalized centro‐symmetric and least squares generalized centro‐symmetric solutions of the matrix equation AYB + CYTD = E

An n×n real matrix P is said to be a symmetric orthogonal matrix if P = P^?1 = P^T. An n × n real matrix Y is called a generalized centro‐symmetric with respect to P, if Y = PYP. It is obvious that every matrix is also a generalized centro‐symmetric matrix with respect to I. In this work by extending the conjugate gradient approach, two iterative methods are proposed for solving the linear matrix equation and the minimum Frobenius norm residual problem over the generalized centro‐symmetric Y, respectively. By the first (second) algorithm for any initial generalized centro‐symmetric matrix, a generalized centro‐symmetric solution (least squares generalized centro‐symmetric solution) can be obtained within a finite number of iterations in the absence of round‐off errors, and the least Frobenius norm generalized centro‐symmetric solution (the minimal Frobenius norm least squares generalized centro‐symmetric solution) can be derived by choosing a special kind of initial generalized centro‐symmetric matrices. We also obtain the optimal approximation generalized centro‐symmetric solution to a given generalized centro‐symmetric matrix Y₀ in the solution set of the matrix equation (minimum Frobenius norm residual problem). Finally, some numerical examples are presented to support the theoretical results of this paper. Copyright © 2011 John Wiley & Sons, Ltd.