Home| Journals | Statistics Online Expert | About Us | Contact Us
    About this Journal  | Table of Contents
Untitled Document

[Abstract] [PDF] [HTML] [Linked References]

 

Gain-function of two-unit non-identical warm standby nuclear power system with failure due to extremely high radiations and failure due to non-availability of heavy water D2O in nuclear power reactors

Ashok Kumar Saini

Associate Professor, Department of Mathematics, B. L. J. S. College, Tosham, Bhiwani, Haryana, INDIA.

Email: [email protected]

Research Article

 

Abstract               Introduction: A pressurized heavy-water reactor (PHWR) is a nuclear power reactor, commonly using enriched natural uranium as its fuel, that uses heavy water (deuterium oxide D2O) as its coolant and moderator. The heavy-water coolant is kept under pressure, allowing it to be heated to higher temperatures without boiling, much as in a PWR. While heavy water is significantly more expensive than ordinary light water, it yields greatly enhanced neutron economy, allowing the reactor to operate without fuel-enrichment facilities (mitigating the additional capital cost of the heavy water) and generally enhancing the ability of the reactor to efficiently make use of alternate fuel cycles. The 1979 Three Mile Island accident and 1986 Chernobyl disaster, along with high construction costs, ended the rapid growth of global nuclear power capacity. A further disastrous release of radioactive materials followed the 2011 Japanese tsunami which damaged the Fukushima I Nuclear Power Plant, resulting in hydrogen gas explosions and partial meltdowns  classified as a Level 7 event. The large-scale release of radioactivity resulted in people being evacuated from a 20 km exclusion zone set up around the power plant, similar to the 30 km radius Chernobyl Exclusion Zone still in effect. In Nuclear Reactor the leakage in form of radiations becomes highly dangerous to the lives of living beings. The radiations from the nuclear reactor are always under serious consideration due to fatal and miserable results to human race. Every precautions and extra care is taken to avoid any miss-happening due to radiations. But still due carelessness or due to failure of some equipment in Nuclear Reactors there occur leakage of radiations causing a major casualty. In the present paper we have taken two-dissimilar warm standby nuclear power system with failure due to extremely high radiations which we abbreviated as FEHR and failure due to Non-availability of heavy water D2O in nuclear power reactor abbreviated as FNAHW. When there are radiations of extremely high magnitude the working of unit stops automatically to avoid excessive damage of the units and when the unit comes in no normal position the repair of the units’ starts immediately. The failure time distribution is taken as exponential and repair time distribution as general. Using Markov regenerative point technique we have calculated different reliability characteristics such as MTSF, reliability of the system, availability analysis in steady state, busy period analysis of the system under repair, expected number of visits by the repairman in the long run and Gain-function. Special case by taking failure and repair as exponential have been derived and graphs are drawn.

Keywords: warm standby, extremely high radiations, non-availability of heavy water in nuclear power reactor, MTSF, Availability, busy period, Gain - function.

 

INTRODUCTION

Water makes an excellent moderator; the hydrogen atoms in the water molecules are very close in mass to a single neutron, and the collisions thus have a very efficient momentum transfer, similar conceptually to the collision of two billiard balls. However, in addition to being a good moderator, water is relatively effective at absorbing neutrons. Using water as a moderator will absorb enough neutrons that there will be too few left over to react with the small amount of 235 U in the fuel, again precluding criticality in natural uranium. Instead, in order to fuel a light-water reactor, first the amount of 235 U in the uranium must be increased, producing enriched uranium, which generally contains between 3% and 5% 235U by weight (the waste from this process is known as depleted uranium, consisting primarily of 238 U). In this enriched form there is enough 235 U to react with the water-moderated neutrons to maintain criticality. An alternative solution to the problem is to use a moderator that does not absorb neutrons as readily as water. In this case potentially all of the neutrons being released can be moderated and used in reactions with the 235 U, in which case there is enough 235 U in natural uranium to sustain criticality. One such moderator is heavy water, or deuterium-oxide. Although it reacts dynamically with the neutrons in a similar fashion to light water (albeit with less energy transfer on average, given that heavy hydrogen, or deuterium, is about twice the mass of hydrogen), it already has the extra neutron that light water would normally tend to absorb. Nuclear power is the fourth-largest source of electricity in India after thermal, hydroelectric and renewable sources of electricity. As of 2013, India has 21 nuclear reactors in operation in 7 nuclear power plants, having an installed capacity of 5308 MW and producing a total of 30,292.91 GW h of electricity while seven other reactors are under construction and are expected to generate an additional 6,100 MW. Despite the opposition, the capacity factor of Indian reactors was at 79% in the year 2011-12 compared to 71% in 2010-11. Nine out of twenty Indian reactors recorded an unprecedented 97% Capacity factor during 2011-12. With the imported uranium from France, the 220 MW Kakrapar 2 PHWR reactors recorded 99% capacity factor during 2011-12. The Availability factor for the year 2011-12 was at 89%. In Nuclear Reactor the leakage in form of radiations becomes highly dangerous to the lives of living beings. The radiations from the nuclear reactor are always under serious consideration due to fatal and miserable results to human race. In the present paper we have taken two-dissimilar warm standby system with failure due to extremely high radiations- FHER and failure due to non-availability of heavy water in nuclear power plant -FNAHW

Assumptions

  1. The failure time distribution is exponential whereas the repair time distribution is arbitrary of two non-identical units.
  2. The repair starts immediately upon failure of units and the repair discipline is FCFS.
  3. The repairs are perfect and start immediately as soon as the extremely high radiations of the system become normal. The radiations of both the units do not go extremely high.
  4. The failure of a unit is detected immediately and perfectly.
  5. The switches are perfect and instantaneous.
  6. All random variables are mutually independent.

 

SYMBOLS FOR STATES OF THE SYSTEM                                                                             

Superscripts: O, WS, SO, FEHR, FNAHW

Operative, Warm Standby, Stops the operation, Failure due to extremely high radiations, failure due to non-availability of heavy water in nuclear power plant respectively

Subscripts: nehr, ehr, nahw, ur, wr, uR

No extremely high radiations. Extremely high radiations, non-availability due to heavy water, under repair, waiting for repair, under repair continued respectively

Up states: 0, 1, 2, 9;

Down states: 3, 4, 5, 6, 7, 8, 10, 11

Regeneration point: 0, 1, 2, 4, 7, 10

States of the System                                                                                                                                                                     

0(Onehr, WSnehr ) One unit is operative and the other unit is warm standby and there is no extremely high radiations in both the units.

1(SOnehr, Onehr)

The operation of the first unit stops automatically due to extremely high radiations and warm standby units starts operating and there is no extremely high radiations.

2(FEHRehr,ur, Onehr)

The first unit fails and undergoes repair after failure due to extremely high radiations are over and the second unit continues to be operative with no extremely high radiations.

3(FEHRehr,uR, SOehr)

The repair of the first unit is continued from state 2 and in the other unit extremely high radiations occur and stops automatically due to extremely high radiations.

4(FEHRehr,ur, SOuehr)

The one unit fails and undergoes repair after the extremely high radiations are over and the other unit also stops automatically due to extremely high radiations.

5(FEHRehr,uR, FEHRehr, wr) The repair of the first unit is continued from state 4 and the other unit is failed due to extremely high radiations in it and is waiting for repair.

6 (Onehr, FEHRehr,ur)

The first unit is operative with no extremely high radiations and the second unit failed due to extremely high radiations is under repair.

7(SOnehr, FNAFnahw,ur)

The operation of the first unit stops automatically due to extremely high radiations and the second unit fails due non-availability of heavy water in nuclear power reactor and undergoes repair.

8(FEHRehr,wr, FNAHWnahw,uR)

The repair of failed switch is continued from state 7 and the first unit is failed after extremely high radiations and waiting for repair.

9(Onehr, SOuehr)

The first unit is operative and the warm standby dissimilar unit is under extremely high radiations

10(SOnehr, FNAHWnahw,ur)                                

The operation of the first unit stops automatically due to extremely high radiations and the second unit fails due to non-availability of heavy water in nuclear power reactor and undergoes repair after the extremely high radiations is over.

11(FEHRehr,wr, FNAHWnahw,uR)

The repair of the second unit is continued from state 10 and the first unit is failed due to extremely high radiations is waiting for repair.

Figure 1: The State Transition Diagram

 

TRANSITION PROBABILITIES

Simple probabilistic considerations yield the following expressions:

p01 = , P07 =  

p09 = , p12 = , p14 =  

P20= G1*( λ1), P22(3) = G1*( λ1)=p23, P72 = G2*( λ4),

P72(8) = G2*( λ4)= P78

We can easily verify that

p01 + p07 + p09 = 1, p12 + p14 = 1, p20 + p23 (=p22(3))= 1, p46(6)= 1 p60 = 1,

p72+ P72(5) + p74 = 1, p9,10 =1, p10,2 + p10,2(11) = 1                                                                                                                       (1)

And mean sojourn time is

µ0 = E(T) =                                                                                                                                                               (2)

Mean Time to System Failure

We can regard the failed state as absorbing

 

,

                                                                                                                                                                               (3-5)

Taking Laplace-Stiltjes transform of eq. (3-5) and solving for

 = N1(s) / D1(s)                                                                                                                                                                    (6)

Where

N1(s) = {  +

D1(s) = 1 -  

Making use of relations (1) and (2) it can be shown that θ0(0) =1, which implies that θ0(t) is a proper distribution.

MTSF = E[T] = d/ds θ0*(0)             = (D1(0) - N1(0)) / D1 (0)

 

     s=0

= ( +p01  + p01 p12  + p09 ) / (1 - p01 p12 p20 )

Where

 +  + ,  + , + (3),

AVAILABILITY ANALYSIS

Let Mi(t) be the probability of the system having started from state I is up at time t without making any other regenerative state. By probabilistic arguments, we have

The value of M0(t), M1(t), M2(t), M4(t) can be found easily.

The point wise availability Ai(t) have the following recursive relations

A0(t) = M0(t) + q01(t)[c]A1(t) + q07(t)[c]A7(t) + q09(t)[c]A9(t)

A1(t) = M1(t) + q12(t)[c]A2(t) + q14(t)[c]A4(t), A2(t) = M2(t) + q20(t)[c]A0(t) + q22(3)(t)[c]A2(t)

A4(t) = q46(3)(t)[c]A6(t), A6(t) = q60(t)[c]A0(t)

A7(t) = (q72(t)+ q72(8)(t)) [c]A2(t) + q74 (t)[c]A4(t)

A9(t) = M9(t) + q9,10(t)[c]A10(t), A10(t) = q10,2(t)[c]A2(t) + q10,2(11)(t)[c]A2(t)                                                                              (7-14)

Taking Laplace Transform of eq. (7-14) and solving for  

 = N2(s) / D2(s)                                                                                                                                                                   (15)

Where

N2(s) = (1 -  22(3)(s)) {  0(s) + 01(s) 1(s) + 09(s) 9(s)}+ 2(s){ 01(s) 42(s) + 07(s)72(s) +  73(8)(s)) +  09

(s) 9,10 (s)( 10,2 (s) + 10,2(11)(s))}

D2(s) = (1 -  22(3)(s)) { 1 -  46(5)(s) 60(s) ( 01(s) 44 (s) + 07(s) 74(s))

 - 20(s)01(s)12(s)07(s)(  72(s)) +  72(8)(s) +  09 (s) 9,10 (s)

( 10,2 (s) + 10,2(11)(s))}

The steady state availability

A0 =  =  =

Using L’ Hospitals rule, we get

A0 =  =                                                                                                                                        (16)

Where

N2(0)= p20(0(0) + p011(0) + p09 9(0) ) + 2(0) (p01p12 + p07 (p72

 + p72(8) + p09 ))

D2(0) = p20{  + p01 + (p01 p14 + p07 p74 )+ p07  + p07  + p09()

                + { 1- ((p01p14 + p07 p74 )}

, ,  

The expected up time of the system in (0, t] is

(t) =  So that                                                                                                        (17)

The expected down time of the system in (0, t] is

 (t) = t- (t) So that                                                                                                                     (18)

The expected busy period of the server for repairing the failed unit under extremely high radiations in (0, t]

R0(t) = S0(t) + q01(t)[c]R1(t) + q07(t)[c]R7(t) + q09(t)[c]R9(t)

R1(t) = S1(t) + q12(t)[c]R2(t) + q14(t)[c]R4(t),

R2(t) = q20(t)[c]R0(t) + q22(3)(t)[c]R2(t)

R4(t) = q46(3)(t)[c]R6(t), R6(t) = q60(t)[c]R0(t)

R7(t) = (q72(t)+ q72(8)(t)) [c]R2(t) + q74 (t)[c]R4(t)

R9(t) = S9(t) + q9,10(t)[c]R10(t), R10(t) = q10,2(t) + q10,2(11)(t)[c]R2(t)                                                                                                (19-26)

Taking Laplace Transform of eq. (19-26) and solving for  

  = N3(s) / D2(s)                                                                                                                                                                   (27)

Where

N2(s) = (1 -  22(3)(s)) {  0(s) + 01(s) 1(s) + 09(s) 9(s)} and D2(s) is already defined.

In the long run, R0 =                                                                                                                                                       (28)

where N3(0)= p20(0(0) + p011(0) + p09 9(0) ) and D2(0) is already defined.

The expected period of the system under extremely high radiations in (0, t] is

(t) =  So that

The expected Busy period of the server for repair of dissimilar units by the repairman in (0, t]

B0(t) = q01(t)[c]B1(t) + q07(t)[c]B7(t) + q09(t)[c]B9(t)

B1(t) = q12(t)[c]B2(t) + q14(t)[c]B4(t), B2(t) = q20(t)[c] B0(t) + q22(3)(t)[c]B2(t)

B4(t) = T4 (t)+ q46(3)(t)[c]B6(t), B6(t) = T6 (t)+ q60(t)[c]B0(t)

B7(t) = (q72(t)+ q72(8)(t)) [c]B2(t) + q74 (t)[c]B4(t)

B9(t) = q9,10(t)[c]B10(t), B10(t) = T10 (t)+ (q10,2(t) + q10,2(11)(t)[c]B2(t)                                                                                              (29- 36)

Taking Laplace Transform of eq. (29-36) and solving for 

 = N4(s) / D2(s)                                                                                                                                                                   (37)

Where

N4(s) = (1 -  22(3)(s)) { 01(s)14(s) 4(s) + 46 (5)(s) 6(s)) + 07(3)(s) 74(s)( 4(s)

 +  46(5)(s) 6(s))+ 09(s)09,10(s)  10(s) )

And D2(s) is already defined.

In steady state, B0 =                                                                                                                                                      (38)

where N4(0)= p20 {( p01 p14 + p07 p74) (4(0) +6(0)) + p09 10(0) } and D2(0) is already defined.

 The expected busy period of the server for repair in (0, t] is

(t) =  So that                                                                                                                        (39)

The expected Busy period of the server for repair of unit for failure due non-availability of heavy water in nuclear power reactors in (o, t]

P0(t) = q01(t)[c]P1(t) + q07(t)[c]P7(t) + q09(t)[c]P9(t)

P1(t) = q12(t)[c]P2(t) + q14(t)[c]P4(t), P2(t) = q20(t)[c]P0(t) + q22(3)(t)[c]P2(t)

P4(t) = q46(3)(t)[c]P6(t), P6(t) = q60(t)[c]P0(t)

P7(t) = L7(t)+ (q72(t)+ q72(8)(t)) [c]P2(t) + q74 (t)[c]P4(t)

P9(t) = q9,10(t)[c]P10(t), P10(t) = (q10,2(t) + q10,2(11)(t))[c]P2(t)                                                                                                           (40-47)

Taking Laplace Transform of eq. (40-47) and solving for

 = N5(s) / D2(s)                                                                                                                                                                    (48)

where N2(s) = 07(s )  7(s) 1 -  22(3)(s)) and D2(s) is defined earlier.

In the long run, P0 =                                                                                                                                                       (49)

where N5(0)= p20 p07 4(0) and D2(0) is already defined.

 The expected busy period of the server for repair of the in (0, t] is

(t) =  So that                                                                                                                          (50)

The expected number of visits by the repairman for repairing the different units in (0, t]

H0(t) = Q01(t)[c]H1(t) + Q07(t)[c]H7(t) + Q09(t)[c]H9(t)

H1(t) = Q12(t)[c][1+H2(t)] + Q14(t)[c][1+H4(t)], H2(t) = Q20(t)[c]H0(t) + Q22(3)(t)[c]H2(t)

H4(t) = Q46(3)(t)[c]H6(t), H6(t) = Q60(t)[c]H0(t)

H7(t) = (Q72(t)+ Q72(8)(t)) [c]H2(t) + Q74 (t)[c]H4(t)

H9(t) = Q9,10(t)[c][1+H10(t)], H10(t) = (Q10,2(t)[c] + Q10,2(11)(t))[c]H2(t)                                                                          (51-58)

Taking Laplace Transform of eq. (51-58) and solving for  

 = N6(s) / D3(s)                                                                                                                                                                    (59)

Where

N6(s) = (1 –  22(3)*(s)) { 01(s)12(s)14(s))  09 (s) 9,10 (s)}

D3(s) = (1 -  22(3)*(s)) { 1 - (01(s) 14 (s) + 07(s) 74(s))46(5)*(s) 60(s)} - 20(s)01(s)12(s)07(s)

(72(s)) +  72(8)(s) +

09 (s)9,10 (s) ( 10,2 (s) +Q 10,2(11)*(s))}

In the long run, H0 =                                                                                                                                                      (60)

where N6(0)= p20 (p01 + p09) and D’3(0) is already defined.

The expected number of visits by the repairman for repairing the unit for failure due non-availability of heavy water in nuclear power reactors in (0, t]

V0(t) = Q01(t)[c]V1(t) + Q07(t)[c]V7(t) + Q09(t)[c]V9(t)

V1(t) = Q12(t)[c]V2(t) + Q14(t)[c]V4(t), V2(t) = Q20(t)[c]V0(t) + Q22(3)(t)[c]V2(t)

V4(t) = Q46(3)(t)[c]V6(t), V6(t) = Q60(t)[c]V0(t)

V7(t) = (Q72(t)[1+V2(t)]+ Q72(8)(t)) [c]V2(t) + Q74 (t)[c]V4(t)

V9(t) = Q9,10(t)[c]V10(t), V10(t) = (Q10,2(t) + Q10,2(11)(t))[c]V2(t)                                                                                       (61-68)

Taking Laplace-Stieltjes transform of eq. (61-68) and solving for  

 = N7(s) / D4(s)                                                                                                                                                                  (69)

where N7(s) =  07 (s) 72 (s) (1 –  22(3)*(s)) and D4(s) is the same as D3(s)

In the long run, V0 =                                                                                                                                                      (70)

where N7(0)= p20 p07 p72 and D’3(0) is already defined.

 

COST BENEFIT ANALYSIS

The cost-benefit function of the system considering mean up-time, expected busy period of the system under extremely high radiations when the units stops automatically, expected busy period of the server for repair of unit for failure due non-availability of heavy water in nuclear power reactors, expected number of visits by the repairman for unit failure, expected number of visits by the repairman for failure due non-availability of heavy water in nuclear power reactors. The expected total cost-benefit incurred in (0, t] is

C (t) = Expected total revenue in (0, t]

  • expected total repair cost for unit failure due non-availability of heavy water in nuclear power reactors in (0,t] - expected total repair cost for repairing the units in (0,t ]
  • expected busy period of the system under extremely high radiations when the units automatically stop in (0,t]
  • expected number of visits by the repairman for repairing the unit failure due non-availability of heavy water in nuclear power reactors in (0,t]
  • expected number of visits by the repairman for repairing of the units in (0,t]

The expected total cost per unit time in steady state is

C = =

 = K1A0 - K2P0 - K3B0 - K4R0 - K5V0 - K6H0

Where

K1: revenue per unit up-time,

K2: cost per unit time for which the system is under repair failure due non-availability of heavy water in nuclear power

K3: cost per unit time for which the system is under unit repair

K4: when units automatically stop cost per unit time for which the system is under extremely high radiations

K5: cost per visit by the repairman for which unit under repair for failure due non-availability of heavy water in nuclear power reactors

K6: cost per visit by the repairman for units repair.

 

CONCLUSION

After studying the system, we have analyzed graphically that when the failure rate due to non-availability of heavy water in nuclear power reactors, failure rate due to extremely high radiations increases, the MTSF and steady state availability decreases and the cost function decreased as the failure increases.

 

REFERENCES

    1. Dhillon, B.S. and Natesen, J, Stochastic Anaysis of outdoor Power Systems in fluctuating environment, Microelectron. Reliab.. 1983; 23, 867-881.
    2. Goel, L.R., Sharma,G.C. and Gupta, Rakesh Cost Analysis of a Two-Unit standby system with different weather conditions, Microelectron. Reliab., 1985; 25, 665-659.
    3. Barlow, R.E. and Proschan, F., Mathematical theory of Reliability, 1965; John Wiley, New York.
    4. Gnedanke, B.V., Belyayar, Yu.K. and Soloyer, A.D., Mathematical Methods of Relability Theory, 1969 ; Academic Press, New York.
    5. Kan, Cheng, Reliability analysis of a system in a randomly changing environment, Acta Math. Appl. Sin. 1985, 2, pp.219-228.
    6. Cao, Jinhua, Stochastic Behavior of a Man Machine System operating under changing environment subject to a Markov Process with two states, Microelectron. Reliab., 1989; 28, pp. 373-378.

     

     

 
 
 
 
 
  Copyrights statperson consultancy www

Copyrights statperson consultancy www.statperson.com  2013. All Rights Reserved.

Developer Details