Next Article in Journal
Potential Benefit of the Quadbar™ on All-Terrain Vehicles
Next Article in Special Issue
Evaluating the Specificity of Community Injury Hospitalization Data over Time
Previous Article in Journal
All-Terrain Vehicle Safety―Potential Effectiveness of the Quadbar as a Crush Prevention Device
Previous Article in Special Issue
Influence of Cognitive Biases in Distorting Decision Making and Leading to Critical Unfavorable Incidents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Damage Survivability of Passenger Ships—Re-Engineering the Safety Factor

by
Jakub Cichowicz
1,*,
Nikolaos Tsakalakis
2,
Dracos Vassalos
3 and
Andrzej Jasionowski
1
1
Brookes Bell Safety at Sea, 280 St Vincent Street, Glasgow G2 5RL, UK
2
Presently at Alpha Marine Services, 55 Kastoros Str., Piraeus 18545, Greece
3
The Ship Stability Research Centre (SSRC), University of Strathclyde, 100 Montrose Street Glasgow G4 0LZ, UK
*
Author to whom correspondence should be addressed.
Submission received: 15 August 2015 / Revised: 25 January 2016 / Accepted: 29 January 2016 / Published: 19 February 2016
(This article belongs to the Special Issue Feature Papers)

Abstract

:
This paper presents a brief summary of the work carried out by SSRC and Brookes Bell Safety at Sea within the EC-funded project GOALDS on the development of a new formulation for assessing the survivability of damaged ships in waves. The proposed formula is meant to be an alternative or replacement to the s-factor in use within the current SOLAS regulations for probabilistic damage stability. The authors briefly discuss concerns related to the current survivability model and present the process of development that led to the re-engineered formulation.

1. Introduction

The current measure of survivability of a damaged ship in a seaway has been shown to be inaccurate and inadequate practically from the moment of its introduction within the probabilistic framework of damage stability, known commonly as SOLAS 2009 [1]. Some research indicated existence of a safety gap between the new SOLAS and SOLAS90 + Stockholm Agreement (SA [2]), e.g., [3] where in conclusions, the authors state that “Applied to two new RoPax designs, it has been revealed that the safety level provided by the new SOLAS 2009 rules is significantly lower compared to the old damage stability rules SOLAS 90 in conjunction with the Stockholm Agreement.”, while other suggested that the new regulations were more demanding (e.g., in terms of GM requirements than the SOLAS90 + Stockholm Agreement) for large RoPax ships while being generally less stringent for small ferries [4]. The latter observation was later “extrapolated” to all passenger vessels and common understanding was that the formula underestimates actual safety levels of large cruise vessels (the HARDER [5], sample of ships included two RoPax ships and one cruise liner) while the safety performance of RoPax ships may be seriously overrated [6].
From this perspective, given the controversies associated with the s-factor it is important to make a clear distinction between findings of and subsequent recommendations by the project HARDER and the formula adopted by the IMO. The brief summary of the outcome of the project HARDER can be found in [7] and there an attentive reader can readily see that the IMO adoption deviated in significant details from the project’s proposals (e.g., the authors recommended use of SEM—static equivalent method [8]—accounting for effects associated with accumulation of floodwater on the vehicle deck of Ro-Ro vessels, either in place or in conjunction with conventional, GZ-based model). Therefore, while considering the shortcomings of the present formula, the reader should not undermine the quality of the HARDER research. On the contrary, the outcome of the project had such gravity that the IMO decided to use it as a basis for the new s-factor formulation, despite the fact that the HARDER’s regression models were derived for a rather modest sample. The historical background behind the development of “harmonised” regulations and concerns related to the adopted survivability model are discussed in detail in [9].
In the course of more recent research activities it has been shown that the issue is much more complex and it is generally a consequence of a few, yet critical, ill-founded decisions made at the time of the s-factor development and its adoption by IMO. As far as the project HARDER is concerned it seems that its most significant shortcoming was the fact that it did not emphasize strongly enough that although the survivability model seemed “satisfactory” at the time of development the evidence supporting it was scarce and insufficient (small, largely homogenous ship sample and incomplete parametric investigation), and that more research was needed in order to conclude the work. Irrespective of this, given that the HARDER’s proposal resulted in a dramatic change in the definition of the factor “s” (as compared with A.265 [10], where the survivability model was based on more intuitive parameters and concepts such as residual GM or freeboard) the decision by IMO to adopt the “conventional” methodology discussed within the project without further studies can be perceived as a serious oversight.
This, combined with some basic notional inconsistencies related to the definition of survivability, led to the situation that the s-factor is biased towards the specific sample. Furthermore, it does not account for certain built-in stability-enhancing features above the deck which practically prohibits harmonisation of the stability framework with relevant goal-based regulations like the Safe Return to Port (SRtP). Finally, the current formulation has been often criticized, e.g., [9] due to the fact that it does not signify key attributes of a damaged ship, such as residual freeboard, accumulation of floodwater on deck [11], etc.
As a result, the s-factor formulation has been perceived as a flawed, inaccurate and confusing instrument that may complicate the design process, misguide the designer and even compromise safety [4].
On one hand (as earlier indications suggested and recent research confirms) the flaws within the s-factor have a minor impact on the attained index of subdivision (the A-Index - an expected (weighted average) probability of survival of any feasible flooding scenario in three distinct loading conditions) of contemporary designs (as the marginal survivability cases, i.e., damage cases with a probability of survival between 0 and 1, constitute only a small percentage of all feasible flooding casualties; this in fact was part of the reason why the current s-factor formulation was adopted in SOLAS 2009). On the other hand, however, the present s-factor formulation is likely to fail precisely where it is particularly needed, that is in those scenarios that are likely to happen and can potentially result in catastrophic losses. Furthermore, these are usually complex cases, difficult to deal with, resulting from design limitations and/or conflicting objectives and in these cases the s-factor fails at least twice. Firstly, this is by failing to provide an accurate assessment of survivability (particularly for the ships of watertight arrangements deviating significantly from that of HARDER’s sample) in dynamic environment. Secondly, it fails because the GZmax and residual range of positive stability cannot be readily linked to the easily measurable design attributes, such as freeboard or compartment length (the relationship between these parameters and the GZ curve exists through several other parameters), which could be used by a designer to resolve the technical problems at hand. Comparison between distinct designs purely on a basis of maximum righting lever and a range of positive stability is neither easy nor complete and in this respect the present s-factor formula may prevent safe and rational design.
The EU-funded project GOALDS (Goal Based Design for Damage Stability, 2009–2012 [12]) set the review of the existing s-factor as one if its core objectives (along with the development of equivalent safety factor for grounding damages and harmonisation of damage stability requirements with the goal-based regulatory framework).
In the following the authors present briefly the work carried within GOALDS and leading to formulation of a new proposal for the factor s.

2. Background Concepts

The three, most fundamental from the perspective of damaged ship survivability, notions are associated with modes of loss, the capsize band and accumulation of water on deck. These concepts were introduced by earlier research and in GOALDS they were examined specifically from the point of view of methodology and generalisation of the results.

2.1. Modes of Ship Loss

In the most general case there are three distinct mechanisms or “modes” of ship loss in a damaged condition: capsize, sinking and transient capsize. The first two are to be understood as gradual processes of diminishing residual stability or residual buoyancy as a direct result of the flooding process. The mechanism of capsize or sinking is a stochastic (but not ergodic) process resulting in quasi-static loss. The third mode is a dynamic process that can be seen as a response of an intact ship to the rapid floodwater ingress. Given that in a transient mode the potential “memory effects” should diminish within a relatively short time (a matter of seconds) for all cases reaching equilibrium post-damage the action of waves can be neglected in this mode of loss.
In the following, the term capsize shall be used to denote any kind of ship loss unless stated otherwise.

2.2. Capsize Band

Capsize band is a concept describing range of sea-states within which a transition from unlikely (lower boundary) to certain (upper boundary) capsizes is observed (in the light of the considerations presented in the previous section, the term “capsize” can be generalized to account for both capsize and sinking).
For a finite observation time, the probability of capsize can be approximated with the help of a sigmoid function [13] of the form:
p f = p f ( H S ) = e H S x 0 Δ x 1 + e H S x 0 Δ x
where: pf is probability of capsizing, HS is the sea-state in question, x0 and Δx are regression parameters—the abscissa of the sigmoid’s inflection point and the bandwidth parameter, respectively.
Alternatively to the above (so called Bolzman’s) representation, earlier publications [14] demonstrated use of normal (Gaussian) representation, where the probability of capsize is approximated by an integral curve of the normal distribution, where x0 would become modal value, μ, of the dpf/dHS distribution and Δx would be replaced by standard deviation, σ.
A significant feature of the capsize band is that as the time of observation increases the capsize band contracts towards its lower boundary, becoming a unit step function as time of exposure to action of environment approaches infinity (see Figure 1).
This property is of major importance and, as will be shown in the following, it had become one of the key findings made during the re-engineering of the s-factor.
It is worth noting that the concept of the critical sea state, HScrit, as used for the s-factor in force, is associated with the sea state at which the probability of capsize (pf) is equal to 0.5, based on half hour tests.

2.3. Water on Deck

Accumulation of water on vehicle deck of a RoRo ships (or more generally accumulation of floodwater inside a damaged ship) has always been considered as a prime factor affecting damaged ship survivability (the SEM approach was derived based on the observed correlation between the amount of floodwater accumulated on deck and survivability of the vessel and the loss of m/v Estonia gave the IMO an impetus for introducing regional requirements for RoPax ships operating in NW Europe, known as the Stockholm Agreement [2]). The problem of floodwater accumulation is usually referred to as the Water on Deck (or WoD) problem.
The major difficulty in addressing the survivability from the WoD angle is the stochastic and non-ergodic character of the flooding process, which generally makes analysis very complex and time consuming (necessity of analysing very long time histories). Proper analysis of flooding history is difficult even for numerical tests and physical experiments suffer additionally due to short time realisations, complex set-up and inherent uncertainty of measurements.
In this context properties of the capsize band came to aid by hinting an unconventional approach for processing WoD time series: the analysis of WoD accumulation should be based not on the cases within the capsize band (i.e., capsizing cases) but just outside its lower boundary, towards which the band contracts (i.e., cases likely to survive for “infinite” time). Such an approach has certain advantages—firstly, long but finite observation time around the lower boundary practically guarantees infinite survival time (the lower boundary can be assumed time invariant, as previously shown) and secondly, analysis of surviving cases allows relaxing the limitation with respect to ergodicity of the process—i.e., flooding, similar to ship response, can be assumed (for engineering purposes) to be ergodic.
The above reasoning allowed formulating a procedure for analysing WoD (or floodwater accumulation in the general case) based on an averaged 95th percentile of floodwater mass (derived from numerical simulations) calculated within the cumulative time (i.e., measured at intervals equal to the wave period).
The main shortcoming of such approach was the difficulty of the inclusion of the transient phase of the time history (non-stationary) but the main purpose was the visualisation and qualitative assessment of the flooding processes without the need to withdraw any quantitative conclusions with respect to the critical amount of floodwater.
The analysis, in spite of certain mathematical deficiencies, conveyed very important information. Firstly, it has been noticed that the 95th percentile curves were asymptotic, i.e., all the surviving cases were limited by a horizontal asymptote whilst a slant asymptote indicated underlying progressive flooding; the latter cases would often result in a loss.
Secondly, although the limiting amount of floodwater would vary with HS, such variability proved to be statistically insignificant (i.e., individual realisations would generally fall within confidence bands of the realisations of sea state corresponding to the lower boundary of the capsize band, e.g., Figure 2).
Furthermore, in most of the cases, should the floodwater characteristic of any particular realisation exceeded (by some statistically significant amount) the upper confidence limit calculated for the lower boundary of the capsize band this would be a clear indication of subsequent loss, see Figure 3.
The final outcome of the WoD analysis was that the flooding process of the surviving cases could be characterised by a statistically unique (for any given damage) limiting mass of floodwater, independent of sea-state and time of exposure. It is important to note that the limiting amount of floodwater coincides with the quantity obtained for the highest sea-state at which no losses were observed—the lower boundary of the capsize band. Increasing the sea-state further does not have an impact on the limiting value but it simply increases the probability of exceeding it (i.e., the higher is the sea state the higher is the probability that the limit amount of floodwater will be exceeded and the ship will capsize). This probability is dependent on time of exposure—in the extreme case of infinite observation time, all realisations in sea-states higher than the lower capsize boundary would at some point experience floodwater accumulation in excess of the critical mass and face inevitable loss (this is well in-line with the unit-step representation of the capsize band).

3. The s-Factor Development

The GOALDS project utilised to a large extent the knowledge generated by the past decades’ research. One of the important lessons learnt was that the key problems of the earlier studies were due to insufficiently diversified sample (in terms of ship types, sizes and internal arrangements) and that the parameters used for regression models were often selected in an arbitrary fashion. In order to avoid similar issues a decision was made to base the development primarily on the numerical results and use the physical experiments mainly for the model verification. This combined with reasonably diversified sample allowed for generating dataset large enough for the systematic parametric investigation.

3.1. Physical and Numerical Survivability Tests

In order to generate the dataset large enough for parametric investigation a significant number of physical tests and numerical simulations were performed with use of several ship models operating at different loading conditions. Although all the ships were subdivided according to SOLAS ’90 (pre-SOLAS2009 vessels) special attention was paid to diversification of the sample—ships varied in size and the subdivisions included stability enhancing features above the traditional bulkhead deck (e.g., side-casings) as well as large unsubdivided cargo spaces below the deck (Long Lower Hold). The GOALDS sample was further supplemented by the ships investigated within the HARDER and EMSA projects.
In total the sample consisted of over twenty various RoPax designs and three cruise ships (small, medium and large).

3.2. Probability of Surviving Collision Damage

Before proceeding further with the results of parametric modelling it might be worthwhile to look into the underlying concepts (though not explicitly referenced to in the SOLAS convention) making the s-factor a measure of ship survivability in waves.
Formally the probability of survival following collision, hull breach and flooding is given as [15]:
s = 0 d H S f H S | c o l l ( H S ) F s u r v ( H S )
where: f H S | c o l l ( H S ) is the probability density distribution of sea states expected to be encountered during collision and F s u r v ( H S ) is the probability of survival in that sea state when exposed to a specific flooding case.
It should be noted here that, given that all the tests performed during the s-factor development were limited to 30 min, the probability of survival is in fact a conditional probability:
F s u r v ( H S ) F s u r v ( t = 30 min | H S )
Therefore it yields:
s ( t = 30 min ) = 0 d H S f H S | c o l l ( H S ) F s u r v ( t = 30 min | H S )
Assuming that the probability of survival, F s u r v ( H S ) , can be approximated by a unit step function centered on the sea state the Hscrit constitutes the 50th percentile significant wave height the vessel subjected to particular damage scenario may survive for 30 min (this corresponds to the abscissa of the inflection point of the sigmoid obtained for t = 30 min and is in line with the previous research):
F s u r v ( H S ) = { 1 H S H S c r i t 0 H S > H S c r i t
It should be noted that although replacing the pf distribution with the step function is in this particular case supported with little evidence, it “does the trick” and allows the integration to be avoided with negligible impact on accuracy of the prediction (for as long as the bandwidth of the capsize band is narrow).
Eventually, the final formulation becomes:
s = 0 H S c r i t d H S f H S | c o l l ( H S ) = exp ( exp ( 0.16 1.2 H S c r i t ) )
where the HS crit is given as:
H S c r i t | t = 30 min = 4 ( min ( G Z max , 0.12 ) 0.12 min ( R a n g e , 16 ) 16 ) = 4 s ( t = 30 min ) 4
In the essence, the approach adopted within the GOALDS project is consistent with other research projects on ship survivability, particularly HARDER and EMSA. To this end the only major difference stems from the adoption of the lower capsize boundary as HS crit and subsequent time invariance of the probability of capsize. At this stage it becomes apparent that the postulated flaw in s-factor (i.e., its inaccuracy) must have been related to the improper or incomplete parametric model.

3.3. Study on Dominant Parameters

In spite of the extensive use of numerical simulations the sample size available within the GOALDS project did not allow for performing the sensitivity analysis based solely on advanced techniques like Design of Experiments (DoE). Due to limited resources not all desired parameters (such as damage extent, loading conditions, etc.) could have been varied to the extent required by, for example, the full factorial design of experiments hence the test matrix (i.e., selection of design variables) had been designed based mainly on past experience and critical reasoning on the outcome of the previous research with only minor amendments introduced during the course of tests whereas the DoE methodology was primarily used to derive response surface models of HS crit as a function of design/damage attributes.
The initial step involved investigation of the present SOLAS s-factor in order to conclude whether the GZmax and Range are sufficient to predict accurately the critical significant sea-state for the generated sample. In order to verify this, the HS crit was calculated with use of the underlying SOLAS model (as given by Equation (7)) and with full quadratic response surface model (denoted in the following as DoE) and given as
H S c r i t = f ( x 1 , x 2 ) = a 0 + a 1 x 1 + a 2 x 2 + a 3 x 1 x 2 + a 4 x 1 2 + a 5 x 2 2
where, a i are the regression coefficients and x 1 , x 2 stand for GZmax and Range, respectively.
The results (see Figure 4) showed that both SOLAS and DoE methods failed to accurately predict the HScrit. The use of a relatively complex model slightly improved the quality of fit (the correlation coefficient increased from the initial value of 0.64 to 0.79) but the exercise only proved that GZmax and Range alone are insufficient to estimate the critical sea state.

3.4. Effect of Scale

The analysis of results pertaining to small and large vessels (from the GOALDS sample), revealed that the significant scatter in results can be attributed to the effect of scale (i.e., dimensions of the actual ships). Indeed, one of the major concerns related to the SOLAS’s s-factor was that it does not account for the ship’s size, which was thought toresult in inaccuracies in predicting survivability of vessels deviating significantly in size from the PRR01, used as a basis (68% of HARDER data points [9]) for its derivation.
In addition to this, the SOLAS 2009 s-factor formulation (based solely on characteristics of the residual GZ curve) is limited to relatively small range and maximum GZ values and therefore it tfails to account for the contribution from watertight volume distributed high enough to be “unseen” by the formulation. This deviates from the good Naval Architecture practice previously expressed through the explicit demand for provision of the residual freeboard.
In order to take the size into consideration the next step of the development involved adding to the initial set of parameters a factor that can best be described as a measure of centroid of the residual volume given, as follows:
z V = ( V C B T V C B D T )
where, VCB T and VCB D, are vertical centres of buoyancy of the intact and damaged compartments (intact and residual watertight envelopes), respectively, and T stands for draught in intact condition. Based on this, the experimental data was approximated by the following regression formulae:
f 1 ( G Z max , z V ) = 18.4 ( G Z max z V ) 0.5 f 2 ( R a n g e , z V ) = 1.5 ( R a n g e z V ) 0.9 H S c r i t = f 1 f 2
As a result the formula offered very accurate prediction with 0.93 correlation with the experimental data. Furthermore, applying DoE to the new parameter set led to a correlation coefficient of 0.97—a good indication that the size of the ship was indeed the missing parameter (see Figure 5).
The major concern of using a regression-based formulation is the possibility of over-reliance on the particular dataset used for its derivation and possible unreliability while applying it for a completely different sample (e.g., the innovative ships of the future).
On the other hand, the formula proved that GZ curve related parameters complemented by the geometrical characteristics of a damage (residual volume in some form) are sufficient for predicting the survival sea-state. Furthermore, during the analysis of the experimental and numerical data it had been found that one of the parameters to which the HS crit is particularly sensitive is the residual metacentric height, GMf (GM of the flooded ship). Hence in the final stage an attempt was made to bring stability characteristic (GZ curve properties) and size of the ship into a formula of possible physical significance.

3.5. Final Formulation

Summarising the foregoing, all investigations led to the conclusion that the parameters to be included in the s-factor formulation should include as minimum the GZmax, Range, some measure of the residual volume (VR) and GMf. Moreover, the ensuing formulation, i.e., H S c r i t = f ( G Z max , R a n g e , V R , G M f ) should have some degree of physical significance in order to minimise overreliance on a specific data set. On this basis, an analytical expression was sought using a combination of the above parameters, as presented in the following.
In order to simplify the modelling process it was assumed that a GZ curve can be, within the range of positive stability, approximated by an isosceles triangle of basis equal to Range and height equal to GZmax. In such a representation the following relations would hold:
A G Z = 1 2 G Z max R a n g e G Z max = 1 2 G M f R a n g e 1 2 R a n g e = G Z max G M f
On this basis the following expression was tested as a first step, deriving from simple dimensional analysis:
( H S c r i t V R 1 / 3 ) = ( G Z max G M f ) H S c r i t = G Z max G M f V R 1 / 3
As this expression produced slightly conservative results, the additional parameters were introduced through various algebraic manipulations and tested against the experimental results. Firstly, multiplying the numerator and denominator of the Equation (12) by half-Range, as shown in the Equation (13) resolved the issue of biased (conservative) prediction, which was due to approximating the GZ curve with an isosceles triangle. Furthermore, it was observed that the numerator of the Equation (13) is nothing else as an area of a triangle of basis equal to Range and height equal to GZmax. Hence, in the final expression, the Equation (14), it was replaced by the explicit reference to the GZ area and introducing Range and AGZ, respectively, led to significantly improved accuracy.
H S c r i t = 1 2 G Z max R a n g e 1 2 G M f R a n g e V R 1 / 3
H S c r i t = A G Z 1 2 G M f R a n g e V R 1 / 3
Deriving from the above, the survival factor, s, is given by the following expression:
( A G Z , V R , R a n g e > 0 ) s = exp ( exp ( 0.16 1.2 H S c r i t ) )

4. Validation and Additional Considerations

4.1. Physical Significance

The following figure shows a graphical interpretation of a ship’s step-response to some external moment, represented on the GZ curve as the heeling lever, l c . The underlying assumption is that the ship is floating in calm water (i.e., impact of waves on restoring moment is not considered) and that the external moment is given as M c ( t ) = { 0 , t = 0 ρ g l c , t > 0 . Furthermore, it is that φ ( t = 0 ) = φ 0 , φ ˙ ( t = 0 ) = 0 , l c = c o n s t and ρ , g and stand for water density, gravity acceleration and displacement (volumetric) of the ship.
The angles shown in the Figure 6 below are respectively
Angle of static equilibrium, φ 0
Angle of heel due to static action of the external moment(stable dynamic equilibrium), φ 1
Angle of dynamic heel, φ 2
Angle of unstable equilibrium (angle of capsize), φ c
Angle of vanishing stability, φ max
The condition “a” = “b” indicates equilibrium condition in which work of the heeling moment is equal to the change of ship’s potential energy.
l c φ 0 φ 2 d φ = φ 0 φ 2 l d φ
From the above it follows that the heeling lever l c corresponds to the mean value of the righting lever l within the interval φ 0 , φ 2 , i.e.,:
l c = φ 0 φ 2 l d φ ( φ 2 φ 0 )
Furthermore, the definition of GM implies that:
l c = l ( φ 1 ) l ( φ 0 ) + d l d φ | φ = φ 0 ( φ 1 φ 0 ) = G M f ( φ 1 φ 0 )
Therefore, combining Equations (17) and (18) yields
( φ 1 φ 0 ) = φ 0 φ 2 l d φ G M f ( φ 2 φ 0 )
The above equation shows that the angle of dynamic equilibrium can be approximately given as a ratio of mean righting lever within the φ 0 , φ 2 range and the residual metacentric height, G M f . In GOALDS formulation the RHS of the above equation is calculated with upper limit corresponding to the point of vanishing stability (point of no return), i.e., φ 2 = φ max . In such a case the heeling lever l c equals the mean value of the righting lever within the interval φ 0 , φ max , that is l c = l ¯ = φ 0 φ max l d φ φ max φ 0 and the Equation (19) becomes:
( φ 1 φ 0 ) = φ c r i t = φ 0 φ max l d φ G M f ( φ max φ 0 ) = A G Z G M f R a n g e
The above implies that the GOALDS’ H S c r i t is proportional to work done by the external heeling moment over the interval φ 0 , φ 1 , i.e., H S c r i t l ¯ φ 0 φ 1 d φ :
H S c r i t = 2 3 V R φ 0 φ max l ( φ ) d φ G M ( φ max φ 0 ) = 2 ( φ 1 φ 0 ) 3 V R = 2 3 V R φ 0 φ 1 d φ = α l ¯ φ 0 φ 1 d φ
where α can be calculated as
α = 2 3 V R l ¯ = 2 3 V R ( φ max φ 0 ) φ 0 φ max l d φ = 2 3 V R R a n g e A G Z
Finally, in order to satisfy the condition that for the specific l ¯ the resultant angle of dynamic heel φ 2 must not exceed the angle of unstable equilibrium, φ c (angle of capsize), i.e., it must be that b max a and therefore:
φ 1 φ c l d φ l ¯ φ 1 φ c d φ l ¯ φ 0 φ 1 d φ φ 1 φ 1 l d φ
After rearranging the condition, φ 2 φ c can be written as
l ¯ φ 0 φ c l d φ φ c φ 0
That is, the average value of the righting lever over the interval φ 0 , φ max , l ¯ must not exceed the average value of the righting lever within the range φ 0 , φ c . This implies that if the condition is not met the critical significant wave height, H S c r i t , could be calculated as follows:
H S c r i t = φ 0 φ c l d φ 1 2 G M f ( φ c φ 0 ) V R 1 3

4.2. Accounting for the Presence of Openings

In case of present unprotected openings (e.g., semi-watertight doors, down-flooding points etc.) within the range of the residual GZ curve, the critical significant wave height should be derived from the following expression:
H S c r i t = 2 ( φ D φ 0 ) V R 1 / 3
where AGZ E is an effective area under the GZ curve taken up to the heel angle, φ D corresponding to the submersion of the opening in question.

4.3. Transient Capsize

As already mentioned, the evidence collected during the initial stages of the development suggested clearly that both modes of the gradual loss, i.e., capsize and sinking, can be handled with the use of the same formulation but in case of a rapid capsize resulting from sudden floodwater ingress, a different formulation is warranted [16]. In this respect, it is proposed to assess vulnerability to transient capsize with the use of the following expression, based on well-known energy balance:
k i = { 1 , i f φ i < E φ D m f i y f i 0 , o t h e r w i s e
where, Eϕ is a dynamic righting lever (integral of the residual GZ curve) up to the heel angle of submersion of the unprotected opening, D is displacement of intact ship, and mfi and yfi are mass and transverse centre of mass of floodwater respectively. The vulnerability check should be performed for all the flooding stages, and the final k is simply a product of ki. This results in the following expression:
s = { e ( e ( 0.16 1.2 H S c r i t ) ) , ( A G Z , V R , R a n g e , k > 0 ) 0 , o t h e r w i s e
The binary nature of k is dictated mainly by the need to highlight those damage cases that require special attention and unless a detailed investigation is performed they should not contribute to the A-Index summation.

4.4. Time to Capsize

In sea states exceeding the critical significant wave height, the probability of survival and time to capsize decrease, the first following a sigmoid pattern [13,17], the latter according to a hyperbolic manner as graphically explained in Figure 7 [17].
The hyperbolic behaviour of the TTC can be captured by the following formulae.
T T C = a H S H S c r i t [ min ]
a = 3 H S c r i t 1.4
Following a systematic investigation of the acquired data, the observation was made that parameter “a” appearing in formulation (29) can also be linked with HsCRIT as shown in Figure 8, where curves for various ships and damages have been plotted together (each point in the Figure 8 corresponds to the average TTC of at least 10 realisations of each sea-state). Figure 9 depicts this relation for various experimental data, derived by means of numerical simulations. The continuous line is the regression model which can be found in the Equation (30).
Figure 10 below shows impact of observation time on the behaviour of the model (single damage case).

4.5. Validation of the Proposed S-Factor

The first stage of validation aimed at verifying whether the chosen parameter set is sufficient for exhaustive description of the dataset. This was done by applying design of experiments (DoE) to the experimental data. The outcome from fully quadratic model constructed accounting for A G Z , G M f and R a n g e resulted in 0.99 correlation with the experimental data, which was considered a satisfactory test concerning completeness of the parameter set (the correlation coefficient calculated for the GOALDS prediction was about 0.90, see Figure 11).
Another important aspect of the verification process concerned impact of ship’s size on accuracy of the prediction and as the Figure 12 shows clearly, GOALDS formula accounts well for size of the ship with the prediction errors of similar magnitudes across all categories. On the other hand, the HARDER formula demonstrates significant variations in prediction errors depending on ship size. In case of GOALDS prediction, 100% of sample points fell within ±1.0 m from the measured H S c r i t compared with much broader interval for harder data (−2.0/+2.5). The response surface model resulted in all the sample points falling within the ±0.5 m interval
In the further stage, A-Index calculated by means of SOLAS and GOALDS survival factors have been compared against each other and against numerical simulations using Monte Carlo sampling (hereafter referred to as performance based assessment, PBA). These results are shown in the figure below:
The impact of the s-factor formulations on the A-index (i.e., average of the probability of survival) was expected to be small due to the limited (in terms of probability) share of the marginal cases of s = (0;1). As it can be seen from the above figure, Figure 13, the discrepancies are general within a few percent, with the new s-factor giving generally higher (than SOLAS) values for all but the small RoPax (R2). The highest difference between current and new s-factors is shown in the case of the very large cruise vessel (C3) for which PBA predicts an average survivability of about 98%. It should be borne in mind, however, that in spite of the small discrepancies between SOLAS 2009 and GOALDS predictions, the results are not conclusive and should be approached with caution. Specifically, it has been demonstrated in the foregoing that the SOLAS s-factor exhibits much lower accuracy than its GOALDS counterpart, which as shown for the tested RoPax vessels resulted in 4.5 m confidence interval for H S c r i t i.e., larger than the entire band of sea states relevant in the s-factor calculations (sea states of H S 0–4 m, as determined by the probability of specific sea state in the event of collision).
The PBA results are generally within a few percentage points from the A-index calculations with the larger discrepancies in the two extreme cases of small RoPax (R2) and the largest cruise ship (C3). Once again, given the inherent uncertainty of the PBA (due to mathematical modelling, duration and number of performed tests) these results cannot be used to draw conclusions on the impact of the s-factor on the overall survivability as measured by the A-index.
Finally, neither the SOLAS nor GOALDS formulation could be sufficiently verified against experimental data for cruise ship models. Due to the complexity of the physical model tests of survivability, the GOALDS data appropriate for comparison are very limited and consist of only two valid experimental points for C1. These points comprise two runs with semi-watertight doors open and closed performed in waves, starting from the damage equilibrium position. These are presented in the Figure 14 below.

5. Concluding Remarks

This paper details all the steps taken during the GOALDS project and aiming at re-formulating the s-factor. The new formulation has all the characteristics that one would intuitively expect, in particular:
The formulation is simple, readily calculable and rational with sound and intuitive physical significance, well in line with other stability concepts (e.g., weather criterion).
It is very much along the lines of the current formulation but uses as its basis the Hs critical for infinite survivability, consistent with the SRtP philosophy. This way water on deck, though not included explicitly in the formulation, is accounted for.
It accounts for scale, as we all suspected it should.
It demonstrates a high degree of correlation with all available results (particularly experimental), using HARDER, GOALDS and EMSA projects data.
It accounts for any level and mode of subdivision, including the watertight envelope above the traditional bulkhead deck.
Finally, and most importantly, it is a generic formulation deriving from basic principles and hence expected to be applicable to all passenger ships, from current to future innovations.
Finally, although the formulation has been only partly validated, it is believed that given all the evidence presented in the foregoing, the new proposal should be taken as a serious alternative to the present s-factor as a design tool and a measure of survivability within the regulatory framework.

Acknowledgments

The authors would like to express their sincere thanks and appreciation to all their colleagues at SSRC and in Project GOALDS who contributed to this breakthrough for their valuable suggestions, discussions and contributions. The financial support from the ECV, EMSA and UK MCA is gratefully acknowledged. Yu Bi, a researcher at the Ship Stability Research Centre and a PhD candidate at the University of Strathclyde, is also acknowledged for her kind support in reviewing and preparing this manuscript for publication.

Disclaimer

The views expressed in this paper are those of the authors and do not necessarily reflect those of the organizations represented by the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IMO. International Convention for the Safety of Life at Sea (SOLAS); IMO: London, UK, 1974. [Google Scholar]
  2. IMO. Resolution 14 Agreement Concerning Specific Stability Requirements for Ro-Ro Passenger Ships Undertaking Regular Scheduled International Voyages Between or To or From Designated Ports in North West Europe and the Baltic Sea; IMO: London, UK, 1996. [Google Scholar]
  3. Dankowski, H.; Krüger, S. On the safety levels of the SOLAS 2009 damage stability rules for RoPax vessels. In Proceedings of the 11th International Symposium on Practical Design of Ships (PRADS), Rio de Janeiro, Brasil, 19–24 September 2010.
  4. Sweden and UK. IMO SLF 52/11/1: Damage Stability Regulations for Ro-Ro Passenger Ships Report of the SDS Correspondece Group; IMO: London, UK, 2009. [Google Scholar]
  5. HARDER (1999–2003): “Harmonisation of Rules and Design Rationale”. EC funded project, DG XII-BRITE. 2000–2003.
  6. Papanikolaou, A.; Mains, C.; Rusaas, S.; Szalek, R.; Tsakalakis, N.; Vassalos, D.; Zaraphonitis, G. GOALDS—Goal Based Damage Stability. In Proceedings of the 11th International Workshop on Stability of Ships and Ocean Vehicles, Wageningen, Netherlands, 21–23 June 2010.
  7. Tagg, R.; Tuzcu, C. A Performance-based Assessment of the Survival of Damaged Ships—Final Outcome of the EU Research Project HARDER. In Proceedings of the 6th International Ship Stability Workshop; Webb Institute: New York, NY, USA, 2002. [Google Scholar]
  8. Vassalos, D.; Pawłowski, M.; Turan, O. A Theoretical Investigation on the Capsizal Resistance of Passenger/RoRo Vessels and Proposal of Survival Criteria; Joint Northwest European Project; Marine and Coastguard Agency: Southampton, UK, 1996. [Google Scholar]
  9. Scott, A.L. Damage Stability of Ro-Pax Ships with Water-on-Deck. In Proceedings of the 11th International Workshop on Stability of Ships and Ocean Vehicles, Wageningen, The Netherlands, 21–23 June 2010.
  10. IMO. Resolution A.265(VIII): Regulations on Subdivision and Stability of Passenger Ships as an Equivalent to Part B of Chapter II of the International Convention for the Safety of Life at Sea; IMO: London, UK, 1973. [Google Scholar]
  11. Tagg, R. Comparison of survivability between SOLAS 90/95 and SOLAS 2009 ships—A retrospective view 10 years on from project HARDER. In Proceedings of the 14th International Ship Stability Workshop, Brest, France, 29 September–01 October 2013.
  12. GOALDS (2009-2012): “GOAL Based Damage Stability”. EC funded project, 233876. 2009–2012.
  13. Tsakalakis, N.; Cichowicz, J.; Vassalos, D. The concept of the capsize band revisited. In Proceedings of the 11th International Workshop on Stability of Ships and Ocean Vehicles, Wageningen, The Netherlands, 2010.
  14. Jasionowski, A.; Vassalos, D.; Scott, A. Ship Vulnerability to Flooding. In Proceedings of the 3rd International Maritime Conference on Design for Safety, Berkeley, CA, USA, 26–28 September 2007.
  15. Jasionowski, A. Study on the Specific Damage Stability Parameters of Ro-Ro Passenger Vessels according to SOLAS 2009 Including Water on Deck Calculation, EMSA/OP/08/2009; European Commission: Brussels, Belgium, 2009. [Google Scholar]
  16. Vassalos, D.; Ikeda, Y.; Jasionowski, A.; Kuroda, T. Transient Flooding on Large Passenger Ships. In Proceedings of the 7th International Ship Stability Workshop, Shanghai, China, 1–3 November 2004.
  17. Jasionowski, A.; Dodworth, K.; Vassalos, D. Proposal for Passenger Survival-Based Criteria for Ro-Ro Vessels, International Shipbuilding Progress; IOS Press: Amsterdam, The Netherlands, 1999; Volume 46, No. 448; pp. 473–487. [Google Scholar]
Figure 1. Assumed (top) and observed (bottom) behaviour of capsize band as function of the observation time.
Figure 1. Assumed (top) and observed (bottom) behaviour of capsize band as function of the observation time.
Safety 02 00004 g001aSafety 02 00004 g001b
Figure 2. Examples of accumulation of floodwater (total) for small RoPax (top) and cruise vessel (bottom). In the latter case lack of asymptote indicates progressive flooding—in fact, as longer time simulations proved, this particular damage case resulted in loss even in calm water.
Figure 2. Examples of accumulation of floodwater (total) for small RoPax (top) and cruise vessel (bottom). In the latter case lack of asymptote indicates progressive flooding—in fact, as longer time simulations proved, this particular damage case resulted in loss even in calm water.
Safety 02 00004 g002
Figure 3. Capsizing realisations plotted against averaged floodwater accumulation characteristics for surviving cases (PRR01).
Figure 3. Capsizing realisations plotted against averaged floodwater accumulation characteristics for surviving cases (PRR01).
Safety 02 00004 g003
Figure 4. Comparison of SOLAS prediction against the DoE regression. Note that SOLAS is meant to predict median sea state and therefore it should generally “over-predict” the experimental HS corresponding to the lower limit of the capsize band.
Figure 4. Comparison of SOLAS prediction against the DoE regression. Note that SOLAS is meant to predict median sea state and therefore it should generally “over-predict” the experimental HS corresponding to the lower limit of the capsize band.
Safety 02 00004 g004
Figure 5. Comparison of the predicted and experimental HS crit.
Figure 5. Comparison of the predicted and experimental HS crit.
Safety 02 00004 g005
Figure 6. Static and dynamic equilibiria in a step-response to an external heeling moment.
Figure 6. Static and dynamic equilibiria in a step-response to an external heeling moment.
Safety 02 00004 g006
Figure 7. The concept of time to capsize.
Figure 7. The concept of time to capsize.
Safety 02 00004 g007
Figure 8. TTC curves for 3 different KG values for a specific damage scenario.
Figure 8. TTC curves for 3 different KG values for a specific damage scenario.
Safety 02 00004 g008
Figure 9. Parameter “a” vs. HsCRIT.
Figure 9. Parameter “a” vs. HsCRIT.
Safety 02 00004 g009
Figure 10. TTC for different simulation time.
Figure 10. TTC for different simulation time.
Safety 02 00004 g010
Figure 11. GOALDS and response surface model predictions plotted against measured data.
Figure 11. GOALDS and response surface model predictions plotted against measured data.
Safety 02 00004 g011
Figure 12. Impact of ship size on accuracy of GOALDS and HARDER predictions.
Figure 12. Impact of ship size on accuracy of GOALDS and HARDER predictions.
Safety 02 00004 g012
Figure 13. Impact of GOALDS survival factor on A-Index—comparison of partial index for DS. The letter C denotes cruise vessel while R stands for RoPax. PBA is the performance-based assessment using numerical simulations and damages generated with Monte Carlo sampling.
Figure 13. Impact of GOALDS survival factor on A-Index—comparison of partial index for DS. The letter C denotes cruise vessel while R stands for RoPax. PBA is the performance-based assessment using numerical simulations and damages generated with Monte Carlo sampling.
Safety 02 00004 g013
Figure 14. Comparison of predicted and measured HS crit for the cruise vessel. Note that the tests with SWD closed had been terminated at HS 4.0 m.
Figure 14. Comparison of predicted and measured HS crit for the cruise vessel. Note that the tests with SWD closed had been terminated at HS 4.0 m.
Safety 02 00004 g014

Share and Cite

MDPI and ACS Style

Cichowicz, J.; Tsakalakis, N.; Vassalos, D.; Jasionowski, A. Damage Survivability of Passenger Ships—Re-Engineering the Safety Factor. Safety 2016, 2, 4. https://0-doi-org.brum.beds.ac.uk/10.3390/safety2010004

AMA Style

Cichowicz J, Tsakalakis N, Vassalos D, Jasionowski A. Damage Survivability of Passenger Ships—Re-Engineering the Safety Factor. Safety. 2016; 2(1):4. https://0-doi-org.brum.beds.ac.uk/10.3390/safety2010004

Chicago/Turabian Style

Cichowicz, Jakub, Nikolaos Tsakalakis, Dracos Vassalos, and Andrzej Jasionowski. 2016. "Damage Survivability of Passenger Ships—Re-Engineering the Safety Factor" Safety 2, no. 1: 4. https://0-doi-org.brum.beds.ac.uk/10.3390/safety2010004

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop