[pruthvi1991]

Hello everybody.

I'm trying to simulate a flapping airfoil with plunge amplitudes as high as 2 chords. Mesh deformation seems to be infeasible for this purpose. Therefore I want to implement the method used by Kinsey and Dumas (2008). They have added a time-varying momentum source term to the Ueqn and also time-varying velocity inlets to the top,bottom and inlet patches. The body force term is independent of space. This is basically a solver in a non-inertial frame of reference. The pitching motion is taken care of by a rotating AMI.

I want to validate the pimpleFoam solver in the non-inertial frame before I proceed to pimpleDyMFoam. I modified the solver as follows by following some previous posts in the forum.

I added an explicit source term to the right hand side of Ueqn as follows

solve(UEqn() == -g -fvc::grad(p));

g is the time varying body force defined as follows

 g=gunits*Foam::cos(runTime.value()*freQ)*vector(0,1,0);

I ran a simulation using the modified solver which is available here. The flow field seems to be reasonable but the value of forceCoeffs is much smaller than the actual value and also very erratic. Check out the forceCoeff plots below.

Cy_non_inertial.png // Non-inertial plot

forceCoeffs_flatplate.jpg // Dynamic Mesh plot

I want to know if there is any mistake in this implementation.

Here is a flow field of the same case using a dynamic mesh.

Do I need to add the body force term to the fvSchemes file? I didn't because 'g' is just a known constant.

Did anybody implement this before? I will be glad to hear your suggestions.

The case file is available here.

Thanks,
Pruthvi.

[kebsiali]

I dont know if what im saying makes sense or not but since you want any suggestion

When changing the momentum equation, shouldn't the pressure correction and flux correction be changed as well?

[pruthvi1991]

Hi Ali,

I don't know if I should or shouldn't. Can you explain why I should?

[JSR]

Do you solve that issue? I am also working on the same type of problem.

[Msure]

Originally Posted by JSR

Do you solve that issue? I am also working on the same type of problem.

Hello,

Have you solved the problem?
Thanks.