[baila12]

I am currently using OpenFOAM 2.0. I generated a mesh of a NACA0012 airfoil using gmsh, then ran the LTSReactingParcelFoam solver on the airfoil. The solver seemed to run correctly as far as the values for p, U, T, etc. Then I ran the foamToVTK command and opened paraview. I next opened the vtk main file, then opened the reacting cloud vtk file within the lagrangian directory and applied the glyph filter to it.

The issue I am running into is that the particle trajectories appear to be going through the airfoil rather than around it. Does anyone have any idea what might be causing this to happen?

[Dommy]

Hello baila12!

Pretty sure your problem is caused by your settings in the constant/reactingCloud1Properties file!

You may change the patchInteractionModel to localInteraction

patchInteractionModel     localInteraction;

and define the interaction by hand - for instance:

 localInteractionCoeffs
    {
        patches
        (
            "WALL"
            {
                type        rebound;
                e                 1;
                mu              0;
            }
            "(INLET|OUTLET)"
            {
                type        escape;
            }
        );
    }

You have to choose between rebound, escape and stick!

Best,
Thomas

[baila12]

Thank you for your reply!

I tried changing the patchInteractionModel to localInteraction and then defining the interaction, but it seems that the particles are still going through the airfoil.
I have the airfoil set to stick and the rest of the patches set to escape. I also tried setting top, bottom, front, and back to rebound and inlet and outlet to escape.
I have attached the reactingCloud1Properties file.

[Dommy]

Can you upload the whole case including mesh? I never did 2D with Euler/Langrange before. Therefore its hard to see only from this file.

What do you want to acchieve? If you are only interested in the motion of the particles and find yourself in the mass load region of one way coupling (which I think is the case) you might switch solvers to a steady state solution for the flow + icoUncoupledParcelFoam.

Its a solver for passive transport of particles by interacting with the flow - but does not solve equations for the flow -> 1-way-coupling

Best,
Thomas

[baila12]

I would like to find the mass flux of water on the airfoil using the lagrangian information.
Even in compressed format the case file is too large to attach here, is there somewhere else I can send it?
Thank You!

[Dommy]

I just need the case basic, cleaned structure (0/ system/ constant/ and mesh or constant with polymesh folder) - I can run it alone. Or you upload it somewhere else (Dropbox, UbuntuOne,...) and post the link here.

I would solve the case with a steady state solver (eg simpleFoam) and use the field for icoUncoupledPracelFoam with a passive transport. By applying the stick interaction you know the mass added to your system and the mass leaving the system - the difference should stick on your airfoil.

In my opinion there is no need to fire with reactingParcelFoam because you are just interested in the distribution of the particles (all reaction models,... are turned off) in front of the airfoil. So a steady state solution might be sufficient.

(Of course reactingParcelFoam might come to the same results but will use much more CPU time..)

Best,
Thomas

[baila12]

Here is the link:
https://drive.google.com/folderview?...G8&usp=sharing
I would like to stick with trying to get the reactingFoam solver to work correctly on the airfoil if possible because of future development, but if it isn't possible I may try that approach!

[baila12]

I discovered that the issue may be that the solver is set to run on porous material, but am unsure of how to fix this.