Sharif Digital Repository

In this work, the implementation of a high-order compact finite-difference lattice Boltzmann method (CFDLBM) is performed in the generalized curvilinear coordinates to improve the computational efficiency of the solution algorithm to handle curved geometries with non-uniform grids. The incompressible form of the discrete Boltzmann equation with the Bhatnagar-Gross-Krook (BGK) approximation with the pressure as the independent dynamic variable is transformed into the generalized curvilinear coordinates. Herein, the spatial derivatives in the resulting lattice Boltzmann (LB) equation in the computational plane are discretized by using the fourth-order compact finite-difference scheme and the... 

Despite significant progress in unstructured grid generation and employment, the robust employment of body-conforming coordinate systems promotes the workers to consider it as a major alternative to treat relatively complex flow fields in irregular geometries. Contrary to the laminar flow treatment, there are many ambiguities around treating turbulent flow on body-fitted coordinate system. The ambiguities are mainly originated from the curvature on the boundaries. Considering the past taken efforts in improving the results of treating turbulent flow on curvilinear coordinate system, we have extended a new formulation on body-fitted coordinate system using physical covariant velocities as the... 

In this paper, Mindlin's second strain gradient theory is formulated and presented in an arbitrary orthogonal curvilinear coordinate system. Equilibrium equations, generalized stress-strain constitutive relations, components of the strain tensor and their first and second gradients, and the expressions for three different types of traction boundary conditions are derived in any orthogonal curvilinear coordinate system. Subsequently, for demonstration, Mindlin's second strain gradient theory is represented in the spherical coordinate system as a highly-practical coordinate system in nanomechanics. Second strain gradient elasticity have been developed mainly for its ability to capture the... 

In the present study, the preconditioned incompressible Navier-Stokes equations with the artificial compressibility method formulated in the generalized curvilinear coordinates are numerically solved by using a high-order compact finite-difference scheme for accurately and efficiently computing the incompressible flows in a wide range of Reynolds numbers. A fourth-order compact finite-difference scheme is utilized to accurately discretize the spatial derivative terms of the governing equations, and the time integration is carried out based on the dual time-stepping method. The capability of the proposed solution methodology for the computations of the steady and unsteady incompressible... 

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to... 

The characteristic boundary conditions are applied and assessed for the solution of incompressible inviscid flows. The two-dimensional incompressible Euler equations based on the artificial compressibility method are considered and then the characteristic boundary conditions are formulated in the generalized curvilinear coordinates and implemented on both the far-field and wall boundaries. A fourth-order compact finite-difference scheme is used to discretize the resulting system of equations. The solution methodology adopted is more suitable for this assessment because the Euler equations and the high-accurate numerical scheme applied are quite sensitive to the treatment of boundary... 

The buoyancy-induced turbulent flows in complex geometries with internal heat source have important application in enclosure fire dynamics, for example, smoke movement in the situation of fire in tunnels. The majority of computational works studying the fire-induced hot air movement in such cases use a Cartesian coordinate system. Hence, to simulate curved geometries, the Cartesian grid should be refined so that it comparatively matches the geometry. In the present work, a three-dimensional numerical method using nonorthogonal curvilinear coordinate system was developed. Numerical simulations were performed to study the flow structures and the heat transfer characteristics of fire-induced... 

The unsteady preconditioned characteristic boundary conditions (UPCBCs) based on the artificial compressibility (AC) method are formulated and applied at artificial boundaries for the direct numerical simulation (DNS) of incompressible flows. The compatibility equations including the unsteady terms are mathematically derived in the generalized curvilinear coordinates and then incorporated as boundary conditions (BCs) in a high-order accurate incompressible flowsolver. The spatial derivative terms of the systemof equations are discretized using the fourth-order compact finite difference (FD) scheme, consistent with the high-order accuracy required for the DNS. The time integration is carried... 

The numerical solution of the parabolized Navier-Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth-order compact finite-difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock fitting procedure is utilized in both the compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are... 

The present study focuses on the numerical simulation of cavitation around the NACA 0015. The unsteady behaviors of cavitation which have worthwhile applications are investigated. The cavitation patterns, velocity fields and frequency of the cavitating flow around hydrofoil is obtained. For multi phase simulation, single-fluid Navier-Stokes equations, along with a volume fraction transport equation, are employed. The bubble dynamics model is utilized to simulate phase change. SIMPLE algorithm is used for velocity and pressure computations. For discretization of equations the finite-volume approach written in body fitted curvilinear coordinates, on collocated grid, is used. In this study,... 

A high-order compact finite-difference scheme is applied and assessed for the numerical simulation of structural dynamics. The two-dimensional elastic stress-strain equations are considered in the generalized curvilinear coordinates and the spatial derivatives in the resulting equations are discretized by a fourth-order compact finite-difference scheme. For the time integration, an implicit second-order dual time-stepping method is utilized in which a fourth-order Runge-Kutta scheme is used to integrate in the pseudo-time level. The accuracy and robustness of the solution procedure proposed are investigated through simulating different two-dimensional benchmark test cases in structural... 

An experimental and numerical study has been carried out to investigate the heat transfer characteristics of a horizontal circular cylinder exposed to a slot jet impingement of air. A square-edged nozzle is mounted parallel with the cylinder axis and jet flow impinges on the bottom of the cylinder. The study is focused on low Reynolds numbers ranging from 120 to 1,210, Grashof numbers up to Gr = 10Re2 and slot-to-cylinder spacing from 2 to 8 of the slot width. The flow field is greatly influenced by the slot exit velocity and the buoyancy force due to density change. A Mach-Zehnder Interferometer is used for measurement of local Nusselt number around the cylinder at 10° interval. It is...