rational curves and rational singularities we study rational curves on algebraic varieties especially on normal affine varieties endowed with a action for varieties with an isolated singularity covered by a family of rational curves with a general member not passing through the singular point we show that this singularity is rational in particular this provides an explanation of classical results due to h a schwartz and g h halphen on polynomial solutions of the generalized fermat equation

a semi implicit numerical scheme for reacting flow an additive semi implicit projection scheme for the simulation of unsteady combustion in two dimensions is constructed the scheme relies on a zero mach number formulation of the compressible conservation equations with detailed chemistry the governing equations are discretized in space using second order differences and integrated in time using a semi implicit approach time integration of the evolution equations for species mass fraction thermodynamic pressure and density is performed using a semi implicit nonsplit scheme that combines a second order predictor corrector treatment of convection and diffusion terms and a stiff integrator for the reaction source terms meanwhile the momentum equations are integrated using a second order projection scheme the projection scheme is based on a predictor corrector approach that couples the evolution of the velocity and density fields in order to stabilize computations of reacting flows with large density variations a pressure poisson equation is inverted following both the predictor and corrector steps using a fast solver the advantages of the stiff integration of reaction source terms are analyzed by comparing the performance of the scheme to that of a predictor corrector scheme in which reaction and diffusion are integrated in a similar nonstiff fashion the comparison in based on both one dimensional 1d unsteady tests of a premixed methane air flame and unsteady two dimensional tests of the same flame interacting with a counterrotating vortex pair in both cases the grimech1 2 reaction mechanism with 32 species and 177 elementary reactions is used computed results show that the stiff reaction scheme enables selection of larger time steps and thus leads to substantial improvement in the performance of the computations for the present reaction mechanism and flame conditions speedup factors of about 10 are achieved in the 1d tests and about five in two dimensions possible extensions of the present scheme to further improve efficiency are also discussed

a new finite variable difference method with application to locally exact numerical scheme a new finite variable difference method fvdm on finite differencing is presented the essence of this method consists in determining the optimum spatial difference such that the total variance of the solution is a minimum under the condition that characteristic roots of the resulting difference equation are always nonnegative to ensure the numerical stability the present fvdm is applied to the locally exact numerical scheme lens the optimum spatial difference of the lens is derived in terms of local mesh reynolds numbers by using this optimum spatial difference the numerical accuracy of the lens for the linear convection diffusion equations is increased without numerical oscillations for all mesh reynolds numbers the present study suggests that an optimum spatial difference from the viewpoint of numerical stability and accuracy exists according to the numerical schemes

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

examples of non algebraic classes in the brown peterson tower for every n ge 0 we construct classes in the brown peterson cohomology bp langle n rangle of smooth projective complex algebraic varieties which are not in the image of the cycle map from the corresponding motivic brown peterson cohomology this generalizes the examples of atiyah and hirzebruch to all finite levels in the brown peterson tower

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

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

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

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

a fast computation technique for the direct numerical simulation of rigid particulate flows in this paper we present a computation technique for the direct numerical simulation of freely moving rigid bodies in fluids we solve three dimensional laminar flow problems using a control volume approach the key feature of this approach is that the computational overhead relative to a pure fluid solver to solve for the motion of rigid particle is very small the formulation is convenient for handling irregular geometries we present results for the sedimentation of particles of different shapes convergence tests are presented to assess the order of accuracy of the numerical scheme various test cases are considered and the numerical results are compared with experimental values to validate the code due to the ability to perform fast computations this method has been used for animations and its application to the direct numerical simulation of turbulent particulate flows merits investigation the technique is not restricted to any constitutive model of the suspending fluid hence it may potentially be used in large eddy simulations les or reynolds averaged navier stokes rans type simulations

meshless schemes for unsteady navier stokes equations in vorticity formulation using radial basis functions in this paper we provide a new scheme for unsteady incompressible flows in vorticity stream function formulation combined with the radial basis functions method it is an efficient meshless method optimal accuracy can be achieved using this method the efficiency and accuracy are demonstrated by numerical examples

a divergence free condition compensated method for incompressible navier stokes equations the present study aims to develop a new formulation to effectively calculate the incompressible navier stokes solutions in non staggered grids the distinguished feature of the proposed method which avoids directly solving the divergence free equation is to add a rigorously derived source term to the momentum equation to ensure satisfaction of the fluid incompressibility for the sake of numerical accuracy dispersion relation preserving upwind scheme developed within the two dimensional context was employed to approximate the convection terms the validity of the proposed mass preserving navier stokes method is justified by solving two benchmark problems at high reynolds and rayleigh numbers based on the simulated navier stokes solutions the proposed formulation is shown to outperform the conventional segregated method in terms of the reduction of cpu time 2007 elsevier b v all rights reserved

torsion cohomology classes and algebraic cycles on complex projective manifolds abstract atiyah and hirzebruch gave examples of even degree torsion classes in the singular cohomology of a smooth complex projective manifold which are not poincare dual to an algebraic cycle we notice that the order of these classes must be small compared to the dimension of the manifold however building upon a construction of kollar one can provide such examples with arbitrary high prime order the dimension being fixed this method also provides examples of torsion algebraic cycles which are non trivial in the griffiths group and lie in a arbitrary high level of the h saito filtration on chow groups

higher entropy conservation and numerical stability of compressible turbulence simulations we present a numerical formulation for the treatment of nonlinear instabilities in shock free compressible turbulence simulations the formulation is high order and contains no artificial dissipation numerical stability is enhanced through semi discrete satisfaction of global conservation properties stemming from the second law of thermodynamics and the entropy equation the numerical implementation is achieved using a conservative skew symmetric splitting of the nonlinear terms the robustness of the method is demonstrated by performing unresolved numerical simulations and large eddy simulations of compressible isotropic turbulence at a very high reynolds number results show the scheme is capable of capturing the statistical equilibrium of low mach number compressible turbulent fluctuations at infinite reynolds number comparisons with the entropy splitting technique j comput phys 162 2000 33 j comput phys 178 2002 307 staggered method j comput phys 191 2 2003 392 and skew symmetric like schemes j comput phys 161 2000 114 confirm the superiority of the current approach we also discuss a flaw in the skew symmetric splitting implemented in the literature very good results are obtained based on the proper splitting

high re solutions for incompressible flow using the navier stokes equations and a multigrid method the vorticity stream function formulation of the two dimensional incompressible navierstokes equations is used to study the effectiveness of the coupled strongly implicit multigrid csi mg method in the determination of high re fine mesh flow solutions the driven flow in a square cavity is used as the model problem solutions are obtained for configurations with reynolds number as high as 10 000 and meshes consisting of as many as 257 x 257 points for re 1000 the 129 x 129 grid solution required 1 5 minutes of cpu time on the amdahl 470 v 6 computer because of the appearance of one or more secondary vortices in the flow field uniform mesh refinement was preferred to the use of one dimensional gridclustering coordinate transformations the past decade has witnessed a great deal of progress in the area of computational fluid dynamics developments in computer technology hardware as well as in advanced numerical algorithms have enabled attempts to be made towards analysis and numerical solution of highly complex flow problems for some of these applications the use of simple iterative techniques to solve the navier stokes equations leads to a rather slow convergence rate for the solutions the solution convergence rate can be seriously affected if the coupling among the various governing differential equations is not properly honored either in the interior of the solution domain or at its boundaries the rate of convergence is also generally strongly dependent on such problem parameters as the reynolds number the mesh size and the total number of computational points this has led several researchers to examine carefully the recently emerging multigrid mg technique as a useful means for enhancing the convergence rate of iterative numerical methods for solving discretized equations at a number of computational grid points so large as to be considered impractical previously

driven cavity flows by efficient numerical techniques abstract efficient and reliable numerical techniques of high order accuracy are presented for solving problems of steady viscous incompressible flow in the plane and are used to obtain accurate solutions for the driven cavity a solution is obtained at reynolds number 10 000 on a 180 180 grid the numerical methods combine an efficient linear system solver an adaptive newton like method for nonlinear systems and a continuation procedure for following a branch of solutions over a range of reynolds numbers

torsion algebraic cycles and etale cobordism we prove that the classical integral cycle class map from algebraic cycles to etale cohomology factors through a quotient of l adic etale cobordism over an algebraically closed field of positive characteristic this shows that there is a strong topological obstruction for cohomology classes to be algebraic and that examples of atiyah hirzebruch and totaro also work in positive characteristic

a monotonicity preserving conservative sharp interface flow solver for high density ratio two phase flows this paper presents a novel approach for solving the conservative form of the incompressible two phase navier stokes equations in order to overcome the numerical instability induced by the potentially large density ratio encountered across the interface the proposed method includes a volume of fluid type integration of the convective momentum transport a monotonicity preserving momentum rescaling and a consistent and conservative ghost fluid projection that includes surface tension effects the numerical dissipation inherent in the volume of fluid treatment of the convective transport is localized in the interface vicinity enabling the use of a kinetic energy conserving discretization away from the singularity two and three dimensional tests are presented and the solutions shown to remain accurate at arbitrary density ratios the proposed method is then successfully used to perform the detailed simulation of a round water jet emerging in quiescent air therefore suggesting the applicability of the proposed algorithm to the computation of realistic turbulent atomization

etale realization on the a1 homotopy theory of schemes abstract we compare friedlander s definition of the etale topological type for simplicial schemes to another definition involving realizations of pro simplicial sets this can be expressed as a notion of hypercover descent for etale homotopy we use this result to construct a homotopy invariant functor from the category of simplicial presheaves on the etale site of schemes over s to the category of pro spaces after completing away from the characteristics of the residue fields of s we get a functor from the morel voevodsky a 1 homotopy category of schemes to the homotopy category of pro spaces

a new numerical approach to solve an elliptic equation a new numerical approach to solve an elliptic partial differential equation that originates from the governing equations of steady state fluid flow and heat transfer is presented the elliptic partial differential equation is transformed by introducing an exponential function to eliminate the convection terms in the equation a fourth order central differencing scheme and a second order central differencing scheme are used to numerically solve the transformed elliptic partial differential equation analytical solutions of this equation are also given comparisons are made between the analytical solutions the numerical results using the present schemes and those using the four classical differencing schemes namely the first order upwind scheme hybrid scheme power law scheme and exponential scheme the comparisons illustrate that the proposed algorithm performs better than the four classical differencing schemes

a weighted meshfree collocation method for incompressible flows using radial basis functions abstract a weighted strong form collocation method using radial basis functions and explicit time integration is proposed to solve the incompressible navier stokes equations the velocities and pressure are solved directly at the same time step and the continuity equation is satisfied at each time step which improve the solution accuracy and stability no artificial compressibility coefficient needs to be introduced for modeling the incompressible flows and no pressure oscillation arises in the numerical solutions optimal convergence can be achieved by imposing the derived proper weights on the boundaries and the continuity equation radial basis collocation method in a lagrangian form is quite easy to capture the moving boundary or free surface in flow problems moreover solid boundary conditions can be enforced directly and no special treatments are required further critical time step for the explicit time integration is predicted in the stability analysis and the influences on the stability are evaluated numerical studies validate the high accuracy as well as good stability of the presented method

numerical simulation of cavity flows based on transformed equations a new numerical algorithm is applied to simulate two dimensional lid driven cavity flows in this new algorithm the momentum equations are first transformed using an exponential function to eliminate the convection terms in the equations then a central differencing scheme is employed to discretize the transformed equations the cavity flows studied in this work include those with non zero velocity component in the y direction on the upper and lower boundaries the results for the velocity components along the geometric centerline stream function patterns and vorticity contours are presented and discussed the predicted results are in excellent agreement with benchmark solutions

closed form numerical solutions of variable coefficient linear second order elliptic problems abstract in this work we develop an alternative numerical technique which allows to construct a numerical solution in closed form of variable coefficient linear second order elliptic problems with dirichlet boundary conditions the elliptic partial differential equation is approximated by a consistent explicit difference scheme and using a discrete separation of the variables method we determine a closed form solution of the two resulting discrete boundary value problems with the separated variables avoiding to have to solve large algebraic systems one of these boundary value problems is a discrete sturm liouville problem which guarantees the qualitative properties of the exact solution of elliptic problem a constructive procedure for the computation of the numerical solution is given and an illustrative example is included

numerical solution of a transformed parabolic equation a numerical analysis of a parabolic partial differential equation pde which originates from the governing equations of transient fluid flow and heat transfer is presented the parabolic pde is transformed by introducing an exponential function to eliminate the convection terms in the equation a fourth order central differencing scheme and a second order central differencing scheme are used to solve the transformed parabolic pde numerically the analytical solutions of this equation are also given comparisons against the analytical solutions are made for the numerical results using the present schemes and those using the four classical differencing schemes namely the first order upwind scheme hybrid scheme power law scheme and exponential scheme

immersed boundary methods for general finite difference and finite volume navier stokes solvers we present an immersed boundary algorithm for incompressible flows with complex boundaries suitable for cartesian or curvilinear grid system the key stages of any immersed boundary technique are the interpolation of a velocity field given on a mesh onto a general boundary a line in 2d a surface in 3d and the spreading of a force field from the immersed boundary to the neighboring mesh points to enforce the desired boundary conditions on the immersed boundary points we propose a technique that uses the reproducing kernel particle method w k liu s jun y f zhang reproducing kernel particle methods int j numer methods fluids 20 8 1995 1081 1106 for the interpolation and spreading unlike other methods presented in the literature the one proposed here has the property that the integrals of the force field and of its moment on the grid are conserved independent of the grid topology uniform or non uniform cartesian or curvilinear the technique is easy to implement and is able to maintain the order of the original underlying spatial discretization applications to two and three dimensional flows in cartesian and non cartesian grid system with uniform and non uniform meshes are presented

stable etale realization and etale cobordism abstract we show that there is a stable homotopy theory of profinite spaces and use it for two main applications on the one hand we construct an etale topological realization of the stable a 1 homotopy theory of smooth schemes over a base field of arbitrary characteristic in analogy to the complex realization functor for fields of characteristic zero on the other hand we get a natural setting for etale cohomology theories in particular we define and discuss an etale topological cobordism theory for schemes it is equipped with an atiyah hirzebruch spectral sequence starting from etale cohomology finally we construct maps from algebraic to etale cobordism and discuss algebraic cobordism with finite coefficients over an algebraically closed field after inverting a bott element

linear stability of incompressible flow using a mixed finite element method we analyzed the linear stability for incompressible flow using the mixed finite element method first by means of the fractional step finite element method the two dimensional base flows were numerically computed over a range of reynolds numbers and then they were perturbed with three dimensional disturbances by using linear stability and normal mode analysis we obtained the partial differential equations governing the evolution of perturbation from linearized navier stokes equations with slight compressibility in terms of the mixed finite element discretization in which six node quadratic lagrange triangular elements with quadratic interpolation for velocities p2 and three node linear lagrange triangular elements for pressure p1 were used a nonsingular generalized eigenproblem was formulated from these equations whose solution gives the dispersion relationship between the complex growth rate and the wave number we examined the stability of lid driven cavity flow in this study the stability results obtained from thep2p1elements were compared with those fromp1p1elements the spatial distribution of the pressure disturbance can be smoothed by using ap2p1element such as that in the simulation of incompressible flows by the mixed finite element method however the choice of elemental type has little influence on the stability results this study presents stability curves to identify the critical reynolds number and the critical wavelength of the neutral mode and discusses the instability mechanism the taylor goertler like vortices in the cavity were obtained by means of the reconstruction of three dimensional flows the stability results and the reconstructed three dimensional flows were in agreement with the observed ones and the other numerical stability results respectively

b spline method and zonal grids for simulations of complex turbulent flows a numerical technique for computations of turbulent flows is described the technique is based on b splines and allows grid embedding in physically significant flow regions numerical tests which include solutions of nonlinear advection diffusion equations and computations of flow over a circular cylinder at reynolds numbers up to 300 indicate that the method is accurate and efficient in computations of flow over a cylinder the lift drag and base suction coefficients agree well with existing experimental data and previous numerical simulations

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

etale endomorphisms of algebraic surfaces with g_m actions let x be an algebraic surface defined over the complex field vec c and endowed with an vec a 1 fibration as such a surface we have a platonic vec a 1 fiber space a weighted hypersurface with its singular point deleted off and more generally an affine algebraic surface with an unmixed g m action and its fixpoint deleted off we consider an etale endomorphism varphi x to x and show that varphi is an automorphism in most cases of particular interest is the case of a platonic vec a 1 fiber space for which varphi being an automorphism is closely related to the jacobian problem for the affine plane vec a 2 we also investigate the automorphism group of such surfaces

a pseudospectral method for solution of the three dimensional incompressible navier stokes equations abstract a new chebyshev pseudospectral technique based on the projection method that was previously applied by the authors to the solution of two dimensional incompressible navier stokes equations in primitive variables for nonperiodic boundary conditions is extended to solve the three dimensional navier stokes equations the crucial point of the method is the requirement that the continuity equation be satisfied everywhere in the domain on the boundaries as well as in the interior the key feature of the work presented in this paper is that the computer storage requirements of the full matrix inversion resulting from direct solution of the pressure poisson equation in three dimensions is greatly reduced by considering an eigenfunction decomposition the method was tested on a two dimensional driven cavity flow and the results were compared with those of the most accurate finite difference calculation the three dimensional driven cavity flow was then calculated at the same reynolds numbers as the two dimensional cases i e re 100 400 and 1000 in the calculated results three dimensional boundary effects were observed in all cases and became more apparent with increasing reynolds number

a geometric approach to the jacobian conjecture in dimension two the jacobian conjecture in dimension two holds true if a conjecture which we call below the generalized sard property gsp holds true for the affine plane and an a1 fibration on a2 the observations will be made for affine pseudo planes with the jacobian conjecture enlarged to the generalized jacobian conjecture

a numerical scheme for transport equations with spatially distributed coefficients based on locally exact difference method a numerical difference scheme to solve stationary transport equations with spatially distributed coefficients is presented the spatial distribution of the coefficients in the transport equations is taken into consideration based on a four region model among three adjacent control volumes in which continuous conditions for solutions are imposed on the boundary between two adjacent regions the coefficients in the difference scheme are determined so that it will be satisfied exactly by any local solution of the continuous equations with piecewise constant coefficients in each region the present scheme is examined through numerical experiments for one dimensional convection diffusion equations with spatially distributed coefficients and source term and a two dimensional cavity flow problem the present scheme shows good solutions

a novel concept for the design of immersed interface methods the objective of this paper is to present a novel robust high order accurate immersed interface method iim for advection diffusion type equations in contrast to other immersed methods that were designed for consistency and accuracy with a posteriori check of the numerical stability we combine local taylor series expansion at irregular grid points with a local stability constraint as part of the design process stability investigations of the iim are employed to demonstrate that the local stability constraint is sufficient for obtaining a globally stable method as long as the neumann number is less than its limiting value one of the key aspects of this iim is that the irregular finite difference stencils can be isolated from the rest of the computational domain to validate our novel immersed interface approach two dimensional and three dimensional test cases for model equations are presented in addition this method is applied to the incompressible navier stokes equations to conduct stability investigations of a boundary layer flow over a rough surface and for investigations of pulsatile stenotic flows stability investigations of wall bounded flows are challenging for immersed methods because the near wall accuracy is important for correctly capturing the characteristics of the hydrodynamic instability mechanisms in particular regarding the wave relation between the wave velocity components close to the wall

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 stochastic projection method for fluid flow i basic formulation we describe the construction and implementation of a stochastic navier stokes solver the solver combines a spectral stochastic uncertainty representation scheme with a finite difference projection method for flow simulation the uncertainty quantification scheme is adapted from the spectral stochastic finite element method ssfem which is based on regarding uncertainty as generating a new dimension and the solution as being dependent on this dimension in the ssfem formalism the stochastic dependence is represented in terms of the polynomial chaos system and the coefficients in the corresponding spectral representation are obtained using a galerkin approach it is shown that incorporation of the spectral uncertainty representation scheme into the projection method results in a coupled system of advection diffusion equations for the various uncertainty fields and in a decoupled system of pressure projection steps this leads to a very efficient stochastic solver whose advantages are illustrated using steady and transient simulations of transport and mixing in a microchannel

an implicit forcing immersed boundary method for simulating viscous flows in irregular domains abstract we present a method for solving the incompressible navier stokes equations in irregular domains these equations are discretized using finite difference method in a uniform cartesian grid stationary rigid boundaries are embedded in the cartesian grid and singular forces are applied at the rigid boundaries to impose the no slip conditions the singular forces are then distributed to the nearby cartesian grid points using the immersed boundary method in the present work the singular forces are computed implicitly by solving a small system of equations at each time step this system of equations is derived from a second order projection method the main advantage of this method is that it imposes the no slip boundary condition exactly and avoids the need for small time step to maintain stability the ability of the method to simulate viscous flows in irregular domains is demonstrated by applying to 2 dimensional flows past a circular cylinder multiple rigid obstacles and 3 dimensional flow past a sphere

augmented lagrangian preconditioner for large scale hydrodynamic stability analysis abstract hydrodynamic linear stability analysis of large scale three dimensional configurations is usually performed with a time stepping approach based on the adaptation of existing solvers for the unsteady incompressible navier stokes equations we propose instead to solve the nonlinear steady equations with the newton method and to determine the largest growth rate eigenmodes of the linearized equations using a shift and invert spectral transformation and a krylov schur algorithm the solution of the shifted linearized navier stokes problem which is the bottleneck of this approach is computed via an iterative krylov subspace solver preconditioned by the modified augmented lagrangian mal preconditioner benzi et al 2011 the well known efficiency of this preconditioned iterative strategy for solving the real linearized steady state equations is assessed here for the complex shifted linearized equations the effect of various numerical and physical parameters is investigated numerically on a two dimensional flow configuration confirming the reduced number of iterations over state of the art steady state and time stepping based preconditioners a parallel implementation of the steady navier stokes and eigenvalue solvers developed in the freefem language suitably interfaced with the petsc slepc libraries is described and made openly available to tackle three dimensional flow configurations its application on a small scale three dimensional problem shows the good performance of this iterative approach over a direct l u factorization strategy in regards of memory and computational time on a large scale three dimensional problem with 75 million unknowns a 80 parallel efficiency on 256 up to 2048 processes is obtained

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 stable and accurate convective modelling procedure based on quadratic upstream interpolation a convective modelling procedure is presented which avoids the stability problems of central differencing while remaining free of the inaccuracies of numerical diffusion associated with upstream differencing for combined convection and diffusion the number of operations at each grid point is comparable to that of standard upstream pluscentral differencing however highly accurate solutions can be obtained with a grid spacing much larger than that required by conventional methods for comparable accuracy with obvious practical advantaged in terms of both speed and storage the algorithm is based on a conservative control volume formulation with cell wall values of each field variable written in terms of a quadratic interpolation using in any one coordinate direction the two adjacent nodal values together with the value at the next upstream node this results in a convective differencing scheme with greater formal accuracy than central differencing while retaining the basic stable convective sensitivity property of upstream weighted schemes the consistent treatment of diffusion terms is equivalent to central differencing with careful modelling numerical boundary conditions are not troublesome some idealized problems are studied showing the practical advantages of the method over other schemes in comparison with exact solutions an application to a complex unsteady two dimensional flow is briefly discussed

a fixed grid sharp interface method for bubble dynamics and phase change abstract a numerical method has been developed for direct simulation of bubble dynamics with large liquid to vapor density ratio and phase change the numerical techniques are based on a fixed grid finite volume method capable of treating the interface as a sharp discontinuity the unsteady axisymmetric navier stokes equations and energy equation in both liquid and vapor phases are computed the mass momentum and energy conditions are explicitly matched at the phase boundary to determine the interface shape and movement the cubic b spline is used in conjunction with a fairing algorithm to yield smooth and accurate information of curvatures nondimensional parameters including reynolds weber and jakob numbers are varied to offer insight into the physical and numerical characteristics of the bubble dynamics based on the present sharp interface approach bubble dynamics for density ratio of 1600 or higher with and without phase change can be successfully computed

torsion algebraic cycles and complex cobordism atiyah and hirzebruch gave the first counterexamples to the hodge conjecture with integral coefficients 3 that conjecture predicted that every integral cohomology class of hodge type p p on a smooth projective variety should be the class of an algebraic cycle but atiyah and hirzebruch found additional topological properties which must be satisfied by the integral cohomology class of an algebraic cycle here we provide a more systematic explanation for their results by showing that the classical cycle map from algebraic cycles modulo algebraic equivalence to integral cohomology factors naturally through a topologically defined ring which is richer than integral cohomology the new ring is based on complex cobordism a well developed topological theory which has been used only rarely in algebraic geometry since hirzebruch used it to prove the riemann roch theorem 17 this factorization of the classical cycle map implies the topological restrictions on algebraic cycles found by atiyah and hirzebruch it goes beyond their work by giving a topological method to show that the classical cycle map can be noninjective as well as nonsurjective the kernel of the classical cycle map is called the griffiths group and the topological proof here that the griffiths group can be nonzero is the first proof of this fact which does not use hodge theory the proof here gives nonzero torsion elements in the griffiths group whereas griffiths s hodge theoretic proof gives nontorsion elements 13 this topological argument also gives examples of algebraic cycles in the kernel of various related cycle maps where few or no examples were known before thus answering some questions posed by colliot thelene and schoen 8 p 14 37 p 13 colliot thelene asked in particular whether the map ch2 x n h4 x z n is injective for all smooth complex projective varieties x here ch x is the group of codimension i algebraic cycles modulo rational equivalence the first examples where colliot thelene s map is not injective were found by kollar and van geemen see 4 p 135 very recently bloch and esnault found examples defined over number fields 7 over nonalgebraically closed fields k there are other examples of smooth projective varieties xk with ch2 xk n hg4t xk z n not injective due to colliot thelene and sansuc as reinterpreted by salberger see 9 and 8 remark 7 6 1 and parimala and suresh 31 these elements of ch2 xk n are not shown to remain nonzero in ch2 xc n however here our topological method

a new numerical algorithm based on transformed equations and its applications to very low re fluid flows a new numerical algorithm is presented in this study to solve the partial differential equations that govern fluid flows this new algorithm is based on first transforming the partial differential equations by introducing an exponential function to eliminate the convection terms a fourth order central differencing scheme and a second order central differencing scheme are used to discretize the transformed equations the algorithm is then applied to simulate fluid flows with exact solutions to validate this new algorithm the fluid flows used in this study are a self designed quasi fluid flow problem stagnation in plane flow hiemenz flow and flow between two concentric cylinders comparisons against the exact solution are made for the results obtained using the new numerical algorithm as well as the power law scheme the comparisons indicate that the present fourth order scheme performs the best and the present second order scheme is the next most accurate