CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Implementation of Temperature Eqn in InterPhaseChangeFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 18, 2013, 08:22
Default Implementation of Temperature Eqn in InterPhaseChangeFoam
  #1
New Member
 
Jibran
Join Date: Oct 2012
Location: Swansea, UK
Posts: 25
Rep Power: 4
Jibran is on a distinguished road
Dear All,

Has anyone sucessfully implemented the temperature equation in interphaseChangeFoam?? I followed 'A interphasechangeFoam tutorial by Martin Andersen' but the results seem to be quite inaccurate. I have attached the relevant files below.

ALPHA1
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 0 0 0 0 0 0];
internalField   uniform 0;
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    
    right_wall
    {
        type            constantAlphaContactAngle;
        theta0          10;
        limit  gradient;
        value           uniform 0;
    }
    
    lower_wall
    {
        type            zeroGradient;
    }
    
    atmosphere
    {
        type  zeroGradient;    
    }
    
    default_faces
    {
        type            empty;
    }
}


P_RGH
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [1 -1 -2 0 0];
internalField   uniform 100000;
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    
    right_wall
    {
        type            fixedFluxPressure;
        adjoint         no;
    }
    
    lower_wall
    {
        type            buoyantPressure;
    }
    
    atmosphere
    { 
        type            fixedValue;
        value           uniform 100000;
    }
    
    default_faces
    {
        type            empty;
    }
}

T
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 0 0 1 0]; // [K]
internalField   uniform 280.15;
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    
    right_wall
    {
        type            fixedValue;
        value         uniform 300.15;
    }
    
    lower_wall
    {
        type            fixedValue;
        value         uniform 300.15;
    }
    
    atmosphere
    {
        type            fixedValue;   
        value         uniform 280.15;
    }
    default_faces
    {
        type            empty;
    }
}

U
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 1 -1 0 0];
internalField   uniform (0 0 0);
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    right_wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    lower_wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    atmosphere
    {
        type            fixedValue;        
        value           uniform (0 0 0);
    }
    
    default_faces
    {
        type            empty;
    }
}

TRANSPORT PROPERTIES
Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
phaseChange on;
 
phaseChangeTwoPhaseMixture Kunz;
 
sigma           sigma  [1 0 -2 0 0 0 0] 0.07;
 
DT   DT  [0 2 -1 0 0 0 0] 0.02; // [m^2/s]
 
phase1
{
    transportModel Newtonian;
    nu              nu [0 2 -1 0 0 0 0] 5.5e-07;
    rho             rho [1 -3 0 0 0 0 0] 988;
}
 
phase2
{
    transportModel Newtonian;
    nu              nu [0 2 -1 0 0 0 0] 4.273e-04;
    rho             rho [1 -3 0 0 0 0 0] 0.02308;
}
KunzCoeffs
{
    UInf            UInf   [0 1 -1 0 0 0 0]     5.0;
    tInf            tInf   [0 0 1 0 0 0 0]       0.016; // L = 0.08 m
    Cc              Cc     [0 0 0 0 0 0 0]      100;
    Cv              Cv     [0 0 0 0 0 0 0]      100;
}
MerkleCoeffs
{ 
    UInf            UInf   [0 1 -1 0 0 0 0]     5.0;
    tInf            tInf   [0 0 1 0 0 0 0]       0.016;  // L = 0.08 m
    Cc              Cc     [0 0 0 0 0 0 0]      1;
    Cv              Cv     [0 0 0 0 0 0 0]      1;
}
SchnerrSauerCoeffs
{
    n                n      [0 -3 0 0 0 0 0]     1.6e+13;
    dNuc           dNuc  [0 1 0 0 0 0 0]     2.0e-06;
    Cc              Cc     [0 0 0 0 0 0 0]      1;
    Cv              Cv     [0 0 0 0 0 0 0]      1;
}
In my results the temperature is going to a maximum of 370K when a value of 300.15K is specified on the right and bottom walls. Any idea where the error is???

Last edited by Jibran; January 18, 2013 at 08:39.
Jibran is offline   Reply With Quote

Old   April 8, 2013, 07:38
Default
  #2
abe
Member
 
ABE
Join Date: Jul 2012
Posts: 37
Rep Power: 5
abe is on a distinguished road
Which version of openfoam do you use?
Have you been able to get good results from new versions of OF (i.e. 2.1 or 2.2)?

ABE
abe is offline   Reply With Quote

Old   April 10, 2013, 05:52
Default
  #3
New Member
 
Jibran
Join Date: Oct 2012
Location: Swansea, UK
Posts: 25
Rep Power: 4
Jibran is on a distinguished road
I am using 2.1 however compressibleInterFoam of 2.2 might be a good starting point as it solves for energy equation.
Jibran is offline   Reply With Quote

Old   April 10, 2013, 05:56
Default
  #4
abe
Member
 
ABE
Join Date: Jul 2012
Posts: 37
Rep Power: 5
abe is on a distinguished road
Hi Jibran,

do you simulate cavitation?
abe is offline   Reply With Quote

Old   April 10, 2013, 06:36
Default
  #5
New Member
 
Jibran
Join Date: Oct 2012
Location: Swansea, UK
Posts: 25
Rep Power: 4
Jibran is on a distinguished road
Hi. I am working on interfacial evaporation and condensation. Are you doing something similar?
Jibran is offline   Reply With Quote

Old   April 10, 2013, 07:37
Default
  #6
abe
Member
 
ABE
Join Date: Jul 2012
Posts: 37
Rep Power: 5
abe is on a distinguished road
I think, the method used in interPhasechange solver is a relatively different. In this case evaporation or condensation can occur anywhere which has pressure lower than saturation pressure. While in the interfacial method, phase change just occur in the interface of the two phases (am I right?)
abe is offline   Reply With Quote

Old   April 10, 2013, 09:17
Default
  #7
New Member
 
Jibran
Join Date: Oct 2012
Location: Swansea, UK
Posts: 25
Rep Power: 4
Jibran is on a distinguished road
Yes that is correct. Therefore I am using a modified phase change model. By including terms such as "alpha1*(1-alpha1)" in the phase change equations one can ensure that the evaporation and condensation takes place just at the interface.
Jibran is offline   Reply With Quote

Old   April 10, 2013, 09:26
Default
  #8
abe
Member
 
ABE
Join Date: Jul 2012
Posts: 37
Rep Power: 5
abe is on a distinguished road
Do you have get any results?
are the both phases considered compressible or just the vapor phase?
abe is offline   Reply With Quote

Old   April 15, 2013, 05:29
Default
  #9
New Member
 
Jibran
Join Date: Oct 2012
Location: Swansea, UK
Posts: 25
Rep Power: 4
Jibran is on a distinguished road
I have some results which seem to be reasonable.
However, I need to validate them.
Yes both phases are compressible.
Jibran is offline   Reply With Quote

Old   September 19, 2013, 00:48
Default
  #10
Member
 
Mohammad Bahreini
Join Date: Dec 2012
Posts: 35
Rep Power: 4
mecman is on a distinguished road
Quote:
Originally Posted by Jibran View Post
Dear All,

Has anyone sucessfully implemented the temperature equation in interphaseChangeFoam?? I followed 'A interphasechangeFoam tutorial by Martin Andersen' but the results seem to be quite inaccurate. I have attached the relevant files below.

ALPHA1
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 0 0 0 0 0 0];
internalField   uniform 0;
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    
    right_wall
    {
        type            constantAlphaContactAngle;
        theta0          10;
        limit  gradient;
        value           uniform 0;
    }
    
    lower_wall
    {
        type            zeroGradient;
    }
    
    atmosphere
    {
        type  zeroGradient;    
    }
    
    default_faces
    {
        type            empty;
    }
}


P_RGH
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [1 -1 -2 0 0];
internalField   uniform 100000;
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    
    right_wall
    {
        type            fixedFluxPressure;
        adjoint         no;
    }
    
    lower_wall
    {
        type            buoyantPressure;
    }
    
    atmosphere
    { 
        type            fixedValue;
        value           uniform 100000;
    }
    
    default_faces
    {
        type            empty;
    }
}

T
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 0 0 1 0]; // [K]
internalField   uniform 280.15;
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    
    right_wall
    {
        type            fixedValue;
        value         uniform 300.15;
    }
    
    lower_wall
    {
        type            fixedValue;
        value         uniform 300.15;
    }
    
    atmosphere
    {
        type            fixedValue;   
        value         uniform 280.15;
    }
    default_faces
    {
        type            empty;
    }
}

U
Code:
 
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions      [0 1 -1 0 0];
internalField   uniform (0 0 0);
boundaryField
{
    axis
    {
        type            symmetryPlane;
    }
    right_wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    lower_wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    atmosphere
    {
        type            fixedValue;        
        value           uniform (0 0 0);
    }
    
    default_faces
    {
        type            empty;
    }
}

TRANSPORT PROPERTIES
Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
phaseChange on;
 
phaseChangeTwoPhaseMixture Kunz;
 
sigma           sigma  [1 0 -2 0 0 0 0] 0.07;
 
DT   DT  [0 2 -1 0 0 0 0] 0.02; // [m^2/s]
 
phase1
{
    transportModel Newtonian;
    nu              nu [0 2 -1 0 0 0 0] 5.5e-07;
    rho             rho [1 -3 0 0 0 0 0] 988;
}
 
phase2
{
    transportModel Newtonian;
    nu              nu [0 2 -1 0 0 0 0] 4.273e-04;
    rho             rho [1 -3 0 0 0 0 0] 0.02308;
}
KunzCoeffs
{
    UInf            UInf   [0 1 -1 0 0 0 0]     5.0;
    tInf            tInf   [0 0 1 0 0 0 0]       0.016; // L = 0.08 m
    Cc              Cc     [0 0 0 0 0 0 0]      100;
    Cv              Cv     [0 0 0 0 0 0 0]      100;
}
MerkleCoeffs
{ 
    UInf            UInf   [0 1 -1 0 0 0 0]     5.0;
    tInf            tInf   [0 0 1 0 0 0 0]       0.016;  // L = 0.08 m
    Cc              Cc     [0 0 0 0 0 0 0]      1;
    Cv              Cv     [0 0 0 0 0 0 0]      1;
}
SchnerrSauerCoeffs
{
    n                n      [0 -3 0 0 0 0 0]     1.6e+13;
    dNuc           dNuc  [0 1 0 0 0 0 0]     2.0e-06;
    Cc              Cc     [0 0 0 0 0 0 0]      1;
    Cv              Cv     [0 0 0 0 0 0 0]      1;
}
In my results the temperature is going to a maximum of 370K when a value of 300.15K is specified on the right and bottom walls. Any idea where the error is???
Hi Jibran and abe
Do you modified youre solver (interPhasechangrFoam) by adding T.eqn and source term?i want simulate evaporation and condensation and i use interPhasechangeFoam solver.i implemented the temperature equation( A interphasechangeFoam tutorial by Martin Andersen) but i dnt know how modified the solver by source term and T.eqn for phase change with increasing Temperature.
mecman is offline   Reply With Quote

Old   June 2, 2015, 07:52
Default
  #11
Senior Member
 
rkhr
Join Date: May 2011
Posts: 211
Rep Power: 7
Kanarya is on a distinguished road
Hi All,

I am trying to simulate condensation with interPhaseChangeFoam by adding scalar T. But it gives me the following error:
PHP Code:
time step continuity errors sum local 1.15855e-05, global = -1.19144e-16cumulative = -1.19144e-16
GAMG
:  Solving for pcorrInitial residual 1, Final residual 3.48266e+19No Iterations 50
GAMG
:  Solving for pcorrInitial residual 0.965555, Final residual 1.08509e+21No Iterations 50
time step continuity errors 
sum local 4.51124e+35, global = -1.41141e+28cumulative = -1.41141e+28
Courant Number mean
3.09991e+35 max5.41686e+38

Starting time loop

Courant Number mean
0.00057227 max1
deltaT 
1.84609e-47
Time 
1.84609e-47

#0  Foam::error::printStack(Foam::Ostream&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  
 
at sigaction.c:0
#3  
 
in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#4  
 
in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#5  
 
in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#6  __libc_start_main in "/lib64/libc.so.6"
#7  
 
in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
Floating point exception 
I am using following BC:
for p_rhg:
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
location    "0";
    
object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [--2 0 0 0 0];

internalField   uniform 1e5;

boundaryField
{
    
inlet
    
        
{
         
type            zeroGradient;
    }
    

   
outlet
    
{
        
type            fixedValue;
        
value           $internalField;
    }

    
  
hot_wall
    
{
         
type           fixedFluxPressure;
adjoint no;
    }
    
ocold_wall
   
{
         
type           fixedFluxPressure;
adjoint no;
    }
    
icold_wall
    

        
type           fixedFluxPressure;
adjoint no;
    }

}

// ************************************************************************* // 
for alpha1:
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
location    "0";
    
object      alpha1;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0];

internalField   uniform 0.00001;

boundaryField
{
    
inlet
    
{
         
type            fixedValue;
        
value           $internalField;
    }

    
outlet
    
{
        
type            inletOutlet;
        
inletValue      $internalField;
    }

    
    
hot_wall
    
{
        
type            zeroGradient;
    }
    
icold_wall
    
{
        
type            zeroGradient;
    }
    
ocold_wall
    
{
        
type            zeroGradient;
    }

}

// ************************************************************************* // 
and for U
PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volVectorField;
    
object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 5.199);

boundaryField
{
    
inlet
    
{
        
type            fixedValue;
        
value           $internalField;
    }

    
outlet
{
        
type            pressureInletOutletVelocity;
        
phi             phi;
        
value           $internalField;
}
    
hot_wall
    
{
        
type            fixedValue;
        
value           uniform (0 0 0);
    }

 
ocold_wall
    
{
        
type            fixedValue;
        
value           uniform (0 0 0);
    }
icold_wall
    
{
        
type            fixedValue;
        
value           uniform (0 0 0);
    }

}

// ************************************************************************* // 
and fr T

PHP Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    
version     2.0;
    
format      ascii;
    class       
volScalarField;
    
object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 366.85;

boundaryField
{
    
hot_wall
    
{
         
type            fixedValue;
        
value           uniform 293.15;
    }
   
ocold_wall
    
{
         
type            zeroGradient;
        
    }
icold_wall
    
{
         
type            zeroGradient;

    }

    
inlet
    
{
        
type            fixedValue;
        
value           uniform 366.85;
    }
    
outlet
    
{
        
type            zeroGradient;
    }
}

// ************************************************************************* // 
I am using version 2.1.0 for it...
can you hep me to fix the problem?

thanks!
Kanarya is offline   Reply With Quote

Old   June 2, 2015, 09:33
Default
  #12
Senior Member
 
rkhr
Join Date: May 2011
Posts: 211
Rep Power: 7
Kanarya is on a distinguished road
Hi Jibran,

Did you add this term (alpha1*(1-alpha1)) to alphaEqn.H?
Did you get the results?

Thanks

Quote:
Originally Posted by Jibran View Post
Yes that is correct. Therefore I am using a modified phase change model. By including terms such as "alpha1*(1-alpha1)" in the phase change equations one can ensure that the evaporation and condensation takes place just at the interface.
Kanarya is offline   Reply With Quote

Old   June 17, 2015, 05:10
Default
  #13
Member
 
Przemek
Join Date: Jun 2011
Posts: 42
Rep Power: 6
gaza is on a distinguished road
Hi Kanaraya
It would be better to use OF2.3.
You can try to decrease Courant number. If you use model from Andersen tutorial it give you wrong results. I try to do use his model for condensation on vertical wall and solver always crashed. It is due to bad model for mass flux.
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Old   June 17, 2015, 05:19
Default
  #14
Senior Member
 
rkhr
Join Date: May 2011
Posts: 211
Rep Power: 7
Kanarya is on a distinguished road
Hi Gaza,
thanks for your reply!
Did you manage to get another model?
Or how did you solve your problem?

thanks!

Quote:
Originally Posted by gaza View Post
Hi Kanaraya
It would be better to use OF2.3.
You can try to decrease Courant number. If you use model from Andersen tutorial it give you wrong results. I try to do use his model for condensation on vertical wall and solver always crashed. It is due to bad model for mass flux.
Kanarya is offline   Reply With Quote

Old   June 17, 2015, 10:47
Default
  #15
Member
 
Przemek
Join Date: Jun 2011
Posts: 42
Rep Power: 6
gaza is on a distinguished road
Hi Kanarya,
I did not solve my problem. I only know that interPhaseChangeFoam is good for cavitation mass transfer models. For heat transfer one has to change the main code a bit and write own functions
for mDotAlphal and mDotP.
Solver is constructed in that way that in alpha equation (alphaEqn.H) condensation mass flux is multiplied by (1 - alpha) and evaporation mass flux is multiplied by alpha. So you have to: 1) adjust your implementation for mass fluxes to the solver or 2) implement them as they should be and change the solver (alpha eqn. and lines in MULES where vDot variables are used).
The same is true about mDotP function.

I think that it would be better to implement model based on TSat not for pSat.
__________________
best regards
pblasiak
gaza is offline   Reply With Quote

Reply

Tags
energy equation, interphasechangefoam, temperature

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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Eliminating exchange energy between phase in Granular Temperature diff eqn. cfdbusiness FLUENT 0 February 2, 2010 23:30
Please explain the implementation of species transport Eqn in reactingFoam kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 07:29
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 13, 2000 23:03
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 11, 2000 03:43


All times are GMT -4. The time now is 19:30.