|March 1, 2012, 11:07||
Rotor-Stator interface boundary condition
Join Date: Feb 2011
Posts: 10Rep Power: 8
I don't know if this is the correct place to solve implementation issues, but let's try it
1) The way it is implemented in Edge now is that the interface boundary averaging is done in the radial direction (at least there must be some nontrivial variation in the radial component of the boundary coordinates).
Do you think that it would be useful to have the possibility of having the interface boundary averaging done in axial direction (i.e., boundaries with constant radial coordinate)? It is not difficult and I can do it - just want to know if it is useful for somebody.
2) I am almost sure (at least have some doubts) that conservation as it is implemented in econsbc_m.f90 and gcorrotstat_m.f90 is not working properly in some cases, at least if rotor and stator are rotational periodic (e.g., one blade modeled).
I have run some tests and the results were different for different relative positions of rotor and stator (which corresponds for instance to the fact that I pick different period). I used the trunk version but used the area weighting instead of mass flux weighting in econsbc_m.f90.
If I force FLSUM(2:4)=0 (don't use conservation for momentum equations) the results were ok - the same radial distribution of variables for both configurations.
I think that
FLSUM(1) = Integral RO * REL.VELOCITY DOT NORMAL
is independent of the position (this number remains the same if you rotate the mesh and velocities accordingly), and so is the conservation term for energy equation and probably also for turbulence (at least k, omega)
FLSUM(2) = Integral RO * REL.VELOCITY DOT NORMAL * U + P*N1
FLSUM(3) = Integral RO * REL.VELOCITY DOT NORMAL * V + P*N2
FLSUM(4) = Integral RO * REL.VELOCITY DOT NORMAL * W + P*N3
Somebody having problem with this? Do you think this could be some bug?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Domain Imbalance||HMR||CFX||5||October 10, 2016 05:57|
|Can anyone give me some hint on how to make traction free boundary condition?||poplar||OpenFOAM||3||January 14, 2015 03:37|
|Boundary Conditions||Thomas P. Abraham||Main CFD Forum||20||July 7, 2013 05:05|
|Setting outlet Pressure boundary condition using CAFFA code||Mukund Pondkule||Main CFD Forum||0||March 16, 2011 04:23|
|How to set boundary condition in Fluent for the fo||Peiyong||FLUENT||1||November 10, 2006 12:44|