|
[Sponsors] |
July 21, 2015, 04:32 |
rhoPimpleFoam for Inviscid Simulation
|
#1 |
Senior Member
KGN
Join Date: Oct 2009
Location: Chennai, India
Posts: 121
Rep Power: 16 |
I am trying to simulate a Combustion Chamber (only flow simulation; no combustion) with 10 outlets.
I wrote a custom BC using groovyBC for inlet (which is burn surface) - Mass generation which changes with time. 10 Outlets - Mass Outflow profile. All the walls are inviscid. First I am trying to get a steady state solution of certain chamber pressure - Mass generation is equal to net mass outflow of 10 outlets.. Since the flow velocity is low subsonic (<0.2), i am using rhoPimpleFoam. As the simulation progresses, the rhoMax keeps on increases and rhoMin keeps in decreasing and also the velocity. The simulation blows saying Maximum number of iteration exceeded. Whether rhoPimpleFoam can be used to simulated Inviscid simulation? I have attached the boundary condition and solver setup Thanks KGN fvSchemes.txt fvSolution.txt 0.zip Last edited by mecbe2002; July 21, 2015 at 06:16. |
|
July 21, 2015, 09:24 |
|
#2 |
Senior Member
KGN
Join Date: Oct 2009
Location: Chennai, India
Posts: 121
Rep Power: 16 |
1. I am using Mass flow Rate boundary condition at Inlet and Outlet. Is it a problem?
2. Whether rhoPimpleFoam can solve Inviscid Simulation? 3. I am able to solve the simulation using rhoCentralFoam, but it takes a long time. Max Co i can use is 0.55. Is there a way to increase the Co number? Also rhoCentralFoam is true transient solver? I guess the problem is with the Boundary condition. Both at inlet and outlet mass flow rate is specified. For pressure zero gradient is used. So this is creating the problem? time step continuity errors : sum local = 6.16642e-08, global = -6.09083e-08, cumulative = -4.11783e-08 rho max/min : 25 8.04696 PIMPLE: iteration 2 smoothSolver: Solving for Ux, Initial residual = 0.000154991, Final residual = 8.7961e-08, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.000494649, Final residual = 4.37157e-07, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.00055955, Final residual = 1.41411e-06, No Iterations 2 smoothSolver: Solving for e, Initial residual = 0.0107836, Final residual = 3.66265e-07, No Iterations 2 [0] [0] [0] --> FOAM FATAL ERROR: [0] Maximum number of iterations exceeded [0] [0] From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const [0] in file /home/rcdt/OpenFOAM/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76. [0] FOAM parallel run aborting [0] [0] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/rcdt/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #1 Foam::error::abort() in "/home/rcdt/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #2 Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double) const) const in "/home/rcdt/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" [0] #3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/home/rcdt/OpenFOAM/OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" Last edited by mecbe2002; July 22, 2015 at 11:41. |
|
July 21, 2015, 22:15 |
|
#3 |
Senior Member
KGN
Join Date: Oct 2009
Location: Chennai, India
Posts: 121
Rep Power: 16 |
Is there any method to set relaxation factor and relaxation tolerance ?
Last edited by mecbe2002; July 22, 2015 at 11:42. |
|
July 22, 2015, 11:45 |
|
#4 |
Senior Member
KGN
Join Date: Oct 2009
Location: Chennai, India
Posts: 121
Rep Power: 16 |
Any reply?
I tried with different settings of relaxation factor and relaxation residual, but still the simulation blows after few iterations. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Waterbell Simulation | o.b.m | Fluent Multiphase | 1 | November 5, 2015 05:55 |
Simulation of high pressure diesel injector - all phases compressible with cavitation | fivos | CFX | 4 | July 30, 2015 06:48 |
mass flow rate issue in supersonic nozzle simulation | xkang | FLUENT | 0 | July 31, 2014 16:06 |
Simulation of a complex wing in solidworks flow simulation | niels1900 | FloEFD, FloWorks & FloTHERM | 6 | April 20, 2011 10:44 |
Continuous vs interrupted simulation | sega | OpenFOAM Running, Solving & CFD | 4 | November 3, 2008 14:29 |