|
[Sponsors] |
Please help with implemetation of boundaries conditions for flow in a pipe |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Member
Michail
Join Date: Apr 2009
Location: Lithuania
Posts: 41
Rep Power: 17 ![]() |
Hello all. Please help with implementation of boundary conditions for fully developed flow in a pipe. I use SIMPLE-algorithm on collocated grid. The enumerating of CV is below. I implemented dp/dx = constant and dU/dx=0 as:
------------------------- For inlet : P(1,J)= P(2,J)-( P(3,J) - P(2,J) )*0.5 U(1,J) = U_inlet ------------------------- For outlet: P(NXmaxC,J)= P(NXmaxC,J)-( P(NXmaxC-3,J) - P(NXmaxC-2,J) )*0.5 U(NXmaxC,j)=U(NXmaxC-1,j) ![]() But the solution diverged. May be I am wrong? |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Emre
Join Date: Oct 2009
Location: Ann Arbor, MI
Posts: 12
Rep Power: 16 ![]() |
dp/dx=0 is not a good idea for pipe flow. Fix the pressure at the outlet and velocity or mass flow rate (or pressure) at the inlet.
------------------------- For inlet : U(1,J) = U_inlet ------------------------- For outlet: P(NXmaxC,J)= P_outlet |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Michail
Join Date: Apr 2009
Location: Lithuania
Posts: 41
Rep Power: 17 ![]() |
Thank You very much for Your suggestion. I'll try immediately and shall answer you today.
But what to do with pressure at inlet? |
|
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
Emre
Join Date: Oct 2009
Location: Ann Arbor, MI
Posts: 12
Rep Power: 16 ![]() |
I suggest linearly extrapolating the inlet pressure from inside
|
|
![]() |
![]() |
![]() |
![]() |
#5 | |
New Member
sadegh
Join Date: Jul 2011
Posts: 2
Rep Power: 0 ![]() |
Quote:
i have problem i want a solve turbulence flows of air on a flat plate in two dimensional with matlab but i dont know what starting from where i request help me,can u help me? my method is control volume and my algorithm is simplec. |
||
![]() |
![]() |
![]() |
![]() |
#6 |
Member
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 30
Rep Power: 14 ![]() |
Hello everyone, I am trying to solve for steady Poiseulle flow on a backward staggered grid using FVM and SIMPLE, these are the initial boundary conditions that I set...
%% Wall 1,The inlet u(2, ![]() v(1, ![]() p(1:2, ![]() %% Wall 2, The top plate %For Poiseulle flow; No slip boundary condition: u(:,nYp)=0; v(:,nY)=0; p(:,nYp)=0; %% Wall 3, Outlet u(nX,2:nYp-1)=u(nX-1,2:nYp-1); % First extrapolate values from the domain u(nX,2:nYp-1)=u(nX-1,2:nYp-1)*sum(u(2,2:nYp-1))/sum(u(nX,2:nYp-1)); % This multiplier is the mass in divided by the mass out, which ensures continuity. v(nXp, ![]() ![]() p(nXp,2:nYp-1)=0; % Constant static pressure value of '0' set at the outlet (I = NI) %% Wall 4, Plate %For Poiseulle flow; No slip boundary condition: u(:,1)=0; v(:,2)=0; p(:,1)=0; From here, I solve for the momentum equations, noting that the source term is the pressure gradient, as well as the source terms in appendix C of Versteeg. The wall shear stress also enters the u-momentum equations in the central coefficient, Versteeg: Sp = -(mu/delta yp)*Acell ... eqn 9.11 pg 275. Since I am using a backward staggered grid, there is no need to incorporate the inlet and outlet boundary conditions as source terms (right hand side of equation), since the hybrid differencing scheme ustilizes the already stored inlet value of uin at i=2, and the u-momentum equations are calculated from i=3. Once u* and v* are calculated, solve for the pressure correction equaiton, being sure to cut the coefficients to the inlet (aW), outlet (aE), and the top and bottom plates (aN and aS), noting also that the coefficients are simply the central coeficients of the u- & v- momentum equations, multiplied by the density and the cellface area (rho*d*a)i,J etc... and that the source term is the continuity imbalance accross a pressure volume, i.e West face - East face + south face -north face. I am using the MATLAB function BICGSTAB to solve the sets of equations... can anyone point out something wrong? REgards, 0565113@gmail.com |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
pipe with buoyantFoam buoyancy, boundary conditions | Thomas Baumann | OpenFOAM | 0 | June 15, 2009 08:58 |
Boundary Conditions for Stratified Pipe Flow | Zhang | FLUENT | 0 | June 30, 2008 13:33 |
Terrible Mistake In Fluid Dynamics History | Abhi | Main CFD Forum | 12 | July 8, 2002 09:11 |
Pipe Flow Boundary conditions | Kunal Jain | Main CFD Forum | 0 | September 4, 2001 20:52 |
fluid flow fundas | ram | Main CFD Forum | 5 | June 17, 2000 21:31 |