We present a new stabilized finite element method for Navier-Stokes equations with friction slip boundary conditions based on Brezzi-Pitkäranta stabilized method. The stability and error estimates of numerical solutions in some norms are derived for standard one-level method. Combining the techniques of two-level discretization method, we propose two-level Newton iteration method and show the stability and error estimate. Finally, the numerical experiments are given to support the theoretical results and to check the efficiency of this two-level iteration method.
Starting from the Boltzmann equation, new boundary conditions are derived to be matched with the Navier—Stokes equations, that are supposed to hold in the main body of a gas. The idea upon which this method is based goes back to Maxwell and Langmuir. Since the distribution function is supposed to be completely determined by the Navier—Stokes equations, this new set of boundary conditions extends in some sense the validity of the macroscopic equations to the transition and free molecular régimes. In fact, it is shown that the free molecular and slip flow régimes are correctly described by this method; the latter is also supposed to give a reasonable approximation for the complete range of Knudsen numbers. The new procedure is applied to different problems such as plane Couette flow, plane and cylindrical Poiseuile flow, heat transfer between parallel plates and concentric cylinders. Results are obtained and compared with the exact numerical solutions for the above-mentioned problems.
A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.
AbstractBased on full domain partition, three parallel iterative finite-element algorithms are proposed and analyzed for the Navier–Stokes equations with nonlinear slip boundary conditions. Since the nonlinear slip boundary conditions include the subdifferential property, the variational formulation of these equations is variational inequalities of the second kind. In these parallel algorithms, each subproblem is defined on a global composite mesh that is fine with size h on its subdomain and coarse with size H (H ≫ h) far away from the subdomain, and then we can solve it in parallel with other subproblems by using an existing sequential solver without extensive recoding. All of the subproblems are nonlinear and are independently solved by three kinds of iterative methods. Compared with the corresponding serial iterative finite-element algorithms, the parallel algorithms proposed in this paper can yield an approximate solution with a comparable accuracy and a substantial decrease in computational time. Contributions of this paper are as follows: (1) new parallel algorithms based on full domain partition are proposed for the Navier–Stokes equations with nonlinear slip boundary conditions; (2) nonlinear iterative methods are studied in the parallel algorithms; (3) new theoretical results about the stability, convergence and error estimates of the developed algorithms are obtained; (4) some numerical results are given to illustrate the promise of the developed algorithms.
Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method.
Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.
Two-Level Brezzi-Pitkäranta Discretization Method Based on Newton Iteration for Navier-Stokes Equations with Friction Boundary Conditions