Monday, 3 April 2023

Process Capability & Performance Indices (Cp, Cpk, Pp, Ppk)

Hello Everyone ....!!!

Hope everything is fine with all. Today i'm trying to post about something new, which will have certain value addition to a chemical engineers profile. The topic is related to statistical analysis and specifically about understanding the Process Capability & Performance Indices. 

What does these mean ?
How to interpret these things ?
What actions can we make based on these indices ?

This is a query that i've received on my mail from one of a process engineer. Today here lets see how to calculate and understand these things in a process engineer's perspective and correlate them with our regular work.

Before going to the calculation part, lets understand the basic terminology.

What is capability ?

Capability is an ability to a certain task and the capability varies from task to task.

What is capability analysis ?

Analysing the capability of a specific process is nothing but capability analysis.

How can we classify capability ?

Capability shall be classified into two types i.e., Short term capability and Long term capability.
Short term capability shall be interpreted using the indices Cp & Cpk,
Long term capability shall be interpreted using the indices Pp & Ppk.

What does Cp / Pp indicate ?

Cp / Pp indicates Process potential index i.e., whether the data points can fit inside the specification limit that is laid by the customer or not.

What does Cpk / Ppk indicate ?

Cpk / Ppk indicates Centering index i.e., whether the data points were centrally located or not.

What is the differnece between Cp & Pp ?

Cp indicates the process potential based on the sample data, whereas the Pp indicates the same based on population data.

Similarly Cpk & Ppk indicates the centering index based on the sample data and population data.

If you still have confusion, all i can say is just remember the indices Cp & Cpk for analysing process capability when sample data is available and Pp & Ppk for analysisng the process capability when population data is available.

What is sample and population ?

Population is a collection of all data points related to a specific process, whereas sample is a set of data points which are picked from a population depending on specific methodology.

What are the various sampling methodologies ?

There are four basic sampling methodologies i.e.,
Random Sampling
Sequential Sampling
Stratified Sampling
Systematic Sampling

How to estimate / calculate the capability indices ?

We can use either manual calculations / Excel based calculations / Calculations using softwares like Minitab.

What are the formulae for manual calculation of these indices ?

Cp = (USL - LSL)/(6 x S)

Pp = (USL - LSL)/(6 x 𝜎)

Cpk = Min (Cpu, Cpl)

Cpu = (USL - x)/(3 x S); Cpl = (x - LSL)/(3 x S)

Ppk = Min (Ppu, Ppl)

Ppu = (USL - X)/(3 x 𝜎); Ppl = (X - LSL)/(3 x 𝜎)

where,
USL - Upper Specification Limit; LSL - Lower Specification Limit
S - Sample standard deviation; 𝜎 - Population Standard deviation
x - Sample mean; X - Population mean

Pro's & Con's of Long term & Short term capability analysis?

Analysing process capability using long term capability is difficult as we need to analyse huge data, whereas the analysis going to be accurate.

Analysing process capability using short term capability is easy as the data is set of sample points, but the results are not accurate as long term capability. 

Lets try to calculate the above said indices with some random data, but before that we should see about the graphs used to represent data distribution and their interpretation using the indices.

Below is the normal distribution curve:

The above shown bell shaped curve indicates the distribution of data points related to a specific process and the process width indicates voice of process which is nothing but how the process delivers.

The graph shall be plotted with below steps:

i. Arrange the data points in ascending order

ii. Form certain classes for data i.e., if the specification limits are 70 - 120 Kg, then classes are going to be like 71 - 80, 81 - 90, 91 - 100, 101 - 110, 111 - 120.

iii. Now identify the number of data points (Frequency) that fall within these classes, for example there are 50 data points that we have collected and the total data points that falls within 71 - 80 is 5, 81 - 90 is 10, 91 - 100 is 15, 101 - 110 is 10 & 111 - 120 is 5.

iv. Plot a bar chart for classes on X - axis and frequency on Y - axis

v. Connect the center points of bars and it will form a bell shaped curve, just like below:
Here in the above plot, the bell shaped curve looks like having sharp edges which is due to the uniform data that i've considered. However, if the data is from real time, then the curve will take a different form.

Example:
Now, lets try to calculate the capability indices from a sample data which are the yields of a pharmaceutical product and the voice of customer is 70 - 120 Kg
Data: 80, 83, 97, 101, 104, 91, 72, 95, 111, 99, 89, 78, 107, 106, 90, 92, 99, 93, 80, 105, 95, 99, 100, 91, 85, 79, 81, 95, 95, 98, 92

Solution:
Voice of customer = Specification limits
Therefore, Upper specification limit = 120 & Lower specification limit = 80;

Here, we got only sample data, hence we need to estimate the Cp & Cpk,
If population data is available, then we need to estimate Pp & Ppk

Sample mean (x) = Average of data points = (Sum of all data points / No. of data points) = 92.97

Sample standard deviation (S) = Σ(xi - x)^2 / (n-1) = 9.62

Cp = (USL - LSL) / (6 x S) = (120 - 70) / (6 x 9.62) = 0.87

Cpk = Min(Cpu, Cpl)

Cpu = (USL - x) / (3 x S) = (120 - 92.97) / (3 x 9.62) = 0.94
Cpl = (x - LSL) / (3 x S) = (92.97 - 70) / (3 x 9.62) = 0.80

Therefore, Cpk = 0.80.

Now lets come to interpretation,
Process potential (Cp) itself is less than 1

that indicates process is not capable of delivering the output within the specification limits.

If Cp is between 1 and 1.33, process is capable, however there are few rejected batches,
If Cp is between 1.33 and 1.67, process is capable and reproducible,
If Cp is between 1.67 and 2, process is highly capable.

Now lets check the interpretation of Cpk,
In the above example, the Cp is less than 1, hence no need to discuss much about Cpk,
because if the data points are not within the specification limits, then automatically the points are not centered, hence no need to discuss about Cpk

If Cp is greater than, then we need to check the Cpk

Below is the representation of Cpk:

If Cpk = 1, the data points were spread across the specification limits and the centering is low
If Cpk = 1.33, the data points spread is reduced and the centering is increased,
If Cpk = 1.67, the centering is further increased
If Cpk = 2, the process is world class

Finally, i would like provide you a correlation between Cpk & Sigma level i.e.,
Sigma level = 3 x Cpk.

For Cpk = 2, Sigma level is 6, so its a Six Sigma process.

This is how, we need to calculate and interpret the process capability using capability indices.


Hope the post is clear to all...

For any queries, please write us at pharmacalc823@gmail.com

Comments are most appreciated.....!!!

SuperSurvey About The Author


Hi! I am Ajay Kumar Kalva, Currently serving as the CEO of this site, a tech geek by passion, and a chemical process engineer by profession, i'm interested in writing articles regarding technology, hacking and pharma technology.
Follow Me on Twitter AjaySpectator & Computer Innovations

Saturday, 1 April 2023

Stossels Criticality Index - Classification of reactions

G
ood day everyone ......!!

Hope everyone is safe and fine. Today i've made sometime to post on the classification of reactions using Stossel's Criticality Index (SCI).

Stossels Criticality Index (SCI) uses the heat of reaction data and thermal stability data to define a reaction class.

Before going to the post in depth, lets see certain nomenclature to have better understanding in the further steps.

What is Tprocess ?

Tprocess (Tp) refers to the operating temperature at which desired reaction occurs.

What is ΔT ad ?

ΔT ad refers to adiabatic temperature raise i.e., the temperature raise which happens when the energy cannot be transferred out of the system. The energy generates here is due to the reaction heat.

What is MTT ?

MTT refers to the Maximum Technical Temperature (MTT), which is the boiling point of solvent available in reaction mass

What is TMR ?

TMR refers to Time to Maximum Rate i.e., the time for the reaction to consume before reaching the onset of decomposition.

What is onset of decomposition ?

Onset of decomposition is the temperature where exothermic decomposition initiates.

What is decomposition ?

Decomposition is an undesired scenario during reaction maintenance which evolves due to sudden shoot up of energy of reaction mass either due to thermal accumulation or due to external supply of energy.

What is thermal accumulation ?

Thermal accumulation is the percentage of energy which is being available with reaction mass after complete dosing of reactant.

Thermal accumulation + Thermal conversion = 100%.

What is TMR, ad24hrs ?

TMR, ad24hrs refers to the onset of decomposition when the reaction mass is maintained for 24 hours under adiabatic conditions i.e., with failure of utility circulation. As the time passes, the onset of decomposition reduces.

What is MTSR ?

MTSR is the Maximum Temperature to Synthesis Reaction i.e., the ability of reaction mass to reach its maximum temperature which is attributed by the heat generated during reaction.

How is the MTSR calculated ?

MTSR shall be calculated as below,

MTSR due to accumulation = Tprocess + (% accumulation x ΔT ad)

MTSR in worst case = Tprocess + ΔT ad


Hope the terms are clear for everyone.

Now lets get into the topic, there shall be five classes in SCI.

The classes are being defined based on the parameters which were described above i.e., Tp, ΔT, TMR, ad24hrs, MTT.

Lets see the classes and their interpretation below.

Here is the representation of SCI:


Class - 1
Observation: In the subject class, MTSR is less than that of MTT and TMR, ad24hrs

Interpretation: Class - 1 is a safe scenario w.r.t. reaction heat, however if energy is supplied through utilities and the total solvent available in the reaction mass is vaporized due to inadequate utility supply to condensers, there shall be risk.

Mitigation plans required: Ensure proper utility to condensers


Class - 2
Observation: In the subject class, MTSR is less than that of TMR, ad24hrs and MTT is beyond TMR, ad24 hrs.

Interpretation: Class - 2 is safe w.r.t. heat released due to reaction. However, if there is any energy in from external sources i.e., the utilities like steam / hot water / hot oil systems, then there is a risk involved.

Mitigation plans required: Control the usage of utilities during reaction.


Class - 3
Observation: In the subject class, the MTSR is greater than the MTT and less than that of TMR, ad24hrs.

Interpretation: The reaction mass potential to reach high temperatures is being controlled by the solvent boiling, as giving more energy leads to reflux of solvent. Hence the subject class is safe till there is adequate solvent available for reflux.

Mitigation plans required: Ensure the process contains enough solvent volumes to compensate the losses due the energy release and ensure adequate utility circulation to condenser.


Class - 4
Observation: MTSR is higher than TMR, ad 24hrs and the MTT is less than TMR, ad24hrs.

Interpretation: The scenario is having risk, however as the MTT is less the risk can be mitigated by consuming more solvent in the process

Mitigation plans required: Ensure adequate utility circulation at condenser and avoid solvent losses.

Class - 5
Observation: MTSR and MTT are above TMR, ad24hrs

Interpretation: Decomposition cannot be detected and its a risky scenario at any point of time.

Mitigation plans required: Don't execute the process on plant scale. Take the process back to R&D and re-engineer it either by changing the solvent or by changing the raw materials to control the adiabatic temperature raise.

So, this is how we need to read the SCI and interpret the classes and provide certain mitigation plans to reduce the risk involved.

If you have any queries in the post, pl reach us at pharmacalc823@gmail.com 

Comments are most appreciated .....!!

for better understanding on the above post we have hosted a video on youtube, you can access the content for free


For more content, please keep on visiting us.
QuizMaker








A
bout The Author


Hi! I am Ajay Kumar Kalva, Currently serving as the CEO of this site, a tech geek by passion, and a chemical process engineer by profession, i'm interested in writing articles regarding technology, hacking and pharma technology.
Follow Me on Twitter AjaySpectator & Computer Innovations