Compression instead of expansion
Hi everybody, I'm simulating an axial turbine with mixing plane and real gas. As you know, a turbine expand the gas, but in my simulation there are zones in which the gas is higly compressed ( so non-physical results). Do you have any hint on how to solve this or how to control? Maybe, i've put wrong boundary condition. Here are the 0/ files:
alphat dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { Inflow_stator { type calculated; value uniform 0; } Outflow_stator { // type mixingPlaneFvPatchField: type mixingPlane; value $internalField; } Wall_stator { type alphatWallFunction; Prt 0.741; value uniform 0; } Perio_stator { type cyclic; value uniform 0; } Inflow_rotor { // type mixingPlaneFvPatchField: type mixingPlane; value $internalField; } Outflow_rotor { type calculated; value uniform 0; } Wall_rotor { type alphatWallFunction; Prt 0.741; value uniform 0; } Perio_rotor1 { type cyclic; value uniform 0; } Perio_rotor2 { type cyclic; value uniform 0; } } epsilon dimensions [0 2 -3 0 0 0 0]; internalField uniform 10000; boundaryField { Inflow_stator { type fixedValue; value uniform 2000; } Outflow_stator { type mixingPlane; value $internalField; } Wall_stator { type compressible::epsilonWallFunction; refValue uniform 0; value $internalField; Cmu 0.09; kappa 0.41; E 9.8; } Perio_stator { type cyclic; value $internalField; } Inflow_rotor { type mixingPlane; value $internalField; } Outflow_rotor { type inletOutlet; inletValue $internalField; value 13450; } Wall_rotor { type compressible::epsilonWallFunction; refValue uniform 0; value $internalField; Cmu 0.09; kappa 0.41; E 9.8; } Perio_rotor1 { type cyclic; value $internalField; } Perio_rotor2 { type cyclic; value $internalField; } } i dimensions [0 2 -2 0 0 0 0]; internalField uniform 27000; boundaryField { Inflow_stator { type fixedValue; value $internalField; } Outflow_stator { // type mixingPlaneEnthalpyJump; // patchType mixingPlane; type mixingPlane; rotating false; value $internalField; } Perio_stator { type cyclic; value $internalField; } Wall_stator { type gradientEnthalpy; gradient uniform 0; value uniform 463864; } Inflow_rotor { // type mixingPlaneEnthalpyJump; // patchType mixingPlane; type mixingPlane; rotating true; value $internalField; } Wall_rotor { type gradientEnthalpy; gradient uniform 0; value $internalField; } Outflow_rotor { type inletOutlet; inletValue $internalField; value $internalField; } Perio_rotor1 { type cyclic; value $internalField; } Perio_rotor2 { type cyclic; value $internalField; } } k dimensions [0 2 -2 0 0 0 0]; internalField uniform 5; boundaryField { Inflow_stator { type fixedValue; value 0.4; } Outflow_stator { type mixingPlane; value $internalField; } Wall_stator { type compressible::kqRWallFunction; value $internalField; } Perio_stator { type cyclic; value $internalField; } Inflow_rotor { type mixingPlane; value $internalField; } Outflow_rotor { type inletOutlet; inletValue $internalField; value $internalField; } Wall_rotor { type compressible::kqRWallFunction; value $internalField; } Perio_rotor1 { type cyclic; value $internalField; } Perio_rotor2 { type cyclic; value $internalField; } } p dimensions [1 -1 -2 0 0 0 0]; internalField uniform 400000; boundaryField { Inflow_stator { type totalPressure; phi phi; rho none; psi psi; U U; gamma 1.4; p0 uniform 5.1e5; value $internalField; } Outflow_stator { type mixingPlane; rotating false; value $internalField; } Wall_stator { type zeroGradient; } Perio_stator { type cyclic; value $internalField; } Inflow_rotor { type mixingPlane; rotating true; value $internalField; } Outflow_rotor { type fixedValue; value uniform 300000; } Wall_rotor { type zeroGradient; } Perio_rotor1 { type cyclic; value $internalField; } Perio_rotor2 { type cyclic; value $internalField; } } T dimensions [0 0 0 1 0 0 0]; internalField uniform 330; boundaryField { Inflow_stator { type totalTemperature; phi phi; rho rho; psi psi; U U; gamma 1.4; T0 uniform 338; value $internalField; } Outflow_stator { // type mixingPlaneEnthalpyJump; // patchType mixingPlane; type mixingPlane; rotating false; value $internalField; } Wall_stator { type zeroGradient; } Perio_stator { type cyclic; value $internalField; } Inflow_rotor { // type mixingPlaneEnthalpyJump; // patchType mixingPlane; type mixingPlane; rotating true; value $internalField; } Outflow_rotor { type zeroGradient; value $internalField; } Wall_rotor { type zeroGradient; } Perio_rotor1 { type cyclic; value $internalField; } Perio_rotor2 { type cyclic; value $internalField; } } U dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 100); boundaryField { Inflow_stator { type pressureDirectedInletVelocity; inletDirection uniform (0 0 1); value uniform (0 0 0); } Outflow_stator { type mixingPlane; rotating false; value $internalField; } Wall_stator { type fixedValue; value uniform (0 0 0); } Perio_stator { type cyclic; value $internalField; } Inflow_rotor { type mixingPlane; rotating true; value $internalField; } Outflow_rotor { type pressureInletOutletVelocity; inletValue $internalField; value $internalField; } Wall_rotor { type fixedValue; value uniform (0 0 0); } Perio_rotor1 { type cyclic; value $internalField; } Perio_rotor2 { type cyclic; value $internalField; } } These are the results after 5000 iterations: Flux divergence min = 4.28614548132e-05 max = 24180860.2673 average: 15406.0985031 Writing divFlux field Mach number min = 0.00198878622956 max = 4.35597144753 Writing Mach number field Field U magnitude min = 0 max = 1000 Field p min = 49202.4992467 max = 1055827.15616 Field rho min = 1.52104950445 max = 27.9890057068 Field T min = 337.498873263 max = 1163.72657359 Mixing plane pair (Outflow_stator, Inflow_rotor) : 0.756677749244 -0.757134457367 Diff = -0.00045670812309 or 0.0603570177062 % End max U should be around 160 m/s and max T and p should be around the ones at the inlet. Any hint would be appreciated. Thanks in advance. |
Hi, please edit your post and use code tags. No one will go through a not formatted text, at least I donīt
|
All times are GMT -4. The time now is 10:13. |