CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Simulation flow through Propeller by SRFSimpleFoam (http://www.cfd-online.com/Forums/openfoam-solving/123182-simulation-flow-through-propeller-srfsimplefoam.html)

baoaero September 6, 2013 05:38

Simulation flow through Propeller by SRFSimpleFoam
 
Dear All
I have problem with simulation flow over propeller by OpenFoam. I used OpenFoam with solver SRFSimpleFoam. I set intinial conditions and parameters follow tutorial Mixer but i dont have result or value of P very big so i don what were the reasons?. I used turbulence model Komega SST. I dont know calculate force and torque of Propeller. Please help me solver this problem? I very confused...:confused:

Thank you very much .....


epsilon
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [ 0 2 -3 0 0 0 0 ];

internalField uniform 2.546;

boundaryField
{
inlet
{
type fixedValue;
value uniform 2.546;
}
outlet
{
type zeroGradient;
}
blade
{
type epsilonWallFunction;
U Urel;
value uniform 2.546;
}
cylinder
{
type zeroGradient;
}
K
}
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [ 0 2 -2 0 0 0 0 ];

internalField uniform 0.84375;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0.84375;
}
outlet
{
type zeroGradient;
}
blade
{
type kqRWallFunction;
value uniform 0.84375;
}
cylinder
{
type zeroGradient;
}

}

OMEGA
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [ 0 0 -1 0 0 0 0 ];

internalField uniform 3.014;

boundaryField
{
inlet
{
type fixedValue;
value uniform 3.014;
}
outlet
{
type zeroGradient;
}
blade
{
type omegaWallFunction;
U Urel;
value uniform 3.014;
}
cylinder
{
type zeroGradient;
}

}


// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [ 0 2 -2 0 0 0 0 ];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
blade
{
type zeroGradient;
}
cylinder
{
type zeroGradient;
}


}

UREL
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object Urel;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [ 0 1 -1 0 0 0 0 ];

internalField uniform ( 0 0 0 );

boundaryField
{
inlet
{
type SRFVelocity;
inletValue uniform ( 0 0 -15 );
relative no;
value uniform ( 0 0 0 );
}
outlet
{
type zeroGradient;
}
blade
{
type fixedValue;
value uniform ( 0 0 0 );
}
cylinder
{
type SRFVelocity;
inletValue uniform ( 0 0 0 );
relative no;
value uniform ( 0 0 0 );
}

}


// ************************************************** *********************** //************************************************** *********** //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object SRFProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

SRFModel rpm;

axis ( 0 0 -1 );

rpmCoeffs
{
rpm 5000;
}

sivakumar September 7, 2013 07:28

Hi,
Actually your questions are not clear, however I will try my best.
Have a look of the following thread for calculating initial value for k and epsilon.

http://www.cfd-online.com/Forums/ope...urbulence.html


To get the forces and torque, just copy and paste the following function underneath your controlDict file.

functions
(
forces
{
type forces;
functionObjectLibs ("libforces.so");
outputControl timeStep;
outputInterval 10;
patches (fan);
// pname p;
// Uname U;
rhoName rhoInf;
log true;
rhoInf 1.20;
CofR (0 0 0);
}
)


Thanks,
Sivakumar

baoaero September 8, 2013 11:39

Dear Sivakumar
I am grateful to you for your reply. I will describe my problem.
I have propeller with diameter 9 inchs and picth 6 inch (Master propeller 9x6 for UAV). I had done experiment in windtunnel with this propeller and determined their thrust and torque. Now, my problem is how determine thrust and torque of propeller by numerrical method using OpenFOAM. I have 3D model of propeller and generated mesh. Parameter of propeller: inlet velocity (speed flight): 10 - 15 m/s; RPM of Propeller 4000 - 8000 RPM; chord of propeller 0.75R = 17 mm.
I used tool estimate yplus of CFD online and calculate yplus. My mesh model has the first node from surface blade yplus = 20.
For simulation flow through this propeller, i used solver SRFSimpleFoam version 2.1.0 and turbulence model Komega SST (flow steady and incompressible). I set value of boundary conditions k, omega, nut, p, Urel as in last post (please see last post). I calculate 1000 time but I have result different pressure very large and force very large (> 10000 N). I think wrong because force and pressyre very large.

Would you mind help me solver this proplem?

1) With my problem, I used solver SRFSimpleFoam is right or wrong? Which is solver of OpenFoam i can calculate thrust and torque of propeller (steady and incompressible) ? If i use another solver of openfoam, do i change mesh mode of propeller?

2) How to calculate value of boundary condition of Problem? which is turbulence model i can choose for solver this problem?

3) Do you have tutorials OpenFoam relation simulation flow through propeller? please share me...?

Please help me...

Thank you very much.
Baoaero.

sivakumar September 9, 2013 04:17

Hi,
There are two options in OF to solve the problem, one is SRF and another one is MRF.

Getting torque and force from MRFSimpleFoam is very simple, just you need to copy and paste the function in your controlDict file. (the function is from the forum only), but I am not sure how to get torque and force from SRF solver, I dont think the function will work with SRF solver.

I can help you for MRFSimpleFoam solver, if you like give me your email id, I will send you the tutorial.

Regarding the turbulence model, its up to your need. If you are interested to calculate stall, the low-Re model is appropriate, then your Y+ value should be 1, otherwise you can use hi-Re model, which the Y+ range should be 30-300.


Thanks,
Sivakumar

baoaero September 9, 2013 04:40

Dear Sivakumar
Thank you for your advances. My email; baodangacv@yahoo.com or quocbao@vietnamairport.vn. I will use MRFSimplefoam for simulations propeller. I need calculate force and torque of propeller. How will i refine mesh model of propeller for using solver MRFSimplefoam?.
Please help me...Thanks you very much
Baoaero.

baoaero September 9, 2013 21:13

Dear Sivakumar
I received the file fan tutorial from you. Thank you for your sharing. Can you guide me for generating Mesh model of Propeller when i use MRFSimplefoam. How many calculated domain of Mesh model?
Please help me...
Thank you very much.
Baoaero

sivakumar September 10, 2013 03:33

Hi,
Have you gone through the tutorial? I have explained everything in that.
The problem is I dont know how you want to model your propeller, in my case I am using 3 volumes.
Inlet volume, fan volume , and outlet volume.

go through the tutorial once again then you will come to know.

Thanks,
Sivakumar

Porisel October 1, 2013 10:44

hi, i hope i can use this thread for my question. i want to simulate a kind of propeller, too. i desgined the propeller in solidworks and put it with snappyhexmesh in a mesh created up by blockmesh. the propeller should rotate in still water. so there is no inflow. my question ist, who to set the boundery settings? i user to solve the problem the SRFSimpleFOAM solver. i set the inlet:
Code:

  inlet
    {
    type  SRFVelocity;
    inletValue uniform (0 0 0);
    relative yes;
    value uniform (0 0 0);

and the outlet as type "slip".
the shape of the propeller is:
Code:

      type            fixedValue;
        value          uniform (0 0 0);

are these bounderysettings ok or are they false?

TReviol October 2, 2013 02:51

Hi,

"slip" for free surface and "fixedValue" for propeller shape are ok. And I never set an inlet condition, too.

Are the boundaries of your fluid domain open for in- and outflow? In this case you could use an "inletOutlet" condition. When your fluid domain is a vessel from process industry, you need neither an inlet nor an outlet. In this case all your boundaries are walls.

I hope I understood your question correctly and I could help.

Regards

Thomas

Porisel October 2, 2013 02:57

hi Thomas, thank you for your answer. the bounderies are patches, except the propeller-boundery, which is a wall.
what do you mean with "open for in- and outlet"? with fluid domain, you mean the initial conditions? the propeller should be in a "open ocean", no borders etc.

in german:
Hallo Thomas, danke für deine Antwort. Eigentlich hab ich nur patches, außer dem Propeller, welche eine wall ist.
Was meinst du mit "open for in- and outlet" und "fluid domain"? der propeller soll im offenen wasser sein, also keine grenzen oder hindernisse.

crusen mind June 24, 2016 06:53

Quote:

Originally Posted by sivakumar (Post 450841)
Hi,
Have you gone through the tutorial? I have explained everything in that.
The problem is I dont know how you want to model your propeller, in my case I am using 3 volumes.
Inlet volume, fan volume , and outlet volume.

go through the tutorial once again then you will come to know.

Thanks,
Sivakumar

can You send the tutorial to my email ID ?
my mail Id is velan45@gmail.com
thanks in advance


All times are GMT -4. The time now is 04:29.