Hyperasymptotics and hyperterminants: Exceptional cases

A new method is introduced for the computation of hyperterminants. It is based on recurrence relations, and can also be used to compute the parameter derivatives of the hyperterminants. These parameter derivatives are needed in hyperasymptotic expansions in exceptional cases. Numerical illustrations and an application are included.     

An efficient transient Navier–Stokes solver on compact nonuniform space grids

In this paper, we propose an implicit higher-order compact (HOC) finite difference scheme for solving the two-dimensional (2D) unsteady Navier–Stokes (N–S) equations on nonuniform space grids. This temporally second-order accurate scheme which requires no transformation from the physical to the computational plane is at least third-order accurate in space, which has been demonstrated with numerical experiments. It efficiently captures both transient and steady-state solutions of the N–S equations with Dirichlet as well as Neumann boundary conditions. The proposed scheme is likely to be very useful for the computation of transient viscous flows involving free and wall bounded shear layers which invariably contain spatial scale variation. Numerical results are presented and compared with analytical as well as established numerical data. Excellent comparison is obtained in all the cases.     

Effect of boundary vorticity discretization on explicit stream‐function vorticity calculations

The numerical solution of the time‐dependent Navier–Stokes equations in terms of the vorticity and a stream function is a well tested process to describe two‐dimensional incompressible flows, both for fluid mixing applications and for studies in theoretical fluid mechanics. In this paper, we consider the interaction between the unsteady advection–diffusion equation for the vorticity, the Poisson equation linking vorticity and stream function and the approximation of the boundary vorticity, examining from a practical viewpoint, global iteration stability and error. Our results show that most schemes have very similar global stability constraints although there may be small stability gains from the choice of method to determine boundary vorticity. Concerning accuracy, for one model problem we observe that there were cases where the boundary vorticity discretization did not propagate to the interior, but for the usual cavity flow all the schemes tested had error close to second order. Copyright © 2005 John Wiley & Sons, Ltd.     

Remarks on the links between low‐order DG methods and some finite‐difference schemes for the Stokes problem

In this paper we demonstrate that some well‐known finite‐difference schemes can be interpreted within the framework of the local discontinuous Galerkin (LDG) methods using the low‐order piecewise solenoidal discrete spaces introduced in (SIAM J. Numer. Anal. 1990; 27 (6): 1466–1485). In particular, it appears that it is possible to derive the well‐known MAC scheme using a first‐order Nédélec approximation on rectangular cells. It has been recently interpreted within the framework of the Raviart–Thomas approximation by Kanschat (Int. J. Numer. Meth. Fluids 2007; published online). The two approximations are algebraically equivalent to the MAC scheme, however, they have to be applied on grids that are staggered on a distance h/2 in each direction. This paper also demonstrates that both discretizations allow for the construction of a divergence‐free basis, which yields a linear system with a ‘biharmonic’ conditioning. Both this paper and Kanschat (Int. J. Numer. Meth. Fluids 2007; published online) demonstrate that the LDG framework can be used to generalize some popular finite‐difference schemes to grids that are not parallel to the coordinate axes or that are unstructured. Copyright © 2008 John Wiley & Sons, Ltd.     

Stability of some low‐order approximations for the Stokes problem

Two‐level low‐order finite element approximations are considered for the inhomogeneous Stokes equations. The elements introduced are attractive because of their simplicity and computational efficiency. In this paper, the stability of a Q₁(h)–Q₁(2h) approximation is analysed for general geometries. Using the macroelement technique, we prove the stability condition for both two‐ and three‐dimensional problems. As a result, optimal rates of convergence are found for the velocity and pressure approximations. Numerical results for three test problems are presented. We observe that for the computed examples, the accuracy of the two‐level bilinear approximation is compared favourably with some standard finite elements. Copyright © 2007 John Wiley & Sons, Ltd.     

Solution to transient Navier–Stokes equations by the coupling of differential quadrature time integration scheme with dual reciprocity boundary element method

The two‐dimensional time‐dependent Navier–Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equation, which links stream function and vorticity with an initial vorticity guess, produces velocity components in turn for the solution to vorticity transport equation. The DRBEM formulation of the vorticity transport equation results in an initial value problem represented by a system of first‐order ordinary differential equations in time. When the DQM discretizes this system in time direction, we obtain a system of linear algebraic equations, which gives the solution vector for vorticity at any required time level. The procedure outlined here is also applied to solve the problem of two‐dimensional natural convection in a cavity by utilizing an iteration among the stream function, the vorticity transport and the energy equations as well. The test problems include two‐dimensional flow in a cavity when a force is present, the lid‐driven cavity and the natural convection in a square cavity. The numerical results are visualized in terms of stream function, vorticity and temperature contours for several values of Reynolds (Re) and Rayleigh (Ra) numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

The steady Navier–Stokes/energy system with temperature‐dependent viscosity—Part 1: Analysis of the continuous problem

In this first part we propose and analyse a model for the study of two‐dimensional incompressible Navier–Stokes equations with a temperature‐dependent viscosity. The flow is supposed in a mixed convection regime and considers an outflow region, leading to a strongly coupled problem between the Navier–Stokes and energy equations, which will be justified theoretically. The coupling in the continuous problem is treated by an outer temperature fixed point strategy. Existence results for a particular variational formulation follows from this study. Further, a particular uniqueness result for small data is also obtained. Copyright © 2007 John Wiley & Sons, Ltd.     

A note on the second problem of Stokes for Maxwell fluids

New exact solutions corresponding to the second problem of Stokes for Maxwell fluids have been established by means of Laplace transforms. For large times, these solutions reduce to the well-known steady-state solutions which are periodic in time and independent of the initial conditions. Furthermore, the transient solutions are in accordance with the previous solutions obtained using the Fourier sine transform. The required time to get the steady-state is determined by graphical illustrations. This time decreases if the frequency of the velocity increases. The effects of the material parameters on the decay of the transients in time are also investigated by graphs.     

Shape reconstruction of an inverse Stokes problem

This paper deals with the shape reconstruction of a viscous incompressible fluid driven by the Stokes flow. For the approximate solution of the ill-posed and nonlinear problem we propose a regularized Newton method. A theoretical foundation for the Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the interior boundary curve in the sense of a domain derivative. The numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible.