CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Open channel flow with interFoam

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

Like Tree1Likes
  • 1 Post By sorockin

Reply
 
LinkBack Thread Tools Display Modes
Old   August 2, 2013, 11:20
Default Open channel flow with interFoam
  #1
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 4
sorockin is on a distinguished road
Dear Foamers, I need some help.

I try to simulate flow in open channel with interFoam. I use setting from wigleyHull tutorial. I only changed the geometry of channel and inlet velocity (0.25 m/s). But I see instability in simulated results. Tel me please what need to change to get flow without instabilities. Screenshots of alpha, U and case setup are attached.

Thanks a lot!
Maxim
Attached Images
File Type: jpg alpha1.jpg (17.9 KB, 78 views)
File Type: jpg U.jpg (27.4 KB, 73 views)
Attached Files
File Type: gz simpleChannel.tar.gz (3.1 KB, 8 views)
sorockin is offline   Reply With Quote

Old   August 3, 2013, 11:36
Default New settings
  #2
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 4
sorockin is on a distinguished road
I have moved form LTSInterFoam to interFoam and changed turbulence values setting. I think that problem with stability was in time constraints or wrong turbulent settings.
Now I use interFoam with maxCo=0.5 and maxAlphaCo=0.5.
I set internalField and inlet value for k from expression k=1.5*(U*I)^2. U=0.25, turbulence intensity I put 1%. So k=9.375e-6.
For mega - from expression omega=Cm^(-0.25)*k^0.5/L, Cm=0.09, turbulence length I put 7% of channel depth. So omega=5.324e-3.
I deleted nut file in 0/ folder to OpenFOAM calculate nut values.
But I still have flow instability.
Help me please I am new in OpenFOAM and don't know what needs to do!
New case setting are attached.

Maxim
Attached Images
File Type: jpg alpha1.jpg (14.4 KB, 30 views)
File Type: jpg Ux.jpg (18.5 KB, 33 views)
Attached Files
File Type: gz simpleChannel.tar.gz (3.8 KB, 10 views)
JFM likes this.
sorockin is offline   Reply With Quote

Old   August 5, 2013, 08:35
Default laminar
  #3
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 4
sorockin is on a distinguished road
If I turn off turbulence and simulate laminar flow I still obtain instability. What is the cause? Solution scheme or boundary condition? But I took boundary conditions from wigleyHull tutorial. I am confused.

Maxim
Attached Images
File Type: jpg Ux.jpg (20.2 KB, 36 views)
File Type: jpg alpha1.jpg (14.7 KB, 31 views)
sorockin is offline   Reply With Quote

Old   August 6, 2013, 08:34
Default
  #4
Member
 
David Gisen
Join Date: Jul 2009
Location: Germany
Posts: 36
Rep Power: 7
David* is on a distinguished road
The second case doesn't look bad in my opinion. You set turbulence at the inlet, and you get turbulence.

In the third case the reason might be some BC.

If you say instability one thinks of numerical instability. You can check it via the residuals. In the forum there is a script for gnuplot that processes the values given in your .log-file (if you write a .log-file).
David* is offline   Reply With Quote

Old   August 6, 2013, 09:41
Default
  #5
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 4
sorockin is on a distinguished road
Thank you for reply David!

The second case doesn't look bad in my opinion. You set turbulence at the inlet, and you get turbulence.

I used initial setting from "wigleyHull" tutorial. "wigleyHull" used turbulence model too but flow is much more uniform. I simulated this case with flow3D (k-epsilon turbulence model, the same mesh) and had uniform flow.

In the third case the reason might be some BC.

How nice it would be to know exactly what it is - BC or solution scheme.

If you say instability one thinks of numerical instability. You can check it via the residuals. In the forum there is a script for gnuplot that processes the values given in your .log-file (if you write a .log-file).

It is good idea, thank you!

Maxim
sorockin is offline   Reply With Quote

Old   August 6, 2013, 10:30
Default laminar, changed setup
  #6
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 4
sorockin is on a distinguished road
I changed setup for laminar simulation. I moved from 3D to 2D (xz plane) and divided inlet plane. Onto inlet and wall.

So now I have channel 85x1x17 m with mesh 340x1x170. Water depth is 15m.
With BCs:

inlet
U: fixedValue; value uniform (0.25 0 0);
alpha1: fixedValue; value uniform 1;
p_rgh: zeroGradient;

inletWall
U: symmetryPlane;
alpha1: symmetryPlane;
p_rgh: symmetryPlane;

outlet
U: zeroGradient;
alpha1: zeroGradient;
p_rgh: zeroGradient;

atmosphere
U: pressureInletOutletVelocity; phi phi; value uniform (0 0 0);
alpha1: inletOutlet; inletValue uniform 0; value uniform 0;
p_rgh: totalPressure;

bottom
U: symmetryPlane;
alpha1: symmetryPlane;
p_rgh: symmetryPlane;

I still have "instability" at inlet. It may be connected with wrong inlet BC. I will be glad any advice. The case is attached.

Maxim
Attached Images
File Type: jpg alpha1.jpg (14.0 KB, 39 views)
File Type: jpg p.jpg (18.5 KB, 44 views)
File Type: jpg Ux.jpg (13.2 KB, 37 views)
Attached Files
File Type: gz simpleChannel.tar.gz (3.6 KB, 5 views)
sorockin is offline   Reply With Quote

Old   August 7, 2013, 11:24
Default Mesh, BC or schems
  #7
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 4
sorockin is on a distinguished road
I made another changes.
My mesh had dx/dz ratio 2.5. I made dx=dz to remove mesh as one of cause that may induce instability.

Then I changed BCs, accordingly Matthew Denno
interFoam simulation blowing up

inlet
U: fixedValue; value uniform (0.25 0 0);
alpha1: fixedValue; value uniform 1;
p_rgh: buoyantPressure; value uniform 0;

inletWall
U: fixedValue; value uniform (0 0 0);
alpha1: zeroGradient;
p_rgh: buoyantPressure; value uniform 0;

outlet
U: inletOutlet; inletValue uniform (0 0 0); value uniform (0 0 0);
alpha1: inletOutlet; inletValue uniform 0; value uniform 0;
p_rgh: inletOutlet; inletValue uniform 0; value uniform 0;

atmosphere
U: pressureInletOutletVelocity; value uniform (0 0 0);
alpha1: inletOutlet; inletValue uniform 0; value uniform 0;
p_rgh: totalPressure;

bottom
U: fixedValue; value uniform (0 0 0);
alpha1: zeroGradient;
p_rgh: buoyantPressure; value uniform 0;

And I still have "instability" that arise at inlet and then result strongly nonuniform flow.
I tried different "working" BC set from another threads but "instability" didn't disappeared.
May be deal in solution schemes. Can anybody help me? Please!

Maxim
Attached Images
File Type: jpg Ux_6s.jpg (16.2 KB, 32 views)
File Type: jpg Ux_14.5s.jpg (26.8 KB, 45 views)
Attached Files
File Type: gz simpleChannel.tar.gz (3.5 KB, 13 views)
sorockin is offline   Reply With Quote

Old   August 9, 2013, 10:14
Default fvSolution
  #8
New Member
 
Maxim Sorockin
Join Date: Oct 2012
Posts: 9
Rep Power: 4
sorockin is on a distinguished road
I took fvSolution setting from Matthew case and looks like "instability" disappeared.
Thank you Matthew for helpful post!

Maxim
Attached Images
File Type: jpg Ux_23s.jpg (14.7 KB, 61 views)
sorockin is offline   Reply With Quote

Old   October 20, 2013, 06:50
Default Mattew Case
  #9
JFM
New Member
 
JFM's Avatar
 
John Anonymous
Join Date: Jan 2011
Location: Melbourne Australia
Posts: 18
Rep Power: 5
JFM is on a distinguished road
Hi Maxim

You refer to Matthew's case in the post above, can you provide a link to this.

Thanks
JFM
JFM is offline   Reply With Quote

Reply

Tags
channel flow, instability, interfoam

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
No fluid in open channel flow VOF flow_CH Fluent Multiphase 2 July 2, 2013 01:05
interFoam vs. simpleFoam channel flow comparison DanM OpenFOAM Running, Solving & CFD 11 January 5, 2013 06:21
OpenFOAM open channel flow downsteam boundary condition Brickman OpenFOAM Running, Solving & CFD 2 November 5, 2012 20:14
Boundary condition problem for open channel flow Andy CFX 8 October 1, 2012 07:45
Open Channel Flow forsumit FLUENT 0 October 1, 2009 02:01


All times are GMT -4. The time now is 12:55.