accurate projection methods for the incompressible navier stokes equations abstract this paper considers the accuracy of projection method approximations to the initial boundary value problem for the incompressible navier stokes equations the issue of how to correctly specify numerical boundary conditions for these methods has been outstanding since the birth of the second order methodology a decade and a half ago it has been observed that while the velocity can be reliably computed to second order accuracy in time and space the pressure is typically only first order accurate in the l norm this paper identifies the source of this problem in the interplay of the global pressure update formula with the numerical boundary conditions and presents an improved projection algorithm which is fully second order accurate as demonstrated by a normal mode analysis and numerical experiments in addition a numerical method based on a gauge variable formulation of the incompressible navier stokes equations which provides another option for obtaining fully second order convergence in both velocity and pressure is discussed the connection between the boundary conditions for projection methods and the gauge method is explained in detail

zeroes of zeta functions and symmetry hilbert and polya suggested that there might be a natural spectral interpretation of the zeroes of the riemann zeta function while at the time there was little evidence for this today the evidence is quite convincing firstly there are the function field analogues that is zeta functions of curves over finite fields and their generalizations for these a spectral interpretation for their zeroes exists in terms of eigenvalues of frobenius on cohomology secondly the developments both theoretical and numerical on the local spacing distributions between the high zeroes of the zeta function and its generalizations give striking evidence for such a spectral connection moreover the low lying zeroes of various families of zeta functions follow laws for the eigenvalue distributions of members of the classical groups in this paper we review these developments in order to present the material fluently we do not proceed in chronological order of discovery also in concentrating entirely on the subject matter of the title we are ignoring the standard body of important work that has been done on the zeta function and l functions 1 the montgomery odlyzko law we begin with the riemann zeta function and some phenomenology associated with it

fast and high accuracy multigrid solution of the three dimensional poisson equation we employ a fourth order compact finite difference scheme fos with the multigrid algorithm to solve the three dimensional poisson equation we test the influence of different orderings of the grid space and different grid transfer operators on the convergence and efficiency of our high accuracy algorithm fourier smoothing analysis is conducted to show that fos has a smaller smoothing factor than the traditional second order central difference scheme cds a new method of fourier smoothing analysis is proposed for the partially decoupled red black gauss seidel relaxation with fos numerical results are given to compare the computed accuracy and the computational efficiency of fos with multigrid against cds with multigrid

a finite difference galerkin formulation for the incompressible navier stokes equations abstract the development of a new computational method for solving the incompressible navier stokes equations in primitive variable form is presented it is found that certain finite difference approximations for these equations can be transformed into an equivalent system which efficiently determines the discrete velocity field and which completely eliminates the pressure two such difference schemes for two dimensional problems are examined and some preliminary numerical results are discussed for the steady driven cavity problem

a coarse grid projection method for accelerating incompressible flow computations we present a coarse grid projection cgp method for accelerating incompressible flow computations which is applicable to methods involving poisson equations as incompressibility constraints the cgp methodology is a modular approach that facilitates data transfer with simple interpolations and uses black box solvers for the poisson and advection diffusion equations in the flow solver after solving the poisson equation on a coarsened grid an interpolation scheme is used to obtain the fine data for subsequent time stepping on the full grid a particular version of the method is applied here to the vorticity stream function primitive variable and vorticity velocity formulations of incompressible navier stokes equations we compute several benchmark flow problems on two dimensional cartesian and non cartesian grids as well as a three dimensional flow problem the method is found to accelerate these computations while retaining a level of accuracy close to that of the fine resolution field which is significantly better than the accuracy obtained for a similar computation performed solely using a coarse grid a linear acceleration rate is obtained for all the cases we consider due to the linear cost elliptic poisson solver used with reduction factors in computational time between 2 and 42 the computational savings are larger when a suboptimal poisson solver is used we also find that the computational savings increase with increasing distortion ratio on non cartesian grids making the cgp method a useful tool for accelerating generalized curvilinear incompressible flow solvers

effects of the computational time step on numerical solutions of turbulent flow effects of large computational time steps on the computed turbulence were investigated using a fully implicit method in turbulent channel flow computations the largest computational time step in wall units which led to accurate prediction of turbulence statistics was determined turbulence fluctuations could not be sustained if the computational time step was near or larger than the kolmogorov time scale

entropy of quantum limits in this paper we show that any measure arising as a weak limit of microlocal lifts of eigenfunctions of the laplacian on certain arithmetic manifolds have dimension at least 11 9 and in particular all ergodic components of this measure with respect to the geodesic flow have positive entropy

application of a fractional step method to incompressible navier stokes equations a numerical method for computing three dimensional time dependent incompressible flows is presented the method is based on a fractional step or time splitting scheme in conjunction with the approximate factorization technique it is shown that the use of velocity boundary conditions for the intermediate velocity field can lead to inconsistent numerical solutions appropriate boundary conditions for the intermediate velocity field are derived and tested numerical solutions for flows inside a driven cavity and over a backward facing step are presented and compared with experimental data and other numerical results

the accuracy of the fractional step method we analyze the accuracy of the fractional step method of kim and moin j comput phys 59 1985 pp 308 323 for the incompressible navier stokes equations we show that the boundary conditions cannot be exactly satisfied in the projection step and that this limits the accuracy of the method we also show that the pressure in any projection method can be at best first order accurate our analysis is simpler and more direct than the previous analyses of this method we also show that there is no numerical boundary layer for velocity or pressure but there is one for the auxiliary pressure variable

stable and accurate pressure approximation for unsteady incompressible viscous flow how to properly specify boundary conditions for pressure is a longstanding problem for the incompressible navier stokes equations with no slip boundary conditions an analytical resolution of this issue stems from a recently developed formula for the pressure in terms of the commutator of the laplacian and leray projection operators here we make use of this formula to a improve the accuracy of computing pressure in two kinds of existing time discrete projection methods implicit in viscosity only and b devise new higher order accurate time discrete projection methods that extend a slip correction idea behind the well known finite difference scheme of kim and moin we test these schemes for stability and accuracy using various combinations of c 0 finite elements for all three kinds of time discretization one can obtain third order accuracy for both pressure and velocity without a time step stability restriction of diffusive type furthermore two kinds of projection methods are found stable using piecewise linear elements for both velocity and pressure

an unsplit higher order godunov method for scalar conservation laws in multiple dimensions abstract in this paper we develop an unsplit higher order godunov method for scalar conservation laws in two dimensions the method represents an extension for the special case being considered of methods developed by colella and by van leer in our method we begin with a piecewise bilinear representation of the solution on each grid cell a new piecewise bilinear representation at the next time level is then obtained from a two step procedure in the first step a conservative predictor corrector scheme derived from the integral form of the differential equation and characteristic considerations is used to obtain average values of the solution over grid cells at the new time level next these new average values are then used to construct a limited piecewise bilinear profile for each cell at the new time level the resulting method is shown to satisfy a maximum principle for constant coefficient linear advection computational results are presented comparing the new method to colella s method for linear advection the method is also applied to two model problems from porous media flow miscible and immiscible displacement the new scheme provides accurate resolution of sharp fronts without any significant distortion

a quadratic spline based interface quasi reconstruction algorithm for accurate tracking of two phase flows a new quadratic spline based interface quasi reconstruction algorithm is presented which provides an accurate and continuous representation of the interface in a multiphase domain and facilitates the direct estimation of local interfacial curvature the fluid interface in each of the mixed cells is represented by piecewise parabolic curves and an initial discontinuous plic approximation of the interface is progressively converted into a smooth quadratic spline made of these parabolic curves the conversion is achieved by a sequence of predictor corrector operations enforcing function c 0 and derivative c 1 continuity at the cell boundaries using simple analytical expressions for the continuity requirements the efficacy and accuracy of the current algorithm has been demonstrated using standard test cases involving reconstruction of known static interface shapes and dynamically evolving interfaces in prescribed flow situations these benchmark studies illustrate that the present algorithm performs excellently as compared to the other interface reconstruction methods available in literature quadratic rate of error reduction with respect to grid size has been observed in all the cases with curved interface shapes only in situations where the interface geometry is primarily flat the rate of convergence becomes linear with the mesh size the flow algorithm implemented in the current work is designed to accurately balance the pressure gradients with the surface tension force at any location as a consequence it is able to minimize spurious flow currents arising from imperfect normal stress balance at the interface this has been demonstrated through the standard test problem of an inviscid droplet placed in a quiescent medium finally the direct curvature estimation ability of the current algorithm is illustrated through the coupled multiphase flow problem of a deformable air bubble rising through a column of water

application of time iterative schemes to incompressible flow the computation of steady incompressible flows by an euler implicit algorithm is studied using both the incompressible equations and the low mach number compressible equations the incompressible equations are handled by adding an artificial time derivative to the continuity equation this allows both the pressure and velocity to be obtained implicitly in one dimensional problems both systems converge rapidly even at low mach numbers where the eigenvalues are very stiff in two dimensions where approximate factorization is required the presence of stiff eigenvalues is highly detrimental stiffness can be avoided in the incompressible equations by selecting an appropriate pseudo mach number this insures reliable convergence and results in an efficient incompressible flow algorithm in the case of compressible equations the mach number cannot be chosen arbitrarily and the contamination introduced by approximate factorization must be removed a matrix preconditioning factor that accomplishes this is developed and demonstrated with this modification the convergence rate is the same as in the incompressible case and is independent of mach number rapid convergence is observed at mach numbers as low as 6 05

a cartesian grid embedded boundary method for hyperbolic conservation laws we present a second order godunov algorithm to solve time dependent hyperbolic systems of conservation laws on irregular domains our approach is based on a formally consistent discretization of the conservation laws on a finite volume grid obtained from intersecting the domain with a cartesian grid we address the small cell stability problem associated with such methods by hybridizing our conservative discretization with a stable nonconservative discretization at irregular control volumes and redistributing the difference in the mass increments to nearby cells in a way that preserves stability and local conservation the resulting method is second order accurate in l 1 for smooth problems and is robust in the presence of large amplitude discontinuities intersecting the irregular boundary

a gradient augmented level set method for unstructured grids we present and test a gradient augmented level set method for unstructured grids which has been designed to allow for easy integration in parallel flow solvers the means to this end is a hermite formulation where the gradient is stored and advected as an independent variable along with the level set function the method uses narrow band storage in conjunction with a fast marching method which is also designed to take advantage of the gradient augmentation we demonstrate through standard passive advection test cases that the proposed method is able to compute interface motion with performance comparable to popular high performing methods on structured grids

on convergence of block centered finite differences for elliptic problems we consider linear selfadjoint elliptic problems with neumann boundary conditions in rectangular domains we demonstrate that with sufficiently smooth data the discrete l 2 norm errors for tensor product block centered finite differences in both the approximate solution and its first derivatives are second order for all nonuniform grids extensions to nonselfadjoint and parabolic problems are discussed

pressure stability in fractional step finite element methods for incompressible flows abstract the objective of this paper is to analyze the pressure stability of fractional step finite element methods for incompressible flows that use a pressure poisson equation for the classical first order projection method it is shown that there is a pressure control which depends on the time step size and therefore there is a lower bound for this time step for stability reasons the situation is much worse for a second order scheme in which part of the pressure gradient is kept in the momentum equation the pressure stability in this case is extremely weak to overcome these shortcomings a stabilized fractional step finite element method is also considered and its stability is analyzed some simple numerical examples are presented to support the theoretical results

on the number of nodal domains of random spherical harmonics let n f be a number of nodal domains of a random gaussian spherical harmonic f of degree n we prove that as n grows to infinity the mean of n f n 2 tends to a positive constant a and that n f n 2 exponentially concentrates around a this result is consistent with predictions made by bogomolny and schmit using a percolation like model for nodal domains of random gaussian plane waves

on error estimates of the projection methods for the navier stokes equations second order schemes we present in this paper a rigorous error analysis of several projection schemes for the approximation of the unsteady incompressible navier stokes equations the error analysis is accomplished by interpreting the respective projection schemes as second order time discretizations of a perturbed system which approximates the navier stokes equations numerical results in agreement with the error analysis are also presented

a new family of 5 5 cc 4oc schemes applicable for unsteady navier stokes equations in this paper a new family of implicit compact finite difference schemes for computation of unsteady convection diffusion equation with variable convection coefficient is proposed the schemes which are fourth order accurate in space and second or lower order accurate in time depending on the choice of weighted time average parameter are then applied to unsteady navier stokes system the proposed schemes where transport variable and its first derivatives are carried as the unknowns combine virtues of compact discretization and pade scheme for spatial derivative these schemes which are based on a five point stencil with constant coefficients named as 5 5 constant coefficient 4th order compact 5 5 cc 4oc give rise to a diagonally dominant system of equations and shows higher accuracy and better phase and amplitude error characteristics than some of the standard methods these schemes are capable of using a grid aspect ratio other than unity and are unconditionally stable they efficiently capture both transient and steady solutions of linear and nonlinear convection diffusion equations with dirichlet as well as neumann boundary conditions subsequently the proposed schemes are applied to problems governed by the incompressible navier stokes equations the results obtained are in excellent agreement with analytical and available numerical results in all cases establishing efficiency and accuracy of the proposed schemes

l infty norms of eigenfunctions of arithmetic surfaces here and elsewhere a l1o0 as a xc is a basic problem of quantum chaos s in any case almost nothing beyond 0 2 is known about 110 joc when the curvature is negative one can push the standard techniques and replace a1 4 by a1 4 log a in this case in this paper we use arithmetic techniques in particular modular correspondences to obtain the first improvement in the exponent 0 2 we also obtain lower bounds on the l norms in these arithmetic cases in more detail let a b be a quaternion division algebra over q a is linearly generated by 1 w q wq over q and w2 a q2 b wq qw 0 here a b e z are square free and we will assume that a 0 as usual the norm and trace are defined by n a a d tr a a oy where if a xo x1w x2q x3wq

an accurate second order approximation factorization method for time dependent incompressible navier stokes equations in spherical polar coordinates a finite difference method for solving three dimensional time dependent incompressible navier stokes equations in spherical polar coordinates is presented in detail a new algorithm which is second order accurate in time and space is considered and decoupling between the velocity and the pressure is achieved by this algorithm further the numerical method is tested by computing the spherical couette flow between two concentric spheres with the inner one rotating a comparison of the numerical solutions with available numerical results and experimental measurements is made it is demonstrated that the numerical code is valid for solving three dimensional unsteady incompressible navier stokes equations in spherical polar coordinates

an accurate adaptive solver for surface tension driven interfacial flows a method combining an adaptive quad octree spatial discretisation geometrical volume of fluid interface representation balanced force continuum surface force surface tension formulation and height function curvature estimation is presented the extension of these methods to the quad octree discretisation allows adaptive variable resolution along the interface and is described in detail the method is shown to recover exact equilibrium to machine accuracy between surface tension and pressure gradient in the case of a stationary droplet irrespective of viscosity and spatial resolution accurate solutions are obtained for the classical test case of capillary wave oscillations an application to the capillary breakup of a jet of water in air further illustrates the accuracy and efficiency of the method the source code of the implementation is freely available as part of the gerris flow solver

efficient implementation of thinc scheme in this paper we improve the multi dimensional thinc tangent of hyperbola for interface capturing scheme f xiao y honma t kono a simple algebraic interface capturing scheme using hyperbolic tangent function int j numer meth fluid 48 2005 1023 the thinc scheme is a vof volume of fluid type method in the original thinc scheme one dimensional thinc scheme was straightforwardly used for multi dimensional cases in this paper we propose the wlic weighed line interface calculation method to extend the thinc scheme to multi dimension in the wlic method the interface is reconstructed by taking an average of line interfaces along x y and z coordinates with the weights calculated from surface normal the wlic method can extend the thinc scheme to multi dimension while maintaining simplicity of implementation and achieve a higher accuracy than the original thinc scheme the wlic method can readily extend the thinc scheme to three dimension

syzygy modules for quasi k gorenstein rings abstract let be a right quasi k gorenstein ring for each dth syzygy module m in mod where 0 d k 1 we obtain an exact sequence 0 b m p c 0 in mod with the properties that it is dual exact p is projective c is a d 1 st syzygy module b is a dth syzygy of ext op d 1 tr m and the right projective dimension of b is less than or equal to d 1 we then give some applications of such an exact sequence as follows 1 we obtain a chain of epimorphisms concerning m and by dualizing it we then get the spherical filtration of auslander and bridger for m 2 we get auslander and bridger s approximation theorem for each reflexive module in mod op 3 we show that for any 0 d k 1 each dth syzygy module in mod has an evans griffith presentation as an immediate consequence of 3 we have that if is a commutative noetherian ring with finite self injective dimension then for any non negative integer d each dth syzygy module in mod has an evans griffith presentation which generalizes an evans and griffith s result to much more general setting

arithmetic and equidistribution of measures on the sphere motivated by problems of mathematical physics quantum chaos questions of equidistribution of eigenfunctions of the laplace operator on a riemannian manifold have been studied by several authors we consider here in analogy with arithmetic hyperbolic surfaces orthonormal bases of eigenfunctions of the laplace operator on the two dimensional unit sphere which are also eigenfunctions of an algebra of hecke operators which act on these spherical harmonics we formulate an analogue of the equidistribution of mass conjecture for these eigenfunctions as well as of the conjecture that their moments tend to moments of the gaussian as the eigenvalue increases for such orthonormal bases we show that these conjectures are related to the analytic properties of degree eight arithmetic l functions associated to triples of eigenfunctions moreover we establish the conjecture for the third moments and give a conditional on standard analytic conjectures about these arithmetic l functions proof of the equidistribution of mass conjecture

a high order discontinuous galerkin method for 2d incompressible flows in this paper we introduce a high order discontinuous galerkin method for two dimensional incompressible flow in the vorticity stream function formulation the momentum equation is treated explicitly utilizing the efficiency of the discontinuous galerkin method the stream function is obtained by a standard poisson solver using continuous finite elements there is a natural matching between these two finite element spaces since the normal component of the velocity field is continuous across element boundaries this allows for a correct upwinding gluing in the discontinuous galerkin framework while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability the method is efficient for inviscid or high reynolds number flows optimal error estimates are proved and verified by numerical experiments

comparison of second and fourth order discretizations for multigrid poisson solvers we combine a compact high order difference approximation with multigrid v cycle algorithm to solve the two dimensional poisson equation with dirichlet boundary conditions this scheme along with several different orderings of grid space and projection operators is compared with the five point formula to show the dramatic improvement in computed accuracy on serial and vector machines

the immersed boundary method for advection electrodiffusion with implicit timestepping and local mesh refinement we describe an immersed boundary method for problems of fluid solute structure interaction the numerical scheme employs linearly implicit timestepping allowing for the stable use of timesteps that are substantially larger than those permitted by an explicit method and local mesh refinement making it feasible to resolve the steep gradients associated with the space charge layers as well as the chemical potential which is used in our formulation to control the permeability of the membrane to the possibly charged solute low reynolds number fluid dynamics are described by the time dependent incompressible stokes equations which are solved by a cell centered approximate projection method the dynamics of the chemical species are governed by the advection electrodiffusion equations and our semi implicit treatment of these equations results in a linear system which we solve by gmres preconditioned via a fast adaptive composite grid fac solver numerical examples demonstrate the capabilities of this methodology as well as its convergence properties

nodal domains statistics a criterion for quantum chaos we consider the distribution of the properly normalized numbers of nodal domains of wave functions in 2d quantum billiards we show that these distributions distinguish clearly between systems with integrable separable or chaotic underlying classical dynamics and for each case the limiting distribution is universal system independent thus a new criterion for quantum chaos is provided by the statistics of the wave functions which complements the well established criterion based on spectral statistics

a stable and convergent scheme for viscoelastic flow in contraction channels we present a new algorithm to simulate unsteady viscoelastic flows in abrupt contraction channels in our approach we split the viscoelastic terms of the oldroyd b constitutive equation using duhamel s formula and discretize the resulting pdes using a semi implicit finite difference method based on a lax wendroff method for hyperbolic terms in particular we leave a small residual elastic term in the viscous limit by design to make the hyperbolic piece well posed a projection method is used to impose the incompressibility constraint we are able to compute the full range of unsteady elastic flows in an abrupt contraction channel from the viscous limit to the elastic limit in a stable and convergent manner we demonstrate the range of our method for unsteady flow of a maxwell fluid with and without viscosity in planar contraction channels we also demonstrate stable and convergent results for benchmark high weissenberg number problems at we 1 and we 10

a cartesian grid method for modeling multiple moving objects in 2d incompressible viscous flow we present an efficient method for solving 2d incompressible viscous flows around multiple moving objects our method employs an underlying regular cartesian grid to solve the system using a streamfunction vorticity formulation and with discontinuities representing the embedded objects the no penentration condition for the moving geometry is satisfied by superposing a homogenous solution to the poisson s equation for the streamfunction the no slip condition is satisfied by generating vorticity on the surfaces of the objects both the initial poisson solution and the evaluation of the homogenous solution require embedding irregular discontinuities in a fast poisson solver computation time is dictated by the time required to do a fast poisson solution plus solve an integral form of laplace s equation there is no significant increase in computational cost if the geometry of the embedded objects is variable and moving relative to the underlying grid we test the method against the canonical example of flow past a cylinder and obtained new results on the flow and forces of two cylinders moving relative to each other

a high order immersed interface method for simulating unsteady incompressible flows on irregular domains immersed boundary methods and immersed interface methods are becoming increasingly popular for the computation of unsteady flows around complex geometries using a cartesian grid while good results both qualitative and quantitative have been obtained most of the methods rely on low order corrections to account for the immersed boundary the objective of the present work is to present as an alternative a high order modified immersed interface method for the 2d unsteady incompressible navier stokes equations in stream function vorticity formulation the method employs an explicit fourth order runge kutta time integration scheme fourth order compact finite differences for computation of spatial derivatives and a nine point fourth order compact discretization of the poisson equation for computation of the stream function corrections to the finite difference schemes are used to maintain high formal accuracy at the immersed boundary as confirmed by analytical tests to validate the method in its application to incompressible flows several physically relevant test cases are computed including uniform flow past a circular cylinder and tollmien schlichting waves in a boundary layer

discrete exterior calculus discretization of incompressible navier stokes equations over surface simplicial meshes a conservative discretization of incompressible navier stokes equations is developed based on discrete exterior calculus dec a distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product the governing equations are first rewritten using the exterior calculus notation replacing vector calculus differential operators by the exterior derivative hodge star and wedge product operators the discretization is then carried out by substituting with the corresponding discrete operators based on the dec framework numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured triangular meshes and first order accuracy for otherwise unstructured meshes by construction the method is conservative in that both mass and vorticity are conserved up to machine precision the relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step

efficient algorithm on a nonstaggered mesh for simulating rayleigh benard convection in a box an efficient semi implicit second order accurate finite difference method is described for studying incompressible rayleigh benard convection in a box with sidewalls that are periodic thermally insulated or thermally conducting operator splitting and a projection method reduce the algorithm at each time step to the solution of four helmholtz equations and one poisson equation and these are are solved by fast direct methods the method is numeri cally stable even though all field values are placed on a single non staggered mesh commensurate with the boundaries the efficiency and accuracy of the method are characterized for several representative convection problems

unsteady solution of incompressible navier stokes equations abstract a new time accurate finite difference numerical method for solving incompressible navier stokes equations is presented with primitive variables as the unknowns the numerical scheme is a crank nicolson implicit treatment of all terms of the equations with central differencing for space derivatives the pressure and the nonlinearities in convection terms are not lagged to obtain the incompressible solution at the advanced time level a continuous auxiliary system is introduced in pseudo time with artificial compressibility time accurate solutions are presented for two dimensional fluid flows in a square cavity with an impulsively starting lid and with an oscillating lid

physically motivated domain decomposition for singularly perturbed equations a domain decomposition algorithm suitable for the efficient and accurate solution of a parabolic reaction convection diffusion equation with small parameter multiplying the diffusion term is presented convergence is established via maximum principle arguments the equation arises in the modeling of laminar transonic flow decomposition into subdomains is accomplished via singular perturbation analysis which dictates regions where certain reduced equations may be solved in place of the full equation effectively preconditioning the problem this paper concentrates on the theoretical basis of the method establishing local and global a priori error bounds

a second order projection method for the incompressible navier stokes equations in this paper we describe a second order projection method for the time dependent incom pressible navier stokes equations as in the original projection method developed by chorin we first solve diffusion convection equations to predict intermediate velocities which are then projected onto the space of divergence free vector fields by introducing more coupling between the diffusion onvection step and the projection step we obtain a temporal discretiza tion that is second order accurate our treatment of the diffusion convection step uses a specialized higher order godunov method for differencing the nonlinear convective terms that provides a robust treatment of these terms at high reynolds number the godunov procedure is second order accurate for smooth flow and remains stable for discontinuous initial data even in the zero viscosity limit we approximate the projection directly using a galerkin procedure that uses a local basis for discretely divergence free vector fields numerical results are presented validating the convergence properties of the method we also apply the method to doubly periodic shear layers to assess the performance of the method on more difficult

a multigrid ghost point level set method for incompressible navier stokes equations on moving domains with curved boundaries abstract in this paper we present a numerical approach to solve the navier stokes equations on moving domains with second order accuracy the space discretization is based on the ghost point method which falls under the category of unfitted boundary methods since the mesh does not adapt to the moving boundary the equations are advanced in time by using crank nicholson the momentum and continuity equations are solved simultaneously for the velocity and the pressure by adopting a proper multigrid approach to avoid the checkerboard instability for the pressure a staggered grid is adopted where velocities are defined at the sides of the cell and the pressure is defined at the centre the lack of uniqueness for the pressure is circumvented by the inclusion of an additional scalar unknown representing the average divergence of the velocity and an additional equation to set the average pressure to zero several tests are performed to simulate the motion of an incompressible fluid around a moving object as well as the lid driven cavity tests around steady objects the object is implicitly defined by a level set approach that allows a natural computation of geometrical properties such as distance from the boundary normal directions and curvature different shapes are tested circle ellipse and flower numerical results show the second order accuracy for the velocity and the divergence that decays to zero with second order and the efficiency of the multigrid that is comparable with the tests available in literature for rectangular domains without objects showing that the presence of a complex shaped object does not degrade the performance

a computational stream function method for two dimensional incompressible viscous flows this work concerns the development of a numerical method based on the stream function formulation of the navier stokes equations to simulate two dimensional plane or axisymmetric viscous flows the main features of the proposed method are the use of the high order finite difference compact method for the discretization of the stream function equation the implicit pseudo transient newton krylov multigrid matrix free method for the stationary stream function equation and the fourth order runge kutta method for the integration of non stationary flows copyright 2005 john wiley sons ltd

nodal intersections for random eigenfunctions on the torus we investigate the number of nodal intersections of random gaussian laplace eigenfunctions on the standard two dimensional flat torus arithmetic random waves with a fixed smooth reference curve with nonvanishing curvature the expected intersection number is universally proportional to the length of the reference curve times the wavenumber independent of the geometry our main result prescribes the asymptotic behavior of the nodal intersections variance for smooth curves in the high energy limit remarkably it is dependent on both the angular distribution of lattice points lying on the circle with radius corresponding to the given wavenumber and the geometry of the given curve in particular this implies that the nodal intersection number admits a universal asymptotic law with arbitrarily high probability

computational design for long term numerical integration of the equations of fluid motion the integral constraints on quadratic quantities of physical importance such as conservation of mean kinetic energy and mean square vorticity will not be maintained in finite difference analogues of the equation of motion for two dimensional incompressible flow unless the finite difference jacobian expression for the advection term is restricted to a form which properly represents the interaction between grid points as derived in this paper it is shown that the derived form of the finite difference jacobian prevents nonlinear computational instability and thereby permits long term numerical integrations

implicit solutions of incompressible navier stokes equations using the pressure gradient method a modified pressure gradient method is developed for solving the incompressible twoand three dimensional navier stokes equations in primitive variables velocity and pressure gradient vectors rather than the velocity and pressure are considered as dependent variables this choice of the dependent variables is consistent with the flow physics as the velocity field of incompressible flows is not dependent on the pressure but rather the pressure gradient this introduces one and two additional dependent variables for two and three dimensions respectively to close the system of equations additional equations are obtained from the identity that the curl of the pressure gradient is zero this choice of dependent variables also leads to dirichlet boundary conditions for the pressure gradient which improves numerical stability and convergence rates the modified pressure gradient method includes a pressure gradient averaging term in the momentum equations to stabilize the solution this term is extended here to general curvilinear coordinates the system of equations is discritized using a conservative finite volume formulation time integration is achieved with a fully implicit scheme including implicit boundary conditions the implicit scheme guarantees that continuity is satisfied to machine zero numerical results compare well with solutions available in the literature for the cavity flow problem and flow about a circular cylinder at reynolds number 40

a non staggered grid fractional step method for time dependent incompressible navier stokes equations in curvilinear coordinates a numerical method for solving three dimensional time dependent incompressible navier stokes equations in curvilinear coordinates is presented the non staggered grid method originally developed by c m rhie and w l chow aiaaj 21 1525 1983 for steady state problems is extended to compute unsteady flows in the computational space the cartesian velocity components and the pressure are defined at the center of a control volume while the volume fluxes are defined at the mid point on their corresponding cell faces the momentum equations are integrated semi implicitly by the approximate factorization technique the intermediate velocities are interpolated onto the faces of the control volume to form the source terms of the pressure poisson equation which is solved iteratively with a multigrid method the compatibility condition of the pressure poisson equation is satisfied in the same manner as in a staggered grid method mass conservation can be satisfied to machine accuracy the pressure boundary condition is derived from the momentum equations solutions of both steady and unsteady problems including the large eddy simulation of a rotating and stratified upwelling flow in an irregular container established the favorable accuracy and efficiency of the present method

short note an interface capturing method with a continuous function the thinc method on unstructured triangular and tetrahedral meshes a novel interface capturing method is proposed to compute moving interfaces on unstructured grids with triangular 2d and tetrahedral 3d elements different from the conventional vof volume of fluid method which involves geometric reconstructions of the interface the present method is based on the algebraic reconstruction approach originally developed in the thinc tangent of hyperbola interface capturing scheme by xiao et al 2005 17 a continuous multidimensional hyperbolic tangent function is employed for retrieving the jump like distribution of the indicator function which avoids the explicit geometric representation of the interface and thus substantially reduces the algorithmic complexity in unstructured grids numerical diffusion and smearing are effectively eliminated and the compact thickness of the jump transition layer in the volume fraction is retained throughout the computation even for largely deformed interface the solution quality of the present scheme is comparable to the vof method with plic piecewise linear interface calculation algorithm

a geometrical area preserving volume of fluid advection method a new class of algorithms that preserve mass exactly for incompressible flows on a cartesian mesh are presented they amount to piecewise linear area preserving mappings of tessellations of the plane they are equivalent to volume of fluid vof advection methods which are decomposed into an eulerian implicit scheme in one direction followed by a lagrangian explicit step in the other one it is demonstrated that mass conservation is exact for incompressible flows and that there are no undershoots or overshoots of the volume fraction which thus always remains constrained between 0 and 1

a cartesian grid projection method for the incompressible euler equations in complex geometries many problems in fluid dynamics require the representation of complicated internal or external boundaries of the flow here we present a method for calculating time dependent incompressible inviscid flow which combines a projection method with a cartesian grid approach for representing geometry in this approach the body is represented as an interface embedded in a regular cartesian mesh the advection step is based on a cartesian grid algorithm for compressible flow in which the discretization of the body near the flow uses a volume of fluid representation a redistribution procedure is used to eliminate time step restrictions due to small cells where the boundary intersects the mesh the projection step uses an approximate projection based on a cartesian grid method for potential flow the method incorporates knowledge of the body through volume and area fractions along with certain other integrals over the mixed cells convergence results are given for the projection itself and for the time dependent algorithm in two dimensions the method is also demonstrated on flow past a half cylinder with vortex shedding

a novel dynamic framework for subgrid scale parametrization of mesoscale eddies in quasigeostrophic turbulent flows abstract this paper puts forth a modular dynamic subgrid scale modeling framework for large eddy simulation of quasigeostrophic turbulence based on minimizing the errors between structural and functional subgrid scale models the approximate deconvolution ad procedure is used to estimate the free modeling parameters for the eddy viscosity coefficient parameterized in space and time using the smagorinsky and leith models the novel idea here is to estimate the modeling parameters using the ad method rather than the traditionally used test filtering approach first a priori and a posteriori analyses are presented for solving a canonical homogeneous isotropic decaying turbulence problem the proposed model is then applied to a wind driven quasigeostrophic four gyre ocean circulation problem which is a standard prototype of more realistic ocean dynamics results show that the proposed model captures the quasi stationary ocean dynamics and provides the time averaged four gyre circulation patterns taking locally resolved flow characteristics into account the model dynamically provides higher eddy viscosity values for lower resolutions

accurate stable and efficient navier stokes solvers based on explicit treatment of the pressure term we present numerical schemes for the incompressible navier stokes equations based on a primitive variable formulation in which the incompressibility constraint has been replaced by a pressure poisson equation the pressure is treated explicitly in time completely decoupling the computation of the momentum and kinematic equations the result is a class of extremely efficient navier stokes solvers full time accuracy is achieved for all flow variables the key to the schemes is a neumann boundary condition for the pressure poisson equation which enforces the incompressibility condition for the velocity field irrespective of explicit or implicit time discretization of the viscous term in the momentum equation the explicit time discretization of the pressure term does not affect the time step constraint indeed we prove unconditional stability of the new formulation for the stokes equation with explicit treatment of the pressure term and first or second order implicit treatment of the viscous term systematic numerical experiments for the full navier stokes equations indicate that a second order implicit time discretization of the viscous term with the pressure and convective terms treated explicitly is stable under the standard cfl condition additionally various numerical examples are presented including both implicit and explicit time discretizations using spectral and finite difference spatial discretizations demonstrating the accuracy flexibility and efficiency of this class of schemes in particular a galerkin formulation is presented requiring only c0 elements to implement

numerical investigation of the three dimensional development in boundary layer transition a numerical method for solving the complete navier stokes equations for incompressible flows is introduced that is applicable for investigating three dimensional transition phenomena in a spatially growing boundary layer results are discussed for a test case with small three dimensional disturbances for which detailed comparison to linear stability theory is possible the validity of our numerical model for investigating nonlinear transition phenomena is demonstrated by realistic spatial simulations of the experiments by kachanov and levchenko1 for a subharmonic resonance breakdown and of the experiments of klebanoff et al 2 for a fundamental resonance breakdown

a cell centered adaptive projection method for the incompressible euler equations we present an algorithm to compute adaptive solutions for incompressible flows using block structured local refinement in both space and time this method uses a projection formulation based on a cell centered approximate projection which allows the use of a single set of cell centered solvers because of refinement in time additional steps are taken to accurately discretize the advection and projection operators at grid refinement boundaries using composite operators which span the coarse and refined grids this ensures that the method is approximately freestream preserving and satisfies an appropriate form of the divergence constraint

arithmetic zeros and nodal domains on the sphere we obtain lower bounds for the number of nodal domains of hecke eigenfunctions on the sphere assuming the generalized lindelof hypothesis we prove that the number of nodal domains of any hecke eigenfunction grows with the eigenvalue of the laplacian by a very different method we show unconditionally that the average number of nodal domains of degree l hecke eigenfunctions grows significantly faster than the uniform growth obtained under lindelof

an overview of projection methods for incompressible flows a series of numerical issues related to the analysis and implementation of fractional step methods for incompressible flows are addressed in this paper these methods are often referred to in the literature as projection methods and can be classified into three classes namely the pressure correction methods the velocity correction methods and the consistent splitting methods for each class of schemes theoretical and numerical convergence results available in the literature are reviewed and open questions are discussed the essential results are summarized in a table which could serve as a useful reference to numerical analysts and practitioners published by elsevier b v msc 65n35 65n22 65f05 35j05

a time integration method using artificial compressibility for unsteady viscous flows abstract a time accurate integration method for the study of unsteady incompressible viscous flows is presented based on a three point backward implicit discretization of the navier stokes equations in real time and a pseudo time relaxation with artificial compressibility to advance the solution between consecutive real time levels the spatial discretization uses central differencing based on a stretched staggered grid a factored adi scheme with lagged non linearities is used to reduce the problem to the solution of scalar tridiagonal systems this method displays very good computing efficiency and accuracy when applied for validation to several typical steady and unsteady flow problems

gerris a tree based adaptive solver for the incompressible euler equations in complex geometries an adaptive mesh projection method for the time dependent incompressible euler equations is presented the domain is spatially discretised using quad octrees and a multilevel poisson solver is used to obtain the pressure complex solid boundaries are represented using a volume of fluid approach second order convergence in space and time is demonstrated on regular statically and dynamically refined grids the quad octree discretisation proves to be very flexible and allows accurate and efficient tracking of flow features the source code of the method implementation is freely available

improving the speed and accuracy of projection type incompressible flow solvers abstract superseding so called first generation incompressible flow solvers of the projection type based on taylor galerkin advection second order pressure damping and element based data structures the current second generation solvers based on high order upwind advection fourth order pressure damping and edge based data structures have now been in use for half a decade and have proven remarkably robust and efficient for many large scale problems in order to achieve higher accuracy and speed these solvers have recently been enhanced in a variety of ways a substepping for advection b implicit treatment of advective terms via sgs and gmres lu sgs iterative solvers c fully implicit time accurate advancement of pressure and velocities and d linelet preconditioning for the pressure poisson equation the combined effect of these third generation improvements leads to speedups of the order of o 1 5 1 10 with similar or even better temporal accuracy as demonstrated on a variety of academic and industrial problems

a new algorithm for solving some mechanical problems this paper explores the utility of a discrete singular convolution algorithm for solving certain mechanical problems benchmark mechanical systems including plate vibrations and incompressible flows are employed to illustrate the robustness and to test accuracy of the present algorithm numerical results indicate that the present approach is very accurate efficient and reliable for solving the aforementioned problems

how many zeros of a random polynomial are real we provide an elementary geometric derivation of the kac integral formula for the expected number of real zeros of a random polynomial with independent standard normally distributed coefficients we show that the expected number of real zeros is simply the length of the moment curve 1 t tn projected onto the surface of the unit sphere divided by the probability density of the real zeros is proportional to how fast this curve is traced out we then relax kac s assumptions by considering a variety of random sums series and distributions and we also illustrate such ideas as integral geometry and the fubini study metric

analysis of a one dimensional model for the immersed boundary method numerical methods are studied for the one dimensional heat equation with a singular forcing term u t u xx c t delta x alpha t the delta function delta x is replaced by a discrete approximation d h x and the resulting equation is solved by a crank nicolson method on a uniform grid the accuracy of this method is analyzed for various choices of d h the case where c t is specified and also the case where c is determined implicitly by a constraint on the solution at the point a are studied these problems serve as a model for the immersed boundary method of peskin for incompressible flow problems in irregular regions some insight is gained into the accuracy that can be achieved and the importance of choosing appropriate discrete delta functions

an improved particle correction procedure for the particle level set method the particle level set method d enright r fedkiw j ferziger i mitchell a hybrid particle level set method for improved interface capturing j comput phys 183 2002 83 116 can substantially improve the mass conservation property of the level set method by using lagrangian marker particles to correct the level set function in the under resolved regions in this study the limitations of the particle level set method due to the errors introduced in the particle correction process are analyzed and an improved particle correction procedure is developed based on a new interface reconstruction scheme moreover the zero level set is anchored as the level set functions are reinitialized hence the additional particle correction after the level set reinitialization is avoided with this new scheme a well defined zero level set can be obtained and the disturbances to the interface are significantly reduced consequently the particle reseeding operation will barely result in the loss of interface characteristics and can be applied as frequently as necessary to demonstrate the accuracy and robustness of the proposed method two extreme particle reseeding strategies one without reseeding and the other with reseeding every time step are applied in several benchmark advection tests and the results are compared with each other three interfacial flow cases a 2d surface tension driven oscillating droplet a 2d gas bubble rising in a quiescent liquid and a 3d drop impact onto a liquid pool are simulated to illustrate the advantages of the current method over the level set and the original particle level set methods with regard to the smoothness of geometric properties and mass conservation in real physical applications

uncoupled temporally second order accurate implicit solver of incompressible navier stokes equations an efficient method with temporal second order accuracy for solving the incompressible navier stokes equations in generalized coordinate systems is developed efficiency is obtained by extending an existing fractional step semi implicit solver into an implicit one with implicit approximation of the convection terms a novel three level linearization scheme is proposed to decouple the momentum equations from one another in addition to the fractional step approach that decouples the continuity equation from the momentum equations consequently all of the governing equations are uncoupled without degrading temporal accuracy or stability and without the need for iterations to solve the nonlinear system of equations at each time level the proposed decoupling technique can be used for other systems of nonlinear and time dependent partial differential equations moreover only minor modifications are required to apply the method to existing two level schemes several test cases confirm that the proposed decoupled scheme has temporal second order accuracy and allows the use of large time steps

a well conditioned augmented system for solving navier stokes equations in irregular domains an augmented method based on a cartesian grid is proposed for the incompressible navier stokes equations in irregular domains the irregular domain is embedded into a rectangular one so that a fast poisson solver can be utilized in the projection method unlike several methods suggested in the literature that set the force strengths as unknowns which often results in an ill conditioned linear system we set the jump in the normal derivative of the velocity as the augmented variable the new approach improves the condition number of the system for the augmented variable significantly using the immersed interface method we are able to achieve second order accuracy for the velocity numerical results and comparisons to benchmark tests are given to validate the new method a lid driven cavity flow with multiple obstacles and different geometries are also presented

a conservative adaptive projection method for the variable density incompressible navier stokes equations in this paper we present a method for solving the equations governing time dependent variable density incompressible flow in two or three dimensions on an adaptive hierarchy of grids the method is based on a projection formulation in which we first solve advection diffusion equations to predict intermediate velocities and then project these velocities onto a space of approximately divergence free vector fields our treatment of the first step uses a specialized second order upwind method for differencing the nonlinear convection terms that provides a robust treatment of these terms suitable for inviscid and high reynolds number flow density and other scalars are advected in such a way as to maintain conservation if appropriate and free stream preservation our approach to adaptive refinement uses a nested hierarchy of logically rectangular girds with simultaneous refinement of the girds in both space and time the integration algorithm on the grid hierarchy is a recursive procedure in which coarse grids are advanced in time fine grids are advanced multiple steps to reach the same time as the coarse grids and the data at different levels are then synchronized the single grid algorithm is described briefly but the emphasis here is on the time stepping procedure for the adaptive hierarchy numerical examples are presented to demonstrate the algorithms s accuracy and convergence properties and illustrate the behavior of the method an additional example demonstrates the performance of the method on a more realistic problem namely a three dimensional variable density shear layer

extreme values of zeta and l functions we introduce a resonance method to produce large values of the riemann zeta function on the critical line and large and small central values of l functions

a coupled level set and volume of fluid method for computing 3d and axisymmetric incompressible two phase flows we present a coupled level set volume of fluid clsvof method for computing 3d and axisymmetric incompressible two phase flows this method combines some of the advantages of the volume of fluid method with the level set method to obtain a method which is generally superior to either method alone we present direct comparisons between computations made with the clsvof method and computations made with the level set method the volume of fluid method and the boundary integral method we also compare our computations to the exact solution for an oscillating ellipse due to lamb and experimental results obtained for a rising gas bubble in liquid obtained by hnat and buckmaster our computational examples focus on flows in which surface tension forces and changes in topology are dominant features of the flow

a diffuse interface model for axisymmetric immiscible two phase flow a diffuse interface model is considered for solving axisymmetric immiscible two phase flow with surface tension the navier stokes ns equations are modified by the addition of a continuum forcing the interface between the two fluids is considered as the half level set of a mass concentration c which is governed by the cahn hilliard ch equation a fourth order degenerate nonlinear parabolic diffusion equation in this work we develop a nonlinear multigrid method to solve the ch equation with degenerate mobility and couple this to a projection method for the incompressible ns equations the diffuse interface method can deal with topological transitions such as breakup and coalescence smoothly without ad hoc cut and connect or other artificial procedures we present results for rayleigh s capillary instability up to forming satellite drops the results agree well with the linear stability theory

a hybrid particle level set method for improved interface capturing in this paper we propose a new numerical method for improving the mass conservation properties of the level set method when the interface is passively advected in a flow field our method uses lagrangian marker particles to rebuild the level set in regions which are underresolved this is often the case for flows undergoing stretching and tearing the overall method maintains a smooth geometrical description of the interface and the implementation simplicity characteristic of the level set method our method compares favorably with volume of fluid methods in the conservation of mass and purely lagrangian schemes for interface resolution the method is presented in three spatial dimensions