Expanding on a wavelet basis the solution of an inverse problem provides several advantages. First of all, wavelet bases yield a natural and efficient multiresolution analysis which allows defining clear optimization strategies on nested subspaces of the solution space. Besides, the continuous representation of the solution with wavelets enables analytical calculation of regularization integrals over the spatial domain. By choosing differentiable wavelets, accurate high-order derivative regularizers can be efficiently designed via the basis's mass and stiffness matrices. More importantly, differential constraints on vector solutions, such as the divergencefree constraint in physics, can be nicely handled with biorthogonal wavelet bases. This paper illustrates these advantages in the particular case of fluid flow motion estimation. Numerical results on synthetic and real images of incompressible turbulence show that divergencefree wavelets and high-order regularizers are particularly relevant in this context.
We explore the potential of a formulation of the Navier-Stokes equations incorporating a random description of the small-scale velocity component. This model, established from a version of the Reynolds transport theorem adapted to a stochastic representation of the flow, gives rise to a large-scale description of the flow dynamics in which emerges an anisotropic subgrid tensor, reminiscent to the Reynolds stress tensor, together with a drift correction due to an inhomogeneous turbulence. The corresponding subgrid model, which depends on the small scales velocity variance, generalizes the Boussinesq eddy viscosity assumption. However, it is not anymore obtained from an analogy with molecular dissipation but ensues rigorously from the random modeling of the flow. This principle allows us to propose several subgrid models defined directly on the resolved flow component. We assess and compare numerically those models on a standard Green-Taylor vortex flow at Reynolds 1600. The numerical simulations, carried out with an accurate divergence-free scheme, outperform classical large-eddies formulations and provides a simple demonstration of the pertinence of the proposed large-scale modeling.
This paper deals with the Helmholtz-Hodge decomposition of a vector field in bounded domain. We present a practical algorithm to compute this decomposition in the context of divergence-free and curlfree wavelets satisfying suitable boundary conditions. The method requires the inversion of divergence-free and curl-free wavelet Gram matrices. We propose an optimal preconditioning which allows to solve the systems with a small number of iterations. Finally, numerical examples prove the accuracy and the efficiency of the method.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.