CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   sequential injection stimulation using the multiphaseInterFoam (https://www.cfd-online.com/Forums/openfoam-solving/123189-sequential-injection-stimulation-using-multiphaseinterfoam.html)

salame September 6, 2013 09:17

sequential injection stimulation using the multiphaseInterFoam
 
Hi, Formers:
these days, i used the multiphaseInterFoam solver for stimulating the sequential injection. but the results not OK. i think the problems emerge in the alpha setting. I want to stimulate first inject O(fluid) in the cavity, some time later, stop O and inject the M (another fluid) with the same inlet.
Relevant files below:
First, inject O.......
1 alphaair
dimensions [0 0 0 0 0 0 0];

internalField uniform 1;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

2 alphaM
dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

3 alphaO
dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

4 setFields
efaultFieldValues
(
volScalarFieldValue alphaair 1
volScalarFieldValue alphaO 0
volScalarFieldValue alphaM 0

);

regions
(
boxToCell
{
box (-0.01 -0.01 0.001) (0.001 0.002 0.003);
fieldValues
(

volScalarFieldValue alphaO 1
volScalarFieldValue alphaM 0
volScalarFieldValue alphaair 0
);
}

Second, after injecion 0.6s, stop O and inject M, I change the alphaO and alphaM in the 0.6 result file without changing others. like:
alphaM:
dimensions [0 0 0 0 0 0 0];

internalField nonuniform List<scalar>
3620
(
-1.34342e-09
-1.3543e-09
...............
)
;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

alphaO:
dimensions [0 0 0 0 0 0 0];

internalField nonuniform List<scalar>
3620
(
-3.24112e-10
-3.36625e-10
.....................
)
;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

the result is (writeInterval 0.2s): 1. 0~0.8s, have no change in the alpha (O, M, air)
http://xiangce.baidu.com/picture/det...dfef5441cc5425
2, from 1s, the alpha changed
1s:http://xiangce.baidu.com/picture/det...6a113f59070475
1.2s:http://xiangce.baidu.com/picture/det...68e8d92babf075
So, can u kindly help me solve the problem?

salame September 11, 2013 09:38

the pictures
1.http://xiangce.baidu.com/picture/det...dfef5441cc5425
2.http://xiangce.baidu.com/picture/det...6a113f59070475
3.http://xiangce.baidu.com/picture/det...68e8d92babf075


Hi, Formers:
these days, i used the multiphaseInterFoam solver for stimulating the sequential injection. but the results not OK. i think the problems emerge in the alpha setting. I want to stimulate first inject O(fluid) in the cavity, some time later, stop O and inject the M (another fluid) with the same inlet.
Relevant files below:
First, inject O.......
1 alphaair
dimensions [0 0 0 0 0 0 0];

internalField uniform 1;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

2 alphaM
dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

3 alphaO
dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

4 setFields
efaultFieldValues
(
volScalarFieldValue alphaair 1
volScalarFieldValue alphaO 0
volScalarFieldValue alphaM 0

);

regions
(
boxToCell
{
box (-0.01 -0.01 0.001) (0.001 0.002 0.003);
fieldValues
(

volScalarFieldValue alphaO 1
volScalarFieldValue alphaM 0
volScalarFieldValue alphaair 0
);
}

Second, after injecion 0.6s, stop O and inject M, I change the alphaO and alphaM in the 0.6 result file without changing others. like:
alphaM:
dimensions [0 0 0 0 0 0 0];

internalField nonuniform List<scalar>
3620
(
-1.34342e-09
-1.3543e-09
...............
)
;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

alphaO:
dimensions [0 0 0 0 0 0 0];

internalField nonuniform List<scalar>
3620
(
-3.24112e-10
-3.36625e-10
.....................
)
;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
down
{
type symmetryPlane;
}
wall
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}

the result is (writeInterval 0.2s): 1. 0~0.8s, have no change in the alpha (O, M, air)
http://xiangce.baidu.com/picture/det...dfef5441cc5425
2, from 1s, the alpha changed
1s:http://xiangce.baidu.com/picture/det...6a113f59070475
1.2s:http://xiangce.baidu.com/picture/det...68e8d92babf075
So, can u kindly help me solve the problem?[/QUOTE]


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