Nagalinga Rajan, Soumyendu Raha¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 91-116, 2008, DOI:10.3970/cmes.2008.023.091

Abstract The article describes a numerical method for time domain integration of noisy dynamical systems originating from engineering applications. The models are second order stochastic differential equations (SDE). The stochastic process forcing the dynamics is treated mainly as multiplicative noise involving a Wiener Process in the Itô sense. The developed numerical integration method is a drift implicit strong order 2.0 method. The method has user-selectable numerical dissipation properties that can be useful in dealing with both multiplicative noise and stiffness in a computationally efficient way. A generalized analysis of the method including the multiplicative noise is presented. Strong order convergence, user-selectable… More >

R.L. Huston¹, C.-Q. Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.2, pp. 143-152, 2005, DOI:10.3970/cmes.2005.010.143

Abstract This paper presents a summary of recent developments in computational methods for multibody dynamics analyses. The developments are presented within the context of an automated numerical analysis. The intent of the paper is to provide a basis for the easy development of computational algorithms. The principal concepts discussed are: differentiation algorithms, partial velocities and partial angular velocities, generalized speeds, Euler parameters, Kane's equations, orthogonal complement arrays, lower body arrays and accuracy testing functions. More >

F. Cosmi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.3, pp. 365-372, 2001, DOI:10.3970/cmes.2001.002.365

Abstract The aim of this paper is to present a methodology for solving the plane elasticity problem using the Cell Method. It is shown that with the use of a parabolic interpolation in a vectorial problem, a convergence rate of 3.5 is obtained. Such a convergence rate compares with, or is even better than, the one obtained with FEM with the same interpolation – depending on the integration technique used by the FEM application. The accuracy of the solution is also comparable or better. More >

P.M.J. Trevelyan¹, L. Elliott¹, D.B. Ingham¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 117-122, 2000, DOI:10.3970/cmes.2000.001.419

Abstract The potential flow in a semi-infinite channel with multiple semi-infinite oblique sub-channels is determined using the Schwarz-Christoffel transformation and complex potential theory. The standard iterative technique, i.e. the Newton-Raphson method with the Jacobian matrix approximated by a finite-difference quotient matrix, was employed with an alternative integration region to that found elsewhere in the literature is employed after integrating across the boundaries to determine the Schwarz-Christoffel transformation parameters which solely depend on the dimensions of the region being considered. Each semi-infinte channel permits integration at infinity perpendicularly across the channel and sub-channels, yielding some analytical relationships between these parameters. The remainder… More >

V. Cristini¹, Y. Renardy²

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.2, pp. 77-94, 2006, DOI:10.3970/fdmp.2006.002.077

Abstract We review studies of a drop of viscous liquid, suspended in another liquid, and undergoing breakup in an impulsively started shear flow. Stokes flow conditions as well as the effects of inertia are reported. They reveal a universal scaling for the fragments, which allows one to use sheared emulsions to produce monodispersity as an alternative to microfluidic devices. More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.6, No.2, pp. 117-128, 2007, DOI:10.3970/cmc.2007.006.117

Abstract In this paper the Landau-Lifshitz equation is subjected to a circularly polarized field in the plane, as well as both a dc field and an anisotropy field along the vertical easy axis perpendicular to the plane. The representation of Landau-Lifshitz equation in the Minkowski space is a Lie-type system. By performing a computation through the Lie-group solvers we can develop a consistent numerical method, which satisfies the consistency condition exactly, and thus can retain the invariant behavior. Then, we use the consistent numerical method to investigate the magnetization reversal, whose switching criterion is displayed through the minimum curve of the… More >

Chia-Cheng Tsai ¹

CMC-Computers, Materials & Continua, Vol.22, No.3, pp. 197-218, 2011, DOI:10.3970/cmc.2011.022.197

Abstract This paper describes the combination of the method of fundamental solutions (MFS) and the dual reciprocity method (DRM) as a meshless numerical method to solve problems of Kirchhoff plates under arbitrary loadings. In the solution procedure, a arbitrary distributed loading is first approximated by either the multiquadrics (MQ) or the augmented polyharmonic splines (APS), which are constructed by splines and monomials. The particular solutions of multiquadrics, splines and monomials are all derived analytically and explicitly. Then, the complementary solutions are solved formally by the MFS. Furthermore, the boundary conditions of lateral displacement, slope, normal moment, and effective shear force are… More >

Gu, M. H.¹, Young, D. L.^1,2, Fan, C. M.¹

CMC-Computers, Materials & Continua, Vol.11, No.3, pp. 185-208, 2009, DOI:10.3970/cmc.2009.011.185

Abstract A meshless numerical algorithm is developed for the solutions of one-dimensional wave equations in this paper. The proposed numerical scheme is constructed by the Eulerian-Lagrangian method of fundamental solutions (ELMFS) together with the D'Alembert formulation. The D'Alembert formulation is used to avoid the difficulty to constitute the linear algebraic system by using the ELMFS in dealing with the initial conditions and time-evolution. Moreover the ELMFS based on the Eulerian-Lagrangian method (ELM) and the method of fundamental solutions (MFS) is a truly meshless and quadrature-free numerical method. In this proposed wave model, the one-dimensional wave equation is reduced to an implicit… More >

J.H. Tian^1,2, X. Han², S.Y. Long³, G.Q. Xie⁴

CMC-Computers, Materials & Continua, Vol.10, No.3, pp. 229-242, 2009, DOI:10.3970/cmc.2009.010.229

Abstract A heat conduction analysis of the functionally graded material (FGM) plates has been investigated based on the hybrid numerical method (HNM). HNM combines the layer element method with the method of Fourier transforms and proves to be efficient and reliable. The FGM plates are infinite large and the material properties vary continuously through thickness. The heat source continually acted one the FGM plates. The temperature distribution of the FGM plates is obtained in different time and different position. Some useful results for heat conduction problems are shown in figures. This article applies HNM to heat conduction firstly and provides us… More >