Numerical solution of transonic and subsonic flow of compressible inviscid and viscous fluid

The work presents two numerical solutions of compressible flows problems with high and very low Mach numbers. Both problems are numerically solved by finite volume method and the explicit MacCormack scheme using a grid of quadrilateral cells. Moved grid of quadrilateral cells is considered in the form of conservation laws using Arbitrary Lagrangian–Eulerian method. In the first case, inviscid transonic flow through cascade DCA 8% is presented and the numerical results are compared to experimental data. The second case, numerical solution of unsteady viscous flow in the channel for upstream Mach number M_∞=0.012 and frequency of the wall motions 100 Hz is presented. The unsteady case can represent a simplified model of airflow coming from the trachea, through the glottal region with periodically vibrating vocal folds to the human vocal tract. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A differential quadrature analysis of unsteady open channel flow

A rapid, convergent and accurate differential quadrature method (DQM) is employed for numerical simulation of unsteady open channel flow. To the best of authors’ knowledge, this is the first attempt to use the DQM in open channel hydraulics. The Saint-Venant equations and the related nonhomogenous, time dependent boundary conditions are discretized in spatial and temporal domain by DQ rules. The unknowns in the entire domain are computed by satisfying governing equations, boundary and initial conditions simultaneously. By employing DQM, accurate results can be obtained using dramatically less grid points in spatial and time domain. The stability of DQM solution is not sensitive to choosing time step or Courant number unlike other methods. Although numerical problems such as instability, oscillation and underestimation near critical depth can be seen by using other methods but DQM solution is smooth and accurate in this case. The results are sensitive to grid distribution in time domain. In light of this, Chebyshev–Gauss–Lobatto distribution performance is excellent. To validate the DQM solutions, the obtained results are compared with those of the characteristic method. In conclusion, DQM is a potential powerful method with minimum computational effort for unsteady flow simulation.     

Uniform grid approximation of nonsmooth solutions with improved convergence for a singularly perturbed convection-diffusion equation with characteristic layers

A mixed boundary value problem for a singularly perturbed elliptic convection-diffusion equation with constant coefficients in a square domain is considered. Dirichlet conditions are specified on two sides orthogonal to the flow, and Neumann conditions are set on the other two sides. The right-hand side and the boundary functions are assumed to be sufficiently smooth, which ensures the required smoothness of the desired solution in the domain, except for neighborhoods of the corner points. Only zero-order compatibility conditions are assumed to hold at the corner points. The problem is solved numerically by applying an inhomogeneous monotone difference scheme on a rectangular piecewise uniform Shishkin mesh. The inhomogeneity of the scheme lies in that the approximating difference equations are not identical at different grid nodes but depend on the perturbation parameter. Under the assumptions made, the numerical solution is proved to converge ?-uniformly to the exact solution in a discrete uniform metric at an O(N ^?3/2ln² N) rate, where N is the number of grid nodes in each coordinate direction.     

A variable boundary method for modelling two dimensional free surface flows with moving boundaries

A coordinate transformation method is proposed for modelling unsteady, depth-averaged shallow water equations for a open channel flow with moving lateral boundaries. The transformation technique which maps the changing domain onto a fixed domain and solves the governing equations in the mapped domain, facilitates the numerical treatment of an irregular boundary configuration. The transformed equations are solved on a staggered grid with a conditionally stable, explicit finite difference scheme. Several numerical experiments are carried out corresponding to different situations, viz., flow with constant discharge, flow with constant discharge and a closed boundary at the downstream, flow in a converging channel with constant discharge and finally flow with varying discharge. The experiments are used to verify the model ability to predict free surface elevation, circulatory pattern and displacement of the boundaries. The simulated results such as displaced area, depth, displacement–time and flow-field are used to evaluate the effects of excess discharge at the upstream on the movement of lateral boundaries.     

Direct simulation of the turbulent boundary layer on a plate

A numerical method for the integration of three-dimensional Navier–Stokes equations for compressible fluid as applied to direct numerical simulation is proposed. By way of example, the boundary layer on a plate is simulated. The computations were carried out for Re_θ = 1500. The computational grid consisted of a half billion nodes. The flow region includes the laminar, transitional, and turbulent zones. The numerically obtained distributions of average velocity, friction, and pulsations are compared with experimental data and available numerical solutions.     

Laminar flow separation in an axi-symmetric sudden smooth expanded circular tube

This paper concerns with the investigation of laminar flow separation and its consequences in a tube over a smooth expansion under the axi-symmetric approximations. A co-ordinate stretching has been made to map the expanded tube into a straight tube. The two-dimensional unsteady Navier-Stokes equations are solved approximately by using primitive variables in staggered grid. A thorough quantitative analysis is performed through numerical simulations of the desired quantities such as wall shear stress, axial velocity, pressure distribution etc. These quantities are presented graphically and their consequences in the flow field are analysed in details. The dependence of the flow field on the physical parameter like expansion height d and on the Reynolds number has been investigated in details. It is interesting to note that the peak value of wall shear stress decreases with increasing height of expansion and also with the increasing Reynolds number.     

Compact Difference Scheme with High Accuracy for One-Dimensional Unsteady Quasi-Linear Biharmonic Problem of Second Kind: Application to Physical Problems

The present paper uses a new two-level implicit difference formula for the numerical study of one-dimensional unsteady biharmonic equation with appropriate initial and boundary conditions. The proposed difference scheme is second-order accurate in time and third-order accurate in space on non-uniform grid and in case of uniform mesh, it is of order two in time and four in space. The approximate solutions are computed without using any transformation and linearization. The simplicity of the proposed scheme lies in its three-point spatial discretization that yields block tri-diagonal matrix structure without the use of any fictitious nodes for handling the boundary conditions. The proposed scheme is directly applicable to singular problems, which is the main utility of our work. The method is shown to be unconditionally stable for model linear problem for uniform mesh. The efficacy of the proposed approach has been tested on several physical problems, including the complex fourth-order nonlinear equations like Kuramoto–Sivashinsky equation and extended Fisher–Kolmogorov equation, where comparison is done with some earlier work. It is clear from numerical experiments that the obtained results are not only in good agreement with the exact solutions but also competent with the solutions derived in earlier research studies.     

A semi‐discrete defect correction finite element method for unsteady incompressible magnetohydrodynamics equations

In this report, we give a semi‐discrete defect correction finite element method for the unsteady incompressible magnetohydrodynamics equations. The defect correction method is an iterative improvement technique for increasing the accuracy of a numerical solution without applying a grid refinement. Firstly, the nonlinear magnetohydrodynamics equations is solved with an artificial viscosity term. Then, the numerical solutions are improved on the same grid by a linearized defect‐correction technique. Then, we give the numerical analysis including stability analysis and error analysis. The numerical analysis proves that our method is stable and has an optimal convergence rate. In order to show the effect of our method, some numerical results are shown. Copyright © 2017 John Wiley & Sons, Ltd.     

Matrix-differential methods of solving boundary-value problems for the nonlinear one-dimensional heat-conduction equation

Matrix numerical differentiation algorithms are applied to construct numerical-analytical methods for approximate solution of boundary-value problems for the nonlinear one-dimensional equation of heat conduction. The problems are reduced to a system of differential equations for the values of the sought approximate solution in the interior grid nodes and also to numerical formulas for the solution values at the boundary nodes. A numerical experiment is conducted. The error relative to grid spacing is established.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 60, pp. 37–43, 1986.     

Least squares moving finite elements

The moving finite element (MFE) method, when applied to purelyhyperbolic partial differential equation, moves nodes with approximatelycharacteristic speeds, which makes the method useless for steady-stateproblems. We introduce the least squares MFE method (LSMFE)for steady-state pure convection problems which corrects thisdefect. We show results for a steady-state pure convection problemin one dimension in which the nodes are no longer swept downstreamas in MFE. The method is then extended to two dimensions andthe grid aligns automatically with the flow, thereby yieldingfar greater accuracy than the corresponding fixed node leastsquares results, as is shown in two-dimensional numerical trials.     

Unsteady convective boundary layer flow of a viscous fluid at a vertical surface with variable fluid properties

In this paper we present numerical solutions to the unsteady convective boundary layer flow of a viscous fluid at a vertical stretching surface with variable transport properties and thermal radiation. Both assisting and opposing buoyant flow situations are considered. Using a similarity transformation, the governing time-dependent partial differential equations are first transformed into coupled, non-linear ordinary differential equations with variable coefficients. Numerical solutions to these equations subject to appropriate boundary conditions are obtained by a second order finite difference scheme known as the Keller-Box method. The numerical results thus obtained are analyzed for the effects of the pertinent parameters namely, the unsteady parameter, the free convection parameter, the suction/injection parameter, the Prandtl number, the thermal conductivity parameter and the thermal radiation parameter on the flow and heat transfer characteristics. It is worth mentioning that the momentum and thermal boundary layer thicknesses decrease with an increase in the unsteady parameter.     

Numerical solution of steady and unsteady flow over a vibrating airfoil in a channel

The work deals with a numerical solution of 2D steady and unsteady inviscid incompressible flow over the profile NACA 0012 in a channel. The finite volume method (FVM) in a form of cell-centered explicit schemes at quadrilateral C-mesh is used. Governing system of equations is the system of incompressible Euler equations. The method of artificial compressibility and time dependent method is applied to steady computations. The small disturbance theory (SDT) applied to a numerical solution of flow over a rotated profile by a small angle only is mentioned. Brief introduction is given to the Arbitrary (Semi) Lagrangian-Eulerian (ALE) method used for unsteady computations. Some numerical results of unsteady flow over a vibrating profile achieved by both SDT and ALE method are presented. Unsteady flow is caused by prescribed oscillations of the profile (one degree of freedom) fixed in an elastic axis. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

Streamwise oscillations of a cylinder beneath a free surface: Free surface effects on vortex formation modes

A computational study of a viscous incompressible two-fluid model with an oscillating cylinder is investigated at a Reynolds number of 200 and at a dimensionless displacement amplitude of A=0.13 and for the dimensionless forcing cylinder oscillation frequency-to-natural vortex shedding frequency ratios, f/f₀=1.5,2.5,3.5. Specifically, two-dimensional flow past a circular cylinder subject to forced in-line oscillations beneath a free surface is considered. The method is based on a finite volume discretization of the two-dimensional continuity and unsteady Navier-Stokes equations (when a solid body is present) on a fixed Cartesian grid. Two-fluid model based on improved volume-of-fluid method is used to discretize the free surface interface. The study focuses on the laminar asymmetric flow structure in the near wake region and lock-on phenomena at a Froude number of 0.2 and for the dimensionless cylinder submergence depths, h=0.25, 0.5 and 0.75. The equivorticity patterns and pressure distribution contours are used for the numerical flow visualization. The code validations in special cases show good comparisons with previous numerical results.     

A numerical method to compute solidification and melting processes

This paper deals with thermodynamically consistent numerical predictions of solidification and melting processes of pure materials using moving grids. Till date, enthalpy-porosity-based formulations of numerical codes have been generally the popular choice, although because of an artificial numerical smearing of the interface, it is virtually impossible to reproduce a sharp melting/solidification front that is supposed to exist for phase changes of pure substances. Numerical techniques based on moving grid methods have been relatively less used as they rely on complex and time-consuming adaptive grid generations. Using the moving grid approach, the authors present a method to solve solidification and melting problems. A simple linear interpolation is used to slide grid nodes along the interface to handle the otherwise obtained grid skewness near the interface. The numerical approach employed is validated with standard test cases available in the literature. The capability of capturing very complex flow field structures and the superiority of the present approach over enthalpy-porosity-based formulations is discussed. The authors also demonstrate the ability of the set-up computer code to solve complex thermofluid processes such as occur during crystal growth in Czochralski reactors.     

Error estimate of a fully discrete defect correction finite element method for unsteady incompressible Magnetohydrodynamics equations

In this study, a fully discrete defect correction finite element method for the unsteady incompressible Magnetohydrodynamics (MHD) equations, which is leaded by combining the Back Euler time discretization with the two-step defect correction in space, is presented. It is a continuous work of our formal paper [Math Method Appl Sci. 2017. DOI:10.1002/mma.4296]. The defect correction method is an iterative improvement technique for increasing the accuracy of a numerical solution without applying a grid refinement. Firstly, the nonlinear MHD equation is solved with an artificial viscosity term. Then, the numerical solutions are improved on the same grid by a linearized defect-correction technique. Then, we introduce the numerical analysis including stability analysis and error analysis. The numerical analysis proves that our method is stable and has an optimal convergence rate. Some numerical results [see Math Method Appl Sci. 2017. DOI:10.1002/mma.4296] show that this method is highly efficient for the unsteady incompressible MHD problems.     

The stream function coordinate system for solution of one-dimensional unsteady compressible flow problems

A one-dimensional unsteady compressible isentropic flow problem is solved using a floating grid finite difference method. GENMIX, a computer program for solution of two-dimensional parabolic flow problems, has been adopted for this purpose. The advantages of utilizing the staggered floating grid method are demonstrated through solution of flow in a shocktube. The present method is able to precisely locate the discontinuities in temperature and density profiles.     

Stability and convergence of two‐grid Crank‐Nicolson extrapolation scheme for the time‐dependent natural convection equations

In this paper, we consider the Crank‐Nicolson extrapolation scheme for the 2D/3D unsteady natural convection problem. Our numerical scheme includes the implicit Crank‐Nicolson scheme for linear terms and the recursive linear method for nonlinear terms. Standard Galerkin finite element method is used to approximate the spatial discretization. Stability and optimal error estimates are provided for the numerical solutions. Furthermore, a fully discrete two‐grid Crank‐Nicolson extrapolation scheme is developed, the corresponding stability and convergence results are derived for the approximate solutions. Comparison from aspects of the theoretical results and computational efficiency, the two‐grid Crank‐Nicolson extrapolation scheme has the same order as the one grid method for velocity and temperature in H¹‐norm and for pressure in L²‐norm. However, the two‐grid scheme involves much less work than one grid method. Finally, some numerical examples are provided to verify the established theoretical results and illustrate the performances of the developed numerical schemes.     

Capabilities of a quasi-gasdynamic algorithm as applied to inviscid gas flow simulation

Ten well-known one-dimensional test problems reflecting the characteristic features of unsteady inviscid gas flows are successfully solved by a unified numerical algorithm based on the quasigasdynamic system of equations. In all the cases, the numerical solution converges to a self-similar one as the spatial grid is refined.     

Numerical simulation of generalized newtonian blood flow past a couple of irregular arterial stenoses

This paper looked at the numerical investigations of the generalized Newtonian blood flow through a couple of irregular arterial stenoses. The flow is treated to be axisymmetric, with an outline of the stenoses obtained from a three dimensional casting of a mild stenosed artery, so that the flow effectively becomes two‐dimensional. The Marker and Cell (MAC) method is developed for the governing unsteady generalized Newtonian equations in staggered grid for viscous incompressible flow in the cylindrical polar co‐ordinates system. The derived pressure‐Poisson equation was solved using Successive‐Over‐Relaxation (S.O.R.) method and the pressure‐velocity correction formulae have been derived. Computations are performed for the pressure drop, the wall shear stress distribution and the separation region. The presented computations show that in comparison to the corresponding Newtonian model the generalized Newtonian fluid experiences higher pressure drop, lower peak wall shear stress and smaller separation region. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 960–981, 2011