Differences

This shows you the differences between two versions of the page.

--- tailorcrete:examples:v-funnel [2012/07/25 11:25] – kolarfil
+++ tailorcrete:examples:v-funnel [2012/09/13 09:14] (current) – bp
@@ Line 1: / Line 1: @@
-==== V-funnel simulation ====
+===== V-funnel simulation =====
-=== Test setup and Geometry ===
+==== Test setup and Geometry ====
-V-funnel problem is one of the first problem, that was modeled. It is focused on boundary condition influence. In it's easiest form is modeled in 2D. The geometry is shown at figure below. Characteristic time can be defined for V-funnel simulation. It is the time between opening the gate and the moment, when it is posible to look through the gate.
+V-funnel problem simulates the emtying of SCC from V-shaped reservoir due to gravity force. The thickness of the specimen is rather small and the effect of the thickness would be negligible, so that the analysis is performed in 2D. This test can be conviniently used to evaluate the influence of different boundary conditions applied on slip walls. The geometry of the test setup is shown in Figure 1a. The basic parameter, compared to the experiments is the time between opening the gate and the moment, when it is posible to look through the gate.
-{{:tailorcrete:examples:v_funnel_geometry.png?200|}}
+|{{:tailorcrete:examples:v_funnel_geometry.png?300|}} | {{ :tailorcrete:examples:v_funnel_mesh_bc.png?300|}} |
+|  Fig. 1.a: Geometry of the test  |  Fig. 1.b: Boundary conditions  |
-=== Computational Model ===
+\\
-Scheme of two-dimensional computational model is shown on the picture below. On the picture is also specification of boundary condition. The whole domain is decomposed into two regions. Initially, in the first region is air modeled as a newtonian fluid. In the second region is concrete, considered as a Bingham fluid. Parameters are therefore yield stress 60 Pa and plastic viscosity 20 Pa*s. So called "Slip with friction" boundary condition was used on walls, when tangential component of the velocity is proportional to the tangential component of the stress tensor.
-{{:tailorcrete:examples:v_funnel_mesh_bc.png?200|}}
+==== Computational Model ====
+{{ :tailorcrete:v-funnel_geometry.png?500|Fig.2: Computational mesh}}
+The setup of 2D computational model is illustrated in Fig. 1b. The problem is modeled as a two-phase flow problem, consisting of two inmissible fluids: one representing the self compacting concrete and the other representing air. The evolving interface between the two fluids is tracked using the Level Set Method, see reference [1] for details. The Figure 1.b also illustrates the applied boundary boundary conditions. At the beginning of the simulation, the problem domain is divided into two regions, one initially filled with SCC (modelled as non-Newtonian fluid using Bingham model) and another one on the top of SCC representing air (modelled as newtonian fuid). The following constitutive parameters for Bingham model of SCC were assumed in the simulation:
+  *yield stress 60 [Pa]
+  *plastic viscosity 20 [Pa*s]
+So called "Slip with friction" boundary condition was used on walls, when tangential component of the velocity is proportional to the tangential component of the stress tensor.
-On the next picture is mesh used for simulation. It is devided into two parts. Lower part is finer to improve accuracy of the flow near the neck. The mesh contains 2907 nodes and 5567 elements.
+The next picture ilustrates the computational mesh. The lower part is refined in order to improve accuracy of predicted flow pattern near the neck. The mesh contains 2907 nodes and 5567 tringular elements with linear interpolation (same approximation used for both velocity and pressure field). Since such element is not satisfying LBB condition, PSPG stabilization is used for preventing oscilations in pressure field. SUPG stabilization improving accuracy in connection with non-linear convective term is also used. For further information, see [2].
-{{:tailorcrete:v-funnel_geometry.png?200|}}
+The link below contains the complete input file for OOFEM code of the V-funnel simulation. Basically, it is a text file containing specifications of geometry (mesh), physics (constitutive model) and numerical solver.
-The link below contains input file of V-funnel simulation for OOFEM program. It is a text file containing specifications about geometry (mesh), physics (constitutive model) and numerical solver.
 {{:tailorcrete:examples:v_funnel_tau60_mu20_oofem.in|}}
-Description of input file can be found here: [[tailorcrete:input_file_description|Input File Description]]
+The general input specification can be found here: [[tailorcrete:input_file_description|Input File Description]].
-=== Results ===
+==== Results ====
 On the video below, motion of concrete-air interface is shown. The characteristic "V-funnel time" (the time from opening the gate at the bottom to the moment when it is possible to look through the gate) is aproximately 2s.
-{{:tailorcrete:examples:v_funnel_tau60_mu20.avi|}}
+{{video>http://www.oofem.org/wiki/lib/exe/fetch.php?media=tailorcrete:examples:v_funnel_tau60_mu20.avi|Animation of V-funnel emptying.}}
+==== References ====
+  - BARTH, T.; SETHIAN, J.A. (2009), Numerical Schemes for the HamiltonJacobi and Level Set Equations on Triangulated Domains. Journal of computational physics, 145 1-40.
+  - TEZDUYAR, T : Stabilized Finite Element Formulations for Incompressible Flow Computations, Advances in Applied Mechanics, Volume 28, 1991, Pages 1-44