We derive an evolution equation for the free-surface dynamics of a thin film of a second-grade fluid over an unsteady stretching sheet using long-wave theory. For the numerical investigation of the viscoelastic effect on the thin-film dynamics, a finite-volume approach on a uniform grid with implicit flux discretization is applied. The present results are in excellent agreement with results available in the literature for a Newtonian fluid. We observe that the fluid thins faster with the rapid stretching rate of the sheet, but the second-grade parameter delays the thinning behaviour of the liquid film.

A hybrid grid based second-order finite volume algorithm has been developed for Detached-Eddy Simulation (DES) of turbulent flows. To alleviate the effect caused by the numerical dissipation of the commonly used second order upwind schemes in implementing DES with unstructured computational fluid dynamics (CFD) algorithms, an improved second-order hybrid scheme is established through modifying the dissipation term of the standard Roe's flux-difference splitting scheme and the numerical dissipation of the scheme can be self-adapted according to the DES flow field information. By Fourier analysis, the dissipative and dispersive features of the new scheme are discussed. To validate the numerical method, DES formulations based on the two most popular background turbulence models, namely, the one equation Spalart-Allmaras (SA) turbulence model and the two equation k–ω Shear Stress Transport model (SST), have been calibrated and tested with three typical numerical examples (decay of isotropic turbulence, NACA0021 airfoil at 60° incidence and 65° swept delta wing). Computational results indicate that the issue of numerical dissipation in implementing DES can be alleviated with the hybrid scheme, the resolution for turbulence structures is significantly improved and the corresponding solutions match the experimental data better. The results demonstrate the potentiality of the present DES solver for complex geometries.

The nonlinear and weakly dispersive Serre equations contain higher-order dispersive terms. These include mixed spatial and temporal derivative flux terms which are difficult to handle numerically. These terms can be replaced by an alternative combination of equivalent temporal and spatial terms, so that the Serre equations can be written in conservation law form. The water depth and new conserved quantities are evolved using a second-order finite-volume scheme. The remaining primitive variable, the depth-averaged horizontal velocity, is obtained by solving a second-order elliptic equation using simple finite differences. Using an analytical solution and simulating the dam-break problem, the proposed scheme is shown to be accurate, simple to implement and stable for a range of problems, including flows with steep gradients. It is only slightly more computationally expensive than solving the shallow water wave equations.

We propose a positivity preserving entropy decreasing finite volume scheme for nonlinear nonlocal equations with a gradient flow structure. These properties allow for accurate computations of stationary states and long-time asymptotics demonstrated by suitably chosen test cases in which these features of the scheme are essential. The proposed scheme is able to cope with non-smooth stationary states, different time scales including metastability, as well as concentrations and self-similar behavior induced by singular nonlocal kernels. We use the scheme to explore properties of these equations beyond their present theoretical knowledge.

In this paper, the numerical study of compressible flow around an airfoil is presented. The flow is analyzed in steady state for subsonic, transonic and even in the supersonic regimes at different angles of attack. In finite-volume method convective fluxes are calculated and compared by two schemes. Modified Jameson flux scheme based on flux averaging with pressure correction is used. Modified Roe scheme which is one of the characteristics-based schemes, with modification in calculation of Jacobian matrix based on Mach number is implemented. Second-order accuracy is used with artificial dissipation to overcome numerical oscillations. The fifth-order Runge–Kutta scheme is used for time discretization. A proper boundary condition based on characteristics is applied. Numerical experiments are performed on the NACA 0012 and also NACA 4412 airfoils. The results confirm the superiority of modified upwind Roe scheme regarding the accuracy, stability and convergence. Results are compared to available results in literature and a good agreement is noticed.