Analysis and numerical analysis of the Helmholtz-Korteweg equation

Technical Summary: Analysis and Numerical Analysis of the Nematic Helmholtz–Korteweg Equation

Problem Statement
This work addresses the mathematical analysis and numerical discretization of the nematic Helmholtz–Korteweg equation, a fourth-order partial differential equation modeling time-harmonic acoustic wave propagation in calamitic fluids (fluids composed of rod-like molecules) exhibiting nematic order. The governing equation is given by:
$$\alpha\Delta^2 u + \beta\nabla\cdot\nabla\left(n^\top(H u)n\right) - \Delta u - k^2 u = f$$
where $u$ represents the density perturbation, $n$ is a piecewise constant unit-vector field describing the local molecular orientation, $H u$ is the Hessian of $u$, and $\alpha \gg \beta \geq 0$ are acoustic susceptibility parameters. The equation arises from modeling nematic liquid crystals as Korteweg fluids with an additional anisotropic stress term dependent on the nematic director $n$.

The problem is posed on a bounded convex domain $\Omega \subset \mathbb{R}^d$ ($d \in \{2,3\}$) with smooth boundary. The authors consider three types of boundary conditions: sound soft, sound hard, and impedance. A primary mathematical challenge identified is that the boundary conditions involve higher-order traces (specifically involving $\Delta u$ and $\nabla(n^\top(H u)n)$) which are not well-defined pointwise for standard $H^2(\Omega)$ solutions, requiring a weak definition via duality.

Methodology
The paper employs a two-pronged approach: continuous functional analysis and discrete finite element analysis.

1. Continuous Analysis (Well-posedness):
   The authors establish the existence and uniqueness of solutions for non-resonant wave numbers $k$ using the T-coercivity technique. This involves:

   • Defining an auxiliary self-adjoint eigenvalue problem by dropping the anisotropic nematic term ($\beta=0$).
   • Constructing a bijective operator $T$ that swaps the sign of eigenfunctions associated with eigenvalues smaller than $k^2$.
   • Proving that the full operator $A$ (including the $\beta$ term) inherits T-coercivity provided $\beta$ is sufficiently small.
   • Addressing the impedance boundary conditions, which pose difficulties due to unbounded terms in the continuous $H^2$ setting, the authors note that a detailed continuous analysis for this case is deferred, though the numerical method accommodates it.

2. Numerical Discretization:
   The authors propose a convergent discretization using $H^2$-conforming finite elements (specifically Argyris and Hsieh–Clough–Tocher elements).

   • Boundary Conditions: To avoid the practical difficulty of imposing essential boundary conditions ($B_0 u = 0$) directly on $C^1$-conforming spaces, Nitsche's method is employed. This allows for the weak enforcement of boundary conditions via penalty and consistency terms added to the sesquilinear form.
   • Stability: Similar to the continuous case, the discrete stability is proven by constructing a discrete T-operator based on the eigenbasis of a discrete auxiliary problem (with $\beta=0$). The authors demonstrate that the discrete operator is uniformly T-coercive for small $\beta$ and sufficiently large stabilization parameters.
   • Convergence: Under the assumption of additional regularity for the continuous solution, a best-approximation result is derived, showing that the error in the discrete solution is bounded by the best approximation error in the finite element space.

Key Contributions and Results

• Well-posedness Proof: The paper provides the first rigorous proof of existence and uniqueness for the nematic Helmholtz–Korteweg equation in 2D and 3D for sound soft and sound hard boundary conditions, valid for small anisotropy parameters $\beta$.
• Convergent Numerical Scheme: It introduces the first convergent numerical method for this equation. The method successfully handles the high regularity requirements ($H^2$) and the non-standard, higher-order boundary conditions using Nitsche's method.
• Validation of Anisotropy: Numerical simulations in 2D confirm that the equation correctly reproduces anisotropic wave propagation. Specifically, the wave speed is shown to be greater in the direction of the nematic director $n$.
• Experimental Reproduction: The authors numerically replicate the Mullen–Lüthi–Stephen experiment, demonstrating that the model captures the anisotropic speed of sound observed in nematic liquid crystals under magnetic fields.
• Tunable Resonators: Simulations of acoustic resonance in a cavity filled with nematic liquid crystals illustrate that the system's eigenfrequencies can be tuned by altering the orientation of the nematic director $n$ (e.g., via external electromagnetic fields). This suggests the potential for designing tunable acoustic resonators.

Significance and Claims
The authors position this work as a foundational step in the mathematical understanding of acoustic wave propagation in nematic fluids. They claim that the derived equation and its numerical analysis provide a valid model for acoustic propagation in nematic calamitic fluids, provided the nematic director field is undistorted ($\nabla n \equiv 0$).

The significance of the work lies in:

1. Theoretical Rigor: Establishing the well-posedness of a novel, higher-order PDE with complex boundary conditions.
2. Numerical Innovation: Overcoming the challenges of $H^2$-conforming discretization and non-standard boundary conditions to provide a practical tool for simulation.
3. Physical Validation: Partially validating the theoretical model proposed in prior work [43] by successfully reproducing known experimental phenomena (anisotropic propagation) and demonstrating the feasibility of simulating tunable acoustic devices.

The paper remains modest regarding the impedance boundary conditions, acknowledging that the continuous analysis for this case requires a different approach (potentially a mixed formulation) and is left for future work, even though the proposed numerical method can handle them. Similarly, the model assumes a fixed director field, noting that larger length scales would require coupling with elastic effects related to director distortions.