Viscous two‐fluid flows in perturbed unbounded domains

Two stationary plane free boundary value problems for the Navier‐Stokes equations are studied. The first problem models the viscous two‐fluid flow down a perturbed or slightly distorted inclined plane. The second one describes the viscous two‐fluid flow in a perturbed or slightly distorted channel. For sufficiently small data and under certain conditions on parameters the solvability and uniqueness results are proved for both problems. The asymptotic behaviour of the solutions is investigated. For the second problem an example of nonuniqueness is constructed. Computational results of flow problems that are very close to the above problems are presented. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Decoupled two‐grid finite element method for the time‐dependent natural convection problem I: Spatial discretization

In this article, a decoupled two grid finite element method (FEM) is proposed and analyzed for the nonsteady natural convection problem using the coarse grid numerical solutions to decouple the nonlinear coupled terms, and the corresponding optimal error estimates are derived. Compared with the standard Galerkin FEM and the usual two‐grid FEM, our algorithm not only keeps good accuracy but also saves a lot of computational cost. Some numerical examples are provided to verify the performances of the decoupled two‐grid FEM. Both theoretical analysis and numerical experiments show the efficiency and effectiveness of the decoupled two‐grid FEM for the nonsteady natural convection problem. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2135–2168, 2015     

WEB‐spline based mesh‐free finite element analysis for the heat equation and the time‐dependent Navier–Stokes equation: A survey

In this article, we give some numerical techniques and error estimates using web‐spline based mesh‐free finite element method for the heat equation and the time‐dependent Navier–Stokes equations on bounded domains. The web‐spline method uses weighted extended B‐splines on a regular grid as basis functions and does not require any grid generation. We demonstrate the method by providing numerical results for the Poisson's and stationary Stokes equation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A discrete duality finite volume discretization of the vorticity‐velocity‐pressure stokes problem on almost arbitrary two‐dimensional grids

We present an application of the discrete duality finite volume method to the numerical approximation of the vorticity‐velocity‐pressure formulation of the two‐dimensional Stokes equations, associated to various nonstandard boundary conditions. The finite volume method is based on the use of discrete differential operators obeying some discrete duality principles. The scheme may be seen as an extension of the classical Marker and Cell scheme to almost arbitrary meshes, thanks to an appropriate choice of degrees of freedom. The efficiency of the scheme is illustrated by numerical examples over unstructured triangular and locally refined nonconforming meshes, which confirm the theoretical convergence analysis led in the article. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1–30, 2015     

Error estimates for two‐level penalty finite volume method for the stationary Navier–Stokes equations

Two‐level penalty finite volume method for the stationary Navier–Stokes equations based on the P₁ ? P₀ element is considered in this paper. The method involves solving one small penalty Navier–Stokes problem on a coarse mesh with mesh size H = ?^{1 / 4}h^{1 / 2}, a large penalty Stokes problem on a fine mesh with mesh size h, where 0 < ? < 1 is a penalty parameter. The method we study provides an approximate solution with the convergence rate of same order as the penalty finite volume solution (u_?h,p_?h), which involves solving one large penalty Navier–Stokes problem on a fine mesh with the same mesh size h. However, our method can save a large amount of computational time. Copyright © 2013 John Wiley & Sons, Ltd.     

A two‐level finite element method for the stationary Navier‐Stokes equations based on a stabilized local projection

In this article, we propose a two‐level finite element method to analyze the approximate solutions of the stationary Navier‐Stokes equations based on a stabilized local projection. The local projection allows to circumvent the Babuska‐Brezzi condition by using equal‐order finite element pairs. The local projection can be used to stabilize high equal‐order finite element pairs. The proposed method combines the local projection stabilization method and the two‐level method under the assumption of the uniqueness condition. The two‐level method consists of solving a nonlinear equation on the coarse mesh and solving a linear equation on fine mesh. The nonlinear equation is solved by the one‐step Newtonian iteration method. In the rest of this article, we show the error analysis of the lowest equal‐order finite element pair and provide convergence rate of approximate solutions. Furthermore, the numerical illustrations coincide with the theoretical analysis expectations. From the view of computational time, the results show that the two‐level method is effective to solve the stationary Navier‐Stokes equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Superconvergence for a time‐discretization procedure for the mixed finite element approximation of miscible displacement in porous media

A time‐stepping procedure is established and analyzed for the problem of miscible displacement in porous medium. The fluid velocity based on mixed element method is post‐processed by interpolation along Gauss lines. Error analysis shows that this procedure has a superconvergence error bound O(h) with respect to the parameter h_p, which is one order higher than the conventional Galerkin or mixed finite element procedures. Compared with the results in the references, the parameter constraint conditions in this article are obviously relaxed. © 2011Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

On pointwise decay rates of time‐periodic solutions to the Navier–Stokes equation

We study the existence of a time‐periodic solution with pointwise decay properties to the Navier–Stokes equation in the whole space. We show that if the time‐periodic external force is sufficiently small in an appropriate sense, then there exists a time‐periodic solution $\{u,p\}$ of the Navier–Stokes equation such that $|{?}^j{u(t,x)|=O(|x|}^{1?n?j})$ and $|{?}^j{p(t,x)|=O(|x|}^{?n?j})(j=0,1,\dots )$ uniformly in $t\in \mathbb{R}$ as $|x|\to \infty$. Our solution decays faster than the time‐periodic Stokes fundamental solution and the faster decay of its spatial derivatives of higher order is also described.     

Numerical solution of two‐dimensional parabolic equation subject to nonstandard boundary specifications using the pseudospectral Legendre method

In this research, the problem of solving the two‐dimensional parabolic equation subject to a given initial condition and nonlocal boundary specifications is considered. A technique based on the pseudospectral Legendre method is proposed for the numerical solution of the studied problem. Several examples are given and the numerical results are shown to demonstrate the efficiently of the newly proposed method. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

An analysis for the compressible Stokes equations by first‐order system of least‐squares finite element method

This article applies the first‐order system least‐squares (fosls) finite element method developed by Cai, Manteuffel and McCormick to the compressible Stokes equations. By introducing a new dependent velocity flux variable, we recast the compressible Stokes equations as a first‐order system. Then it is shown that the ellipticity and continuity hold for the least‐squares functionals employing the mixture of H^?1 and L², so that the fosls finite element methods yield best approximations for the velocity flux and velocity. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:689–699, 2001     

Asymptotic preserving time‐discretization of optimal control problems for the Goldstein–Taylor model

We consider the development of implicit‐explicit time integration schemes for optimal control problems governed by the Goldstein–Taylor model. In the diffusive scaling, this model is a hyperbolic approximation to the heat equation. We investigate the relation of time integration schemes and the formal Chapman–Enskog‐type limiting procedure. For the class of stiffly accurate implicit–explicit Runge–Kutta methods, the discrete optimality system also provides a stable numerical method for optimal control problems governed by the heat equation. Numerical examples illustrate the expected behavior. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1770–1784, 2014     

Analysis of a combined mixed finite element and discontinuous Galerkin method for incompressible two‐phase flow in porous media

We analyze a combined method consisting of the mixed finite element method for pressure equation and the discontinuous Galerkin method for saturation equation for the coupled system of incompressible two‐phase flow in porous media. The existence and uniqueness of numerical solutions are established under proper conditions by using a constructive approach. Optimal error estimates in L²(H¹) for saturation and in L ^∞ (H(div)) for velocity are derived. Copyright © 2013 John Wiley & Sons, Ltd.     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

On the uniqueness of the solution of the two‐dimensional Navier–Stokes equation with a Dirac mass as initial vorticity

We propose two different proofs of the fact that Oseen's vortex is the unique solution of the two‐dimensional Navier–Stokes equation with a Dirac mass as initial vorticity. The first argument, due to C. E. Wayne and the second named author, is based on an entropy estimate for the vorticity equation in self‐similar variables. The second proof is new and relies on symmetrization techniques for parabolic equations. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical solutions of the coupled Burgers’ equation by the Galerkin quadratic B‐spline finite element method

In this paper, a coupled Burgers’ equation has been numerically solved by a Galerkin quadratic B‐spline FEM. The performance of the method has been examined on three test problems. Results obtained by the method have been compared with known exact solution and other numerical results in the literature. A Fourier stability analysis of the method is also investigated. Copyright © 2013 John Wiley & Sons, Ltd.     

Identification of the thermal growth characteristics of coagulated tumor tissue in laser‐induced thermotherapy

We consider an inverse problem arising in laser‐induced thermotherapy, a minimally invasive method for cancer treatment, in which cancer tissues are destroyed by coagulation. For the dosage planning, numerical simulations play an important role. To this end, a crucial problem is to identify the thermal growth kinetics of the coagulated zone. Mathematically, this problem is a nonlinear and nonlocal parabolic heat source inverse problem. The solution to this inverse problem is defined as the minimizer of a nonconvex cost functional in this paper. The existence of the minimizer is proven. We derive the Gateaux derivative of the cost functional, which is based on the adjoint system, and use it for a numerical approximation of the optimal coefficient. Numerical implementations are presented to show the validity of the optimization schemes. Copyright © 2012 John Wiley & Sons, Ltd.     

Single cell discretization of O(kh2 + h4) for the estimates of ${\partial u}\over{\partial n}$ for the two‐space dimensional quasi‐linear parabolic equation

In this article, new stable two‐level explicit difference methods of O(kh² + h⁴) for the estimates of for the two‐space dimensional quasi‐linear parabolic equation are derived, where k > 0 and h > 0 are grid sizes in time and space directions, respectively. We use a single computational cell for the methods, which are applicable to the problems both in cartesian and polar coordinates. The proposed methods have the simplicity in nature and employ the same marching type technique of solution. Numerical results obtained by the proposed methods for several different problems were compared with the exact solutions. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 250–261, 2001     

Existence,blow‐up,and exponential decay estimates for a system of semilinear wave equations associated with the helicalflows of Maxwell fluid

The paper is devoted to the study of a system of semilinear wave equations associated with the helical flows of Maxwell fluid. First, based on Faedo–Galerkin method and standard arguments of density corresponding to the regularity of initial conditions, we establish two local existence theorems of weak solutions. Next, we prove that any weak solutions with negative initial energy will blow up in finite time. Finally, we give a sufficient condition to guarantee the global existence and exponential decay of weak solutions via the construction of a suitable Lyapunov functional. Copyright © 2015 John Wiley & Sons, Ltd.