[13msmemusman]

i entered following equation in interFoam.C
CoM = sum(rho*mesh.V()*mesh.C().dimensionedInternalField ())/sum(rho*mesh.V());
then defined it in creatFields.H as following
dimensionedVector CoM("CoM", dimLength, vector::zero);
then i compiled it works well. but now i want to take its gradient. i wrote
gradCoM = grad (CoM);
and defined it in creatFields.H as following
uniformDimensionedVectorField gradCoM("gradCoM", dimLength, vector::zero);
now it gives following error please help
In file included from myInterFoamDDT.C:60:0:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:2:73: error: no matching function for call to ‘Foam::UniformDimensionedField<Foam::Vector<double > >::UniformDimensionedField(const char [8], const Foam::dimensionSet&, const Foam::Vector<double>&)’
uniformDimensionedVectorField gradCoM("gradCoM", dimLength, vector::zero);
one thing is obvious in my views that problem is in defining the field as UniformDimensionedField...... please help

[bkassar]

Hi Muhammad,
The compiler is complaining about the constructor of uniformDimensionedVectorField.
However, I think your main problem is not that.
The main issue is that you are trying to compute the gradient of a vector. Notice that the variable 'CoM' is not a Field.
cheers,
Bruno

[demichie]

In addition if you have a field A the gradient of A should have the units of A divided by the unit of lenght.

Ciao
Mattia

[13msmemusman]

Sir you are right that its not a field. i am taking gradient of a vector. then what should i use instead of DimensionedVector???? actually i read that gradient of a vector is a tensor thats why i wrote that

[bkassar]

Hi Muhammad,
As far as I see, you computed the center of mass and stored it in 'CoM'.
Could you tell us what you want to compute?
It is not clear to me why you need the gradient...
cheers,
Bruno

[13msmemusman]

Sir i want to take double derivative of every component of this vector as i have to put it in a formulae

[13msmemusman]

Sir i want to take double derivative of every component of this vector as i have to put it in a formulae. i am sorry at first i was wrong. i was thinking that formulae requires double spacial derivative of COM completely. i want to take double derivative of every component of CoM. thank you sir

[bkassar]

I see. So, by double derivative you mean double derivative with respect to time? You want the acceleration of this center of mass?
If the answer is yes, you could store the values of CoM in three different time steps and perform the derivatives by a second-order differencing scheme you may choose.
Hope that helps.

[13msmemusman]

sir i stored components of Varible CoM in 3 different varibles in following way

CoM = sum(rho*mesh.V()*mesh.C().dimensionedInternalField ())/sum(rho*mesh.V());
Info << "Center of Mass " << CoM.dimensions() << " = " << CoM << nl << endl;
vector CoMx = CoM.x();
vector CoMy = CoM.y();
vector CoMz = CoM.z();


but it gives following error

myInterFoamDDT.C: In function ‘int main(int, char**)’:
myInterFoamDDT.C:111:19: error: ‘Foam::dimensionedVector’ has no member named ‘x’
vector CoMx = CoM.x();
^
myInterFoamDDT.C:112:19: error: ‘Foam::dimensionedVector’ has no member named ‘y’
vector CoMy = CoM.y();
^
myInterFoamDDT.C:113:19: error: ‘Foam::dimensionedVector’ has no member named ‘z’


Please help me.

[13msmemusman]

please help

[bkassar]

Hi Muhammad,
The second derivative you want to take is with respect to time?

[13msmemusman]

yes sir its with respect to time

[13msmemusman]

The problem specified above is solved.... now i just want to take its derivative with respect to time. means i need acceleration of center of mass.

[bkassar]

So now, you just need the value of CoM at three consecutive time steps; let's say CoM(i-2), CoM(i-1) and CoM(i). Choose a second-order differencing scheme and use these three values to obtain the acceleration. There are many ways to do it. The easiest one would be the second-order backward scheme. For your problem, it would give something like:
A = (CoM(i)-2*CoM(i-1)+CoM(i-2))/(dt^2)
Hope that helps!
cheers,
Bruno

[13msmemusman]

Thank you very much sir....