Sharif Digital Repository

Hydrodynamics of highly nonlinear cnoidal waves and their subsequent strong plunging breakers are among the least understood and most significant issues in coastal engineering. In this work, a weakly compressible smoothed particle hydrodynamics (SPH) formulation is used for the study of the generation and propagation of cnoidal waves and investigation of the characteristics of the induced strong plunging breakers. Numerical results show the capability of the SPH scheme for properly simulating the cnoidal waves. For the case of strong plunging breakers, dynamic and kinematic features of the flow are computed and compared with certain implementations of other numerical techniques. SPH is shown... 

Vibration suppression of a strain gradient Euler-Bernoulli beam in presence of disturbance and uncertainties is considered in this investigation. Vibration of the system is suppressed by an adaptive boundary controller which has robustness to the environmental and control effort disturbances. The direct Lyapunov stability theorem is used to design the controller and adaptation law. The numerical results are presented to demonstrate the effectiveness of the proposed controller. Copyright © 2016 by ASME  

Surface turbulence during the filling of the mold triggers the entrainment of oxide films, which appears to be detrimental to the soundness of the final casting. Nonpressurized and bottom-gating systems have been employed in order to avoid such casting defects by reducing the surface velocity of the liquid metal. However, recent studies have shown that the melt front velocity in the mold entrance may exceed the critical value in the nonpressurized and bottom-gating systems. Therefore, a study was conducted on numerical simulation melt flow pattern in nonpressurized and bottom-gating systems. It was noted that the liquid metal enters the gate and then mold cavity with a higher velocity by... 

Sloshing, or liquid free surface oscillation, in containers has many important applications in a variety of engineering fields. The modal method can be used to solve linear sloshing problems and is the most efficient reduced order method that has been used during the previous decade. In the present article, the modal method is used to solve a nonlinear sloshing problem. The method is based on a potential flow solution that implements a two-phase analysis on sloshing in a rectangular container. According to this method, the solution to the mass conservation equation, with a nonpenetration condition at the tank walls, results in velocity potential expansion; this is similar to the mode shapes... 

Studies on two-phase flow in small scale pipes have become more important, because of the application of mini-scale devices in several engineering fields including, high heat-flux compact heat exchangers, and cooling systems of various types of equipment. In a mini pipe the behavior of two phase flow is not the same as flow in conventional pipes. The difference is caused by different effective forces; for e. g. inside a mini pipe capillary forces are more important in comparison with gravitational forces. This paper is devoted to numerical simulation of gas-liquid two phase flow in a vertical mini pipe. Prediction of bubble shape and the effects of gas and liquid velocities on flow... 

In this paper, a two-dimensional numerical approach is used to study the effect of micro combustor height, 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 is used. The transient gas phase energy and species conservation equations result in an Advection-Diffusion-Reaction system that leads to two stiff systems of PDEs. In the present work, the computational domain is solved through the Strang splitting method, which is suitable for a nonlinear stiff system of PDEs. A revised boundary condition for the velocity equation is applied and its effect on the flow... 

We theoretically demonstrate the tunabiltiy of terahertz random lasers composed of high temperature superconductor YBCO and ruby layers as active medium. The considered system is a one-dimensional disordered medium made of ruby grain and YBCO. Finite-difference time domain method is used to calculate the emission spectrum and spatial distribution of electric field at different temperatures. Our numerical results reveal that the superconductor based random lasers exhibit large temperature tunability in the terahertz domain. The emission spectrum is significantly temperature dependent, the number of lasing modes and their intensities increase with decreasing temperature. Also, we make some... 

A simple and efficient transmission line model is proposed here to study how the transmission characteristics of photonic crystal waveguides are tailored by introduction of stubs patterned in the photonic crystal lattice. It is shown that band-pass and band-stop optical filters can be easily designed and optimized when stubs of appropriate length are brought in. Since the lengths of the designed stubs are not necessarily integer multiples of the photonic crystal lattice constant, a geometric shift in a portion of the photonic crystal structure is shown to be essential. The proposed model is verified by using a rigorous numerical method. An excellent agreement is observed between the... 

The effect of surface and interface elasticity in analysis of the Saint-Venant torsion problem of an eccentrically two-phase circular nanobar is considered. The problem is formulated in the context of Gurtin's surface elasticity. For a rigorous solution of the proposed problem, conformal mapping together with a Laurent series expansion are employed. At the nanoscales the usual classical theories cease to hold and the corresponding results deteriorate. The numerical results well illustrate that the torsional rigidity of the mentioned nanosized structural elements are significantly affected by the size. Some applications of the given results can be contemplated in the design of micro/nano... 

In this letter, we propose a novel communication scheme for macro-scale molecular timing channels (MTCs), where the information carrier is the release time of the molecules. Next, the symbol constellation, receiver's sampling time, and receiver's decision rule are designed and analyzed. We further propose a computationally simple yet optimal receiver. Our numerical results show that macro-scale MTC is an approach for achieving ultra-reliable molecular communication. © 2012 IEEE  

The objective of the dynamic economic dispatch (DED) problem is to schedule power generation for the online units for a given time horizon economically, satisfying various operational constraints. Due to the effect of valve-point effects and prohibited operating zones (POZs) in the generating units cost functions, DED problem is a highly non-linear and non-convex optimization problem. The DED problem even may be more complicated if transmission losses and ramp-rate constraints are taken into account. This paper presents a novel and heuristic algorithm to solve DED problem of generating units, by employing time varying acceleration coefficients iteration particle swarm optimization... 

Two-phase flow simulations around a body were not studied before and considering these flows can play a significant role in long-term reliability and safety of industrial systems. In this paper, flow regimes, drag coefficient and void fraction around different cross-section prisms were considered. To achieve this aim, main equations of flow have been developed for investigation of drag coefficient in air-water two phase. Our numerical analyses were performed by a designed and written CFD package which is based on Eulerian-Eulerian approach. Geometries, which have been studied in this article, are: circle, rectangle and triangle, for different aspect ratio (length over width) and leading edge... 

In this paper, we evaluate the effect of power randomization on the saturation throughput of a WLAN comprised of four classes of users, corresponding to four access categories (ACs) of IEEE 802.11e EDCA. To this end, we consider a specific number of power levels for all users. By some simplifications, we find the probability distributions for the power levels at each class in order to optimize the saturation throughput of a specific class. We evaluate the effect of such optimization onto the saturation throughput of other ACs as well as the saturation throughput of the whole network. Finally, we present some numerical results to show the effect of our optimization in different conditions  

The urine formation and excretion system have long been of interest for mathematicians and physiologists to elucidate the obscurities within the process happens in renal tissue. In this study, a novel three-dimensional approach is utilized for modeling the urine concentrating mechanism in rat renal outer medulla which is essentially focused on demonstrating the significance of tubule's architecture revealed in anatomic studies and physiological literature. Since nephrons and vasculatures work interdependently through a highly structured arrangement in outer medulla which is dominated by vascular bundles, a detailed functional unit is proposed based on this specific configuration.... 

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

A three-dimensional elasticity solution for the analysis of functionally graded rotating cylinders with variable thickness profile is proposed. The axisymmetric structure has been divided in several divisions in the radial direction. Constant mechanical properties and thickness profile are assumed within each division. The solution is considered for four different thickness profiles, namely constant, linear, concave, and convex. It is shown that the linear, concave, and convex thickness profiles have smaller stress values compared to a constant thickness profile. The effects of various grading indices as well as different boundary conditions, namely solid, free-free hollow and fixed-free... 

An analytical derivation of the elastodynamic fundamental solutions for a transversely isotropic tri-material full-space is presented by means of a complete representation using two displacement potentials. The complete set of three-dimensional point-load, patch-load, and ring-load Green's functions for stresses and displacements are given, for the first time, in the complex-plane line-integral representations. The formulation includes a complete set of transformed stress-potential and displacement-potential relations in the framework of Fourier expansions and Hankel integral transforms, that is useful in a variety of elastodynamic as well as elastostatic problems. For the numerical... 

In this paper, the three-dimensional adaptive finite element modeling is presented for cohesive fracture analysis of non-planer crack growth. The technique is performed based on the Zienkiewicz-Zhu error estimator by employing the modified superconvergent patch recovery procedure for the stress recovery. The Espinosa-Zavattieri bilinear constitutive equation is used to describe the cohesive tractions and displacement jumps. The 3D cohesive fracture element is employed to simulate the crack growth in a non-planar curved pattern. The crack growth criterion is proposed in terms of the principal stress and its direction. Finally, several numerical examples are analyzed to demonstrate the... 

In this paper a saturable absorber medium is employed as an optical limiter to reduce the spot size to the range of several tens of nanometres. The characteristics of a Gaussian beam are theoretically analysed upon propagation through the saturable absorber medium. Based on Maxwell equations a system of coupled nonlinear ordinary differential equations for intensity, beam radius and beam curvature is obtained. Theoretical analyses and numerical results show that the behaviour of a Gaussian beam in a saturable absorber medium strongly depends on the initial characteristics of the laser beam. Numerical results indicate that, depending on the initial conditions and a suitable saturable absorber... 

In this work, we present a level set method to simulate steady and unsteady mass transfer from a single droplet moving in a second phase fluid under buoyant force. We initially use level set to determine the interface between the two phases, where the shape of drop forms. Next, we extend this method to solve a unique mass transfer equation for the entire solution domain without considering the discontinuity appeared at the interface. We use a finite element method incorporated with the characteristic-based split (CBS) algorithm to implement axi-symmetric mass transfer equations on a stationary Eulerian grid. Of course, the convection-diffusion modeling of mass transfer is different from the...