A.R.Skovoroda¹, R.V.Goldstein²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 457-472, 2003, DOI:10.3970/cmes.2003.004.457

Abstract We discuss the inverse problem of the recovery of the distribution of the elastic moduli of a stratified layer, based on measurements of the surface displacement under localized surface loads. A general parametric solution and a numerical procedure for computing the parameters are presented. Examples of numerical results are given. The problem and its solution are related to the monitoring of elastic properties of living tissues. More >

I. J. Rao¹, J.D. Humphrey², K.R. Rajagopal³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 439-456, 2003, DOI:10.3970/cmes.2003.004.439

Abstract Recent discoveries in molecular and cell biology reveal that many cell types sense and respond (via altered gene expression) to changes in their mechanical environment. Such mechanotransduction mechanisms are responsible for many changes in structure and function, including the growth and remodeling process. To understand better, and ultimately to use (e.g., in tissue engineering), biological growth and remodeling, there is a need for mathematical models that have predictive and not just descriptive capability. In contrast to prior models based on reaction-diffusion equations or the concept of volumetric growth, we examine here a newly proposed constrained More >

K.J. Fischer^1,2, J.A. Bastidas³, H.J. Pfaeffle², J.D. Towers²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.3&4, pp. 397-404, 2003, DOI:10.3970/cmes.2003.004.397

Abstract The two bones of the forearm, the radius and the ulna, have been shown to bear different proportions of the overall forearm load at the wrist and the elbow. This biomechanical data suggests load transfer between the bones occurs through the soft tissues of the forearm. Load transfer from radius to ulna through passive soft tissues such as the interosseous ligament (IOL) has been experimentally measured. Ex vivo studies of the forearm, however, cannot account for the effect of internal loads generated by the muscles and, in some cases, external forces acting directly on the… More >

Jan Sladek¹, Vladimir Sladek¹, Chuanzeng Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 637-648, 2003, DOI:10.3970/cmes.2003.004.637

Abstract A new computational method for solving transient elastodynamic initial-boundary value problems in continuously non-homogeneous solids, based on the meshless local Petrov-Galerkin (MLPG) method, is proposed in the present paper. The moving least squares (MLS) is used for interpolation and the modified fundamental solution as the test function. The local Petrov-Galerkin method for unsymmetric weak form in such a way is transformed to the local boundary integral equations (LBIE). The analyzed domain is divided into small subdomains, in which a weak solution is assumed to exist. Nodal points are randomly spread in the analyzed domain and More >

Q. Li¹, S. Shen¹, Z. D. Han¹, S. N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 571-586, 2003, DOI:10.3970/cmes.2003.004.571

Abstract In this paper, a truly meshless method, the Meshless Local Petrov-Galerkin (MLPG) Method, is developed for three-dimensional elasto-statics. The two simplest members of MLPG family of methods, the MLPG type 5 and MLPG type 2, are combined, in order to reduce the computational requirements and to obtain high efficiency. The MLPG5 method is applied at the nodes inside the 3-D domain, so that any domain integration is eliminated altogether, if no body forces are involved. The MLPG 2 method is applied at the nodes that are on the boundaries, and on the interfaces of material More >

A. Aimi¹, M. Diligenti¹, F. Lunardini¹, A. Salvadori²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.1, pp. 31-50, 2003, DOI:10.3970/cmes.2003.004.031

Abstract This paper is devoted to the study of a new application of the Panel Clustering Method [Hackbusch and Sauter (1993); Hackbusch and Nowak (1989)]. By considering a classical 3D Neumann screen problem in its boundary integral formulation discretized with the Galerkin BEM, which requires the evaluation of double integrals with hypersingular kernel, we recall and use some recent results of analytical evaluation of the inner hypersingular integrals. Then we apply the Panel Clustering Method (PCM) for the evaluation of the outer integral. For this approach error estimate is shown. Numerical examples and comparisons with classical More >

Gerhard Klimeck^1,2, Fabiano Oyafuso², Timothy B. Boykin³, R. Chris Bowen², Paul von Allmen⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.5, pp. 601-642, 2002, DOI:10.3970/cmes.2002.003.601

Abstract Material layers with a thickness of a few nanometers are common-place in today's semiconductor devices. Before long, device fabrication methods will reach a point at which the other two device dimensions are scaled down to few tens of nanometers. The total atom count in such deca-nano devices is reduced to a few million. Only a small finite number of "free'' electrons will operate such nano-scale devices due to quantized electron energies and electron charge. This work demonstrates that the simulation of electronic structure and electron transport on these length scales must not only be fundamentally… More >

Amitesh Maiti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.5, pp. 589-600, 2002, DOI:10.3970/cmes.2002.003.589

Abstract Atomistic modeling and simulations are becoming increasingly important in the design of new devices at the nanoscale. In particular, theoretical modeling of carbon nanotubes have provided useful insight and guidance to many experimental efforts. To this end, we report simulation results on the electronic, structural and transport properties for two different applications of carbon nanotube-based devices: (1) effect of adsorbates on field emission; and (2) effect of mechanical deformation on the electronic transport. The reported simulations are based on First Principles Density Functional Theory (DFT), classical molecular mechanics, and tight-binding transport based on the recursive More >

Leonid V. Zhigilei¹, Avinash M. Dongare¹

CMES-Computer Modeling in Engineering & Sciences, Vol.3, No.5, pp. 539-556, 2002, DOI:10.3970/cmes.2002.003.539

Abstract Computational modeling has a potential of making an important contribution to the advancement of laser-driven methods in nanotechnology. In this paper we discuss two computational schemes developed for simulation of laser coupling to organic materials and metals and present a multiscale model for laser ablation and cluster deposition of nanostructured materials. In the multiscale model the initial stage of laser ablation is reproduced by the classical molecular dynamics (MD) method. For organic materials, the breathing sphere model is used to simulate the primary laser excitations and the vibrational relaxation of excited molecules. For metals, the… More >

S.G. Bardenhagen¹, J.E. Guilkey², K.M. Roessig³, J.U. Brackbill⁴, W.M. Witzel⁵, J.C.Foster⁶

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.4, pp. 509-522, 2001, DOI:10.3970/cmes.2001.002.509

Abstract Contact between deformable bodies is a difficult problem in the analysis of engineering systems. A new approach to contact has been implemented using the Material Point Method for solid mechanics, Bardenhagen, Brackbill, and Sulsky (2000a). Here two improvements to the algorithm are described. The first is to include the normal traction in the contact logic to more appropriately determine the free separation criterion. The second is to provide numerical stability by scaling the contact impulse when computational grid information is suspect, a condition which can be expected to occur occasionally as material bodies move through… More >