[kk415]

Hello All


I am new to OpenFoam. I have this small portion of code in my solver. I think right now it is working in a very efficient. Can anyone suggest me any more efficient way of doing the same thing.


CODE


forAll(solidRegions[0].cells(), cellI)
{
vector r=solidRegions[0].C()[cellI]-fluidRegions[0].C()[celladd_to_src[cellI]];
Us[cellI]=Uf[celladd_to_src[cellI]];
Us[cellI].component(0)=Us[cellI].component(0)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(0)*r.component(0);
Us[cellI].component(1)=Us[cellI].component(1)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(1)*r.component(0);
Us[cellI].component(2)=Us[cellI].component(2)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(2)*r.component(0);
Us[cellI].component(0)=Us[cellI].component(0)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(3)*r.component(1);
Us[cellI].component(1)=Us[cellI].component(1)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(4)*r.component(1);
Us[cellI].component(2)=Us[cellI].component(2)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(5)*r.component(1);
Us[cellI].component(0)=Us[cellI].component(0)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(6)*r.component(2);
Us[cellI].component(1)=Us[cellI].component(1)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(7)*r.component(2);
Us[cellI].component(2)=Us[cellI].component(2)+fvc::grad(Uf)()[celladd_to_src[cellI]].component(8)*r.component(2);
Tf[celladd_to_src[cellI]]=0;
}

[dlahaye]

Possibly it helps to look into Chapter 9 of Moukalled, Manga and Darwish, The finite volume method in computational fluid dynamics, Springer 2015 , http://www.springer.com/gp/book/9783319168739

[msaravia]

If I am right you need to update Us as:

Us += fvc::grad(Uf).T() & r

Writing this would probably be faster, in your loop you are evaluating the gradient several times and you should do it only once.