|
[Sponsors] |
September 6, 2019, 11:34 |
Pressure drop - simpleFoam
|
#1 |
New Member
E. Lowry
Join Date: Sep 2019
Posts: 8
Rep Power: 6 |
I am currently working on matching the analytical solution for pressure drop in a pipe. I created a mesh of a pipe with 0.0254m diameter and 0.254m length with an inlet and outlet. The mesh was created using Ansys and refined in the near wall region to help with wall effects. The analytical solution indicates that the pressure drop with water flowing at 2 m/s should be about 1600 Pa/m. The simulation with simpleFoam using k-epsilon are giving an absurdly low value (usually less than 10 Pa/m). Things I have tried:
1. Ensured the viscosity value and units are correct 2. Calculated the correct k and epsilon values for the ICs 3. BCs are all okay, (no slip at boundaries, zero pressure outlet, etc.) 3. Corrected the kinematic pressure values to Pascal for final comparison 5. Ran checkMesh to ensure there are no irregularities in the mesh. Everything is OK according to checkMesh Suggestions on what other things I could try would be greatly appreciated. All run files except for the mesh (too big, ~2M cells) are attached for reference. Here is the last time in the log file with the residuals: Code:
Time = 362 smoothSolver: Solving for Ux, Initial residual = 0.00449619, Final residual = 0.000437133, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.00438297, Final residual = 0.000422318, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00501302, Final residual = 0.00046774, No Iterations 2 GAMG: Solving for p, Initial residual = 0.00122809, Final residual = 7.17839e-05, No Iterations 2 time step continuity errors : sum local = 0.000218834, global = -8.03457e-07, cumulative = -1.74794 smoothSolver: Solving for epsilon, Initial residual = 0.0133315, Final residual = 0.00125709, No Iterations 2 smoothSolver: Solving for k, Initial residual = 0.00878087, Final residual = 0.000864259, No Iterations 2 ExecutionTime = 497.16 s ClockTime = 498 s fieldValueDelta deltaTotalPressure_inOutlet write: subtract(areaIntegrate(inlet,p),areaIntegrate(outlet,p)) = 0.00114703 Many thanks in advance, E.Lowry |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Computed Pressure Drop is lower than experimental data | Ash Kot | FLUENT | 2 | May 17, 2017 09:41 |
lower pressure drop compared to experimental value in Two phase flow calculations. | Ash Kot | Fluent Multiphase | 0 | May 16, 2017 16:29 |
Pressure drop over throttle disappears (simpleFoam --> pipe simulation) | highpressuretube | OpenFOAM Running, Solving & CFD | 18 | December 21, 2016 12:07 |
How to plot pressure drop with Periodic BC? | bigfans | FLUENT | 7 | November 8, 2016 11:28 |
Pressure drop using Fan type BC | Alexis Sack | OpenFOAM Running, Solving & CFD | 2 | September 22, 2014 09:18 |