Sharif Digital Repository

The growth of slant cracks by fretting fatigue of a half plane in contact with a flat rounded pad was studied. The mode I and mode II stress intensity factors for cracks of various lengths and directions were calculated using the semi-analytical method of the distribution of dislocations, and their cumulative effect on the crack growth was investigated using the strain energy density criterion. The results showed dominance of mode I fracture on crack growth, and based on the observation of strain energy density factor versus crack orientation, the results also suggest that depending on the crack length, the most critical orientations of cracks are at 0° to 20° outward the contact zone. Good... 

The lateral movement of a liquefiable soil layer on gentle slopes is the most visible and devastating type of liquefaction-induced ground failure. Recent earthquakes have shown that this phenomenon causes severe damages to coastal structures, pier of the bridges and life-lines by exerting large lateral forces on the structures. In this paper coupled dynamic field equations of extended Biot's theory with u-p formulation are used for simulating the phenomenon and the soil behavior is modeled by a critical state two-surface plasticity model for sands. Furthermore, in this study variation of permeability coefficient during liquefaction is taken into account. The permeability coefficient is... 

Thermal cracking is the prevalent type of distress experienced by asphalt pavements in cold regions. It is widely assumed that when thermal stresses induced in the pavement exceed the tensile strength of the asphalt surface layer, cracking occurs; however, the role of thermal fatigue should not be ignored. To better describe the low-temperature (LT) performance properties of modified asphalt binders, a new parameter, normalised tensile stress (NTS), was defined in this research. NTS values were compared with bending beam rheometer (BBR), direct tension (DT) and semi-circular bending (SCB) fracture test results. The m-value parameter of the BBR test specifies the type of asphalt binder. This... 

In this work, numerical investigation of a gas turbine model combustor (GTMC) was carried out using two different turbulence-chemistry interaction models: the EDC (eddy dissipation concept) and TPDF (transported probability density function). GTMC with good optical access for laser measurements provides a useful database for swirling CH4/Air diffusion flames at atmospheric pressure. Modeling was performed by solving Reynolds-averaged Navier-Stokes (RANS) and Reynolds stress model (RSM) equations for a two-dimensional (2D) axisymmetric computational domain accompanied by swirl and the combustion chamber was investigated for both reacting and nonreacting conditions. A detailed reduced... 

Recently, the mechanisms of natural undulatory locomotion of aquatic animal swimming have become one of the most significant issues for the researchers, swimmers and engineers. This study aims to elucidate and compare the propulsive vortical signature and performance of backward (negative undulation) and forward (positive undulation) travelling waves through a typical fishlike propulsor by a systematic numerical study. The numerical approach uses a pressure-based finite volume method solver to solve Navier–Stokes equations in an arbitrary Lagrangian–Eulerian framework domain containing a two-dimensional NACA 0012 foil moving with prescribed kinematics. Some of the important findings are: (1)... 

In this paper, a one-dimensional numerical approach is used to study the effect of various parameters such as micro combustor diameter, mass flow rate and external convection heat transfer coefficient on the temperature and species mass fraction profiles. A premixed mixture of H2-Air with a multi-step chemistry (9 species and 19 reactions) is used and thermal conductivity of the mixture is considered as a function of species thermal conductivity and temperature by using a set of new relations. The transient gas phase energy and species conservation equations result in an Advection-Diffusion-Reaction system (A-D-R) that leads to two stiff systems of PDEs, which can not be solved by... 

Thermal cracking is the prevalent type of distress experienced by asphalt pavements in cold regions. It is widely assumed that when thermal stresses induced in the pavement exceed the tensile strength of the asphalt surface layer, cracking occurs; however, the role of thermal fatigue should not be ignored. To better describe the low-temperature (LT) performance properties of modified asphalt binders, a new parameter, normalised tensile stress (NTS), was defined in this research. NTS values were compared with bending beam rheometer (BBR), direct tension (DT) and semi-circular bending (SCB) fracture test results. The m-value parameter of the BBR test specifies the type of asphalt binder. This... 

This article is devoted to the numerical investigation of the shock wave/boundary layer interaction (SWBLI) as the main factor influencing the aerodynamic performance of transonic bumped airfoils and wings. The numerical analysis is conducted for the ONERA-M6 wing through a shock control bump (SCB) shape optimisation process using the adjoint optimisation method. SWBLI is analyzed for both clean and bumped airfoils and wings, and it is shown how the modified wave structure originating from upstream of the SCB reduces the wave drag, by improving the boundary layer velocity profile downstream of the shock wave. The numerical simulation of the turbulent viscous flow and a gradient-based adjoint... 

The present work is a numerical investigation into the waves generated by a pressure source moving in straight channels with a nonrectangular cross-channel depth profile. Wave fields generated by the moving pressure source are described, and the effects of channel bathymetry on the generated wave characteristics of wave height, wave-breaking location, wave-breaking type and intensity, and peel angle are presented. Four different channel cross-section shapes were investigated, and the results were analyzed. It was concluded that the maximum wave height depended primarily on pressure-source parameters and that the investigated channel shapes did not have a significant effect on the maximum... 

The aim of the present work is to assess the effects of wettability, heterogeneity, and viscosity differences on water-oil displacement process in micromodel porous media through numerical modelling. The two-phase flow was simulated by Cahn-Hilliard phase field method (PFM) using a finite element package. The micromodel was initially saturated with oil (wetting phase) and oil was produced through invasion of the displacing phase into the matrix. The computed oil and water saturations were in good agreement with those obtained by the visual flooding experiment. Using the validated model, sensitivity analysis was performed to investigate the effects of different wettability states,... 

This paper addresses numerically-derived tsunami wave loads on bridge superstructures using smoothed particle hydrodynamics (SPH), which is a type of mesh-free methods. Although there exist some relationships for the case of impinged loads on bridges exerted by regular (sinusoidal) waves, for the case of solitary waves such as tsunamis, no relation has yet been proposed in the literature. This shortcoming is partly due to the lack of understanding the mechanism of wave action on the bridge superstructures. In this study, three water depths, three wave amplitudes and four submergence depths of the deck are considered for the process of numerical investigation of tsunami-induced loads on... 

In this article, the behavior of gas bubbles in tapered fluidized beds is investigated with the use of a two-fluid model incorporating kinetic theory of granular flow. The effects of various parameters such as apex angle, particle size, and particle density on the size distribution and the rise velocity of gas bubbles were examined. In addition, the simulation results for the bubble fraction and axial velocity of gas bubbles were compared with experimental data reported in the literature and good agreement was observed. As the apex angle was increased, the fraction of gas bubbles with large sizes increased and the fraction of bubbles with small sizes decreased. As the particle size... 

In this article, ferrohydrodynamic forced-convection heat transfer from a heated sphere embedded in a ferrofluid in the presence of the uniform external magnetic field has been studied numerically for the first time over a wide range of Reynolds number value, nanoparticle diameter, particle volume fraction, and magnetic field intensity. Despite the uniform external magnetic field applied, the internal magnetic field near the sphere could be nonuniform due to the considerable difference between the relative magnetic permeability of the sphere and the surrounding medium. Kelvin body force arises from this nonuniformity and induces vortexes near the sphere. These vortexes disturb the boundary... 

This paper presents a new procedure for determining the fracture toughness of rock-like specimens using the diametric compression test with the center-cracked horseshoe disk (CCHD) method. Using finite element analysis, a dimensionless stress intensity factor was obtained and a formula was rendered for determining mode I fracture toughness. To evaluate the accuracy of the measurement results produced by the CCHD method, fracture toughness experiments were conducted on the same rock-like material using the notched Brazilian disk (NBD) method. The CCHD tests were simulated using a two-dimensional particle flow code for validation of the experimental results, and a great agreement between the... 

Low-crested breakwaters have become more attractive because they do not hinder the beautiful coastal landscape and are more environmentally friendly comparing to traditional breakwaters. The available formulae for design of these types of structures mainly focus on wave transmission, reflection and dissipation aspects. However, stability of the armour blocks are not directly taken into account in the design process. In this study the Lagrangian meshfree method of SPH along with a discrete element method i.e. Polyhedral DEM are used to simulate the interaction between the sea water waves with different characteristics and low-crested breakwaters constructed using cubic armour blocks. SPH is... 

In this article, the behavior of gas bubbles in tapered fluidized beds is investigated with the use of a two-fluid model incorporating kinetic theory of granular flow. The effects of various parameters such as apex angle, particle size, and particle density on the size distribution and the rise velocity of gas bubbles were examined. In addition, the simulation results for the bubble fraction and axial velocity of gas bubbles were compared with experimental data reported in the literature and good agreement was observed. As the apex angle was increased, the fraction of gas bubbles with large sizes increased and the fraction of bubbles with small sizes decreased. As the particle size... 

The purpose of the present study is to investigate the thermal performance of a photovoltaic/thermal system, integrated with phase change materials in porous medium. For this purpose, a metal foam was employed as porous medium and the performance of five different PCMs, as organic and inorganic, were examined as well. Moreover, the effects of different key parameters such as the mass flow rate, solar irradiance, inlet water temperature and inclination were studied. Finally, the simulation results were compared with a water-cooled photovoltaic/thermal without incorporating PCMs and porous medium, and thermal performance of the three PV/T cases were reported. The highest thermal efficiency of... 

This study investigates GH2/LOX coaxial jet flame at trans- and supercritical conditions using the Reynolds averaged Navier-Stokes approach. Four two-equation-turbulence models, three real equation of states, two chemical mechanisms, and three different chamber pressures are examined. Predictions show good agreement with measurements qualitatively and quantitatively. Based on the results, the predictions of the Soave-Redlich-Kwong equation of state (EOS) are closer to the experiment, while the Aungier-Redlich-Kwong EOS has more deviation than the others. Moreover, the k-ω shear stress transport model has better performance than the other turbulence models. It is also found that the flow... 

Cavitating flow is investigated around marine propellers, experimentally and numerically. Two different types of conventional model propellers are used for the study. The first one is a four bladed model propeller, so called model A, and the second one is a three bladed propeller, model B. Model A is tested in different cavitation regimes in a K23 cavitation tunnel. The results are presented in characteristic curves and related pictures. Finally, the results are discussed. Model B is investigated based on existing experimental results. In addition, model B is used for validation of the numerical solution prior to the testing of model A. The cavitation phenomenon is predicted numerically on a... 

In this research the fluid dynamics characteristics of a stellar turbulent jet flow is studied numerically and the results of three dimensional jet issued from a stellar nozzle are presented. A numerical method based on control volume approach with collocated grid arrangement is employed. The turbulent stresses are approximated using k-ε and k-ω models with four different inlet conditions. The velocity field is presented and the rate of decay at jet centerline is noted. Special attention is drawn on the influence of corner angle and number of wings on mixing in stellar cross section jets. Stellar jets with three; four and five wings and 15-65° corner angles are studied. Also the effect of...