
[Sponsors] 
November 11, 2005, 21:48 
question on bounday layer modeling

#1 
Guest
Posts: n/a

Hi, friends:
I'd like to ask a simple question on bounday layer simulation. The problem is : given free stream velocity U(t)=A*cos(wt), where w is angular frequency , how to solve a linear laminar boundary layer numerically (my problem is rough wall turbulence, but I simplify it to laminar case here for simplicity). Seems that it's a quite simple homework problem. We have governing equation: Equation(1): d(u)/dt = nu* d^(u)/dz^2 look for solution u(z,t), where z is distance from wall and nu is kinematic viscosity. Boundary conditions are : Equation(2): a)u=0 on z =0; Equation(3): b)u=U on z=h where h is is the thickness of the boundary or Equation(4): c)tau=nu*d(u)/dz =0 on z=h. Normally h(t) is not a constat, if we choose h large enough, we should be able to set h=constant. My question is really: how can we satisfy both (3) and (4)?? Seems that if we choose to use (3), then we won't be able to satisfy (4), or if we want to satisfy (4), we won't be able to satisfy (3). Because this is a 2 point BVP, only two boundary conditions are needed, and (2) can't be neglected. My problem is that : on the top of the bounday layer, if we want shear stress tau(h)=0, then we won't have solution u(h) =U; if we want u(h)=U, then tau(h)\=0. This is puzzling me for sometime, and I want to know your opinion. I also know that there is analytical solution to this problem (Stokes 1847, Lamb 1932) which says tau=0 at infinity, u=U at infinity. But I want to use finite difference for a far more complicated application, so I have to choose a finite domain 0<= z <=h. Now, I'm kind of biased to using (3) instead of (4), because it's important to keep u(h)=U(t), otherwise after some time of time integration, the solution u(t,z) will be out of whack of the free stream velocity U(t), but I don't know whether this is a standard practice for oscillatory boundary layer modeling. Best, Wen 

November 12, 2005, 05:50 
Re: question on bounday layer modeling

#2 
Guest
Posts: n/a

You don't try to satisfy both; e.g. If you take the full nonlinear problem
u_t + uu_x + vu_y = P_x + u_yy, u_x + v_y = 0, P_y = 0, (You can linearize these if you want it won't change what I'm saying) with the boundary condition u=v=0 on y=0. For the upper boundary condition, as y> infinity, you have two equivalent options: (1) set u>U(x,t) => P_x = U_t + UU_x (Bernoulli), (2) set P_x using Bernoulli and U_y > 0. Basically you choose one of these and, provided the lid is highing enough, the other condition is automatically satisfied. 

November 12, 2005, 18:29 
Re: question on bounday layer modeling

#3 
Guest
Posts: n/a

Tom:
Thank you for the tip. Wen 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
To thin boundary layer?  Luk  CFX  3  February 27, 2009 04:22 
Question, modeling airflow through radiator  Roland  CFX  6  August 7, 2006 07:25 
Question about layer additionremove  liu  OpenFOAM Running, Solving & CFD  6  February 21, 2006 17:26 
errors  Fahad  Main CFD Forum  0  March 23, 2004 14:20 
Boundary Layer Flow Paradox  Wen Long  Main CFD Forum  3  September 24, 2002 09:47 