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.

Scalable sequential alternating proximal methods for sparse structural SVMs and CRFs

Structural Support Vector Machines (SSVMs) and Conditional Random Fields (CRFs) are popular discriminative methods used for classifying structured and complex objects like parse trees, image segments and part-of-speech tags. The datasets involved are very large dimensional, and the models designed using typical training algorithms for SSVMs and CRFs are non-sparse. This non-sparse nature of models results in slow inference. Thus, there is a need to devise new algorithms for sparse SSVM and CRF classifier design. Use of elastic net and L1-regularizer has already been explored for solving primal CRF and SSVM problems, respectively, to design sparse classifiers. In this work, we focus on dual elastic net regularized SSVM and CRF. By exploiting the weakly coupled structure of these convex programming problems, we propose a new sequential alternating proximal (SAP) algorithm to solve these dual problems. This algorithm works by sequentially visiting each training set example and solving a simple subproblem restricted to a small subset of variables associated with that example. Numerical experiments on various benchmark sequence labeling datasets demonstrate that the proposed algorithm scales well. Further, the classifiers designed are sparser than those designed by solving the respective primal problems and demonstrate comparable generalization performance. Thus, the proposed SAP algorithm is a useful alternative for sparse SSVM and CRF classifier design.

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.

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.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Efficient simulation of unitary operators by combining two numerical algorithms: An NMR simulation of the mirror-inversion propagator of an XY spin chain

Precise experimental implementation of unitary operators is one of the most important tasks for quantum information processing. Numerical optimization techniques are widely used to find optimized control fields to realize a desired unitary operator. However, finding high-fidelity control pulses to realize an arbitrary unitary operator in larger spin systems is still a difficult task. In this work, we demonstrate that a combination of the GRAPE algorithm, which is a numerical pulse optimization technique, and a unitary operator decomposition algorithm Ajoy et al., Phys. Rev. A 85, 030303 (2012)] can realize unitary operators with high experimental fidelity. This is illustrated by simulating the mirror-inversion propagator of an XY spin chain in a five-spin dipolar coupled nuclear spin system. Further, this simulation has been used to demonstrate the transfer of entangled states from one end of the spin chain to the other end.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

A detailed numerical study of NOx kinetics in low calorific value H-2/CO syngas flames

The present study provides an extensive and detailed numerical analysis of NO chemical kinetics in low calorific value H-2/CO syngas flames utilizing predictions by five chemical kinetic mechanisms available out of which four deal with H-2/CO while the fifth mechanism (GRI 3.0) additionally accounts for hydrocarbon chemistry. Comparison of predicted axial NO profiles in premixed flat flames with measurements at 1 bar, 3.05 bar and 9.15 bar shows considerably large quantitative differences among the various mechanisms. However, at each pressure, the quantitative reaction path diagrams show similar NO formation pathways for most of the mechanisms. Interestingly, in counterflow diffusion flames, the quantitative reaction path diagrams and sensitivity analyses using the various mechanisms reveal major differences in the NO formation pathways and reaction rates of important reactions. The NNH and N2O intermediate pathways are found to be the major contributors for NO formation in all the reaction mechanisms except GRI 3.0 in syngas diffusion flames. The GRI 3.0 mechanism is observed to predict prompt NO pathway as the major contributing pathway to NO formation. This is attributed to prediction of a large concentration of CH radical by the GRI 3.0 as opposed to a relatively negligible value predicted by all other mechanisms. Also, the back-conversion of NNH into N2O at lower pressures (2-4 bar) was uniquely observed for one of the five mechanisms. The net reaction rates and peak flame temperatures are used to correlate and explain the differences observed in the peak NO] at different pressures. This study identifies key reactions needing assessment and also highlights the need for experimental data in syngas diffusion flames in order to assess and optimize H-2/CO and nitrogen chemistry. Copyright (C) 2014, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Compiler/Runtime Framework for Dynamic Dataflow Parallelization of Tiled Programs

Task-parallel languages are increasingly popular. Many of them provide expressive mechanisms for intertask synchronization. For example, OpenMP 4.0 will integrate data-driven execution semantics derived from the StarSs research language. Compared to the more restrictive data-parallel and fork-join concurrency models, the advanced features being introduced into task-parallelmodels in turn enable improved scalability through load balancing, memory latency hiding, mitigation of the pressure on memory bandwidth, and, as a side effect, reduced power consumption. In this article, we develop a systematic approach to compile loop nests into concurrent, dynamically constructed graphs of dependent tasks. We propose a simple and effective heuristic that selects the most profitable parallelization idiom for every dependence type and communication pattern. This heuristic enables the extraction of interband parallelism (cross-barrier parallelism) in a number of numerical computations that range from linear algebra to structured grids and image processing. The proposed static analysis and code generation alleviates the burden of a full-blown dependence resolver to track the readiness of tasks at runtime. We evaluate our approach and algorithms in the PPCG compiler, targeting OpenStream, a representative dataflow task-parallel language with explicit intertask dependences and a lightweight runtime. Experimental results demonstrate the effectiveness of the approach.

Faster computation of the Karhunen-Loeve expansion using its domain independence property

The goal of this work is to reduce the cost of computing the coefficients in the Karhunen-Loeve (KL) expansion. The KL expansion serves as a useful and efficient tool for discretizing second-order stochastic processes with known covariance function. Its applications in engineering mechanics include discretizing random field models for elastic moduli, fluid properties, and structural response. The main computational cost of finding the coefficients of this expansion arises from numerically solving an integral eigenvalue problem with the covariance function as the integration kernel. Mathematically this is a homogeneous Fredholm equation of second type. One widely used method for solving this integral eigenvalue problem is to use finite element (FE) bases for discretizing the eigenfunctions, followed by a Galerkin projection. This method is computationally expensive. In the current work it is first shown that the shape of the physical domain in a random field does not affect the realizations of the field estimated using KL expansion, although the individual KL terms are affected. Based on this domain independence property, a numerical integration based scheme accompanied by a modification of the domain, is proposed. In addition to presenting mathematical arguments to establish the domain independence, numerical studies are also conducted to demonstrate and test the proposed method. Numerically it is demonstrated that compared to the Galerkin method the computational speed gain in the proposed method is of three to four orders of magnitude for a two dimensional example, and of one to two orders of magnitude for a three dimensional example, while retaining the same level of accuracy. It is also shown that for separable covariance kernels a further cost reduction of three to four orders of magnitude can be achieved. Both normal and lognormal fields are considered in the numerical studies. (c) 2014 Elsevier B.V. All rights reserved.

NUMERICAL SIMULATIONS OF LIQUID JET BREAK UP IN A CROSSFLOW

Atomization is the process of disintegration of a liquid jet into ligaments and subsequently into smaller droplets. A liquid jet injected from a circular orifice into cross flow of air undergoes atomization primarily due to the interaction of the two phases rather than an intrinsic break up. Direct numerical simulation of this process resolving the finest droplets is computationally very expensive and impractical. In the present study, we resort to multiscale modelling to reduce the computational cost. The primary break up of the liquid jet is simulated using Gerris, an open source code, which employs Volume-of-Fluid (VOF) algorithm. The smallest droplets formed during primary atomization are modeled as Lagrangian particles. This one-way coupling approach is validated with the help of the simple test case of tracking a particle in a Taylor-Green vortex. The temporal evolution of the liquid jet forming the spray is captured and the flattening of the cylindrical liquid column prior to breakup is observed. The size distribution of the resultant droplets is presented at different distances downstream from the location of injection and their spatial evolution is analyzed.

Interphase anisotropy effects on lamellar eutectics: A numerical study

In directional solidification of binary eutectics, it is often observed that two-phase lamellar growth patterns grow tilted with respect to the direction z of the imposed temperature gradient. This crystallographic effect depends on the orientation of the two crystal phases alpha and beta with respect to z. Recently, an approximate theory was formulated that predicts the lamellar tilt angle as a function of the anisotropy of the free energy of the solid(alpha)-solid(beta) interphase boundary. We use two different numerical methods-phase field (PF) and dynamic boundary integral (BI)-to simulate the growth of steady periodic patterns in two dimensions as a function of the angle theta(R) between z and a reference crystallographic axis for a fixed relative orientation of alpha and beta crystals, that is, for a given anisotropy function (Wulff plot) of the interphase boundary. For Wulff plots without unstable interphase-boundary orientations, the two simulation methods are in excellent agreement with each other and confirm the general validity of the previously proposed theory. In addition, a crystallographic ``locking'' of the lamellae onto a facet plane is well reproduced in the simulations. When unstable orientations are present in the Wulff plot, it is expected that two distinct values of the tilt angle can appear for the same crystal orientation over a finite theta(R) range. This bistable behavior, which has been observed experimentally, is well reproduced by BI simulations but not by the PF model. Possible reasons for this discrepancy are discussed.

NUMERICAL STUDIES OF DYNAMO ACTION IN A TURBULENT SHEAR FLOW. I.

We perform numerical experiments to study the shear dynamo problem where we look for the growth of a large-scale magnetic field due to non-helical stirring at small scales in a background linear shear flow in previously unexplored parameter regimes. We demonstrate the large-scale dynamo action in the limit where the fluid Reynolds number (Re) is below unity while the magnetic Reynolds number (Rm) is above unity; the exponential growth rate scales linearly with shear, which is consistent with earlier numerical works. The limit of low Re is particularly interesting, as seeing the dynamo action in this limit would provide enough motivation for further theoretical investigations, which may focus attention on this analytically more tractable limit of Re < 1 compared to the more formidable limit of Re > 1. We also perform simulations in the regimes where (i) both (Re, Rm) < 1, and (ii) Re > 1 and Rm < 1, and compute all of the components of the turbulent transport coefficients (alpha(ij) and alpha(ij)) using the test-field method. A reasonably good agreement is observed between our results and the results of earlier analytical works in similar parameter regimes.

Numerical study of solidification of A356 aluminum alloy flowing on an oblique plate with experimental validation

Molten A356 aluminum alloy flowing on an oblique plate is water cooled from underneath. The melt partially solidifies on plate wall with continuous formation of columnar dendrites. These dendrites are continuously sheared off into equiaxed/fragmented grains and carried away with the melt by producing semisolid slurry collected at plate exit. Melt pouring temperature provides required solidification whereas plate inclination enables necessary shear for producing slurry of desired solid fraction. A numerical model concerning transport equations of mass, momentum, energy and species is developed for predicting velocity, temperature, macrosegregation and solid fraction. The model uses FVM with phase change algorithm, VOF and variable viscosity. The model introduces solid phase movement with gravity effect as well. Effects of melt pouring temperature and plate inclination on hydrodynamic and thermo-solutal behaviors are studied subsequently. Slurry solid fractions at plate exit are 27%, 22%, 16%, and 10% for pouring temperatures of 620 degrees C, 625 degrees C, 630 degrees C, and 635 degrees C, respectively. And, are 27%, 25%, 22%, and 18% for plate inclinations of 30, 45, 60, and 75, respectively. Melt pouring temperature of 625 degrees C with plate inclination of 60 generates appropriate quality of slurry and is the optimum. Both numerical and experimental results are in good agreement with each other. (C) 2015 Taiwan Institute of Chemical Engineers. Published by Elsevier B.V. All rights reserved.