A fixed-grid finite element based enthalpy formulation for generalized phase change problems: role of superficial mushy region

The enthalpy method is primarily developed for studying phase change in a multicomponent material, characterized by a continuous liquid volume fraction (phi(1)) vs temperature (T) relationship. Using the Galerkin finite element method we obtain solutions to the enthalpy formulation for phase change in 1D slabs of pure material, by assuming a superficial phase change region (linear (phi(1) vs T) around the discontinuity at the melting point. Errors between the computed and analytical solutions are evaluated for the fluxes at, and positions of, the freezing front, for different widths of the superficial phase change region and spatial discretizations with linear and quadratic basis functions. For Stefan number (St) varying between 0.1 and 10 the method is relatively insensitive to spatial discretization and widths of the superficial phase change region. Greater sensitivity is observed at St = 0.01, where the variation in the enthalpy is large. In general the width of the superficial phase change region should span at least 2-3 Gauss quadrature points for the enthalpy to be computed accurately. The method is applied to study conventional melting of slabs of frozen brine and ice. Regardless of the forms for the phi(1) vs T relationships, the thawing times were found to scale as the square of the slab thickness. The ability of the method to efficiently capture multiple thawing fronts which may originate at any spatial location within the sample, is illustrated with the microwave thawing of slabs and 2D cylinders. (C) 2002 Elsevier Science Ltd. All rights reserved.

Automatic decomposition of unstructured meshes employing genetic algorithms for parallel FEM computations

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

Effect of contact at the interface on the compressive properties of fly ash-epoxy composites

Particulate composites based on polymer matrices generally contain fillers, especially those that are abundantly available and are cheaper. The inclusion of these, besides improving the properties, makes the system costwise viable, In the present study, fly ash was tried as a filler in epoxy. The filler particle surfaces were modified using three chemical surface treatment techniques in order to elicit the effect of adhesion at the interface on the mechanical properties of these composites. The compatibilizing of the filler with the use of a silane coupling agent yielded the best compression strength values. Scanning Electron Microscopy (SEM) has been used to characterize and supplement the mechanical test data.

Modeling of concrete cracking-A hybrid technique of using displacement discontinuity element method and direct boundary element method

This paper presents a new approach by making use of a hybrid method of using the displacement discontinuity element method and direct boundary element method to model concrete cracking by incorporating fictitious crack model. Fracture mechanics approach is followed using the Hillerborg's fictitious crack model. A boundary element based substructure method and a hybrid technique of using displacement discontinuity element method and direct boundary element method are compared in this paper. In order to represent the process zone ahead of the crack, closing forces are assumed to act in such a way that they obey a linear normal stress-crack opening displacement law. Plain concrete beams with and without initial crack under three-point loading were analyzed by both the methods. The numerical results obtained were shown to agree well with the results from existing finite element method. The model is capable of reproducing the whole range of load-deflection response including strain-softening and snap-back behavior as illustrated in the numerical examples. (C) 2011 Elsevier Ltd. All rights reserved.

Influence of the roughness of aggregate surface on the interface bond strength

An experimental investigation on the bond strength of the interface between mortar and aggregate is reported. Composite compact specimens were used for applying Mode I and Mode 11 loading effects. The influence of the type of mortar and type of aggregate and its roughness on the bond strength of the interface has been studied. It has been observed that the bond strength of the interface in tension is significantly low, though the mortars exhibited higher strength. The highest tensile bond strength values have been observed with rough concrete surface with M-13 mortar. The bond strength of the interface in Mode I load depends on the type of aggregate surface and its roughness, and the type of mortar, The bond strength of the interface between mortar M-13 cast against rough concrete in direct tension seems to be about one third of the strength of the mortar. However, it is about 1/20th to 1/10th with the mortar M-12 in sandwiched composite specimens. The bond strength of the interface in shear (Mode IT) significantly increases as the roughness and the phase angle of the aggregate surface increase. The strength of mortar on the interface bond strength has been very significant. The sandwiched composite specimens show relatively low bond strength in Mode I loading. The behavior of the interface in both Mode I and Mode 11 loading effects has been brittle, indicating catastrophic failure. (C) 2002 Elsevier Science Ltd. All rights reserved.

Modelling, simulation, and graphical user interface for industrial gas carburising process

A mathematical model has been developed for the gas carburising (diffusion) process using finite volume method. The computer simulation has been carried out for an industrial gas carburising process. The model's predictions are in good agreement with industrial experimental data and with data collected from the literature. A study of various mass transfer and diffusion coefficients has been carried out in order to suggest which correlations should be used for the gas carburising process. The model has been interfaced in a Windows environment using a graphical user interface. In this way, the model is extremely user friendly. The sensitivity analysis of various parameters such as initial carbon concentration in the specimen, carbon potential of the atmosphere, temperature of the process, etc. has been carried out using the model.

Systematic trends in structural and thermodynamic properties of ternary oxides in the systems Ln-Pd-O (Ln = lanthanide element)

Isothermal sections of the phase diagrams for the systems Ln-Pd-O (Ln = lanthanide element) at 1223 K indicate the presence of two inter-oxide compounds Ln(4)PdO(7) and Ln(2)Pd(2)O(5) for Ln = La, Pr, Nd, Sm, three compounds Ln(4)PdO(7), Ln(2)PdO(4) and Ln(2)Pd(2)O(5) for Ln = Eu, Gd and only one compound of Ln(2)Pd(2)O(5) for Ln = Tb to Ho. The lattice parameters of the compounds Ln(4)PdO(7), Ln(2)PdO(4) and Ln(2)Pd(2)O(5) show systematic nonlinear variation with atomic number. The unit cell volumes decrease with increasing atomic number. The standard Gibbs energies, enthalpies and entropies of formation of the ternary oxides from their component binary oxides (Ln(2)O(3) and PdO) have been measured recently using an advanced version of the solid-state electrochemical cell. The Gibbs energies and enthalpies of formation become less negative with increasing atomic number of Ln. For all the three compounds, the variation in Gibbs energy and enthalpy of formation with atomic number is markedly non-linear. The decrease in stability with atomic number is most pronounced for Ln(2)Pd(2)O(5), followed by Ln(4)PdO(7) and Ln(2)PdO(4). This is probably related to the repulsion between Pd2+ ions on the opposite phases Of O-8 cubes in Ln(2)Pd(2)O(5), and the presence of Ln-filled O-8 cubes that share three faces with each other in Ln4PdO7. The values for entropy of formation of all the ternary oxides from their component binary oxides are relatively small. Although the entropies of formation show some scatter, the average value for Ln = La, Pr, Nd is more negative than the average value for the other lanthanide elements. From this difference, an average value for the structure transformation entropy of Ln(2)O(3) from C-type to A-type is estimated as 0.87 J.mol(-1).K-1. The standard Gibbs energies of formation of these ternary oxides from elements at 1223 K are presented as a function of lanthanide atomic number. By invoking the Neumann-Kopp rule for heat capacity, thermodynamic properties of the inter-oxide compounds at 298.15 K are estimated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Nonlinear Finite-Element Modeling of Batter Piles under Lateral Load

In this paper, a finite-element model is developed in which the nonlinear soil behavior is represented by a hyperbolic relation for static load condition and modified hyperbolic relation, which includes both degradation and gap for a cyclic load condition. Although batter piles are subjected to lateral load, the soil resistance is also governed by axial load, which is incorporated by considering the P-Δ moment and geometric stiffness matrix. By adopting the developed numerical model, static and cyclic load analyses are performed adopting an incremental-iterative procedure where the pile is idealized as beam elements and the soil as elastoplastic spring elements. The proposed numerical model is validated with published laboratory and field pile test results under both static and cyclic load conditions. This paper highlights the importance of the degradation factor and its influence on the soil resistance-displacement (p-y) curve, number of cycles of loading, and cyclic load response.

Novel phase-transition behavior near liquid/liquid critical points of aqueous solutions: Formation of a third phase at the interface

We report a novel phase behavior in aqueous solutions of simple organic solutes near their liquid/liquid critical points, where a solid-like third phase appears at the liquid/liquid interface. The phenomenon has been found in three different laboratories. It appears in many aqueous systems of organic solutes and becomes enhanced upon the addition of salt to these solutions.

Water ingression and the role of interface on the compressive properties of epoxy filled with treated fly ash

Polymer composites are generally filled with either fibrous or particulate materials to improve the mechanical properties. In choosing the fillers one looks for materials that are inexpensive and available in abundance, in order to realize a cost reduction also. Also, often these fibres/fillers are treated to improve the matrix adhesion and thereby mechanical properties. The present study is focussed on the influence of water ingression in such filler-modified composites and the attendant changes in the compressive properties. The changes in property effected following exposure to aqueous media and the influence interface modification has on the scenario is emphasized in the work. It is seen that for plain epoxy and fly ash filled systems the strengths are increased following exposure to aqueous media. The composites with surface-treated ash particles, on the other hand, record a drop in the values. Modulus values show are increased to varying degree in unfilled and filled systems. The study also includes a fractographic analysis of the tested samples with and without exposure to water.

A new beam finite element for the analysis of functionally graded materials

A new beam element is developed to study the thermoelastic behavior of functionally graded beam structures. The element is based on the first-order shear deformation theory and it accounts for varying elastic and thermal properties along its thickness. The exact solution of static part of the governing differential equations is used to construct interpolating polynomials for the element formulation. Consequently, the stiffness matrix has super-convergent property and the element is free of shear locking. Both exponential and power-law variations of material property distribution are used to examine different stress variations. Static, free vibration and wave propagation problems are considered to highlight the behavioral difference of functionally graded material beam with pure metal or pure ceramic beams. (C) 2003 Elsevier Science Ltd. All rights reserved.

Reconstructing Solid Model from 2D Scanned Images of Biological Organs for Finite Element Simulation

This work presents a methodology to reconstruct 3D biological organs from image sequences or other scan data using readily available free softwares with the final goal of using the organs (3D solids) for finite element analysis. The methodology deals with issues such as segmentation, conversion to polygonal surface meshes, and finally conversion of these meshes to 3D solids. The user is able to control the detail or the level of complexity of the solid constructed. The methodology is illustrated using 3D reconstruction of a porcine liver as an example. Finally, the reconstructed liver is imported into the commercial software ANSYS, and together with a cyst inside the liver, a nonlinear analysis performed. The results confirm that the methodology can be used for obtaining 3D geometry of biological organs. The results also demonstrate that the geometry obtained by following this methodology can be used for the nonlinear finite element analysis of organs. The methodology (or the procedure) would be of use in surgery planning and surgery simulation since both of these extensively use finite elements for numerical simulations and it is better if these simulations are carried out on patient specific organ geometries. Instead of following the present methodology, it would cost a lot to buy a commercial software which can reconstruct 3D biological organs from scanned image sequences.

A nonlinear spectral finite element model for analysis of wave propagation in solid with internal friction and dissipation

A geometrically non-linear Spectral Finite Flement Model (SFEM) including hysteresis, internal friction and viscous dissipation in the material is developed and is used to study non-linear dissipative wave propagation in elementary rod under high amplitude pulse loading. The solution to non-linear dispersive dissipative equation constitutes one of the most difficult problems in contemporary mathematical physics. Although intensive research towards analytical developments are on, a general purpose cumputational discretization technique for complex applications, such as finite element, but with all the features of travelling wave (TW) solutions is not available. The present effort is aimed towards development of such computational framework. Fast Fourier Transform (FFT) is used for transformation between temporal and frequency domain. SFEM for the associated linear system is used as initial state for vector iteration. General purpose procedure involving matrix computation and frequency domain convolution operators are used and implemented in a finite element code. Convergnence of the spectral residual force vector ensures the solution accuracy. Important conclusions are drawn from the numerical simulations. Future course of developments are highlighted.