Next Article in Journal
Acoustic Roughness Measurement of Railhead Surface Using an Optimal Sensor Batch Algorithm
Next Article in Special Issue
The Transient Flow behind an Instantaneously Started Circular Cylinder with Two Symmetrical Strips
Previous Article in Journal
Application of Raman Spectroscopy for Characterizing Synthetic Non-Metallic Inclusions Consisting of Calcium Sulphide and Oxides
Previous Article in Special Issue
Numerical Investigations on the Blade Tip Clearance Excitation Forces in an Unshrouded Turbine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CFD-Based Metamodeling of the Propagation Distribution of Styrene Spilled from a Ship

1
School of Mechanical Engineering, Chung-Ang University, Seoul 06974, Korea
2
LG Electronics, Changwon 51554, Korea
3
Korea Research Institute of Ships & Ocean Engineering (KRISO), Daejeon 34103, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(6), 2109; https://doi.org/10.3390/app10062109
Submission received: 13 January 2020 / Revised: 9 March 2020 / Accepted: 15 March 2020 / Published: 20 March 2020
(This article belongs to the Special Issue Computational Fluid Mechanics and Heat Transfer)

Abstract

:
The present study aimed to numerically establish a new metamodel for predicting the propagation distribution of styrene, which is one of the hazardous and noxious substances (HNSs) spilled from ships. Three-dimensional computational fluid dynamics (CFD) simulations were conducted for 80 different scenarios to gather large amounts of data on the spatial distribution of the change in concentration over time. We used the commercial code of ANSYS Fluent (V.17.2) to solve the Reynolds-averaged Navier–Stokes equations, together with the scalar transport equation. Based on the CFD results, we adopted the well-known kriging model to create a metamodel that estimated the propagation velocity and spatial distributions by considering the effect of the current surface velocity, deep current velocity, surface layer depth, and crack position. The results show that the metamodel accurately predicted the changes in the local distribution of styrene over time. This model was also evaluated using the hidden-point test.

1. Introduction

Hazardous and noxious substances (HNSs), often transported by ships, are harmful to living resources, marine life, and even human health when spilled into the sea [1,2,3]. Because most HNSs, unlike oil, are transparent, the HNSs are difficult to recognize directly after leakage. In general, the HNS can be classified into gases (G), evaporates (E), floaters (F), dissolvers (D), and sinkers (S), as they float, evaporate, or dissolve [4]. For instance, the Ievoli Sun, a chemical tanker, sank in the English Channel and released approximately 4000 tons of styrene [5]. In this accident, high-intensity currents caused significant dilution, spreading styrene in the seawater, and styrene was detected in the gill tissues of crabs in the vicinity of the wreck. More than 2000 HNSs are being transported by ships currently; thus, a robust and quick response strategy for spill accidents should be developed.
According to the literature, there are many reports on oil spill problems. Can et al. [6] performed a numerical analysis of the leakage from an oil container ship in the Istanbul Strait and provided information on the oil propagation. Elhakeem et al. [7] developed a combined model using the hydrodynamic and oil-spill models. They also conducted a numerical simulation of the oil-spill propagation in the case of real accidents in the Arabian Gulf region, showing that the simulated results agreed well with the observed spill trajectories. As mentioned above, the oil spill was not difficult to visualize in the experiment because the oil was not transparent. However, because HNSs are optically transparent and very dangerous it is very difficult to conduct experiments in a real marine environment. Hence, many experimental works have been conducted in restricted areas and conditions. For example, Fuhrer et al. [8] released a floating cell containing a small amount of styrene into the sea and monitored the residual concentration over time. They reported that only 50% of the initial amount of spilled styrene remained just after one hour. Furthermore, experimental data have been reported in the Standard European Behavior Classification (SEBC) code. However, these data have limitations because the experiments were mostly conducted for pure products under atmospheric conditions in a laboratory. It turns out that these conditions significantly differ from those in the case of incidents at sea. Consequently, actual propagation behavior is different from the experimental results reported previously.
Another issue in carrying out the experiments is about the safety caused by the lethal toxicity of the HNSs. When the real experiment is conducted in the sea, the marine ecosystem and living resources will be at risk and eventually destroyed. Thus, doing a real experiment is strictly limited. Computational works are a promising alternative for the detailed investigation of HNS propagation behavior under various conditions. Even if the computational fluid dynamics (CFD) approach yields a lot of information on HNS propagation for various scenarios, it requires a huge amount of time to obtain the final solutions. Practically, we should make specific responses as quickly as possible just after HNS spill accidents. In this regard, CFD simulation cannot be used practically; however, a direct prediction tool, such as a metamodel, can be a powerful way to suggest how to respond to the HNS accident quickly.
Indeed, the metamodel is a surrogate model, which includes the mathematical relationship between inputs and outputs [9,10,11]. To develop the metamodel, we needed to gather a lot of information on the related physics. Hence, the present study extensively conducted CFD simulations for different cases and it suggested important parameters that can be used to mathematically connect inputs and outputs. Previously, a regression model was reported and tested with the CFD simulation results [12,13]. The current study aimed to construct the metamodel with the simulation results for 80 scenarios. Styrene was used as the HNS material and the propagation behavior of styrene in the sea was predicted in detail over a 360° degree range with 10° degree interval. This study also considered the representative current-velocity profile expressed in terms of surface current velocity (vs), deep current velocity (vd), and depth of the surface layer (ds). Moreover, the hidden point tests were conducted to obtain more accurate solutions.

2. Mathematical Representations

2.1. Kriging Model

The metamodel can be efficient and practical for the quick prediction of HNS propagation without further CFD simulations, contributing to a rapid response after an accident [14]. The present study simulated a total of 80 cases and used the well-known kriging model to develop the meta-model [14,15]. According to the kriging model, we first set the unknown function y(x) consisting of a regression model and the stochastic process:
y ( x ) = f ( x ) + Z ( x )
where f(x) represents the approximation function, and Z(x) represents the localized deviation interpolating the n sampled values by considering a Gaussian distribution. A correlation between Z(xi) and Z(xj) is given by a covariance matrix of the samples as follows:
C o v [ Z ( x i ) , Z ( x j ) ] = σ 2 R ( [ R ( x i , x j ) ] ) ,
where R is the correlation matrix, which is an (n × n) diagonal-symmetric matrix, and σ2 is the process variance. R(xi,xj) represents the Gaussian correlation function, which is expressed as:
R ( x k i , x k j ) = exp [ i , j = 1 n θ k | x k i x k j | 2 ] .
The predictor can be described at an arbitrary point x using:
y = β + r T ( x ) R 1 ( y f β )
where y is the CFD simulation results (a column vector of length n), and f is an n × 1 unit column vector. rT(x) represents a column vector {x1, ..., xn} of arbitrary points x. This column vector is expressed as follows:
r T ( x ) = [ R ( x k , x k 1 ) , R ( x k , x k 2 ) , , R ( x k , x k n ) ] T
In Equation (4), β is mathematically given by ( f T R 1 f ) 1 f T R 1 y . The correlation parameter θk in Equation (3) is also determined using maximum likelihood estimation as follows:
max Φ θ k > 0 ( θ k ) = [ n ln ( σ 2 ) + ln | R | ] / 2 ,
where σ 2 is estimated using:
σ 2 = [ ( y f β ) T R 1 ( y f β ) ] / n .

2.2. CFD Model

We needed a lot of data to develop the metamodel of the styrene propagation in the sea. Hence, the present study extensively conducted three-dimensional CFD simulations to obtain the predictions of the styrene concentration distribution by solving the Reynolds-averaged Navier–Stokes equations and the scalar transport equations at the same time. We used the commercial code of ANSYS Fluent (V.17.2) for simulations. Since the characteristic length is in the order of several hundred meters, the flow can be regarded as a turbulent flow around the ship. The present study adopted the standard k-ε model for simplicity. Numerically, we solved the mass and momentum conservation equations, which are expressed as follows:
ρ t + ( ρ u ̄ i ) x i = 0 ,
t ( ρ u ̄ i ) + x i ( ρ u ̄ i u ̄ j ) = p x i + x i [ μ ( u ̄ i x j + u ̄ j x i 2 3 δ i j u ̄ l x l ) ] + x j ( ρ u i u j ) ,
where ρ denotes density, ui denotes the velocity vector, and μ denotes viscosity. The present simulation used the standard k-ε model, and the wall function to predict turbulent stresses at the wall in the simulation [16]. The standard k-ε model can be expressed as:
t ( ρ k ) + x i ( ρ u ̄ i k ) = x j [ ( μ + μ t σ k ) k x j ] + G k + G b ρ ε ,
t ( ρ ε ) + x i ( ρ u ̄ i ε ) = x j [ ( μ + μ t σ ε ) ε x j ] + C 1 ε k ( G k + G b ) C 2 ρ ε 2 k ,
μ t = C μ ρ k 2 ε ,
where k and ε represent the turbulent kinetic energy and the dissipation rate, respectively. Furthermore, μt denotes the turbulent viscosity and Gk denotes the production of turbulent kinetic energy due to the mean shear rates. Here, Gb indicates the turbulent kinetic energy production due to buoyancy. The model constants C1, C2, and Cμ were set to C1 = 1.22, C2 = 1.92, and Cμ = 0.09. The turbulent Prandtl numbers were set to Prk = 1.0 and Prε= 1.3. The species transport equation was adopted to describe the convection–diffusion for the kth species as follows:
t ( ρ Y k ) + ( ρ u ̄ i Y k ) x i = J j k x j ,
where Jjk and Yk are the diffusion flux and the mass fraction of species k, respectively.

2.3. Numerical Details

Among various HNSs, the present study used a styrene material for the simulation based on knowledge of detailed scenarios with various current–velocity conditions. Figure 1 shows the schematic of the model that represents the current-velocity profile. It is known that the current-velocity profile consists of a tide-induced current, a density-driven current, and a wind-driven current [17,18]. Because of the complexity of atmospheric conditions, we modeled the current-velocity profile in terms of the surface current velocity (vs), deep current velocity (vd), and depth of surface layer (ds). The present study considered four parameters, as follows: i) the surface velocity (vs), ii) the ratio of the deep current velocity to surface current velocity (vd/vs), iii) the depth of the surface layer (ds), and iv) the crack positions as the side (collision) and bottom (grounding). The surface velocity was set in the range from 0.1 to 1 m/s, and the deep current velocity was taken to be 30% to 100% of the surface current velocity (vs). Furthermore, the depth of the surface layer was in the range of 10 to 100 m. The simulation conditions were finally set for 80 cases. The present simulation ignored the interfacial waves between the seawater and air for simplicity because it would have required tremendous computational resources and complicated data about atmospheric conditions. Hence, the slip condition at sea level was taken as the boundary condition. Furthermore, the pressure out condition was used for all surfaces except for the bottom surface, upon which the wall boundary condition was applied. The density, viscosity, and diffusivity properties of styrene were set to 906 (kg/m3), 0.76 × 10−3 (kg/(m∙s)) and 2.0 × 10−9 (m2/s), respectively. It was also important to establish the accident scenario. This simulation was conducted under the scenarios considering two cracks from the side (collision) and bottom (grounding) [19,20]. Based on information provided by KRISO (Korea Research Institute of Ships & Ocean Engineering), a mass flow rate from a crack was set as 500 kg/s. With a time step of 0.05 s, the transient simulation was performed for a total simulation time of 1000 s. The corresponding averaged Courant number was estimated to be 0.1231, and the maximum Courant number of 0.9232 was found near the crack position. The estimated y+ values were in the range of 100 to 220, showing an acceptable range for predicting turbulent stresses appropriately.
As shown in Figure 2, the computational domain had a size of 1000 m × 1000 m × 100 m. The detailed configuration of the ship was provided by KRISO. The accident ship had a length of 160 m, a width of 30 m, and a draft of 14.5 m. The crack with its area of 0.25 m2 was positioned at 80 m from the front of the ship. Furthermore, another scenario considered a side crack that could form owing to collision-induced damage. The side-crack was located at 10 m below sea level, and the bottom crack due to grounding was located at the center of the bottom surface of the ship. We used the tetrahedral grid system and built much finer grids around the ship, as shown in Figure 2b. For the grid convergence test, we compared the velocity profile behind the ship (x = 290, z = 0, 450 m < y < 550 m) in the case of vs = 1, vd/vs = 0.5, and ds = 30 m. As shown Figure 3a, the maximum deviation between grid numbers of 1,770,958 and 2,482,851 was estimated to be 0.7% at y = 500 m. Therefore, the present study adopted 1,770,958 elements and 322,925 nodes after the grid convergence tests. We set the convergence criterion for the continuity as 10−4, for the energy as 10−6, and for others as 10−3.

3. Results and Discussion

3.1. Propagation Characteristics

Figure 4a,b presents the predicted propagation behavior of styrene with different current velocities at a flow time 1000 s after leaking began. As expected, the surface current velocity (vs) substantially affects the diffusive characteristics of styrene propagation. As the surface current velocity increased, styrene propagated much faster owing to the advection effect, whereas the styrene propagation was barely affected by the crack positions. The latter was because styrene, as a floating material, rapidly moved upward after being spilled from the crack, indicating the propagation behaviors at the side and bottom cracks were almost the same. Unlike the bottom crack, the styrene distribution in the side crack case showed a slightly asymmetric shape. After considering three factors, namely the surface current velocity, deep current velocity, and surface layer depth, it was found that the surface current velocity significantly affected the propagation behavior after a styrene spill.

3.2. Estimation of Styrene Distribution

As stated above, saving time costs and getting accurate solutions are the ultimate objectives of metamodeling. Hence, we numerically simulated 80 cases and stored the acquired data for different variables, such as local velocity, kinetic energy, pressure, concentration, and so on. A metamodel connects inputs (scenario) and outputs (prediction) mathematically. Thus, we need to find out new parameters that are useful for responding to an accident. The present study suggested two parameters, as illustrated in Figure 5, which were the change in propagation distribution with time and the propagation velocity of the styrene interface.
First, the two parameters were estimated from the CFD results, and then those parameters were implemented in the kriging model. In Figure 5, a criterion regarding styrene concentration was taken as 100 ppm, corresponding to the exposure limit of styrene specified by the Occupational Safety and Health Administration [21]. Based on this criterion, the propagation velocity of styrene interface could be calculated as follows: first, the center point was selected as a middle point of the styrene interface at t = 300 s after the spill, and the interface area was then divided into intervals of 10° based on the center point, as shown in Figure 5. At a certain time, a two-dimensional interface contour line and the lines formed every 10° intersected each other at some points; the propagation velocity of the styrene interface and the propagation distribution were then estimated at these points. Finally, all data successively stored with time were used to make the metamodel.
Figure 6 shows the predicted styrene interface location at t = 1000 s according to the surface current velocity, which was predicted using the present metamodel in the bottom crack case with vd/vs = 0.5 and ds = 30 m. The blue dot represents the styrene interface. It shows that at a low surface current velocity, the diffusion effect became dominant, and the styrene interface had an almost circular shape. As the surface current velocity increased, the shape of the styrene interface was closer to an ellipse due to the large advection effect. The goodness of fit was evaluated using R2 and RMSE (root mean square error). As listed in Table 1, the R2 values were greater than 0.9, which means that the metamodel provided a good fit.

3.3. Hidden-Point Tests for the Evaluation

Mostly, comparison with experimental data has been widely adopted for model evaluation in many relevant numerical studies. Unfortunately, there is no available data for direct comparison of styrene under the specified conditions. Instead, we used another method, the so-called hidden point (or layer) method, which has been mainly introduced in machine learning and neural networks. More clearly, the hidden point (or condition) is not included in 80 different CFD cases that are used in making the current metamodel. Hidden points are new scenarios. Based on them, the CFD simulation was first conducted and then the metamodel ran separately, and finally, the two results were compared. We also selected three additional parameters, the first of which was the boundary arrival time tb, defined as a time at which the styrene interface reached the outmost boundary face of the computational domain. The two other parameters were the x and y coordinates of the styrene interface at t = 1000 s after the spill started. Figure 7 shows the contours calculated by the metamodel at the given Ds of 30 m in the case of a bottom crack. The surface represents the response of the metamodel for the inputs, the black dot represents the CFD simulation result from 80 cases taken in the present study, and the red dot indicates the metamodel prediction based on the hidden point.
The comparison result is summarized in Table 2. The CFD simulations were performed for eight additional cases, which were regarded as the hidden points. Those eight cases are listed in Table 2 were not included in the 80 cases used for making the metamodel. For measuring the relative error between the CFD result and the metamodel prediction, the quantitative error of boundary arrival time could be estimated using:
δ 1 = V V V × 100 ( % )
where V’ denotes the metamodel prediction at the hidden point and V represents the CFD result at the same point. For the x- and y-coordinates of the styrene interface, the relative error was obtained using:
δ 2 = i = 1 3 j = 1 36 k = 1 42 ( V V V × 100 ) i j k ( % )
where i denotes the coordinate, and j and k represent the angle and the time, respectively. The estimated error of the styrene interface was less than 9.2%, which shows a relatively good performance. As shown in Table 2, the created metamodel showed relatively good for qualitatively predicting how the spatial distribution changed with time, though it provided a big difference and showed a low accuracy.
The kriging method is a kind of inverse distance weighted (IDW) method. It uses the interpolation by giving larger weights to the adjacent points. Therefore, the accuracy of the metamodel can be increased by extending the database. To improve the accuracy of the metamodel for predicting the boundary arrival time, four cases (3, 4, 6, and 7) were selected and put into the previous database for 80 cases. Then, with the newly updated database with 84 cases, the metamodel could be updated by calculating the coefficients again with the use of Equation (4). For the evaluation, we compared the results of the modified metamodel with the CFD results for the hidden point cases 1, 2, 5, and 6. Table 3 shows the relative errors for the modified metamodel, which shows a better performance at predicting the boundary arrival time. For instance, the 37.4% error for case 6 was reduced to 11.8%. This model will be updated continuously by performing further simulations for additional cases to enhance the model’s performance.

4. Conclusions

The present study developed a new metamodel for predicting the styrene propagation distribution and boundary arrival time. Two important parameters were introduced and the extensive CFD simulations were carried out. A mathematical relationship between inputs and outputs was established using the kriging model. To feed the big data that were necessary to make the metamodel, the present study extensively conducted three-dimensional transient CFD simulations for 80 different cases from the taken scenarios. The following conclusions were drawn.
The present study presented the CFD simulation results and analyzed the influence of current profiles on the propagation characteristics. Due to the floating behavior of styrene, the concentration fields were concentrated near the sea level and the interface of styrene moved much faster in the streamwise direction as the surface velocity increased. This was clearly because the advection effect became dominant and most of the styrene propagated near the sea’s surface.
  • A new metamodel was provided to estimate the propagation distribution and the boundary arrival time. Based on CFD results, the main parameters were calculated and implemented in the metamodel. A comparison was made between CFD results and the metamodel prediction, showing good agreement between them. This verified that the current model could predict the transient characteristics of styrene propagation well. Thus, the use of the metamodel would be a powerful tool for the quick estimation of HNS propagation that can be used to formulate a fast response in the early stages after accidents.
  • This metamodel was evaluated using the hidden point tests. By adding data from eight additional cases, the performance of the metamodel was improved. For instance, a 37.4% error was reduced to 11.8% due to the modification. Thus, the current model will be updated continuously to achieve better accuracy via modification with additional data from further CFD simulations.

Author Contributions

Conceptualization and supervision, S.H.L. and M.L.; investigation, C.H.J. and M.K.K.; methodology, C.H.J. and M.K.K.; supervision, S.H.L.; writing of the original draft, C.H.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was fundded by the Development of Human Resources for Design Specialists in Applications of Functional Materials (Brain Korea (BK) 21) and the Chung-Ang University Graduate Research Scholarship (2016).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IMO. Protocol on preparedness, Response and co-operation to pollution incidents by hazardous and noxious substances. 2000 (OPRC-HNS Protocol); IMO: London, UK, 2000. [Google Scholar]
  2. Neuparth, T.; Moreira, S.; Santos, M.M.; Reis-Henriques, M.A. Hazardous and noxious substances (HNS) in the marine environment: Prioritizing HNS that pose major risk in a european context. Mar. Pollut. Bull. 2011, 62, 21–28. [Google Scholar] [CrossRef] [PubMed]
  3. Harold, P.D.; De Souza, A.S.; Louchart, P.; Russell, D.; Brunt, H. Development of a risk-based prioritisation methodology to in-form public health emergency planning and preparedness in case of accidental spill at sea of hazardous and noxious substances (HNS). Environ. Int. 2014, 72, 157–163. [Google Scholar] [CrossRef] [PubMed]
  4. Bonn Agreement, Counter pollution manual. Bonn Agreemment Counter Pollution Manual. 1991. Available online: https://www.bonnagreement.org/publications.
  5. Law, R.J.; Kelly, C.; Matthiessen, P.; Aldridge, J. The loss of the chemical tanker Ievoli Sun in the English Channel, October 2000. Mar. Pollut. Bull. 2003, 46, 254–257. [Google Scholar] [CrossRef]
  6. Can, S.; Celik, F.; Yilmaz, H.; Bak, O.A. Oil spill simulation: A case study in the strait of Istanbul. Fresenius Environ. Bull. 2007, 16, 1517–1522. [Google Scholar]
  7. Elhakeem, A.A.; Elshorbagy, W.; Chebbi, R. Oil spill simulation and validation in the Arabian (Persian) Gulf with special reference to the UAE coast. Water Air Soil Pollut. 2007, 184, 243–254. [Google Scholar] [CrossRef]
  8. Fuhrer, M.; Peron, O.; Hofer, T.; Morrissette, M.; Le Floch, S. Offshore experiments on styrene spillage in marine waters for risk assessment. Mar. Pollut. Bull. 2012, 64, 1367–1374. [Google Scholar] [CrossRef]
  9. Wen, J.; Yang, H.Z.; Jian, G.P.; Tong, X.; Li, K.; Wang, S.M. Energy and cost optimization of shell and tube heat exchanger with helical baffles using kriging metamodel based on MOGA. Int. J. Heat Mass Tranf. 2016, 98, 29–39. [Google Scholar] [CrossRef]
  10. Rulik, S.; Wroblewski, W.; Fraczek, D. Metamodel-based optimization of the Labyrinth Seal. Arch. Mech. Eng. 2017, 64, 75–91. [Google Scholar] [CrossRef]
  11. Mao, X.F.; Zhang, W.X.; Qian, J.K.; Wu, M.H. CFD based multi-disciplinary optimization design of high performance deep sea seismic vessel. In Proceedings of the ASME 36th International Conference on Ocean, Offshore and Arctic Engineering, Trondheim, Norway, 25–30 June 2017; p. 1. [Google Scholar]
  12. Jeong, C.H.; Ko, M.K.; Lee, M.; Lee, S.H. Numerical simulation of propagation characteristics of hazardous noxious substances spilled from transport ships. Appl. Sci. 2018, 8, 2409. [Google Scholar] [CrossRef] [Green Version]
  13. Ko, M.K.; Jeong, C.H.; Lee, M.; Lee, S.H. Development of a metamodel for predicting near-field propagation of hazardous and noxious substances spilled from a ship. Appl. Sci. 2019, 9, 3838. [Google Scholar] [CrossRef] [Green Version]
  14. Simpson, T.W.; Mistree, F. Kriging models for global approximation in simulation-based multidisciplinary design optimization. AIAA J. 2001, 39, 2233–2241. [Google Scholar] [CrossRef] [Green Version]
  15. Wahba, G. Spline Models for Observational Data (CBMS-NSF Regional Conference Series in Applied Mathematics); Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1990; p. 169. [Google Scholar]
  16. Launder, B.E.; Spalding, D.B. The numerical computation of turbulent flows. Comput. Methods Appl. Mech. Eng. 1974, 3, 269–289. [Google Scholar] [CrossRef]
  17. Sellar, B.G.; Wakelam, G.; Sutherland, D.R.J.; Ingram, D.M.; Venugopal, V. Characterisation of tidal flows at the European marine energy center in the absence of ocean waves. Energies 2018, 11, 176. [Google Scholar] [CrossRef] [Green Version]
  18. Samolyubov, B.I.; Ivanova, I.N. The evolution of velocity profiles and turbulent viscosity in a system of currents with wind-induced and density flows. Mosc. Univ. Phys. Bull. 2014, 69, 426–432. [Google Scholar] [CrossRef]
  19. Ventikos, N.P.; Papanikolaou, A.D.; Louzis, K.; Koimtzoglou, A. Statistical analysis and critical review of navigational accidents in adverse weather conditions. Ocean Eng. 2018, 163, 502–517. [Google Scholar] [CrossRef]
  20. Friis-Hansen, P.; Simonsen, B.C. GRACAT: Software for grounding and collision risk analysis. Mar. Struct. 2002, 15, 383–401. [Google Scholar] [CrossRef]
  21. Occupational Safety and Health Administration, OSHA Occupational Chemical Database. 2018. Available online: https://www.osha.gov/chemicaldata/chemResult.html?recNo=14 (accessed on 18 December 2018).
Figure 1. The schematic of the representative current profile.
Figure 1. The schematic of the representative current profile.
Applsci 10 02109 g001
Figure 2. Computational domain and configuration of the ship (provided by the Korea Research Institute of Ships & Ocean Engineering).
Figure 2. Computational domain and configuration of the ship (provided by the Korea Research Institute of Ships & Ocean Engineering).
Applsci 10 02109 g002
Figure 3. (a) Grid convergence test results for the case of vs = 1, vd/vs = 0.5, and ds = 30 m at x = 290 m and z = 0 m, (b) along the black solid line A–A`.
Figure 3. (a) Grid convergence test results for the case of vs = 1, vd/vs = 0.5, and ds = 30 m at x = 290 m and z = 0 m, (b) along the black solid line A–A`.
Applsci 10 02109 g003
Figure 4. Propagation behavior of the styrene with different surface current velocities (vs), spilled from the (a) side crack and (b) bottom crack (t = 1000 s, vd/vs = 0.5, and ds = 30 m).
Figure 4. Propagation behavior of the styrene with different surface current velocities (vs), spilled from the (a) side crack and (b) bottom crack (t = 1000 s, vd/vs = 0.5, and ds = 30 m).
Applsci 10 02109 g004
Figure 5. Schematic showing the interface of styrene and the interface propagation velocity.
Figure 5. Schematic showing the interface of styrene and the interface propagation velocity.
Applsci 10 02109 g005
Figure 6. Interface of the spilled styrene predicted using the metamodel (black dots) and computational fluid dynamics (CFD) simulation (blue dots) in the case of t = 1000 s, vd/vs = 0.5, ds = 30 m, and Lc = bottom crack. The criterion for the interface was defined as a mass fraction corresponding to the exposure limit (100 ppm).
Figure 6. Interface of the spilled styrene predicted using the metamodel (black dots) and computational fluid dynamics (CFD) simulation (blue dots) in the case of t = 1000 s, vd/vs = 0.5, ds = 30 m, and Lc = bottom crack. The criterion for the interface was defined as a mass fraction corresponding to the exposure limit (100 ppm).
Applsci 10 02109 g006
Figure 7. Contours of the (a) boundary arrival time, as well as the (b) x- and (c) y-coordinates of the styrene interface facing 50° from the middle point at t = 1000 s after spill starts.
Figure 7. Contours of the (a) boundary arrival time, as well as the (b) x- and (c) y-coordinates of the styrene interface facing 50° from the middle point at t = 1000 s after spill starts.
Applsci 10 02109 g007
Table 1. The summary of fit with using R2 and RMSE.
Table 1. The summary of fit with using R2 and RMSE.
Vs (m/s)Vd/VsDs (m)LcR2RMSE
0.10.530Bottom0.97732.17
0.20.530Bottom0.93356.43
0.50.530Bottom0.95756.91
1.00.530Bottom0.949412.83
Table 2. Eight additional cases for the hidden-point test and the corresponding errors of the metamodel for the boundary arrival time (tb).
Table 2. Eight additional cases for the hidden-point test and the corresponding errors of the metamodel for the boundary arrival time (tb).
Case No.Vs (m/s)Vd/VsDs (m)Lcδ1 (%)δ2 (%)
Case 10.300.530Side12.39.2
Case 20.300.530Bottom12.17.3
Case 30.400.530Side11.06.4
Case 40.400.530Bottom12.15.1
Case 50.700.530Side19.64.8
Case 60.700.530Bottom37.45.6
Case 70.800.530Side37.83.7
Case 80.800.530Bottom46.64.7
Table 3. Error analysis for the modified metamodel.
Table 3. Error analysis for the modified metamodel.
Case No.δ1 (%)
Case 15.0
Case 24.0
Case 59.7
Case 611.8

Share and Cite

MDPI and ACS Style

Jeong, C.H.; Ko, M.K.; Lee, M.; Lee, S.H. CFD-Based Metamodeling of the Propagation Distribution of Styrene Spilled from a Ship. Appl. Sci. 2020, 10, 2109. https://doi.org/10.3390/app10062109

AMA Style

Jeong CH, Ko MK, Lee M, Lee SH. CFD-Based Metamodeling of the Propagation Distribution of Styrene Spilled from a Ship. Applied Sciences. 2020; 10(6):2109. https://doi.org/10.3390/app10062109

Chicago/Turabian Style

Jeong, Chan Ho, Min Kyu Ko, Moonjin Lee, and Seong Hyuk Lee. 2020. "CFD-Based Metamodeling of the Propagation Distribution of Styrene Spilled from a Ship" Applied Sciences 10, no. 6: 2109. https://doi.org/10.3390/app10062109

APA Style

Jeong, C. H., Ko, M. K., Lee, M., & Lee, S. H. (2020). CFD-Based Metamodeling of the Propagation Distribution of Styrene Spilled from a Ship. Applied Sciences, 10(6), 2109. https://doi.org/10.3390/app10062109

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