共查询到20条相似文献，搜索用时 0 毫秒
1.
In this paper, we describe an implicit hybrid finite volume (FV)/element (FE) incompressible Navier–Stokes solver for turbulent flows based on the Spalart–Allmaras detached eddy simulation (SA‐DES). The hybrid FV/FE solver is based on the segregated pressure correction or projection method. The intermediate velocity field is first obtained by solving the original momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centers and the auxiliary variable at vertices, making the current solver a staggered‐mesh scheme. The SA‐DES turbulence equation is solved after the velocity and the pressure fields have been updated at the end of each time step. The same matrix‐free FV method as the one used for momentum equations is used to solve the turbulence equation. The turbulence equation provides the eddy viscosity, which is added to the molecular viscosity when solving the momentum equation. In our implementation, we focus on the accuracy, efficiency and robustness of the SA‐DES model in a hybrid flow solver. This paper will address important implementation issues for high‐Reynolds number flows where highly stretched elements are typically used. In addition, some aspects of implementing the SA‐DES model will be described to ensure the robustness of the turbulence model. Several numerical examples including a turbulent flow past a flat plate and a high‐Reynolds number flow around a high angle‐of‐attack NACA0015 airfoil will be presented to demonstrate the accuracy and efficiency of our current implementation. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
2.
We present a Lagrangian formulation for finite element analysis of quasi‐incompressible fluids that has excellent mass preservation features. The success of the formulation lays on a new residual‐based stabilized expression of the mass balance equation obtained using the finite calculus method. The governing equations are discretized with the FEM using simplicial elements with equal linear interpolation for the velocities and the pressure. The merits of the formulation in terms of reduced mass loss and overall accuracy are verified in the solution of 2D and 3D quasi‐incompressible free‐surface flow problems using the particle FEM ( www.cimne.com/pfem ). Examples include the sloshing of water in a tank, the collapse of one and two water columns in rectangular and prismatic tanks, and the falling of a water sphere into a cylindrical tank containing water. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
3.
A hybrid particle‐mesh method was developed for efficient and accurate simulations of two‐phase flows. In this method, the main component of the flow is solved using the constrained interpolated profile/multi‐moment finite volumemethod; the two‐phase interface is rendered using the finite volume particle (FVP) method. The effect of surface tension is evaluated using the continuum surface force model. Numerical particles in the FVP method are distributed only on the surface of the liquid in simulating the interface between liquid and gas; these particles are used to determine the density of each mesh grid. An artificial term was also introduced to mitigate particle clustering in the direction of maximum compression and sparse discretization errors in the stretched direction. This enables accurate interface tracking without diminishing numerical efficiency. Two benchmark simulations are used to demonstrate the validity of the method developed and its numerical stability. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
4.
Shahrouz Aliabadi Christopher Bigler Erdal Yilmaz Sridhar Palle Bela Soni 《International Journal of Computational Fluid Dynamics》2013,27(4):175-189
An implicit hybrid finite element (FE)/volume solver has been extended to incompressible flows coupled with the energy equation. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centred finite volume (FV) method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centres and the auxiliary variable at vertices. The Generalized Minimal Residual (GMRES) matrix-free strategy is adapted to solve the governing equations in both FE and FV methods. The presented 2D and 3D numerical examples show the robustness and accuracy of the numerical method. 相似文献
5.
Numerical issues arising in computations of viscous flows in corners formed by a liquid–fluid free surface and a solid boundary are considered. It is shown that on the solid a Dirichlet boundary condition, which removes multivaluedness of velocity in the ‘moving contact‐line problem’ and gives rise to a logarithmic singularity of pressure, requires a certain modification of the standard finite‐element method. This modification appears to be insufficient above a certain critical value of the corner angle where the numerical solution becomes mesh‐dependent. As shown, this is due to an eigensolution, which exists for all angles and becomes dominant for the supercritical ones. A method of incorporating the eigensolution into the numerical method is described that makes numerical results mesh‐independent again. Some implications of the unavoidable finiteness of the mesh size in practical applications of the finite‐element method in the context of the present problem are discussed. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
6.
A numerical model has been developed for the 2D simulation of free surface flows or, more generally speaking, moving interface ones. The bulk fluids on both sides of the interface are taken into account in simulating the incompressible laminar flow state. In the case of heat transfer the whole system, i.e. walls as well as possible obstacles, is considered. This model is based on finite element analysis with an Eulerian approach and an unstructured fixed mesh. A special technique to localize the interface allows its temporal evolution through this mesh. Several numerical examples are presented to demonstrate the capabilities of the model. © 1997 by John Wiley & Sons, Ltd. 相似文献
7.
In this work we present a numerical method for solving the incompressible Navier–Stokes equations in an environmental fluid mechanics context. The method is designed for the study of environmental flows that are multiscale, incompressible, variable‐density, and within arbitrarily complex and possibly anisotropic domains. The method is new because in this context we couple the embedded‐boundary (or cut‐cell) method for complex geometry with block‐structured adaptive mesh refinement (AMR) while maintaining conservation and second‐order accuracy. The accurate simulation of variable‐density fluids necessitates special care in formulating projection methods. This variable‐density formulation is well known for incompressible flows in unit‐aspect ratio domains, without AMR, and without complex geometry, but here we carefully present a new method that addresses the intersection of these issues. The methodology is based on a second‐order‐accurate projection method with high‐order‐accurate Godunov finite‐differencing, including slope limiting and a stable differencing of the nonlinear convection terms. The finite‐volume AMR discretizations are based on two‐way flux matching at refinement boundaries to obtain a conservative method that is second‐order accurate in solution error. The control volumes are formed by the intersection of the irregular embedded boundary with Cartesian grid cells. Unlike typical discretization methods, these control volumes naturally fit within parallelizable, disjoint‐block data structures, and permit dynamic AMR coarsening and refinement as the simulation progresses. We present two‐ and three‐dimensional numerical examples to illustrate the accuracy of the method. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
8.
In the present work a finite‐difference technique is developed for the implementation of a new method proposed by Aristov and Pukhnachev (Doklady Phys. 2004; 49 (2):112–115) for modeling of the axisymmetric viscous incompressible fluid flows. A new function is introduced that is related to the pressure and a system similar to the vorticity/stream function formulation is derived for the cross‐flow. This system is coupled to an equation for the azimuthal velocity component. The scheme and the algorithm treat the equations for the cross‐flow as an inextricably coupled system, which allows one to satisfy two conditions for the stream function with no condition on the auxiliary function. The issue of singularity of the matrix is tackled by adding a small parameter in the boundary conditions. The scheme is thoroughly validated on grids with different resolutions. The new numerical tool is applied to the Taylor flow between concentric rotating cylinders when the upper and lower lids are allowed to rotate independently from the inner cylinder, while the outer cylinder is held at rest. The phenomenology of this flow is adequately represented by the numerical model, including the hysteresis that takes place near certain specific values of the Reynolds number. Thus, the present results can be construed to demonstrate the viability of the new model. The success can be attributed to the adequate physical nature of the auxiliary function. The proposed technique can be used in the future for in‐depth investigations of the bifurcation phenomena in rotating flows. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
9.
Robert Jurjevi&#x; 《国际流体数值方法杂志》1999,31(3):601-626
In this paper, a Galerkin weighted residual finite element numerical solution method, with velocity material time derivative discretisation, is applied to solve for a classical fluid mechanics system of partial differential equations modelling two‐dimensional stationary incompressible Newtonian fluid flow. Classical examples of driven cavity laminar flow and laminar flow past a cylinder are presented. Numerical results are compared with data found in the literature. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
10.
In this paper we present a problem we have encountered using a stabilized finite element method on fixed grids for flows with interfaces modelled with the level set approach. We propose a solution based on enriching the pressure shape functions on the elements cut by the interface. The enrichment is used to enable the pressure gradient to be discontinuous at the interface, thus improving the ability to simulate the behaviour of fluids with different density under a gravitational force. The additional shape function used is local to each element and the corresponding degree of freedom can therefore be condensed prior to assembly, making the implementation quite simple on any existing finite element code. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
11.
Samuel P. Schofield Mark A. Christon Vadim Dyadechko Rao V. Garimella Robert B. Lowrie Blair K. Swartz 《国际流体数值方法杂志》2010,63(8):931-952
This paper compares the numerical performance of the moment‐of‐fluid (MOF) interface reconstruction technique with Youngs, LVIRA, power diagram (PD), and Swartz interface reconstruction techniques in the context of a volume‐of‐fluid (VOF) based finite element projection method for the numerical simulation of variable‐density incompressible viscous flows. In pure advection tests with multiple materials MOF shows dramatic improvements in accuracy compared with the other methods. In incompressible flows where density differences determine the flow evolution, all the methods perform similarly for two material flows on structured grids. On unstructured grids, the second‐order MOF, LVIRA, and Swartz methods perform similarly and show improvement over the first‐order Youngs' and PD methods. For flow simulations with more than two materials, MOF shows increased accuracy in interface positions on coarse meshes. In most cases, the convergence and accuracy of the computed flow solution was not strongly affected by interface reconstruction method. Published in 2009 by John Wiley & Sons, Ltd. 相似文献
12.
In previous studies, the moment‐of‐fluid interface reconstruction method showed dramatic accuracy improvements in static and pure advection tests over existing methods, but this did not translate into an equivalent improvement in volume‐tracked multimaterial incompressible flow simulation using low‐order finite elements. In this work, the combined effects of the spatial discretization and interface reconstruction in flow simulation are examined. The mixed finite element pairs, Q1Q0 (with pressure stabilization) and Q2P ? 1 are compared. Material order‐dependent and material order‐independent first and second‐order accurate interface reconstruction methods are used. The Q2P ? 1 elements show significant improvements in computed flow solution accuracy for single material flows but show reduced convergence using element‐average piecewise constant density and viscosity in volume‐tracked simulations. In general, a refined Q1Q0 grid, with better material interface resolution, provided an accuracy similar to the Q2P ? 1 element grid with a comparable number of degrees of freedom. Moment‐of‐fluid shows more benefit from the higher‐order accurate flow simulation than the LVIRA, Youngs', and power diagram interface reconstruction methods, especially on unstructured grids, but does not recover the dramatic accuracy improvements it has shown in advection tests. Published 2012. This article is a US Government work and is in the public domain in the USA. 相似文献
13.
A non‐linear method, PREC, for computation of the movement of a free surface is proposed here. The method is composed of three steps: identifying the free surface by using a non‐linear function from the volume fraction matrix, updating the volume fraction matrix using a volume projection method with error correction, and treatment of the results using overshooting or undershooting. Identification of the free surface includes using a polynomial function with 2, 4, or 8 coefficients for one‐, two‐, or three‐dimensional problems, respectively. The polynomial reconstruction involves non‐negligible numerical error. The second advection step includes a linear projection method in space and time. Advection of the volume fraction matrix is computed from the occupying volume of the mesh at the previous time step. At the new time step, the error at each grid point is assumed to be similar to the error at the previous time step and is used for correction. Overshooting or undershooting develops around the free surface mesh points due to the solution's finite time increment. The third step includes truncating the numerical overshooting or undershooting volumes, i.e. isotropic spreading of the excess fluid volumes. The PREC method is evaluated for a one‐dimensional flow case and several two‐dimensional simple flow cases with circular sections (cases include transition parallel to a coordinate, transition with an intersection angle to a coordinate, and rotation). The results from the present method are compared with analytical solutions and results from a donor‐cell VOF method. As a result of these comparisons, the PREC method is validated. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
14.
The paper's leitmotiv is condensed in one word: robustness. This is a real hindrance for the successful implementation of any multigrid scheme for solving the Navier–Stokes set of equations. In this paper, many hints are given to improve this issue. Instead of looking for the best possible speed‐up rate for a particular set of problems, at a given regime and in a given condition, the authors propose some ideas pursuing reasonable speed‐up rates in any situation. In a previous paper, the authors presented a multigrid method for solving the incompressible turbulent RANS equations, with particular care in the robustness and flexibility of the solution scheme. Here, these concepts are further developed and extended to compressible laminar and turbulent flows. This goal is achieved by introducing a non‐linear multigrid scheme for compressible laminar (NS equations) and turbulent flow (RANS equations), taking benefit of a convenient master–slave implementation strategy. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
15.
In this paper, we report our development of an implicit hybrid flow solver for the incompressible Navier–Stokes equations. The methodology is based on the pressure correction or projection method. A fractional step approach is used to obtain an intermediate velocity field by solving the original momentum equations with the matrix‐free implicit cell‐centred finite volume method. The Poisson equation derived from the fractional step approach is solved by the node‐based Galerkin finite element method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centres and the auxiliary variable at cell vertices, making the current solver a staggered‐mesh scheme. Numerical examples demonstrate the performance of the resulting hybrid scheme, such as the correct temporal convergence rates for both velocity and pressure, absence of unphysical pressure boundary layer, good convergence in steady‐state simulations and capability in predicting accurate drag, lift and Strouhal number in the flow around a circular cylinder. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
16.
By treating it as a contact discontinuity in the density field, a free surface between two immiscible fluids can be automatically ‘captured’ by the enforcement of conservation laws. A surface‐capturing method of this kind requires no special tracking or fitting treatment for the free surface, thereby offering the advantage of algorithm simplicity over the surface‐tracking or the surface‐fitting method. A surface‐capturing method based on a new multi‐fluid incompressible Navier–Stokes formulation is developed. It is applied to a variety of free‐surface flows, including the Rayleigh–Taylor instability problem, the ship waves around a Wigley hull and a model bubble‐rising problem to demonstrate the validity and versatility of the present method. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
17.
Two‐phase flows around fluid particles are often considered to be in infinite domains, to avoid influence of the domain walls. Numerical simulations, however, must be modeled with a bounded domain, thus introducing artificial boundaries. Modeling of fluid flow in a domain with such artificial boundaries requires a careful choice of suitable boundary conditions. Slip boundary conditions for example can have a large impact on the computational results if the domain is chosen to be too small, because they model impermeable walls. This paper introduces an artificial boundary condition for simulations of the flow around single rising or settling fluid particles based on the approximated decay behavior of the velocity and the pressure field in the surrounding liquid. This is applied to the simulation of rising gas bubbles in systems with a Reynolds number of up to 50, and the outcome is compared with experimental results and simulations with slip boundary condition. It is found that domain size can be reduced by a factor of about two compared with slip boundary conditions without loss of accuracy. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
18.
In this paper, we present a new method for simulating the motion of a disperse particle phase in a carrier gas through porous media. We assume a sufficiently dilute particle‐laden flow and compute, independently of the disperse phase, the steady laminar fluid velocity using the immersed boundary method. Given the velocity of the carrier gas, the equations of motion for the particles experiencing the Stokes drag force are solved to determine their trajectories. The ‘no‐slip consistent’ particle tracking algorithm avoids possible numerical filtration of very small particles due to the nonzero velocity field at the solid–fluid interface introduced by the immersed boundary method. This physically consistent tracking allows a reliable estimation of the filtration efficiency of porous filters due to inertial impaction. We illustrate and test our new approach for model porous media consisting of a structured array of aligned rectangular fibers, arranged in line and staggered. In the staggered geometry, the effect of the residual velocity at the solid–fluid interface is significant for particles with low inertia. Without adopting the developed no‐slip consistent numerical method, an artificial numerical filtration is observed, which becomes dominant for small enough particles. For both the in line and the staggered geometries, the filtration rate depends quite strongly and non monotonically on the particle inertia. This is expressed most clearly in the staggered arrangement in which a very strong increase in the filtration efficiency is observed at a well‐defined critical droplet size, corresponding to a qualitative change in the dominant particle paths in the porous medium. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
19.
The coupling between the equations governing the free‐surface flows, the six degrees of freedom non‐linear rigid body dynamics, the linear elasticity equations for mesh‐moving and the cables has resulted in a fluid‐structure interaction technology capable of simulating mooring forces on floating objects. The finite element solution strategy is based on a combination approach derived from fixed‐mesh and moving‐mesh techniques. Here, the free‐surface flow simulations are based on the Navier–Stokes equations written for two incompressible fluids where the impact of one fluid on the other one is extremely small. An interface function with two distinct values is used to locate the position of the free‐surface. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian–Eulerian domain. This allows us to handle the motion of the time dependent geometries. Forces and momentums exerted on the floating object by both water and hawsers are calculated and used to update the position of the floating object in time. In the mesh moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The non‐linear rigid body dynamics equations are coupled with the governing equations of fluid flow and are solved simultaneously to update the position of the floating object. The numerical examples includes a 3D simulation of water waves impacting on a moored floating box and a model boat and simulation of floating object under water constrained with a cable. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
