BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Endre Süli

Endre Süli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath.NAmath-phmath.MPNumerical Analysiscond-mat.mtrl-scicond-mat.softmath.FAnlin.PS

Catalog footprint

What is connected

21works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

Existence of Large-Data Global Weak Solutions to Kinetic Models of Nonhomogeneous Dilute Polymeric Fluids

We prove the existence of large-data global-in-time weak solutions to a general class of coupled bead-spring chain models with finitely extensible nonlinear elastic (FENE) type spring potentials for nonhomogeneous incompressible dilute polymeric fluids in a bounded domain in $\mathbb{R}^d$, $d=2$ or $3$. The class of models under consideration involves the Navier--Stokes system with variable density, where the viscosity coefficient depends on both the density and the polymer number density, coupled to a Fokker--Planck equation with a density-dependent drag coefficient. The proof is based on combining a truncation of the probability density function with a two-stage Galerkin approximation and weak compactness and compensated compactness techniques to pass to the limits in the sequence of Galerkin approximations and in the truncation level.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical analysis of a topology optimization problem for Stokes flow

T. Borrvall and J. Petersson [Topology optimization of fluids in Stokes flow, International Journal for Numerical Methods in Fluids 41 (1) (2003) 77--107] developed the first model for topology optimization of fluids in Stokes flow. They proved the existence of minimizers in the infinite-dimensional setting and showed that a suitably chosen finite element method will converge in a weak(-*) sense to an unspecified solution. In this work, we prove novel regularity results and extend their numerical analysis. In particular, given an isolated local minimizer to the infinite-dimensional problem, we show that there exists a sequence of finite element solutions, satisfying necessary first-order optimality conditions, that strongly converges to it. We also provide the first numerical investigation into convergence rates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Random vortex dynamics via functional stochastic differential equations

In this paper we present a novel, closed three-dimensional (3D) random vortex dynamics system, which is equivalent to the Navier--Stokes equations for incompressible viscous fluid flows. The new random vortex dynamics system consists of a stochastic differential equation which is, in contrast with the two-dimensional random vortex dynamics equations, coupled with a finite-dimensional ordinary functional differential equation. This new random vortex system paves the way for devising new numerical schemes (random vortex methods) for solving three-dimensional incompressible fluid flow equations by Monte Carlo simulations. In order to derive the 3D random vortex dynamics equations, we have developed two powerful tools: the first is the duality of the conditional distributions of a couple of Taylor diffusions, which provides a path space version of integration by parts; the second is a forward type Feynman--Kac formula representing solutions to nonlinear parabolic equations in terms of functional integration. These technical tools and the underlying ideas are likely to be useful in treating other nonlinear problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite Element Approximation of a Strain-Limiting Elastic Model

We construct a finite element approximation of a strain-limiting elastic model on a bounded open domain in $\mathbb{R}^d$, $d \in \{2,3\}$. The sequence of finite element approximations is shown to exhibit strong convergence to the unique weak solution of the model. Assuming that the material parameters featuring in the model are Lipschitz-continuous, and assuming that the weak solution has additional regularity, the sequence of finite element approximations is shown to converge with a rate. An iterative algorithm is constructed for the solution of the system of nonlinear algebraic equations that arises from the finite element approximation. An appealing feature of the iterative algorithm is that it decouples the monotone and linear elastic parts of the nonlinearity in the model. In particular, our choice of piecewise constant approximation for the stress tensor (and continuous piecewise linear approximation for the displacement) allows us to compute the monotone part of the nonlinearity by solving an algebraic system with $d(d+1)/2$ unknowns independently on each element in the subdivision of the computational domain. The theoretical results are illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On incompressible heat-conducting viscoelastic rate-type fluids with stress-diffusion and purely spherical elastic response

We prove the existence of large-data global-in-time weak solutions to an evolutionary PDE system describing flows of incompressible \emph{heat-conducting} viscoelastic rate-type fluids with stress-diffusion, subject to a stick-slip boundary condition for the velocity and a homogeneous Neumann boundary condition for the extra stress tensor. In the introductory section we develop the thermodynamic foundations of the proposed model, and we document the role of thermodynamics in obtaining critical structural relations between the quantities of interest. These structural relations are then exploited in the mathematical analysis of the governing equations. In particular, the definition of weak solution is motivated by the thermodynamic basis of the model. The extra stress tensor describing the elastic response of the fluid is in our case purely spherical, which is a simplification from the physical point of view. The model nevertheless exhibits features that require novel mathematical ideas in order to deal with the technically complex structure of the associated internal energy and the more complicated forms of the corresponding entropy and energy fluxes. The paper provides the first rigorous proof of the existence of large-data global-in-time weak solutions to the governing equations for \emph{coupled thermo-mechanical processes} in viscoelastic rate-type fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Finite Element Approximation of Elliptic Homogenization Problems in Nondivergence-Form

We use uniform $W^{2,p}$ estimates to obtain corrector results for periodic homogenization problems of the form $A(x/\varepsilon):D^2 u_{\varepsilon} = f$ subject to a homogeneous Dirichlet boundary condition. We propose and rigorously analyze a numerical scheme based on finite element approximations for such nondivergence-form homogenization problems. The second part of the paper focuses on the approximation of the corrector and numerical homogenization for the case of nonuniformly oscillating coefficients. Numerical experiments demonstrate the performance of the scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Well-posedness of the fractional Zener wave equation for heterogenous viscoelastic materials

We explore the well-posedness of the fractional version of Zener's wave equation for viscoelastic solids, which is based on a constitutive law relating the stress tensor $\boldsymbolσ$ to the strain tensor $\boldsymbol\varepsilon(\bf u)$, with $\bf u$ being the displacement vector, defined by: $(1+τD_t^α) {\boldsymbolσ}=(1+ρD_t^α)[2μ{\boldsymbol\varepsilon}({\bf u})+λ\text{tr}(\boldsymbol\varepsilon(\bf u)) \bf ]$. Here $μ,λ\in\mathrm{L}^\infty(Ω)$, $μ$ is the shear modulus bounded below by a positive constant, and $λ\geq 0$ is first Lamé coefficient, $D_t^α$, with $α\in (0,1)$, is the Caputo time-derivative, $τ>0$ is the characteristic relaxation time and $ρ\geqτ$ is the characteristic retardation time. We show that, when coupled with the equation of motion $\varrho \ddot{\bf u} = \text{Div}{\boldsymbolσ} + \bf F$, considered in a bounded open Lipschitz domain $Ω$ in $\mathbb{R}^3$ and over a time interval $(0,T]$, where $\varrho\in \mathrm{L}^\infty(Ω)$ is the density of the material, bounded below by a positive constant, and $\bf F$ is a specified load vector, the resulting model is well-posed in the sense that the associated initial-boundary-value problem, with initial conditions ${\bf u}(0,\mathbf{x}) = {\bf g}(\mathbf{x})$, $\dot{\bf u}(0,\mathbf{x}) = \bf h(\mathbf{x})$, ${\boldsymbolσ}(0,\mathbf{x}) = {\bf s}(\mathbf{x})$, for $\mathbf{x} \in Ω$, and a homogeneous Dirichlet boundary condition, possesses a unique weak solution for any choice of ${\bf g }\in [\mathrm{H}^1_0(Ω)]^3$, ${\bf h}\in [\mathrm{L}^2(Ω)]^3$, and ${\bf S} = {\bf S}^{\rm T} \in [\mathrm{L}^2(Ω)]^{3 \times 3}$, and any load vector ${\bf F} \in\mathrm{L}^2(0,T;[\mathrm{L}^2(Ω)]^3)$, and that this unique weak solution depends continuously on the initial data and the load vector.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A partial Fourier transform method for a class of hypoelliptic Kolmogorov equations

We consider hypoelliptic Kolmogorov equations in $n+1$ spatial dimensions, with $n\geq 1$, where the differential operator in the first $n$ spatial variables featuring in the equation is second-order elliptic, and with respect to the $(n+1)$st spatial variable the equation contains a pure transport term only and is therefore first-order hyperbolic. If the two differential operators, in the first $n$ and in the $(n+1)$st co-ordinate directions, do not commute, we benefit from hypoelliptic regularization in time, and the solution for $t>0$ is smooth even for a Dirac initial datum prescribed at $t=0$. We study specifically the case where the coefficients depend only on the first $n$ variables. In that case, a Fourier transform in the last variable and standard central finite difference approximation in the other variables can be applied for the numerical solution. We prove second-order convergence in the spatial mesh size for the model hypoelliptic equation $\frac{\partial u}{\partial t} + x \frac{\partial u}{\partial y} = \frac{\partial^2 u}{\partial x^2}$ subject to the initial condition $u(x,y,0) = δ(x) δ(y)$, with $(x,y) \in \mathbb{R} \times\mathbb{R}$ and $t>0$, proposed by Kolmogorov, and for an extension with $n=2$. We also demonstrate exponential convergence of an approximation of the inverse Fourier transform based on the trapezium rule. Lastly, we apply the method to a PDE arising in mathematical finance, which models the distribution of the hedging error under a mis-specified derivative pricing model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Analysis of a viscosity model for concentrated polymers

The paper is concerned with a class of mathematical models for polymeric fluids, which involves the coupling of the Navier-Stokes equations for a viscous, incompressible, constant-density fluid with a parabolic-hyperbolic integro-differential equation describing the evolution of the polymer distribution function in the solvent, and a parabolic integro-differential equation for the evolution of the monomer density function in the solvent. The viscosity coefficient appearing in the balance of linear momentum equation in the Navier-Stokes system includes dependence on the shear-rate as well as on the weight-averaged polymer chain length. The system of partial differential equations under consideration captures the impact of polymerization and depolymerization effects on the viscosity of the fluid. We prove the existence of global-in-time, large-data weak solutions under fairly general hypotheses.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dissipative weak solutions to compressible Navier-Stokes-Fokker-Planck systems with variable viscosity coefficients

Motivated by a recent paper by Barrett and Süli [J.W. Barrett & E. Süli: Existence of global weak solutions to compressible isentropic finitely extensible bead-spring chain models for dilute polymers, Math. Models Methods Appl. Sci., 26 (2016)], we consider the compressible Navier--Stokes system coupled with a Fokker--Planck type equation describing the motion of polymer molecules in a viscous compressible fluid occupying a bounded spatial domain, with polymer-number-density-dependent viscosity coefficients. The model arises in the kinetic theory of dilute solutions of nonhomogeneous polymeric liquids, where the polymer molecules are idealized as bead-spring chains with finitely extensible nonlinear elastic (FENE) type spring potentials. The motion of the solvent is governed by the unsteady, compressible, barotropic Navier--Stokes system, where the viscosity coefficients in the Newtonian stress tensor depend on the polymer number density. Our goal is to show that the existence theory developed in the case of constant viscosity coefficients can be extended to the case of polymer-number-density-dependent viscosities, provided that certain technical restrictions are imposed, relating the behavior of the viscosity coefficients and the pressure for large values of the solvent density. As a first step in this direction, we prove here the weak sequential stability of the family of dissipative (finite-energy) weak solutions to the system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Existence of large-data finite-energy global weak solutions to a compressible Oldroyd-B model

A compressible Oldroyd--B type model with stress diffusion is derived from a compressible Navier--Stokes--Fokker--Planck system arising in the kinetic theory of dilute polymeric fluids, where polymer chains immersed in a barotropic, compressible, isothermal, viscous Newtonian solvent, are idealized as pairs of massless beads connected with Hookean springs. We develop a-priori bounds for the model, including a logarithmic bound, which guarantee the nonnegativity of the elastic extra stress tensor, and we prove the existence of large data global-in-time finite-energy weak solutions in two space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Degenerate Mobilities in Phase Field Models are Insufficient to Capture Surface Diffusion

Phase field models frequently provide insight to phase transitions, and are robust numerical tools to solve free boundary problems corresponding to the motion of interfaces. A body of prior literature suggests that interface motion via surface diffusion is the long-time, sharp interface limit of microscopic phase field models such as the Cahn-Hilliard equation with a degenerate mobility function. Contrary to this conventional wisdom, we show that the long-time behaviour of degenerate Cahn-Hilliard equation with a polynomial free energy undergoes coarsening, reflecting the presence of bulk diffusion, rather than pure surface diffusion. This reveals an important limitation of phase field models that are frequently used to model surface diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Existence of global weak solutions to compressible isentropic finitely extensible nonlinear bead-spring chain models for dilute polymers

We prove the existence of global-in-time weak solutions to a general class of models that arise from the kinetic theory of dilute solutions of nonhomogeneous polymeric liquids, where the polymer molecules are idealized as bead-spring chains with finitely extensible nonlinear elastic (FENE) type spring potentials. The class of models under consideration involves the unsteady, compressible, isentropic, isothermal Navier-Stokes system in a bounded domain $Ω$ in $\mathbb{R}^d$, $d = 2$ or $3$, for the density, the velocity and the pressure of the fluid. The right-hand side of the Navier-Stokes momentum equation includes an elastic extra-stress tensor, which is the sum of the classical Kramers expression and a quadratic interaction term. The elastic extra-stress tensor stems from the random movement of the polymer chains and is defined through the associated probability density function that satisfies a Fokker-Planck-type parabolic equation, a crucial feature of which is the presence of a centre-of-mass diffusion term. We require no structural assumptions on the drag term in the Fokker-Planck equation; in particular, the drag term need not be corotational. With a nonnegative initial density for the continuity equation; a square-integrable initial velocity datum for the Navier-Stokes momentum equation; and a nonnegative initial probability density function for the Fokker-Planck equation, which has finite relative entropy with respect to the Maxwellian associated with the spring potential in the model, we prove, via a limiting procedure on certain discretization and regularization parameters, the existence of a global-in-time bounded-energy weak solution to the coupled Navier-Stokes-Fokker-Planck system, satisfying the prescribed initial condition.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sharp Interface Limits of the Cahn-Hilliard Equation with Degenerate Mobility

In this work, the sharp interface limit of the degenerate Cahn-Hilliard equation (in two space dimensions) with a polynomial double well free energy and a quadratic mobility is derived via a matched asymptotic analysis involving exponentially large and small terms and multiple inner layers. In contrast to some results found in the literature, our analysis reveals that the interface motion is driven by a combination of surface diffusion flux proportional to the surface Laplacian of the interface curvature and an additional contribution from nonlinear, porous-medium type bulk diffusion, For higher degenerate mobilities, bulk diffusion is subdominant. The sharp interface models are corroborated by comparing relaxation rates of perturbations to a radially symmetric stationary state with those obtained by the phase field model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Tensor-Sparsity of Solutions to High-Dimensional Elliptic Partial Differential Equations

A recurring theme in attempts to break the curse of dimensionality in the numerical approximations of solutions to high-dimensional partial differential equations (PDEs) is to employ some form of sparse tensor approximation. Unfortunately, there are only a few results that quantify the possible advantages of such an approach. This paper introduces a class $Σ_n$ of functions, which can be written as a sum of rank-one tensors using a total of at most $n$ parameters and then uses this notion of sparsity to prove a regularity theorem for certain high-dimensional elliptic PDEs. It is shown, among other results, that whenever the right-hand side $f$ of the elliptic PDE can be approximated with a certain rate $\mathcal{O}(n^{-r})$ in the norm of ${\mathrm H}^{-1}$ by elements of $Σ_n$, then the solution $u$ can be approximated in ${\mathrm H}^1$ from $Σ_n$ to accuracy $\mathcal{O}(n^{-r'})$ for any $r'\in (0,r)$. Since these results require knowledge of the eigenbasis of the elliptic operator considered, we propose a second "basis-free" model of tensor sparsity and prove a regularity theorem for this second sparsity model as well. We then proceed to address the important question of the extent such regularity theorems translate into results on computational complexity. It is shown how this second model can be used to derive computational algorithms with performance that breaks the curse of dimensionality on certain model high-dimensional elliptic PDEs with tensor-sparse data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite element approximation of steady flows of incompressible fluids with implicit power-law-like rheology

We develop the analysis of finite element approximations of implicit power-law-like models for viscous incompressible fluids. The Cauchy stress and the symmetric part of the velocity gradient in the class of models under consideration are related by a, possibly multi--valued, maximal monotone $r$-graph, with $1<r<\infty$. Using a variety of weak compactness techniques, including Chacon's biting lemma and Young measures, we show that a subsequence of the sequence of finite element solutions converges to a weak solution of the problem as the finite element discretization parameter $h$ tends to 0. A key new technical tool in our analysis is a finite element counterpart of the Acerbi--Fusco Lipschitz truncation of Sobolev functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Greedy approximation of high-dimensional Ornstein-Uhlenbeck operators

We investigate the convergence of a nonlinear approximation method introduced by Ammar et al. (J. Non-Newtonian Fluid Mech. 139:153-176, 2006) for the numerical solution of high-dimensional Fokker-Planck equations featuring in Navier-Stokes-Fokker-Planck systems that arise in kinetic models of dilute polymers. In the case of Poisson's equation on a rectangular domain in R^2, subject to a homogeneous Dirichlet boundary condition, the mathematical analysis of the algorithm was carried out recently by Le Bris, Lelièvre and Maday (Const. Approx. 30:621-651, 2009), by exploiting its connection to greedy algorithms from nonlinear approximation theory, explored, for example, by DeVore and Temlyakov (Adv. Comput. Math. 5:173-187, 1996); hence, the variational version of the algorithm, based on the minimization of a sequence of Dirichlet energies, was shown to converge. Here, we extend the convergence analysis of the pure greedy and orthogonal greedy algorithms considered by Le Bris et al. to a technically more complicated situation, where the Laplace operator is replaced by an Ornstein-Uhlenbeck operator of the kind that appears in Fokker-Planck equations that arise in bead-spring chain type kinetic polymer models with finitely extensible nonlinear elastic potentials, posed on a high-dimensional Cartesian product configuration space D = D_1 x ... x D_N contained in R^(N d), where each set D_i, i = 1, ..., N, is a bounded open ball in R^d, d = 2, 3.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Existence and equilibration of global weak solutions to finitely extensible nonlinear bead-spring chain models for dilute polymers

We show the existence of global-in-time weak solutions to a general class of coupled FENE-type bead-spring chain models that arise from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The class of models involves the unsteady incompressible Navier-Stokes equations in a bounded domain in two or three space dimensions for the velocity and the pressure of the fluid, with an elastic extra-stress tensor appearing on the right-hand side in the momentum equation. The extra-stress tensor stems from the random movement of the polymer chains and is defined by the Kramers expression through the associated probability density function that satisfies a Fokker-Planck-type parabolic equation, a crucial feature of which is the presence of a center-of-mass diffusion term. We require no structural assumptions on the drag term in the Fokker-Planck equation; in particular, the drag term need not be corotational. With a square-integrable and divergence-free initial velocity datum for the Navier-Stokes equation and a nonnegative initial probability density function for the Fokker-Planck equation, which has finite relative entropy with respect to the Maxwellian of the model, we prove the existence of a global-in-time weak solution to the coupled Navier-Stokes-Fokker-Planck system. It is also shown that in the absence of a body force, the weak solution decays exponentially in time to the equilibrium solution, at a rate that is independent of the choice of the initial datum and of the centre-of-mass diffusion coefficient.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Existence of global weak solutions to finitely extensible nonlinear bead-spring chain models for dilute polymers with variable density and viscosity

We show the existence of global-in-time weak solutions to a general class of coupled bead-spring chain models that arise from the kinetic theory of dilute solutions of nonhomogeneous polymeric liquids with noninteracting polymer chains, with finitely extensible nonlinear elastic (FENE) spring potentials. The class of models under consideration involves the unsteady incompressible Navier-Stokes equations with variable density and density-dependent dynamic viscosity in a bounded domain in two and three space dimensions, for the density, the velocity and the pressure of the fluid, with an elastic extra-stress tensor appearing on the right-hand side in the momentum equation. The extra-stress tensor stems from the random movement of the polymer chains and is defined by the Kramers expression through the associated probability density that satisfies a Fokker-Planck-type parabolic equation, a crucial feature of which is the presence of a centre-of-mass diffusion term and a nonlinear density-dependent drag coefficient. With a bounded and positive initial density for the continuity equation; a square-integrable and divergence-free initial velocity datum for the Navier-Stokes equation; and a nonnegative initial probability density function for the Fokker-Planck equation, which has finite relative entropy with respect to the Maxwellian associated with the spring potential in the model, we prove, via a limiting procedure on certain regularization parameters, the existence of a global-in-time weak solution to the coupled Navier-Stokes-Fokker-Planck system, satisfying the given initial condition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reflections on Dubinskii's nonlinear compact embedding theorem

We present an overview of a result by Ju. A. Dubinskii [Mat. Sb. 67 (109) (1965); translated in Amer. Math. Soc. Transl. (2) 67 (1968)], concerning the compact embedding of a seminormed set in $L^p(0,T; \mathcal{A}_0)$, where $\mathcal{A}_0$ is a Banach space and $p \in [1,\infty]$; we establish a variant of Dubinskii's theorem, where a seminormed nonnegative cone is used instead of a seminormed set; and we explore the connections of these results with a nonlinear compact embedding theorem due to E. Maitre [Int. J. Math. Math. Sci. 27 (2003)].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Existence and equilibration of global weak solutions to Hookean-type bead-spring chain models for dilute polymers

We show the existence of global-in-time weak solutions to a general class of coupled Hookean-type bead-spring chain models that arise from the kinetic theory of dilute solutions of polymeric liquids with noninteracting polymer chains. The class of models involves the unsteady incompressible Navier-Stokes equations in a bounded domain in two or three space dimensions for the velocity and the pressure of the fluid, with an elastic extra-stress tensor appearing on the right-hand side in the momentum equation. The extra-stress tensor stems from the random movement of the polymer chains and is defined by the Kramers expression through the associated probability density function that satisfies a Fokker-Planck-type parabolic equation, a crucial feature of which is the presence of a center-of-mass diffusion term. We require no structural assumptions on the drag term in the Fokker-Planck equation; in particular, the drag term need not be corotational. With a square-integrable and divergence-free initial velocity datum for the Navier-Stokes equation and a nonnegative initial probability density function for the Fokker-Planck equation, which has finite relative entropy with respect to the Maxwellian of the model, we prove the existence of a global-in-time weak solution to the coupled Navier-Stokes-Fokker-Planck system. It is also shown that in the absence of a body force, the weak solution decays exponentially in time to the equilibrium solution, at a rate that is independent of the choice of the initial datum and of the centre-of-mass diffusion coefficient.

0 citations0 reviewsNo code linkedOpen access