An enthalpy-based model of dendritic growth in a convecting binary alloy melt

Purpose-In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues. Design/methodology/approach-The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed. Findings-Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated. Originality/value-As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.

Energy Loss in the Impact of Elastic Spheres on a Rigid Half-Space in Presence of Adhesion

Adhesion can cause energy losses in asperities or particles coming into dynamic contact resulting in frictional dissipation, even if the deformation occurring is purely elastic. Such losses are of special significance in impact of nanoparticles and friction between surfaces under low contact pressure to hardness ratio. The objective of this work is to study the effect of adhesion during the normal impact of elastic spheres on a rigid half-space, with an emphasis on understanding the mechanism of energy loss. We use finite element method for modeling the impact phenomenon, with the adhesion due to van der Waals force and the short-range repulsion included as body forces distributed over the volume of the sphere. This approach, in contrast with commonly used surface force approximation, helps to model the interactions in a more precise way. We find that the energy loss in impact of elastic spheres is negligible unless there are adhesion-induced instabilities. Significant energy loss through elastic stress waves occurs due to jump-to-contact and jump-out-of-contact instabilities and can even result in capture of the elastic sphere on the half-space.

A parametric study on the buckling of functionally graded material plates with internal discontinuities using the partition of unity method

In this paper, the effect of local defects, viz., cracks and cutouts on the buckling behaviour of functionally graded material plates subjected to mechanical and thermal load is numerically studied. The internal discontinuities, viz., cracks and cutouts are represented independent of the mesh within the framework of the extended finite element method and an enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The properties are assumed to vary only in the thickness direction and the effective properties are estimated using the Mori-Tanaka homogenization scheme. The plate kinematics is based on the first order shear deformation theory. The influence of various parameters, viz., the crack length and its location, the cutout radius and its position, the plate aspect ratio and the plate thickness on the critical buckling load is studied. The effect of various boundary conditions is also studied. The numerical results obtained reveal that the critical buckling load decreases with increase in the crack length, the cutout radius and the material gradient index. This is attributed to the degradation in the stiffness either due to the presence of local defects or due to the change in the material composition. (C) 2013 Elsevier Masson SAS. All rights reserved.

Automatic finite element formulation and assembly of hyperelastic higher order structural models

The formulation of higher order structural models and their discretization using the finite element method is difficult owing to their complexity, especially in the presence of non-linearities. In this work a new algorithm for automating the formulation and assembly of hyperelastic higher-order structural finite elements is developed. A hierarchic series of kinematic models is proposed for modeling structures with special geometries and the algorithm is formulated to automate the study of this class of higher order structural models. The algorithm developed in this work sidesteps the need for an explicit derivation of the governing equations for the individual kinematic modes. Using a novel procedure involving a nodal degree-of-freedom based automatic assembly algorithm, automatic differentiation and higher dimensional quadrature, the relevant finite element matrices are directly computed from the variational statement of elasticity and the higher order kinematic model. Another significant feature of the proposed algorithm is that natural boundary conditions are implicitly handled for arbitrary higher order kinematic models. The validity algorithm is illustrated with examples involving linear elasticity and hyperelasticity. (C) 2013 Elsevier Inc. All rights reserved.

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

Analysis of Entropy Generation During Mixed Convective Heat Transfer of Nanofluids Past a Rotating Circular Cylinder

The entropy generation due to mixed convective heat transfer of nanofluids past a rotating circular cylinder placed in a uniform cross stream is investigated via streamline upwind Petrov-Galerkin based finite element method. Nanosized copper (Cu) particles suspended in water are used with Prandtl number (Pr)=6.9. The computations are carried out at a representative Reynolds number (Re) of 100. The dimensionless cylinder rotation rate, a, is varied between 0 and 2. The range of nanoparticle volume fractions (phi) considered is 0 <= phi <= 5%. Effect of aiding buoyancy is brought about by considering two fixed values of the Richardson number (Ri) as 0.5 and 1.0. A new model for predicting the effective viscosity and thermal conductivity of dilute suspensions of nanoscale colloidal particles is presented. The model addresses the details of the agglomeration-deagglomeration in tune with the pertinent variations in the effective particulate dimensions, volume fractions, as well as the aggregate structure of the particulate system. The total entropy generation is found to decrease sharply with cylinder rotation rates and nanoparticle volume fractions. Increase in nanoparticle agglomeration shows decrease in heat transfer irreversibility. The Bejan number falls sharply with increase in alpha and phi.

A robust layer adapted difference method for singularly perturbed two-parameter parabolic problems

A finite difference method for a time-dependent singularly perturbed convection-diffusion-reaction problem involving two small parameters in one space dimension is considered. We use the classical implicit Euler method for time discretization and upwind scheme on the Shishkin-Bakhvalov mesh for spatial discretization. The method is analysed for convergence and is shown to be uniform with respect to both the perturbation parameters. The use of the Shishkin-Bakhvalov mesh gives first-order convergence unlike the Shishkin mesh where convergence is deteriorated due to the presence of a logarithmic factor. Numerical results are presented to validate the theoretical estimates obtained.

A Smooth Discretization Bridging Finite Element and Mesh-free Methods Using Polynomial Reproducing Simplex Splines

This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.

Performance Analysis of Heat Sinks With Phase-Change Materials Subjected to Transient and Cyclic Heating

Phase-change cooling technique is a suitable method for thermal management of electronic equipment subjected to transient or cyclic heat loads. The thermal performance of a phase-change based heat sink under cyclic heat load depends on several design parameters, namely, applied heat flux, cooling heat transfer coefficient, thermophysical properties of phase-change materials (PCMs), and physical dimensions of phase-change storage system during melting and freezing processes. A one-dimensional conduction heat transfer model is formulated to evaluate the effectiveness of preliminary design of practical PCM-based energy storage units. In this model, the phase-change process of the PCM is divided into melting and solidification subprocesses, for which separate equations are written. The equations are solved sequentially and an explicit closed-form solution is obtained. The efficacy of analytical model is estimated by comparing with a finite-volume-based numerical solution for both transient and cyclic heat loads.

Implicit scheme for meshless compressible Euler solver

In this paper, an implicit scheme is presented for a meshless compressible Euler solver based on the Least Square Kinetic Upwind Method (LSKUM). The Jameson and Yoon's split flux Jacobians formulation is very popular in finite volume methodology, which leads to a scalar diagonal dominant matrix for an efficient implicit procedure (Jameson & Yoon, 1987). However, this approach leads to a block diagonal matrix when applied to the LSKUM meshless method. The above split flux Jacobian formulation, along with a matrix-free approach, has been adopted to obtain a diagonally dominant, robust and cheap implicit time integration scheme. The efficacy of the scheme is demonstrated by computing 2D flow past a NACA 0012 airfoil under subsonic, transonic and supersonic flow conditions. The results obtained are compared with available experiments and other reliable computational fluid dynamics (CFD) results. The present implicit formulation shows good convergence acceleration over the RK4 explicit procedure. Further, the accuracy and robustness of the scheme in 3D is demonstrated by computing the flow past an ONERA M6 wing and a clipped delta wing with aileron deflection. The computed results show good agreement with wind tunnel experiments and other CFD computations.

A hybrid method for stochastic response analysis of a vibrating structure

Response analysis of a linear structure with uncertainties in both structural parameters and external excitation is considered here. When such an analysis is carried out using the spectral stochastic finite element method (SSFEM), often the computational cost tends to be prohibitive due to the rapid growth of the number of spectral bases with the number of random variables and the order of expansion. For instance, if the excitation contains a random frequency, or if it is a general random process, then a good approximation of these excitations using polynomial chaos expansion (PCE) involves a large number of terms, which leads to very high cost. To address this issue of high computational cost, a hybrid method is proposed in this work. In this method, first the random eigenvalue problem is solved using the weak formulation of SSFEM, which involves solving a system of deterministic nonlinear algebraic equations to estimate the PCE coefficients of the random eigenvalues and eigenvectors. Then the response is estimated using a Monte Carlo (MC) simulation, where the modal bases are sampled from the PCE of the random eigenvectors estimated in the previous step, followed by a numerical time integration. It is observed through numerical studies that this proposed method successfully reduces the computational burden compared with either a pure SSFEM of a pure MC simulation and more accurate than a perturbation method. The computational gain improves as the problem size in terms of degrees of freedom grows. It also improves as the timespan of interest reduces.

A rapid method to assess reactive oxygen species in yeast using H2DCF-DA

A protocol to efficiently assess Reactive Oxygen Species (ROS) levels in yeast cells using H2DCF-DA is described here. This method employs lithium acetate to permeate the cell wall, and thus, augments the release of the fluorescent product, dichlorofluorescein from the cells. This protocol obviates the need for both physical and enzymatic lysis methods that are arduous and time consuming. This method is simple, less time consuming and reproducible, especially while dealing with a large sample size. The lithium acetate method gave significantly reproducible and linear results (P < 0.0001), as compared with direct measurement (P = 0.0005), sonication (P = 0.1466) and bead beating (P = 0.0028).

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.