Diffusive approximation of a time-fractional Burger's equation in nonlinear acoustics

A fractional time derivative is introduced into the Burger's equation to model losses of nonlinear waves. This term amounts to a time convolution product, which greatly penalizes the numerical modeling. A diffusive representation of the fractional derivative is adopted here, replacing this nonlocal operator by a continuum of memory variables that satisfy local-in-time ordinary differential equations. Then a quadrature formula yields a system of local partial differential equations, well-suited to numerical integration. The determination of the quadrature coefficients is crucial to ensure both the well-posedness of the system and the computational efficiency of the diffusive approximation. For this purpose, optimization with constraint is shown to be a very efficient strategy. Strang splitting is used to solve successively the hyperbolic part by a shock-capturing scheme, and the diffusive part exactly. Numerical experiments are proposed to assess the efficiency of the numerical modeling, and to illustrate the effect of the fractional attenuation on the wave propagation.Comment: submitted to Siam SIA

Modelling and numerical analysis of energy-dissipating systems with nonlocal free energy

The broad objective of this thesis is to design finite-volume schemes for a family of energy-dissipating systems. All the systems studied in this thesis share a common property: they are driven by an energy that decreases as the system evolves. Such decrease is produced by a dissipation mechanism, which ensures that the system eventually reaches a steady state where the energy is minimised. The numerical schemes presented here are designed to discretely preserve the dissipation of the energy, leading to more accurate and cost-effective simulations. Most of the material in this thesis is based on the publications [16, 54, 65, 66, 243]. The research content is structured in three parts. First, Part II presents well-balanced first-, second- and high-order finite-volume schemes for a general class of hydrodynamic systems with linear and nonlinear damping. These well-balanced schemes preserve stationary states at machine precision, while discretely preserving the dissipation of the discrete free energy for first- and second-order accuracy. Second, Part III focuses on finite-volume schemes for the Cahn-Hilliard equation that unconditionally and discretely satisfy the boundedness of the phase eld and the free-energy dissipation. In addition, our Cahn-Hilliard scheme is employed as an image inpainting filter before passing damaged images into a classification neural network, leading to a significant improvement of damaged-image prediction. Third, Part IV introduces nite-volume schemes to solve stochastic gradient-flow equations. Such equations are of crucial importance within the framework of fluctuating hydrodynamics and dynamic density functional theory. The main advantages of these schemes are the preservation of non-negative densities in the presence of noise and the accurate reproduction of the statistical properties of the physical systems. All these fi nite-volume schemes are complemented with prototypical examples from relevant applications, which highlight the bene fit of our algorithms to elucidate some of the unknown analytical results.Open Acces