ActuatorDiskExplicitForce in OF2.1. Help
Are there a guide to now how to change a solver code prepared for OF1.5dev to a newer OF versions(OF2.1)?
I have been checking the solver windSimpleFoam of the version 2.0.1 and due to the fact than it also takes into account the thrust to modify the momentum equation.
I investigate and I found the following report:
It seems very interesting because it can be modelled the trust and the swirl with a radial distribution based on a analytical approach. It can be fixed the thrust and the torque of your wind turbine.
Although it was done based on OF 1.5dev, I have been trying to implement the same approach for OF 2.1 but I am not able to do it. Are there a guide to now how to change a solver between 1.5dev a newer OF versions?
When I try to compile, I receive the following error:
"readSIMPLEControls.H","initConvergenceCheck.H " and "convergenceCheck.H"
This is my first/second approach to the code of OF, so I will try to continue with the problem and If a find the mistake, I will post the solution.
With this new Ueqn.H code, I have received the next error:
I have compiled the new code, but....
But it seems that it does not work properly.
It computes the forces (axial and tangential), but the volumeForce is always zero. There is not a source.
Please, I think that I am near to solve the problem but I do not find where it is.
i am very interested in this case of actuator dis of Erik Svenning
but i found the same problem of you :( ???
are you solved??
please contact me
hello drsattar and be_inspired,
about the Erik Svenning code there's a couple of things you ahve to take into account,
first that code was made for OpenFoam 1.5 dev, and the coding the solvers has changed from that version to 2.1.x. so that's why the simplefoam solver that comes with the code won't compile in more recent versions of openfoam. if you insist on using it you would have to cut and paste the parts of the code into the newer version of simplefoam (it can be a bit tricky - it took me a while)
second, that code was made on a openfoam workshop, and not for a phd or masters thesis, i've personally updated the code and made it run in the new version of openfoam and i've encounter a couple of mistakes in the code:
- first in the code he sets the axial and tangential force to be in the same direction, which is not true, even though openfoam works in cartesian coordiantes you have to decompose the tagential force in Z and Y.
-also the way the code selects the cells to insert the bodyforce form points and faces leads to a erratic distribution of the force.
i finally gave up on that code :s
i would refer you to this thread:
here they show a very simple and practical way to generate and insert a bodyforce into a solver :),
and depending on you C++ it turns out to be very flexible to insert mathematical expresions to calculate the force.
i hope this helps and good luck!
i am very thanful for your reply, if you like send me your email to send you the final case which i created it of actuator disk model. i wish if you able to check it...
also the case of Fan (body force mode) is not work with me??? :o
my email (email@example.com)
thank again :)
it's not necessary for you to send me an e-mail :),
my recomendation is not to use the actuador disk model form Erik Svenning, as i told you it has mistakes that i was not able to correct,
you should go with the bodyforce model, i'm going to attach a sample solver and a sample case, (it's only set to add an axial force, you're gonna have to modify it to include the tangential force),
to make a good case,
the first thing you should do is to validate your code against the analytical solution, here's and article that talks about that, and gives the formulas for the anlytical solution,
after you can set up an actual real case and compare it against experimental data, or against other already validated model.
i wish i could help you more, but i also have do my own simulations, i have 2 semesters before i finish my masters :) and i have to concentrate on that.
if you have any other questions ask, but it might take me a while to answer,
the best of luck on your research, i know it's not easy
actuator disk (bodyForce)
I am very grateful to you about your answer on my problem, now i will try to use and modify the attachment files ......:)
i wish to be contact with me in the future.....
With your solver, pressure wiggles are not fixed:
This is quite quite important for me.
The report that you have linked is my reference point but I do not know how to modify simpleFoam and rotorDiskSource to manage body forces without pressure wiggles.
The problem happens when body force is directly applied over the UEqn and velocity and pressure are not coupled.
It seems like you have not implemented the correction proposed by your link
Do you have a modified version of your simpleFoamForce to solve that?
|All times are GMT -4. The time now is 11:49.|