Search results for: Numerical Schlieren
Commenced in January 2007
Frequency: Monthly
Edition: International
Paper Count: 3550

Search results for: Numerical Schlieren

3370 Numerical Simulation of Bio-Chemical Diffusion in Bone Scaffolds

Authors: Masoud Madadelahi, Amir Shamloo, Seyedeh Sara Salehi

Abstract:

Previously, some materials like solid metals and their alloys have been used as implants in human’s body. In order to amend fixation of these artificial hard human tissues, some porous structures have been introduced. In this way, tissues in vicinity of the porous structure can be attached more easily to the inserted implant. In particular, the porous bone scaffolds are useful since they can deliver important biomolecules like growth factors and proteins. This study focuses on the properties of the degradable porous hard tissues using a three-dimensional numerical Finite Element Method (FEM). The most important studied properties of these structures are diffusivity flux and concentration of different species like glucose, oxygen, and lactate. The process of cells migration into the scaffold is considered as a diffusion process, and related parameters are studied for different values of production/consumption rates.

Keywords: bone scaffolds, diffusivity, numerical simulation, tissue engineering

Procedia PDF Downloads 384
3369 Numerical Study on the Heat Transfer Characteristics of Composite Phase Change Materials

Authors: Gui Yewei, Du Yanxia, Xiao Guangming, Liu Lei, Wei Dong, Yang Xiaofeng

Abstract:

A phase change material (PCM) is a substance which absorbs a large amount of energy when undergoing a change of solid-liquid phase. The good physical and chemical properties of C or SiC foam reveal the possibility of using them as a thermal conductivity enhancer for the PCM. C or SiC foam composite PCM has a high effective conductivity and becomes one of the most interesting thermal storage techniques due to its advantage of simplicity and reliability. The paper developed a numerical method to simulate the heat transfer of SiC and C foam composite PCM, a finite volume technique was used to discretize the heat diffusion equation while the phase change process was modeled using the equivalent specific heat method. The effects of the porosity were investigated based on the numerical method, and the effects of the geometric model of the microstructure on the equivalent thermal conductivity was studies.

Keywords: SiC foam, composite, phase change material, heat transfer

Procedia PDF Downloads 510
3368 Numerical Investigation of Static and Dynamic Responses of Fiber Reinforced Sand

Authors: Sandeep Kumar, Mahesh Kumar Jat, Rajib Sarkar

Abstract:

Soil reinforced with randomly distributed fibers is an attractive means to improve the performance of soil in a cost effective manner. Static and dynamic characterization of fiber reinforced soil have become important to evaluate adequate performance for all classes of geotechnical engineering problems. Present study investigates the behaviour of fiber reinforced cohesionless soil through numerical simulation of triaxial specimen. The numerical model has been validated with the existing literature of laboratory triaxial compression testing. A parametric study has been done to find out optimum fiber content for shear resistance. Cyclic triaxial testing has been simulated and the stress-strain response of fiber-reinforced sand has been examined considering different combination of fiber contents. Shear modulus values and damping values of fiber-reinforced sand are evaluated. It has been observed from results that for 1.0 percent fiber content shear modulus increased 2.28 times and damping ratio decreased 4.6 times. The influence of amplitude of cyclic strain, confining pressure and frequency of loading on the dynamic properties of fiber reinforced sand has been investigated and presented.

Keywords: damping, fiber reinforced soil, numerical modelling, shear modulus

Procedia PDF Downloads 278
3367 Perforation Analysis of the Aluminum Alloy Sheets Subjected to High Rate of Loading and Heated Using Thermal Chamber: Experimental and Numerical Approach

Authors: A. Bendarma, T. Jankowiak, A. Rusinek, T. Lodygowski, M. Klósak, S. Bouslikhane

Abstract:

The analysis of the mechanical characteristics and dynamic behavior of aluminum alloy sheet due to perforation tests based on the experimental tests coupled with the numerical simulation is presented. The impact problems (penetration and perforation) of the metallic plates have been of interest for a long time. Experimental, analytical as well as numerical studies have been carried out to analyze in details the perforation process. Based on these approaches, the ballistic properties of the material have been studied. The initial and residual velocities laser sensor is used during experiments to obtain the ballistic curve and the ballistic limit. The energy balance is also reported together with the energy absorbed by the aluminum including the ballistic curve and ballistic limit. The high speed camera helps to estimate the failure time and to calculate the impact force. A wide range of initial impact velocities from 40 up to 180 m/s has been covered during the tests. The mass of the conical nose shaped projectile is 28 g, its diameter is 12 mm, and the thickness of the aluminum sheet is equal to 1.0 mm. The ABAQUS/Explicit finite element code has been used to simulate the perforation processes. The comparison of the ballistic curve was obtained numerically and was verified experimentally, and the failure patterns are presented using the optimal mesh densities which provide the stability of the results. A good agreement of the numerical and experimental results is observed.

Keywords: aluminum alloy, ballistic behavior, failure criterion, numerical simulation

Procedia PDF Downloads 312
3366 A Counter-flow Vortex Tube With Energy Separation: An Experimental Study and CFD Analysis

Authors: Li̇zan Mahmood Khorsheed Zangana

Abstract:

Experimental and numerical investigations have been carried out to study the mechanism of separation energy and flow phenomena in the counter-flow vortex tube. This manuscript presents a complete comparison between the experimental investigation and CFD analysis. The experimental model tested under different inlet pressures. Three-dimensional numerical modelling using the k-ε model. The results show any increase in both cold mass fraction and inlet pressure caused to increase ΔTc, and the maximum ΔTc value occurs at P = 6 bar. The coefficient of performance (COP) of two important factors in the vortex tube have been evaluated, which ranged from 0.25 to 0.74. The maximum axial velocity is 93, where it occurs at the tube axis close the inlet exit (Z/L=0.2). The results showed a good agreement for experimental and numerical analysis.

Keywords: counter flow, vortex tube, computational fluid dynamics analysis, energy separation, experimental study

Procedia PDF Downloads 79
3365 Numerical Simulation of the Flow around Wing-In-Ground Effect (WIG) Craft

Authors: A. Elbatran, Y. Ahmed, A. Radwan, M. Ishak

Abstract:

The use of WIG craft is representing an ambitious technology that will support in reducing time, effort, and money of the conventional marine transportation in the future. This paper investigates the aerodynamic characteristic of compound wing-in-ground effect (WIG) craft model. Drag coefficient, lift coefficient and Lift and drag ratio were studied numerically with respect to the ground clearance and the wing angle of attack. The modifications of the wing has been done in order to investigate the most suitable wing configuration that can increase the wing lift-to-drag ratio at low ground clearance. A numerical investigation was carried out in this research work using finite volume Reynolds-Averaged Navier-Stokes Equations (RANSE) code ANSYS CFX, Validation was carried out by using experiments. The experimental and the numerical results concluded that the lift to drag ratio decreased with the increasing of the ground clearance.

Keywords: drag Coefficient, ground clearance, navier-stokes, WIG

Procedia PDF Downloads 380
3364 Numerical Modeling of Turbulent Natural Convection in a Square Cavity

Authors: Mohammadreza Sedighi, Mohammad Said Saidi, Hesamoddin Salarian

Abstract:

A numerical study has been performed to investigate the effect of using different turbulent models on natural convection flow field and temperature distributions in partially heated square cavity compare to benchmark. The temperature of the right vertical wall is lower than that of heater while other walls are insulated. The commercial CFD codes are used to model. Standard k-w model provided good agreement with the experimental data.

Keywords: Buoyancy, Cavity, CFD, Heat Transfer, Natural Convection, Turbulence

Procedia PDF Downloads 341
3363 Numerical Simulation of Phase Transfer during Cryosurgery for an Irregular Tumor Using Hybrid Approach

Authors: Rama Bhargava

Abstract:

In the current paper, numerical simulation has been performed for the two-dimensional time dependent Pennes’ heat transfer model which is solved for irregular diseased tumor cells. An elliptic cryoprobe of varying sizes is taken at the center of the computational domain in such a manner that the location of the probe is fixed throughout the computation. The phase transition occurs due to the effect of probe with infusion of different nanoparticles Au, Al₂O₃, Fe₃O₄. The cooling performance of these nanoparticles injected at very low temperature, has been studied by implementing a hybrid FEM/EFGM method in which the whole domain is decomposed into two subdomains. The results are shown in terms of temperature profile inside the computational domain. Rate of cooling is obtained for various nanoparticles and it is observed that infusion of Au nanoparticles is very much efficient in increasing the heating rate than other nanoparticles. Such numerical scheme has direct applications where the domain is irregular.

Keywords: cryosurgery, hybrid EFGM/FEM, nanoparticles, simulation

Procedia PDF Downloads 242
3362 Study of the Electromagnetic Resonances of a Cavity with an Aperture Using Numerical Method and Equivalent Circuit Method

Authors: Ming-Chu Yin, Ping-An Du

Abstract:

The shielding ability of a shielding cavity is affected greatly by its resonances, which include resonance modes and frequencies. The equivalent circuit method and numerical method of transmission line matrix (TLM) are used to analyze the effect of aperture-cavity coupling on electromagnetic resonances of a cavity with an aperture in this paper. Both theoretical and numerical results show that the resonance modes of a shielding cavity with an aperture can be considered as the combination of cavity and aperture inherent resonance modes with resonance frequencies shifting, and the reason of this shift is aperture-cavity coupling. Because aperture sizes are important parameters to aperture-cavity coupling, variation rules of electromagnetic resonances of a shielding cavity with its aperture sizes are given, which will be useful for the design of shielding cavities.

Keywords: aperture-cavity coupling, equivalent circuit method, resonances, shielding equipment

Procedia PDF Downloads 444
3361 Performance Comparison and Visualization of COMSOL Multiphysics, Matlab, and Fortran for Predicting the Reservoir Pressure on Oil Production in a Multiple Leases Reservoir with Boundary Element Method

Authors: N. Alias, W. Z. W. Muhammad, M. N. M. Ibrahim, M. Mohamed, H. F. S. Saipol, U. N. Z. Ariffin, N. A. Zakaria, M. S. Z. Suardi

Abstract:

This paper presents the performance comparison of some computation software for solving the boundary element method (BEM). BEM formulation is the numerical technique and high potential for solving the advance mathematical modeling to predict the production of oil well in arbitrarily shaped based on multiple leases reservoir. The limitation of data validation for ensuring that a program meets the accuracy of the mathematical modeling is considered as the research motivation of this paper. Thus, based on this limitation, there are three steps involved to validate the accuracy of the oil production simulation process. In the first step, identify the mathematical modeling based on partial differential equation (PDE) with Poisson-elliptic type to perform the BEM discretization. In the second step, implement the simulation of the 2D BEM discretization using COMSOL Multiphysic and MATLAB programming languages. In the last step, analyze the numerical performance indicators for both programming languages by using the validation of Fortran programming. The performance comparisons of numerical analysis are investigated in terms of percentage error, comparison graph and 2D visualization of pressure on oil production of multiple leases reservoir. According to the performance comparison, the structured programming in Fortran programming is the alternative software for implementing the accurate numerical simulation of BEM. As a conclusion, high-level language for numerical computation and numerical performance evaluation are satisfied to prove that Fortran is well suited for capturing the visualization of the production of oil well in arbitrarily shaped.

Keywords: performance comparison, 2D visualization, COMSOL multiphysic, MATLAB, Fortran, modelling and simulation, boundary element method, reservoir pressure

Procedia PDF Downloads 491
3360 Buckling Analysis of Composite Shells under Compression and Torsional Loads: Numerical and Analytical Study

Authors: Güneş Aydın, Razi Kalantari Osgouei, Murat Emre Öztürk, Ahmad Partovi Meran, Ekrem Tüfekçi

Abstract:

Advanced lightweight laminated composite shells are increasingly being used in all types of modern structures, for enhancing their structural efficiency and performance. Such thin-walled structures are susceptible to buckling when subjected to various loading. This paper focuses on the buckling of cylindrical shells under axial compression and torsional loads. Effects of fiber orientation on the maximum buckling load of carbon fiber reinforced polymer (CFRP) shells are optimized. Optimum fiber angles have been calculated analytically by using MATLAB program. Numerical models have been carried out by using Finite Element Method program ABAQUS. Results from analytical and numerical analyses are also compared.

Keywords: buckling, composite, cylindrical shell, finite element, compression, torsion, MATLAB, optimization

Procedia PDF Downloads 587
3359 Numerical Modelling of Effective Diffusivity in Bone Tissue Engineering

Authors: Ayesha Sohail, Khadija Maqbool, Anila Asif, Haroon Ahmad

Abstract:

The field of tissue engineering is an active area of research. Bone tissue engineering helps to resolve the clinical problems of critical size and non-healing defects by the creation of man-made bone tissue. We will design and validate an efficient numerical model, which will simulate the effective diffusivity in bone tissue engineering. Our numerical model will be based on the finite element analysis of the diffusion-reaction equations. It will have the ability to optimize the diffusivity, even at multi-scale, with the variation of time. It will also have a special feature, with which we will not only be able to predict the oxygen, glucose and cell density dynamics, more accurately, but will also sort the issues arising due to anisotropy. We will fix these problems with the help of modifying the governing equations, by selecting appropriate spatio-temporal finite element schemes, by adaptive grid refinement strategy and by transient analysis.

Keywords: scaffolds, porosity, diffusion, transient analysis

Procedia PDF Downloads 541
3358 A Three-Dimensional (3D) Numerical Study of Roofs Shape Impact on Air Quality in Urban Street Canyons with Tree Planting

Authors: Bouabdellah Abed, Mohamed Bouzit, Lakhdar Bouarbi

Abstract:

The objective of this study is to investigate numerically the effect of roof shaped on wind flow and pollutant dispersion in a street canyon with one row of trees of pore volume, Pvol = 96%. A three-dimensional computational fluid dynamics (CFD) model for evaluating air flow and pollutant dispersion within an urban street canyon using Reynolds-averaged Navier–Stokes (RANS) equations and the k-Epsilon EARSM turbulence model as close of the equation system. The numerical model is performed with ANSYS-CFX code. Vehicle emissions were simulated as double line sources along the street. The numerical model was validated against the wind tunnel experiment. Having established this, the wind flow and pollutant dispersion in urban street canyons of six roof shapes are simulated. The numerical simulation agrees reasonably with the wind tunnel data. The results obtained in this work, indicate that the flow in 3D domain is more complicated, this complexity is increased with presence of tree and variability of the roof shapes. The results also indicated that the largest pollutant concentration level for two walls (leeward and windward wall) is observed with the upwind wedge-shaped roof. But the smallest pollutant concentration level is observed with the dome roof-shaped. The results also indicated that the corners eddies provide additional ventilation and lead to lower traffic pollutant concentrations at the street canyon ends.

Keywords: street canyon, pollutant dispersion, trees, building configuration, numerical simulation, k-Epsilon EARSM

Procedia PDF Downloads 365
3357 Numerical Modeling to Validate Theoretical Models of Toppling Failure in Rock Slopes

Authors: Hooman Dabirmanesh, Attila M. Zsaki

Abstract:

Traditionally, rock slope stability is carried out using limit equilibrium analysis when investigating toppling failure. In these equilibrium methods, internal forces exerted between columns are not clearly defined, and to the authors’ best knowledge, there is no consensus in literature with respect to the results of analysis. A discrete element method-based numerical model was developed and applied to simulate the behavior of rock layers subjected to toppling failure. Based on this calibrated numerical model, a study of the location and distribution of internal forces that result in equilibrium was carried out. The sum of side forces was applied at a point on a block which properly represents the force to determine the inter-column force distribution. In terms of the side force distribution coefficient, the result was compared to those obtained from laboratory centrifuge tests. The results of the simulation show the suitable criteria to select the correct position for the internal exerted force between rock layers. In addition, the numerical method demonstrates how a theoretical method could be reliable by considering the interaction between the rock layers.

Keywords: contact bond, discrete element, force distribution, limit equilibrium, tensile stress

Procedia PDF Downloads 143
3356 Numerical Investigation of Flow and Heat Transfer Characteristics of a Natural Refrigerant within a Vortex Tube

Authors: Mirza Popovac

Abstract:

This paper investigates the application of the vortex tubes towards increasing the efficiency of high temperature heat pumps based on natural refrigerants, by recovering a part of the expansion work within the refrigerant cycle. To this purpose the 3D Navier-Stokes solver is used to perform a set of numerical simulations, investigating the vortex tube performance. Firstly, the fluid flow and heat transfer characteristics are analyzed for standard configurations of vortex tubes, and the obtained results are validated against the experimental and numerical data available in literature. Subsequently, different geometry specifications are analyzed, as well as the interplay between relevant heat pump operating conditions and the properties of natural refrigerants. Finally, the characteristic curve of performance will be derived for investigated vortex tubes specifications when used within high temperature heat pumps.

Keywords: heat pump, vortex tube, CFD, natural refrigerant

Procedia PDF Downloads 140
3355 An Entropy Stable Three Dimensional Ideal MHD Solver with Guaranteed Positive Pressure

Authors: Andrew R. Winters, Gregor J. Gassner

Abstract:

A high-order numerical magentohydrodynamics (MHD) solver built upon a non-linear entropy stable numerical flux function that supports eight traveling wave solutions will be described. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver is especially well-suited for flows involving strong discontinuities due to its strong stability without the need to enforce artificial low density or energy limits. Furthermore, a new formulation of the numerical algorithm to guarantee positivity of the pressure during the simulation is described and presented. By construction, the solver conserves mass, momentum, and energy and is entropy stable. High spatial order is obtained through the use of a third order limiting technique. High temporal order is achieved by utilizing the family of strong stability preserving (SSP) Runge-Kutta methods. Main attributes of the solver are presented as well as details on an implementation of the new solver into the multi-physics, multi-scale simulation code FLASH. The accuracy, robustness, and computational efficiency is demonstrated with a variety of numerical tests. Comparisons are also made between the new solver and existing methods already present in FLASH framework.

Keywords: entropy stability, finite volume scheme, magnetohydrodynamics, pressure positivity

Procedia PDF Downloads 343
3354 Experimental and Numerical Processes of Open Die Forging of Multimetallic Materials with the Usage of Different Lubricants

Authors: Isik Cetintav, Cenk Misirli, Yilmaz Can, Damla Gunel

Abstract:

This work investigates experimental and numerical analysis of open die forging of multimetallic materials. Multimetallic material production has recently become an interesting research field. The mechanical properties of the materials to be used for the formation of multimetallic materials and the mechanical properties of the multimetallic materials produced will be compared and the material flows of the use of different lubricants will be examined. Furthermore, in this work, the mechanical properties of multimetallic metallic materials produced using different materials will be examined by using different lubricants. The advantages and disadvantages of different lubricants will be approached with the bi-metallic material to be produced. Cylindrical specimens consisting of two different materials were used in the experiments. Specimens were prepared as aluminum sleeve and copper core and upset at different reduction. This metal combination present a material model of which chemical composition is different. ABAQUS software was used for the simulations. Simulation and experimental results have also shown reasonable agreement.

Keywords: multimetallic, forging, experimental, numerical

Procedia PDF Downloads 278
3353 Development of Numerical Method for Mass Transfer across the Moving Membrane with Selective Permeability: Approximation of the Membrane Shape by Level Set Method for Numerical Integral

Authors: Suguru Miyauchi, Toshiyuki Hayase

Abstract:

Biological membranes have selective permeability, and the capsules or cells enclosed by the membrane show the deformation by the osmotic flow. This mass transport phenomenon is observed everywhere in a living body. For the understanding of the mass transfer in a body, it is necessary to consider the mass transfer phenomenon across the membrane as well as the deformation of the membrane by a flow. To our knowledge, in the numerical analysis, the method for mass transfer across the moving membrane has not been established due to the difficulty of the treating of the mass flux permeating through the moving membrane with selective permeability. In the existing methods for the mass transfer across the membrane, the approximate delta function is used to communicate the quantities on the interface. The methods can reproduce the permeation of the solute, but cannot reproduce the non-permeation. Moreover, the computational accuracy decreases with decreasing of the permeable coefficient of the membrane. This study aims to develop the numerical method capable of treating three-dimensional problems of mass transfer across the moving flexible membrane. One of the authors developed the numerical method with high accuracy based on the finite element method. This method can capture the discontinuity on the membrane sharply due to the consideration of the jumps in concentration and concentration gradient in the finite element discretization. The formulation of the method takes into account the membrane movement, and both permeable and non-permeable membranes can be treated. However, searching the cross points of the membrane and fluid element boundaries and splitting the fluid element into sub-elements are needed for the numerical integral. Therefore, cumbersome operation is required for a three-dimensional problem. In this paper, we proposed an improved method to avoid the search and split operations, and confirmed its effectiveness. The membrane shape was treated implicitly by introducing the level set function. As the construction of the level set function, the membrane shape in one fluid element was expressed by the shape function of the finite element method. By the numerical experiment, it was found that the shape function with third order appropriately reproduces the membrane shapes. The same level of accuracy compared with the previous method using search and split operations was achieved by using a number of sampling points of the numerical integral. The effectiveness of the method was confirmed by solving several model problems.

Keywords: finite element method, level set method, mass transfer, membrane permeability

Procedia PDF Downloads 250
3352 Numerical Study of Natural Convection in a Triangular Enclosure as an Attic for Different Geometries and Boundary Conditions

Authors: H. Golchoobian, S. Saedodin, M. H. Taheri, A. Sarafraz

Abstract:

In this paper, natural convection in an attic is numerically investigated. The geometry of the problem is considered to be a triangular enclosure. ANSYS Fluent software is used for modeling and numerical solution. This study is for steady state. Four right-angled triangles with height to base ratios of 2, 1, 0.5 and 0.25 are considered. The behavior of various parameters related to its performance, including temperature distribution and velocity vectors are evaluated, and graphs for the Nusselt number have been drawn. Also, in this study, the effect of geometric shape of enclosure with different height-to-base ratios has been evaluated for three types of boundary conditions of winter, summer day and one another state. It can be concluded that as the bottom side temperature and ratio of base to height of the enclosure increases, the convective effects become more prominent and circulation happened.

Keywords: enclosure, natural convection, numerical solution, Nusselt number, triangular

Procedia PDF Downloads 197
3351 A Numerical Solution Based on Operational Matrix of Differentiation of Shifted Second Kind Chebyshev Wavelets for a Stefan Problem

Authors: Rajeev, N. K. Raigar

Abstract:

In this study, one dimensional phase change problem (a Stefan problem) is considered and a numerical solution of this problem is discussed. First, we use similarity transformation to convert the governing equations into ordinary differential equations with its boundary conditions. The solutions of ordinary differential equation with the associated boundary conditions and interface condition (Stefan condition) are obtained by using a numerical approach based on operational matrix of differentiation of shifted second kind Chebyshev wavelets. The obtained results are compared with existing exact solution which is sufficiently accurate.

Keywords: operational matrix of differentiation, similarity transformation, shifted second kind chebyshev wavelets, stefan problem

Procedia PDF Downloads 402
3350 Biases in Numerically Invariant Joint Signatures

Authors: Reza Aghayan

Abstract:

This paper illustrates that numerically invariant joint signatures suffer biases in the resulting signatures. Next, we classify the arising biases as Bias Type 1 and Bias Type 2 and show how they can be removed.

Keywords: Euclidean and affine geometries, differential invariant signature curves, numerically invariant joint signatures, numerical analysis, numerical bias, curve analysis

Procedia PDF Downloads 597
3349 Large Eddy Simulation of Hydrogen Deflagration in Open Space and Vented Enclosure

Authors: T. Nozu, K. Hibi, T. Nishiie

Abstract:

This paper discusses the applicability of the numerical model for a damage prediction method of the accidental hydrogen explosion occurring in a hydrogen facility. The numerical model was based on an unstructured finite volume method (FVM) code “NuFD/FrontFlowRed”. For simulating unsteady turbulent combustion of leaked hydrogen gas, a combination of Large Eddy Simulation (LES) and a combustion model were used. The combustion model was based on a two scalar flamelet approach, where a G-equation model and a conserved scalar model expressed a propagation of premixed flame surface and a diffusion combustion process, respectively. For validation of this numerical model, we have simulated the previous two types of hydrogen explosion tests. One is open-space explosion test, and the source was a prismatic 5.27 m3 volume with 30% of hydrogen-air mixture. A reinforced concrete wall was set 4 m away from the front surface of the source. The source was ignited at the bottom center by a spark. The other is vented enclosure explosion test, and the chamber was 4.6 m × 4.6 m × 3.0 m with a vent opening on one side. Vent area of 5.4 m2 was used. Test was performed with ignition at the center of the wall opposite the vent. Hydrogen-air mixtures with hydrogen concentrations close to 18% vol. were used in the tests. The results from the numerical simulations are compared with the previous experimental data for the accuracy of the numerical model, and we have verified that the simulated overpressures and flame time-of-arrival data were in good agreement with the results of the previous two explosion tests.

Keywords: deflagration, large eddy simulation, turbulent combustion, vented enclosure

Procedia PDF Downloads 244
3348 A Comparative Study between FEM and Meshless Methods

Authors: Jay N. Vyas, Sachin Daxini

Abstract:

Numerical simulation techniques are widely used now in product development and testing instead of expensive, time-consuming and sometimes dangerous laboratory experiments. Numerous numerical methods are available for performing simulation of physical problems of different engineering fields. Grid based methods, like Finite Element Method, are extensively used in performing various kinds of static, dynamic, structural and non-structural analysis during product development phase. Drawbacks of grid based methods in terms of discontinuous secondary field variable, dealing fracture mechanics and large deformation problems led to development of a relatively a new class of numerical simulation techniques in last few years, which are popular as Meshless methods or Meshfree Methods. Meshless Methods are expected to be more adaptive and flexible than Finite Element Method because domain descretization in Meshless Method requires only nodes. Present paper introduces Meshless Methods and differentiates it with Finite Element Method in terms of following aspects: Shape functions used, role of weight function, techniques to impose essential boundary conditions, integration techniques for discrete system equations, convergence rate, accuracy of solution and computational effort. Capabilities, benefits and limitations of Meshless Methods are discussed and concluded at the end of paper.

Keywords: numerical simulation, Grid-based methods, Finite Element Method, Meshless Methods

Procedia PDF Downloads 389
3347 Control of a Plane Jet Spread by Tabs at the Nozzle Exit

Authors: Makito Sakai, Takahiro Kiwata, Takumi Awa, Hiroshi Teramoto, Takaaki Kono, Kuniaki Toyoda

Abstract:

Using experimental and numerical results, this paper describes the effects of tabs on the flow characteristics of a plane jet at comparatively low Reynolds numbers while focusing on the velocity field and the vortical structure. The flow visualization and velocity measurements were respectively carried out using laser Doppler velocimetry (LDV) and particle image velocimetry (PIV). In addition, three-dimensional (3D) plane jet numerical simulations were performed using ANSYS Fluent, a commercially available computational fluid dynamics (CFD) software application. We found that the spreads of jets perturbed by large delta tabs and round tabs were larger than those produced by the other tabs tested. Additionally, it was determined that a plane jet with square tabs had the smallest jet spread downstream, and the jet’s centerline velocity was larger than those of jets perturbed by the other tabs tested. It was also observed that the spanwise vortical structure of a plane jet with tabs disappeared completely. Good agreement was found between the experimental and numerical simulation velocity profiles in the area near the nozzle exit when the laminar flow model was used. However, we also found that large eddy simulation (LES) is better at predicting the developing flow field of a plane jet than the laminar and the standard k-ε turbulent models.

Keywords: plane jet, flow control, tab, flow measurement, numerical simulation

Procedia PDF Downloads 334
3346 Geomechanical Technologies for Assessing Three-Dimensional Stability of Underground Excavations Utilizing Remote-Sensing, Finite Element Analysis, and Scientific Visualization

Authors: Kwang Chun, John Kemeny

Abstract:

Light detection and ranging (LiDAR) has been a prevalent remote-sensing technology applied in the geological fields due to its high precision and ease of use. One of the major applications is to use the detailed geometrical information of underground structures as a basis for the generation of a three-dimensional numerical model that can be used in a geotechnical stability analysis such as FEM or DEM. To date, however, straightforward techniques in reconstructing the numerical model from the scanned data of the underground structures have not been well established or tested. In this paper, we propose a comprehensive approach integrating all the various processes, from LiDAR scanning to finite element numerical analysis. The study focuses on converting LiDAR 3D point clouds of geologic structures containing complex surface geometries into a finite element model. This methodology has been applied to Kartchner Caverns in Arizona, where detailed underground and surface point clouds can be used for the analysis of underground stability. Numerical simulations were performed using the finite element code Abaqus and presented by 3D computing visualization solution, ParaView. The results are useful in studying the stability of all types of underground excavations including underground mining and tunneling.

Keywords: finite element analysis, LiDAR, remote-sensing, scientific visualization, underground stability

Procedia PDF Downloads 173
3345 The Development of a New Block Method for Solving Stiff ODEs

Authors: Khairil I. Othman, Mahfuzah Mahayaddin, Zarina Bibi Ibrahim

Abstract:

We develop and demonstrate a computationally efficient numerical technique to solve first order stiff differential equations. This technique is based on block method whereby three approximate points are calculated. The Cholistani of varied step sizes are presented in divided difference form. Stability regions of the formulae are briefly discussed in this paper. Numerical results show that this block method perform very well compared to existing methods.

Keywords: block method, divided difference, stiff, computational

Procedia PDF Downloads 428
3344 Seismic Retrofit of Rectangular Columns Using Fiber Reinforced Polymers

Authors: E. L. Elghazy, A. M. Sanad, M. G. Ghoneim

Abstract:

Over the past two decades research has shown that fiber reinforced polymers can be efficiently, economically and safely used for strengthening and rehabilitation of reinforced concrete (RC) structures. Designing FRP confined concrete columns requires reliable analytical tools that predict the level of performance and ductility enhancement. A numerical procedure is developed aiming at determining the type and thickness of FRP jacket needed to achieve a certain level of ductility enhancement. The procedure starts with defining the stress strain curve, which is used to obtain moment curvature relationship then displacement ductility ratio of reinforced concrete cross-sections subjected to bending moment and axial force. Three sets of published experimental tests were used to validate the numerical procedure. Comparisons between predicted results obtained by using the proposed procedure and actual results of experimental tests proved the reliability of the proposed procedure.

Keywords: columns, confinement, ductility, FRP, numerical

Procedia PDF Downloads 448
3343 Numerical Investigation of Pressure and Velocity Field Contours of Dynamics of Drop Formation

Authors: Pardeep Bishnoi, Mayank Srivastava, Mrityunjay Kumar Sinha

Abstract:

This article represents the numerical investigation of the pressure and velocity field variation of the dynamics of pendant drop formation through a capillary tube. Numerical simulations are executed using volume of fluid (VOF) method in the computational fluid dynamics (CFD). In this problem, Non Newtonian fluid is considered as dispersed fluid whereas air is considered as a continuous fluid. Pressure contours at various time steps expose that pressure varies nearly hydrostatically at each step of the dynamics of drop formation. A result also shows the pressure variation of the liquid droplet during free fall in the computational domain. The evacuation of the fluid from the necking region is also shown by the contour of the velocity field. The role of surface tension in the Pressure contour of the dynamics of drop formation is also studied.

Keywords: pressure contour, surface tension, volume of fluid, velocity field

Procedia PDF Downloads 405
3342 Numerical Study on Jatropha Oil Pool Fire Behavior in a Compartment

Authors: Avinash Chaudhary, Akhilesh Gupta, Surendra Kumar, Ravi Kumar

Abstract:

This paper presents the numerical study on Jatropha oil pool fire in a compartment. A fire experiment with jatropha oil was conducted in a compartment of size 4 m x 4 m x m to study the fire development and temperature distribution. Fuel is burned in the center of the compartment in a pool diameter of 0.5 m with an initial fuel depth of 0.045 m. Corner temperature in the compartment, doorway temperature and hot gas layer temperature at various locations are measured. Numerical simulations were carried out using Fire Dynamics Simulator (FDS) software at grid size of 0.05 m, 0.12 m and for performing simulation heat release rate of jatropha oil measured using mass loss method were inputted into FDS. Experimental results shows that like other fuel fires, the whole combustion process can be divided into four stages: initial stage, growth stage, steady profile or developed phase and decay stage. The fire behavior shows two zone profile where upper zone consists of mainly hot gases while lower zone is relatively at colder side. In this study, predicted temperatures from simulation are in good agreement in upper zone of compartment. Near the interface of hot and cold zone, deviations were reported between the simulated and experimental results which is probably due to the difference between the predictions of smoke layer height by FDS. Also, changing the grid size from 0.12 m to 0.05 m does not show any effect in temperatures at upper zone while in lower zone, grid size of 0.05 m showed satisfactory agreement with experimental results. Numerical results showed that calculated temperatures at various locations matched well with the experimental results. On the whole, an effective method is provided with reasonable results to study the burning characteristics of jatropha oil with numerical simulations.

Keywords: jatropha oil, compartment fire, heat release rate, FDS (fire dynamics simulator), numerical simulation

Procedia PDF Downloads 257
3341 Numerical Study of Fluid Flow and Heat Transfer in the Spongy-Porous Media

Authors: Zeinab Sayed Abdel Rehim, M. A. Ziada, H. Salwa El-Deeb

Abstract:

Numerical study of fluid flow, heat transfer and thermal energy storing or released in/from spongy-porous media to predict the thermal performance and characteristics of the porous media as packed bed system is presented in this work. This system is cylindrical channel filled with porous media (carbon foam). The system consists of working fluid (air) and spongy-porous medium; they act as the heat exchanger (heating or cooling modes) where thermal interaction occurs between the working fluid and the porous medium. The spongy-porous media are defined by the different type of porous medium employed in the storing or cooling modes. Two different porous media are considered in this study: Carbon foam, and Silicon rubber. The flow of the working fluid (air) is one dimensional in the axial direction from the top to downward and steady state conditions. The numerical results of transient temperature distribution for both working fluid and the spongy-porous medium phases and the amount of stored/realized heat inside/from the porous medium for each case with respect to the operating parameters and the spongy-porous media characteristics are illustrated.

Keywords: fluid flow, heat transfer, numerical analysis, spongy-porous media, thermal performance, transient conditions

Procedia PDF Downloads 545