edge location to subpixel values in digital imagery a new method for locating edges in digital data to subpixel values and which is invariant to additive and multiplicative changes in the data is presented for one dimensional edge patterns an ideal edge is fit to the data by matching moments it is shown that the edge location is related to the so called christoffel numbers also presented is the study of the effect of additive noise on edge location the method is extended to include two dimensional edge patterns where a line equation is derived to locate an edge this in turn is compared with the standard hueckel edge operator an application of the new edge operator as an edge detector is also provided and is compared with sobel and hueckel edge detectors in presence and absence of noise

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 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 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

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

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

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

gauge transformation between the first order nonisospectral and isospectral heisenberg hierarchies abstract by means of the ablowitz kaup newell segur system with symbolic computation the gauge transformation between the first order nonisospectral and isospectral heisenberg spectral problems is obtained correspondingly the equivalence relation between the first order nonisospectral and isospectral heisenberg hierarchies is given in addition further applications and potential research directions of our method are discussed

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

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

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

multidimensional upwind methods for hyperbolic conservation laws abstract we present a class of second order conservative finite difference algorithms for solving numerically time dependent problems for hyperbolic conservation laws in several space variables these methods are upwind and multidimensional in that the numerical fluxes are obtained by solving the characteristic form of the full multidimensional equations at the zone edge and that all fluxes are evaluated and differenced at the same time in particular operator splitting is not used correct behavior at discontinuities is obtained by the use of solutions to the riemann problem and by limiting some of the second order terms numerical results are presented which show that the methods described here yield the same high resolution as the corresponding operator split methods

darboux transformation and hamiltonian structure for the jaulent miodek hierarchy abstract starting from the jaulent miodek spectral problem we derive the associated hierarchy of nonlinear evolution equations in this paper it is shown that this hierarchy is completely integrable in the liouville sense and possesses the hamiltonian structure moreover by virtue of symbolic computation two types of darboux transformations for the whole hierarchy are explicitly constructed which enables us to find the new soliton like shock and anti shock solutions for the jaulent miodek hierarchy figures are presented to discuss the properties of the new soliton like solutions

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

gauge transformations of constrained kp flows new integrable hierarchies integrable systems in 1 1 dimensions arise from the kp hierarchy as symmetry reductions involving square eigenfunctions exploiting the residual gauge freedom in these constraints new integrable systems are derived they include generalizations of the hierarchy of the kundu eckhaus equation and higher order extensions of the yajima oikawa and melnikov hierarchies constrained modified kp flows yield further integrable equations such as the hierarchies of the derivative nls equation the gerdjikov ivanov equation and the chen lee liu equation

on the nonlocal equations and nonlocal charges associated with the harry dym hierarchy a large class of nonlocal equations and nonlocal charges for the harry dym hierarchy is exhibited they are obtained from nonlocal casimirs associated with its bi hamiltonian structure the lax representation for some of these equations is also given

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

gauge transformations and symmetries of integrable systems we analyse several integrable systems in zero curvature form within the framework of invariant gauge theory in the drinfeld sokolov gauge we derive a two parameter family of nonlinear evolution equations which as special cases include the korteweg de vries kdv and harry dym hd equations we find residual gauge transformations which lead to infinitesimal symmetries of this family of equations for kdv and hd equations we find an infinite hierarchy of such symmetry transformations and we investigate their relation with local conservation laws constants of the motion and the bi hamiltonian structure of the equations applying successive gauge transformations of miura type we obtain a sequence of gauge equivalent integrable systems among them the modified kdv and calogero kdv equations

on the birkhoff factorization problem for the heisenberg magnet and nonlinear schrodinger equations a geometrical description of the heisenberg magnet hm equation with classical spins is given in terms of flows on the quotient space g h where g is an infinite dimensional lie group and h is a subgroup of g it is shown that the hm flows are induced by an action of r2 on g h and that the hm equation can be integrated by solving a birkhoff factorization problem for g for the hm flows that are laurent polynomials in the spectral variable we derive an algebraic transformation between solutions of the nonlinear schrodinger nls and heisenberg magnet equation the birkhoff factorization problem for g is treated in terms of the geometry of the segal wilson grassmannian gr h the solution of the problem is given in terms of a pair of baker functions for special subspaces in gr h the baker functions are constructed explicitly for subspaces that yield multisoliton solutions of nls and hm equations

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 wave propagation algorithm for viscoelastic fluids with spatially and temporally varying properties an algorithm is presented for simulating the fluid and stress equations that arise in a continuum model of platelet aggregation a l fogelson r d guy platelet wall interactions in continuum models of platelet thrombosis formulation and numerical solution math med biol 21 2004 293 334 the model is equivalent to models of viscoelastic fluids with spatially and temporally varying material parameters a subsystem containing the elastic terms is handled using a wave propagation algorithm two different methods for propagating waves are compared in computational tests one method linearizes the equations about cell edges and the other is based on the propagation of waves in a heterogeneous elastic medium when the viscoelastic fluid is in contact with a newtonian fluid the method based on linearization about the edges gives more accurate results no high weissenberg number instabilities were observed as long as the time step is chosen to obey a cfl condition the algorithm is stable and no unphysical oscillations appear in the solution 2007 elsevier b v all rights reserved

derivation and implicit solution of the harry dym equation and its connections with the korteweg de vries equation the harry dym equation which is related to the classical string problem is derived in three different ways an implicit cusp solitary wave solution is constructed via a simple direct method the existing connections between the harry dym and the korteweg de vries equations are uniformised and simplified and transformations between their respective solutions are carried out explicitly whenever possible physical insights are provided

stability of a viscoelastic burgers flow the system of equations proposed by burgers to model turbulent flow in a channel is extended to include viscoelastic affects the stability and bifurcation properties are examined in the neighborhood of the critical reynolds number for highly elastic fluids the bifurcated state is periodic with a shift in frequency

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

the k m n equation with generalized evolution term studied by symmetry reductions and qualitative analysis abstract in this paper we obtain symmetry reductions of the k m n equation with generalized evolution term the reduction to ordinary differential equations comes from an optimal system of subalgebras some of these equations admit symmetries which lead to further reductions and one of them comes out suitable for qualitative analysis its dynamical behavior is fully described and conservative quantities are stated

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

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 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

traveling waves and shocks in a viscoelastic generalization of burgers equation we consider traveling wave phenomena for a viscoelastic generalization of burgers equation for asymptotically constant velocity profiles we find three classes of solutions corresponding to smooth traveling waves piecewise smooth waves and piecewise constant shock solutions each solution type is possible for a given pair of asymptotic limits and we characterize the dynamics in terms of the relaxation time and viscosity

a two phase area level line edge detector abstract in this paper a new feature detector is proposed to extract both edge and line points existing in a same image it takes an n n block as a processing unit and locates feature points through two phases classification and extraction in the classification phase based on some local and global properties of the feature points each block is classified as homogeneous a line or an edge block in the extraction phase for each line or edge block according to its block type the proposed approach takes a corresponding feature extractor to locate all feature points as the block type is known the extraction work can be speeded up and the extraction result will be more accurate furthermore since the type of extracted feature point is known this will make segmentation easier experimental results are also given to show the effectiveness of the proposed method

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

on detecting edges an edge in an image corresponds to a discontinuity in the intensity surface of the underlying scene it can be approximated by a piecewise straight curve composed of edgels i e short linear edge elements each characterized by a direction and a position the approach to edgel detection here is to fit a series of one dimensional surfaces to each window kernel of the operator and accept the surface description which is adequate in the least squares sense and has the fewest parameters a one dimensional surface is one which is constant along some direction the tanh is an adequate basis for the stepedge and its combinations are adequate for the roofedge and the line edge the proposed method of step edgel detection is robust with respect to noise for step size noise 2 5 it has subpixel position localization position and an angular localization better than 10 further it is designed to be insensitive to smooth shading these results are demonstrated by some simple analysis statistical data and edgelimages also included is a comparison of performance on a real image with a typical operator difference of gaussians the results indicate that the proposed operator is superior with respect to detection localization and resolution

reduction technique for matrix nonlinear evolution equations solvable by the spectral transform the main purpose of this paper is to describe a technique of reduction whereby from the class of evolution equations for matrices of order n solvable via the spectral transform associated to the matrix linear schrodinger eigenvalue problem one derives subclasses of nonlinear evolution equations involving less than n2 fields to illustrate the method from the equations for matrices of order 2 two subclasses of equations for 2 fields rather than 4 are obtained the first class coincides or rather includes that solvable via the spectral transform associated to the generalized zakharov shabat spectral problem further reduction to nonlinear evolution equations for a single field reproduces a number of well known equations but also yields a novel one highly nonlinear the second class also yields highly nonlinear equations some examples are given including another novel evolution equation for a single field

a computational approach to edge detection this paper describes a computational approach to edge detection the success of the approach depends on the definition of a comprehensive set of goals for the computation of edge points these goals must be precise enough to delimit the desired behavior of the detector while making minimal assumptions about the form of the solution we define detection and localization criteria for a class of edges and present mathematical forms for these criteria as functionals on the operator impulse response a third criterion is then added to ensure that the detector has only one response to a single edge we use the criteria in numerical optimization to derive detectors for several common image features including step edges on specializing the analysis to step edges we find that there is a natural uncertainty principle between detection and localization performance which are the two main goals with this principle we derive a single operator shape which is optimal at any scale the optimal detector has a simple approximate implementation in which edges are marked at maxima in gradient magnitude of a gaussian smoothed image we extend this simple detector using operators of several widths to cope with different signal to noise ratios in the image we present a general method called feature synthesis for the fine to coarse integration of information from operators at different scales finally we show that step edge detector performance improves considerably as the operator point spread function is extended along the edge

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

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

on modified nls kaup and nlbq equations differential transformations and bilinearization new transformations between the nonlinear schrodinger kaup and non local boussinesq equations as well as their modified counterparts are found and analysed the bilinear representations of these equations including an alternative bilinear form of the chen lee liu equation are obtained by a direct method based on the bell s exponential polynomials explicit wronskian solutions to these equations are also presented

miura and generalized backlund transformation for kdv hierarchy using the fact that miura transformation can be expressed in the form of gauge transformation connecting the kdv and mkdv equations we discuss the derivation of the backlund transformation and its miura gauge transformation connecting both hierarchies

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

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