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... 

A new shock-detecting sensor for properly switching between a second-order and a higher-order filter is developed and assessed. The sensor is designed based on an order analysis. The nonlinear filter with the proposed sensor ensures damping of the high-frequency waves in smooth regions and at the same time removes the Gibbs oscillations around the discontinuities when using high-order compact finite difference schemes. In addition, a suitable scaling is proposed to have dissipation proportional to the shock strength and also to minimize the effects of the second-order filter on the very small scales. Several numerical experiments are carried out and the accuracy of the nonlinear filter with... 

The preconditioned characteristic boundary conditions based on the artificial compressibility (AC) method are implemented at artificial boundaries for the solution of two- and three-dimensional incompressible viscous flows in the generalized curvilinear coordinates. The compatibility equations and the corresponding characteristic variables (or the Riemann invariants) are mathematically derived and then applied as suitable boundary conditions in a high-order accurate incompressible flow solver. The spatial discretization of the resulting system of equations is carried out by the fourth-order compact finite-difference (FD) scheme. In the preconditioning applied here, the value of AC parameter... 

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is extended and applied to accurately simulate two-phase liquid-vapor flows with high density ratios. Herein, the He-Shan-Doolen-type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting equations are discretized by using the fourth-order compact finite-difference scheme and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient two-phase flow solver. A high-order spectral-type low-pass compact nonlinear filter is used to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions... 

A high-order compact finite-difference total Lagrangian method (CFDTLM) is developed and applied to nonlinear structural dynamic analysis. The two-dimensional simulation of thermo-elastodynamics is numerically performed in generalized curvilinear coordinates by taking into account the geometric and material nonlinearities. The spatial discretization is carried out by a fourth-order compact finite-difference scheme and an implicit second-order accurate dual time-stepping method is applied for the time integration. The accuracy and capability of the proposed solution methodology for the nonlinear structural analysis is investigated through simulating different static and dynamic benchmark... 

This paper investigates implementation of upwind compact implicit and explicit high-order finite difference schemes for solution of the level set equation. The upwind compact implicit and explicit high-order finite difference schemes are well-known techniques to descritize spatial derivatives for convection term in hyperbolic equations. Applying of upwind high-order schemes on the level set equation leads to less error and CPU time reduction compared to essential non-oscillatory (ENO), weighted essential non-oscillatory schemes (WENO), and even different particle level set methods. The results indicate the error based on area loss decreases drastically with applying high-order upwind,... 

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... 

In this work, the capability and performance of the vorticity confinement (VC) implemented in a high-order accurate flow solver in predicting two-dimensional (2D) compressible mixing layer flows on coarse grids are investigated. Here, the system of governing equations with incorporation of the VC in the formulation is numerically solved by the fourth-order compact finite difference scheme. To stabilize the numerical solution, a low-pass high-order filter is applied, and the nonreflective boundary conditions are used at the farfield and outflow boundaries to minimize the reflections. At first, the numerical results without applying the VC are validated by available direct numerical... 

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... 

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 objective of this work is to numerically study the fluid flow and acoustic field of a supersonic impinging jet by applying the vorticity confinement (VC) method. For this aim, the three-dimensional compressible Navier-Stokes equations with the incorporation of the VC method are considered and the resulting system of equations is solved by using the sixth-order compact finite-difference scheme. To eliminate the numerical instability, a low-pass high-order filter is used. The nonreflective boundary conditions are applied for all the free boundaries and the radiated sound field is obtained by the Kirchhoff surface integration. Comparisons of the present results with the experimental data... 

The objective of this work is to numerically study the fluid flow and acoustic field of a supersonic impinging jet by applying the vorticity confinement (VC) method. For this aim, the three-dimensional compressible Navier-Stokes equations with the incorporation of the VC method are considered and the resulting system of equations is solved by using the sixth-order compact finite-difference scheme. To eliminate the numerical instability, a low-pass high-order filter is used. The nonreflective boundary conditions are applied for all the free boundaries and the radiated sound field is obtained by the Kirchhoff surface integration. Comparisons of the present results with the experimental data...