 |
 |
 |
|
|
|
[Sponsors] | ||||
Water level is raising as simulation progresses |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 31, 2022, 09:31
|
Water level is raising as simulation progresses
|
#1 |
|
New Member
Akshay
Join Date: Jan 2020
Posts: 28
Rep Power: 6  |
I am conducting two-phase flow simulations in a separator with 3 inlets and 2 outlets using twoPhaseEulerFoam in OF version 8.
The inlet velocity BC is specified at 2 inlets and pressure outlet BC is specified at the 2 outlets. Initially, half of the separator volume is filled with water using setFieldDict. I am expecting the water level to remain constant but strangely the water level increases. I was expecting flow in = flow out, the velocity vector at the water outlet is not correctly calculating the water velocity at the outlet. I have attached my BC for U, p_rgh, alpha fields below. Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet1
{
type fixedValue;
value uniform 0.698;
}
inlet2
{
type fixedValue;
value uniform 0.698;
}
inlet3
{
type fixedValue;
value uniform 0.698;
}
walls
{
type zeroGradient;
}
baffle1
{
type zeroGradient;
}
baffle2
{
type zeroGradient;
}
baffle3
{
type zeroGradient;
}
outlet1
{
type inletOutlet;
phi phi.water;
inletValue uniform 1;
value uniform 1;
}
outlet2
{
type inletOutlet;
phi phi.water;
inletValue uniform 0;
value uniform 0;
}
}
// ************************************************************************* //
Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alpha.oxygen;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 1;
boundaryField
{
inlet1
{
type fixedValue;
value uniform 0.302;
}
inlet2
{
type fixedValue;
value uniform 0.302;
}
inlet3
{
type fixedValue;
value uniform 0.302;
}
walls
{
type zeroGradient;
}
baffle1
{
type zeroGradient;
}
baffle2
{
type zeroGradient;
}
baffle3
{
type zeroGradient;
}
outlet1
{
type inletOutlet;
phi phi.oxygen;
inletValue uniform 0;
value uniform 0;
}
outlet2
{
type inletOutlet;
phi phi.oxygen;
inletValue uniform 1;
value uniform 1;
}
}
// ************************************************************************* //
Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 4.01325e5;
boundaryField
{
inlet1
{
type fixedFluxPressure;
value $internalField;
}
inlet2
{
type fixedFluxPressure;
value $internalField;
}
inlet3
{
type fixedFluxPressure;
value $internalField;
}
walls
{
type fixedFluxPressure;
value $internalField;
}
baffle1
{
type fixedFluxPressure;
value $internalField;
}
baffle2
{
type fixedFluxPressure;
value $internalField;
}
baffle3
{
type fixedFluxPressure;
value $internalField;
}
outlet1 //water outlet
{
type prghPressure;
p $internalField;
value $internalField;
}
outlet2 //oxygen outlet
{
type fixedValue;
value $internalField;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Code:
FoamFile
{
version 2.0;
format binary;
class volVectorField;
object U.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet1
{
type fixedValue;
value uniform (0 0 3.83);
}
inlet2
{
type fixedValue;
value uniform (0 0 -3.83);
}
inlet3
{
type fixedValue;
value uniform (0 0 -3.83);
}
walls
{
type noSlip;
}
baffle1
{
type noSlip;
}
baffle2
{
type noSlip;
}
baffle3
{
type noSlip;
}
outlet1
{
type pressureInletOutletVelocity;
phi phi.water;
value $internalField;
}
outlet2
{
type pressureInletOutletVelocity;
phi phi.water;
value $internalField;
}
}
// ************************************************************************* //
Code:
FoamFile
{
version 2.0;
format binary;
class volVectorField;
object U.oxygen;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet1
{
type fixedValue;
value uniform (0 0 1.66);
}
inlet2
{
type fixedValue;
value uniform (0 0 -1.66);
}
inlet3
{
type fixedValue;
value uniform (0 0 -1.66);
}
walls
{
type fixedValue;
value uniform (0 0 0);
}
baffle1
{
type fixedValue;
value uniform (0 0 0);
}
baffle2
{
type fixedValue;
value uniform (0 0 0);
}
baffle3
{
type fixedValue;
value uniform (0 0 0);
}
outlet1
{
type pressureInletOutletVelocity;
phi phi.oxygen;
value $internalField;
}
outlet2
{
type pressureInletOutletVelocity;
phi phi.oxygen;
value $internalField;
}
}
// ************************************************************************* //
out.0000.png out.0399.png Thank you. |
|
|
|
|
|
June 23, 2022, 02:40
|
I have modifies the boundary conditions
|
#2 |
|
New Member
Akshay
Join Date: Jan 2020
Posts: 28
Rep Power: 6 |
It is really difficult to maintain a water level inside the separator, I am playing with certain boundary conditions. Please let me know if my present adopted boundary condition is good.
Code:
{
version 2.0;
format binary;
class volVectorField;
object U.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet1
{
type fixedValue;
value uniform (0 0 3.83);
}
inlet2
{
type fixedValue;
value uniform (0 0 -3.83);
}
inlet3
{
type fixedValue;
value uniform (0 0 -3.83);
}
walls
{
type noSlip;
}
pipes
{
type noSlip;
}
baffle1
{
type noSlip;
}
baffle2
{
type noSlip;
}
baffle3
{
type noSlip;
}
outlet1
{
type outletPhaseMeanVelocity;
UnMean 3.07;
alpha alpha.water;
value uniform (0 0 0);
}
outlet2
{
type zeroGradient;
}
}
Code:
FoamFile
{
version 2.0;
format binary;
class volVectorField;
object U.oxygen;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet1
{
type fixedValue;
value uniform (0 0 1.66);
}
inlet2
{
type fixedValue;
value uniform (0 0 -1.66);
}
inlet3
{
type fixedValue;
value uniform (0 0 -1.66);
}
walls
{
type noSlip;
}
pipes
{
type noSlip;
}
baffle1
{
type noSlip;
}
baffle2
{
type noSlip;
}
baffle3
{
type noSlip;
}
outlet1
{
type zeroGradient;
}
outlet2
{
type pressureInletOutletVelocity;
phi phi.oxygen;
value $internalField;
}
}
// ************************************************************************* //
Code:
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 4.01325e5;
boundaryField
{
inlet1
{
type fixedFluxPressure;
value $internalField;
}
inlet2
{
type fixedFluxPressure;
value $internalField;
}
inlet3
{
type fixedFluxPressure;
value $internalField;
}
walls
{
type fixedFluxPressure;
value $internalField;
}
pipes
{
type fixedFluxPressure;
value $internalField;
}
baffle1
{
type fixedFluxPressure;
value $internalField;
}
baffle2
{
type fixedFluxPressure;
value $internalField;
}
baffle3
{
type fixedFluxPressure;
value $internalField;
}
outlet1
{
type fixedFluxPressure;
value $internalField;
}
outlet2
{
type fixedValue;
value $internalField;
}
/*
{
type totalPressure;
p0 $internalField;
value $internalField;
// Optional entries
U U.oxygen;
phi phi;
rho rho;
}
*/
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Code:
Mesh stats
points: 188392
faces: 539938
internal faces: 516140
cells: 176193
faces per cell: 5.99387
boundary patches: 10
point zones: 0
face zones: 7
cell zones: 0
Overall number of cells of each type:
hexahedra: 164588
prisms: 7831
wedges: 0
pyramids: 0
tet wedges: 67
tetrahedra: 0
polyhedra: 3707
Breakdown of polyhedra by number of faces:
faces number of cells
4 467
5 413
6 836
7 45
8 62
9 1230
10 15
11 11
12 475
13 5
14 4
15 127
16 3
18 14
Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
baffle1 1686 1713 ok (non-closed singly connected)
baffle2 1670 1702 ok (non-closed singly connected)
baffle3 1684 1715 ok (non-closed singly connected)
inlet1 59 67 ok (non-closed singly connected)
inlet2 60 68 ok (non-closed singly connected)
inlet3 57 63 ok (non-closed singly connected)
outlet1 218 232 ok (non-closed singly connected)
outlet2 16 18 ok (non-closed singly connected)
pipes 1771 1985 ok (non-closed singly connected)
walls 16577 17438 ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (-0.108 -0.348498 -0.4122) (2.8772 0.348494 0.4122)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (-2.41681e-17 5.65559e-16 3.60538e-16) OK.
Max cell openness = 2.84095e-16 OK.
Max aspect ratio = 6.71026 OK.
Minimum face area = 6.42656e-07. Maximum face area = 0.000725531. Face area magnitudes OK.
Min volume = 5.69092e-08. Max volume = 1.11177e-05. Total volume = 1.04276. Cell volumes OK.
Mesh non-orthogonality Max: 44.9781 average: 5.05001
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.95597 OK.
Coupled point location match (average 0) OK.
Mesh OK.
End
|
|
|
|
|
|
|
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Solid object falling in water tank simulation | cassini83 | CFX | 3 | August 9, 2021 17:05 |
| postProcessing to visualize areas touched by water during a sprinkling simulation | TheBigLebowski | OpenFOAM | 0 | August 1, 2021 10:20 |
| Water level doesn't rise in CFD-DEM simulation of air-water-solid. What should I chec | gsiitd | OpenFOAM | 0 | January 23, 2019 00:25 |
| For a single phase flow, is there any way to check water surface level? | Tanjina | FLUENT | 0 | February 23, 2014 22:56 |
| Water subcooled boiling | Attesz | CFX | 7 | January 5, 2013 03:32 |
Linear Mode
