[matin]

Dear Foamers,

I'm currently using DPMFoam and MagneticFoam to simulate flow containing particles in a magnetic field. To do this I have to first create the volVectorField HdotGradH, which I already did. I also have the particles in the channel. In order to apply the magnetic field I need to add the paramagnetic force to the kinematicCloudProperties directory, which Includes magneticSusceptibility & volvectorField HdotGradH.
I created my submodel as follow:

particleForces
{
WenYuDrag
{
alphac alpha.water;
}
gravity;

paramagnetic
{
magneticSusceptibility 6e-4; //m^3/kg
HdotGradH HdotGradH; //magnetic force
}

when I run the case, it gives the following error:

request for volVectorField HdotGradH from objectRegistry region0 failed
available objects of type volVectorField are
1(U.water)

From function objectRegistry::lookupObject<Type>(const word&) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.

FOAM aborting

#0 Foam::error:: srintStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const& Foam:: objectRegistry::lookupObject<Foam::GeometricField< Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#3 Foam::ParamagneticForce<Foam::KinematicCloud<Foam: :Cloud<Foam::CollidingParcel<Foam::KinematicParcel <Foam:: particle> > > > >::cacheFields(bool) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/liblagrangianIntermediate.so"
#4
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/DPMFoam"
.
.
My problem is that Openfoam does not recognize the HdotGradH file in the Zero directory and only sees the U.water as the volvectorfield !!

I tried every possible format, like the following:

HdotGradH HdotGradH.force;
HdotGradHName HdotGradH;
.
.
.
but the only case that program could run is:
HdotGradH U.water;

but it's not correct, cause I need to have magnetic field, not velocity field !

I would really appreciate if anyone could help me with this problem to see why OpenFoam does not see the HdotGradH file in the Zero directory.

Any help or reference regarding DPMFoam/MPPICFoam would be highly appreciated.

Best regards,

Matin

[vonboett]

although a late reply I'd like to note that it is helpful to look into createFields.H of the solver to see how fields are added to the object registry...