Sharif Digital Repository

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 study, the nonlinear partial differential equations governing two phase flow through porous media are solved using two different methods, namely, finite difference and finite element. The capillary pressure term is considered in the mathematical model. The numerical results on a 2-D test case are then compared with the experimental drainage process and water flooding performed on a glass type micromodel. Based on the obtained results, finite difference technique needs less computational time for solving governing equations of two phase flow, but findings of this method show less agreement with the experimental data. The finite element scheme was found to be more adequate and its... 

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

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

Flow movement in unsaturated soil can be expressed by a partial differential equation, named Richards equation. The objective of this study is the finding of an appropriate implicit numerical solution for head based Richards equation. Some of the well known finite difference schemes (fully implicit, Crank Nicolson and Runge-Kutta) have been utilized in this study. In addition, the effects of different approximations of moisture capacity function, convergence criteria and time stepping methods were evaluated. Two different infiltration problems were solved to investigate the performance of different schemes. These problems include of vertical water flow in a wet and very dry soils. The... 

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

Mesh point positions in Finite Difference Method (FDM) of discretization for the neutron diffusion equation can remarkably affect the averaged neutron fluxes as well as the effective multiplication factor. In this study, by aid of improving the mesh point positions, an enhanced finite difference scheme for the neutron diffusion equation is proposed based on the neutron importance function. In order to determine the neutron importance function, the adjoint (backward) neutron diffusion calculations are performed in the same procedure as for the forward calculations. Considering the neutron importance function, the mesh points can be improved through the entire reactor core. Accordingly, in... 

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

This paper uses a fourth-order compact finite-difference scheme for solving steady incompressible flows. The high-order compact method applied is an alternating direction implicit operator scheme, which has been used by Ekaterinaris for computing two-dimensional compressible flows. Herein, this numerical scheme is efficiently implemented to solve the incompressible Navier-Stokes equations in the primitive variables formulation using the artificial compressibility method. For space discretizing the convective fluxes, fourth-order centered spatial accuracy of the implicit operators is efficiently obtained by performing compact space differentiation in which the method uses block-tridiagonal... 

The main purpose of this paper is modeling and simulation of in-situ releasing of smart nano-sized core-shell particles at the water-oil interface during polymer flooding. During the polymer flooding process, when these nano-particles reach the water-oil interface, migrate to the oil phase and the hydrophobic layer of them dissolves in this phase. After dissolution of this protective nano-sized layer, the hydrophilic core containing a water-soluble ultra high molecular weight polymer diffuses back into the water phase and with dissolving in this phase, dramatically increases viscosity of flooding water in the neighborhood of the water-oil interface. In this study, two different... 

The current study investigates the conjugate heat transfer characteristics for laminar flow in backward facing step channel. All of the channel walls are insulated except the lower thick wall under a constant temperature. The upper wall includes a insulated obstacle perpendicular to flow direction. The effect of obstacle height and location on the fluid flow and heat transfer are numerically explored for the Reynolds number in the range of 10 ≤ Re ≤ 300. Incompressible Navier-Stokes and thermal energy equations are solved simultaneously in fluid region by the upwind compact finite difference scheme based on flux-difference splitting in conjunction with artificial compressibility method. In... 

The current study investigates the conjugate heat transfer characteristics for laminar flow in backward facing step channel. All of the channel walls are insulated except the lower thick wall under a constant temperature. The upper wall includes a insulated obstacle perpendicular to flow direction. The effect of obstacle height and location on the fluid flow and heat transfer are numerically explored for the Reynolds number in the range of 10 ≤ Re ≤ 300. Incompressible Navier-Stokes and thermal energy equations are solved simultaneously in fluid region by the upwind compact finite difference scheme based on flux-difference splitting in conjunction with artificial compressibility method. In... 

This paper has developed an axisymmetric laminar and turbulent two-phase flow solver to simulate pressure-swirl atomizers. Equations include the explicit algebraic Reynolds stress model, the Reynolds-averaged Navier-Stokes, and the level set equation. Applying a high-order compact upwind finite difference scheme with the level set equation being culminated to capture the interface between air-liquid two-phase flow and decreasing the mass conservation error in the level set equation. The results show that some recirculation zones are observed close to the wall in the swirl chamber and to the axis. This model can predict converting the Rankin vortex in the swirl chamber to the forced vortex in... 

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

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

This study reports on the simulation of temperature distribution of human tooth under a laser beam based on non-Fourier models. The temperature in the tooth depth that directly results from the conduction heat transfer process is caused by the lengthy thermal relaxation time in the tooth layers. A detailed tooth composed of enamel, dentin, and pulp with unstructured shape, uneven boundaries, and realistic thicknesses was considered. A finite difference scheme was separately adopted to solve time-dependent equations in solid layers and soft tissue of the tooth. In this study, a dual-phase-lag non-Fourier heat conduction model was applied to evaluate temperature distribution induced by laser... 

The main goal of this study is to assess the application of high-order compact finite-difference schemes for the solution of the Euler equations in conjunction with the compressible vorticity confinement method on both uniform Cartesian and curvilinear grids. Here, the spatial discretization of the governing equations is performed by the fourth-order compact finite-difference scheme and the temporal term is discretized by the fourth-order Runge-Kutta method. To stabilize the numerical solution, appropriate dissipation terms are applied and a detail assessment is performed to study the effects of the values of confinement and dissipation coefficients on the solution to reasonably preserve the... 

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 this work, a preconditioned high-order weighted essentially non-oscillatory (WENO) finite-difference lattice Boltzmann method (WENO-LBM) is applied to deal with the incompressible turbulent flows. Two different turbulence models namely, the Spalart–Allmaras (SA) and k−ωSST models are used and applied in the solution method for this aim. The spatial derivatives of the two-dimensional (2D) preconditioned LB equation in the generalized curvilinear coordinates are discretized by using the fifth-order WENO finite-difference scheme and an implicit–explicit Runge–Kutta scheme is adopted for the time discretization. For the convective and diffusive terms of the turbulence transport equations, the...