modelMulticomponentDistillationFUG

Extends from DistillationBase.

Information

It is a non standard implementation of the Fenske-Underwood-Gilliland shortcut method for multicomponent distillation.
The following simplifications are made:
- The distribution of components between the distillate and bottom, obtained by the Fenske equation for total reflux, is considered valid for the actual and minimum reflux.
- It is considered just one pinch point, close to the feed composition. As the distillate composition is taken from the Fenske's equation, only one Underwood set of equations is needed in order to calculate the minimum reflux. This may give poor calculation if there is an important distribution of any component between distillate and bottom.
The feed must be specified as flow, in kmol/h or kg/h. The top light, and the bottom heavy, can be specified as flow or concentration. Flow is recommended for best convergence.
The model is lacking one equation. You can configure the needed equation by selecting the required parameter. It is done at the graphical menu under the "Operating conditions" tab. You can select one of the options as true: fix the reflux to a constant value, or as a multiplier of the minimum reflux. Or fix the heat exchange at the condenser or the reboiler. Or you can add the extra equation (for example the value of the vpor flow at the stripping section) by hand.  


Parameters

TypeNameDefaultDescription
Modelica.Units.SI.MolarMassmwL (from DistillationBalance)0light component molecular weight
Modelica.Units.SI.MolarMassmwH (from DistillationBalance)0heavy component molecular weight
BooleanmolarBasis (from DistillationBalance)trueif true, use molar fraction and kmol/h, otherwise mass fraction and kg/h
Modelica.Units.SI.MolarMass[nNonKey]mwIfill(0, nNonKey)nonkey components molecular weight
Feed data
BooleanisBinaryDistillation (from DistillationBalance)trueif true only two components are considered
Realffl (from DistillationBalance)fixed feed flow rate of light component in kmol/h or kg/h
Realffh (from DistillationBalance)fixed feed flow rate of heavy component in kmol/h or kg/h
Modelica.Units.SI.TemperatureDifferencefSuperHeating (from DistillationBalance)0feed superheating, subcooling if negative
FractionfGfSat (from DistillationBalance)0gas fraction of the feed, if in saturated conditions
Modelica.Units.SI.MolarEnthalpyhvF (from DistillationBalance)41000feed molar vaporization enthalpy
Modelica.Units.SI.MolarHeatCapacitycpFl (from DistillationBalance)0feed liquid heat capacity, if subcooled
Modelica.Units.SI.MolarHeatCapacitycpFv (from DistillationBalance)0feed vapor heat capacity, if superheated
StringfLight""key light component name
StringfHeavy""key heavy component name
StringfNonKeys""feed nonkey components
IntegernNonKey1number of nonkey components
Real[nNonKey]ffifill(0, nNonKey)fixed feed flow of nonkey components in kmol/h or kg/h
Top data
FreeFluids.Types.PerformanceSpecificationtopFixed (from DistillationBalance)FreeFluids.Types.PerformanceSpecification.FixedCompositioncomposition, flow, recovery or none are fixed
Fractiondyl (from DistillationBalance)0.0fixed distillate light fraction
Realdfl (from DistillationBalance)0.0fixed distillate light flow. kmol/h or kg/h
Fractiondlr (from DistillationBalance)0.95fixed light recovery in distillate
Modelica.Units.SI.TemperatureDifferencedSubCooling (from DistillationBalance)0distillate subcooling, if total condenser used
Modelica.Units.SI.MolarEnthalpyhvD (from DistillationBalance)41000distillate molar vaporization enthalpy
Modelica.Units.SI.MolarHeatCapacitycpDl (from DistillationBalance)0distillate liquid heat capacity
Bottom data
FreeFluids.Types.PerformanceSpecificationbottomFixed (from DistillationBalance)FreeFluids.Types.PerformanceSpecification.FixedCompositioncomposition, flow, recovery or none are fixed
Fractionbxh (from DistillationBalance)0.0fixed bottom heavy fraction
Realbfh (from DistillationBalance)0.0fixed bottom heavy flow. kmol/h or kg/h
Fractionbrh (from DistillationBalance)0.95fixed heavy recovery in bottom
Modelica.Units.SI.MolarEnthalpyhvB (from DistillationBalance)41000bottom molar vaporization enthalpy
Operating conditions
BooleanfixRatioToMinReflux (from DistillationBalance)falseif true, the ratio reflux/minimum reflux is fixed
RealminRratio (from DistillationBalance)1.3ratio used, if previous is true
BooleanfixRefluxRatio (from DistillationBalance)falseif true, the ratio reflux/distillate is fixed
Realrr (from DistillationBalance)0fixed internal reflux/distillate ratio
BooleanfixCondenserPower (from DistillationBalance)false
Modelica.Units.SI.Powerwc (from DistillationBalance)0fixed power removed at condenser
BooleanfixReboilerPower (from DistillationBalance)false
Modelica.Units.SI.Powerwb (from DistillationBalance)0fixed power added at reboiler
BooleanisPartialCondenser (from DistillationBalance)falseif true the distillate leaves as gas
FreeFluids.Types.ReboilerTypereboilerType (from DistillationBase)FreeFluids.Types.ReboilerType.Kettle
Equilibrium data
Types.EquilibriumDataTypedataType (from DistillationBase)Types.EquilibriumDataType.EquilibriumConstant
Realfkl (from DistillationBase)0Feed equilibrium constant of light component
Realfkh (from DistillationBase)0Feed equilibrium constant of heavy component
Realdkl (from DistillationBase)0Distillate equilibrium constant of light component
Realdkh (from DistillationBase)0Distillate equilibrium constant of heavy component
Realbkl (from DistillationBase)0Bottom equilibrium constant of light component
Realbkh (from DistillationBase)0Bottom equilibrium constant of heavy component
Realfrv (from DistillationBase)0Feed relative volatility light/heavy
Realdrv (from DistillationBase)0Distillate relative volatility light/heavy
Realbrv (from DistillationBase)0Bottom relative volatility light/heavy
Real[nNonKey]fkifill(0, nNonKey)Feed equilibrium constant of nonkey components
Real[nNonKey]dkifill(0, nNonKey)Distillate equilibrium constant of nonkey components
Real[nNonKey]bkifill(0, nNonKey)Bottom equilibrium constant of nonkey components
Real[nNonKey]frvi_hfill(0, nNonKey)feed relative volatility of nonkey components
Real[nNonKey]drvi_hfill(0, nNonKey)distillate relative volatility of nonkey components
Real[nNonKey]brvi_hfill(0, nNonKey)bottom relative volatility of nonkey components

Components

TypeNameDefaultDescription
Modelica.Units.SI.MolarFlowRateF (from DistillationBalance)feed molar flow rate
Modelica.Units.SI.MolarFlowRateD (from DistillationBalance)distillate molar flow rate
Modelica.Units.SI.MolarFlowRateB (from DistillationBalance)bottom purge molar flow rate
Modelica.Units.SI.MolarMassFmw (from DistillationBalance)feed molar mass
Modelica.Units.SI.MassFlowRateFmass (from DistillationBalance)feed mass flow rate
Modelica.Units.SI.MassFlowRateFflMass (from DistillationBalance)feed light mass flow rate
Modelica.Units.SI.MassFlowRateFfhMass (from DistillationBalance)feed heavy mass flow rate
Modelica.Units.SI.MolarMassDmw (from DistillationBalance)distillate molar mass
Modelica.Units.SI.MassFlowRateDmass (from DistillationBalance)distillate mass flow rate
Modelica.Units.SI.MassFlowRateDflMass (from DistillationBalance)distillate light mass flow rate
Modelica.Units.SI.MassFlowRateDfhMass (from DistillationBalance)distillate heavy mass flow rate
Modelica.Units.SI.MolarMassBmw (from DistillationBalance)bottom molar mass
Modelica.Units.SI.MassFlowRateBmass (from DistillationBalance)bottom mass flow rate
Modelica.Units.SI.MassFlowRateBfhMass (from DistillationBalance)bottom heavy mass flow rate
Modelica.Units.SI.MassFlowRateBflMass (from DistillationBalance)bottom light mass flow rate
FractionFZl (from DistillationBalance)feed molar fraction of light component
FractionFZlMass (from DistillationBalance)feed mass fraction of light component
FractionFZh (from DistillationBalance)feed molar fraction of heavy component
FractionFZhMass (from DistillationBalance)feed mass fraction of heavy component
Modelica.Units.SI.MolarFlowRateFfl (from DistillationBalance)feed molar flow rate of light component
Modelica.Units.SI.MolarFlowRateFfh (from DistillationBalance)feed molar flow rate of heavy component
FractionDYl (from DistillationBalance)distillate molar fraction of light comp.
FractionDYlMass (from DistillationBalance)distillate mass fraction of light comp.
FractionDYh (from DistillationBalance)distillate molar fraction of heavy comp.
FractionDYhMass (from DistillationBalance)distillate mass fraction of heavy comp.
Modelica.Units.SI.MolarFlowRateDfl (from DistillationBalance)distillate molar flow rate of light component
Modelica.Units.SI.MolarFlowRateDfh (from DistillationBalance)distillate molar flow rate of heavy component
FractionBXl (from DistillationBalance)bottom molar fraction of light comp.
FractionBXlMass (from DistillationBalance)bottom mass fraction of light comp.
FractionBXh (from DistillationBalance)bottom molar fraction of heavy comp.
FractionBXhMass (from DistillationBalance)bottom mass fraction of heavy comp.
Modelica.Units.SI.MolarFlowRateBfl (from DistillationBalance)bottom molar flow rate of light component
Modelica.Units.SI.MolarFlowRateBfh (from DistillationBalance)bottom molar flow rate of heavy component
FractionFgf (from DistillationBalance)feed gas fraction
RealQ (from DistillationBalance)feed enthalpy parameter
RealRR (from DistillationBalance)internal reflux ratio, after heating of subcooled reflux
RealRRext (from DistillationBalance)external reflux ratio, as observed at condenser
RealBR (from DistillationBalance)boil ratio
Modelica.Units.SI.PowerWc (from DistillationBalance)heat removed at condenser
Modelica.Units.SI.PowerWb (from DistillationBalance)heat added at the reboiler
Modelica.Units.SI.MolarFlowRateLe (from DistillationBalance)liquid molar flow rate in the enrichment section
Modelica.Units.SI.MolarFlowRateVe (from DistillationBalance)gas molar flow rate in the enrichment section
Modelica.Units.SI.MolarFlowRateLs (from DistillationBalance)liquid molar flow rate in the stripping section
Modelica.Units.SI.MolarFlowRateVs (from DistillationBalance)gas molar flow rate in the stripping section
RealRmin (from DistillationBase)minimum reflux ratio
FractionRYl (from DistillationBase)operation lines intersection Y value
FractionRXl (from DistillationBase)operation lines intersection X value
RealFrv (from DistillationBase)relative volatility of light/heavy at feed conditions
RealDrv (from DistillationBase)top of column relative volatility l/h
RealBrv (from DistillationBase)bottom of column relative volatility l/h
RealRVavg (from DistillationBase)average relative volatility l/h
RealNstagesMin (from DistillationBase)number of states at total reflux
RealNstages (from DistillationBase)number of stages at working reflux
RealNplates (from DistillationBase)number of plates needed at working reflux
RealFstage (from DistillationBase)optimum stage for fee, from top
Fraction[nNonKey]FZimolar fraction of nonkey components in feed
Modelica.Units.SI.MolarFlowRate[nNonKey]Ffimolar fraction of nonkey components in feed
Fraction[nNonKey]DYimolar fraction of nonkey components in distillate
Modelica.Units.SI.MolarFlowRate[nNonKey]Dfimolar fraction of nonkey components in distillate
Fraction[nNonKey]BXimolar fraction of nonkey components in bottom
Modelica.Units.SI.MolarFlowRate[nNonKey]Bfimolar fraction of nonkey components in distillate
Real[nNonKey]Frvi_htop of column relative volatility nonkey/heavy
Real[nNonKey]Drvi_htop of column relative volatility nonkey/heavy
Real[nNonKey]Brvi_hbottom of column relative volatility nonkey/heavy
Real[nNonKey]RVavgi_haverage relative volatility nonkey/heavy
RealSigma1
Real[nNonKey]Sigma1i
RealSigma2
Real[nNonKey]Sigma2i
RealThetaUnderwood parameter
RealRmin1minimum reflux for type 1 pinch point