CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Phoenics (http://www.cfd-online.com/Forums/phoenics/)
-   -   Ground Coding Using PLANT (http://www.cfd-online.com/Forums/phoenics/52315-ground-coding-using-plant.html)

Alex January 23, 2009 06:34

Ground Coding Using PLANT
 
Dear friends

I have currently considered a simple 2D flow in a duct (X-Z coordinate system, CARTES=T), and manage to get a converged solution. However I tried to compute the overall momentum balance using PLANT eventhough the result file have provide the information. I have manage to obtain the shear stress contibution in (Pa/area).

The problem is that to get the distance between the cell centers in x-direction is unrealistic when using either XG2D or DXG2D (when printed out it give the following values). I have the correct result using ZGNZ, I would be appreciate for any suggestion

cell IX=1 -10101.00 cell IX=2 -10101.00 -10101.00 -10101.00 -10101.00 -10101.00 -10101.00 -10101.00 -10101.00

1.000000

3.000000

7.000000

145.0000

146.0000

147.0000

148.0000

149.0000

150.0000

3798.000 cell IX=20 3818.000


Alex January 23, 2009 08:52

Re: Ground Coding Using PLANT
 
Dear friend

Thanks again for the silence, as this make me to strunggle. I discover that as satlit is processed the ground coding, it defined the DXG2D as LFDXG2=ABS(ANYZ(DXG2D ,IZZ)) as shown below. I have to change this to LFDXG2=L0F(DXG2D ) to make its effect. Thanks for motivating me.

IF(ISTEP.GE.1 .AND.ISTEP.LE.LSTEP )THEN

DO 19701 IZZ=1 ,NZ

LFWTAU =ABS(ANYZ(INAME('WTAU '),IZZ))

LFDXG2=ABS(ANYZ(DXG2D ,IZZ))

DO 19701 IX=1 ,NX

IADD=NY*(IX-1)

DO 19701 IY=1 ,NY

I=IY+IADD

L0WTAU=LFWTAU+I

L0DXG2=LFDXG2+I 19701 F(L0WTAU )=F(L0DXG2)

ENDIF


All times are GMT -4. The time now is 05:35.