Suscripción a Biblioteca: Guest
Portal Digitalde Biblioteca Digital eLibros Revistas Referencias y Libros de Ponencias Colecciones
International Journal for Uncertainty Quantification
Factor de Impacto: 4.911 Factor de Impacto de 5 años: 3.179 SJR: 1.008 SNIP: 0.983 CiteScore™: 5.2

ISSN Imprimir: 2152-5080
ISSN En Línea: 2152-5099

Acceso abierto

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2018021021
pages 383-403

DATA ASSIMILATION FOR NAVIER-STOKES USING THE LEAST-SQUARES FINITE-ELEMENT METHOD

Alexander Schwarz
Institut für Mechanik, University of Duisburg-Essen, Universitätsstraße 15, 45141 Essen, Germany
Richard P. Dwight
Aerodynamics Group, Faculty of Aerospace, TU Delft, P.O. Box 5058, 2600GB Delft, The Netherlands

SINOPSIS

We investigate theoretically and numerically the use of the least-squares finite-element method (LSFEM) to approach data-assimilation problems for the steady-state, incompressible Navier-Stokes equations. Our LSFEM discretization is based on a stress-velocity-pressure (S-V-P) first-order formulation, using discrete counterparts of the Sobolev spaces H(div)×H1×L2 for the variables respectively. In general, S-V-P formulations are promising when the stresses are of special interest, e.g., for non-Newtonian, multiphase or turbulent flows. Resolution of the system is via minimization of a least-squares functional representing the magnitude of the residual of the equations. A simple and immediate approach to extend this solver to data assimilation is to add a data-discrepancy term to the functional. Whereas most data assimilation techniques require a large number of evaluations of the forward simulation and are therefore very expensive, the approach proposed in this work uniquely has the same cost as a single forward run. However, the question arises: what is the statistical model implied by this choice? We answer this within the Bayesian framework, establishing the latent background covariance model and the likelihood. Further we demonstrate that−in the linear case−the method is equivalent to application of the Kalman filter, and derive the posterior covariance. We practically demonstrate the capabilities of our method on a backward-facing step case. Our LSFEM formulation (without data) is shown to have good approximation quality, even on relatively coarse meshes−in particular with respect to mass conservation and reattachment location. Adding limited velocity measurements from experiment, we show that the method is able to correct for discretization error on very coarse meshes, as well as correct for the influence of unknown and uncertain boundary conditions.


Articles with similar content:

OPTIMIZATION-BASED SAMPLING IN ENSEMBLE KALMAN FILTERING
International Journal for Uncertainty Quantification, Vol.4, 2014, issue 4
Alexander Bibov, Heikki Haario, Antti Solonen, Johnathan M. Bardsley
AN ENSEMBLE KALMAN FILTER USING THE CONJUGATE GRADIENT SAMPLER
International Journal for Uncertainty Quantification, Vol.3, 2013, issue 4
Heikki Haario, Antti Solonen, Albert Parker, Marylesa Howard, Johnathan M. Bardsley
AN ADAPTIVE TWO-DIMENSIONAL MESH REFINEMENT METHOD BASED ON THE LAW OF MASS CONSERVATION
Journal of Flow Visualization and Image Processing, Vol.15, 2008, issue 1
Zhenquan Li
NONLINEAR NONLOCAL MULTICONTINUA UPSCALING FRAMEWORK AND ITS APPLICATIONS
International Journal for Multiscale Computational Engineering, Vol.16, 2018, issue 5
Eric T. Chung, Yalchin Efendiev, Wing T. Leung, Mary Wheeler
A POSTERIORI ERROR ESTIMATION FOR A CUT CELL FINITE VOLUME METHOD WITH UNCERTAIN INTERFACE LOCATION
International Journal for Uncertainty Quantification, Vol.5, 2015, issue 5
Simon Tavener, J. B. Collins, Donald Estep