|
[Sponsors] |
BC for counter-current sieve-plate separation column |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 18, 2015, 03:15 |
BC for counter-current sieve-plate separation column
|
#1 |
New Member
Gian
Join Date: Feb 2015
Posts: 6
Rep Power: 11 |
I am trying to run a simulation of a countercurrent sieve-plate liquid-liquid separation column through twoPhaseEuler but the simulations keeps on crashing with the following message:
--> FOAM FATAL ERROR: Maximum number of iterations exceeded 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 in file /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76. FOAM aborting #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::error::abort() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 Foam::species::thermo<Foam::hConstThermo<Foam:er fectFluid<Foam::specie> >, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectFluid<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectFluid<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:e rfectFluid<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double) const) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #3 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectFluid<F oam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #4 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectFluid<F oam::specie> >, Foam::sensibleEnthalpy> > > >::correct() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so" #5 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam" #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/twoPhaseEulerFoam" Aborted (core dumped) ------------------------------------------------------------- Physical System: Water is entering through a topInlet and should exit through the bottomOutlet. Toulene is entering through bottomInlet and should exit through topOutlet. thecolumn is filled Initiallywith waterThere is also a pressureOutlet at the top of thecolumn __________________________________________________ _- BC 0/ alpha.toulene: internalField uniform 0; boundaryField { bottomInlet { type fixedValue; value uniform 1; } topInlet { type fixedValue; value uniform 0; } topOutlet { type inletOutlet; phi phi.toulene; inletValue uniform 1; value uniform 1; } bottomOutlet { type inletOutlet; phi phi.toulene; inletValue uniform 0; value uniform 0; } walls { type zeroGradient; } pressureOutlet_1 { type inletOutlet; phi phi.air; inletValue uniform 1; value uniform 1; } defaultFaces { type empty; } } _________________________________________________- 0/p: internalField uniform 1e5; boundaryField { bottomInlet { type fixedFluxPressure; value $internalField; } topInlet { type fixedFluxPressure; value $internalField; } topOutlet { type fixedFluxPressure; value $internalField; } bottomOutlet { type fixedFluxPressure; value $internalField; } walls { type fixedFluxPressure; value $internalField; } pressureOutlet_1 { type fixedValue; value uniform 1e5; } defaultFaces { type empty; } } __________________________________________________ ______ 0/U.toulene: internalField uniform (0 0 0); boundaryField { bottomInlet { type fixedValue; value uniform (0.0175 0 0); } topInlet { type fixedValue; value uniform (0 0 0); } topOutlet { type fixedValue; value uniform (0.0175 0 0); } bottomOutlet { // type zeroGradient; type fixedValue; value uniform (0 0 0); } walls { type fixedValue; value uniform (0 0 0); } pressureOutlet_1 { type pressureInletOutletVelocity; phi phi.air; value uniform (0 0 0); } defaultFaces { type empty; } } --------------------------------------------------------------------------------------- 0/U.water internalField uniform (0 0 0); boundaryField { bottomInlet { type fixedValue; value uniform (0 0 0); } topInlet { type fixedValue; value uniform (0.0175 0 0); } topOutlet { type fixedValue; value uniform (0 0 0); } bottomOutlet { type fixedValue; value uniform (0.0175 0 0); } walls { type fixedValue; value uniform (0 0 0); } pressureOutlet_1 { type pressureInletOutletVelocity; phi phi.water; value uniform (0 0 0); } defaultFaces { type empty; } } ----------------------------------------------------- The initial (water-filled) configuration is shown as attached. On left lower side, toulene bottomInlet is seen with alpha.air=1 ALso a picture of the alpha,air (=alpha.toulene) at 3 s is attached. So, are my BC's reasonable?. Flow is taken to be laminar. Thanks for any answers/suggestions. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Current density visualisation (PEM fuel cell add-on module) | pchoopanya | FLUENT | 10 | August 21, 2023 15:33 |
[ICEM] Meshing/Blocking a sieve tray of a distillation column | Vdelmonego | ANSYS Meshing & Geometry | 1 | February 23, 2016 04:23 |
perforated plate boundary - Air injector bubble column | Jonas Ansoni | OpenFOAM Pre-Processing | 2 | August 23, 2014 11:29 |
[General] 2 datas on one plot | Akuji | ParaView | 46 | December 1, 2013 15:06 |
heat transfer on dimpled plate | mech | FLUENT | 4 | February 6, 2007 16:15 |