噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D...
Transcript of 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D...
![Page 1: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/1.jpg)
噴煙柱数値モデルのベンチマークテスト
S. Barsotti (INGV, Italy; IMO)C. Bonadonna (Univ. Geneva, Switzerland)M. Bursik (Univ. Buffalo, USA)G. Carazzo (IPG, France)M. Cerminara (INGV, Italy)W. J. Degruyter (Georgia Tech, USA)M. de’Michieli Vitturi (INGV, Italy)L. C. Denby (Univ. Cambridge, UK)B. Devenish (Met. Office, UK)S. Engwell (INGV, Italy)A. Folch (Barcelona Supercomputing Center, Spain)
F. Girault (IPG, France)M. Herzog (Univ. Cambridge, UK)A. J. Hogg (Univ. Bristol, UK)E. Kaminski (IPG, France)G. Macedonio (INGV, Italy) L. G. Mastin (USGS, USA)A. Neri (INGV, Italy)J. C. Phillips (Univ. Bristol, UK)S. Tait (IPG, France)A. van Eaton (USGS, USA)C. Witham (Met. Office, UK)M. Woodhouse (Univ. Bristol, UK)
鈴木雄治郎(地震研)・Antonio Costa (INGV)
![Page 2: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/2.jpg)
噴煙高度
Modified from Fig. 1 in Mastin et al. (2009)
/ (Morton et al., 1956)
log(MER [kg/s])Co
lumn he
ight [km]
大気中での火山灰拡散分布地表での堆積分布
噴煙高度
噴出率 (MER [kg/s])
次元解析モデル
![Page 3: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/3.jpg)
噴煙柱モデル1D model*Morton et al. (1956)のBuoyant Plume Theoryに基づく
*定常噴煙*中心軸に沿った保存則
3D model*非定常*Navier‐Stokes 方程式
計算量:少エントレインメント仮説に基づく‐> 実用的
計算量:大流れ場の直接計算‐> ダイナミクス理解に有用
異なるモデルから得られる結果を比較モデル間の不確定性を評価
![Page 4: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/4.jpg)
様々な噴煙柱モデルLabel Name Dimension Air Entrainment Corr. Author
1 Puffin 1D = 0.15, = 1.0 M. Bursik
2 Degruyter&Bonadonna 1D = 0.10, = 0.5 W. Degruyter
3 PlumeMoM 1D = 0.09, = 0.6 M. de’Michieli Vitturi
4 Devenish 1D = 0.10, = 0.5 B. Devenish
5 FPluMe 1D = f(Ri), = g(Ri) A. Folch
6 PPM 1D = f(Ri), = 0.5 F. Girault
7 Plumeria 1D = 0.09, = 0.5 L. Mastin
8 PlumeRise 1D = 0.09, = 0.9 M. Woodhouse
9 Cerminara1D 1D = 0.10, = 0.0 M. Cerminara
10 ATHAM 3D LES M. Herzog
11 SK‐3D 3D no‐LES Y. J. Suzuki
12 ASHEE 3D LES M. Cerminara
13 PDAC 3D LES T. Esposti Ongaro
14 Mastin et al. (2009) 0D L. G. Mastin
15 Degruyter&Bonadonna(2012) 0D = 0.10, = 0.5 W. Degruyter
16 Woodhouse et al. (2013) 0D M. Woodhouse
![Page 5: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/5.jpg)
計算条件# 噴火強度 風 Fixed parameter Used models
WP2 Weak ‐ Column height (6,000 m) 1D, 0D
WP4 Weak ✔1 Column height (6,000 m) 1D, 0D
SP2 Strong ‐ Column height (37,000 m) 1D, 0D
SP4 Strong ✔2 Column height (37,000 m) 1D, 0D
WP1 Weak ‐ MER (1.5x106 kg/s) 1D, 3D, 0D
WP3 Weak ✔1 MER (1.5x106 kg/s) 1D, 3D, 0D
SP1 Strong ‐ MER (1.5x109 kg/s) 1D, 3D, 0D
SP3 Strong ✔2 MER (1.5x109 kg/s) 1D, 3D, 0D
1 2
Z (km)
W‐>E wind velocity (m/s)
Z (km)
S‐>N wind velocity (m/s) W‐>E wind velocity (m/s) S‐>N wind velocity (m/s)
Ser.A
Ser.B
![Page 6: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/6.jpg)
1Dモデルの計算例 (Plumeria)
Weak Plume without Wind
Strong Plume without Wind
![Page 7: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/7.jpg)
3Dモデルの計算例 (SK‐3D)Weak Plume
Strong Plume
![Page 8: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/8.jpg)
MERs for Fixed Column HeightsWeak plume, No wind Weak plume, With wind
Strong plume, No wind Strong plume, With wind
SD=55%
SD=47%
SD=97%
SD=46%
SD: Standard Deviation
1D
1D
1D0D 0D
0D 0D1D
Mean(1D)
Mean(1D)
Mean(1D)
Mean(1D)
MER
(kg/s)
MER
(kg/s)
MER
(kg/s)
MER
(kg/s)
![Page 9: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/9.jpg)
Column heights for Fixed MERsWeak plume, No wind Weak plume, With wind
Strong plume, No wind Strong plume, With wind
SD=9%SD=8%
SD=27%SD=21%
1D
1D 1D
1D3D
3D
3D
3D
0D 0D
0D 0DMean(1D+3D)
Mean(1D+3D)
Mean(1D+3D)
Mean(1D+3D)
Column He
ight (m
)Co
lumn He
ight (m
)
Column He
ight (m
)Co
lumn He
ight (m
)
![Page 10: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/10.jpg)
Strong plumeWindless
Windy
![Page 11: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/11.jpg)
Weak plumeWindless
Windy
![Page 12: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/12.jpg)
Inter‐comparison of 3D models (1)
SK‐3DATHAM
ASHEE PDAC
Strong plume without wind
Time averaging between 900 to 960 sec.
![Page 13: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/13.jpg)
Inter‐comparison of 3D models (2)SK‐3DATHAM
ASHEE PDAC
Weak plumewithout wind
![Page 14: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/14.jpg)
Inter‐comparison of 3D models (3)
ATHAM SK‐3DWeakplume
Strongplume
Windy cases
![Page 15: 噴煙柱数値モデルの ベンチマークテスト10 ATHAM 3D LES M. Herzog 11 SK‐3D 3D no‐LES Y.J. Suzuki 12 ASHEE 3D LES M. Cerminara 13 PDAC 3D LES T. Esposti Ongaro 14](https://reader033.fdocument.pub/reader033/viewer/2022052003/6015e0130c432f58b31cb61f/html5/thumbnails/15.jpg)
まとめ*異なるモデルで「噴出率-噴煙高度」関係を求めた・噴出率固定:噴煙高度はモデル間で良い一致・噴煙高度固定:噴出率はモデルに敏感・小規模噴煙,風速大の条件で違いは顕著
*小規模噴煙では,1Dと3Dモデルで鉛直プロファイルが類似.ただし噴煙高度は異なる
*3Dシミュレーションを用いたエントレイメント係数決定が必要
*火山灰拡散や降灰予測で噴煙柱モデルを用いる場合にはモデルによる不確定性があることを認識しなければならない.