CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

BC for counter-current sieve-plate separation column

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 18, 2015, 03:15
Default BC for counter-current sieve-plate separation column
  #1
New Member
 
Gian
Join Date: Feb 2015
Posts: 5
Rep Power: 4
vishmehr24 is on a distinguished road
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.
Attached Images
File Type: jpg init.jpg (14.0 KB, 6 views)
File Type: jpg alpha-3s.jpg (11.4 KB, 6 views)
vishmehr24 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Current density visualisation (PEM fuel cell add-on module) pchoopanya FLUENT 6 July 14, 2016 15:26
[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 10:29
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


All times are GMT -4. The time now is 06:17.