A second-order cut-cell method for the numerical simulation of 2D flows past obstacles

Pergamon-Elsevier Science Ltd, 2012
Article de conférence

We present a new second-order method, based on the MAC scheme on cartesian grids, for the numerical simulation of two-dimensional incompressible flows past obstacles. In this approach, the solid boundary is embedded in the cartesian computational mesh. Discretizations of the viscous and convective terms are formulated in the context of finite volume methods ensuring local conservation properties of the scheme. Classical second-order centered schemes are applied in mesh cells which are sufficiently far from the obstacle. In the mesh cells cut by the obstacle, first-order approximations are proposed. The resulting linear system is nonsymmetric but the stencil remains local as in the classical MAC scheme on cartesian grids. The linear systems are solved by a fast direct method based on the capacitance matrix method. The time integration is achieved with a second-order projection scheme. While in cut-cells the scheme is locally first-order, a global second-order accuracy is recovered. This property is assessed by computing analytical solutions for a Taylor-Couette problem. The efficiency and robustness of the method is supported by numerical simulations of 2D flows past a circular cylinder at Reynolds number up to 9500. Good agreement with experimental and published numerical results are obtained. (c) 2012 Elsevier Ltd. All rights reserved.

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

Pergamon-Elsevier Science Ltd, 2012
Article de conférence

Résumé

Official source

https://infoscience.epfl.ch/record/184528?ln=fr

À propos de ce résultat

Concepts associés (14)

Concepts associés

Chargement

Publications associées (30)

Publications associées

Chargement

Simulation informatique

vignette|upright=1|Une simulation informatique, sur une étendue de , de l'évolution du typhon Mawar produite par le Modèle météorologique Weather Research and Forecasting La simulation informatique

Analyse numérique

L’analyse numérique est une discipline à l'interface des mathématiques et de l'informatique. Elle s’intéresse tant aux fondements qu’à la mise en pratique des méthodes permettant de résoudre, par des

Direct numerical simulation

A direct numerical simulation (DNS) is a simulation in computational fluid dynamics (CFD) in which the Navier–Stokes equations are numerically solved without any turbulence model. This means that the

, ,

A numerical model is presented for the simulation of complex fluid flows with free surfaces. The unknowns are the velocity and pressure fields in the liquid region, together with a function defining the volume fraction of liquid. Although the mathematical formulation of the model is similar to the volume of fluid (VOF) method, the numerical schemes used to solve the problem are different. A splitting method is used for the time discretization. At each time step, two advection problems and a generalized Stokes problem are to be solved. Two different grids are used for the space discretization. The two advection problems are solved on a fixed, structured grid made out of small rectangular cells, using a forward characteristic method. The generalized Stokes problem is solved using a finite element method on a fixed, unstructured mesh. Numerical results are presented for several test cases: the filling of an S-shaped channel, the filling of a disk with core, the broken dam in a confined domain. (C) 1999 Academic Press.

1999

Chargement

Numerical simulation of three-dimensional free surface flows

Vincent Maronnier, Marco Picasso, Jacques Rappaz

A numerical model is presented for the simulation of complex fluid flows with free surfaces in three space dimensions. The model described in Maronnier et al. (J. Comput. Phys. 1999; 155(2):439) is extended to three dimensional situations. The mathematical formulation of the model is similar to that of the volume of fluid (VOF) method, but the numerical procedures are different. A splitting method is used for the time discretization. At each time step, two advection problems-one for the predicted velocity field and the other for the volume fraction of liquid-are to be solved. Then, a generalized Stokes problem is solved and the velocity field is corrected. Two different grids are used for the space discretization. The two advection problems are solved on a fixed, structured grid made out of small cubic cells, using a forward characteristic method. The generalized Stokes problem is solved using continuous, piecewise linear stabilized finite elements on a fixed, unstructured mesh of tetrahedrons. The three-dimensional implementation is discussed. Efficient postprocessing algorithms enhance the quality of the numerical solution. A hierarchical data structure reduces memory requirements. Numerical results are presented for complex geometries arising in mold filling. Copyright (C) 2003 John Wiley Sons, Ltd.

2003

Chargement

Numerical Simulation of Macrosegregation : a Comparison between FVM and FEM Predictions and a Confrontation with Experiments

Jean-Luc Desbiolles, Jacques Rappaz, Michel Rappaz

Micro-macrosegregation calculations have been performed for a rectangular cavity containing either a Pb-48 wt pct Sn alloy or a Sn-5 wt pct Pb alloy. The numerical results calculated with a finite volume method (FVM) and a finite element method (FEM) are compared with experimental results previously obtained by Hebditch and Hunt.([1]) The two methods are based on the same average conservation equations governing heat and mass transfer and the same assumptions: lever rule, equal and constant density of the solid and liquid phases (except in the buoyancy term), permeability of the mushy zone given by the Carman-Kozeny relation, and no transport of the solid phase. Although the same parameters are used in both calculations, small differences are observed as a result of the different formulations. In particular, the instabilities appearing in the mushy zone (channels) of the Sn-5 wt pct Pb alloy are more pronounced with the FVM formulation as compared with FEM, whereas the opposite trend is observed for the Pb-48 wt pct Sn alloy. Nevertheless, the final segregation maps at the end of solidification compare fairly well with the experimental findings.

1998

Chargement