Application of Extended Lattice Gas Automata Model in Ship Evacuation Simulation

Lattice gas automaton is a mathematical model that is used to simulate the horizontal uniform evacuation behaviour of a group. However, extended lattice gas automata model is proposed to examine marine evacuation behaviour, which is subject to deck heeling. ,e application of distance accumulation algorithm and the conversion probability mostly make the extended model, while the approach deals with the most complicated ship evacuation. Moreover, the suggested model is expected to enhance the safety and efficiency of evacuation. ,e distance accumulation lattice gas automata model considers multiple movement behaviours, flow density, deck heeling, counterflow, and congestion. Movement behaviour will be severely affected in deck heeling process where people may walk normally, walk while bent over, or crawl. To verify the proposed model, 11 test scenarios and several emergency evacuation scenarios are demonstrated. ,e simulation results explain the validity of another experimental model. ,e number of people in counterflow, deck heeling, and difference in movement have a direct effect on evacuation, which is as discussed in results.,is research article provides a brief study on ship design and crew response behaviour in case of mishap/ accident.


Introduction
Ships are having complex internal structure, narrow space, and limited line of sight which are the same in conditions of highrise buildings [1,2]. However, fire accidents happen in cabin where people evacuate by walking, bending over, or crawling to avoid toxic smoke [3,4]. e process of people evacuation can be made possible in different postures according to individual factors. is happens in framework heels at a very large angle. For example, at 8 : 48 on April 16, 2014, the Motor Vessel (MV) Sewol heeled with sharp turn, while at 8 : 52, the heeling angle was estimated up to 30°to the port side. During evacuation, most of the crew usually grab nearby items while moving out through bending over, walking, and crawling [5][6][7]. In addition, the mentioned movements are considered as common evacuation behaviour on ships.
e IMO developed requirements for travel time, outset time, and launching time for newly built and existing ships [8,9]. e essentials of IMO need to be verified for newly built and existing ships. Simulation tools are employed for deck heeling conditions to evacuate people. Maritime evacuation software were used in the early 2000s, such as EVI and EXODUS. EVI focuses on overall evacuation where every agent must have the same velocity and attenuates the density flow [7].
EXODUS software is also used for evacuation simulation analysis, which considers change in distance. Moreover, when agent moves from starting point to the arrival point requires evacuation time [10]. e most common evacuation simulation software for land is Pathfinder, which was developed by the US Fire Technology Society. e use of this approach helps Fire Prevention Engineers to simulate the evacuation of people from buildings. erefore, such evacuation software reduces the velocity and increases the cases of heeling hull. Behaviour such as bending over, walking, and crawling cannot simulate the evacuation behaviour, but in comparison to normal walking, people have collision problem.
DMS-series Marine Engine Room Simulator is a commercial device that is widely used in marine accidents and training [11,12]. However, algorithm scripts are written in Java, which simulates evacuations in the case of accidents, operation, and driving of marine equipment for training. While unity engine is very useful in the whole process [13,14]. e major contribution points of this research article are given below.
Distance accumulation lattice gas automata model is designed. Evacuation simulation if DMS-2013 is utilized for marine engine room. Different evacuation behaviours are studied with various deck heeling scenarios. IMO-based several evacuation methods are considered in experiments. Results are formulated from DMS-2013 marine engine room simulator.
Finally, a few scenes before the MV Sewol capsized are demonstrated, which is compared with other evacuation software. Hull heeling increases evacuation time, recommendations to save lives, and properties.

Update Rule for the Distance Accumulation Lattice Gas
Automata Model. Lattice gas automata are the concrete form of cellular automata, where in between grid points there is one particle which moves alone with grid lines at the same time. In 1973, Hardy et al. proposed HPP model for lattice gas automata [15]. en, in 1987, Reynolds utilized Boids model, which gives optimal results of evacuation behaviour [16]. Nowadays, different novel versions of lattice gas automata model are extensively used for evacuation. e distance acquisition lattice gas automata model is a similar kind of approach but quite changed with the lattice gas automata model. However, the updated rule allows cells to move in proper mobility pattern. erefore, 3D grids with unit length d with nine-dimensional array (x, y, z, h, b, v, s, n, c) represent the state of node that can be occupied on deck, where x, y, z represent the coordinates of node and b represents the occupied state at the coordinate position: b � 0. For no cell, b � 1, which specifies cell in the standing state, b � 2 indicates an individual who is bent over, walking, or crawling and b � 3 designates an individual who is bent over, walking, or crawling, also representing trunk. However, v is the instantaneous rate of the cell at the coordinate s having distance of the cell. Moving nodes have initial value of zero, and n is the number of an individual (cells for b � 2 and b � 3, where the same n represents the head and trunk of the exact same individual). is mechanism positions an individual's shoulder width to grid unit d, where the length of two grids illustrates front-to-back distance of an individual. Also, bent-over, walking, or crawling need to be calculated easily. However, h is the height coordinate of ceiling with the above node. Moreover, c is evacuation behaviour code (c � 0 means walking, c � 1 means bent-over walking, and c � 3 means crawling). erefore, the approach uses to be convenient which use different colours to represent various evacuation behaviours.
In addition, different from the previous lattice gas automata model, this method attempts to use the discrete time t of distance accumulation as a parameter for simulation. e approach is not having any fixed value where random values within a certain range are used. However, every cell is not having the right to move at each discrete time. But in some intervals, discreate time is set to t 1 , t 2 , t 3 , . . . , t n up to nth. Moreover, the remaining distance after the kth step is given as s � vt 1 + vt 2 + vt 3 , . . . , vt n − dk.
e value of the remaining distance s for all cells is updated accordingly with respect to discrete time. In addition, the values related with distance s need to have cell with b � 1, 2 and unit grid length d are compared. If s < d, then the cell has no right to move in any direction. But let us suppose s ≥ d, where the value of s is reassigned to cell. erefore, in this method, the cell has right to move one step. However, each cell moves one step, where the discrete time t is selected due to the fact that distance accumulated by the fastest cell at each discrete time is not higher than one grid length d. is ensures that V max t ≤ d; hence, the range of the discrete time t should be 0 < t ≤ d/V max .

Movement Rule for the Distance Accumulation Lattice Gas
Automata Model. Similarly, the general lattice gas automata approach has individuals in terms of distance which exists at grid points and move freely up, down, left, or right. Particularly, range in field is considered standard von Neumann neighbourhood that follows the nonbackwards movement rule [17]. Movement rules are examined to be the main difference between bent-over walking and crawling. Figure 1 shows the x-coordinate of centre position, which is symbolled as x 0 and cell on the left side satisfies X ≤ X 0 . As the cell cannot move backwards in −x-direction, therefore upward and right routes are considered the main gateway near exit. e left direction is reviewed bypass having the probability of cell taking to left is given by P −x � P. e probability of cell moving upwards is given as P y � (1 − P) × |y − y 0 |/(|x − x 0 | + |y − y 0 |), where probability of cell moving rightwards is specified as Individuals with b � 1 and 2 have eight possible conditions, as discussed in Figure 1. Movement rules applicable to these conditions are as below: (1) e above movement rules are used to update the position of each cell with b � 1 and 2 in random order. When the head cell of b � 2 moves, the coordinates with the same label n and b � 3 become unmanned. Coordinates of head cell with b � 2 are assigned the value of b � 3.
e new position is formulated for trunk cell to move properly. However, position of head cell with b � 2 is updated with the passage of time. Movement rule diagram of the distance accumulation lattice gas automata model is represented in Figure 1.

Method and Algorithm
e distance accumulation lattice gas automata model can simulate multiple behaviour modes on horizontal ground. In order to deal with complex sea conditions like DMS-2013, Marine Engine Room simulator combines the model with evacuation direction, deck heeling angles, congestion, and counterflow algorithm which do evacuation behaviour of passengers in cabin.

Selection of the Evacuation Direction.
When the DMS-2013 Marine Engine Room Simulator starts modelling, the xaxis is parallel to the centre axis of hull which is pointed towards the bow and stern. e y-axis is pointed to the port and starboard sides of hull where xy-plane is parallel to the deck. e z-axis is perpendicular to the xy-plane and pointed above the hull. e DMS-2013 Marine Engine Room Simulator uses Dijkstra algorithm to calculate the shortest escape path for passengers. After the escape path is determined, each step of a passenger is mainly along the path. Similar to the principle of random deviation without moving backwards is described in Section 2.2, while moving in perpendicular to main direction is considered bypassing behaviour. As shown in Figure 2, congestion and counterflow are not considered as main direction for passenger a which is positive direction of the y-axis, and positive and negative directions along the x-axis are bypass directions. If the total probability of bypassing both sides is represented by D, where Dε[0, 1], therefore movement rules corresponding to the distance accumulation lattice gas automata model are as follows:

Deck Inclination Algorithm.
To consider the impact of natural environment on evacuation, the effect of deck heeling on the velocity of passengers walking on deck while going upstairs and downstairs should be evaluated. Figure 3 shows all possible evacuation scenarios for tilting hull. e first step is used to determine whether passengers are evacuating on flat ground, going upstairs or downstairs. However, walking state is determined as heeling uphill or downhill. e node of passenger k is represented by the coordinates (x n k , y n k , z n k ) at the discrete time n. en, a passenger is determined to be evacuating on flat ground or going upstairs or downstairs as follows: When a ship capsizes, which will heel to port or starboard in most cases. Moreover, hull heels to starboard and the positive direction of the y-axis points which port opposite to the heeling angle of deck. While y n k − y n−1 k > 0 means the passenger is moving upwards, y n k − y n−1 k � 0 Mobile Information Systems 3 means the passenger is moving transversally, and y n k − y n−1 k < 0 conveys the passenger is moving downwards. Because of the frictional force is too small, where normal person cannot move along a slope with an angle greater than 35° [9,18]. erefore, the reduction of the moving velocity with a passenger is not holding a nearby object which varies within the heeling angle of 0°-35°. While the heeling angle is 0°, and the degree of reduction is 1, the heeling angle must be greater than 35°, and degree of reduction is considered 0. All evacuation scenarios for a tilting hull are shown in Figure 3. e evacuation speed varies greatly because of different factors, such as gender and evacuation behaviour. However, the differences in the trends directly represent change in speed, which is caused by an increase in heeling angle. Also, heeling angle is very small for groups of different genders or evacuation behaviours [18,19]. us, the velocity reduction factor is introduced to approximate the influence of heeling angle on the evacuation speed. Research has shown the reduction in moving velocity with an increase in heeling angle. erefore, the angle is usually linear where the minimum change is 5° [19,20]. Heeling angle range is 5°-25°, while a value is set for the velocity reduction on every 5. However, linear interpolation is used to determine the change in velocity reduction within a heeling angle interval of 5°. e default values for the velocity reduction are shown in Figure 4. is research study discusses the experimental data utilized from Cho [18]. e user manually have input or modified the reduction values for simulation with their own experimental data. Relationship between the reduction of moving velocity and heeling angle of deck is shown in Figure 4.

Bypass and Congestion Algorithm.
Congestion occurs not only in reality but also in software simulation [21]. However, in event-based environment congestion often uses to bypass the main path. erefore, for high flow density ahead, use bypass algorithm for evacuation of peoples [7]. e distance accumulation lattice gas automata model considers direction, which is perpendicular to main flow. However, bypassing probability uses to occupy number of nodes on path which directly establishes relationship for solving the problem. Figure 2 represents passenger a, which is away from the exit and selects exit k at the time t for evacuation. erefore, path needs to be calculated by Dijkstra's algorithm, N k (t) which is having occupancy rate of nodes on path at the time t. However, by using N k (t) instead of bypassing probability D, this update rules (d)-(h). e distance accumulation lattice gas automata model stays the same. In addition, rules (a)-(c) are updated and given below to ensure the probability of moving along with the main direction.   Mobile Information Systems N k (t) defaults the occupancy rate of all nodes on specified path, which can be manually set to the node occupancy rate on predetermined path within n steps of the passenger.

Effects of Congestion Evacuation on Velocity.
e moving velocity of passengers decreases as the surrounding flow density increases [22]. When the flow density is 6.0 persons/m 2 or more, then moving velocity becomes 0.1 m/s, where 1/12 is calculated the normal speed [23], while, in this study, raycasting is used to reflect the effect of the flow density on the moving velocity related to passengers [24]. Although the area occupancy method is specific ray-casting algorithm from simulated environment, area occupancy can be obtained from passenger using the mentioned algorithm. Table 1 presents the default relationship in between occupancy of surrounding nodes and velocity reduction factor. In addition, data processing can be obtained from experiment of Togawa [23] with binomial. A user manually input or modify the reduction values for a simulation using their own experimental data.

Effects of Counterflow on Evacuation Direction.
Counterflow is very commonly utilized in reality, while evacuees try to choose different exits and rescue workers with different objectives to overcome counterflow [25,26]. Road traffic must have pedestrian pathways to establish appropriate rules for avoiding counterflow and congestion. However, most commonly move in one direction having one side. erefore, similar rules are adopted in the evacuation simulation to solve counterflow. e system by default walk on the right side in case of counterflow. e user may manually choose to walk on the left side.
In general, if a node is occupied with the travelling opposite to other node, counterflow behaviour is suggested [27]. In Figure 1, the nodes along with main direction are occupied in (d), (f ), (g), and (h), which may lead to counterflow. When individual i is evacuated, the nodes along with main direction are occupied by individual j. i ⇀ is the unit vector along with main travelling direction of individual i, and j ⇀ is the unit vector along main travelling direction of individual j. erefore, the counterflow is determined as follows: Two persons' counterflow case is discussed with the rule of walking on the right side, which applies the main travelling direction. e positive direction of y-axis must apply movement rules in Figure 1. However, the rules are transformed as follows: In two persons' counterflow where the method should change the direction, the passengers behind them will behave the same which will find a person avoiding the counterflow.

Experiment and Simulation e evacuation function in DMS-2013 Marine Engine Room
Simulator is used to implement about 11 scenarios in MSC/ CIRC 1238 which are proposed by IMO. Related experiments are carried out for some scenarios, where simulation results must be compared other models. However, several scenarios of 2014 try to capsize the MV Sewol which is also simulated. All simulation parameters are set to have default values for the different simulation software.

IMO Test.
e 11 scenarios are directly taken from Annex 3 of MSC/CIRC 1238, which are presented in Table 2.
Test 6 involves increasing evacuation around the corners. However, to simulate walking evacuation experiments are performed for six additional scenarios which include combination of bent-over, walking, and crawling, Mobile Information Systems experimental results are compared to verify the reliability of reflected evacuation with walking, bent-over walking, and crawling. Previously, about 20 undergraduates of nearly the age 18-29 years were tested for the velocities where they are evacuated by walking, bent-over, and crawling. In addition, Table 3 presents the data of about 20 experiments which was imported into DMS-2013 Marine Engine Room Simulator where three different scenarios are simulated. erefore, in each experimentation, about 20 subjects are properly arranged in 5 × 4 array. Figure 5 shows simulation scenarios and distribution of passengers. Moreover, for simulation (a), the 20 subjects walked normally; in experimentation (b), 20 subjects walked bent-over. Although in simulation (c), about 20 subjects crawled. Also, at the end simulation, (a) must need to require Test 6.
e simulation results are shown in Figure 6. Guardrails are used to build the Test 6 scenarios for experiment. Also, infrared timer is utilized at the exit. However, the total evacuation time must be recorded for at least three experiments. Table 4 compares the simulation and experimental results.
Test 8 scenarios are used to verify the counterflow condition. While two rooms with dimensions of 10 × 10 are connected by 10 × 2 channels. However, one hundred occupants in cabin 1 need to move cabin 2, and counterflow from cabin 2 to cabin 1 is directly set to 0, 10, 50, and 100 people.
e test results are shown in Figure 7. Table 5 compares the simulation results of EVIs, Cho, and DMS. Figure 8 represents the scenarios which need to be required for Test 10; 23 people are assigned to at least 12 cabins, while the passengers in Cabins 1-4 and 7-10 need to escape from the main exit. erefore, in other cabins there should be required to escape from secondary exit. In Figure 9, simulation results are properly discussed.

MV Sewol Evacuation Simulation.
In turn, passenger safety is directly related to the size of the MV Sewol heel tipping port. On time 8 : 48 at April 16, 2014, accident occurred where ship was carrying about 476 people. However, 17 on the Nav Deck (5th), 336 on A Deck (4th), and 123 on B Deck (3rd). In addition, at 8 : 50, the heeling angle was 30°at sharp 10 : 09 the ship sunk at heeling angle of 74°. Clearly, the passengers need to perceive heeling angle of 30°and the accident happened during the daytime. Passengers must need to board the rescue boat 123 Vessel through the middle zone near with port side of Nav Deck, A Deck, and B Deck. However, in simulation environment, passengers escaped by walking along the accommodation ladder in the middle of   Maintaining set walking speed in corridor ✓ 7 2 Maintaining set walking speed up staircase ✓ 7 3 Maintaining set walking speed down staircase ✓ 7 4 Exit flow rate ✓ 7 5 Response time Rounding corners ✓ 7 7 Assignment of population demographics parameters ✓ 7 8 Counterflow: two rooms connected via a corridor 7 ✓ 9 Exit flow: crowd dissipation from a large public room 7 ✓ 10 Exit route allocation 7 ✓ 11 Staircase 7 ✓  zone. e port side of each deck before boarding rescue boat needs to be in use. In different situations, passengers need to be evacuated, and there is no specific ratio that needs to be utilized while evacuation. Except of all possible scenarios deck heeling is considered for walking, bent-over walking, and crawling. Also, mixed mode is presented as 1 : 1:1. However, a large proportion of passengers on MV Sewol are having 30 years of age. Different velocities are used for walking, bent-over walking, and crawling which must have rules set for male and female passengers. erefore,    passengers. Figure 10 shows the distribution of passengers evacuated in mixed behaviour mode. Other scenarios had the same coordinates but different evacuation behaviours. Figure 11 shows the evacuation conditions for passengers on Deck A after 60 s of walking at heeling angles of 0°a nd 30°. More passengers reached the exit when the deck was not heeling. Figure 12 compares the scenarios after 60 s of walking, bent-over walking, crawling, and mixed evacuation at a heeling angle of 30°. Together with Figure 11(b), the results showed that walking had the highest evacuation efficiency, while crawling had the lowest efficiency [28][29][30][31]. Table 6 compares the evacuation simulation results for MV Sewol. e number of evacuated passengers, velocity range, and location of gathering points are the same. But the simulation models are completely different. e distance accumulation lattice gas automata model is used to increase the evacuation condition of passengers at constant deck heeling angle. While the simulation was performed with a velocity-based fine network model, dynamic changes in heeling angle need to be considered. Table 5 presents EVI software and Sol which have simulation duration with respect to counterflow density, while counterflow has direct impact on evacuation where passengers from cabin 1 to cabin 2. e consistency has been observed in results especially for test 6 in MSC/CIRC 1238. Figure 11 and Table 6 describe that whether the evacuation is increased with deck heeling angle. e total evacuation duration must be greater in comparison with deck heeling angle of 30°. ese analysis results are observed in simulation from both distance accumulation lattice gas automata model and velocity-based fine network model. e deck heeling angle directly affects the evacuation velocity of the passengers, which determine the evacuation efficiency and total evacuation duration. Figures 11(b), 12 and Table 6 show walking has highest evacuation efficiency while crawling is considered lowest efficiency while deck heeling angle is hardly 30°. erefore, bent-over walking and crawling have lower velocities than walking and larger collision volume, which cause congestion. Although bent-over walking and crawling have the same collision volume and lower velocity, due to that lowest evacuation efficiency is estimated.

Discussion.
e simulation results give insight overview on distance accumulation lattice gas automata model accurately. Various common marine evacuation scenarios with      Figure 11: Scenarios on A Deck after 60 s of walking at heeling angles of (a) 0°and (b) 30°. 0.5 s, but the velocity of an agent cannot be distinguished. Although velocity-based network model discriminates agent speed which refines grid, it cannot increase evacuation modes other than walking, bent-over walking, and crawling. erefore, distance accumulation lattice gas automata model has advantages to increase different evacuation behaviours of passengers. e distance accumulation lattice gas automata model uses von Neumann neighbourhood, where agents mostly move up, down, left, and right. Less moving directions are compared with Moore neighbourhood due to that move upper left, upper right, lower left, and lower right. Cabin passages are either vertical or parallel to the ship's central axis. Due to this approach, von Neumann is considered more suitable. However, drawback can be modified to change the square grid into equilateral triangle. In addition, the evacuation can be increased by using the function of DMS-2013 Marine Engine Room which simulate and fix deck heeling angle. Also, the effects of panic factor and social repulsive forces help in evacuation behaviour.

Conclusion
Distance accumulation lattice gas automata model tries to increase 11 tests which deals with IMO. Evacuation scenarios are based on capsizing of MV Sewol where reliability model needs to be verified. However, simulation results give us the following awareness: (1) Increase in counterflow extend evacuation time, signs which must be provided in cabins to indicate evacuation direction. (2) Reduction of evacuation velocity directly affects evacuation efficiency. erefore, instruction should be issued when ship capsizes. (3) Walking has higher velocity than bent-over walking, crawling, and mixed behaviour. us, the frictional force of ground material should be increased. However, guardrails should be provided on both sides of important evacuation passages to increase the possibility of walking which reduce the evacuation time as well.  Simulating marine evacuation behaviour is a complex process. e scenarios considered in this study are relatively complex. In future work, dynamic changes in evacuation environment and effects of psychological factors on evacuation need to be considered.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.