CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

LTSReactingParcelFoam Solver on an Airfoil

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 18, 2013, 12:43
Default LTSReactingParcelFoam Solver on an Airfoil
  #1
New Member
 
Cay Myers
Join Date: Jun 2013
Posts: 14
Rep Power: 5
baila12 is on a distinguished road
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?
baila12 is offline   Reply With Quote

Old   July 18, 2013, 16:54
Default
  #2
New Member
 
Thomas F.
Join Date: Jun 2013
Location: Graz, Austria
Posts: 14
Rep Power: 5
Dommy is on a distinguished road
Hello baila12!

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

You may change the patchInteractionModel to localInteraction
Code:
patchInteractionModel     localInteraction;
and define the interaction by hand - for instance:

Code:
 localInteractionCoeffs
    {
        patches
        (
            "WALL"
            {
                type        rebound;
                e                 1;
                mu              0;
            }
            "(INLET|OUTLET)"
            {
                type        escape;
            }
        );
    }
You have to choose between rebound, escape and stick!

Best,
Thomas
Dommy is offline   Reply With Quote

Old   July 19, 2013, 10:00
Default
  #3
New Member
 
Cay Myers
Join Date: Jun 2013
Posts: 14
Rep Power: 5
baila12 is on a distinguished road
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.
Attached Files
File Type: txt reactingCloud1Properties.txt (5.2 KB, 32 views)
baila12 is offline   Reply With Quote

Old   July 19, 2013, 11:11
Default
  #4
New Member
 
Thomas F.
Join Date: Jun 2013
Location: Graz, Austria
Posts: 14
Rep Power: 5
Dommy is on a distinguished road
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
Dommy is offline   Reply With Quote

Old   July 19, 2013, 11:42
Default
  #5
New Member
 
Cay Myers
Join Date: Jun 2013
Posts: 14
Rep Power: 5
baila12 is on a distinguished road
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!
baila12 is offline   Reply With Quote

Old   July 19, 2013, 16:48
Default
  #6
New Member
 
Thomas F.
Join Date: Jun 2013
Location: Graz, Austria
Posts: 14
Rep Power: 5
Dommy is on a distinguished road
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
Dommy is offline   Reply With Quote

Old   July 22, 2013, 08:15
Default
  #7
New Member
 
Cay Myers
Join Date: Jun 2013
Posts: 14
Rep Power: 5
baila12 is on a distinguished road
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 is offline   Reply With Quote

Old   July 31, 2013, 13:54
Default
  #8
New Member
 
Cay Myers
Join Date: Jun 2013
Posts: 14
Rep Power: 5
baila12 is on a distinguished road
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.
baila12 is offline   Reply With Quote

Reply

Tags
airfoil, ltsreactingparcelfoam, paraview

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with airfoil shape optimization robyTKD SU2 Shape Design 6 June 13, 2013 19:13
Unexplained Error during Solver Runs cfb CFX 6 November 9, 2012 16:42
Strange residuals of the Density Based Solver Pat84 FLUENT 0 October 22, 2012 15:59
Quarter Burner mesh with periosic condition SamCanuck FLUENT 2 August 31, 2011 11:34
Oscillating Airfoil - Dynamic meshing or user-defined velocity profile. DarrenC Main CFD Forum 5 July 19, 2010 22:33


All times are GMT -4. The time now is 22:48.