0000000000010111
AUTHOR
Maria Lukacova-medvidova
Stochastic Galerkin method for cloud simulation
AbstractWe develop a stochastic Galerkin method for a coupled Navier-Stokes-cloud system that models dynamics of warm clouds. Our goal is to explicitly describe the evolution of uncertainties that arise due to unknown input data, such as model parameters and initial or boundary conditions. The developed stochastic Galerkin method combines the space-time approximation obtained by a suitable finite volume method with a spectral-type approximation based on the generalized polynomial chaos expansion in the stochastic space. The resulting numerical scheme yields a second-order accurate approximation in both space and time and exponential convergence in the stochastic space. Our numerical results…
IMEX Finite Volume Methods for Cloud Simulation
We present new implicit-explicit (IMEX) finite volume schemes for numerical simulation of cloud dynamics. We use weakly compressible equations to describe fluid dynamics and a system of advection-diffusion-reaction equations to model cloud dynamics. In order to efficiently resolve slow dynamics we split the whole nonlinear system in a stiff linear part governing the acoustic and gravitational waves as well as diffusive effects and a non-stiff nonlinear part that models nonlinear advection effects. We use a stiffly accurate second order IMEX scheme for time discretization to approximate the stiff linear operator implicitly and the non-stiff nonlinear operator explicitly. Fast microscale clou…
A multiscale approach to liquid flows in pipes I: The single pipe
Abstract In the present paper we study the propagation of pressure waves in a barotropic flow through a pipe, with a possibly varying cross-sectional area. The basic model is the Saint–Venant system. We derive two multiscale models for the cases of weak and strong damping, respectively, which describe the time evolution of the piezometric head and the velocity. If the damping is weak, then the corresponding first-order hyperbolic system is linear but contains an additional integro-differential equation that takes into account the damping. In the case of strong damping, the system is nonlinear. The full and multiscale models are compared numerically; we also discuss results obtained by a lar…
Numerical Simulation of a Contractivity Based Multiscale Cancer Invasion Model
We present a problem-suited numerical method for a particularly challenging cancer invasion model. This model is a multiscale haptotaxis advection-reaction-diffusion system that describes the macroscopic dynamics of two types of cancer cells coupled with microscopic dynamics of the cells adhesion on the extracellular matrix. The difficulties to overcome arise from the non-constant advection and diffusion coefficients, a time delay term, as well as stiff reaction terms.
Existence of global weak solutions to the kinetic Peterlin model
Abstract We consider a class of kinetic models for polymeric fluids motivated by the Peterlin dumbbell theories for dilute polymer solutions with a nonlinear spring law for an infinitely extensible spring. The polymer molecules are suspended in an incompressible viscous Newtonian fluid confined to a bounded domain in two or three space dimensions. The unsteady motion of the solvent is described by the incompressible Navier–Stokes equations with the elastic extra stress tensor appearing as a forcing term in the momentum equation. The elastic stress tensor is defined by Kramer’s expression through the probability density function that satisfies the corresponding Fokker–Planck equation. In thi…
A Second Order Accurate Kinetic Relaxation Scheme for Inviscid Compressible Flows
In this paper we present a kinetic relaxation scheme for the Euler equations of gas dynamics in one space dimension. The method is easily applicable to solve any complex system of conservation laws. The numerical scheme is based on a relaxation approximation for conservation laws viewed as a discrete velocity model of the Boltzmann equation of kinetic theory. The discrete kinetic equation is solved by a splitting method consisting of a convection phase and a collision phase. The convection phase involves only the solution of linear transport equations and the collision phase instantaneously relaxes the distribution function to an equilibrium distribution. We prove that the first order accur…
Convergence of a finite volume scheme for the compressible Navier–Stokes system
We study convergence of a finite volume scheme for the compressible (barotropic) Navier–Stokes system. First we prove the energy stability and consistency of the scheme and show that the numerical solutions generate a dissipative measure-valued solution of the system. Then by the dissipative measure-valued-strong uniqueness principle, we conclude the convergence of the numerical solution to the strong solution as long as the latter exists. Numerical experiments for standard benchmark tests support our theoretical results.
A multi-scale method for complex flows of non-Newtonian fluids
We introduce a new heterogeneous multi-scale method for the simulation of flows of non-Newtonian fluids in general geometries and present its application to paradigmatic two-dimensional flows of polymeric fluids. Our method combines micro-scale data from non-equilibrium molecular dynamics (NEMD) with macro-scale continuum equations to achieve a data-driven prediction of complex flows. At the continuum level, the method is model-free, since the Cauchy stress tensor is determined locally in space and time from NEMD data. The modelling effort is thus limited to the identification of suitable interaction potentials at the micro-scale. Compared to previous proposals, our approach takes into acco…
Systematic derivation of hydrodynamic equations for viscoelastic phase separation
(abridged) We present a detailed derivation of a simple hydrodynamic two-fluid model, which aims at the description of the phase separation of non-entangled polymer solutions, where viscoelastic effects play a role. It is directly based upon the coarse-graining of a well-defined molecular model, such that all degrees of freedom have a clear and unambiguous molecular interpretation. The considerations are based upon a free-energy functional, and the dynamics is split into a conservative and a dissipative part, where the latter satisfies the Onsager relations and the Second Law of thermodynamics. The model is therefore fully consistent with both equilibrium and non-equilibrium thermodynamics.…
Adaptive discontinuous evolution Galerkin method for dry atmospheric flow
We present a new adaptive genuinely multidimensional method within the framework of the discontinuous Galerkin method. The discontinuous evolution Galerkin (DEG) method couples a discontinuous Galerkin formulation with approximate evolution operators. The latter are constructed using the bicharacteristics of multidimensional hyperbolic systems, such that all of the infinitely many directions of wave propagation are considered explicitly. In order to take into account multiscale phenomena that typically appear in atmospheric flows nonlinear fluxes are split into a linear part governing the acoustic and gravitational waves and a nonlinear part that models advection. Time integration is realiz…
Mathematical and Numerical Analysis of Some FSI Problems
In this chapter we deal with some specific existence and numerical results applied to a 2D/1D fluid–structure coupled model, for an incompressible fluid and a thin elastic structure. We will try to underline some of the mathematical and numerical difficulties that one may face when studying this kind of problems such as the geometrical nonlinearities or the added mass effect. In particular we will point out the link between the strategies of proof of weak or strong solutions and the possible algorithms to discretize these type of coupled problems.
Energy-Stable Numerical Schemes for Multiscale Simulations of Polymer–Solvent Mixtures
We present a new second-order energy dissipative numerical scheme to treat macroscopic equations aiming at the modeling of the dynamics of complex polymer–solvent mixtures. These partial differential equations are the Cahn-Hilliard equation for diffuse interface phase fields and the Oldroyd-B equations for the hydrodynamics of the polymeric mixture. A second-order combined finite volume/finite difference method is applied for the spatial discretization. A complementary approach to study the same physical system is realized by simulations of a microscopic model based on a hybrid Lattice Boltzmann/Molecular Dynamics scheme. These latter simulations provide initial conditions for the numerical…
Molecular dynamics simulations in hybrid particle-continuum schemes: Pitfalls and caveats
Heterogeneous multiscale methods (HMM) combine molecular accuracy of particle-based simulations with the computational efficiency of continuum descriptions to model flow in soft matter liquids. In these schemes, molecular simulations typically pose a computational bottleneck, which we investigate in detail in this study. We find that it is preferable to simulate many small systems as opposed to a few large systems, and that a choice of a simple isokinetic thermostat is typically sufficient while thermostats such as Lowe-Andersen allow for simulations at elevated viscosity. We discuss suitable choices for time steps and finite-size effects which arise in the limit of very small simulation bo…
$$\mathscr {K}$$-Convergence of Finite Volume Solutions of the Euler Equations
We review our recent results on the convergence of invariant domain-preserving finite volume solutions to the Euler equations of gas dynamics. If the classical solution exists we obtain strong convergence of numerical solutions to the classical one applying the weak-strong uniqueness principle. On the other hand, if the classical solution does not exist we adapt the well-known Prokhorov compactness theorem to space-time probability measures that are generated by the sequences of finite volume solutions and show how to obtain the strong convergence in space and time of observable quantities. This can be achieved even in the case of ill-posed Euler equations having possibly many oscillatory s…
Comparison of continuous and discontinuous Galerkin approaches for variable-viscosity Stokes flow
We describe a Discontinuous Galerkin (DG) scheme for variable-viscosity Stokes flow which is a crucial aspect of many geophysical modelling applications and conduct numerical experiments with different elements comparing the DG approach to the standard Finite Element Method (FEM). We compare the divergence-conforming lowest-order Raviart-Thomas (RT0P0) and Brezzi-Douglas-Marini (BDM1P0) element in the DG scheme with the bilinear Q1P0 and biquadratic Q2P1 elements for velocity and their matching piecewise constant/linear elements for pressure in the standard continuous Galerkin (CG) scheme with respect to accuracy and memory usage in 2D benchmark setups. We find that for the chosen geodynami…
Well-balanced bicharacteristic-based scheme for multilayer shallow water flows including wet/dry fronts
The aim of this paper is to present a new well-balanced finite volume scheme for two-dimensional multilayer shallow water flows including wet/dry fronts. The ideas, presented here for the two-layer model, can be generalized to a multilayer case in a straightforward way. The method developed here is constructed in the framework of the Finite Volume Evolution Galerkin (FVEG) schemes. The FVEG methods couple a finite volume formulation with evolution operators. The latter are constructed using the bicharacteristics of multidimensional hyperbolic systems. However, in the case of multilayer shallow water flows the required eigenstructure of the underlying equations is not readily available. Thus…
Visualization of Parameter Sensitivity of 2D Time-Dependent Flow
In this paper, we present an approach to analyze 1D parameter spaces of time-dependent flow simulation ensembles. By extending the concept of the finite-time Lyapunov exponent to the ensemble domain, i.e., to the parameter that gives rise to the ensemble, we obtain a tool for quantitative analysis of parameter sensitivity both in space and time. We exemplify our approach using 2D synthetic examples and computational fluid dynamics ensembles.
Analysis of a viscoelastic phase separation model
A new model for viscoelastic phase separation is proposed, based on a systematically derived conservative two-fluid model. Dissipative effects are included by phenomenological viscoelastic terms. By construction, the model is consistent with the second law of thermodynamics, and we study well-posedness of the model, i.e., existence of weak solutions, a weak-strong uniqueness principle, and stability with respect to perturbations, which are proven by means of relative energy estimates. A good qualitative agreement with mesoscopic simulations is observed in numerical tests.
Numerical analysis of the Oseen-type Peterlin viscoelastic model by the stabilized Lagrange-Galerkin method, Part II: A linear scheme
This is the second part of our error analysis of the stabilized Lagrange-Galerkin scheme applied to the Oseen-type Peterlin viscoelastic model. Our scheme is a combination of the method of characteristics and Brezzi-Pitk\"aranta's stabilization method for the conforming linear elements, which leads to an efficient computation with a small number of degrees of freedom especially in three space dimensions. In this paper, Part II, we apply a semi-implicit time discretization which yields the linear scheme. We concentrate on the diffusive viscoelastic model, i.e. in the constitutive equation for time evolution of the conformation tensor a diffusive effect is included. Under mild stability condi…
Energy-stable linear schemes for polymer-solvent phase field models
We present new linear energy-stable numerical schemes for numerical simulation of complex polymer-solvent mixtures. The mathematical model proposed by Zhou, Zhang and E (Physical Review E 73, 2006) consists of the Cahn-Hilliard equation which describes dynamics of the interface that separates polymer and solvent and the Oldroyd-B equations for the hydrodynamics of polymeric mixtures. The model is thermodynamically consistent and dissipates free energy. Our main goal in this paper is to derive numerical schemes for the polymer-solvent mixture model that are energy dissipative and efficient in time. To this end we will propose several problem-suited time discretizations yielding linear scheme…
Shear-Thinning in Oligomer Melts—Molecular Origins and Applications
We investigate the molecular origin of shear-thinning in melts of flexible, semiflexible and rigid oligomers with coarse-grained simulations of a sheared melt. Entanglements, alignment, stretching and tumbling modes or suppression of the latter all contribute to understanding how macroscopic flow properties emerge from the molecular level. In particular, we identify the rise and decline of entanglements with increasing chain stiffness as the major cause for the non-monotonic behaviour of the viscosity in equilibrium and at low shear rates, even for rather small oligomeric systems. At higher shear rates, chains align and disentangle, contributing to shear-thinning. By performing simulations …
A study on time discretization and adaptive mesh refinement methods for the simulation of cancer invasion: The urokinase model
In the present work we investigate a model that describes the chemotactically and proteolytically driven tissue invasion by cancer cells. The model is a system of advection-reaction-diffusion equations that takes into account the role of the serine protease urokinase-type plasminogen activator. The analytical and numerical study of such a system constitutes a challenge due to the merging, emerging, and traveling concentrations that the solutions exhibit. Classical numerical methods applied to this system necessitate very fine discretization grids to resolve these dynamics in an accurate way. To reduce the computational cost without sacrificing the accuracy of the solution, we apply adaptive…
Existence and uniqueness of global classical solutions to a two species cancer invasion haptotaxis model
We consider a haptotaxis cancer invasion model that includes two families of cancer cells. Both families, migrate on the extracellular matrix and proliferate. Moreover the model describes an epithelial-to-mesenchymal-like transition between the two families, as well as a degradation and a self-reconstruction process of the extracellular matrix. We prove positivity and conditional global existence and uniqueness of the classical solutions of the problem for large initial data.
Kinematic splitting algorithm for fluid–structure interaction in hemodynamics
Abstract In this paper we study a kinematic splitting algorithm for fluid–structure interaction problems. This algorithm belongs to the class of loosely-coupled fluid–structure interaction schemes. We will present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Fluid flow is described by the conservation laws with nonlinearities in convective and diffusive terms. For simplicity of presentation the structure is modelled by the generalized string equation, but the results presented in the paper may be generalized to more complex structure models. The arbitrary Lagrangian–Eulerian approach is used in order to take…
Entropy dissipation of moving mesh adaptation
Non-uniform grids and mesh adaptation have been a growing part of numerical simulation over the past years. It has been experimentally noted that mesh adaptation leads not only to locally improved solution but also to numerical stability of the underlying method. There have been though only few results on the mathematical analysis of these schemes due to the lack of proper tools that incorporate both the time evolution and the mesh adaptation step of the overall algorithm. In this paper we provide a method to perform the analysis of the mesh adaptation method, including both the mesh reconstruction and evolution of the solution. We moreover employ this method to extract sufficient condition…