Next Article in Journal
The 2022 Seismic Sequence in the Northern Adriatic Sea and Its Long-Term Simulation
Next Article in Special Issue
Performance of Base-Isolated RC School Building under Blast Loading
Previous Article in Journal
Study on Dynamic Impact Mechanical Properties of UHPC with High-Content and Directional Reinforced Steel Fiber
Previous Article in Special Issue
Research on the Mechanism and Safe Thickness of Karst Tunnel-Induced Water Inrush under the Coupling Action of Blasting Load and Water Pressure
 
 
Article
Peer-Review Record

Behavior of Barrier Wall under Hydrogen Storage Tank Explosion with Simulation and TNT Equivalent Weight Method

Appl. Sci. 2023, 13(6), 3744; https://doi.org/10.3390/app13063744
by Seungwon Kim, Taejin Jang, Topendra Oli * and Cheolwoo Park *
Reviewer 1:
Reviewer 2: Anonymous
Reviewer 3: Anonymous
Reviewer 4:
Appl. Sci. 2023, 13(6), 3744; https://doi.org/10.3390/app13063744
Submission received: 3 February 2023 / Revised: 11 March 2023 / Accepted: 13 March 2023 / Published: 15 March 2023
(This article belongs to the Special Issue Blast Loading and Blast Effect on Structures)

Round 1

Reviewer 1 Report

Major revision has to be performed in order to make the paper easier to read and to clarify certain points.

Main points:

1)     Explain why for this study a 75% ratio of hydrogen with air was taken.

2)     Why do you take explosive energy of TNT = 4.6 MJ/kg instead of generally accepted 4.2 MJ/kg?

3)     Could you please give and explanation as to why the authors of [27] consider a value of 0.01 for the explosion efficiency (Eq. 1). Was it based on experimental data or numerical analysis?

4)     For the LS-DYNA material models it would be better to provide corresponding formulae for the MAT_JHONSON_COOK and other models. For a reader who is not familiar with these models, the Table 4 and 5 are not informative. Same for the HIGH_EXPLOSIVE_BURN – need to present the analytical version of the model.

5)     What is the justification for the choice of the mesh sizes (100 mm in air and 20 mm inside the structure). Have you performed grid sensitivity studies? If Yes, please provide the results.

Minor remarks:

a)      The phrase (lines 11-12) is not clear. Please revise.

b)     Line 43-44: the lower and upper limits for combustion (not laminar flame).

c)      Line 194: hydrogem..

d)     Line 211: In figure 2 i_s is….

e)     Line 220: “The model was modelled…” Please revise.

Line 315: …table 7

Author Response

Response to Reviewer 1 Comments

Thank you for review

 

Point 1

Explain why for this study a 75% ratio of hydrogen with air was taken

Response 1

Thank you for your comment.

In present of 4 % to 75 % ratio of hydrogen with air, explode can occur. If we take less than 75 %, the equivalent weight of TNT will be less. To create worse case scenario 75% was selected.

 

Point 2

Why do you take explosive energy of TNT = 4.6 MJ/kg instead of generally accepted 4.2 MJ/kg?

Response 2

Thank you for your comment.

Explosive energy was taken from this article. This article calculate the explosion energy of TNT.

Kwon, S.; Park, J. A Review of TNT Equivalent Method for Evaluating Explosion Energy Due to Gas Explosion.

 

Point 3

Could you please give and explanation as to why the authors of [27] consider a value of 0.01 for the explosion efficiency (Eq. 1). Was it based on experimental data or numerical analysis?

Response 3

Thank you for your comments.

It has been updated in manuscript.

 

Point 4

For the LS-DYNA material models it would be better to provide corresponding formulae for the MAT_JHONSON_COOK and other models. For a reader who is not familiar with these models, the Table 4 and 5 are not informative. Same for the HIGH_EXPLOSIVE_BURN – need to present the analytical version of the model.

Response 4

Thank you for your comment.

References has been added.

 

Point 5

What is the justification for the choice of the mesh sizes (100 mm in air and 20 mm inside the structure). Have you performed grid sensitivity studies? If Yes, please provide the results.

Response 5

Thank you for your comment.

Grid sensitivity analysis had not been studied but the mesh size was taken from the verification of previous studies. 

 

Point 6

Minor remarks:

a) The phrase (lines 11-12) is not clear. Please revise.

Response 6 a

Thank you for your comment.

It has been corrected.

b) Line 43-44: the lower and upper limits for combustion (not laminar flame).

Response 6 b

Thank you for your comment.

It has been corrected.

c) Line 194: hydrogem..

Response 6 c

Thank you for your comment.

It has been corrected.

d) Line 211: In figure 2 i_s is….

Response 6 d

Thank you for your comment.

It has been corrected.

e) Line 220: “The model was modelled…” Please revise.

Response 6 e

Thank you for your comment.

It has been corrected

f) Line 315: …table 7

Response 6 f

Thank you for your comment.

It has been corrected.

Author Response File: Author Response.pdf

Reviewer 2 Report

This manuscript numerically investigated the dynamic response of H-beam walls under blast loading. A total of 36 numerical models were developed to examine the effects of structural layout, wall thickness, blast distance and TNT mass factors. As the ALE method is very expensive and the numerical models in the manuscript had large dimensions, a very large and costly computational effort was performed by the authors. The manuscript is of high quality and the research can contribute to the study of the impact dynamics of RC structures, within the scope of Applied Sciences, and the reviewer recommends publication of the manuscript after responding the following comments.

-        L76, Nassr et al., 2012, citation error.

-        In the introductory section, a literature review on the numerical simulation of explosions in RC structures was missing.

-        Eq.(1), L296, “X” should be “×”.

-        In section 2.2, an original reference to the Friedlander equation should be added to the presentation of the overpressure relation for free air explosion waves for the convenience of the reader. e.g., (10.1098/rspa.1946.0046)

-        L217, the names of a,b,c in figure 4 were duplicated.

-        It would be better to present more details about the finite element model, such as the mesh discretization, the cross section of the RC wall for different concrete thickness, etc.

-        L249, The KCC and CSC models were applied to simulate concrete, but an introduction to the models, or references to them, were missing. References to material constitutive models should to be added. For the KCC model: (10.1061/(ASCE)EM.1943-7889.0000952), (10.1016/j.cemconcomp.2022.104911); For the CSC model: (FHWA-HRT-05-062), (10.1016/j.compstruct.2022.116424).

-        L255, the original work of Malvar should be cited. (10.1016/s0734-743x(97)00023-7)

-        How was the fluid-solid coupling in the manuscript achieved? Was it by using the traditional *CONSTRAINED_LAGRANGE_IN_SOLID, or did it use *ALE_STRUCTURED_FSI in the structured ALE approach. please clarify.

-        The reviewer did not find a description of the contact between the reinforcement, H-beam and concrete, please check.

-        It is necessary to add the error (%) between the numerical analysis and the experimental results in Table 6.

-        The ALE method of calculating blast wave propagation requires a very fine mesh, and the underestimation of peak pressure at 1600 mm might be the reason for this. The "KPa" in Table 7 should be "kPa"

-        It is advisable to revise Figure 8 to make the graph easier to understand and the axis ranges need to be modified to fit the data.

-        From Figure 10, it seems that the element erosion method was used, but the reviewer did not find a description in the manuscript, please explain.

-        If the element erosion method was used, then the volume of the failed element can be calculated in LS-prepost, which would provide a quantitative analysis for this study.

-        The effective plastic strain in Figure 10 was essentially an index of damage, with strain hardening between 0 and 1 and strain softening between 1 and 2. Please explain this historical variable in the revision to avoid misinterpretation.

-        Scaled distances for different load cases should be calculated and provided.

-        Is it possible to obtain advice on which structural forms and dimensions are safe at certain blast distances under the most adverse blast loading scenario?

-        The reinforcement arrangement in the RC wall needs to be supplemented.

Author Response

Response to Reviewer 2 Comments

This manuscript numerically investigated the dynamic response of H-beam walls under blast loading. A total of 36 numerical models were developed to examine the effects of structural layout, wall thickness, blast distance and TNT mass factors. As the ALE method is very expensive and the numerical models in the manuscript had large dimensions, a very large and costly computational effort was performed by the authors. The manuscript is of high quality and the research can contribute to the study of the impact dynamics of RC structures, within the scope of Applied Sciences, and the reviewer recommends publication of the manuscript after responding the following comments.

Thank you for review

 

Point 1

L76, Nassr et al., 2012, citation error.

Response 1

Thank you for your comment.

It has been corrected.

 

Point 2

In the introductory section, a literature review on the numerical simulation of explosions in RC structures was missing.

Response 2

Thank you for your comment.

Literature section has been updated.

 

Point 3

Eq.(1), L296, “X” should be “×”.

Response 3

Thank you for your comments.

It has been corrected.

 

Point 4

In section 2.2, an original reference to the Friedlander equation should be added to the presentation of the overpressure relation for free air explosion waves for the convenience of the reader. e.g., (10.1098/rspa.1946.0046)

Response 4

Thank you for your comment.

It has been corrected.

 

Point 5

L217, the names of a,b,c in figure 4 were duplicated.

Response 5

Thank you for your comment.

It has been corrected.

Titile of the article had been corrected.

 

Point 6

It would be better to present more details about the finite element model, such as the mesh discretization, the cross section of the RC wall for different concrete thickness, etc

Response 6

Thank you for your comment.

It has been added and mesh discretization of the simulation had been missing for this study.

 

Point 7

L249, The KCC and CSC models were applied to simulate concrete, but an introduction to the models, or references to them, were missing. References to material constitutive models should to be added. For the KCC model: (10.1061/(ASCE)EM.1943-7889.0000952), (10.1016/j.cemconcomp.2022.104911); For the CSC model: (FHWA-HRT-05-062), (10.1016/j.compstruct.2022.116424).

Response 7

Thank you for your comment.

It has been corrected

 

Point 8

L255, the original work of Malvar should be cited. (10.1016/s0734-743x(97)00023-7)

Response 8

Thank you for your comment.

It has been corrected

 

Point 9

How was the fluid-solid coupling in the manuscript achieved? Was it by using the traditional *CONSTRAINED_LAGRANGE_IN_SOLID, or did it use *ALE_STRUCTURED_FSI in the structured ALE approach. please clarify.

Response 9

Thank you for your comment.

It has been updated

 

Point 10

The reviewer did not find a description of the contact between the reinforcement, H-beam and concrete, please check.

Response 10

Thank you for your comment.

It has been updated

 

Point 11

It is necessary to add the error (%) between the numerical analysis and the experimental results in Table 6.

Response 11

Thank you for your comment.

It has been updated

 

Point 12

The ALE method of calculating blast wave propagation requires a very fine mesh, and the underestimation of peak pressure at 1600 mm might be the reason for this. The "KPa" in Table 7 should be "kPa"

Response 12

Thank you for your comment.

It has been corrected

 

Point 13

It is advisable to revise Figure 8 to make the graph easier to understand and the axis ranges need to be modified to fit the data.

Response 13

Thank you for your advise

 

Point 14

From Figure 10, it seems that the element erosion method was used, but the reviewer did not find a description in the manuscript, please explain.

Response 14

Thank you for your comment.

It has been added

 

Point 15

If the element erosion method was used, then the volume of the failed element can be calculated in LS-prepost, which would provide a quantitative analysis for this study.

Response 15

Thank you for your comment.

For this study, erode volume of elements after the simulation was not studied.

 

Point 16

The effective plastic strain in Figure 10 was essentially an index of damage, with strain hardening between 0 and 1 and strain softening between 1 and 2. Please explain this historical variable in the revision to avoid misinterpretation.

Response 16

Thank you for your comment.

It has been upadated

 

Point 17

Scaled distances for different load cases should be calculated and provided.

Response 17

Thank you for your comment.

It has been added in section 4.1.

 

Point 18

Is it possible to obtain advice on which structural forms and dimensions are safe at certain blast distances under the most adverse blast loading scenario?

Response 18

Thank you for your comment.

From this study, it is confirmed that when the thickness of the protection wall increased to 200 mm the structure gets low damage for both reinforced and composite structure at the distance of 3.6 m. Future study should be continue.

 

Point 19

The reinforcement arrangement in the RC wall needs to be supplemented.

Response 19

Thank you for your comment.

It has been corrected.

Author Response File: Author Response.pdf

Reviewer 3 Report

The manuscript investigates the behavior of barrier walls subject to explosions. The research is consistent with the aims and scopes of the journal.

The manuscript is very well written and I did not find any particular issue that should be addressed.

Possibly, authors could consider a better representation of the displacements in Figure 8 and reduce the font size of colorbars in Figure 10 (if possible). Nevertheless, such minor issues do not represent real drawbacks and could be implemented during the publication phase.

In conclusion, I recommend it to be accepted in the present form.

 

Author Response

Response to Reviewer 3 Comments

The manuscript investigates the behavior of barrier walls subject to explosions. The research is consistent with the aims and scopes of the journal.

The manuscript is very well written and I did not find any particular issue that should be addressed.

Possibly, authors could consider a better representation of the displacements in Figure 8 and reduce the font size of colorbars in Figure 10 (if possible). Nevertheless, such minor issues do not represent real drawbacks and could be implemented during the publication phase.

In conclusion, I recommend it to be accepted in the present form

Thank you for review

 

Point 1

Possibly, authors could consider a better representation of the displacements in Figure 8 and reduce the font size of colorbars in Figure 10 (if possible).

Response 1

Thank you for your comment.

Author Response File: Author Response.pdf

Reviewer 4 Report

1. There was no text describing Table 1. Please add a paragraph discussing the differences in the requirement of barrier wall in different countries. Also please add references for safety code requirement in Japan and in France.

2. The heading in Table 1 " Minimum requirement of barrier wall" would be better understood if it is changed to "Minimum requirement of barrier wall thickness."

3. Line 98 to 103: It is quite confusing the way it is now in the manuscript because most of the description is about the accident at hydrogen station in Gangneung, but in the middle of this, but another hydrogen at Samcheok is briefly described. It is suggested that the sentence is split into two sentences to separate the two hydrogen stations as follows: Figure 1 (a)  shows the hydrogen tank explosion that occurred in Gangneung, South Korea, on May 2019. Figure 1 (b) shows the site visit at the hydrogen station at Samcheok, South Korea.

3.1 Then move the following sentences that describe Figure 1(b) to line 108 by inserting them in front of "it is necessary... " These sentences are"Figure 1 (b) shows the site visit at the hydrogen station at Samcheok, South Korea. The thickness of the reinforced concrete is 0.35 m in figure 1 (b), but the distance between the hydrogen gas trailer and the protection wall is 1.53 m. The distance between 102 the protection wall and hydrogen gas storage follows the International Fire Code (IFC) 2009 (States, 2007) [17]."

3.2 Then please describe the purpose of showing the picture of  the Samcheok station next to the Gangneung station. Is it because they are similar in design, or something else?

4. In the text starting from Lines 207 to 212, Figure 3 was mistyped as Figure 2.

5. Line 190-191: What is X in Expression 1? How is its value determined?
6. Line 195: Please provide a justification for choosing the value of The explosion efficiency of 0.01.
7.The sentence is broken from Line 209 to 210, "
After the positive
phase the pressure becomes smaller than the ambient pressure and at
again re
turns to it called negative phase." Please rewrite with two or more sentences, it will be easier to understand.

8. Line 189: For better terminology, it is suggested that the word "ratio" is changed to " volumetric concentration" so that it would read "
a 75 % volumetric concentration of hydrogen with air.
9. Line 295: Remove the word "of."
10. Line 255: What is Malvar? Is there a reference for this?
11. Table 6, third column: Change "form" to "from".
12. Line 315: Change "Table 6" to "Table 7."
13. Line 342-351: Figure 7 appears several times in the paragraph. Should the correct figure be "8" instead of "7."
14. Is it correct to say that Section 4.2 provides displacement results without considering failure from a detonation of 69.7 kg of TNT. And, the failures are discussed in section 4.3 based on the same simulation results as in Section 4.2.
If that is the case, please provide the failure criteria used to produce results in Section 4.3.
15. In figure 10, plots of plastic strain are shown according to color. Please add description on the meaning of cells that have no color.
16. Line 355-357, 376-377: The position of TNT was 0.55 m, but the detonation point was 0.75 m. Why the two position are different? Please describe the physical reason behind this.

Author Response

Response to Reviewer 4 Comments

 

Point 1

There was no text describing Table 1. Please add a paragraph discussing the differences in the requirement of barrier wall in different countries. Also please add references for safety code requirement in Japan and in France

Response 1

Thank you for your comment.

It has been corrected.

 

Point 2

The heading in Table 1 " Minimum requirement of barrier wall" would be better understood if it is changed to "Minimum requirement of barrier wall thickness."

Response 2

Thank you for your comment.

It has been corrected.

 

Point 3

Line 98 to 103: It is quite confusing the way it is now in the manuscript because most of the description is about the accident at hydrogen station in Gangneung, but in the middle of this, but another hydrogen at Samcheok is briefly described. It is suggested that the sentence is split into two sentences to separate the two hydrogen stations as follows: Figure 1 (a)  shows the hydrogen tank explosion that occurred in Gangneung, South Korea, on May 2019. Figure 1 (b) shows the site visit at the hydrogen station at Samcheok, South Korea.

Response 3

Thank you for your comments.

It has been corrected.

 

Point 3.1

Then move the following sentences that describe Figure 1(b) to line 108 by inserting them in front of "it is necessary... " These sentences are"Figure 1 (b) shows the site visit at the hydrogen station at Samcheok, South Korea. The thickness of the reinforced concrete is 0.35 m in figure 1 (b), but the distance between the hydrogen gas trailer and the protection wall is 1.53 m. The distance between 102 the protection wall and hydrogen gas storage follows the International Fire Code (IFC) 2009 (States, 2007) [17]."

Response 3

Thank you for your comments.

It has been corrected.

 

Point 3.2

Then please describe the purpose of showing the picture of  the Samcheok station next to the Gangneung station. Is it because they are similar in design, or something else?

Response 3

Thank you for your comments.

It has been added.

 

Point 4

In the text starting from Lines 207 to 212, Figure 3 was mistyped as Figure 2.

Response 4

Thank you for your comment.

It has been corrected.

 

Point 5

Line 190-191: What is X in Expression 1? How is its value determined?

Response 5

Thank you for your comment.

It has been corrected.

 

Point 6

Line 195: Please provide a justification for choosing the value of The explosion efficiency of 0.01.

Response 6

Thank you for your comment.

It has been corrected.

 

Point 7

The sentence is broken from Line 209 to 210, "After the positive
phase the pressure becomes smaller than the ambient pressure and at
again returns to it called negative phase." Please rewrite with two or more sentences, it will be easier to understand.

Response 7

Thank you for your comment.

It has been corrected

 

Point 8

Line 189: For better terminology, it is suggested that the word "ratio" is changed to " volumetric

Response 8

Thank you for your comment.

It has been corrected.

 

Point 9

Line 295: Remove the word "of."

Response 9

Thank you for your comment.

It has been corrected.

 

Point 10

Line 255: What is Malvar? Is there a reference for this?.

Response 10

Thank you for your comment.

Reference has been added.

 

Point 11

Table 6, third column: Change "form" to "from".

Response 11

Thank you for your comment.

It has been corrected.

 

Point 12

Line 315: Change "Table 6" to "Table 7."

Response 12

Thank you for your comment.

It has been corrected

 

Point 13

Line 342-351: Figure 7 appears several times in the paragraph. Should the correct figure be "8" instead of "7."

Response 13

Thank you for your comment.

It has been corrected.

 

Point 14

Is it correct to say that Section 4.2 provides displacement results without considering failure from a detonation of 69.7 kg of TNT. And, the failures are discussed in section 4.3 based on the same simulation results as in Section 4.2.

Response 14

Thank you for your comment.

Yes, section 4.2 and section 4.3 are based on the same simulation results.

 

Point 15

In figure 10, plots of plastic strain are shown according to color. Please add description on the meaning of cells that have no color.

Response 15

Thank you for your comment.

It has been corrected.

 

Point 16

Line 355-357, 376-377: The position of TNT was 0.55 m, but the detonation point was 0.75 m. Why the two position are different? Please describe the physical reason behind this.

Response 16

Thank you for your comment.

It has been corrected.

Author Response File: Author Response.pdf

Round 2

Reviewer 1 Report

Thank you for the answers.

I still disagree with your answer to the point 1. I undrestand  that you

would like to maximise the amount af TNT by increasing the hydrogen content to 75 %. The problem is that a significant amount of hydrogen gas will be left after explosion due to a very rich mixture. 

If your domain is open, than you would need an additional time for the remaining hydrogen be mixed with air, following by secondary explosion.

Please clarify this point.

Author Response

Thank you for review

 

I still disagree with your answer to the point 1. I undrestand  that you would like to maximise the amount af TNT by increasing the hydrogen content to 75 %. The problem is that a significant amount of hydrogen gas will be left after explosion due to a very rich mixture. 

If your domain is open, than you would need an additional time for the  remaining hydrogen be mixed with air, following by secondary explosion.

Please clarify this point.

Thank you for review

 

Point 1

Response 1

Thank you for your comment.

In present of 4 % to 75 % ratio of hydrogen with air, explode can occur (from  references). In this study maximum percent was taken for all cases. TNT  equivalent weight method was used for this study. TNT explosions are  detonations, while vapor cloud explosions are deflagrations i.e., the shape and the velocity of their respective blast waves are different. TNT equivalency method is straightforward, it remains widely used. It has been added in manuscript with references.

Figure 8 represents the pressure vs time curve and for this study 10 milliseconds simulation is sufficient for the simulations. Increasing the simulation time can erode the few concrete elements after 10 milliseconds but the process is too  slow and displacements on structure has not a large difference.

For secondary explosion current study did not consider it in this analysis. Future analysis and experiment study should be done. 

Author Response File: Author Response.pdf

Back to TopTop