Spectral method for the unsteady incompressible Navier-Stokes equations in gauge formulation

A spectral method which uses a gauge method, as opposed to a projection method, to decouple the computation of velocity and pressure in the unsteady incompressible Navier-Stokes equations, is presented. Gauge methods decompose velocity into the sum of an auxilary field and the gradient of a gauge variable, which may, in principle, be assigned arbitrary boundary conditions, thus overcoming the issue of artificial pressure boundary conditions in projection methods. A lid-driven cavity flow is used as a test problem. A subtraction method is used to reduce the pollution effect of singularities at the top corners of the cavity. A Chebyshev spectral collocation method is used to discretize spatially. An exponential time differencing method is used to discretize temporally. Matrix diagonalization procedures are used to compute solutions directly and efficiently. Numerical results for the flow at Reynolds number Re = 1000 are presented, and compared to benchmark results. It is shown that the method, called the spectral gauge method, is straightforward to implement, and yields accurate solutions if Neumann boundary conditions are imposed on the gauge variable, but suffers from reduced convergence rates if Dirichlet boundary conditions are imposed on the gauge variable

Drop Formation in a One-Dimensional Approximation of the Navier-Stokes Equation

We consider the viscous motion of a thin, axisymmetric column of fluid with a free surface. A one-dimensional equation of motion for the velocity and the radius is derived from the Navier-Stokes equation. We compare with recent experiments on the breakup of a liquid jet and on the bifurcation of a drop suspended from an orifice. The equations form singularities as the fluid neck is pinching off. The nature of the singularities is investigated in detail

Adjoint recovery of superconvergent functionals from PDE approximations

Motivated by applications in computational fluid dynamics, a method is presented for obtaining estimates of integral functionals, such as lift or drag, that have twice the order of accuracy of the computed flow solution on which they are based. This is achieved through error analysis that uses an adjoint PDE to relate the local errors in approximating the flow solution to the corresponding global errors in the functional of interest. Numerical evaluation of the local residual error together with an approximate solution to the adjoint equations may thus be combined to produce a correction for the computed functional value that yields the desired improvement in accuracy. Numerical results are presented for the Poisson equation in one and two dimensions and for the nonlinear quasi-one-dimensional Euler equations. The theory is equally applicable to nonlinear equations in complex multi-dimensional domains and holds great promise for use in a range of engineering disciplines in which a few integral quantities are a key output of numerical approximations

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement distributions in various bodies containing cracks. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The results obtained show a considerable potential for using this method in the three-dimensional analysis of finite geometry solids and suggest a possible extension of this technique to nonlinear material behavior

A fast high-order solver for problems of scattering by heterogeneous bodies

A new high-order integral algorithm for the solution of scattering problems by heterogeneous bodies is presented. Here, a scatterer is described by a (continuously or discontinuously) varying refractive index n(x) within a two-dimensional (2D) bounded region; solutions of the associated Helmholtz equation under given incident fields are then obtained by high-order inversion of the Lippmann-Schwinger integral equation. The algorithm runs in O(Nlog(N)) operations where N is the number of discretization points. A wide variety of numerical examples provided include applications to highly singular geometries, high-contrast configurations, as well as acoustically/electrically large problems for which supercomputing resources have been used recently. Our method provides highly accurate solutions for such problems on small desktop computers in CPU times of the order of seconds

Numerical Relativity: A review

Computer simulations are enabling researchers to investigate systems which are extremely difficult to handle analytically. In the particular case of General Relativity, numerical models have proved extremely valuable for investigations of strong field scenarios and been crucial to reveal unexpected phenomena. Considerable efforts are being spent to simulate astrophysically relevant simulations, understand different aspects of the theory and even provide insights in the search for a quantum theory of gravity. In the present article I review the present status of the field of Numerical Relativity, describe the techniques most commonly used and discuss open problems and (some) future prospects.Comment: 2 References added; 1 corrected. 67 pages. To appear in Classical and Quantum Gravity. (uses iopart.cls