|
[Sponsors] |
Density Modelling Problems in Compressible Case |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 21, 2016, 10:17 |
Density Modelling Problems in Compressible Case
|
#1 |
Member
Charles K.
Join Date: Feb 2016
Location: Germany
Posts: 37
Rep Power: 10 |
Hi Foamers,
I am trying to simulate a compressible case with rhoPimpleFoam and rhoSimpleFoam. The problem is that the density isn't calculated correctly as it is almost constant trough the entire domain except for the INLET and the 2 OUTLETS, where I configured a specific density as I am using flowRateInletVelocity (mass flow) as Inlet/outlet BC. There should be strong density variations in certain areas of the Model, e.g. the flow is accelerated through a nozzle. There is probably something wrong with my BCs. I would be very thankful if someone could take a look at my settings below. Let me know if you need more specific informations. Cheers! Charles U Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { PLENUM { type fixedValue; value uniform (0 0 0); } ROTOR { type rotatingWallVelocity; origin (0 0 0); axis (1 0 0); omega constant 838; value uniform (0 0 0); } STATOR { type fixedValue; value uniform (0 0 0); } OUTLET_LABY { type flowRateInletVelocity; massFlowRate constant -0.00031042; rhoInlet 2.776162; } INLET { type flowRateInletVelocity; massFlowRate constant 0.002989; rhoInlet 3.09; value uniform (0 0 0); } OUTLET { type pressureInletOutletVelocity; value uniform (0 0 0); } CUT { type slip; } SPALT_PERIODIC_1 { type cyclicAMI; } SPALT_PERIODIC_2 { type cyclicAMI; } PLENUM_PERIODIC_1 { type cyclicAMI; } PLENUM_PERIODIC_2 { type cyclicAMI; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 300000; boundaryField { PLENUM { type zeroGradient; } ROTOR { type zeroGradient; } STATOR { type zeroGradient; } OUTLET_LABY { type zeroGradient; } INLET { type zeroGradient; } OUTLET { type fixedValue; value uniform 250304; } CUT { type zeroGradient; } SPALT_PERIODIC_1 { type cyclicAMI; } SPALT_PERIODIC_2 { type cyclicAMI; } PLENUM_PERIODIC_1 { type cyclicAMI; } PLENUM_PERIODIC_2 { type cyclicAMI; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 390; boundaryField { PLENUM { type zeroGradient; } ROTOR { type zeroGradient; } STATOR { type fixedValue; value uniform 389; } OUTLET_LABY { type inletOutlet; inletValue uniform 369.989; value uniform 369.989; } INLET { type fixedValue; value uniform 400; } OUTLET { type inletOutlet; inletValue uniform 390; value uniform 390; } CUT { type zeroGradient; } SPALT_PERIODIC_1 { type cyclicAMI; } SPALT_PERIODIC_2 { type cyclicAMI; } PLENUM_PERIODIC_1 { type cyclicAMI; } PLENUM_PERIODIC_2 { type cyclicAMI; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object thermophysicalProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // thermoType { type hePsiThermo; mixture pureMixture; transport sutherland; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } mixture { specie { nMoles 1; //no combustion = 1 molWeight 28.9; // für trockene Luft in g/mol } thermodynamics { Cp 1007; Hf 2.544e+06; } transport { As 1.4584e-06; Ts 110.33; } } // ************************************************************************* // |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Setting the height of the stream in the free channel | kevinmccartin | CFX | 12 | October 13, 2022 21:43 |
compressible pipe flow with varying density | nandiganavishal | OpenFOAM | 0 | July 23, 2011 19:20 |
The Density of Compressible Flow | panou | Main CFD Forum | 0 | April 27, 2011 04:58 |
Free surface boudary conditions with SOLA-VOF | Fan | Main CFD Forum | 10 | September 9, 2006 12:24 |
2-D Compressible Axisymmetric modelling | jehazneb | FLUENT | 0 | December 20, 2004 05:12 |