CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM

channelFoam for LES of Flow over a Wavy Wall

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 6, 2012, 01:50
Angry channelFoam for LES of Flow over a Wavy Wall
  #1
Member
 
Nick
Join Date: Nov 2011
Location: Tongji University,Shanghai,China
Posts: 33
Blog Entries: 6
Rep Power: 5
sunliming is on a distinguished road
Hi Foamers,
I adopted channelFoam to simulate flow over a wavy wall which is sinusoidal shaped,
seeing that merely the geometrical configurations are different from tutorial channel395.

----------------------------------------------------------------------------------------
Sizes of the domain are:
(Lx,Ly,Lz)=(4Lambda,Lambda,2Lambda),whereLambda=1 is wave length;
X:streamwise;
Y:vertical or channel depth direction;
Z:spanwise.
Refer to Case 76 of http://cfd.mace.manchester.ac.uk/ercoftac/ for details!

----------------------------------------------------------------------------------------
Grids:192 X 96 X 96;
Vertcial expansion ratio 1.04 with y+ equals about 0.8;
Bulk velocity Um=1m/s;
Re=Um*Lambda/nu=6760;

----------------------------------------------------------------------------------------
Boundary Conditions:
All are the same as channel395 but the topWall,here I removed the topWall exposing
the channel to the atmosphere,making it an open channel .Hence I set BC of slip for both
U and p at the topWall,i.e.,
topWall patch;
bottomWall: Wall;
inlet,outlet,frontWall backWall: cyclic;

----------------------------------------------------------------------------------------
Sub-grid model:homogenousdynSmagorinsky

----------------------------------------------------------------------------------------
Unluckily,I encountered several problems:

*First, after about 4 flow-times through the channel I got some results,and I'd like to get
a field average in spanwise direction Z ONLY,but the postChannel application exerts averaging
throuth all over the field including streamwise direction which is undesired in my simulation,
cause my homogenous direction is ONLY Z! So how can I revise the postChannel code and recompile
to meet my command? Or are there already some existing codes that can help?

*Second,after about 7 flow-times,the calculation exploded due to the gradually growing Courant
number which started at a mean value of 0.2(maximum 0.6) and ended up to an amazing value of several
thousands!Meanwhile, the pressure gradient ‘gradP’is dropping all the way from 0.00886 to below
0.0001,yielding friction velocity u_tau dropped form 0.09 to 0.01!
I'v no idea how this happened!So I wanna know how is the gradP calculated in the straight
channel,is it a mean value of the whole field which no longer suit my wavy case?Can I still use
this channleFoam application for my simulation od should I change for another application e.g.
pisoFoam?

----------------------------------------------------------------------------------------
Any advise is appreciated!
Best Wishes!
sunliming
Tongji University
Shanghai, China
Email:sunliming59045@hotmail.com
sunliming is offline   Reply With Quote

Old   March 6, 2012, 02:07
Default
  #2
Member
 
Nick
Join Date: Nov 2011
Location: Tongji University,Shanghai,China
Posts: 33
Blog Entries: 6
Rep Power: 5
sunliming is on a distinguished road
Attached is the geometrical configuration of the open channel copied from Hudson's paper,but the coordinate system should be adjusted.
Attached Images
File Type: jpg 123.jpg (56.5 KB, 48 views)
sunliming is offline   Reply With Quote

Old   June 7, 2012, 00:35
Default moving wavy channel
  #3
Member
 
Nick
Join Date: Nov 2011
Location: Tongji University,Shanghai,China
Posts: 33
Blog Entries: 6
Rep Power: 5
sunliming is on a distinguished road
moving wavy channel
Attached Images
File Type: jpg moving boundary.jpg (39.6 KB, 236 views)
sunliming is offline   Reply With Quote

Old   June 7, 2012, 00:57
Default moving wavy channel
  #4
Member
 
Nick
Join Date: Nov 2011
Location: Tongji University,Shanghai,China
Posts: 33
Blog Entries: 6
Rep Power: 5
sunliming is on a distinguished road
Hello all,
I am doing a LES simulation of flow over a flexible wall undergoing a streamwise travelling wave motion, the lower boundary of the domain is shown in the picture below. The elevation of the lower wall is a function of both time and space, i.e. y_wall=amplitude*sin[k(x-ct)], where c is the phase speed and k is frequency, which means every point of the lower boundary oscillates up and down its own equilibrium position with different phase angle.
So how can I deal with such kind of boundary motion and corresponding dynamic mesh? Can I use codeFixedValue BC? How?
Thank you in advance!
Beat regards!
Sunliming
Tongji University,
Shanghai, China,
Email: sunliming59045@hotmail.com
sunliming is offline   Reply With Quote

Old   May 1, 2013, 16:59
Default
  #5
Member
 
pooyan
Join Date: Nov 2011
Posts: 61
Rep Power: 5
sam1364 is on a distinguished road
How did you fix Re for this flow geometry?


Quote:
Originally Posted by sunliming View Post
Hello all,
I am doing a LES simulation of flow over a flexible wall undergoing a streamwise travelling wave motion, the lower boundary of the domain is shown in the picture below. The elevation of the lower wall is a function of both time and space, i.e. y_wall=amplitude*sin[k(x-ct)], where c is the phase speed and k is frequency, which means every point of the lower boundary oscillates up and down its own equilibrium position with different phase angle.
So how can I deal with such kind of boundary motion and corresponding dynamic mesh? Can I use codeFixedValue BC? How?
Thank you in advance!
Beat regards!
Sunliming
Tongji University,
Shanghai, China,
Email: sunliming59045@hotmail.com
sam1364 is offline   Reply With Quote

Reply

Tags
channel flow, courant number increasing, postchannel, pressure gradient, wavy boundary

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
Slip flow boundary condition on wall using UDFs mahdi haddad FLUENT 1 November 11, 2014 07:13
3 phase flow (wall boundary condition) Mohsin FLUENT 1 October 6, 2011 08:22
problems for the flow over flat wall iilw1314 Main CFD Forum 0 March 19, 2010 06:53
Influece of wall velocity in the main flow marvin CFX 0 March 22, 2008 04:05
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19


All times are GMT -4. The time now is 10:22.