Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDerivation and numerical approximation of hyperbolic viscoelastic flow systems: Saint-Venant 2D equations for Maxwell fluids
328
0
0.0
(
0
)
Authors
Sebastien Boyaval
Ask ChatGPT about the research
No Arabic abstract
We pursue here the development of models for complex (viscoelastic) fluids in shallow free-surface gravity flows which was initiated by [Bouchut-Boyaval, M3AS (23) 2013] for 1D (translation invariant) cases. The models we propose are hyperbolic quasilinear systems that generalize Saint-Venant shallow-water equations to incompressible Maxwell fluids. The models are compatible with a formulation of the thermo-dynamics second principle. In comparison with Saint-Venant standard shallow-water model, the momentum balance includes extra-stresses associated with an elastic potential energy in addition to a hydrostatic pressure. The extra-stresses are determined by an additional tensor variable solution to a differential equation with various possible time rates. For the numerical evaluation of solutions to Cauchy problems, we also propose explicit schemes discretizing our generalized Saint-Venant systems with Finite-Volume approximations that are entropy-consistent (under a CFL constraint) in addition to satisfy exact (discrete) mass and momentum conservation laws. In comparison with most standard viscoelastic numerical models, our discrete models can be used for any retardation-time values (i.e. in the vanishing solvent-viscosity limit). We finally illustrate our hyperbolic viscoelastic flow models numerically using computer simulations in benchmark test cases. On extending to Maxwell fluids some free-shear flow testcases that are standard benchmarks for Newtonian fluids, we first show that our (numerical) models reproduce well the viscoelastic physics, phenomenologically at least, with zero retardation-time. Moreover, with a view to quantitative evaluations, numerical results in the lid-driven cavity testcase show that, in fact, our models can be compared with standard viscoelastic flow models in sheared-flow benchmarks on adequately choosing the physical parameters of our models. Analyzing our models asymptotics should therefore shed new light on the famous High-Weissenberg Number Problem (HWNP), which is a limit for all the existing viscoelastic numerical models.
rate research
Read More
The standard multilayer Saint-Venant system consists in introducing fluid layers that are advected by the interfacial velocities. As a consequence there is no mass exchanges between these layers and each layer is described by its height and its average velocity. Here we introduce another multilayer system with mass exchanges between the neighborhing layers where the unknowns are a total height of water and an average velocity per layer. We derive it from Navier-Stokes system with an hydrostatic pressure and prove energy and hyperbolicity properties of the model. We also give a kinetic interpretation leading to effective numerical schemes with positivity and energy properties. Numerical tests show the versatility of the approach and its ability to compute recirculation cases with wind forcing.
In this work, we use the finite differences in time domain (FDTD) numerical method to compute and assess the validity of Hopf solutions, or hopfions, for the electromagnetic field equations. In these solutions, field lines form closed loops characterized by different knot topologies which are preserved during their time evolution. Hopfions have been studied extensively in the past from an analytical perspective but never, to the best of our knowledge, from a numerical approach. The implementation and validation of this technique eases the study of more complex cases of this phenomena; e.g. how these fields could interact with materials (e.g. anisotropic or non-linear), their coupling with other physical systems (e.g. plasmas), and also opens the path on their artificial generation by different means (e.g. antenna arrays or lasers).
Using equations of motion with the anisotropic dissipative term for quantum particle and quantum-mechanical commutation rules, the general Maxwell-type differential equations are derived. The direct modifications of the well-known Maxwell equations due to the medium effects (openness of the system) are discussed.
We propose a new reduced model for gravity-driven free-surface flows of shallow elastic fluids. It is obtained by an asymptotic expansion of the upper-convected Maxwell model for elastic fluids. The viscosity is assumed small (of order epsilon, the aspect ratio of the thin layer of fluid), but the relaxation time is kept finite. Additionally to the classical layer depth and velocity in shallow models, our system describes also the evolution of two scalar stresses. It has an intrinsic energy equation. The mathematical properties of the model are established, an important feature being the non-convexity of the physically relevant energy with respect to conservative variables, but the convexity with respect to the physically relevant pseudo-conservative variables. Numerical illustrations are given, based on a suitable well-balanced finite-volume discretization involving an approximate Riemann solver.
The efficient numerical integration of large-scale matrix differential equations is a topical problem in numerical analysis and of great importance in many applications. Standard numerical methods applied to such problems require an unduly amount of computing time and memory, in general. Based on a dynamical low-rank approximation of the solution, a new splitting integrator is proposed for a quite general class of stiff matrix differential equations. This class comprises differential Lyapunov and differential Riccati equations that arise from spatial discretizations of partial differential equations. The proposed integrator handles stiffness in an efficient way, and it preserves the symmetry and positive semidefiniteness of solutions of differential Lyapunov equations. Numerical examples that illustrate the benefits of this new method are given. In particular, numerical results for the efficient simulation of the weather phenomenon El Ni~no are presented.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
