Interaction between a notch and cylindrical voids in aluminum single crystals: Experimental observations and numerical simulations

A combined 3D finite element simulation and experimental study of interaction between a notch and cylindrical voids ahead of it in single edge notch (tension) aluminum single crystal specimens is undertaken in this work. Two lattice orientations are considered in which the notch front is parallel to the crystallographic 10 (1) over bar] direction. The flat surface of the notch coincides with the (010) plane in one orientation and with the (1 (1) over bar1) plane in the other. Three equally spaced cylindrical voids are placed directly ahead of the notch tip. The predicted load-displacement curves, slip traces, lattice rotation and void growth from the finite element analysis are found to be in good agreement with the experimental observations for both the orientations. Finite element results show considerable through-thickness variation in both hydrostatic stress and equivalent plastic slip which, however, depends additionally on the lattice orientation. The through-thickness variation in the above quantities affects the void growth rate and causes it to differ from the center-plane to the free surface of the specimen. (c) 2012 Elsevier Ltd. All rights reserved.

Tiling stencil computations to maximize parallelism

Most stencil computations allow tile-wise concurrent start, i.e., there always exists a face of the iteration space and a set of tiling hyperplanes such that all tiles along that face can be started concurrently. This provides load balance and maximizes parallelism. However, existing automatic tiling frameworks often choose hyperplanes that lead to pipelined start-up and load imbalance. We address this issue with a new tiling technique that ensures concurrent start-up as well as perfect load-balance whenever possible. We first provide necessary and sufficient conditions on tiling hyperplanes to enable concurrent start for programs with affine data accesses. We then provide an approach to find such hyperplanes. Experimental evaluation on a 12-core Intel Westmere shows that our code is able to outperform a tuned domain-specific stencil code generator by 4% to 27%, and previous compiler techniques by a factor of 2x to 10.14x.

Numerical simulations of CO2 migration during charnockite genesis

Charnockite is considered to be generated either through the dehydration of granitic magma by CO2 purging or by solid-state dehydration through CO2 metasomatism during granulite facies metamorphism. To understand the extent of dehydration, CO2 migration is quantitatively modeled in silicate melt and metasomatic fluid as a function of temperature, H2O wt%, pressure, basal CO2 flux and dynamic viscosity. Numerical simulations show that CO2 advection through porous and permeable high-grade metamorphic rocks can generate dehydrated patches close to the CO2 flow path, as illustrated by the occurrences of ``incipient charnockites.'' CO2 reaction-front velocity constrained by field observations is 0.69 km/m.y., a reasonable value, which matches well with other studies. On the other hand, temperature, rate of cooling, and basal CO2 flux are the critical parameters affecting CO2 diffusion through a silicate melt. CO2 diffusion through silicate melt can only occur at temperature greater than 840 degrees C and during slow cooling (<= 3.7 x 10(-5) degrees C/yr), features that are typical of magma emplacement in the lower crust. Stalling of CO2 fluxing at similar to 840 degrees C explains why some deep-level plutons contain both hydrous and anhydrous (charnockitic) mineral assemblages. CO2 diffusion through silicate melt is virtually insensitive to pressure. Addition of CO2 basal flux facilitates episodic dehydrated melt migration by generating fracture pathways.

Numerical analysis of geocell reinforced ballast overlying soft clay subgrade

Geotextiles and geogrids have been in use for several decades in variety of geo-structure applications including foundation of embankments, retaining walls, pavements. Geocells is one such variant in geosynthetic reinforcement of recent years, which provides a three dimensional confinement to the infill material. Although extensive research has been carried on geocell reinforced sand, clay and layered soil subgrades, limited research has been reported on the aggregates/ballast reinforced with geocells. This paper presents the behavior of a railway sleeper subjected to monotonic loading on geocell reinforced aggregates, of size ranging from 20 to 75 mm, overlying soft clay subgrades. Series of tests were conducted in a steel test tank of dimensions 700 mm x 300 mm x 700 mm. In addition to the laboratory model tests, numerical simulations were performed using a finite difference code to predict the behavior of geocell reinforced ballast. The results from numerical simulations were compared with the experimental data. The numerical and experimental results manifested the importance that the geocell reinforcement has a significant effect on the ballast behaviour. The results depicted that the stiffness of underlying soft clay subgrade has a significant influence on the behavior of the geocell-aggregate composite material in redistributing the loading system.

Vorticity moments in four numerical simulations of the 3D Navier-Stokes equations

The issue of intermittency in numerical solutions of the 3D Navier-Stokes equations on a periodic box 0, L](3) is addressed through four sets of numerical simulations that calculate a new set of variables defined by D-m(t) = (pi(-1)(0) Omega(m))(alpha m) for 1 <= m <= infinity where alpha(m) = 2m/(4m - 3) and Omega(m)(t)](2m) = L-3 integral(v) vertical bar omega vertical bar(2m) dV with pi(0) = vL(-2). All four simulations unexpectedly show that the D-m are ordered for m = 1,..., 9 such that Dm+1 < D-m. Moreover, the D-m squeeze together such that Dm+1/D-m NE arrow 1 as m increases. The values of D-1 lie far above the values of the rest of the D-m, giving rise to a suggestion that a depletion of nonlinearity is occurring which could be the cause of Navier-Stokes regularity. The first simulation is of very anisotropic decaying turbulence; the second and third are of decaying isotropic turbulence from random initial conditions and forced isotropic turbulence at fixed Grashof number respectively; the fourth is of very-high-Reynolds-number forced, stationary, isotropic turbulence at up to resolutions of 4096(3).

Numerical simulations of bubble formation from submerged needles under non-uniform direct current electric field

In several chemical and space industries, small bubbles are desired for efficient interaction between the liquid and gas phases. In the present study, we show that non-uniform electric field with appropriate electrode configurations can reduce the volume of the bubbles forming at submerged needles by up to three orders of magnitude. We show that localized high electric stresses at the base of the bubbles result in slipping of the contact line on the inner surface of the needle and subsequent bubble formation occurs with contact line inside the needle. We also show that for bubble formation in the presence of highly non-uniform electric field, due to high detachment frequency, the bubbles go through multiple coalescences and thus increase the apparent volume of the detached bubbles. (C) 2013 AIP Publishing LLC.

Thermal Efficiency of LPG and PNG-fired burners: Experimental and numerical studies

In the present work, the thermal efficiency of a conventional domestic burner is studied both experimentally and numerically for liquefied petroleum gas (LPG) and piped natural gas (PNG) fuels. Three-dimensional computational fluid dynamic (CFD) modeling of the steady-state flow, combustion and heat transfer to the vessel is reported for the first time in such burners. Based on the insights from the CFD model concerning the flow and heat transfer, design modifications in the form of a circular insert and a radiant sheet are proposed which are observed to increase thermal efficiency for LPG. For PNG, predictions showed that loading height was a much more important factor affecting efficiency than these design modifications and an optimal loading height could be identified. Experiments confirm these trends by showing an improvement in burner thermal efficiency of 2.5% for LPG with the modified design, and 10% for PNG with the optimal loading height, demonstrating that the CFD modeling approach developed in the present work is a useful tool to study domestic burners. (C) 2013 Elsevier Ltd. All rights reserved.

Numerical and experimental investigation and optimization of an open wheeled race car cooling air duct

In this study the cooling performance due to air flow and aerodynamics of the Formula Student open wheeled race car has been investigated and optimized with the help of CFD simulations and experimental validation. The race car in context previously suffered from overheating problems. Flow analysis was carried out based on the detailed race car 3D model (NITK Racing 2012 formula student race car). Wind tunnel experiments were carried out on the same. The results obtained from the computer simulations are compared with experimental results obtained from wind tunnel testing of the full car. Through this study it was possible to locate the problem areas and hence choose the best configuration for the cooling duct. The CFD analysis helped in calculating the mass flow rate, pressure and velocity distribution for different velocities of the car which is then used to determine the heat dissipated by the radiator. Area of flow separation could be visualized and made sure smooth airflow into the radiator core area. This significantly increased the cooling performance of the car with reduction in drag.

Effect of Inlet Jet injection angle on a calandria based reactor-An investigation with numerical analysis of heat and mass transfer for an optimum reactor design

Heat and mass transfer studies in a calandria based reactor is quite complex both due to geometry and due to the complex mixing flow. It is challenging to devise optimum operating conditions with efficient but safe working range for such a complex configuration. Numerical study known to be very effective is taken up for investigation. In the present study a 3D RANS code with turbulence model has been used to compute the flow fields and to get the heat transfer characteristics to understand certain design parameters of engineering importance. The angle of injection and of the coolant liquid has a large effect on the heat transfer within the reactor.

Numerical study of the impact of inoculant and grain transport on macrosegregation and microstructure formation during solidification of an Al-22%Cu alloy

We investigate the impact of the nucleation law for nucleation on Al-Ti-B inoculant particles, of the motion of inoculant particles and of the motion of grains on the predicted macrosegregation and microstructure in a grain-refined Al-22 wt.% Cu alloy casting. We conduct the study by numerical simulations of a casting experiment in a side-cooled 76×76×254 mm sand mould. Macrosegregation and microstructure formation are studied with a volume-averaged two-phase model accounting for macroscopic heat and solute transport, melt convection, and transport of inoculant particles and equiaxed grains. On the microscopic scale it accounts for nucleation on inoculant particles with a given size distribution (and corresponding activation undercooling distribution)and for the growth of globular solid grains. The growth kinetics is described by accounting for limited solute diffusion in both liquid and solid phases and for convective effects. We show that the consideration of a size distribution of the inoculants has a strong impact on the microstructure(final grain size) prediction. The transport of inoculants significantly increases the microstructure heterogeneities and the grain motion refines the microstructure and reduces the microstructure heterogeneities.

An experimental and numerical study of piezoceramic actuator hysteresis in helicopter active vibration control

An aeroelastic analysis is used to investigate the rate dependent hysteresis in piezoceramic actuators and its effect on helicopter vibration control with trailing edge flaps. Hysteresis in piezoceramic materials can cause considerable complications in the use of smart actuators as prime movers in applications such as helicopter active vibration control. Dynamic hysteresis of the piezoelectric stack actuator is investigated for a range of frequencies (5 Hz (1/rev) to 30 Hz (6/rev)) which are of practical importance for helicopter vibration analysis. Bench top tests are conducted on a commercially available piezoelectric stack actuator. Frequency dependent hysteretic behavior is studied experimentally for helicopter operational frequencies. Material hysteresis in the smart actuator is mathematically modeled using the theory of conic sections. Numerical simulations are also performed at an advance ratio of 0.3 for vibration control analysis using a trailing edge flap with an idealized linear and a hysteretic actuator. The results indicate that dynamic hysteresis has a notable effect on the hub vibration levels. It is found that the theory of conic sections offers a straight forward approach for including hysteresis into aeroelastic analysis.

Analysis of composite single lap joints using numerical and experimental approach

Adhesives are widely used to execute the assembly of aerospace and automotive structures due to their ability to join dissimilar materials, reduced stress concentration, and improved fatigue resistance. The mechanical behavior of adhesive joints can be studied either using analytical models or by conducting mechanical tests. However, the complexity owing to multiple interfaces, layers with different properties, material and geometric nonlinearity and its three-dimensional nature combine to increase the difficulty in obtaining an overall system of governing equations to predict the joint behavior. On the other hand, experiments are often time consuming and expensive due to a number of parameters involved. Finite element analysis (FEA) is profoundly used in recent years to overcome these limitations. The work presented in this paper involves the finite element modeling and analysis of a composite single lap joint where the adhesive-adherend interface region was modeled using connector elements. The computed stresses were compared with the experimental stresses obtained using digital image correlation technique. The results showed an agreement. Further, the failure load predicted using FEA was found to be closer to the actual failure load obtained by mechanical tests.

Effect of angle of incidence on mixed convective wake dynamics and heat transfer past a square cylinder in cross flow at Re=100

In this paper, a numerical investigation is performed to study the mixed convective flow and heat transfer characteristics past a square cylinder in cross flow at incidence. Utilizing air (Pr = 0.71) as an operating fluid, computations are carried out at a representative Reynolds number (Re) of 100. Angles of incidences are varied as, 0 degrees <= alpha <= 45 degrees. Effect of superimposed positive and negative cross-flow buoyancy is brought about by varying the Richardson number (RI) in the range -1.0 <= Ri <= 1.0. The detail features of flow topology and heat transport are analyzed critically for different angles of incidences. The thermo fluidic forces acting on the cylinder during mixed convection are captured in terms of the drag (C-D), lift (C-L), and moment (C-M) coefficients. The results show that the lateral width of the cylinder wake reduces with increasing alpha and the isotherms spread out far wide. In the range 0 degrees < alpha < 45 degrees, C-D reduces with increasing Ri. The functional dependence of C-M with Ri reveals a linear relationship. Thermal boundary layer thickness reduces with increasing angle of incidences. The global rate of heat transfer from the cylinder increases with increasing alpha. (C) 2014 Elsevier Ltd. All rights reserved.

Experimental and numerical studies on friction welding of thixocast A356 aluminum alloy

This paper highlights the role of globular microstructure on the weldability of semi-solid processed aluminum alloys via high temperature flow behavior. The investigation was carried out on the joining of thixocast A356 aluminum alloy components by friction welding. A thermomechanical model was developed to predict the temperature and stress distributions, as well as to identify the suitable and safe range of parameters. Good comparisons between numerical and experimental results were observed. In addition, metallographic examinations and hardness and tensile tests of the welded samples were carried out. It was found that the tensile strength of the joint is higher than the tensile strength of the parent material for the optimum set of parameters. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.