scholarly journals Ensemble Kalman Filter for non-conservative moving mesh solvers with a joint physics and mesh location update

Author(s):  
Christian Sampson ◽  
Alberto Carrassi ◽  
Ali Aydogdu ◽  
Chris Jones

<p>Numerical solvers using adaptive meshes can focus computational power on important regions of a model domain capturing important or unresolved physics. The adaptation can be informed by the model state, external information, or made to depend on the model physics. <br> In this latter case, one can think of the mesh configuration  as part of the model state. If observational data is to be assimilated into the model, the question of updating the mesh configuration with the physical values arises. Adaptive meshes present significant challenges when using popular ensemble Data Assimilation (DA) methods. We develop a novel strategy for ensemble-based DA for which the adaptive mesh is updated along with the physical values. This involves including the node locations as a part of the model state itself allowing them to be updated automatically at the analysis step. This poses a number of challenges which we resolve to produce an effective approach that promises to apply with some generality. We evaluate our strategy with two testbed models in 1-d comparing to a strategy that we previously developed that does not update the mesh configuration. We find updating the mesh improves the fidelity and convergence of the filter. An extensive analysis on the performance of our scheme beyond just the RMSE error is also presented.</p>

Atmosphere ◽  
2018 ◽  
Vol 9 (11) ◽  
pp. 444 ◽  
Author(s):  
Jinxi Li ◽  
Jie Zheng ◽  
Jiang Zhu ◽  
Fangxin Fang ◽  
Christopher. Pain ◽  
...  

Advection errors are common in basic terrain-following (TF) coordinates. Numerous methods, including the hybrid TF coordinate and smoothing vertical layers, have been proposed to reduce the advection errors. Advection errors are affected by the directions of velocity fields and the complexity of the terrain. In this study, an unstructured adaptive mesh together with the discontinuous Galerkin finite element method is employed to reduce advection errors over steep terrains. To test the capability of adaptive meshes, five two-dimensional (2D) idealized tests are conducted. Then, the results of adaptive meshes are compared with those of cut-cell and TF meshes. The results show that using adaptive meshes reduces the advection errors by one to two orders of magnitude compared to the cut-cell and TF meshes regardless of variations in velocity directions or terrain complexity. Furthermore, adaptive meshes can reduce the advection errors when the tracer moves tangentially along the terrain surface and allows the terrain to be represented without incurring in severe dispersion. Finally, the computational cost is analyzed. To achieve a given tagging criterion level, the adaptive mesh requires fewer nodes, smaller minimum mesh sizes, less runtime and lower proportion between the node numbers used for resolving the tracer and each wavelength than cut-cell and TF meshes, thus reducing the computational costs.


2015 ◽  
Vol 11 (A29B) ◽  
pp. 776-778
Author(s):  
Xin Wang ◽  

AbstractWe present new emission line identifications and improve the lensing reconstruction of the mass distribution of galaxy cluster Abell 2744 using the Grism Lens-Amplified Survey from Space (GLASS) spectroscopy and the Hubble Frontier Fields (HFF) imaging. We performed blind and targeted searches for faint line emitters on all objects, including the arc sample, within the field of view (FoV) of GLASS prime pointings. We report 55 high quality spectroscopic redshifts, 5 of which are for arc images. We also present an extensive analysis based on the HFF photometry, measuring the colors and photometric redshifts of all objects within the FoV, and comparing the spectroscopic and photometric redshift estimates. In order to improve the lens model of Abell 2744, we develop a rigorous algorithm to screen arc images, based on their colors and morphology, and selecting the most reliable ones to use. As a result, 25 systems (corresponding to 72 images) pass the screening process and are used to reconstruct the gravitational potential of the cluster pixellated on an adaptive mesh. The resulting total mass distribution is compared with a stellar mass map obtained from the Spitzer Frontier Fields data in order to study the relative distribution of stars and dark matter in the cluster.


2012 ◽  
Vol 17 (5) ◽  
pp. 732-748 ◽  
Author(s):  
Andrej Bugajev ◽  
Raimondas Čiegis

We consider a singular second-order boundary value problem. The differential problem is approximated by the Galerkin finite element scheme. The main goal is to compare the well known apriori Bakhvalov and Shishkin meshes with the adaptive mesh based on the aposteriori dual error estimators. Results of numerical experiments are presented.


Author(s):  
Sasuga Ito ◽  
Masato Furukawa ◽  
Yamada Kazutoyo ◽  
Kaito Manabe

Abstract Turbulence is one of the most important phenomena in fluid dynamics. In general, turbulent phenomena can be resolved more clearly with Large Eddy Simulation (LES) compared with Unsteady Reynolds Averaged Navier-Stokes (URANS), and the numerical solution shows good agreements with that based on Direct Numerical Simulation (DNS). However, more time and computational power are needed on LES than those on URANS. If possible, the ideal simulation method is that the method is able to resolve the turbulent phenomena same quality as the results based on DNS and LES with less time and less computational power same as that on URANS. This paper shows an adaptive simulation method based on URANS and Ensemble Kalman Filter (Enkf) to reproduce the flow fields based on LES. In this study, a two-dimensional turbine cascade flow has been solved with URANS and LES. The adaptive simulation method has been also applied to the cascade flow. As the results, in the flow field of URANS with the assimilated turbulence model’s parameters, the separation phenomenon and the boundary layer thickness was close to that of the time averaged LES.


2018 ◽  
Vol 245 ◽  
pp. 08004 ◽  
Author(s):  
Maria Churilova

The article is devoted to comparison of finite element marking criteria for adaptive mesh refinement while solving plane Cosserat elasticity problems. The goal is to compare the resulting adaptive meshes obtained with different marking strategies. Mesh refinement and error control is done using the functional type a posteriori error majorant. Implemented algorithms use the zero-order Raviart-Thomas approximation on triangular meshes. Four widely used marking criteria are utilized for mesh adaptation. The comparative analysis is presented for two plane-strain problems.


2012 ◽  
Vol 23 (05) ◽  
pp. 1250033
Author(s):  
H. W. ZHENG ◽  
X. J. LI ◽  
G. W. YANG ◽  
C. SHU

In this paper, the interaction of shock waves with multi-fluids interfaces is investigated by numerical simulations using unstructured quadrilateral adaptive meshes. In order to obtain a detailed structure of the interface, a solution adaptive method for compressible multi-fluid flows developed by Zheng et al. is employed. Firstly, the method is verified by a planar shock and interface interaction problem, which is compared with the front tracking method for the Richtmyer–Meshkov instability problem. Following the verification, the interaction between a circular shock and a sinusoidally perturbed circular interface in cylinder vessel is firstly investigated in our paper. The results show that the solution adaptive method can be employed to study the compressible multi-fluid cases with relatively complex geometry as well as capturing the fine details of interfacial structures of the interaction.


2017 ◽  
Vol 21 (1) ◽  
pp. 635-650 ◽  
Author(s):  
Chengcheng Huang ◽  
Andrew J. Newman ◽  
Martyn P. Clark ◽  
Andrew W. Wood ◽  
Xiaogu Zheng

Abstract. In this study, we examine the potential of snow water equivalent data assimilation (DA) using the ensemble Kalman filter (EnKF) to improve seasonal streamflow predictions. There are several goals of this study. First, we aim to examine some empirical aspects of the EnKF, namely the observational uncertainty estimates and the observation transformation operator. Second, we use a newly created ensemble forcing dataset to develop ensemble model states that provide an estimate of model state uncertainty. Third, we examine the impact of varying the observation and model state uncertainty on forecast skill. We use basins from the Pacific Northwest, Rocky Mountains, and California in the western United States with the coupled Snow-17 and Sacramento Soil Moisture Accounting (SAC-SMA) models. We find that most EnKF implementation variations result in improved streamflow prediction, but the methodological choices in the examined components impact predictive performance in a non-uniform way across the basins. Finally, basins with relatively higher calibrated model performance (> 0.80 NSE) without DA generally have lesser improvement with DA, while basins with poorer historical model performance show greater improvements.


Sign in / Sign up

Export Citation Format

Share Document