[BSengupta]

Hi,


Let me first try to describe what I wish to do:


I am trying to make a mesh deformation solver. I am working on aircraft icing. After all the other calculations are done, the thickness of the ice is obtained as output. I need of deform or morph the mesh according to thickness of ice accreted on the surface.



I wrote the following code:

volScalarField UNn_X = UN.component(vector::X);
        volScalarField UNn_Y = UN.component(vector::Y);
        volScalarField UNn_Z = UN.component(vector::Z);
        
        UNn = (-UNn_X*i) + (-UNn_Y*j) + (-UNn_Z*k);
        
        pointMesh pMesh(mesh);
        volPointInterpolation vpi(mesh);
         pointVectorField Um = vpi.interpolate(UNn);
         pointScalarField Thp = vpi.interpolate(Thick);
         //~ pointVectorField pointDisplacement = Thp*Um;
         pointDisplacement_ = Thp*Um/lref;
         
        pointDisplacement_.correctBoundaryConditions();
        Info<< "pointDisplacement_ = " << pointDisplacement_ << endl;
        runTime++;
        Info<< "Time = " << runTime.timeName() << nl << endl;

        Info<<"deformationFOAM START"<< endl;
        
        Info<<"Mesh Moving"<<endl;
        mesh.update();
         mesh.write();

Here, UN is the normal vector to boundary.

But the code do not calculate pointDisplacement automatically, but when the displacement is given manually it deforms the mesh as expected. But I need to generate pointDisplacement automatically according to given distance.



In createFields, I defined the quantities as following:

volScalarField Thick
    (
        IOobject
        (
           "Thick",
            runTime.timeName(),
            mesh,
            IOobject::READ_IF_PRESENT,
            IOobject::AUTO_WRITE
         ),
         mesh,
         dimensionedScalar("Thick",dimensionSet(0,1,0,0,0),0.),
         zeroGradientFvPatchScalarField::typeName
    );
    
    pointVectorField pointDisplacement_
    (
    IOobject
        (
            "pointDisplacement",
            mesh.time().timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        pointMesh::New(mesh),
        dimensionedVector("pointDisplacement",dimensionSet(0,1,0,0,0,0,0),vector::zero),
        zeroGradientFvPatchVectorField::typeName
    );

dimensionedScalar lref("A",dimensionSet(0,1,0,0,0,0,0), 1.0);

I am using foam-extend 4.0.


I shall be really grateful and helped if anyone can suggest me necessary changes in the code to achieve what I need.


Thanks a lot,
Bidesh