N.A. Dumont¹, E.Y. Mamani¹

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 109-150, 2011, DOI:10.3970/cmes.2011.078.109

Abstract A particular implementation of the hybrid boundary element method is presented for the two-dimensional analysis of potential and elasticity problems, which, although general in concept, is suited for fracture mechanics applications. Generalized Westergaard stress functions, as proposed by Tada, Ernst and Paris in 1993, are used as the problem's fundamental solution. The proposed formulation leads to displacement-based concepts that resemble those presented by Crouch and Starfield, although in a variational framework that leads to matrix equations with clear mechanical meanings. Problems of general topology, such as in the case of unbounded and multiply-connected domains, may More >

Y.C. Shiah¹, C. L. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.2, pp. 95-108, 2011, DOI:10.3970/cmes.2011.078.095

Abstract By differentiating the Green function of Ting and Lee (1997) for 3D general anisotropic elastotatics in a spherical coordinate system as an intermediate step, and then using the chain rule, derivatives of up to the second order of this fundamental solution are obtained in exact, explicit, algebraic forms. No tensors of order higher than two are present in these derivatives, thereby allowing these quantities to be numerically evaluated quite expeditiously. These derivatives are required for the computation of the internal point displacements and stresses via Somigliana's identity in BEM analysis. Some examples are presented to More >

A. Sellier, N. Ghalia

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.1, pp. 25-50, 2011, DOI:10.3970/cmes.2011.078.025

Abstract The Green tensor complying with anisotropic slip conditions at the surface of a plane, impermeable, motionless and slipping wall is theoretically obtained and an efficient numerical method is proposed to accurately compute at a very reasonable cpu time cost each of its Cartesian component. The accuracy of the advocated numerical strategy is tested against the Maple Software and the employed procedure makes it possible to calculate the Green tensor for a non-isotropic slip condition at a cpu time cost comparable with the one needed for the less complicated isotropic Navier condition. More >

A. Brancati¹, M. H. Aliabadi¹, A. Milazzo²

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.1, pp. 1-24, 2011, DOI:10.3970/cmes.2011.078.001

Abstract This paper presents an improved adaptive cross approximation (ACA) approach developed in conjunction with the Hierarchical format matrix and the GMRES solver. A novel scheme to generate the cluster tree (based upon preliminary considerations of the prescribed boundary conditions) and an improved ACA algorithm (approximating the system matrix for mixed Robin conditions) are described. The asymptotic smoothness property of a kernel generated by a linear combination of two asymptotic smooth kernels is demonstrated. Numerical results show the new approach to be up to 50% faster than the conventional ACA approach. More >

V. Ungvichian¹, P. Kanongchaiyos¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 201-220, 2011, DOI:10.3970/cmes.2011.077.201

Abstract This paper describes an enhancement to Garland and Heckbert's mesh simplification method by using the principal curvatures and directions of each vertex. We calculate the values and directions, before using them to determine the absolute normal curvature in the direction of contraction, and multiplying the curvature with the edge length, the maximum absolute cosine of the angles between the edge and the normals of faces adjacent to either endpoint, and the quadric error of the collapse. We also apply penalties based on compactness and angular and dihedral deviations of the resulting faces. We have implemented… More >

M.Souli¹, F.Erchiqui²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 183-200, 2011, DOI:10.3970/cmes.2011.077.183

Abstract During the design process of membrane structure to resist to high pressure loading, and the characterization of hyperelastic material, a structure made up of thin rubber undergoes large deformation and rotation under high pressure loading out of high pressurized gas. Until recently, to simulate the inflation of the hyperelastic membrane, a uniform pressure based on thermodynamic model or experimental tests is applied to the structure, as boundary conditions. From a computational time point of view, this approach is very fast, since no computational fluid dynamics is involved in the simulation. However, at the late stage… More >

Minghui Wang¹, Musheng Wei², Shanrui Hu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 173-182, 2011, DOI:10.3970/cmes.2011.077.173

Abstract The mapping from the symmetric solution set to its independent parameter space is studied and an iterative method is proposed for the least-squares minimum-norm symmetric solution of AXB = E. Numerical results are reported that show the efficiency of the proposed methods. More >

Hsin-Ping Chu¹, Cheng-Ying Lo²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 161-172, 2011, DOI:10.3970/cmes.2011.077.161

Abstract This paper presents the application of the differential transform method to solve strongly nonlinear equations with cubic nonlinearities and self-excitation terms. First, the equations are transformed by the differential transform method into the algebra equations in terms of the transformed functions. Secondly, the higher-order transformed functions are calculated in terms of other lower-order transformed functions through the iterative procedure. Finally, the solutions are approximated by the n-th partial sum of the infinite series obtained by the inverse differential transform. Two strongly nonlinear equations with different coefficients and initial conditions are given as illustrative examples. More >

Zhou YH^1,2, Wang XM², Wang JZ^1,2 , Liu XJ²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 137-160, 2011, DOI:10.3970/cmes.2011.077.137

Abstract In this paper, we present an efficient wavelet-based algorithm for solving a class of fractional vibration, diffusion and wave equations with strong nonlinearities. For this purpose, we first suggest a wavelet approximation for a function defined on a bounded interval, in which expansion coefficients are just the function samplings at each nodal point. As the fractional differential equations containing strong nonlinear terms and singular integral kernels, we then use Laplace transform to convert them into the second type Voltera integral equations with non-singular kernels. Certain property of the integral kernel and the ability of explicit More >

A.B. Tatomir^1,2, A.Szymkiewicz³, H. Class¹, R. Helmig¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 81-112, 2011, DOI:10.3970/cmes.2011.077.081

Abstract We present a two phase flow conceptual model, the corresponding simulator (2pMINC) and a workflow for large-scale fractured reservoirs, based on a continuum fracture approach which uses the multiple interacting continua (MINC) method complemented with an improved upscaling technique. The complex transient behavior of the flow processes in fractured porous media is captured by subgridding the coarse blocks in nested volume elements which have effective properties calculated from the detailed representation of the fracture system. In this way, we keep a physically based approach, preserve the accuracy of the model, avoid the common use of… More >