Eintrag weiter verarbeiten
Data assimilation using adaptive, non-conservative, moving mesh models
Gespeichert in:
| Zeitschriftentitel: | Nonlinear Processes in Geophysics |
|---|---|
| Personen und Körperschaften: | , , , , |
| In: | Nonlinear Processes in Geophysics, 26, 2019, 3, S. 175-193 |
| Format: | E-Article |
| Sprache: | Englisch |
| veröffentlicht: |
Copernicus GmbH
|
| Schlagwörter: |
| author_facet |
Aydoğdu, Ali Carrassi, Alberto Guider, Colin T. Jones, Chris K. R. T Rampal, Pierre Aydoğdu, Ali Carrassi, Alberto Guider, Colin T. Jones, Chris K. R. T Rampal, Pierre |
|---|---|
| author |
Aydoğdu, Ali Carrassi, Alberto Guider, Colin T. Jones, Chris K. R. T Rampal, Pierre |
| spellingShingle |
Aydoğdu, Ali Carrassi, Alberto Guider, Colin T. Jones, Chris K. R. T Rampal, Pierre Nonlinear Processes in Geophysics Data assimilation using adaptive, non-conservative, moving mesh models General Mathematics |
| author_sort |
aydoğdu, ali |
| spelling |
Aydoğdu, Ali Carrassi, Alberto Guider, Colin T. Jones, Chris K. R. T Rampal, Pierre 1607-7946 Copernicus GmbH General Mathematics http://dx.doi.org/10.5194/npg-26-175-2019 <jats:p>Abstract. Numerical models solved on adaptive moving meshes have become increasingly prevalent in recent years. Motivating problems include the study of fluids in a Lagrangian frame and the presence of highly localized structures such as shock waves or interfaces. In the former case, Lagrangian solvers move the nodes of the mesh with the dynamical flow; in the latter, mesh resolution is increased in the proximity of the localized structure. Mesh adaptation can include remeshing, a procedure that adds or removes mesh nodes according to specific rules reflecting constraints in the numerical solver. In this case, the number of mesh nodes will change during the integration and, as a result, the dimension of the model's state vector will not be conserved. This work presents a novel approach to the formulation of ensemble data assimilation (DA) for models with this underlying computational structure. The challenge lies in the fact that remeshing entails a different state space dimension across members of the ensemble, thus impeding the usual computation of consistent ensemble-based statistics. Our methodology adds one forward and one backward mapping step before and after the ensemble Kalman filter (EnKF) analysis, respectively. This mapping takes all the ensemble members onto a fixed, uniform reference mesh where the EnKF analysis can be performed. We consider a high-resolution (HR) and a low-resolution (LR) fixed uniform reference mesh, whose resolutions are determined by the remeshing tolerances. This way the reference meshes embed the model numerical constraints and are also upper and lower uniform meshes bounding the resolutions of the individual ensemble meshes. Numerical experiments are carried out using 1-D prototypical models: Burgers and Kuramoto–Sivashinsky equations and both Eulerian and Lagrangian synthetic observations. While the HR strategy generally outperforms that of LR, their skill difference can be reduced substantially by an optimal tuning of the data assimilation parameters. The LR case is appealing in high dimensions because of its lower computational burden. Lagrangian observations are shown to be very effective in that fewer of them are able to keep the analysis error at a level comparable to the more numerous observers for the Eulerian case. This study is motivated by the development of suitable EnKF strategies for 2-D models of the sea ice that are numerically solved on a Lagrangian mesh with remeshing.</jats:p> Data assimilation using adaptive, non-conservative, moving mesh models Nonlinear Processes in Geophysics |
| doi_str_mv |
10.5194/npg-26-175-2019 |
| facet_avail |
Online Free |
| format |
ElectronicArticle |
| fullrecord |
blob:ai-49-aHR0cDovL2R4LmRvaS5vcmcvMTAuNTE5NC9ucGctMjYtMTc1LTIwMTk |
| id |
ai-49-aHR0cDovL2R4LmRvaS5vcmcvMTAuNTE5NC9ucGctMjYtMTc1LTIwMTk |
| institution |
DE-Zwi2 DE-D161 DE-Zi4 DE-Gla1 DE-15 DE-Pl11 DE-Rs1 DE-14 DE-105 DE-Ch1 DE-L229 DE-D275 DE-Bn3 DE-Brt1 |
| imprint |
Copernicus GmbH, 2019 |
| imprint_str_mv |
Copernicus GmbH, 2019 |
| issn |
1607-7946 |
| issn_str_mv |
1607-7946 |
| language |
English |
| mega_collection |
Copernicus GmbH (CrossRef) |
| match_str |
aydogdu2019dataassimilationusingadaptivenonconservativemovingmeshmodels |
| publishDateSort |
2019 |
| publisher |
Copernicus GmbH |
| recordtype |
ai |
| record_format |
ai |
| series |
Nonlinear Processes in Geophysics |
| source_id |
49 |
| title |
Data assimilation using adaptive, non-conservative, moving mesh models |
| title_unstemmed |
Data assimilation using adaptive, non-conservative, moving mesh models |
| title_full |
Data assimilation using adaptive, non-conservative, moving mesh models |
| title_fullStr |
Data assimilation using adaptive, non-conservative, moving mesh models |
| title_full_unstemmed |
Data assimilation using adaptive, non-conservative, moving mesh models |
| title_short |
Data assimilation using adaptive, non-conservative, moving mesh models |
| title_sort |
data assimilation using adaptive, non-conservative, moving mesh models |
| topic |
General Mathematics |
| url |
http://dx.doi.org/10.5194/npg-26-175-2019 |
| publishDate |
2019 |
| physical |
175-193 |
| description |
<jats:p>Abstract. Numerical models solved on adaptive moving meshes have become increasingly prevalent in recent years. Motivating problems include the study of fluids in a Lagrangian frame and the presence of highly localized structures such as shock waves or interfaces. In the former case, Lagrangian solvers move the nodes of the mesh with the dynamical flow; in the latter, mesh resolution is increased in the proximity of the localized structure. Mesh adaptation can include remeshing, a procedure that adds or removes mesh nodes according to specific rules reflecting constraints in the numerical solver. In this case, the number of mesh nodes will change during the integration and, as a result, the dimension of the model's state vector will not be conserved. This work presents a novel approach to the formulation of ensemble data assimilation (DA) for models with this underlying computational structure. The challenge lies in the fact that remeshing entails a different state space dimension across members of the ensemble, thus impeding the usual computation of consistent ensemble-based statistics. Our methodology adds one forward and one backward mapping step before and after the ensemble Kalman filter (EnKF) analysis, respectively. This mapping takes all the ensemble members onto a fixed, uniform reference mesh where the EnKF analysis can be performed. We consider a high-resolution (HR) and a low-resolution (LR) fixed uniform reference mesh, whose resolutions are determined by the remeshing tolerances. This way the reference meshes embed the model numerical constraints and are also upper and lower uniform meshes bounding the resolutions of the individual ensemble meshes. Numerical experiments are carried out using 1-D prototypical models: Burgers and Kuramoto–Sivashinsky equations and both Eulerian and Lagrangian synthetic observations. While the HR strategy generally outperforms that of LR, their skill difference can be reduced substantially by an optimal tuning of the data assimilation parameters. The LR case is appealing in high dimensions because of its lower computational burden. Lagrangian observations are shown to be very effective in that fewer of them are able to keep the analysis error at a level comparable to the more numerous observers for the Eulerian case. This study is motivated by the development of suitable EnKF strategies for 2-D models of the sea ice that are numerically solved on a Lagrangian mesh with remeshing.</jats:p> |
| container_issue |
3 |
| container_start_page |
175 |
| container_title |
Nonlinear Processes in Geophysics |
| container_volume |
26 |
| format_de105 |
Article, E-Article |
| format_de14 |
Article, E-Article |
| format_de15 |
Article, E-Article |
| format_de520 |
Article, E-Article |
| format_de540 |
Article, E-Article |
| format_dech1 |
Article, E-Article |
| format_ded117 |
Article, E-Article |
| format_degla1 |
E-Article |
| format_del152 |
Buch |
| format_del189 |
Article, E-Article |
| format_dezi4 |
Article |
| format_dezwi2 |
Article, E-Article |
| format_finc |
Article, E-Article |
| format_nrw |
Article, E-Article |
| _version_ |
1792347573795684360 |
| geogr_code |
not assigned |
| last_indexed |
2024-03-01T17:57:25.006Z |
| geogr_code_person |
not assigned |
| openURL |
url_ver=Z39.88-2004&ctx_ver=Z39.88-2004&ctx_enc=info%3Aofi%2Fenc%3AUTF-8&rfr_id=info%3Asid%2Fvufind.svn.sourceforge.net%3Agenerator&rft.title=Data+assimilation+using+adaptive%2C+non-conservative%2C+moving+mesh+models&rft.date=2019-07-24&genre=article&issn=1607-7946&volume=26&issue=3&spage=175&epage=193&pages=175-193&jtitle=Nonlinear+Processes+in+Geophysics&atitle=Data+assimilation+using+adaptive%2C+non-conservative%2C+moving+mesh+models&aulast=Rampal&aufirst=Pierre&rft_id=info%3Adoi%2F10.5194%2Fnpg-26-175-2019&rft.language%5B0%5D=eng |
| SOLR | |
| _version_ | 1792347573795684360 |
| author | Aydoğdu, Ali, Carrassi, Alberto, Guider, Colin T., Jones, Chris K. R. T, Rampal, Pierre |
| author_facet | Aydoğdu, Ali, Carrassi, Alberto, Guider, Colin T., Jones, Chris K. R. T, Rampal, Pierre, Aydoğdu, Ali, Carrassi, Alberto, Guider, Colin T., Jones, Chris K. R. T, Rampal, Pierre |
| author_sort | aydoğdu, ali |
| container_issue | 3 |
| container_start_page | 175 |
| container_title | Nonlinear Processes in Geophysics |
| container_volume | 26 |
| description | <jats:p>Abstract. Numerical models solved on adaptive moving meshes have become increasingly prevalent in recent years. Motivating problems include the study of fluids in a Lagrangian frame and the presence of highly localized structures such as shock waves or interfaces. In the former case, Lagrangian solvers move the nodes of the mesh with the dynamical flow; in the latter, mesh resolution is increased in the proximity of the localized structure. Mesh adaptation can include remeshing, a procedure that adds or removes mesh nodes according to specific rules reflecting constraints in the numerical solver. In this case, the number of mesh nodes will change during the integration and, as a result, the dimension of the model's state vector will not be conserved. This work presents a novel approach to the formulation of ensemble data assimilation (DA) for models with this underlying computational structure. The challenge lies in the fact that remeshing entails a different state space dimension across members of the ensemble, thus impeding the usual computation of consistent ensemble-based statistics. Our methodology adds one forward and one backward mapping step before and after the ensemble Kalman filter (EnKF) analysis, respectively. This mapping takes all the ensemble members onto a fixed, uniform reference mesh where the EnKF analysis can be performed. We consider a high-resolution (HR) and a low-resolution (LR) fixed uniform reference mesh, whose resolutions are determined by the remeshing tolerances. This way the reference meshes embed the model numerical constraints and are also upper and lower uniform meshes bounding the resolutions of the individual ensemble meshes. Numerical experiments are carried out using 1-D prototypical models: Burgers and Kuramoto–Sivashinsky equations and both Eulerian and Lagrangian synthetic observations. While the HR strategy generally outperforms that of LR, their skill difference can be reduced substantially by an optimal tuning of the data assimilation parameters. The LR case is appealing in high dimensions because of its lower computational burden. Lagrangian observations are shown to be very effective in that fewer of them are able to keep the analysis error at a level comparable to the more numerous observers for the Eulerian case. This study is motivated by the development of suitable EnKF strategies for 2-D models of the sea ice that are numerically solved on a Lagrangian mesh with remeshing.</jats:p> |
| doi_str_mv | 10.5194/npg-26-175-2019 |
| facet_avail | Online, Free |
| format | ElectronicArticle |
| format_de105 | Article, E-Article |
| format_de14 | Article, E-Article |
| format_de15 | Article, E-Article |
| format_de520 | Article, E-Article |
| format_de540 | Article, E-Article |
| format_dech1 | Article, E-Article |
| format_ded117 | Article, E-Article |
| format_degla1 | E-Article |
| format_del152 | Buch |
| format_del189 | Article, E-Article |
| format_dezi4 | Article |
| format_dezwi2 | Article, E-Article |
| format_finc | Article, E-Article |
| format_nrw | Article, E-Article |
| geogr_code | not assigned |
| geogr_code_person | not assigned |
| id | ai-49-aHR0cDovL2R4LmRvaS5vcmcvMTAuNTE5NC9ucGctMjYtMTc1LTIwMTk |
| imprint | Copernicus GmbH, 2019 |
| imprint_str_mv | Copernicus GmbH, 2019 |
| institution | DE-Zwi2, DE-D161, DE-Zi4, DE-Gla1, DE-15, DE-Pl11, DE-Rs1, DE-14, DE-105, DE-Ch1, DE-L229, DE-D275, DE-Bn3, DE-Brt1 |
| issn | 1607-7946 |
| issn_str_mv | 1607-7946 |
| language | English |
| last_indexed | 2024-03-01T17:57:25.006Z |
| match_str | aydogdu2019dataassimilationusingadaptivenonconservativemovingmeshmodels |
| mega_collection | Copernicus GmbH (CrossRef) |
| physical | 175-193 |
| publishDate | 2019 |
| publishDateSort | 2019 |
| publisher | Copernicus GmbH |
| record_format | ai |
| recordtype | ai |
| series | Nonlinear Processes in Geophysics |
| source_id | 49 |
| spelling | Aydoğdu, Ali Carrassi, Alberto Guider, Colin T. Jones, Chris K. R. T Rampal, Pierre 1607-7946 Copernicus GmbH General Mathematics http://dx.doi.org/10.5194/npg-26-175-2019 <jats:p>Abstract. Numerical models solved on adaptive moving meshes have become increasingly prevalent in recent years. Motivating problems include the study of fluids in a Lagrangian frame and the presence of highly localized structures such as shock waves or interfaces. In the former case, Lagrangian solvers move the nodes of the mesh with the dynamical flow; in the latter, mesh resolution is increased in the proximity of the localized structure. Mesh adaptation can include remeshing, a procedure that adds or removes mesh nodes according to specific rules reflecting constraints in the numerical solver. In this case, the number of mesh nodes will change during the integration and, as a result, the dimension of the model's state vector will not be conserved. This work presents a novel approach to the formulation of ensemble data assimilation (DA) for models with this underlying computational structure. The challenge lies in the fact that remeshing entails a different state space dimension across members of the ensemble, thus impeding the usual computation of consistent ensemble-based statistics. Our methodology adds one forward and one backward mapping step before and after the ensemble Kalman filter (EnKF) analysis, respectively. This mapping takes all the ensemble members onto a fixed, uniform reference mesh where the EnKF analysis can be performed. We consider a high-resolution (HR) and a low-resolution (LR) fixed uniform reference mesh, whose resolutions are determined by the remeshing tolerances. This way the reference meshes embed the model numerical constraints and are also upper and lower uniform meshes bounding the resolutions of the individual ensemble meshes. Numerical experiments are carried out using 1-D prototypical models: Burgers and Kuramoto–Sivashinsky equations and both Eulerian and Lagrangian synthetic observations. While the HR strategy generally outperforms that of LR, their skill difference can be reduced substantially by an optimal tuning of the data assimilation parameters. The LR case is appealing in high dimensions because of its lower computational burden. Lagrangian observations are shown to be very effective in that fewer of them are able to keep the analysis error at a level comparable to the more numerous observers for the Eulerian case. This study is motivated by the development of suitable EnKF strategies for 2-D models of the sea ice that are numerically solved on a Lagrangian mesh with remeshing.</jats:p> Data assimilation using adaptive, non-conservative, moving mesh models Nonlinear Processes in Geophysics |
| spellingShingle | Aydoğdu, Ali, Carrassi, Alberto, Guider, Colin T., Jones, Chris K. R. T, Rampal, Pierre, Nonlinear Processes in Geophysics, Data assimilation using adaptive, non-conservative, moving mesh models, General Mathematics |
| title | Data assimilation using adaptive, non-conservative, moving mesh models |
| title_full | Data assimilation using adaptive, non-conservative, moving mesh models |
| title_fullStr | Data assimilation using adaptive, non-conservative, moving mesh models |
| title_full_unstemmed | Data assimilation using adaptive, non-conservative, moving mesh models |
| title_short | Data assimilation using adaptive, non-conservative, moving mesh models |
| title_sort | data assimilation using adaptive, non-conservative, moving mesh models |
| title_unstemmed | Data assimilation using adaptive, non-conservative, moving mesh models |
| topic | General Mathematics |
| url | http://dx.doi.org/10.5194/npg-26-175-2019 |