11 - Simulation of Flow Through a Pipe Part 2 2

This document represents the second part of the project related to the simulation of an incompressible laminar flow through a pipe. In the previous part of the project, a 4 edge degrees simulation was carried out using icoFoam and wedge boundary conditions. The purpose of this part is to employ the symmetry boundary conditions, applied at a geometry with 3 different angles (10,25 and 45) and compare both results. As the transient part is no longer of interest, and for further comparison, the solver simpleFoam was used in this part.

First, all the new meshes were generated with the Matlab program of the first part of the project, and the boundary conditions were changed accordingly. In the following picture, the mesh for the 45 degrees case (symmetry BC) is shown as an example. The results of the computational simulations will be discussed next.



Comparison of the velocity profiles (Fully developed vs Analytical)



Comparison of the Pressure Profile through the Pipe


Average Velocity [m/s]:

Hagen-Poiseuille 4 Degrees 10 Degrees 25 Degrees 45 Degrees
 0.10731 0.1110   0.1096  0.1082  0.1076

Maximum Velocity [m/s]:

Hagen-Poiseuille 4 Degrees 10 Degrees 25 Degrees 45 Degrees
 0.2146 0.2210   0.2188  0.2166  0.2148

 Pressure Drop [Pa]:

Hagen-Poiseuille 4 Degrees 10 Degrees 25 Degrees 45 Degrees
43.4 37.45   37.52 37.62  38.85


The results above show that there is not an important difference between the wedge boundary conditions and the symmetry boundary conditions. The small differences are most likely due to the angle simulated. In all of them, profile is fully developed, as the velocity profile is pretty similar to the analytic solution. The first thing to notice is that, the higher the angle simulated, the higher the fidelity of the results (higher values of mean and maximum velocity, closer to the analytical). Also notice that, as expected, the biggest differences appear near the wall and near the axis. 

In term of the pressure drop, all the simulations approximate to the analytical solution, but none of them reaches the analytical value close enough. The tendency shows that also the pressure drop increases when we increase the angle simulated. This means that the simulation cannot predict with exactitud the pressure drop in the tube, with the simple model used here. Further prove with other solvers, maybe icoFoam, should be carried out in order to choose the correct one.

Projects by Jorge Martinez

The purpose of this project is to show the effect of initial temperature and pressure in the ignition delay in a combustion process. For this purpose, a Constant Volume Reactor and the corresponding network objects are created in Cantera, and the combustion process is s Read more

The purpose of this project, is to setup and analyze a transient shock tube simulation using Converge Studio. The project will analyze how the interaction of the original shock waves and the reflected waves affect the fluid in the tube, and how the computational mesh ha Read more

The purpose of this project is to simulate the conjugate heat transfer in a pipe using Converge Studio. For this, an air flow is simulated through an aluminum pipe, which is receiveing a constant heat flux of 10,000 W/m2 at the outside wall. The objective is to understa Read more

This project objective is to simulate the Prandtl Meyer shock wave phenomena using Converge Studio. First, a quick literature review about shock waves and their boundary conditions is provided. Then , the problem is set up and solved, focusing in how the different param Read more

The purpose of this assignment is to create a Python program including a Multivariate Newton Rhapson Solver, to solve a non-linear coupled differential system. This method is really useful for stiff systems, where the explicit solver are unstable. The system of equatio Read more

The objective of this part of the project, is to simulate the transient flow through an elbow pipe with a throttle valve in the middle. This valve will be rotating from zero to 25 degrees, and then stay steady.  Case Setup. In the steady state part of the project Read more

This project purpose is to simulate the flow over an elbow body that contains a throttle using CONVERGE Studio. In this first part of the project, an steady state analysis will be set up. This means that the throttle will not move. This case helps set up the real analys Read more

In this challenge, the efficiency of using a heat recuperator to pre-heat the air is analyzed. For that, we simulate a furnace burning methane and air mixture, and we analyze first how the pre-heating affects at the AFL, and then the total heat/work that could be extrac Read more

This project main objective is to calculate the Adiabatic Flame Temperature using Python and Cantera. In the first part of the project, the effect of the equivalence ratio for the methane combustion in a constant volume chamber is analyzed, and the python and cantera re Read more

This projects aims to simulate a backward facing step flow using CONVERGE Studio. The flow is creating when a not reacting specie, air (23% O2 and 77% N2 in mass fraction), suffers a pressure gradient between the ends of the channel. The velocity, pressure and recirulat Read more


The End