[jeytsav]

Hello people

I would like to ask what is the way to initialize a spherical farfield domain?

Using the following BCs I get a weird initialization with potentialFoam (OFv4.1) as shown below

view1.PNG view2.PNG view3.PNG

It is about a simple spherical body into a farfield domain. I have only 3 boundaries (see code)



Here is my code


Boundaries "/constant/polyMesh/boundary"

FoamFile
{
	version 2.0;
	format ascii;
	class polyBoundaryMesh;
	location "";
	object boundary;
}

3
(

	farfield
	{
		type patch;
		startFace 77179;
		nFaces 3670;
	}

	body
	{
		type wall;
		startFace 80849;
		nFaces 162;
	}

	symmetry
	{
		type symmetry;
		startFace 81011;
		nFaces 2438;
	}

)

Velocity boundary condition "/0/U"

FoamFile
{
	version 2.0;
	format ascii;
	class volVectorField;
	location "";
	object U;
}

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

internalField uniform ( 2.5 0. 0. );


boundaryField
{
	farfield
	{
		type inletOutlet;
		inletValue uniform (2.5 0. 0.);
		value	 uniform (2.5 0. 0.);
	}

	body
	{
		type fixedValue;
		value	 uniform (0. 0. 0.);
	}

	symmetry
	{
		type symmetry;
	}


}

Pressure boundary conditions "0/p"

FoamFile
{
	version 2.0;
	format ascii;
	class volScalarField;
	location "";
	object p;
}

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.;

boundaryField
{
	farfield
	{
		type outletInlet;
                value uniform 0;
                outletValue uniform 0;

	}

	body
	{
		type zeroGradient;
	}

	symmetry
	{
		type symmetry;
	}


}

The velocity does not advance inside the domain (as would be expected) and just stays in the half boundary (following strange patterns).

Does anyone know what is the correct way to initialize such a field.? I have tried various boundary conditions until now.

Thanks everyone in advance for spending their time for this!

Have a nice day!