CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

fully developed velocity profile in 3D

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 27, 2003, 05:00
Default fully developed velocity profile in 3D
  #1
S.S.
Guest
 
Posts: n/a
I have a laminar, viscous 3D flow in a duct (brick type geometry) with rectangular cross-section. I want to use velocity inlet and pressure outlet boundary conditions. There is a nice expression for fully-developed velocity profile (parabola) in 2D (written in terms of u_max). How is it written in 3D? Any help or reference will be highly appreciated. Thanks in anticipation.
  Reply With Quote

Old   July 27, 2003, 15:12
Default BC at the inlet of the tube
  #2
Jackie
Guest
 
Posts: n/a
Hi all,

I've had a problem of specifying the boundary condition at the inlet of an opened-end tube. I use CFD to simulate flow inside the tube having an embedded momentum source. The momentum source is located at the center of the horizontal tube and both ends (inlet and outlet) of the tube is opened to the atmosphere.

I specify zero static pressure at the outlet (reference to ambient). The flow inside the tube is driven by the given momentum source. However, I don't know what the boundary condition at the inlet of the tube should be.

Could anybody please give me some suggestion?

Thank you very much

Jackie

  Reply With Quote

Old   July 27, 2003, 16:00
Default Re: BC at the inlet of the tube
  #3
xueying
Guest
 
Posts: n/a
B.C. at inlet might be fully developed, i.e. \vec{n} \cdot \nabla \vec{v} = \vec{0} (in latex form).
  Reply With Quote

Old   July 27, 2003, 16:05
Default Re: fully developed velocity profile in 3D
  #4
xueying
Guest
 
Posts: n/a
The fully developed flow velocity profile for a rectangular channel is a serial solution; following is the analytical solution written in latex form:

\[

u(x,y) = \sum_{n=1}^{\infty} u_n(y) \phi_n(x) \] where \[

u_n(y) = \frac{a_n}{\lambda_n} \left( -1+\frac{sinh\sqrt{\lambda_n}y+sinh\sqrt{\lambda_n} (b-y)}{sinh\sqrt{\lambda_n}b}\right) \]

\[

\phi_n(x) = \sqrt{\frac{2}{a}}sin(\frac{n \pi x}{a}), \;\; n=1,2,... \]

where $a$ is the length in the $x$ direction, $b$ is the length in the $y$ direction, $\lambda_n = \frac{n^2\pi^2}{a^2}$, $a_n = \sqrt{\frac{2}{a}}\alpha \frac{a}{n\pi}[(-1)^n-1] $ and $\alpha = - \Delta p/(L\mu)$, $\Delta p$ is pressure drop in the $z$ direction with length $L$, and flow viscosity $\mu$.
  Reply With Quote

Old   July 27, 2003, 16:47
Default Re: fully developed velocity profile in 3D
  #5
S.S.
Guest
 
Posts: n/a
Yes, but can it be expressed in terms of u_max (maximum velocity)? In 2-D, there is a nice relation for u(y) which is given in terms of u_max (maximum velocity which occurs at the centerline of the parabola), and the y-coordinates. This expression is very convenient to use as velocity inlet b.c. profile in CFD codes (provided that u_max is known). Is there anything like that in 3D?

Also, can someone give me a reference where the detailed derivation of this exact 3D solution can be found?

Thanks.
  Reply With Quote

Old   July 27, 2003, 23:07
Default Re: fully developed velocity profile in 3D
  #6
Hall
Guest
 
Posts: n/a
Hi,

If u have a nice expression for fully-developed velocity profile (parabola) in 2D, I think it is easy to expand it to 3d, just let the profile rotate around the longitude axis

let's say if x is the longitude axis, and at the inlet, u = U(y) is a function of y. Now at 3-D, again you set the x' as longitude axis. At inlet, you will have coordinates like (y',z'), take r = sqrt(y'^2 + z'^2), and let u = U(r). Now, u expand the 2-d velocity field to 3-d
  Reply With Quote

Old   July 29, 2003, 12:02
Default Re: fully developed velocity profile in 3D
  #7
xueying
Guest
 
Posts: n/a
The detailed derivation of this exact 3D solution can be found in book titled as Mathematical Methods in Chemical Engineering; the authors are: Varma, A. & Morbidelli, M.; Oxford University Press, New York Oxford, page 494.
  Reply With Quote

Old   July 6, 2013, 15:36
Default
  #8
New Member
 
Ugur Goktolga
Join Date: Jul 2013
Posts: 2
Rep Power: 0
mugurg is on a distinguished road
I know it has been decade since this thread was opened, but just in case somebody googles and finds this page, I wanted to help.

You can check the equations 335-338 from Shah and London. Generating parabolic velocity profile for 3D flows might not be a good idea since it does not represent fully developed flow for especially high aspect ratio channels. So, better check the equations that I have written. Actually, I have written a MATLAB script that generates fully developed velocity profile for rectangular channels and uploaded it to MATLAB file exchange forums, but it is not approved yet. If somebody replies this thread, I can supply the link for the file, or send as an attachment.
mugurg is offline   Reply With Quote

Reply

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
time dependent BC with fully developed velocity profile susheil FLUENT 6 December 20, 2012 06:41
Fully Developed Flow in Star-cd SMM STAR-CD 0 September 5, 2011 22:08
Fully developed flow profile icfd Main CFD Forum 3 May 20, 2009 22:21
fully developed velocity profile KS FLUENT 1 November 12, 2008 09:36
inlet velocity profile for fully developed turbule pankaj saha Main CFD Forum 4 August 24, 2007 16:53


All times are GMT -4. The time now is 13:34.