Abstract:In this article we develop a high order accurate method to solve the incompressible boundary layer equations in a provably stable manner.~We first derive continuous energy estimates,~and then proceed to the discrete setting.~We formulate the discrete approximation using high-order finite difference methods on summation-by-parts form and implement the boundary conditions weakly using the simultaneous approximation term method.~By applying the discrete energy method and imitating the continuous analysis,~the discrete estimate that resembles the continuous counterpart is obtained proving stability.~We also show that these newly derived boundary conditions removes the singularities associated with the null-space of the nonlinear discrete spatial operator.~Numerical experiments that verifies the high-order accuracy of the scheme and coincides with the theoretical results are presented.~The numerical results are compared with the well-known Blasius similarity solution as well as that resulting from the solution of the incompressible Navier Stokes equations.

What problem does this paper attempt to address?

This paper attempts to solve the following problems: 1. **Develop high - order accurate numerical methods**: The author aims to develop a high - order accurate method for incompressible boundary - layer equations to solve these problems stably and efficiently. Traditional low - order methods (such as second - order finite differences) have deficiencies when dealing with boundary - layer flows. For example, they require too many grid points to capture the details within the boundary layer and may introduce non - physical dissipation. 2. **Ensure numerical stability**: By deriving continuous energy estimates and discretizing them, the author hopes to prove that the proposed numerical method is stable. Specifically, they use the Summation - by - Parts (SBP) form in high - order finite - difference methods and combine the Simultaneous Approximation Term (SAT) method to weakly impose boundary conditions, thereby ensuring the stability of the numerical scheme. 3. **Eliminate singularities and improve boundary conditions**: It is mentioned in the paper that the newly derived boundary conditions can eliminate singularities related to the null space of the nonlinear discrete - space operator. This helps to improve the quality and reliability of the numerical solution. 4. **Verify high - order accuracy**: The author verifies the high - order accuracy of this method through numerical experiments and compares the results with theoretical analysis. In addition, comparisons are also made with the classical Blasius similarity solution and the results of solving the incompressible Navier - Stokes equations to confirm the effectiveness and accuracy of the new method. ### Presentation of Mathematical Formulas To ensure the correctness and readability of the formulas, some of the key formulas involved in the paper are presented in Markdown format as follows: #### Incompressible Boundary - Layer Equations \[ u_t+uu_x + vu_y=-p_x+\mu u_{yy} \] \[ 0 = -p_y \] \[ u_x + v_y=0 \] #### Discrete Energy Estimate \[ \frac{d}{dt}\|U\|^2_I + 2\mu\|U_y\|^2_I=-\oint_\Gamma\left(unu^2 + 2unp-2\mu uu_{yn}\right)ds \] #### Boundary Condition Matrix Decomposition \[ M=\bar{X}\Lambda_M\bar{X}^T \] where $\Lambda_M = \text{diag}(\lambda_1,\lambda_2,\lambda_3,\lambda_4)$ and $\bar{X}$ is the eigenvector matrix. #### Energy - Stable Boundary Conditions \[ u = 0, v = 0\quad\text{(South boundary)} \] \[ \alpha\frac{v}{2}u-\mu u_y = 0, p = p_\infty\quad\text{(North boundary)} \] \[ u = U_\infty\quad\text{(West boundary)} \] \[ p = p_\infty\quad\text{(East boundary)} \] These formulas show how the author ensures that the proposed method has both high - order accuracy and numerical stability through strict mathematical derivations and numerical experiments.

A provably stable and high-order accurate finite difference approximation for the incompressible boundary layer equations

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

A Higher Order Compact Finite Difference Algorithm for Solving the Incompressible Navier–Stokes Equations

Accuracy and Stability Enhancements in the Incompressible Finite-Volume-particle Method for Multiphase Flow Simulations.

Compact finite-difference scheme for some Sobolev type equations with Dirichlet boundary conditions

Weighted L2-stability of a Discrete Kinetic Approximation for the Incompressible Navier–Stokes Equations on Bounded Domains

High Order Boundary Conditions for High Order Finite Difference Schemes on Curvilinear Coordinates Solving Compressible Flows

A two-stage fourth order time-accurate discretization for Lax--Wendroff type flow solvers II. High order numerical boundary conditions

Exponential Boundary-Layer Approximation Space for Solving the Compressible Laminar Navier-Stokes Equations.

Application Of High-Order Compact Difference Scheme In The Computation Of Incompressible Wall-Bounded Turbulent Flows

A high-order upwind compact difference scheme for solving the streamfunction-velocity formulation of the unsteady incompressible Navier–Stokes equations

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier-Stokes: Applications to fluid-structure interaction

Asymptotically-consistent analytical solutions for the non-Newtonian Sakiadis boundary layer

A Less Time-Consuming Upwind Compact Difference Method with Adjusted Dissipation Property for Solving the Unsteady Incompressible Navier-Stokes Equations

An Arbitrarily High-Order Fully Well-balanced Hybrid Finite Element-Finite Volume Method for a One-dimensional Blood Flow Model

Penalty method for the Navier-Stokes-Fourier system with Dirichlet boundary conditions: convergence and error estimates

Boundary Scheme for a Discrete Kinetic Approximation of the Navier–Stokes Equations

Locally energy-stable finite element schemes for incompressible flow problems: Design and analysis for equal-order interpolations

Nonlinear Boundary Conditions for Energy and Entropy Stable Initial Boundary Value Problems in Computational Fluid Dynamics

Entropy-stable discontinuous Galerkin approximation with summation-by-parts property for the incompressible Navier-Stokes/Cahn-Hilliard system

An Upwind Compact Difference Scheme for Solving the Streamfunction-Velocity Formulation of the Unsteady Incompressible Navier-Stokes Equation