CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   How to model a fan fixing the mass flow rate? (https://www.cfd-online.com/Forums/openfoam-solving/80628-how-model-fan-fixing-mass-flow-rate.html)

herbert October 5, 2010 08:14

Hi Mad,

sorry it was my fault. The definiton in the boundary file is correct, but in your U file the BC-type must be directMapped (not directMappedFixedValue).

Sorry again,
Stefan

maddalena October 5, 2010 08:47

Quote:

Originally Posted by herbert (Post 277868)
sorry it was my fault. The definition in the boundary file is correct, but in your U file the BC-type must be directMapped (not directMappedFixedValue).

Ah yes! I should have noticed by myself! :o There is an example among the tutorials...
Just as an information: do you know why the BC condition is directMapped and in the /OpenFOAM/OpenFOAM-1.6.x/src/finiteVolume/fields/fvPatchFields/derived folder there are:
  • directMappedFixedValue and
  • directMappedVelocityFluxFixedValue?
What is the difference between the two?

Cheers

mad

herbert October 5, 2010 09:26

I think directMappedVelocityFluxFixedValue recyles the velocity and the flux whereas directMappedFixedValue (which we are using here) only patches velocity in this case.

The biggest disadvantage of the first one if you try to apply it to your case is that you can't set the average value there. My guess is, that you could need this BC for compressible simulations where phi is not only a function of U.

Regards,
Stefan

maddalena October 5, 2010 09:32

Quote:

Originally Posted by herbert (Post 277880)
The biggest disadvantage of the first one if you try to apply it to your case is that you can't set the average value there.

Yes, I saw that is not possible to fix the average value when using the directMappedVelocityFluxFixedValue. This means that I cannot fix the mass flux in my case.
I have the first simulation running. I will report as soon as I have some results.
Thanks for your support.

cheers

mad

robingilbert October 5, 2010 16:49

hello,

I have a heat source (constant heat addition using setFields) very next to the fan. i dont know the temperature at the fan. i just have a heat source befor it. how do i set up the temperature at the fan boundary?? do i use directMapped as well??

so its like:
fanOutlet=domain inlet
fanOutlet
{
type directMapped;
value uniform 0;
setAverage false;
}
what abt fanInlet??

maddalena October 6, 2010 03:45

Hi Robin,
are you using cht? If it is so,
Quote:

Originally Posted by robingilbert (Post 277952)
how do i set up the temperature at the fan boundary?

I think that you can use the fan as an interface coupling, indeed this is what you do when applying the directMapped. Why do not create the two fan faces in two different regions, and then couple them using a solidWallMixedTemperatureCoupled applied to T, while U and p are as described above.
Quote:

Originally Posted by robingilbert (Post 277952)
I have a heat source (constant heat addition using setFields) very next to the fan. i dont know the temperature at the fan. i just have a heat source before it.

maybe you can find this thread interesting: http://www.cfd-online.com/Forums/ope...egionfoam.html I believe the discussion can be applied to a heat source in general (even defined as a setField), since we traced back the problem to the coupling condition. But perhaps you will be luckier than me... :rolleyes:

mad

robingilbert October 6, 2010 04:18

i am using simpleFoam.

maddalena October 6, 2010 09:18

1 Attachment(s)
Hello,
bad news: the following caused velocity instabilities (i guess due to the numerics) close to the fanOutlet,
  • fanOutlet = domain inlet
    • U: directMapped, setAverage: massFlowRate/Area
    • p: zeroGradient
  • fanInlet = domain outlet
    • U: zeroGradient
    • p: fixedValue p = 0
You can see the velocity field in the attached picture, while there is a brief explanation of my geometry in this post. the flow goes from left to right.
Note that at the moment turbulence is switched off to not include k and epsilon.
I used the fixedValue p = 0 BC since I usually fix velocity in one of the boundary and pressure in one of the other. However, Stefan suggested here a zeroGradient + pRefCell close to the outlet. How can I get its number? where should I check it? EDIT: answer here: http://www.cfd-online.com/Forums/ope...-pressure.html

Is it correct to not fix pressure anywhere and use zeroGradient only? May the pressure fixedValue be the reason of these instabilities?

cheers

mad

maddalena October 7, 2010 03:41

No way. Simulation crashed after few time steps using the zeroGradient condition at the inlet and at the outlet. Ideas?

mad

eugene October 8, 2010 05:47

If you have switched off turbulence, what is your Re and viscosity? If you do not have a reasonably high nu and Re is high as well, the solution will of course blow up (especially if you are using 2nd order numerics).

The flow is behaving very strangely downstream of the fan. Do you have any idea why this is? I suggest you look at the velocity values at the inlet in a bit more detail.

maddalena October 8, 2010 06:26

2 Attachment(s)
Hi Eugene,
I used a laminar model only on the first 100 time steps, just to help the solver. Then I switched on the turbulence.
My schemes are the following:
  • gradSchemes faceMDLimited Gauss linear 0.5;
  • divSchemes Gauss linearUpwindV cellLimited Gauss linear 1;
  • laplacianSchemes Gauss linear limited 0.5;
and the strange behaviour of the outlet is due to the BC suggested above (k and epsilon fixedValue at the fanOutlet patch, I guess). Now I changed them in the following way:
  • fanOutlet = domain inlet:
    • U: directMapped, average massflow/Area
    • p: zeroGradient;
    • k, epsilon: inletOutlet
  • fanInlet = domain outlet:
    • U: zeroGradient;
    • p: fixedValue 0; pRefCell close to the outlet;
    • k, epsilon: inletOutlet
As you can see, the velocity field looks better now, even close to the fan.
In order to make the simulation more stable, I also lowered relaxation factor:
  • p 0.15;
  • U 0.5;
  • k 0.3;
  • epsilon 0.3;
This worked. After 550 time steps the velocity field does not shows any particular problems. However, pressure residual convergence is not nice, see the attached file.
What parameters should I change to tune my simulation?

mad

eugene October 10, 2010 19:46

Hi Maddalena,

I'm not sure how to proceed next. You numerics looks stable. You could try mapping your turbulence properties as well and use pure upwind convection for turbulence. I guess you are having stability problems because there is a discontinuity in your pressure field over the fan while your velocity field is continuous. You could try reverse mapping your pressure field with a prescribed mean, but this might just make things worse. You could also try zero gradient p at the outlet, but again this will not necessarily improve stability. There is a "fixedMeanValue" boundary type in OpenFOAM-Extend that might be a better fit for your problem.

maddalena October 11, 2010 03:50

Hi Eugene and thank you for your answer.
Quote:

Originally Posted by eugene (Post 278602)
I guess you are having stability problems because there is a discontinuity in your pressure field over the fan while your velocity field is continuous.

Yes, that sounds a logical explanation. But this is a fan indeed!
Quote:

Originally Posted by eugene (Post 278602)
You could also try zero gradient p at the outlet, but again this will not necessarily improve stability.

Already tried that (some posts above, but it is "hidden"). You are right, stability does not improve.
Quote:

Originally Posted by eugene (Post 278602)
There is a "fixedMeanValue" boundary type in OpenFOAM-Extend that might be a better fit for your problem.

Ok, this will be my following step. I will try this solution as soon as my processors will be free. Thanks once again.
Cheers

mad

maddalena October 12, 2010 10:45

Quote:

Originally Posted by eugene (Post 278602)
There is a "fixedMeanValue" boundary type in OpenFOAM-Extend that might be a better fit for your problem.

Tried that, without success... :(

NicolasB June 27, 2011 11:35

Hi Maddalena,

Sorry for digging out this old post.

I'd like to set the mass flow rate on a fan. But I've a problem with my mesh, and I think you could help me.
I made my mesh with Gambit, setting the fan as "fan". I import it with the "fluent3DMeshToFoam" utility. The fan become a cyclic patch, and when I give the next BC:
type flowRateInletVelocity;
patchType cyclic;
flowRate 8.79;
value uniform (2 0 0);
it's applied on both sides of the fan... with external normal :(
Could you, please, teach me how to use the mesh utilities in order to have a "fan_inlet" and a "fan_outlet"? unless there is an other solution, like giving a vector.
I've successfully tried to give a velocity, but that's not what I need.

Regards,
Nicolas

maddalena June 28, 2011 04:06

Hi Nicolas,
of course this will not work:
Quote:

Originally Posted by NicolasB (Post 313740)
type flowRateInletVelocity;
patchType cyclic;
flowRate 8.79;
value uniform (2 0 0);

You cannot define something as cyclic and than apply to it a standard BC! Cyclic wants something that requires cyclic, i.e. cyclic or fan. Not something else.

As for concluding my discussion above, I used a modified simpleFoam version in such a way that some properties may be given to specific cellSets. robingilbert made this solver, so you should ask to him to have a copy of it.

mad

NicolasB June 30, 2011 09:09

Hi Maddalena,

Thank you very much for your answer.
I gave up this option and dug a bit more the standard "fan" BC. In fact I had a problem since I regrouped 3 fans in the same patch. I fixed it by splitting my patch and setting up the pressure jump with the right sign for each fan.

Now this is working, I'd like to complicate.
I described my case in this thread, in which you wrote a lot :)

Regards,
Nicolas

dhruv October 27, 2011 06:10

same problem
 
Hi Maddalena,

I am trying to simulate a fan in a channel, and I am having the same problem. The pressure at the inlet of the fan is higher than the pressure at the outlet. How did you correct it?

Regards,

Dhruv.

Quote:

Originally Posted by maddalena (Post 277614)
Hi,
Just to say I have solved this:

and now pressure has the expected distribution. However, the simulation is not stable and residual never get low. My fvScheme and fvSolution should not be the problem: they are tuned for similar cases. I guess the BC are not set correctly to simulate what I want to solve. Is there anyone that can help me to understand what my errors can be?
Regards

maddalena


maddalena October 27, 2011 08:41

Hello,
Quote:

Originally Posted by dhruv (Post 329670)
I am having the same problem. The pressure at the inlet of the fan is higher than the pressure at the outlet.

what is the same problem? unstable simulation or flow going on the opposite direction you want to?
What OF version are you using? starting from 2.0.0 there is a new bc to simulate fan properly: http://www.cfd-online.com/Forums/ope...essure-bc.html

mad

dhruv November 2, 2011 04:41

flow in opposite direction
 
Hi Maddalena,

thanks for the reply. The problem was flow going in opposite direction, which I have solved now, by changing the order for the patches in the boundary file.

Regards,
Dhruv

Quote:

Originally Posted by maddalena (Post 329704)
Hello,
what is the same problem? unstable simulation or flow going on the opposite direction you want to?
What OF version are you using? starting from 2.0.0 there is a new bc to simulate fan properly: http://www.cfd-online.com/Forums/ope...essure-bc.html

mad



All times are GMT -4. The time now is 16:42.