|
[Sponsors] |
Flow Between Parallel Plates With Parabolic Entrance Velocity |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 4, 2014, 04:30 |
Flow Between Parallel Plates With Parabolic Entrance Velocity
|
#1 |
New Member
saeb
Join Date: Mar 2013
Posts: 15
Rep Power: 13 |
Dear Friends,
Previously I worked on the steady flow between parallel plates with uniform entrance velocity and I got satisfactory results. Since then I have been working on parabolic entrance velocity but I am not sure about the validity of the results. The following plot is the centerline axial velocity when the mean axial velocity is Um=0.1004708 (which seems alright): but when I decrease the range of the vertical axis it seems that there is an error which I am not sure if this is an acceptable numerical error or not. The following is the same plot with smaller range of vertical axis: As you see this error only shows itself in very small range of vertical axis. Last edited by saaghi; February 4, 2014 at 12:02. |
|
February 7, 2014, 04:06 |
|
#2 |
Member
Ardalan
Join Date: Feb 2010
Location: Baku
Posts: 50
Rep Power: 16 |
what are your outlet boundary conditions?
It seems to be numerical errors. Parabolic velocity profile is the exact solution to this problem. but your solution is numerical and have some errors depending on your order of dicretization. So I suggest that you refine mesh in y direction. If the error decreases as you refine the mesh, be sure that it is just discretization error. |
|
February 7, 2014, 04:48 |
|
#3 |
New Member
saeb
Join Date: Mar 2013
Posts: 15
Rep Power: 13 |
At the outlet we have fully developed boundary condition (zero gradient for velocity components).
Do you suggest to use a finer mesh in Y direction or axial direction(X)? I don't think it is from the discretization scheme, because I checked the results with upwind, QUICK, and CUBISTA schemes and got almost the same results for all. |
|
February 7, 2014, 05:04 |
|
#4 |
Member
Ardalan
Join Date: Feb 2010
Location: Baku
Posts: 50
Rep Power: 16 |
I think your flow is not fully developed at the outlet. increase the length of channel.
if you are using simple algorithm, check the global mass conservation. your mass flow rate is not equal to the integral of inlet parabolic profile. |
|
February 7, 2014, 05:23 |
|
#5 |
New Member
saeb
Join Date: Mar 2013
Posts: 15
Rep Power: 13 |
Since the inlet velocity profile is parabolic, the channel length doesn't matter anymore, the flow is fully developed all along the length.
I am using the simple algorithm and I think the problem is in the mass conservation too. Are you saying I have to check the corrections? |
|
February 7, 2014, 06:38 |
|
#6 |
Member
Ardalan
Join Date: Feb 2010
Location: Baku
Posts: 50
Rep Power: 16 |
Implement boundary condition as follows:
at every time step calculate inlet and outlet mass flow rate (M_in , M_out) then at the outlet use U(imax)=U(imax-1)*M_in/M_out instead of zero boundary condition for pressure correction use zero gradient at the outlet This will solve your problem |
|
February 7, 2014, 07:09 |
|
#7 |
New Member
saeb
Join Date: Mar 2013
Posts: 15
Rep Power: 13 |
I have already done all these, but the problem should be somewhere else.
Thank you anyway. |
|
February 7, 2014, 07:19 |
|
#8 |
Member
Ardalan
Join Date: Feb 2010
Location: Baku
Posts: 50
Rep Power: 16 |
in this case M_in is not equal to U_m*h
M_in=SUM(u_1j*Dy) Have you considered this? |
|
February 7, 2014, 08:06 |
|
#9 |
New Member
saeb
Join Date: Mar 2013
Posts: 15
Rep Power: 13 |
Yes, of course I did.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
How do I prescribe Average Velocity /Total Inlet Flow for a simple channel problem? | skuznet | OpenFOAM Pre-Processing | 4 | February 16, 2022 09:44 |
Laminar steady flow between parallel plates using Fortran | saaghi | Main CFD Forum | 10 | July 2, 2020 09:23 |
Emergency:UDF for a time dependent parabolic velocity | zumaqiong | Fluent UDF and Scheme Programming | 12 | March 25, 2010 12:00 |
Influece of wall velocity in the main flow | marvin | CFX | 0 | March 22, 2008 02:05 |
Variables Definition in CFX Solver 5.6 | R P | CFX | 2 | October 26, 2004 02:13 |