SILProof / White papers / Markov Modelling

White paper · Quantitative methods

The 5% that decides your 2oo3

Markov modelling for functional safety, from failure states to a PFDavg you can defend. One worked 2oo3 model, solved two independent ways, and the single assumption that moves the answer by a factor of 48.

13 October 2026 · SILProof Engineering · ~12 minute read

Download the full paper (PDF, 8 pages) Five connected state cards with a dashed common-cause path, illustrating a Markov model of 2oo3 voting

Three level transmitters. Two-out-of-three voting. Any single transmitter can fail and the trip still works. It feels like the safest sensor arrangement on the P&ID, and against independent failures it nearly is. Then you add one modest assumption, that 5% of dangerous failures hit all three channels together, and the calculated PFDavg gets 48.5 times worse. This paper builds the model that shows why, step by step, and checks it against an independent closed-form solution.

1.339 × 10−52oo3 PFDavg with β = 5%, 60 FIT per channel, annual proof test
2.761 × 10−7The same group if common cause is ignored (β = 0)
48.5×The gap. About 98% of the 2oo3 PFDavg comes from the common-cause path

The paper in 60 seconds

  1. States are engineering decisions. Keep a state only if it changes demand success, diagnostics, repair or test response.
  2. A proof test is scheduled, not random. A plain Markov chain is memoryless, so the test is applied as a separate reset.
  3. Redundancy kills independent failures but leaves a common-cause floor. That floor usually decides the number.
  4. Exact maths can still be wrong for the plant if data boundaries, coverage, repair policy or dependencies are wrong.
  5. PFDavg is one gate of three. Architecture and systematic capability remain separate requirements.

1. Why reach for a Markov model at all?

For a single unrepaired channel, the familiar shortcut is PFDavg ≈ λDU × T / 2. It is good when its assumptions hold and the failure probability is small.

Once voting, diagnostics, repair and common causes enter, the shortcut hides more than it shows. A Markov chain does the opposite. It lists every state the system can be in and tracks how probability moves between them over time, so every assumption sits in plain view as a state or an arrow.

That is not a claim that Markov always beats fault trees, reliability block diagrams or formulas. The rule is simpler: use the simplest method that captures the dependencies and time behaviour your decision depends on.

2. Five states are enough

Track the 2oo3 group by how many channels have failed dangerously and undetected since the last proof test: 0, 1, 2 or 3. Add one more state, C, for a common-cause event that takes out the whole group at once. With two or more channels gone, or in C, the group cannot act on a demand.

Proof test every T = 8,760 h: every state resets to 0 (perfect test) 0 failedAvailable1 failedAvailable2 failedDangerous3 failedDangerous 3λᵢ 2λᵢ λᵢ C · common causeDangerous · absorbing until the test λc from any state Green: at least two channels can still vote. Red: the 2oo3 group cannot act on demand.
Figure 1. The five-state continuous-time Markov chain. Solid arrows are independent channel failures; dashed gold arrows are the common-cause shock, possible from any state.

The worked numbers: λDU = 60 FIT per channel (1 FIT = 10−9 failures per hour) and β = 5%. The β-factor splits each channel's rate into an independent part λi = (1 − β) λDU = 57 FIT and a group-wide common-cause rate λc = β λDU = 3 FIT.

FromToRateWhy
013λiAny of three healthy channels can fail
122λiTwo healthy channels remain
23λiThe last channel fails
0, 1, 2, 3CλcOne shock disables the group

C is absorbing until the next test. Treating common cause as one group-wide shock is deliberately simple; real systems may need partial common-cause groups.

3. The whole method in two lines

Put the state probabilities in a row vector p(t) and the transition rates in a generator matrix Q. Each off-diagonal entry is a rate and each row sums to zero. Then:

dp/dt = p(t) Q   ⇒   p(t) = p(0) exp(Qt)

Mark the dangerous states with a reward vector r = [0, 0, 1, 1, 1]. The instantaneous unavailability is p(t)·r, and the number we report is its average over the test interval:

PFDavg = (1/T) ∫0T p(t) r dt
The subtle part: exponential waiting times are memoryless, and a proof test every 12 months is not. So the chain governs the time between tests, and the test itself is a scheduled reset applied at t = T. This is a finite-horizon average, not a steady-state availability.

4. Solve it twice, from two directions

A result is only as credible as its cross-check. For this model there is an independent route that never mentions a matrix. A surviving channel has failed by time t with probability q(t) = 1 − e−λit, and a common-cause shock has occurred with probability c(t) = 1 − e−λct. The group succeeds only if no shock occurred and at least two channels survive:

PFD(t) = 1 − 3e−(λc+2λi)t + 2e−(λc+3λi)t

Averaged over T, with A(a,T) = (1 − e−aT)/(aT):

PFDavg = 1 − 3A(λc+2λi, T) + 2A(λc+3λi, T)

The state-space answer and the event-based answer must agree for identical assumptions. If they do not, one of them is wrong. (In software, evaluate small aT with expm1 to avoid cancellation.)

5. The result

00.51.01.52.02.53.0 × 10⁻⁵ 03 mo6 mo9 mo12 mo · test PFDavg = 1.339×10⁻⁵ (area ÷ T) 2oo3 with β = 5% 2oo3 with common cause removed (β = 0) Time since the last proof test
Figure 2. Instantaneous PFD across one perfect proof-test interval. The shaded area divided by T is PFDavg. With common cause removed, the curve barely leaves the axis.
Method / scenarioPFDavg
Markov chain (CTMC), β = 5%1.338904489e−5
Independent closed form, β = 5%1.338904489e−5
No common-cause mechanism (β = 0)2.760739324e−7
Common-cause shock alone1.313988489e−5

The two methods agree to numerical precision. Look at the last row: the common-cause shock on its own is 1.314 × 10−5, almost the entire answer. The redundancy is doing its job against independent failures. It just cannot touch the failure that hits all three channels together.

Redundancy suppresses independent failures. It does nothing for the failure all three channels share.The central lesson of the 2oo3 example

6. The 5% that decides everything

Hold each channel at 60 FIT and change only β. The result is nearly a straight line: PFDavg is governed by the common-cause assumption far more than by the voting.

0123456 × 10⁻⁵ 0%5%10%15%20% 1.34×10⁻⁵2.65×10⁻⁵5.27×10⁻⁵ Common-cause β factor (λDU per channel held at 60 FIT)
Figure 3. 2oo3 PFDavg against the assumed β. Illustration only; the curve is not an empirical β estimate.
β0%1%3%5%10%20%
PFDavg2.761e−72.899e−68.144e−61.339e−52.650e−55.273e−5

Now the surprise. With the same data, here is how three architectures compare over one year:

ArchitecturePFDavgWhat has to happen
1oo12.628e−4The one channel fails
1oo21.322e−5Both channels fail, or a group-wide common cause
2oo31.339e−5Two of three fail, or a group-wide common cause

1oo2 and 2oo3 land within 1.3% of each other, because the shared common-cause floor dominates both. That says nothing about spurious trips, diagnostics or operating preference, which need their own models. It does say that adding a third transmitter does not buy dangerous-failure performance while the common cause stays where it is.

Try it: move β and watch the floor

This uses the paper's closed form for a DU-only group with a perfect proof test. It is for building intuition, not for a SIL claim.

2oo3 PFDavg
1oo2 PFDavg
1oo1 PFDavg
vs no common cause

7. From a maths model to a plant model

The worked example is deliberately clean. A real SIF needs more, and each item below changes the model, not just the numbers.

Detection and repair

A dangerous detected failure might trip the plant, call maintenance, or leave the function running impaired under compensating measures. Those are different states. Repair can be a rate μ = 1/MTTR where that approximation is defensible; separate detection time, administrative delay, hands-on repair and restoration testing where they matter.

Imperfect and partial proof tests

Model the test as a state-transfer matrix R at t = T. A perfect test sends everything back to healthy. An imperfect one must keep the probability in failure modes it cannot reveal. Resetting everything would understate PFDavg. Partial tests may need multi-phase models.

Dependence is wider than one β

Shared impulse lines, calibration methods, power supplies, environment, design defects and maintenance practice are all coupling mechanisms. A constant β is a transparent approximation, not a physical explanation, so record the boundary and source behind it.

What a plain Markov chain cannot do

Fixed test dates, fixed repair durations, wear-out and demand damage break the memoryless assumption. Piecewise or phased models, semi-Markov methods, age states or Monte Carlo simulation are the remedies. More states do not fix the wrong stochastic assumptions.

8. A review checklist you can use on Monday

Credibility comes from checks anyone can repeat, not from the complexity of a diagram. Before a Markov result goes into a report, confirm:

  • Probability is conserved: p(t) sums to one and every state stays within [0, 1].
  • The generator is valid: Qij ≥ 0 off the diagonal and each row sums to zero.
  • Boundaries behave: PFD(0) = 0 after a perfect test, and β = 0 reproduces independent voting.
  • A second method agrees: closed form or numerical integration matches the chain.
  • Sensitivities point the right way: a longer test interval or higher β never improves the result.
  • The reset is honest: perfect tests reset covered states; imperfect tests keep uncovered ones.
  • Data is traceable: channel boundary, failure-rate category, source version, units and β meaning are recorded.
  • Numerics are controlled: solver tolerance and integration convergence are stated.
  • The claim is separated: a PFDavg pass is not reported as architectural or systematic compliance.

For this benchmark, the chain and the closed form differ by less than 10−12 in absolute PFDavg. That checks the mathematics of the example; it is not evidence about any particular software implementation.

9. Leave a record another engineer can rerun

"We used a Markov model" is not a calculation record. The minimum is the state dictionary, the transition table, Q, initial conditions, the demand-failure reward, the reset policy, the analysis horizon, the integration method, parameter sources and the numerical checks. With those, anyone can rebuild the result: assemble the 5 × 5 generator, evaluate p(0)exp(Qt), sum states 2, 3 and C, average over 8,760 h, and compare with the closed form.

10. The takeaway

A Markov result is not credible because it uses matrices. It is credible when the physical model, parameter meanings, transitions, test policy and validation all survive review. The 2oo3 example teaches three things: redundancy suppresses independent simultaneous failures; a modest common-cause assumption can dominate what is left; and an independent cross-check can test a model before anyone relies on it.

Model the failure mechanism before calculating the probability. Then show why the model suits the safety function.The central principle

Inside the IEC 61508 / IEC 61511 lifecycle, Markov modelling adds rigour to quantitative SIL verification. It does not replace risk allocation, architectural constraints, systematic integrity, field evidence or independent review. For how those fit together, read The Three Gates to SIL.

Abbreviations and symbols
CCFCommon-cause failure
CTMCContinuous-time Markov chain
DD / DUDangerous detected / dangerous undetected failure
FITFailures in time: 10−9 failures per hour
PFDavgAverage probability of dangerous failure on demand
βBeta factor (simplified common-cause allocation)
λi / λcIndependent channel / group common-cause failure rate
Q / rGenerator matrix / demand-failure reward vector
TProof-test interval
Technical references
  1. IEC 61511-1:2016+AMD1:2017, Functional safety: safety instrumented systems for the process industry sector, Part 1.
  2. IEC 61508-6:2010, Guidelines on the application of IEC 61508-2 and -3.
  3. Innal, F., Lundteigen, M.A., Liu, Y., Barros, A. (2016). PFDavg generalized formulas for SIS subject to partial and full periodic tests based on multi-phase Markov models. Reliability Engineering & System Safety 150, 160–170. doi:10.1016/j.ress.2016.01.022
  4. Ottermo, M.V., Hauge, S., Håbrekke, S. (2021). Reliability Data for Safety Equipment: PDS Data Handbook, 2021 Edition. SINTEF Digital.
  5. Trivedi, K.S., Ciardo, G., Malhotra, M., Sahner, R.A. (1993). Dependability and performability analysis. NASA-CR-191565.
  6. IEC 61508-2:2010, Requirements for E/E/PE safety-related systems.
  7. Butler, R.W. (1985). An Abstract Specification Language for Markov Reliability Models. NASA TM-86423.
About this white paper
An explanatory methods paper built around one controlled mathematical benchmark. It is not an assessment or SIL certificate, does not validate the SILProof calculation engine, and is not a substitute for the licensed IEC standards, project requirements or competent review. Download the PDF.