Numerical Analysis of Flow through Abrasive Water Suspension Jet: The Effect of Garnet, Aluminum Oxide and Silicon Carbide Abrasive on Skin Friction Coefficient Due To Wall Shear and Jet Exit Kinetic Energy

It is well known that the abrasive particles in the abrasive water suspension has significant effect on the erosion characteristics of the inside surface of the nozzle. Abrasive particles moving with the flow cause severe skin friction effect, there by altering the nozzle diameter due to wear which in turn reflects on the life of the nozzle for effective machining. Various commercial abrasives are available for abrasive water jet machining. The erosion characteristic of each abrasive is different. In consideration of this aspect, in the present work, the effect of abrasive materials namely garnet, aluminum oxide and silicon carbide on skin friction coefficient due to wall shear stress and jet kinetic energy has been analyzed. It is found that the abrasive material of lower density produces a relatively higher skin friction effect and higher jet exit kinetic energy.

Study of Natural Convection in a Triangular Cavity Filled with Water: Application of the Lattice Boltzmann Method

The Lattice Boltzmann Method (LBM) with double populations is applied to solve the steady-state laminar natural convective heat transfer in a triangular cavity filled with water. The bottom wall is heated, the vertical wall is cooled, and the inclined wall is kept adiabatic. The buoyancy effect was modeled by applying the Boussinesq approximation to the momentum equation. The fluid velocity is determined by D2Q9 LBM and the energy equation is discritized by D2Q4 LBM to compute the temperature field. Comparisons with previously published work are performed and found to be in excellent agreement. Numerical results are obtained for a wide range of parameters: the Rayleigh number from  to  and the inclination angle from 0° to 360°. Flow and thermal fields were exhibited by means of streamlines and isotherms. It is observed that inclination angle can be used as a relevant parameter to control heat transfer in right-angled triangular enclosures.  

Oil Refineries Emissions: Source and Impact: A Study using AERMOD

The main objectives of this paper are to measure pollutants concentrations in the oil refinery area in Kuwait over three periods during one year, obtain recent emission inventory for the three refineries of Kuwait, use AERMOD and the emission inventory to predict pollutants concentrations and distribution, compare model predictions against measured data, and perform numerical experiments to determine conditions at which emission rates and the resulting pollutant dispersion is below maximum allowable limits.

Performance Evaluation of the Post-Installed Anchor for Sign Structure

Numerous experimental tests for post-installed anchor systems drilled in hardened concrete were conducted in order to estimate pull-out and shear strength accounting for uncertainties such as torque ratios, embedment depths and different diameters in demands. In this study, the strength of the systems was significantly changed by the effect of those three uncertainties during pull-out experimental tests, whereas the shear strength of the systems was not affected by torque ratios. It was also shown that concrete cone failure or damage mechanism was generally investigated during and after pull-out tests and in shear strength tests, mostly the anchor systems were failed prior to failure of primary structural system. Furthermore, 3D finite element model for the anchor systems was created by ABAQUS for the numerical analysis. The verification of finite element model was identical till the failure points to the load-displacement relationship specified by the experimental tests.

Prognostic and Diagnostic Modes of Mathematical Model for the Pre-operation of Suspended Sediment Transport model in Estuaries and Coastal areas

Both prognostic and diagnostic modes of a 3D baroclinic model in hydrodynamic and sediment transport models of the Princeton Ocean Model (POM) were conducted to separate prognose and diagnose effects of different hydrodynamic factors on transport of suspended sediment discharged from the rivers to the Gulf of Thailand (GoT). Both transport modes of suspended sediment distribution in the GoT were numerically simulated. It could be concluded that the suspended sediment discharged from the rivers around the GoT. Most of sediments in estuaries and coastal areas are deposited outside the GoT under the condition of wind-driven current, and very small amount of the sediments of them are transported faraway. On the basis of wind forcing, sediments from the lower GoT to the upper GoT are mainly transported south-northwestward and also continuously moved north-southwestward. An obvious 3D characteristic of suspended sediment transport is produced in the wind-driven current residual circulation condition. In this study, the transport patterns at the third layer are generally consistent with the typhoon-induced strong currents in two case studies of Typhoon Linda 1997. The case studies presented the prognostic and diagnostic modes during 00UTC28OCT1997 to 12UTC06NOV1997 in a short period with the current condition for pre-operation of the suspended sediment transport model in estuaries and coastal areas.

Numerical Study of Natural Convection Effects in Latent Heat Storage using Aluminum Fins and Spiral Fillers

A numerical investigation has carried out to understand the melting characteristics of phase change material (PCM) in a fin type latent heat storage with the addition of embedded aluminum spiral fillers. It is known that melting performance of PCM can be significantly improved by increasing the number of embedded metallic fins in the latent heat storage system but to certain values where only lead to small improvement in heat transfer rate. Hence, adding aluminum spiral fillers within the fin gap can be an option to improve heat transfer internally. This paper presents extensive computational visualizations on the PCM melting patterns of the proposed fin-spiral fillers configuration. The aim of this investigation is to understand the PCM-s melting behaviors by observing the natural convection currents movement and melting fronts formation. Fluent 6.3 simulation software was utilized in producing twodimensional visualizations of melting fractions, temperature distributions and flow fields to illustrate the melting process internally. The results show that adding aluminum spiral fillers in Fin type latent heat storage can promoted small but more active natural convection currents and improve melting of PCM.

A Fully Implicit Finite-Difference Solution to One Dimensional Coupled Nonlinear Burgers’ Equations

A fully implicit finite-difference method has been proposed for the numerical solutions of one dimensional coupled nonlinear Burgers’ equations on the uniform mesh points. The method forms a system of nonlinear difference equations which is to be solved at each iteration. Newton’s iterative method has been implemented to solve this nonlinear assembled system of equations. The linear system has been solved by Gauss elimination method with partial pivoting algorithm at each iteration of Newton’s method. Three test examples have been carried out to illustrate the accuracy of the method. Computed solutions obtained by proposed scheme have been compared with analytical solutions and those already available in the literature by finding L2 and L∞ errors.

The Contraction Point for Phan-Thien/Tanner Model of Tube-Tooling Wire-Coating Flow

The simulation of extrusion process is studied widely in order to both increase products and improve quality, with broad application in wire coating. The annular tube-tooling extrusion was set up by a model that is termed as Navier-Stokes equation in addition to a rheological model of differential form based on singlemode exponential Phan-Thien/Tanner constitutive equation in a twodimensional cylindrical coordinate system for predicting the contraction point of the polymer melt beyond the die. Numerical solutions are sought through semi-implicit Taylor-Galerkin pressurecorrection finite element scheme. The investigation was focused on incompressible creeping flow with long relaxation time in terms of Weissenberg numbers up to 200. The isothermal case was considered with surface tension effect on free surface in extrudate flow and no slip at die wall. The Stream Line Upwind Petrov-Galerkin has been proposed to stabilize solution. The structure of mesh after die exit was adjusted following prediction of both top and bottom free surfaces so as to keep the location of contraction point around one unit length which is close to experimental results. The simulation of extrusion process is studied widely in order to both increase products and improve quality, with broad application in wire coating. The annular tube-tooling extrusion was set up by a model that is termed as Navier-Stokes equation in addition to a rheological model of differential form based on single-mode exponential Phan- Thien/Tanner constitutive equation in a two-dimensional cylindrical coordinate system for predicting the contraction point of the polymer melt beyond the die. Numerical solutions are sought through semiimplicit Taylor-Galerkin pressure-correction finite element scheme. The investigation was focused on incompressible creeping flow with long relaxation time in terms of Weissenberg numbers up to 200. The isothermal case was considered with surface tension effect on free surface in extrudate flow and no slip at die wall. The Stream Line Upwind Petrov-Galerkin has been proposed to stabilize solution. The structure of mesh after die exit was adjusted following prediction of both top and bottom free surfaces so as to keep the location of contraction point around one unit length which is close to experimental results.

A Comparison Study of a Symmetry Solution of Magneto-Elastico-Viscous Fluid along a Semi- Infinite Plate with Homotopy Perturbation Method and4th Order Runge–Kutta Method

The equations governing the flow of an electrically conducting, incompressible viscous fluid over an infinite flat plate in the presence of a magnetic field are investigated using the homotopy perturbation method (HPM) with Padé approximants (PA) and 4th order Runge–Kutta method (4RKM). Approximate analytical and numerical solutions for the velocity field and heat transfer are obtained and compared with each other, showing excellent agreement. The effects of the magnetic parameter and Prandtl number on velocity field, shear stress, temperature and heat transfer are discussed as well.

Hexagonal Honeycomb Sandwich Plate Optimization Using Gravitational Search Algorithm

Honeycomb sandwich panels are increasingly used in the construction of space vehicles because of their outstanding strength, stiffness and light weight properties. However, the use of honeycomb sandwich plates comes with difficulties in the design process as a result of the large number of design variables involved, including composite material design, shape and geometry. Hence, this work deals with the presentation of an optimal design of hexagonal honeycomb sandwich structures subjected to space environment. The optimization process is performed using a set of algorithms including the gravitational search algorithm (GSA). Numerical results are obtained and presented for a set of algorithms. The results obtained by the GSA algorithm are much better compared to other algorithms used in this study.

A Study on the Effect of Variation of the Cross-sectional Area of Spiral Volute Casing for Centrifugal Pump

The impeller and the casing are the key components of a centrifugal pump. Although there have been many studies on the impeller and the volute casing of centrifugal pump, further study of the volute casing to improve the performance of centrifugal pumps is needed. In this paper, the effect of cross-sectional area on the performance of volute casing was investigated using a commercial CFD code. The performance characteristics, not only at the off-design point but also for a full type model are required these days. So we conducted numerical analysis for all operating points by using complete geometry through transient analysis. Transient analysis on the complete geometry of a real product has the advantage of simulating realistic flow. The results of this study show the variation of a performance curve by modifying the above-mentioned design parameter.

Finite Element Analysis for Damped Vibration Properties of Panels Laminated Porous Media

A numerical method is proposed to calculate damping properties for sound-proof structures involving elastic body, viscoelastic body, and porous media. For elastic and viscoelastic body displacement is modeled using conventional finite elements including complex modulus of elasticity. Both effective density and bulk modulus have complex quantities to represent damped sound fields in the porous media. Particle displacement in the porous media is discretised using finite element method. Displacement vectors as common unknown variables are solved under coupled condition between elastic body, viscoelastic body and porous media. Further, explicit expressions of modal loss factor for the mixed structures are derived using asymptotic method. Eigenvalue analysis and frequency responded were calculated for automotive test panel laminated viscoelastic and porous structures using this technique, the results almost agreed with the experimental results.

Computation of the Filtering Properties of Photonic Crystal Waveguide Discontinuities Using the Mode Matching Method

In this paper, the application of the Mode Matching (MM) method in the case of photonic crystal waveguide discontinuities is presented. The structure under consideration is divided into a number of cells, which supports a number of guided and evanescent modes. These modes can be calculated numerically by an alternative formulation of the plane wave expansion method for each frequency. A matrix equation is then formed relating the modal amplitudes at the beginning and at the end of the structure. The theory is highly efficient and accurate and can be applied to study the transmission sensitivity of photonic crystal devices due to fabrication tolerances. The accuracy of the MM method is compared to the Finite Difference Frequency Domain (FDFD) and the Adjoint Variable Method (AVM) and good agreement is observed.

Joint Use of Factor Analysis (FA) and Data Envelopment Analysis (DEA) for Ranking of Data Envelopment Analysis

This article combines two techniques: data envelopment analysis (DEA) and Factor analysis (FA) to data reduction in decision making units (DMU). Data envelopment analysis (DEA), a popular linear programming technique is useful to rate comparatively operational efficiency of decision making units (DMU) based on their deterministic (not necessarily stochastic) input–output data and factor analysis techniques, have been proposed as data reduction and classification technique, which can be applied in data envelopment analysis (DEA) technique for reduction input – output data. Numerical results reveal that the new approach shows a good consistency in ranking with DEA.

Reflection of Plane Waves at Free Surface of an Initially Stressed Dissipative Medium

The paper discuses the effect of initial stresses on the reflection coefficients of plane waves in a dissipative medium. Basic governing equations are formulated in context of Biot's incremental deformation theory. These governing equations are solved analytically to obtain the dimensional phase velocities of plane waves propagating in plane of symmetry. Closed-form expressions for the reflection coefficients of P and SV waves- incident at the free surface of an initially stressed dissipative medium are obtained. Numerical computations, using these expressions, are carried out for a particular model. Computations made with the results predicted in presence and absence of the initial stresses and the results have been shown graphically. The study shows that the presence of compressive initial stresses increases the velocity of longitudinal wave (P-wave) but diminishes that of transverse wave (SV-wave). Also the numerical results presented indicate that initial stresses and dissipation might affect the reflection coefficients significantly.

RBF- based Meshless Method for Free Vibration Analysis of Laminated Composite Plates

The governing differential equations of laminated plate utilizing trigonometric shear deformation theory are derived using energy approach. The governing differential equations discretized by different radial basis functions are used to predict the free vibration behavior of symmetric laminated composite plates. Effect of orthotropy and span to thickness ratio on frequency parameter of simply supported laminated plate is presented. Numerical results show the accuracy and good convergence of radial basis functions.

Verification of K-ω SST Turbulence Model for Supersonic Internal Flows

In this work, we try to find the best setting of Computational Fluid Dynamic solver available for the problems in the field of supersonic internal flows. We used the supersonic air-toair ejector to represent the typical problem in focus. There are multiple oblique shock waves, shear layers, boundary layers and normal shock interacting in the supersonic ejector making this device typical in field of supersonic inner flows. Modeling of shocks in general is demanding on the physical model of fluid, because ordinary conservation equation does not conform to real conditions in the near-shock region as found in many works. From these reasons, we decided to take special care about solver setting in this article by means of experimental approach of color Schlieren pictures and pneumatic measurement. Fast pressure transducers were used to measure unsteady static pressure in regimes with normal shock in mixing chamber. Physical behavior of ejector in several regimes is discussed. Best choice of eddy-viscosity setting is discussed on the theoretical base. The final verification of the k-ω SST is done on the base of comparison between experiment and numerical results.

Evaluation of Geosynthetic Forces in GRSRW under Dynamic Condition

Geosynthetics have proved to be suitable for reinforced soil retaining walls. Based on the increasing uses of geosynthetic reinforced soil systems in the regions, which bear frequent earthquakes, the study of dynamic behavior of structures seems necessary. Determining the reinforcement forces is; therefore, one of the most important and main points of discussions in designing retaining walls, by which we prevent from conservative planning. Thus, this paper intended to investigate the effects of such parameters as wall height, acceleration type, vertical spacing of reinforcement, type of reinforcement and soil type on forces and deformation through numerical modeling of the geosynthetic reinforced soil retaining walls (GRSRW) under dynamic loading with finite difference method by using FLAC. The findings indicate rather positive results with each parameter.

Numerical Modeling and Computer Simulation of Ground Movement above Underground Mine

This paper describes topic of computer simulation with regard to the ground movement above an underground mine. Simulation made with software package ADINA for nonlinear elastic-plastic analysis with finite elements method. The one of representative profiles from Mine 'Stara Jama' in Zenica has been investigated. A collection and selection of both geo-mechanical data and geometric parameters of the mine was necessary for performing these simulations. Results of estimation have been compared with measured values (vertical displacement of surface), and then simulation performed with assumed dynamic and dimensions of excavation, over a period of time. Results are presented with bitmaps and charts.

Capsule-substrate Adhesion in the Presence of Osmosis by the Immersed Interface Method

A two-dimensional thin-walled capsule of a flexible semi-permeable membrane is adhered onto a rigid planar substrate under adhesive forces (derived from a potential function) in the presence of osmosis across the membrane. The capsule is immersed in a hypotonic and diluted binary solution of a non-electrolyte solute. The Stokes flow problem is solved by the immersed interface method (IIM) with equal viscosities for the enclosed and surrounding fluid of the capsule. The numerical results obtained are verified against two simplified theoretical solutions and the agreements are good. The osmotic inflation of the adhered capsule is studied as a function of the solute concentration field, hydraulic conductivity, and the initial capsule shape. Our findings indicate that the contact length shrinks in dimension as capsule inflates in the hypotonic medium, and the equilibrium contact length does not depend on the hydraulic conductivity of the membrane and the initial shape of the capsule.