Sharif Digital Repository

In this study, two finite-difference lattice Boltzmann methods (FDLBM) are applied and assessed for the simulation of two-phase liquid-vapor flows with high density ratios. For this aim, the He-Shan-Doolen type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting system of equations are discretized by using the second-order central difference and modified Lax-Wendroff schemes. Suitable numerical dissipation terms and filters are applied to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions and at the same time to remove the numerical oscillations in the interface region of the two phases.Three... 

In the present study, a high-order finite-difference lattice Boltzmann solver is applied for simulating steady and unsteady three-dimensional incompressible flows. To achieve an accurate and robust flow solver, the incompressible form of the lattice Boltzmann equation in the three-dimensional generalized curvilinear coordinates is discretized spatially based on the fifth-order weighted essentially non-oscillatory (WENO) finite-difference scheme. To ensure the stability and temporal accuracy of the flow solver, the fourth-order Runge-Kutta method is used for the time integration. To examine the accuracy and performance of the flow solver, different three-dimensional incompressible flow... 

Progress aerodynamic and structural applications has lead to solid-fluid interaction (SFI) engineering. One major application of SFI is to improve of aerodynamic characteristics of airplanes. One option is to use elastic shells.this field of study is so novel and a challenge for future researches. This what current thesis is focused on We use numerical simulation in both fluid and solid parts.In other word, we employ the finite –element method to solve the solid mechanics governing equation and the finite volume element method to treat the fluid dynamics governing equations. The fluid governing equations are 2D navier- stokes equations and the solid part is the 1D Euler – Bernoulli equations... 

In the present study, the simulation of incompressible Electro-Magneto-hydrodynamic flows is performed using a finite volume lattice Boltzmann method (FVLBM). The Boltzmann transport equation is solved using a cell-centered finite volume method on structured meshes. A central difference scheme is used to discretize the spatial derivatives and the fourth-order numerical dissipation term is added to stabilize the solution. To discretize the temporal derivative, the fourth-order Runge-Kutta time stepping scheme is applied. The standard collision-streaming lattice Boltzmann method has been used to simulate EMHD flows in the literature, however, it has several deficiencies such as the... 

In this study, a high-order finite-difference lattice Boltzmann method (FDLBM) is used to simulate the two-dimensional incompressible flows. Here, the incompressible form of the lattice Boltzmann (LB) equation in the two-dimensional generalized curvilinear coordinates is considered and the resulting equation is discretized based on both the third- and fifth-order upwind finite-difference schemes. The time integration of the present flow solver is performed by the fourth-order Runge-Kutta method. Several incompressible laminar flow problems are simulated to examine the accuracy and performance of the developed high-order FDLBM solver. The present results are compared with the existing... 

In this study, fluid-solid interaction (FSI) is simulated computationally by using a high-order accurate numerical method. The two-dimensional incompressible viscous flows are considered in the fluid domain. The primary problem with solutions of the incompressible Navier–Stokes equations is the difficulty of coupling changes in the velocity field with changes in the pressure field while satisfying the continuity equation. Herein, the artificial compressibility method is used to overcome this difficulty. Preconditioning is implemented to reduce the stiffness of the system of equations to increase the convergence rate of the solution. Using preconditioning, physical solutions even at low... 

In the present study, the preconditioned incompressible Navier‐Stokes equations with the artificial compressibility (AC) 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. A fourth‐order compact finite‐difference scheme is utilized to discretize the spatial derivative terms of the resulting system of equations and the time integration is carried out based on the dual time‐stepping method. The capability of the proposed solution methodology for computing the steady and unsteady incompressible viscous flows in a wide range of Reynolds... 

In the present study, an incompressible smoothed particle hydrodynamics method based on vorticity-stream function (VSF-SPH) formulation is developed and assessed for simulating steady and unsteady incompressible flows. The vorticity-stream function formulation in the Eulerian reference frame is written in a Lagrangian reference frame to provide an appropriate incompressible SPH algorithm. The advantage of the proposed smoothed particle hydrodynamics method based on the vorticity-stream function (VSF-SPH) formulation over the weakly compressible SPH (WCSPH) is that the VSF-SPH method is a truly incompressible SPH algorithm and it does not involve any approximate enforcement of the... 

In the present thesis, a high-order compact finite-difference lattice Boltzmann method (CFDLBM) is proposed and applied for an accurate and efficient numerical simulation of liquid-vapor two-phase flows. At first, the stability of the fourth-order CFDLBM is performed by using the von Neumann stability analysis for the D2Q7 and D2Q9 lattices. The stability analysis indicates that the CFDLBM proposed is stable and thus suitable for the simulation of high Reynolds number flows. The high-order CFDLBM is then developed and applied to accurately compute 2-D and 3-D incompressible flows in the Cartesian coordinates. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized... 

In the present study, the numerical simulation of incompressible laminar and turbulent flows using a high-order finite difference lattice Boltzmann method is presented. To handle curved geometries with non uniform grids, the incompressible form of lattice Boltzmann equation is transformed into the generalized curvilinear coordinates and the spatial derivatives of the resulting equation are discretized using the fifth-order WENO scheme. The advantage of using the WENO-LBM developed is that it needs less number of grid points and remains stable even at high Reynolds number flows. For the temporal term, the fourth-order explicit Rung-Kutta scheme is adopted for laminar flow calculations and... 

In this study, the numerical simulation of the two-dimensional incompressible flow around an oscillating airfoil is performed. For this aim, the incompressible Navier-Stokes equations based on the artificial compressibility approach written in the arbitrary Lagrangian-Eulerian form are considered. Then, the vorticity confinement method is incorporated in the formulation and the resulting system of equations is solved by a second-order central-difference finite volume method with the controllable dissipation terms. For the time integration, the implicit dual-time stepping scheme is implemented. At first, the numerical solution of the incompressible flow over the oscillating NACA0012 airfoil... 

The purpose of this study is to investigate the distribution of air flow and temperature in electric machines. The fluid flow equations are solved by using the characteristic Galerkin. The unsteady equations of the air flow and heat transfer between the rotor, stator and radial cooling channels of the generator are solved by finite element method. Tetrahedral elements are used, but hexahedral elements can also be used. The upwind method is used to discretize the convection terms, but the central difference for the diffusion parts, that is equal to Petrov-Galerkin method in steady state situation. Due to the large number of channels and specific geometry, to reduce the computation costs, the... 

In this research, an anistropic thin-walled beam composite is used to consider the effects of fiber orientation and lay-up configuration on the aeroelastic stability and response of an aircraft wing. The circum ferentially asymmetric stiffness (CAS) model consists of a group of non-classical effects; such as transverse shear, material anisotropic, warping inhibition, etc. The aerodynamic load has been modeled based on the strip theory and compressible unsteady flow in the finite- state form. In addition, to form the mass, stiffness and damping matrices of non-conservative aeroelastic systems the extended Galerkin’s method (EGM) and the method of separation of variables has been used. After...