Hello all....!!!
Have taken a gap..... because of my profession!!
After many days gone through the mail and comments, and most of the requests were related to vent sizing, ejectors. And even in my circle many were expressing the same query related to vent line sizing.
Also recently have seen a disaster in a chineese chemical plant (Xiangshui), nearly 78 persons were dead and many went missing, which is majorly due to lack of proper safety measures and may be avoidable if adequate planning for relief is available. And especially it can be considered as a black day in every process engineer's career. See Below video:
As per every management, Process safety is everyone's responsibility, but as per me a Process Engineer holds major accountability than others. Hence being an process engineer, its quite necessary to evaluate every hidden risk before going further with a process. As in a manufacturing company process engineer would be like a head light defining right direction to the operators.
As per every management, Process safety is everyone's responsibility, but as per me a Process Engineer holds major accountability than others. Hence being an process engineer, its quite necessary to evaluate every hidden risk before going further with a process. As in a manufacturing company process engineer would be like a head light defining right direction to the operators.
Today here i gonna add a small notes related to that,
Actually its a very simple task if we have all the data which is required for calculation, as based on assumption its quite waste to calculate because the calculation would be something like satisfying ourselves.
Maximum allowable working pressure: Maximum gauge pressure at a given/designated temperature.
Operating pressure: Gauge pressure, usually 10-12% below the MAWP.
Over pressure: Pressure above the set pressure during relief process.
Back pressure: Outlet pressure of the relief device during pressure knocking. Usually it will be equal to atmospheric pressure.
Accumulated pressure: Pressure raise above MAWP during relief process.
You might also like:
[How to] Select a FBD blower ?
[How to] Scale-up process equipment's ?Now lets go into the calculation part.
Usually the vent size shall be calculated from the data obtained from ARC(Accelerated rate calorimetry), RSST(Reactive System Screening Tool), ARSST(Advanced Reactive System Screening Tool).
2. Rate of pressure raise (P, psi/min).
p [1 + (1.98x10^3)/(p^1.75)]^0.285 Cd
where,
V - Volume of gas / reactant (m3),
Cd - discharge coefficient = Actual flowrate / theoretical flowrate,
Now, i'll demonstrate with a case study,
Details of the reaction from the ARC study in run away case:
Rate of pressure raise = 0.05 bar/min = 0.73 psi/min,
Set pressure of the relief device = 4.0 Kg/cm2 = 56.89 psig,
Volume of the gas = 30 m3, (This may seems higher, but it is quite regular as the gas densities were very low),
Then using above formula,
A = ( 3.5 x 10^-3 ) / [49.78 x (1 + (1.98*10^3 / 49.78^1.75)^0.285] x [30*(0.2+0.73)/0.7]
= 0.00203 m2.
Dia, d = ((4 / 3.141) x A)^0.5 = 0.0507 m = 5.07 cms = 1.999 in ~2 in.
This is quite lengthy and requires most of the data from the vendor side,
the details include the constants like discharge coefficient, back pressure correction factor.
And other details include the vapour/gas compressibility factor, heat capacity ratio, Mass flowrate.
Here is the formula for the vent size of a storage vessel,
A = [ w / C x Kd x P1 x Kb ] x ( T x Z / M )^0.5.
w - mass flowrate of the solvent vapour,
C = [ j x g / R ] x [ (2 / ( j + 1 )) ^ (( j+1)/(j-1))],
j - Heat capacity ratio of vapour or gas,
R - universal gas constant,
P1 (psia) - Pmax (psig) + 14.7,
Kd - discharge coefficient (~ 0.7),
Kb - back pressure correction factor (~1.0),
Now, let me demonstrate with a case study by calculating the required vent size for a storage vessel of 60 KL capacity containing toluene,
Heat capacity ration of the vapour be 1.2, Solvent be methanol - mol wt.: 92 g/mole,
Vessel dia be 2.5 m, height of the vessel be 3 m (torispherical vessel),
Wetted surface area of the vessel (Awet) = ( 𝜋 d l ) + ( 1.084 x D^2 )
= 30.33 m2 = 326.5 ft2.
Total heat absorption, Q = 34500 x (Awet)^0.82 = 3973775.31 BTU / hr.
Let the rate of vaporization be 20 Kgs/sec [the value is quite lower than the actual, usually for hydrocarbons it will vary between 30-50 Kg/sec as they possess lower latent heats] ,
P1 = Pmax + 14.7 = (Pmawp x 1.21) + 14.7 = (4 x 1.21) + 1.013
= 5.853 bar abs = 5.853 x 10^5 N/m2.
Also Read:
Reduce drying time-cycles in ANFD effectively
Selection of filtration equipment
C = [ (1.3 x 1 / 8314) x (2/(1.3+1))^(1.3/0.3)]^0.5 = 7.32*10^-3,
A = ( 20 / (7.32*10^-3 x 0.7 x 5.853 x 10^5 x 1)) x ((273.15+110.65) x 1 / 92)^0.5
= 0.0067 x 2.042 = 0.0136 m2.
Dia, d = ( 4 x 0.0136 / 3.141 )^0.5 = 0.132 m = 13.2 cms = 5.185 in ~5.5 / 6 in.
The above demonstrated calculation would work great if relevant authenticated safety reports are available. If not, then as an engineer we should be in a position to provide a temporary way forward without giving a major break in production.
Vent Size Calculation without Safety Data:
In such cases we should depend on the theoretical material balance,
I'll give you a case study with a random chemical reaction,
Yeah, i got a Balz Schiemann reaction, where a Aniline undergoes diazotization with Nitrous acid followed by reaction with fluoroboric acid giving aryl diazonium salts and thermal decomposition of dry aryl diazonium salts gives aryl fluoride and bi-product of the reaction is Nitrogen gas and boron trifluoride.
Below is the scheme of the reaction:
C6H5NH2 + HNO2 + HBF4 ----- > C6H5F + N2 + BF3
Note: Property of Pharma Engineering |
Now lets suppose, we are manufacturing fluorobenzene with a batch size of 100 Kgs.
Below are the raw materials quantities used in the batch manufacturing:
Aniline - 100 Kgs,
HNO2 - 50 Kgs,
Fluoroboric acid - 150 Kgs.
Below is the material balance for the above scheme with the raw material Qty.'s:
So based on the material balance, it can be said that
29.80 Kgs of Nitrogen and 72.13 Kgs of Boron trifluoride is liberating.
As mostly the Diazonium salts are highly unstable at lower temperatures, the preferable temperature of addition would be below 0° C. Hence the first reaction would be formation of Diazo salts, and nitrogen would be liberating in this point and it depends on the addition rate of HNO2, as the Nitrogen formation is attributed by Nitrous acid addition.
Also Read:
[How to] Calculate the product content in liquid without concentrating ?
Bond energy calculations - Alternate to RC1e study
Reaction - 1: Formation of Diazonium salts
Hence in the first step, Nitrogen would be liberating and the Qty.'s is 29.80 Kgs.
Converting into volumes, it would be 29.80 / 1.25 = 23.84 m3, which is quite higher.
and let the addition of nitrous acid took 2 hours,
In 1st case, lets consider Nitrogen is liberating uniformly throughout the addition,
i.e., 23.84 / (2 x 60) = 0.199 m3/min = 199 L /min.
But while performing the calculations through material balance, we need to consider the worst scenario's, as liberation would go in 10 minutes with a peak liberation of 90%.
The worst scenario would be 23.84 x 0.9 / 10 = 2.15 m3/min.
Now to position ourselves in a safe state during addition, we shall make the system ready to knock out 2.15 m3 of nitrogen per minute,
So, now one more assumption we need to consider is velocity of nitrogen, i.e., based on my experience i'll consider as 10 m/sec.
So, the relief device(vent) area shall be (A) = (2.15/60) / 10 = 0.0036 m2.
Dia, d = ( 0.0036 x 4 / 3.141 )^0.5 = 0.0675 m = 6.75 cms = 2.66 in ~3 in.
Hence for the first reaction we should provide a vent size of 3 inch.
Reaction - 2: Formation of aryl fluoride
Now enter the dragon, lets move to next reaction,
i.e., formation of fluorobenzene, where boron trifluoride is a bi-product and the gas is quite toxic and the same can't be left into atmosphere, it need to be scrubbed.
If you need a scrubber design, pl refer below link, in this subject post i'll describe only about vent size.
72.13 Kgs of Boron trifluoride is being liberated,
Lets follow the same procedure like above,
Converting it to volume, it will be 72.13 / 2.76 = 26.13 m3.
Same 90% of peak load in 10 minutes,
it will be 26.13 x 0.9 / 10 = 2.35 m3 /min.
And as this time, i'll be supported by scrubber, hence i'll consider the velocity as 10 m/sec,
Here you might get a doubt that, "why am i considering velocity as 10 m/sec even having a scrubber advantage?".
Answering your query, "dear engineers, if you can see the density above, it is 2.76 Kg/m3 which is more than double of nitrogen(1.25 Kg/m3)". Hope you got it.
Let's get back into calculation, 2.35 m3 of BF3 to be knocked out in one minute.
Relief device(vent) area shall be (A) = (2.35/60) / 10 = 0.0039 m2.
Dia, d = ( 4 x 0.0039 / 3.141 ) ^ 0.5 = 0.071 m = 7.065 cms = 2.78 in ~ 3 inches.
Hence, from the reaction 2 it is observed that a minimum size of 3 inches is required for both reactions.
And here the lucky fact is both are liberating in different reactions, or else it would be a complicate case to solve.
Reaction - 3: Knocking out the Unreacted compounds
And also we need to check the vent sizing during quenching operations, as in the above case, there is an excess reagent used in process i.e., HBF4,
If water is used for quenching the unreacted HBF4, then the reaction would be as below:
So here, in this case, lets suppose 100 lts of water used for quenching the unreacted Fluoroboric acid. products will be boric acid (Solid) and hydrogen fluoride (gas).
And the sane calculation shall be repeated. [I'll leave it and let it be an exercise to the readers].
That's it......!!!
Hope everything is understandable to all.......!!!
Any queries feel free to comment / message....!!!
This post is dedicated to my Dearest brother Govindarajan V.
Finally, I'm expressing my view here, "If you are unable to Knockout a process risk, you'll be knocked out by it".
Related Posts:
[How to] Design a decanter for workups in Pharma ?
Agitated Thin Film Dryer Design calculations
[How to] Calculate vapour column diameter ?
[How to] Perform bond energy calculation ?
About The Author
Thank you so much sir, you are really doing a great work.sir if possible then can you write something on ejector system both water and steam and calculation involved in it's design and selection
ReplyDeleteA = ( 3.5 x 10^-3 ) / [3.5 x (1 + (1.98*10^3 / 49.78^1.75)^0.285] x [30*(0.2+0.73)/0.7]
ReplyDelete= 0.00185 m2.
in above calculation, u r taken pressure in PSI & kg/cm2 in some case, i am not understand how it will come.
Dear Nitin,
DeleteThanks for correcting, modified the calculation, the vent size increased from 1.911 inch to 1.999 inch.
Best Regards,
AJAY K
How to interpret heat of reaction
ReplyDeleteHii prasad,
DeleteOnce go through these posts: https://www.pharmacalculations.com/2018/07/reaction-calorimetry-RC1e.html
https://www.pharmacalculations.com/2017/07/bond-energy-calculations.html
https://www.pharmacalculations.com/2019/06/how-to-perform-energy-balance.html
If any queries feel free to connect with me at pharmacalc823@gmail.com
Regards,
AJAY K
Thankyou sir .
ReplyDeleteThanks for the information. Just wanted to ask where can I find the derivation or explanation of the formula you used above.
ReplyDeleteHii ,
DeleteIt is taken from API 520 guide
Thanks for the info. Can you give me the link for the explanation or derivation of the formula you used above?
ReplyDeleteHii ,
DeleteIt is taken from API 520 guide
Ajay this post is very nice, now i am working on a project size selection of PSV& rupture disc.this post is very helpful for me.
ReplyDeleteI have a one quary regarding W value.how you put w value 20 kg/sec.pls elabrate.
In A = ( 20 / (7.32*10^-3 x 0.7 x 5.853 x 10^5 x 1)) x ((273.15+110.65) x 1 / 92)^0.5
When Q=3973775.31 BTU / hr.
Dear ,
DeleteQ is the vaporation losses, in the above line i've mentioned it, thats a considered value only.
Next time pl comment with your good name.
Dear ajay,
ReplyDeleteSorry for that i didn't mention my name next time it will not done.
I m reffering Crosby engineering handbook
Where in Pg. No 7.23
W = Q / Hvap
As per this book we have to calculate Q also. We can't consider .
Please brief me.
Regards
Raj
Are you referring to the vaporation losses ?
DeleteIf you have any psv /rupture disc sizing
ReplyDeleteCalculation sheet please share with me.
Regards
Raj
Hii Raj,
DeletePl share your id at pharmacalc823@gmail.com.