Sharif Digital Repository

In recent years, there has been a growing interest in the simulation and analysis of piezoelectric transducers with the he help of equivalent electrical circuit simulations (EECS). This paper has been devoted to study of such approach for two designed and fabricated ultrasonic sandwich transducers. By using analytical analysis, the dimensions of components of the two piezoelectric transducers were determined for the assumed resonance frequencies of 30 kHz and 40 kHz. Then, by using a two dimensional finite element model, and regarding two different modeling techniques for the transducers and by application of a current source which was connected directly to the piezoelectric pieces of the... 

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin... 

Numerical modeling of the fully coupled phenomena of solid deformation-fluid flow in partially saturated porous media is of great interest in many branches of science and engineering. In this study, a new formulation based on one of the famous mesh-less methods, called Element-Free Galerkin (EFG), is developed to simulate the water and air movement through variably saturated soils. For this purpose, the governing partial differential equations including the equilibrium equation and mass conservation laws for each fluid phase are discretized in space using the same EFG shape functions. To enforce the essential boundary conditions, penalty method is employed. Temporal discretization is... 

In the past few decades, numerical simulation of multiphase flow systems has received increasing attention because of its importance in various fields of science and engineering. In this paper, a three-dimensional numerical model is developed for the analysis of simultaneous flow of two fluids through porous media. The numerical approach is fairly new based on the element-free Galerkin (EFG) method. The EFG is a type of mesh-less method which has rarely been used in the field of flow in porous media. The weak forms of the governing partial differential equations are derived by applying the weighted residual method and Galerkin technique. The penalty method is utilized for imposition of the... 

Meshless methods are a relatively new type of numerical methods that have attracted the attention of many researchers over the past years. So far, a number of meshless methods have been developed and applied to solve problems in various fields of engineering, including solid mechanics and geotechnical problems. The Element-Free Galerkin (EFG) method is adopted in this study for solving the governing partial differential equations of equilibrium and continuity of pore fluid flow for numerical simulation of coupled hydro-mechanical problems. For this purpose, the weak form of the governing equations is derived by applying the weighted residual method and Galerkin technique. The penalty method... 

The soil permeability coefficient plays a key role in the process of numerical simulation of the liquefaction phenomenon. Liquefaction causes a considerable increase in soil permeability, due to the creation of easier paths for water flow. The work presented in this paper tries to investigate the effects of permeability coefficient on the results of numerical modeling of the liquefaction phenomenon. To do this, a fully coupled (u-P) formulation is employed to analyze soil displacements and pore water pressures. Two different versions of a well-calibrated critical state bounding surface plasticity model, which possesses the capability to utilize a single set of material parameters for a wide... 

Liquefaction phenomenon is usually accompanied by large amounts of settlement owing to disruption of soil structure. In addition to that, large settlement also occurs by a significant increase in soil permeability during seismic excitation. To properly simulate the post-liquefaction settlement, it is important to take the compressibility properties of the liquefied sand as well as the permeability increase into account. Using initial permeability coefficient in the course of simulation of liquefaction leads to underestimation of settlement. In addition to that, using unrealistic values for permeability may cause erroneous predictions of other aspects of soil behavior. Therefore, an accurate... 

An experimental study is conducted to measure small-strain shear modulus of clay-cement mixture using bender element apparatus setup in a triaxial cell. Bender element tests were conducted on cement-treated soils and the results were analyzed to study the variation of shear modulus properties of soil specimens at different cement contents, confining pressures, curing times, and compaction moisture contents. Based on the obtained results increasing the cement ratio has a significant effect on the small-strain shear modulus of the treated soils, and this effect signifies with increasing the moisture content and curing time. Rates of shear modulus enhancements due to cement content, curing... 

In this paper, a fully coupled three-dimensional dynamic analysis is carried out to investigate the dynamic behavior of pile foundations in liquefied ground. A critical state bounding surface plasticity model is used to model soil skeleton, while a fully coupled (u- P) formulation is employed to analyze soil displacements and pore water pressures. Furthermore, in this study, variation of permeability coefficient during liquefaction is taken into account; the permeability coefficient is related to excess pore water pressure ratio. Results of a centrifuge test on pile foundations are used to demonstrate the capability of the model for reliable analysis of piles under dynamic loading. Then, the... 

In this study, a coupled numerical approach based on Lattice Boltzmann Method (LBM) and Discrete Element Method (DEM) is employed for 2D simulation of fluid flow in porous media comprising of movable circular particles. The developed model is used for simulation of sand production which is one of the important problems in petroleum industry. The numerical tool has proved to have the capability of investigating the mechanisms involved in sand production problem. The results show that the rate of sand production is strongly affected by flow rate and confining pressure  

In this paper, Lattice Boltzmann method (LBM) has been used to study the effects of permeability and tortuosity on flow through saturated particulate media and identify the relationships between permeability and tortuosity with other parameters such as particles diameter, grain specific surface, and porosity. LBM is a simple kinematic model that can incorporate the essential physics of microscopic and mesoscopic processes involved in flow through granular soils. The obtained results indicate that the 2D LB model, due to its inherent theoretical advantages, is capable of demonstrating that the porosity and specific surface are the most influential parameters in determining the intrinsic... 

Relative permeability curves have practical implications in petroleum reservoir simulations. Study of the effects of reservoir wettability, pore shape geometry, and viscosity ratio of flowing fluids on the relative permeabilities is of great importance in reservoir modeling. In this paper, lattice Boltzmann method (LBM) is employed for analyzing the two-fluid flow in rigid porous media. The developed LBM code proved to be a robust numerical tool for analyzing the factors that influence the relative permeabilities of two immiscible fluids flowing through porous media. The numerically derived relative permeability curves demonstrate that in neutrally wet reservoirs, the effect of viscosity... 

In this study, a coupled numerical approach based on Lattice Boltzmann Method (LBM) and Discrete Element Method (DEM) is employed for two-dimensional simulation of fluid flow in deformable particulate media comprising of movable circular particles. The developed LB-DE code is validated against the results of a bi-axial shear test as well as two well-known benchmark problems including settling of a circular particle under gravity force inside a viscous fluid, and motion of a neutrally buoyant particle released in a Poiseuille flow. The verified code is then utilized for simulation of "Sand Production" phenomenon which is of importance for oil producing wells in weakly cemented sandstone... 

Abstract Lateral Spreading, which usually occurs as a consequence of liquefaction in gently sloped loose saturated sand layers, is known to be a major source of earthquake-induced damages to structures such as quay walls, bridge piers, pipelines, and highway/railways. Therefore evaluation of the liquefaction potential and using appropriate methods for prediction of the adverse consequences of lateral spreading is of great importance. In this study, numerical modeling has been used to study lateral spreading phenomenon behind rigid waterfront structures. Coupled dynamic field equations of the extended Biot's theory with u-P formulation are used for simulating the phenomenon. A fully coupled... 

Permeability of porous materials is an important characteristic which is extensively used in various engineering disciplines. There are a number of issues that influence the permeability coefficient among which the porosity, size of particles, pore shape, tortuosity, and particle size distribution are of great importance. In this paper a C++ GPU code based on three-dimensional lattice Boltzmann method (LBM) has been developed and used for investigating the effects of the above mentioned factors on the permeability coefficient of granular materials. Multirelaxation time collision scheme of the LBM equations is used in the simulator, which is capable of modeling the exact position of the... 

In this paper, a large plasticity deformation finite element modeling is presented for three-dimensional dynamic analysis of unsaturated soils with special reference to the failure of lower San Fernando dam under the 1971 earthquake. The finite element method is applied to the governing equations for the spatial discretization, followed by a generalized Newmark scheme used for the time domain discretization. Time stepping scheme is used in the fully implicit coupled method and a direct solution procedure is used for the coupled equation system. The framework of generalized plasticity is presented and the numerical results of unsaturated soils are demonstrated based on the Pastor-Zienkiewicz... 

Occurrence of liquefaction in saturated sand deposits underlying foundation of structure can cause a wide range of structural damages starting from minor settlement, and ending to general failure due to loss of bearing capacity. If the bearing capacity failure is not the problem, reliable estimation of the liquefaction-induced settlement will be of prime importance in assessment of the overall performance of the structure. Currently, there are few procedures with limited application in practice for estimation of settlement of foundations on liquefied ground. Therefore, development of a general relationship is important from the practical viewpoint. In this paper, the dynamic response of... 

Dynamic compaction is a widely used method for improvement of loose granular deposits. Its applicability in saturated layers generally considered to be less effective because of the fact that part of the applied energy is absorbed by pore water. Up to now the majority of numerical simulations have focused on the analysis of dynamic compaction in dry/moist soils. In this paper, a fully coupled hydro-mechanical finite element code has been developed and employed to evaluate the dynamic compaction effects on saturated granular soils. After verification of the results by comparing the numerical results with those measured in a real field case of DC treatment in a highway, some sensitivity... 

Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi-analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element-free Galerkin (EFG) mesh-less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity... 

In performance-based geotechnical earthquake engineering, the required degree and spatial extent of ground densification for mitigation of liquefaction beneath a structure should be determined based on the acceptable levels of performance of foundation. Currently, there is no solution for evaluation of the amount of settlement and tilt of footings constructed on a densified ground which is surrounded by a liquefiable soil. This implies the need for numerical procedures for simulation of seismic behavior of shallow foundations supported on both liquefiable and densified subsoil. In this paper, the dynamic response of shallow foundations on a densified ground is studied using a 3D fully...