Paul R. Dawson¹, Donald E. Boyce², Ronald Rogge³

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.2, pp. 123-142, 2005, DOI:10.3970/cmes.2005.010.123

Abstract Computational mechanics provides a powerful environment for modeling the evolution of material structure during deformation processes and for associating that evolution with changes to the mechanical properties. In this paper, we illustrate a two-scale formulation that links the mechanical loading applied at the scale of a component (the continuum scale) to the responses of the material at the scale of the crystals that comprise it (the crystal scale). Employing the capabilities offered by computational mechanics, we can better understand how heterogeneity of deformation arising at both the continuum and crystal scales influences the behaviors observed More >

Zhenyu Xue¹, Ashkan Vaziri¹, John W. Hutchinson¹

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 79-96, 2005, DOI:10.3970/cmes.2005.010.079

Abstract A constitutive model for the elastic-plastic behavior of plastically compressible orthotropic materials is proposed based on an ellipsoidal yield surface with evolving ellipticity to accommodate non-uniform hardening or softening associated with stressing in different directions. The model incorporates rate-dependence arising from material rate-dependence and micro-inertial effects. The basic inputs are the stress-strain responses under the six fundamental stress histories in the orthotropic axes. Special limits of the model include classical isotropic hardening theory, the Hill model for incompressible orthotropic solids, and the Deshpande-Fleck model for highly porous isotropic foam metals. A primary motivation is application… More >

P. A. C. Porto¹, A. B. Jorge¹, G. O. Ribeiro²

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 65-78, 2005, DOI:10.3970/cmes.2005.010.065

Abstract This work deals with a numerical solution technique for the self-regular gradient form of Green's identity, the flux boundary integral equation (flux-BIE). The required C^1,α inter-element continuity conditions for the potential derivatives are imposed in the boundary element method (BEM) code through a non-symmetric variational formulation. In spite of using Lagrangian C⁰ elements, accurate numerical results were obtained when applied to heat transfer problems with singular or quasi-singular conditions, like boundary points and interior points which may be arbitrarily close to the boundary. The numerical examples proposed show that the developed algorithm based on the self-regular More >

C.-W. Chang¹, C.-S. Liu², J.-R. Chang^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 13-38, 2005, DOI:10.3970/cmes.2005.010.013

Abstract In this paper, the inverse heat conduction problem governed by sideways heat equation is investigated numerically. The problem is ill-posed because the solution, if it exists, does not depend continuously on the data. To begin with, this ill-posed problem is analyzed by considering the stability of the semi-discretization numerical schemes. Then the resulting ordinary differential equations at the discretized times are numerically integrated towards the spatial direction by the group preserving scheme, and the stable range of the index r = 1/2ν Δt is investigated. When the numerical results are compared with exact solutions, it More >

Z. D. Han¹, A. M. Rajendran², S.N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 1-12, 2005, DOI:10.3970/cmes.2005.010.001

Abstract A nonlinear formulation of the Meshless Local Petrov-Galerkin (MLPG) finite-volume mixed method is developed for the large deformation analysis of static and dynamic problems. In the present MLPG large deformation formulation, the velocity gradients are interpolated independently, to avoid the time consuming differentiations of the shape functions at all integration points. The nodal values of velocity gradients are expressed in terms of the independently interpolated nodal values of displacements (or velocities), by enforcing the compatibility conditions directly at the nodal points. For validating the present large deformation MLPG formulation, two example problems are considered: 1)… More >

Marcello Lappa¹

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.3, pp. 213-234, 2005, DOI:10.3970/fdmp.2005.001.213

Abstract The physical properties of many emulsions and metal alloys strongly depend on the multiphase morphology which is controlled to a great degree by particle-particle interaction during the related processing. In the present article significant effort is devoted to illustrate the philosophy of modeling for these phenomena and some insights into the physics. Within such a context working numerical techniques that have enjoyed a widespread use over recent years are presented and/or reviewed. Finally a focused and critical comparison of these possible approaches is reported illustrating advantages and disadvantages, strengths and weaknesses, past history and future More >

Marcello Lappa¹

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.2, pp. 171-188, 2005, DOI:10.3970/fdmp.2005.001.171

Abstract The paper presents a comparative and critical analysis of some theoretical/experimental/numerical arguments concerning the possible stabilization of the surface-tension-driven (Marangoni) flow in the Floating Zone technique and in various related fluid-dynamic models. It is conceived as a natural extension of the focused overview published in Cryst. Res. Tech. 40(6), 531, (2005) where much room was devoted to discuss the intrinsic physical mechanisms responsible for three-dimensional and oscillatory flows in a variety of technological processes. Here, a significant effort is provided to illustrate the genesis of possible control strategies (many of which are still in a More >

Cosmina S. Hogea¹, Bruce T. Murray¹, James A. Sethian^2,3

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.2, pp. 109-130, 2005, DOI:10.3970/fdmp.2005.001.109

Abstract A computational framework for simulating growth and transport in biological materials based on continuum models is proposed. The advantages of the finite difference methodology employed are generality and relative simplicity of implementation. The Cartesian mesh/level set method developed here provides a computational tool for the investigation of a host of transport-based tissue/tumor growth models, that are posed as free or moving boundary problems and may exhibit complicated boundary evolution including topological changes. The methodology is tested here on a widely studied "incompressible flow" type tumor growth model with a numerical implementation in two dimensions; comparisons More >

Gustav Amberg¹, Junichiro Shiomi²

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.1, pp. 81-96, 2005, DOI:10.3970/fdmp.2005.001.081

Abstract A few issues in materials science are reviewed with regard to the importance of fluid flows. The effect of convection on generic solidification problems is discussed. One relevant class of flows in melts is those driven by surface tension gradients. In welding this thermo- or solutocapillary flow will determine the penetration depth, and will depend very sensitively on the composition of the material, through the dependence of surface tension on temperature, presence of surfactants, etc. In crystal growth the convective motion in the melt may cause instabilities that are often undesired in practical processes. The More >

A. La Rocca¹, H. Power¹, V. La Rocca², M. Morale²

CMC-Computers, Materials & Continua, Vol.2, No.4, pp. 239-250, 2005, DOI:10.3970/cmc.2005.002.239

Abstract This work presents a meshless numerical scheme for the solution of time dependent non linear heat transfer problems in terms of a radial basis function Hermite collocation approach. The proposed scheme is applied to foodstuff's samples during freezing process; evaluation of the time evolution of the temperature profile along the sample, as well as at the core, is carried out. The moving phase-change zone is identified in the domain and plotted at several timesteps. The robustness of the proposed scheme is tested by a comparison of the obtained numerical results with those found using a More >