Hokyung Shim¹, Vinh Thuy Tran Ho^1,,Semyung Wang^1,2, Daniel A. Tortorelli³

CMES-Computer Modeling in Engineering & Sciences, Vol.37, No.2, pp. 175-202, 2008, DOI:10.3970/cmes.2008.037.175

Abstract This paper presents a topological shape optimization technique for electromagnetic problems using a level set method and radial basis functions. The proposed technique is a level set (LS) based optimization dealing with geometrical shape derivatives and topological design. The shape derivative is computed by an adjoint variable method to avoid numerous sensitivity evaluations. A level set model embedded into the scalar function of higher dimensions is propagated to represent the design boundary of a domain. The level set function interpolated into a fixed initial domain is evolved by using the Hamilton-Jacobi equation. The moving free… More >

A.I. Abreu¹, W.J. Mansur¹, D. Soares Jr^1,2, J.A.M. Carrer³

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 123-132, 2008, DOI:10.3970/cmes.2008.030.123

Abstract This paper is concerned with the numerical computation of space derivatives of a time-domain (TD-) Boundary Element Method (BEM) formulation for the analysis of scalar wave propagation problems. In the present formulation, the Convolution Quadrature Method (CQM) is adopted, i.e., the basic integral equation of the TD-BEM is numerically substituted by a quadrature formula, whose weights are computed using the Laplace transform of the fundamental solution and a linear multi-step method. In order to numerically compute space derivatives, the present work properly transforms the quadrature weights of the CQM-BEM, adopting the so-called Complex-Variable-Differentiation Method (CVDM). More >

K. Hashiguchi¹

The International Conference on Computational & Experimental Engineering and Sciences, Vol.4, No.4, pp. 257-264, 2007, DOI:10.3970/icces.2007.004.257

Abstract General corotational rate of tensors in arbitrary order having the objectivity is shown first, and then it is verified that the material-derivative of yield condition can be replaced generally to the corotational derivative, i.e. the consistency condition. More >

K. Hashiguchi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.1, pp. 55-62, 2007, DOI:10.3970/cmes.2007.017.055

Abstract Constitutive equation describing the mechanical properties of material has to be formulated in an identical form independent of coordinate systems by which it is described even if there exist any mutual configuration and/or mutual rotation between the material and coordinate systems. This mechanical requirement is attained by describing rate variables as corotational rate tensors with objectivity in constitutive equations in rate form. Besides, in order to use the material-time derivative of yield condition as a consistency condition it has to be replaced to the corotational derivative. In this note a general corotational rate for tensors More >

A. Frangi¹, L. Ghezzi, P. Faure-Ragani²

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.1, pp. 41-48, 2006, DOI:10.3970/cmes.2006.015.041

Abstract Three dimensional magneto-mechanical problems at low frequency are addressed by means of a coupled fast Boundary Element - Finite Element approach with total scalar potential and focusing especially on the issue of global force calculation on movable ferromagnetic parts. The differentiation of co-energy in this framework and the use of Maxwell tensor are critically discussed and the intrinsic links are put in evidence. Three examples of academic and industrial applications are employed for validation. More >

A.I.Tolstykh, D.A. Shirobokov¹

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.2, pp. 207-222, 2005, DOI:10.3970/cmes.2005.007.207

Abstract A way of using RBF as the basis for PDE's solvers is presented, its essence being constructing approximate formulas for derivatives discretizations based on RBF interpolants with local supports similar to stencils in finite difference methods. Numerical results for different types of elasticity equations showing reasonable accuracy and good$h$-convergence properties of the technique are presented. Applications of the technique to problems with non-self-adjoint operators (like those for the Navier-Stokes equations) are also considered. More >

Zai You Yan^1,2, Fang Sen Cui², Kin Chew Hung²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.1, pp. 97-106, 2005, DOI:10.3970/cmes.2005.007.097

Abstract Taking the normal derivative of solid angles on the surface into account, a modified Burton and Miller's formulation is derived. From which, a more reasonable expression of the hypersingular operator is obtained. To overcome the hypersingular integral, the regularization scheme developed recently is employed. Plane acoustic wave scattering from a rigid sphere is computed to show the correctness of the modified formulation with the regularization scheme. In the computation, eight-nodded isoparametric element is applied. More >

S. Choi¹, M.D. Marcozzi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.3, pp. 201-212, 2004, DOI:10.3970/cmes.2004.005.201

Abstract The general intractability of derivative security pricing models to numerical techniques arguably remains one of the preeminant problems of mathematical finance. In particular, the valuations of such models may be represented as solutions of variational inequalities of evolutionary type typically characterized by their high number of degrees of freedom, unbounded domains, and asymptotic behavior. We consider the application of Multi-Quadratic Radial Basis Functions (MQ-RBF) to the problem of option pricing. More >

W.M.G.. So¹, K.J. Lau¹, S.W. Ng¹

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.3, pp. 189-200, 2004, DOI:10.3970/cmes.2004.005.189

Abstract A new finite element analysis procedure is implemented for the determination of complex stress intensity factors in interfacial cracks. Only nodal displacements and strain energies of the near-crack-tip elements are involved in this procedure so that element stiffness matrices need not be made available. The method is first tested using a closed form solution for infinite media to obtain a suitable finite element mesh. It is then applied to finite plates and four-point bending specimens containing interfacial cracks. In cases where reference values are available for comparison, good agreement of results can be obtained with More >