Traveling wave solutions to a general incompressible Navier-Stokes-Fourier system with free boundary

1. Problem Formulation

The paper investigates the existence and uniqueness of traveling wave solutions for a generalized Navier-Stokes-Fourier (NSF) system modeling a viscous, incompressible, heat-conducting fluid.

• Physical Setting: The fluid occupies a "layer-like" domain $\Omega_{b+\zeta(\cdot, t)}$, which is a horizontally infinite strip of finite depth bounded below by a rigid flat surface ($\Sigma_0$) and above by a free surface ($\Sigma_{b+\zeta}$) interacting with the atmosphere.
• Governing Equations: The system couples the conservation of momentum, mass, and energy:
  $$\begin{cases} \partial_t v + v \cdot \nabla v + \nabla \cdot S = F & \text{(Momentum)} \\ \nabla \cdot v = 0 & \text{(Incompressibility)} \\ \partial_t \theta + v \cdot \nabla \theta + \nabla \cdot \Psi = 0 & \text{(Energy)} \end{cases}$$
  where $v$ is velocity, $\theta$ is temperature, $S = PI - \Pi$ is the stress tensor, and $\Psi$ is the heat flux.
• Generalization: Unlike standard NSF models, this system allows for temperature-dependent viscosity ($\Pi(\theta, Dv)$) and temperature-dependent capillary coefficients ($\varsigma(\theta)$). The latter induces Marangoni stresses on the free surface.
• Boundary Conditions:
  • Bottom ($\Sigma_0$): No-slip ($v=0$) and fixed temperature ($\theta = \theta_0$).
  • Free Surface ($\Sigma_{b+\zeta}$): Dynamic balance involving gravity, external bulk forces ($F$), external stress ($t_{ext}$), external heat sources ($h_{ext}$), surface tension, and Marangoni effects.
• Traveling Wave Ansatz: The authors seek solutions that are stationary in a frame moving with constant velocity $\gamma e_1$. This transforms the time-dependent free boundary problem into a stationary problem on a fixed domain with an unknown free surface profile $\eta(x')$.

2. Methodology

The authors employ a rigorous nonlinear functional analysis approach, specifically utilizing the Implicit Function Theorem (IFT) in Banach spaces. The proof strategy involves several sophisticated steps:

A. Reformulation in a Fixed Domain

To handle the moving free boundary, the authors introduce a diffeomorphism $F_\eta$ that maps the unknown domain $\Omega_{b+\eta}$ to a fixed reference domain $\Omega = \mathbb{R}^{n-1} \times (0, b)$. This transforms the free boundary problem into a system of PDEs with variable coefficients depending on $\eta$.

B. Linearization and the "e-Stokes" System

The core of the analysis lies in studying the linearized operator around the trivial solution (flat surface, zero velocity, constant temperature). The linearized system decouples into:

1. A generalized Stokes system for velocity and pressure.
2. A heat equation for temperature.
3. Coupling terms arising from the free boundary conditions (Marangoni stress and surface tension).

The authors define a linear operator $\Upsilon$ mapping the solution variables $(u, \psi, p, \eta)$ to the data $(f, g, l, k, h, m)$. To apply the IFT, they must prove that $\Upsilon$ is an isomorphism between appropriate Sobolev spaces.

C. Solvability of the Linearized System

Proving the isomorphism property requires establishing injectivity and surjectivity:

• Injectivity: Proven via energy estimates and uniqueness arguments for the linearized system.
• Surjectivity (The Main Challenge): This is the most technically demanding part. The authors use a duality argument involving the formal adjoint system.
  • They characterize the range of the forward operator as the annihilator of the kernel of the adjoint.
  • They derive a compatibility condition (an integral identity involving the data and the adjoint solution).
  • Fourier Analysis: By applying the horizontal Fourier transform, the compatibility condition is reduced to an algebraic equation involving Fourier symbols ($\omega_v, \omega_q, \omega_\phi$, etc.).
  • Asymptotic Analysis: The authors perform a delicate asymptotic analysis of these symbols for both low frequencies ($|\xi| \to 0$) and high frequencies ($|\xi| \to \infty$). This involves solving ODEs explicitly in the vertical direction to determine the behavior of the symbols.
  • Key Result: They show that the symbol $\rho(\xi)$ in the compatibility condition behaves like $|\xi|^2$ for small $\xi$ and $|\xi|$ for large $\xi$. This specific behavior dictates the regularity of the free surface function $\eta$, requiring it to belong to a specialized Hilbert space $X^{s+5/2}$ (identified in previous works on viscous traveling waves).

D. Nonlinear Analysis

Once the linear isomorphism is established, the authors verify that the nonlinear terms in the transformed system define a $C^1$ map between the chosen Sobolev spaces. This relies on standard algebra properties of Sobolev spaces and composition theorems for functions with bounded derivatives.

3. Key Contributions

1. Thermal Driving of Traveling Waves: A major conceptual contribution is the demonstration that heat sources alone can drive a viscous incompressible fluid into a traveling wave state. The paper proves that even if the bulk force and external stress are zero, a non-zero external heat source is sufficient to generate a unique traveling wave solution.
2. Generalized NSF System: The results extend previous work (which often assumed constant viscosity and no heat coupling) to a general system where viscosity and surface tension depend on temperature. This includes the generation of Marangoni stresses.
3. Small Data Well-Posedness: The paper establishes a local well-posedness theory for small data in high-order Sobolev spaces ($H^s$). It proves the existence of a unique solution that depends continuously (and locally Lipschitz) on the forcing data.
4. Technical Advancement in Fourier Multipliers: The paper provides a detailed analysis of the Fourier symbols associated with the coupled fluid-thermal system. The asymptotic analysis of these symbols, particularly the interaction between the thermal equation and the free boundary, is more complex than in purely hydrodynamic cases and is rigorously resolved here.

4. Main Results

• Theorem 1.1 (Flattened Formulation): Under specific conditions on the physical parameters (viscosity $\mu$, thermal conductivity $\kappa$, and the derivative of surface tension $\sigma'(0)$), specifically $\max(\|Q_1\|^2/4, \|Q_2\|^2/4)|\sigma'(0)|^2 < 2\mu\kappa$, there exist open sets of data and solutions such that the linearized operator is an isomorphism. Consequently, for small data, there exists a unique solution to the nonlinear problem in the flattened domain.
• Theorem 1.2 (Eulerian Formulation): The solutions obtained in the flattened domain can be mapped back to the original physical domain. The theorem guarantees the existence of a unique free surface $\eta$, velocity $w$, temperature $\phi$, and pressure $q$ that solve the original free boundary problem. The solution map is $C^1$ and locally Lipschitz.
• Regularity: The solutions possess high regularity ($C^k$ smoothness) provided the data is sufficiently smooth, determined by the Sobolev embedding theorems.

5. Significance

• Bridging Hydrodynamics and Thermodynamics: This work is significant because it rigorously couples the complex dynamics of free boundaries with thermal effects in a viscous fluid. Previous literature often treated these separately or assumed simplified models.
• Mechanism for Wave Generation: It identifies a new physical mechanism for generating traveling waves: thermal gradients. This has potential implications for understanding natural phenomena where heat flux drives fluid motion (e.g., in atmospheric or oceanic layers, or industrial coating processes).
• Mathematical Robustness: The paper sets a high standard for the analysis of free boundary problems with coupled physics. The techniques developed for handling the Fourier symbols and the specific compatibility conditions for the heat-conducting case provide a blueprint for analyzing other coupled systems (e.g., magnetohydrodynamics or multi-phase flows).
• Foundation for Future Research: By establishing the well-posedness of the linearized problem and the nonlinear mapping properties, this paper lays the groundwork for studying the stability, long-time behavior, and bifurcation of these thermal traveling waves.