Hi all,
I’m trying to simulate a multiphase FSI case (water seeping into soil). The soil part is removed in fluid mesh and it’s 0.1m wide and 0.2m high. For solid part, the cell size in x direction on the interface is 0.01m, which is two times of fluid cell size. The water level is 0.1m above the soil’s top surface. The pressure on the interface should be 981Pa. However, the solid solver didn’t get the right Pressure.
The BC for fluid part and preciceDict are as follows:
alpha.water
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5-dev |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
“wall.”
{
type zeroGradient;
}
seepageTest
{
type zeroGradient;
}
frontAndBack
{
type zeroGradient;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
}
// ************************************************************************* //
pointDisplacement
FoamFile
{
version 2.0;
format ascii;
class pointVectorField;
object pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 0 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type slip;
}
outlet
{
type slip;
}
seepageTest
{
type fixedValue;
value $internalField;
}
atmosphere
{
type slip;
}
"wall."
{
type slip;
}
frontAndBack
{
type slip;
}
}
// ************************************************************************* //
p_rgh
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5-dev |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
“wall.”
{
type fixedFluxPressure;
value uniform 0;
}
seepageTest
{
type fixedFluxPressure;
value uniform 0;
}
frontAndBack
{
type fixedFluxPressure;
value uniform 0;
}
outlet
{
type fixedFluxPressure;
value uniform 0;
}
inlet
{
type fixedFluxPressure;
value uniform 0;
}
atmosphere
{
type totalPressure;
U U;
phi phi;
rho rho;
psi none;
gamma 1;
p0 uniform 0;
value uniform 0;
}
}
// ************************************************************************* //
U
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.7.x |
| \ / A nd | Web: www.OpenFOAM.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location “0”;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 0);
}
atmosphere
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
“wall.”
{
type fixedValue;
value uniform (0 0 0);
}
seepageTest
{
type fixedValue;
value $internalField;
}
frontAndBack
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type fixedValue;
value uniform (0 0 0);
}
}
// ************************************************************************* //
preciceDict
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location “system”;
object preciceDict;
}
preciceConfig “…/precice-config.xml”;
participant Fluid;
modules (FSI);
interfaces
{
Interface1
{
mesh Fluid-Mesh;
patches (seepageTest);
locations faceCenters;
readData
(
);
writeData
(
Pressure
);
};
Interface2
{
mesh Fluid-Mesh;
patches (seepageTest);
locations faceNodes;
readData
(
Displacement
Velocity
);
writeData
(
);
};
};
FSI
{
rho rho [1 -3 0 0 0 0 0] 1;
}