An operator-splitting finite element method for the efficient parallel solution of multidimensional population balance systems

A finite element method for solving multidimensional population balance systems is proposed where the balance of fluid velocity, temperature and solute partial density is considered as a two-dimensional system and the balance of particle size distribution as a three-dimensional one. The method is based on a dimensional splitting into physical space and internal property variables. In addition, the operator splitting allows to decouple the equations for temperature, solute partial density and particle size distribution. Further, a nodal point based parallel finite element algorithm for multi-dimensional population balance systems is presented. The method is applied to study a crystallization process assuming, for simplicity, a size independent growth rate and neglecting agglomeration and breakage of particles. Simulations for different wall temperatures are performed to show the effect of cooling on the crystal growth. Although the method is described in detail only for the case of d=2 space and s=1 internal property variables it has the potential to be extendable to d+s variables, d=2, 3 and s >= 1. (C) 2011 Elsevier Ltd. All rights reserved.

Growth and consumption rates of the phase layers during interdiffusion in a diffusion couple with finite end member

The relations for the growth and consumption rates of a layer with finite thickness as an end member and the product phases in the interdiffusion zone are developed. We have used two different methodologies, the diffusion based and the physico-chemical approach to develop the same relations. We have shown that the diffusion based approach is rather straightforward; however, the physico-chemical approach is much more versatile than the other method. It was found that the position of the marker plane becomes vague in the second stage of the interdiffusion process in pure A thin layer/B couple, where two phases grow simultaneously.

Rotation field in wedge indentation of metals

A study is made of the rotation field in wedge indentation of metals using copper as the model material system. Wedges with apical angles of 60 and 120 are used to indent annealed copper, and the deformation is mapped using image correlation. The indentation of annealed and strain-hardened copper is simulated using finite element analysis. The rotation field, derived from the deformation measurements, provides a clear way of distinguishing between cutting and compressive modes of deformation. Largely unidirectional rotation on one side of the symmetry line with small spatial rotation gradients is characteristic of compression. Bidirectional rotation with neighboring regions of opposing rotations and locally high rotation gradients characterizes cutting. In addition, the rotation demarcates such characteristic regions as the pile-up zone in indentation of a strain-hardened metal. The residual rotation field obtained after unloading is essentially the same as that at full load, indicating that it is a scalar proxy for plastic deformation as a whole.

Fracture analysis of cracked plate panels using NI-MVCCI technique

The objective of this paper is to propose a numerically integrated modified virtual crack closure integral (NI-MVCCI) technique for fracture analysis of cracked plate panels. NI-MVCCI technique is generalized one and the expressions for computing the strain energy release rate (SERR) are independent of the finite element employed. NI-MVCCI technique has been demonstrated for 4-noded, 8-noded (regular and quarter-point) and 9-noded isoparametric finite elements. Numerical studies on fracture analysis of 2-D crack (mode-I and mode-II) problems have been conducted employing these elements. SERR and stress intensity factors (SIF) have been computed for these problems and found to be in good agreement with the respective analytical solutions available in the literature. The appropriate Gauss numerical integration order to be employed for each of these elements for accurate computation of SERR and SIF has been recommended based on the studies.

An embedded grid adaptation strategy for unstructured data based finite volume computations

A grid adaptation strategy for unstructured data based codes, employing a combination of hexahedral and prismatic elements, generalizable to tetrahedral and pyramidal elements has been developed.

Appraisal of Bishop's method of slope stability analysis

The stability of slopes is a major problem in geotechnical engineering. Of the methods available for the analysis of soil slopes such as limit equilibrium methods, limit analysis and numerical methods such as FEM and FDM, limit equilibrium methods are popular and generally used, owing to their simplicity in formulation and in evaluating the overall factor of safety of slope. However limit equilibrium methods possess certain disadvantages. They do not consider whether the slope is an embankment or natural slope or an excavation and ignore the effect of incremental construction, initial stress, stress strain behavior etc. In the work reported in this paper, a comparative study of actual state of stress and actual factor of safety and Bishop's factor of safety is performed. The actual factor of safety is obtained by consideration of contours of mobilised shear strains. Using Bishop's method of slices, the critical slip surfaces of a number of soil slopes with different geometries are determined and both the factors of safety are obtained. The actual normal stresses and shear stresses are determined from finite difference formulation using FLAG (Fast Lagrangian Analysis of Continuaa) with Mohr-Coulomb model. The comparative study is performed in terms of parameter lambda(c phi) (= gamma H tan phi/c). I is shown that actual factor of safety is higher than Bishop's factor of safety depending on slope angle and lambda(c phi).

Analysis of the consolidation of preloaded clays

Precompression, wherein the probable settlements are achieved at an accelerated pace through preloading, well before the construction is take up, has been widely used in areas of ground improvement with respect to soft clays. By applying a temporary surcharge load in excess or less than the permanent load, the soil achieves higher initial effective stress and when the final load is applied, the soil experiences, lower load increment ratio or negative load increment ratio. In this paper, based on the laboratory experiments conducted on cochin marine clays and Mangalore marine clays, attempts have been made to analyse the volume change behaviour of preloaded clays. It has been brought out that for a preloaded clay, the final load increment ratio has an important role in its behaviour. Effective preloading not only reduces the final settlement due to primary, the secondary consolidation settlement also gets reduced.

3D finite element modelling of a composite patch repair

Composite-patching on cracked/weak metallic aircraft structures improves structural integrity. A Boron Epoxy patch employed to repair a cracked Aluminum sheet is modeled employing 3D Finite Element Method (FEM). SIFs extracted using ''displacement extrapolation'' are used to measure the repair effectiveness. Two issues viz., patch taper and symmetry have been looked into.

Exciting forces due to interaction of water waves with a submerged sphere in an ice-covered two-layer fluid of finite depth

Scattering of water waves by a sphere in a two-layer fluid, where the upper layer has an ice-cover modelled as an elastic plate of very small thickness, while the lower one has a rigid horizontal bottom surface, is investigated within the framework of linearized water wave theory. The effects of surface tension at the surface of separation is neglected. There exist two modes of time-harmonic waves - the one with lower wave number propagating along the ice-cover and the one with higher wave number along the interface. Method of multipole expansions is used to find the particular solution for the problem of wave scattering by a submerged sphere placed in either of the layers. The exciting forces for vertical and horizontal directions are derived and plotted against different values of the wave number for different submersion depths of the sphere and flexural rigidity of the ice-cover. When the flexural rigidity and the density of the ice-cover are taken to be zero, the numerical results for the exciting forces for the problem with free surface are recovered as particular cases. (C) 2011 Elsevier Ltd. All rights reserved.

Unitary invariants for Hilbert modules of finite rank

We associate a sheaf model to a class of Hilbert modules satisfying a natural finiteness condition. It is obtained as the dual to a linear system of Hermitian vector spaces (in the sense of Grothendieck). A refined notion of curvature is derived from this construction leading to a new unitary invariant for the Hilbert module. A division problem with bounds, originating in Douady's privilege, is related to this framework. A series of concrete computations illustrate the abstract concepts of the paper.

Effect of strain rate and temperature on the plastic deformation behaviour of a bulk metallic glass composite

The composites consisting of amorphous matrix reinforced with crystalline dendrites offer extraordinary combinations of strength, stiffness, and toughness and can be processed in bulk. Hence, they have been receiving intense research interest, with a primary focus to study their mechanical properties. In this paper, the temperature and strain rate effects on the uniaxial compression response of a tailored bulk metallic glass (BMG) composite has been investigated. Experimental results show that at temperatures ranging between ambient to 500 K and at all strain rates; the onset of plastic deformation in the composite is controlled by that in the dendrites. As the temperature is increased to the glass transition temperature of the matrix and beyond, flow in the amorphous matrix occurs readily and hence it dictates the composite's response. The role of the constituent phases in controlling the deformation mechanism of the composite has been verified by assessing the strain rate sensitivity and the activation volume for deformation. The composite is rate sensitive at room temperature with values of strain rate sensitivity and activation volume being similar to that of the dendrites. At test temperatures near to the glass transition temperature, the composite however becomes rate-insensitive corresponding to that of the matrix phase. At low strain rates, serrated flow akin to that of dynamic strain ageing in crystalline alloys was observed and the serration magnitude decreases with increasing temperature. Initiation of the shear bands at the dendrite/matrix interface and propagation of them through the matrix ligaments until their arrest at another interface is the responsible mechanism for this. (C) 2011 Elsevier B.V. All rights reserved.

Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants

A finite-element scheme based on a coupled arbitrary Lagrangian-Eulerian and Lagrangian approach is developed for the computation of interface flows with soluble surfactants. The numerical scheme is designed to solve the time-dependent Navier-Stokes equations and an evolution equation for the surfactant concentration in the bulk phase, and simultaneously, an evolution equation for the surfactant concentration on the interface. Second-order isoparametric finite elements on moving meshes and second-order isoparametric surface finite elements are used to solve these equations. The interface-resolved moving meshes allow the accurate incorporation of surface forces, Marangoni forces and jumps in the material parameters. The lower-dimensional finite-element meshes for solving the surface evolution equation are part of the interface-resolved moving meshes. The numerical scheme is validated for problems with known analytical solutions. A number of computations to study the influence of the surfactants in 3D-axisymmetric rising bubbles have been performed. The proposed scheme shows excellent conservation of fluid mass and of the total mass of the surfactant. (C) 2012 Elsevier Inc. All rights reserved.

Shear sum rules at finite chemical potential

We derive sum rules which constrain the spectral density corresponding to the retarded propagator of the T-xy component of the stress tensor for three gravitational duals. The shear sum rule is obtained for the gravitational dual of the N = 4 Yang-Mills, theory of the M2-branes and M5-branes all at finite chemical potential. We show that at finite chemical potential there are additional terms in the sum rule which involve the chemical potential. These modifications are shown to be due to the presence of scalars in the operator product expansion of the stress tensor which have non-trivial vacuum expectation values at finite chemical potential.

Quasi-static crack tip fields in rate-sensitive FCC single crystals

In this work, the effects of loading rate, material rate sensitivity and constraint level on quasi-static crack tip fields in a FCC single crystal are studied. Finite element simulations are performed within a mode I, plane strain modified boundary layer framework by prescribing the two term (K-T) elastic crack tip field as remote boundary conditions. The material is assumed to obey a rate-dependent crystal plasticity theory. The orientation of the single crystal is chosen so that the crack surface coincides with the crystallographic (010) plane and the crack front lies along 101] direction. Solutions corresponding to different stress intensity rates K., T-stress values and strain rate exponents m are obtained. The results show that the stress levels ahead of the crack tip increase with K. which is accompanied by gradual shrinking of the plastic zone size. However, the nature of the shear band patterns around the crack tip is not affected by the loading rate. Further, it is found that while positive T-stress enhances the opening and hydrostatic stress levels ahead of crack tip, they are considerably reduced with imposition of negative T-stress. Also, negative T-stress promotes formation of shear bands in the forward sector ahead of the crack tip and suppresses them behind the tip.