[AshwaniAssam]

Hi,


I am trying to modify sonicFOAM. I am trying to create a source term (phi_g) to momentum equation which is d(phi_g)/dx. Where phi_g is function of coordinate x. It is created in the createFields.h in the following manner.

volScalarField phi_g
(
IOobject
(
"phi_g",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
(mesh.C().component(0))*dummy
);

Now, in the pressure equation it will become Laplacian(phi_g).


Now, when I compile the code it compiles fine. But while running a case it gives the following error:

--> FOAM FATAL ERROR:
cannot be called for a calculatedFvPatchField
on patch sides of field phi_g in file "/home/asd/OpenFOAM/asd-v1906/run/shockTube/constant/phi_g"
You are probably trying to solve for a field with a default boundary condition.

From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientIntern alCoeffs() const [with Type = double]
in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 188.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::calculatedFvPatchField<double>::gradientInte rnalCoeffs() const at ??:?
#3 Foam::fv::gaussLaplacianScheme<double, Foam::SymmTensor<double> >::fvmLaplacianUncorrected(Foam::GeometricField<do uble, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#4 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 ? at ??:?
#6 ? at ??:?
#7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8 ? at ??:?
Aborted (core dumped)

My question is that the boundary condition for the new scalar field should be calculated in the same manner as its definition.

Do I need to specify a boundary condition for phi_g in 0 folder with types calculated?