Sharif Digital Repository

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

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

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, two preconditioners proposed by Eriksson, and Choi and Merkel are implemented on a 3D upwind Euler flow solver on unstructured meshes. The mathematical formulations of these preconditioning schemes for the set of primitive variables Q→p 1=[ρ u v , w p] T are drawn and their eigenvalues and eigenvectors are compared with each others. A cell-centered finite volume Roe's method is used for discretization of the 3D preconditioned Euler equations. The accuracy and performance of these preconditioning schemes are examined by computing low Mach number flows over the ONERA M6 wing for different conditions  

A numerical treatment for the prediction of cavitating flows is presented and assessed. The algorithm uses the preconditioned multiphase Euler equations with appropriate mass transfer terms. A central difference finite volume scheme with suitable dissipation terms to account for density jumps across the cavity interface is shown to yield an effective method for solving the multiphase Euler equations. The Euler equations are utilized herein for the cavitation modeling, because some certain characteristics of cavitating flows can be obtained using the solution of this system of equations with relative low computational effort. In addition, the Euler equations are appropriate for the assessment... 

An improved progressive preconditioning method for analyzing steady inviscid and laminar flows around fully wetted and sheet-cavitating hydrofoils is presented. The preconditioning matrix is adapted automatically from the pressure and/or velocity flow-field by a power-law relation. The cavitating calculations are based on a single fluid approach. In this approach, the liquid/vapour mixture is treated as a homogeneous fluid whose density is controlled by a barotropic state law. This physical model is integrated with a numerical resolution derived from the cell-centered Jameson's finite volume algorithm. The stabilization is achieved via the second-and fourth-order artificial dissipation... 

A new method of spatial filtering in high-order finite volume methods is presented and assessed. The base of this method is to filter face-averaged variables (fluxes) and then the recovery of cell-averaged ones. Two kinds of filtering method are proposed. The first kind is highly dissipative and appropriate for the numerical regions that need high dissipation, e.g. sponge zones. The second kind, on the other hand, is a precise method and hence is suitable for applying the high-order finite difference filters to the finite volume methods. Applying high-order finite difference filters directly to the high-order finite volume methods without using the proposed method causes stability problems... 

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

In the present work, two cavitation modeling strategies, namely the barotropic cavitation model and the transport equation-based model are applied and assessed for the numerical simulation of inviscid cavitating flows over two-dimensional and axisymmetric geometries. The algorithm uses the preconditioned Euler equations employing the interface capturing method for both the cavitation models. A same numerical solution procedure is used herein for discretizing the governing equations resulting from these two cavitation modeling strategies for the assessment to be valid and reliable. A central difference finite-volume scheme employing the suitable dissipation terms to account for density jumps... 

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth-order compact FD scheme, and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient incompressible flow solver. A high-order spectral-type low-pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also... 

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

In this study, the Chebyshev collocation spectral lattice Boltzmann method (CCSLBM) is developed and assessed for the computation of low-speed flows. Both steady and unsteady flows are considered here. The discrete Boltzmann equation with the Bhatnagar-Gross-Krook approximation based on the pressure distribution function is considered and the space discretization is performed by the Chebyshev collocation spectral method to achieve a highly accurate flow solver. To provide accurate unsteady solutions, the time integration of the temporal term in the lattice Boltzmann equation is made by the fourth-order Runge-Kutta scheme. To achieve numerical stability and accuracy, physical boundary... 

In this work, a high-order nodal discontinuous Galerkin method is applied and assessed for the simulation of compressible non-cavitating and cavitating flows. The one-fluid approach with the thermal effects is used to properly model the cavitation phenomenon. Here, the spatial and temporal derivatives in the system of governing equations are discretized using the nodal discontinuous Galerkin method and the third-order TVD Runge–Kutta method, respectively. Various numerical fluxes such as the Roe, Rusanov, HLL, HLLC and AUSM+-up and two discontinuity capturing methods, namely, the generalized MUSCL limiter and a generalized exponential filter are implemented in the solution algorithm. At... 

In this work, the Chebyshev collocation spectral lattice Boltzmann method is implemented in the generalized curvilinear coordinates to provide an accurate and efficient low-speed LB-based flow solver to be capable of handling curved geometries with non-uniform grids. The low-speed form of the D2Q9 and D3Q19 lattice Boltzmann equations is transformed into the generalized curvilinear coordinates and then the spatial derivatives in the resulting equations are discretized by using the Chebyshev collocation spectral method and the temporal term is discretized with the fourth-order Runge–Kutta scheme to provide an accurate and efficient low-speed flow solver. All boundary conditions are... 

In the present study, shock-disturbances interaction in high-speed compressible inviscid flow is simulated utilizing the weighted essentially non-oscillatory (WENO) scheme by employing the shock-capturing technique. For this aim, the two-dimensional Euler equations in strong conservative form are discretized by using the explicit third-order TVD Runge–Kutta scheme in time and the fifth-order WENO finite difference scheme in space. The main advantage of using the WENO scheme is its capability for properly solving the discontinuities in the domain without needing any artificial viscosity, limiter function or filter. Hence, this scheme is stable, and thus, it is suitable for simulating very... 

In this study, a theoretical analysis is performed to assess the interaction of freestream disturbances with a plane normal shock considering real gas effects. Such effects are important in a field with high velocities and high temperatures. To perform the theoretical analysis, the downstream disturbances field is expressed as a mathematical function of the upstream one by incorporating real gas effects in the formulation. Here, the linearized one-dimensional perturbed unsteady Euler equations are used for the classification of the downstream/upstream disturbances field and the linearized one-dimensional perturbed Rankine-Hugoniot equations are applied to provide a relationship between the... 

In the present work, the spectral difference lattice Boltzmann method (SDLBM) is implemented on unstructured meshes for the solution methodology to be capable of accurately simulating the compressible flows over complex geometries. Both the inviscid and viscous compressible flows are computed by applying the unstructured SDLBM. The compressible form of the discrete Boltzmann–BGK equation with the Watari model is considered and the solution of the resulting system of equations is obtained by applying the spectral difference method on arbitrary quadrilateral meshes. The accuracy and robustness of the unstructured SDLBM for simulating the compressible flows are demonstrated by simulating four... 

In this study, a theoretical analysis is performed to assess the interaction of freestream disturbances with a plane normal shock considering real gas effects. Such effects are important in a field with high velocities and high temperatures. To perform the theoretical analysis, the downstream disturbances field is expressed as a mathematical function of the upstream one by incorporating real gas effects in the formulation. Here, the linearized one-dimensional perturbed unsteady Euler equations are used for the classification of the downstream/upstream disturbances field and the linearized one-dimensional perturbed Rankine-Hugoniot equations are applied to provide a relationship between the...