alphaFeedback

WorksMathematics

Well-Posedness and Regularity for a Modified Liquid Crystal Model with Transport Noise

AI contributions
Not recorded for this paper
View details

Abstract

We prove well-posedness for a stochastic Lin-Liu model of modified liquid crystals, coupling an incompressible fluid to a director field under transport and multiplicative noise. On smooth bounded domains in two dimensions, we prove global well-posedness for energy-class initial data, provided the leading coefficient of the polynomial director potential is positive. On the torus, we establish local well-posedness in time-weighted LpL^p-LqL^q frameworks with Besov initial data, including regimes of low and negative regularity.

Introduction

In this chapter, we study a stochastic version of a simplified model for the dynamics of nematic liquid crystals. The model couples an incompressible fluid velocity field with a director field describing the preferred orientation of the molecules. Our aim is to formulate the system as an abstract stochastic evolution equation and to prove well-posedness results in both a two-dimensional L2L^{2}-framework and an LpL^{p}-LqL^{q} framework on the torus.

Description of the model

In the deterministic setting, the dynamics of liquid crystals are described by the Ericksen–Leslie theory, derived independently by Ericksen and Leslie from conservation laws and nematic modelling; see [6, 9]. In [10], Lin and Liu introduced a simplified version of the Ericksen–Leslie system, which retains the main coupling between the fluid velocity and the molecular orientation while neglecting some stretching and rotational effects. Throughout this chapter, uu denotes the incompressible fluid velocity, dd denotes the director field, and pp denotes the pressure. For some k∈Nk \in\mathbb{N}, k≥1k \ge1, and coefficients (ai)0≤i≤k∈Rk+1(a_i)_{0 \le i \le k} \in\mathbb{R}^{k+1}, we set

Φ(d):=(∑i=0kai∣d∣2i)d.(1)\Phi(d) := \left(\sum_{i=0}^{k} a_i |d|^{2i}\right)d. \tag*{(1)}

The deterministic Lin–Liu system considered in [10] is then given by

du=(Δu−(u⋅∇)u−∇p−∇⋅(∇d⊙∇d)) dt,dd=(Δd−(u⋅∇)d−Φ(d)) dt,∇⋅u=0,\begin{aligned} \mathrm{d}u &= \left(\Delta u - (u \cdot\nabla)u - \nabla p - \nabla\cdot(\nabla d \odot\nabla d)\right)\,\mathrm{d}t,\\ \mathrm{d}d &= \left(\Delta d - (u \cdot\nabla)d - \Phi(d)\right)\,\mathrm{d}t,\\ \nabla\cdot u &= 0, \end{aligned}

where we use the notation

∇d⊙∇d:=∑i=1d∇di⊗∇di.(2)\nabla d \odot\nabla d := \sum_{i=1}^{d} \nabla d_i \otimes\nabla d_i. \tag*{(2)}

This system can be viewed as a simplified Ericksen-Leslie model in which the stretching and rotation effects generated by the director field are neglected. This simplification makes the model closer to the incompressible Navier-Stokes equations and allows one to develop existence and uniqueness results using related techniques, although it also removes some features of the full physical theory; see, for instance, [7].

The stochastic model we study consists in the aforementioned system (Lin-Liu), where transport and multiplicative noises have been included. More explicitly, let U=ℓ2U = \ell^{2} and let WW be a cylindrical Wiener process on U\mathcal{U}, then define the transport noise terms respectively as,

∑n≥1(b1,n⋅∇)u dWnand∑n≥1(b2,n⋅∇)d dWn,(3)\sum_{n\geq1}(b_{1,n}\cdot\nabla)u\,\mathrm{d}W^{n}\quad\text{and}\quad\sum_{n\geq1}(b_{2,n}\cdot\nabla)d\,\mathrm{d}W^{n}, \tag*{(3)}

where (bi,n)n∈N∈ℓ2(Rd)(b_{i,n})_{n\in\mathbb{N}}\in\ell^{2}(\mathbb{R}^{d}), i∈{1,2}i\in\{1,2\}. Similarly, the multiplicative noise terms read

∑n≥1g1,n(u,d) dWnand∑n≥1g2,n(d) dWn,(4)\sum_{n\geq1}g_{1,n}(u,d)\,\mathrm{d}W^{n}\quad\text{and}\quad\sum_{n\geq1}g_{2,n}(d)\,\mathrm{d}W^{n}, \tag*{(4)}

with (g1,n)n∈N:Rd×Rd→ℓ2(Rd)(g_{1,n})_{n\in\mathbb{N}}:\mathbb{R}^{d}\times\mathbb{R}^{d}\to\ell^{2}(\mathbb{R}^{d}) and (g2,n)n∈N:Rd→ℓ2(Rd)(g_{2,n})_{n\in\mathbb{N}}:\mathbb{R}^{d}\to\ell^{2}(\mathbb{R}^{d}). In addition, a martingale component of pressure increment ∑n≥1πndWn\sum_{n\geq1}\pi_{n}\mathrm{d}W^{n} is introduced in the system. As for the deterministic contribution, it is defined up to a constant, and such that

∇πn=(I−P)[(b1,n⋅∇)u+g1,n(u,d)].(5)\nabla\pi_{n}=(I-\mathbb{P})\left[(b_{1,n}\cdot\nabla)u+g_{1,n}(u,d)\right]. \tag*{(5)}

where P\mathbb{P} is the Helmholtz-Leray projection defined below.

Moreover, the system is studied on a convex bounded smooth domain O⊂Rd\mathcal{O}\subset\mathbb{R}^{d} and is equipped with Dirichlet and Neumann boundary conditions on uu and dd, respectively. Thus, the system of interest is the following, that we label as (MLC).

(MLC)

{du=(Δu−(u⋅∇)u−∇p−∇⋅(∇d⊙∇d))dt+∑n≥1[(b1,n⋅∇)u+g1,n(u,d)−∇πn]dWnon O,dd=(Δd−(u⋅∇)d−Φ(d))dt+∑n≥1[(b2,n⋅∇)d+g2,n(d)]dWnon O,∇⋅u=0on O,u=0 and ∂νd=0on ∂O,u∣t=0=u0 and d∣t=0=d0.\left\{ \begin{aligned} \mathrm{d}u&=\left(\Delta u-(u\cdot\nabla)u-\nabla p-\nabla\cdot(\nabla d\odot\nabla d)\right)\mathrm{d}t+\sum_{n\geq1}\left[(b_{1,n}\cdot\nabla)u+g_{1,n}(u,d)-\nabla\pi_{n}\right]\mathrm{d}W^{n}\quad\text{on }\mathcal{O},\\ \mathrm{d}d&=\left(\Delta d-(u\cdot\nabla)d-\Phi(d)\right)\mathrm{d}t+\sum_{n\geq1}\left[(b_{2,n}\cdot\nabla)d+g_{2,n}(d)\right]\mathrm{d}W^{n}\quad\text{on }\mathcal{O},\\ \nabla\cdot u&=0\quad\text{on }\mathcal{O},\\ u&=0\ \text{and}\ \partial_{\nu}d=0\quad\text{on }\partial\mathcal{O},\\ u|_{t=0}&=u_{0}\ \text{and}\ d|_{t=0}=d_{0}. \end{aligned} \right.

The transport noise terms ∑n≥1(b1,n⋅∇)u dWn\sum_{n\geq1}(b_{1,n}\cdot\nabla)u\,\mathrm{d}W^{n} and ∑n≥1(b2,n⋅∇)d dWn\sum_{n\geq1}(b_{2,n}\cdot\nabla)d\,\mathrm{d}W^{n} can be thought of as perturbations induced by unresolved small scale dynamics, while the multiplicative noise terms ∑n≥1g1,n(u,d) dWn\sum_{n\geq1}g_{1,n}(u,d)\,\mathrm{d}W^{n} and ∑n≥1g2,n(d) dWn\sum_{n\geq1}g_{2,n}(d)\,\mathrm{d}W^{n} can be understood as external forcing contributions. Furthermore, to ensure the energy balance of the system, we assume below that the sequences (bi,n)n≥1(b_{i,n})_{n\geq1}, i∈{1,2}i\in\{1,2\}, fulfil a parabolicity condition – see [1, 4] for example.

The abstract Cauchy problem

First, we denote by

Lσ2(O)={u∈L2(O,Rd)∣div⁡u=0, u⋅ν=0 in H−1/2(∂O;R)},(6)L^{2}_{\sigma}(\mathcal{O})=\{u\in L^{2}(\mathcal{O},\mathbb{R}^{d})\mid\operatorname{div}u=0,\ u\cdot\nu=0\text{ in }H^{-1/2}(\partial\mathcal{O};\mathbb{R})\}, \tag*{(6)}
H01(O)=H01(O,Rd)∩Lσ2(O).(7)H^{1}_{0}(\mathcal{O})=H^{1}_{0}(\mathcal{O},\mathbb{R}^{d})\cap L^{2}_{\sigma}(\mathcal{O}). \tag*{(7)}

Then, we consider the spaces

H=Lσ2(O)×H1(O;Rd)andV=H01(O)×{d∈H2(O;Rd)∣∂νd=0 on ∂O},(8)H=L^{2}_{\sigma}(\mathcal{O})\times H^{1}(\mathcal{O};\mathbb{R}^{d})\quad\text{and}\quad V=H^{1}_{0}(\mathcal{O})\times\{d\in H^{2}(\mathcal{O};\mathbb{R}^{d})\mid\partial_{\nu}d=0\text{ on }\partial\mathcal{O}\}, \tag*{(8)}

equipped with the following inner products: for all (u,d),(u′,d′)∈H(u,d),(u',d')\in H,

((u,d),(u′,d′))H:=∫O(u⋅u′+d⋅d′+∇d:∇d′) dO,(9)((u,d),(u',d'))_{H}:=\int_{\mathcal{O}}(u\cdot u'+d\cdot d'+\nabla d:\nabla d')\,\mathrm{d}\mathcal{O}, \tag*{(9)}

and, for all (u,d),(u′,d′)∈V(u,d),(u',d')\in V,

((u,d),(u′,d′))V:=∫O(u⋅u′+d⋅d′+∇u:∇u′+∇d:∇d′+Δd⋅Δd′) dO.(10)((u,d),(u',d'))_{V}:=\int_{\mathcal{O}}(u\cdot u'+d\cdot d'+\nabla u:\nabla u'+\nabla d:\nabla d'+\Delta d\cdot\Delta d')\,\mathrm{d}\mathcal{O}. \tag*{(10)}

Then, we introduce the following duality pairing ⟨⋅,⋅⟩V×V′\langle\cdot,\cdot\rangle_{V\times V'}: for all (u,d)∈V(u,d)\in V and $(u',d')\in H,

⟨(u,d),(u′,d′)⟩V×V′=((u,d),(u′,d′))H.\langle(u,d),(u',d')\rangle_{V\times V'}=((u,d),(u',d'))_{H}.

As such, V⊂H⊂V′V\subset H\subset V' is a Gelfand triple.

We consider W:=(Wn)n≥1W:=(W^{n})_{n\geq1} to be an ℓ2\ell^{2}-cylindrical Brownian motion. Moreover, we denote by Y=(u,d)Y=(u,d) and write the linear operator formally as

AY=(−PΔu−Δd)AY=\begin{pmatrix}-\mathbb{P}\Delta u\\-\Delta d\end{pmatrix}

where P:L2(O;Rd)→L2(O)\mathbb{P}: L^2(\mathcal{O};\mathbb{R}^d) \to L^2(\mathcal{O}) is the orthogonal Helmholtz projection. In the variational setting, A:V→V′A: V \to V' is defined by

References

  1. [1]A. Agresti and M.C. Veraar. “Nonlinear SPDEs and maximal regularity: An extended survey”. In: Nonlinear Differ. Equ. Appl. 32 (2025), p. 123. DOI: 10.1007/s00030-025-01090-2.DOI
  2. [2]A. Agresti and M.C. Veraar. “Stochastic maximal L^p (L^q)-regularity for second order systems with periodic boundary conditions”. In: Ann. Inst. Henri Poincaré (B) Probab. Stat. 60.1 (2024), pp. 413–430. DOI: 10.1214/22-AIHP1333.DOI
  3. [3]A. Agresti and M.C. Veraar. “Stochastic Navier–Stokes equations for turbulent flows in critical spaces”. In: Commun. Math. Phys. 405.2 (2024), Paper no. 43. DOI: 10.1007/s00220-023-04867-7.DOI
  4. [4]A. Agresti et al. “The stochastic primitive equations with transport noise and turbulent pressure”. In: Stoch. Partial Differ. Equ. Anal. Comput. 12.1 (2024), pp. 53–133. DOI: 10.1007/s40072-022-00277-3.
  5. [5]D. Chae. “On the well-posedness of the Euler equations in the Triebel-Lizorkin spaces”. In: Comm. Pure Appl. Math. 55.5 (2002), pp. 654–678. DOI: 10.1002/cpa.10029.DOI
  6. [6]J.L. Ericksen. “Conservation laws for liquid crystals”. In: Trans. Soc. Rheol. 5.1 (1961), pp. 23–34. DOI: 10.1122/1.548883.DOI
  7. [7]E. Feireisl et al. “A new approach to non-isothermal models for nematic liquid crystals”. In: Arch. Ration. Mech. Anal. 205.2 (2012), pp. 651–672. DOI: 10.1007/s00205-012-0517-4.DOI
  8. [8]T. Hytönen et al. Analysis in Banach Spaces. Vol. II: Probabilistic Methods and Operator Theory. Cham: Springer, 2017. DOI: 10.1007/978-3-319-69808-3.
  9. [9]F.M. Leslie. “Some constitutive equations for liquid crystals”. In: Arch. Ration. Mech. Anal. 28.4 (1968), pp. 265–283. DOI: 10.1007/BF00251810.DOI
  10. [10]F. Lin and C. Liu. “Nonparabolic dissipative systems modeling the flow of liquid crystals”. In: Comm. Pure Appl. Math. 48.5 (1995), pp. 501–537. DOI: 10.1002/cpa.3160480503.DOI
  11. [11]M.E. Taylor. Partial differential equations. 1, Basic theory. Springer, 1996.

Paper details

vlnc:W_v47a7r22
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

Not recorded for this paper · View statement →

Looking for related papers…
Contents