
[Sponsors] 
July 27, 2003, 05:00 
fully developed velocity profile in 3D

#1 
Guest
Posts: n/a

I have a laminar, viscous 3D flow in a duct (brick type geometry) with rectangular crosssection. I want to use velocity inlet and pressure outlet boundary conditions. There is a nice expression for fullydeveloped 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.


July 27, 2003, 15:12 
BC at the inlet of the tube

#2 
Guest
Posts: n/a

Hi all,
I've had a problem of specifying the boundary condition at the inlet of an openedend 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 

July 27, 2003, 16:00 
Re: BC at the inlet of the tube

#3 
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).


July 27, 2003, 16:05 
Re: fully developed velocity profile in 3D

#4 
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} (by)}{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)^n1] $ and $\alpha =  \Delta p/(L\mu)$, $\Delta p$ is pressure drop in the $z$ direction with length $L$, and flow viscosity $\mu$. 

July 27, 2003, 16:47 
Re: fully developed velocity profile in 3D

#5 
Guest
Posts: n/a

Yes, but can it be expressed in terms of u_max (maximum velocity)? In 2D, 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 ycoordinates. 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. 

July 27, 2003, 23:07 
Re: fully developed velocity profile in 3D

#6 
Guest
Posts: n/a

Hi,
If u have a nice expression for fullydeveloped 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 3D, 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 2d velocity field to 3d 

July 29, 2003, 12:02 
Re: fully developed velocity profile in 3D

#7 
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.


July 6, 2013, 15:36 

#8 
New Member
Ugur Goktolga
Join Date: Jul 2013
Posts: 2
Rep Power: 0 
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 335338 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. 

January 10, 2017, 01:48 
fully developed profile at the inlet of a duct of rectangular section

#9  
Member
vidyadhar
Join Date: Jul 2016
Posts: 68
Rep Power: 3 
Quote:
I have seen this thread more than a decade after you have posted it. I want to have fully developed velocity profile at the inlet of a duct with rectangular cross section. I request you to help me in executing the same. Thank you, vidyadhar 

January 10, 2017, 02:23 

#10  
Senior Member
Lane Carasik
Join Date: Aug 2014
Posts: 588
Rep Power: 7 
Quote:


January 11, 2017, 13:39 
Equation used to generate fully developed profile in the rectangular duct

#11  
Member
vidyadhar
Join Date: Jul 2016
Posts: 68
Rep Power: 3 
Quote:
Thank you for the quick reply. I have seen the code provided in the link. May I request you to please let me know the equation used in the code.Also, what is um in the equation given. Though the reference mentioned in the code is a book by Shah and London, I am unable to access it. Thanks & Regards, Vidyadhar 

January 11, 2017, 13:44 

#12  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,492
Rep Power: 40 
Quote:
the solution for a laminar steady developed flow in a rectangular duct satisfies the Poisson equation Lap u = Re*dp/dx with u=0 on the boundaries. 

January 11, 2017, 13:49 

#13  
Member
vidyadhar
Join Date: Jul 2016
Posts: 68
Rep Power: 3 
Quote:
Thank you for the reply. Since there are many analytical expressions available in various sources for the Poisson equation, I would like to know which form was used by the coder to develop the code. I request you to please let me know the same. Thanks & Regards, vidyadhar 

January 11, 2017, 14:00 

#14 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,492
Rep Power: 40 

January 11, 2017, 14:06 

#15  
Member
vidyadhar
Join Date: Jul 2016
Posts: 68
Rep Power: 3 
Quote:
I will try to find the equation in the book, if it is not hidden. Regards, vidyadhar 

Thread Tools  
Display Modes  


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 Starcd  SMM  STARCD  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 