I noticed that in simulations
I noticed that in simulations with porous zones, the velocity oscillates at the exit of the porous zone. Moreover something unexpected happens also close to the inlet, as shown in the picture.
The system is a 2D channel, 1m long (diameter 0.07), with a porous zone in the first 0.5m (porosity = 0.5, d = 10^7, f = 0). The inlet velocity is 1m/s. All solver settings are taken from the rhoSimplePorousFoam implicit tutorial case.
The result shown here is obtained with a coarse grid. The behaviour improves significantly with a very grid with cells of some millimeter, but some issue is still present at the end of the porous zone, as shown below (cell 0.001 x 0.001).
Is this a known issue or are there particular requirements for the mesh around the porous zone?
With kind regards,
Sorry for the typo in "The
Sorry for the typo in
"The behaviour improves significantly with a very grid with cells..."
which should be:
"The behaviour improves significantly with a very fine grid with cells..."
Hi, here is just a guess, when
Hi, here is just a guess, when the density changes significantly for the inlet and outlet, a shock wave case probably occurs, to overcome this, some type of artificial viscosity should be added to the momentum equation, see Von Neumann and Richtmyer (1950) and Jameson et al
Thanks Roy, I agree because
I agree because the porous zone introduces a discontinuity in the resistance at the borders of the porous zone. However I think the right way to solve this is to work on the interpolation schemes/solution algorithms to manage the shock.
What gives me something to think about is the too high sensitivity to the mesh refinement. Even on a very refined mesh, the velocity profile changes significantly.
Grid with squared cells of 3mm:
Grid with squared cells of 1mm:
|All times are GMT -4. The time now is 10:26.|