The fictitious domain method with L2‐penalty for the Stokes problem with the Dirichlet boundary condition

We consider the fictitious domain method with L²‐penalty for the Stokes problem with the Dirichlet boundary condition. First, we investigate the error estimates for the penalty method at the continuous level. We obtain the convergence of order in H¹‐norm for the velocity and in L²‐norm for the pressure, where is the penalty parameter. The L²‐norm error estimate for the velocity is upgraded to . Moreover, we derive the a priori estimates depending on for the solution of the penalty problem. Next, we apply the finite element approximation to the penalty problem using the P1/P1 element with stabilization. For the discrete penalty problem, we prove the error estimate in H¹‐norm for the velocity and in L²‐norm for the pressure, where h denotes the discretization parameter. For the velocity in L²‐norm, the convergence rate is improved to . The theoretical results are verified by the numerical experiments.     

An optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows

This article is concerned about an optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows. Using the method, an classical domain decomposition problem is transformed into a constrained minimization problem for which the objective functional is chosen to measure the jump in the dependent variables across the common interfaces between subdomains. The Lagrange multiplier rule is used to transform the constrained optimization problem into an unconstrained one and that rule is applied to derive an optimality system from which optimal solutions may be obtained. The optimality system is also derived using “sensitivity” derivatives instead of the Lagrange multiplier rule. We consider a gradient‐type approach to the solution of domain decomposition problem. The results of some numerical experiments are presented to demonstrate the feasibility and applicability of the algorithm developed in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Composite penalty method of a low order anisotropic nonconforming finite element for the Stokes problem

Composite penalty method of a low order anisotropic nonconforming quadrilateral finite element for the Stokes problem is presented. This method with a large penalty parameter can achieve the same accuracy as the stand method with a small penalty parameter and the convergence rate of this method is two times as that of the standard method under the condition of the same order penalty parameter. The superconvergence for velocity is established as well. The results of this paper are also valid to the most of the known nonconforming finite element methods.     

Two-level iteration penalty and variational multiscale method for steady incompressible flows

In this paper, we study two-level iteration penalty and variational multiscale method for the approximation of steady Navier-Stokes equations at high Reynolds number. Comparing with classical penalty method, this new method does not require very small penalty parameter $\varepsilon$. Moreover, two-level mesh method can save a large amount of CPU time. The error estimates in $H^1$ norm for velocity and in $L^2$ norm for pressure are derived. Finally, two numerical experiments are shown to support the efficiency of this new method.     

On approximation in the fictitious domain method for a bidimensional transport equation

In this article a numerical method for solving a two‐dimensional transport equation in the stationary case is presented. Using the techniques of the variational calculus, we find the approximate solution for a homogeneous boundary‐value problem that corresponds to a square domain D₂. Then, using the method of the fictitious domain, we extend our algorithm to a boundary value problem for a set D that has an arbitrary shape. In this approach, the initial computation domain D (called physical domain) is immersed in a square domain D₂. We prove that the solution obtained by this method is a good approximation of the exact solution. The theoretical results are verified with the help of a numerical example. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A three-step Oseen-linearized finite element method for incompressible flows with damping

This paper introduces a three-step Oseen-linearized finite element method for the 2D/3D steady incompressible Navier–Stokes equations with nonlinear damping term. Within this method, we solve a nonlinear problem over a coarse grid followed by solving two Oseen-linearized problems over a fine grid, which possess the same stiffness matrices with only various right-hand sides. We theoretically analyze the stability of the present method, and derive optimal error estimates of the finite element solutions. We conduct a series of numerical experiments which support the theoretical analysis and test the effectiveness of the proposed method. We demonstrate numerically that there is a significant improvement in the accuracy of the approximate solutions over those for the standard two-level method.     

On the accuracy of the viscous form in simulations of incompressible flow problems

We present a numerical study of drag/lift and flux estimates using two forms of Navier‐Stokes equations (NSE) that are equivalent in the continuum formulation but not in the discrete finite element formulation. The two investigated forms of the NSE differ in the viscous term that is represented in one form by νΔ u with ν being the viscosity and 2ν?·?^S u in the other form where ?^S represents the deformation tensor. The study consists of numerical analysis of the two forms and computations of drag/lift, pressure drop on the cylinder problem and computations of flux for the Poiseuille flow. The main objective is to provide a clear comparison of the reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time for the two forms of NSEs. Our computational results of the reference values do not differ significantly between the two forms, but the differences are there. For the Poiseuille flow, the differences in the flux computations were much smaller, and this agreed with the computationally obtained results of the divergence of the velocity field. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 523–541, 2012     

A stabilized finite element method for a fictitious domain problem allowing small inclusions

The purpose of this work is to approximate numerically an elliptic partial differential equation posed on domains with small perforations (or inclusions). The approach is based on the fictitious domain method, and as the method's interest lies in the case in which the geometrical features are not resolved by the mesh, we propose a stabilized finite element method. The stabilization term is a simple, non‐consistent penalization that can be linked to the Barbosa‐Hughes approach. Stability and convergence are proved, and numerical results confirm the theory.     

The Method of Fictitious Domains in the Signorini Problem

We justify the method of fictitious domains for an elliptic equation with nonlinear Signorini boundary conditions. The method makes it possible to construct a family of auxiliary problems defined in a wider domain and possessing the property that their solutions converge in an appropriate sense to a solution of the original problem.     

A balancing domain decomposition by constraints preconditioner for a weakly over‐penalized symmetric interior penalty method

We develop a balancing domain decomposition by constraints preconditioner for a weakly over‐penalized symmetric interior penalty method for second‐order elliptic problems. We show that the condition number of the preconditioned system satisfies similar estimates as those for conforming finite element methods. Corroborating numerical results are also presented. Copyright © 2012 John Wiley & Sons, Ltd.     

Study of the mixed finite volume method for Stokes and Navier‐Stokes equations

We present finite volume schemes for Stokes and Navier‐Stokes equations. These schemes are based on the mixed finite volume introduced in (Droniou and Eymard, Numer Math 105 (2006), 35‐71), and can be applied to any type of grid (without “orthogonality” assumptions as for classical finite volume methods) and in any space dimension. We present numerical results on some irregular grids, and we prove, for both Stokes and Navier‐Stokes equations, the convergence of the scheme toward a solution of the continuous problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A two‐level method in space and time for the Navier‐Stokes equations

A two‐level method in space and time for the time‐dependent Navier‐Stokes equations is considered in this article. The approximate solution u_M∈H_M is decomposed into the large eddy component v∈H_m(m < M) and the small eddy component w∈H. We obtain the large eddy component v by solving a standard Galerkin equation in a coarse‐level subspace H_m with a time step length k, whereas the small eddy component w is derived by solving a linear equation in an orthogonal complement subspace H with a time step length pk, where p is a positive integer. The analysis shows that our two‐level scheme has long‐time stability and can reach the same accuracy as the standard Galerkin method in fine‐level subspace H_M for an appropriate configuration of p and m. Moreover, some numerical examples are provided to complement our theoretical analysis. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Numerical simulations for two‐dimensional stochastic incompressible Navier‐Stokes equations

This article mainly concerns modeling the stochastic input and its propagation in incompressible Navier‐Stokes(N‐S) flow simulations. The stochastic input is represented spectrally by employing orthogonal polynomial functionals from the Askey scheme as trial basis to represent the random space. A standard Galerkin projection is applied in the random dimension to derive the equations in the weak form. The resulting set of deterministic equations is then solved with standard methods to obtain the mean solution. In this article, the main method employs the Hermite polynomial as the basis in random space. Numerical examples are given and the error analysis is demonstrated for a model problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Robust penalty function method for an uncertain multi-time control optimization problems

This paper gives some new results on multi-time first-order PDE constrained control optimization problem in the face of data uncertainty (MCOPU). We obtain the robust sufficient optimality conditions for (MCOPU). Further, we construct an unconstrained multi-time control optimization problem (MCOPU)_? corresponding to (MCOPU) via absolute value penalty function method. Then, we show that the robust optimal solution to the constrained problem and a robust minimizer to the unconstrained problem are equivalent under suitable hypotheses. Moreover, we give some non-trivial examples to validate the results established in this paper.     

A fast finite difference method for biharmonic equations on irregular domains and its application to an incompressible Stokes flow

Biharmonic equations have many applications, especially in fluid and solid mechanics, but is difficult to solve due to the fourth order derivatives in the differential equation. In this paper a fast second order accurate algorithm based on a finite difference discretization and a Cartesian grid is developed for two dimensional biharmonic equations on irregular domains with essential boundary conditions. The irregular domain is embedded into a rectangular region and the biharmonic equation is decoupled to two Poisson equations. An auxiliary unknown quantity Δu along the boundary is introduced so that fast Poisson solvers on irregular domains can be used. Non-trivial numerical examples show the efficiency of the proposed method. The number of iterations of the method is independent of the mesh size. Another key to the method is a new interpolation scheme to evaluate the residual of the Schur complement system. The new biharmonic solver has been applied to solve the incompressible Stokes flow on an irregular domain.      

NONCONFORMING FINITE ELEMENT PENALTY METHOD FOR STOKES EQUATION

a special penalty method is presented to improve the accuracy of the standard penaltymethod (or solving Stokes equation with nonconforming finite element, It is shown that thismethod with a larger penalty parameter can achieve the same accuracy as the staodaxd methodwith a smaller penalty parameter. The convergence rate of the standard method is just hall order of this penalty method when using the same penalty parameter, while the extrapolationmethod proposed by Faik et al can not yield so high accuracy of convergence. At last, we alsoget the super-convergence estimates for total flux.     

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.