CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Verification & Validation

SimpleFoam: velocity profile does not match analytical solution

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   April 7, 2020, 01:58
Default SimpleFoam: velocity profile does not match analytical solution
New Member
Deutschland (DEU)
Join Date: Jan 2020
Posts: 3
Rep Power: 5
Abhinav_Nagarajan is on a distinguished road
Dear Foamers,
I am simulating a reactor with a weird geometry (please find picture attached), with a diffuser section connected to a rectangular duct.
Simulation conditions:
Inlet velocity: 0.55m/s
Reynold's number(at inlet): 1161

The issue: The centerline velocity decreases along the length of the duct, basically the centreline decelerates as it passes through the rectangular duct, whereas it should increase. Therefore, the results are not validated as per existing analytical models.

I am unable to figure where the problem lies, since I have a mesh with no errors. Could you please help me out with this issue. I have attached my fvschemes, fvsolution and the reactor snip along with the velocity vs length profile.

Thank you!
Attached Images
File Type: png reactor snip.png (59.5 KB, 43 views)
Attached Files
File Type: txt fvSchemes.txt (1.3 KB, 5 views)
File Type: txt fvSolution.txt (1.4 KB, 4 views)
File Type: txt checkMesh.txt (3.7 KB, 4 views)
File Type: txt simpleFoam.txt (179.0 KB, 9 views)
Abhinav_Nagarajan is offline   Reply With Quote

Old   May 12, 2022, 14:00
Join Date: Aug 2017
Posts: 61
Rep Power: 7
SHUBHAM9595 is on a distinguished road
Hi Abhinav,

Unfortunately, I've to disagree with you
I am unable to figure where the problem lies, since I have a mesh with no errors
As you can clearly see in the output of checkMesh, we do have some issues with the grid.


    Cell determinant (wellposedness) : minimum: 0.108244 average: 7.7354
    Cell determinant check OK.
 ***Concave cells (using face planes) found, number of cells: 32
  <<Writing 32 concave cells to set concaveCells
    Face interpolation weight : minimum: 0.114561 average: 0.498237
    Face interpolation weight check OK.
    Face volume ratio : minimum: 0.0876268 average: 0.992505
    Face volume ratio check OK.

Failed 1 mesh checks.

U can visualize this SPECIFIC set (concaveCells) in paraview with something likw ....
foamToVTK -cellSets concaveCells -time 0
and then get an idea where exactly u need to make changes within the domain...
SHUBHAM9595 is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
kindly help me .. i have and error at line number 147.. m zubair Fluent UDF and Scheme Programming 0 February 10, 2019 12:25
Defining Womersley or laminar velocity profile grkmgcl CFX 1 March 27, 2018 08:32
Fixed velocity profile: exploding solution CeesH FLUENT 1 August 31, 2017 11:36
Velocity profile disturbance due to loss coefficient rks171 Main CFD Forum 3 May 25, 2012 18:30
Wall functions Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16

All times are GMT -4. The time now is 07:47.