Home About us Contact
|
Home About us Contact | ||
|
|
||
Pressure Poisson Equation (pressure + poisson_equation)
Selected AbstractsA monolithic approach for interaction of incompressible viscous fluid and an elastic body based on fluid pressure Poisson equationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 2 2005Daisuke Ishihara Abstract This paper describes a new monolithic approach based on the fluid pressure Poisson equation (PPE) to solve an interaction problem of incompressible viscous fluid and an elastic body. The PPE is derived so as to be consistent with the coupled equation system for the fluid-structure interaction (FSI). Based on this approach, we develop two kinds of efficient monolithic methods. In both methods, the fluid pressure is derived implicitly so as to satisfy the incompressibility constraint, and all other unknown variables are derived fully explicitly or partially explicitly. The coefficient matrix of the PPE for the FSI becomes symmetric and positive definite and its condition is insensitive to inhomogeneity of material properties. The arbitrary Lagrangian,Eulerian (ALE) method is employed for the fluid part in order to take into account the deformable fluid-structure interface. To demonstrate fundamental performances of the proposed approach, the developed two monolithic methods are applied to evaluate the added mass and the added damping of a circular cylinder as well as to simulate the vibration of a rectangular cylinder induced by vortex shedding. Copyright © 2005 John Wiley & Sons, Ltd. [source] A collocated, iterative fractional-step method for incompressible large eddy simulationINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2008Giridhar Jothiprasad Abstract Fractional-step methods are commonly used for the time-accurate solution of incompressible Navier,Stokes (NS) equations. In this paper, a popular fractional-step method that uses pressure corrections in the projection step and its iterative variants are investigated using block-matrix analysis and an improved algorithm with reduced computational cost is developed. Since the governing equations for large eddy simulation (LES) using linear eddy-viscosity-based sub-grid models are similar in form to the incompressible NS equations, the improved algorithm is implemented in a parallel LES solver. A collocated grid layout is preferred for ease of extension to curvilinear grids. The analyzed fractional-step methods are viewed as an iterative approximation to a temporally second-order discretization. At each iteration, a linear system that has an easier block-LU decomposition compared with the original system is inverted. In order to improve the numerical efficiency and parallel performance, modified ADI sub-iterations are used in the velocity step of each iteration. Block-matrix analysis is first used to determine the number of iterations required to reduce the iterative error to the discretization error of. Next, the computational cost is reduced through the use of a reduced stencil for the pressure Poisson equation (PPE). Energy-conserving, spatially fourth-order discretizations result in a 7-point stencil in each direction for the PPE. A smaller 5-point stencil is achieved by using a second-order spatial discretization for the pressure gradient operator correcting the volume fluxes. This is shown not to reduce the spatial accuracy of the scheme, and a fourth-order continuity equation is still satisfied to machine precision. The above results are verified in three flow problems including LES of a temporal mixing layer. Copyright © 2008 John Wiley & Sons, Ltd. [source] Pressure segregation methods based on a discrete pressure Poisson equation.INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2008An algebraic approach Abstract In this paper, we introduce some pressure segregation methods obtained from a non-standard version of the discrete monolithic system, where the continuity equation has been replaced by a pressure Poisson equation obtained at the discrete level. In these methods it is the velocity instead of the pressure the extrapolated unknown. Moreover, predictor,corrector schemes are suggested, again motivated by the new monolithic system. Key implementation aspects are discussed, and a complete stability analysis is performed. We end with a set of numerical examples in order to compare these methods with classical pressure-correction schemes. Copyright © 2007 John Wiley & Sons, Ltd. [source] Pressure boundary condition for the time-dependent incompressible Navier,Stokes equationsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2006R. L. Sani Abstract In Gresho and Sani (Int. J. Numer. Methods Fluids 1987; 7:1111,1145; Incompressible Flow and the Finite Element Method, vol. 2. Wiley: New York, 2000) was proposed an important hypothesis regarding the pressure Poisson equation (PPE) for incompressible flow: Stated there but not proven was a so-called equivalence theorem (assertion) that stated/asserted that if the Navier,Stokes momentum equation is solved simultaneously with the PPE whose boundary condition (BC) is the Neumann condition obtained by applying the normal component of the momentum equation on the boundary on which the normal component of velocity is specified as a Dirichlet BC, the solution (u, p) would be exactly the same as if the ,primitive' equations, in which the PPE plus Neumann BC is replaced by the usual divergence-free constraint (, · u = 0), were solved instead. This issue is explored in sufficient detail in this paper so as to actually prove the theorem for at least some situations. Additionally, like the original/primitive equations that require no BC for the pressure, the new results establish the same thing when the PPE approach is employed. Copyright © 2005 John Wiley & Sons, Ltd. [source] Computation of unsteady viscous incompressible flows in generalized non-inertial co-ordinate system using Godunov-projection method and overlapping meshesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2002H. Pan Abstract Time-dependent incompressible Navier,Stokes equations are formulated in generalized non-inertial co-ordinate system and numerically solved by using a modified second-order Godunov-projection method on a system of overlapped body-fitted structured grids. The projection method uses a second-order fractional step scheme in which the momentum equation is solved to obtain the intermediate velocity field which is then projected on to the space of divergence-free vector fields. The second-order Godunov method is applied for numerically approximating the non-linear convection terms in order to provide a robust discretization for simulating flows at high Reynolds number. In order to obtain the pressure field, the pressure Poisson equation is solved. Overlapping grids are used to discretize the flow domain so that the moving-boundary problem can be solved economically. Numerical results are then presented to demonstrate the performance of this projection method for a variety of unsteady two- and three-dimensional flow problems formulated in the non-inertial co-ordinate systems. Copyright © 2002 John Wiley & Sons, Ltd. [source] An approximate projection method for incompressible flowINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2002David E. Stevens This paper presents an approximate projection method for incompressible flows. This method is derived from Galerkin orthogonality conditions using equal-order piecewise linear elements for both velocity and pressure, hereafter Q1Q1. By combining an approximate projection for the velocities with a variational discretization of the continuum pressure Poisson equation, one eliminates the need to filter either the velocity or pressure fields as is often needed with equal-order element formulations. This variational approach extends to multiple types of elements; examples and results for triangular and quadrilateral elements are provided. This method is related to the method of Almgren et al. (SIAM J. Sci. Comput. 2000; 22: 1139,1159) and the PISO method of Issa (J. Comput. Phys. 1985; 62: 40,65). These methods use a combination of two elliptic solves, one to reduce the divergence of the velocities and another to approximate the pressure Poisson equation. Both Q1Q1 and the method of Almgren et al. solve the second Poisson equation with a weak error tolerance to achieve more computational efficiency. A Fourier analysis of Q1Q1 shows that a consistent mass matrix has a positive effect on both accuracy and mass conservation. A numerical comparison with the widely used Q1Q0 (piecewise linear velocities, piecewise constant pressures) on a periodic test case with an analytic solution verifies this analysis. Q1Q1 is shown to have comparable accuracy as Q1Q0 and good agreement with experiment for flow over an isolated cubic obstacle and dispersion of a point source in its wake. Copyright © 2002 John Wiley & Sons, Ltd. [source] Three-dimensional numerical modelling of free surface flows with non-hydrostatic pressureINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2002Musteyde B. Koçyigit Abstract A three-dimensional numerical model is developed for incompressible free surface flows. The model is based on the unsteady Reynolds-averaged Navier,Stokes equations with a non-hydrostatic pressure distribution being incorporated in the model. The governing equations are solved in the conventional sigma co-ordinate system, with a semi-implicit time discretization. A fractional step method is used to enable the pressure to be decomposed into its hydrostatic and hydrodynamic components. At every time step one five-diagonal system of equations is solved to compute the water elevations and then the hydrodynamic pressure is determined from a pressure Poisson equation. The model is applied to three examples to simulate unsteady free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field. Emphasis is focused on applying the model to wave problems. Two of the examples are about modelling small amplitude waves where the hydrostatic approximation and long wave theory are not valid. The other example is the wind-induced circulation in a closed basin. The numerical solutions are compared with the available analytical solutions for small amplitude wave theory and very good agreement is obtained. Copyright © 2002 John Wiley & Sons, Ltd. [source] Numerical simulation of high-Reynolds number flow around circular cylinders by a three-step FEM,BEM modelINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2001D. L. Young Abstract An innovative computational model, developed to simulate high-Reynolds number flow past circular cylinders in two-dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier,Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three-step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high-Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd. [source] Interface tracking finite volume method for complex solid,fluid interactions on fixed meshesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 2 2002H. S. Udaykumar Abstract We present a numerical technique for computing flowfields around moving solid boundaries immersed in fixed meshes. The mixed Eulerian,Lagrangian framework treats the immersed boundaries as sharp solid,fluid interfaces and a conservative finite volume formulation allows boundary conditions at the moving surfaces to be exactly applied. A semi-implicit second-order accurate spatial and temporal discretization is employed with a fractional-step scheme for solving the flow equations. A multigrid accelerator for the pressure Poisson equations has been developed to apply in the presence of multiple embedded solid regions on the mesh. We present applications of the method to two types of problems: (a) solidification in the presence of flows and particles, (b) fluid,structure interactions in flow control. In both these problems, the sharp interface method presents advantages by being able to track arbitrary interface motions, while capturing the full viscous, unsteady dynamics. Copyright © 2001 John Wiley & Sons, Ltd. [source] | ||||||||||