Yunpeng Ma¹, Lifeng Wang^{1, *}, Mingxu Yi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.113, No.1, pp. 89-108, 2017, DOI:10.3970/cmes.2017.113.086

Abstract In this paper, an analytical time domain formulation based on Ffowcs Williams-Hawkings (FW-H) equation is derived for the prediction of the acoustic velocity field generated by moving bodies. This provides the imposition of the Neumann boundary condition on a rigid scattering surface. In order to calculate the scattering sound pressure of the duct, a thin-body boundary element method (BEM) has been proposed. The radiate sound pressure is calculated using the acoustic analogy FW-H equation. The scattering effect of the duct wall on the propagation of the sound wave is presented using the thin-body BEM. Computational More >

Liu Liqi¹, Wang Haitao^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.109-110, No.4, pp. 303-324, 2015, DOI:10.3970/cmes.2015.109.303

Abstract This paper presents a three-dimensional (3-D) boundary element method (BEM) scheme based on the Radial Integration Method (RIM) for wave propagation analysis of continuously non-homogeneous problems. The Kelvin fundamental solutions are adopted to derive the boundary-domain integral equation (BDIE). The RIM proposed by Gao (Engineering Analysis with Boundary Elements 2002; 26(10):905-916) is implemented to treat the domain integrals in the BDIE so that only boundary discretization is required. After boundary discretization, a set of second-order ordinary differential equations with respect to time variable are derived, which are solved using the Wilson-q method. Main advantages of More >

Y.C. Shiah¹, C.L. Tan^2,3, Li-Ding Chan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.109-110, No.1, pp. 15-33, 2015, DOI:10.3970/cmes.2015.109.015

Abstract In the conventional boundary element method (BEM), the presence of singular kernels in the boundary integral equation or integral identities causes serious inaccuracy of the numerical solutions when the source and field points are very close to each other. This situation occurs commonly in elastostatic analysis of thin structures. The numerical inaccuracy issue can be resolved by some regularization process. Very recently, the self-regularization scheme originally proposed by Cruse and Richardson (1996) for 2D stress analysis has been extended and modified by He and Tan (2013) to 3D elastostatics analysis of isotropic bodies. This paper More >

Jun Huang^1,2, Chaopu Zhang¹, Song Xiang², Liu Yang¹, Mingxu Yi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.108, No.5, pp. 315-330, 2015, DOI:10.3970/cmes.2015.108.315

Abstract In order to accurately predict the aerodynamic noise of the propeller, a hybrid method combining Computational Fluid Dynamics (CFD) method with Boundary Element Method (BEM) is developed in this paper. The calculation includes two steps: firstly, the unsteady viscous flow around the propeller is calculated using the CFD method to acquire the noise source information; secondly, the radiated sound pressure is calculated using BEM method in the frequency domain. In comparison with the experimental results from wind tunnel, the calculated results of aerodynamic performance are rather desirable. The simulation and experimental results of aerodynamic noise More >

T. Gortsas¹, S.V. Tsinopoulos², D. Polyzos^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.107, No.4, pp. 321-343, 2015, DOI:10.3970/cmes.2015.107.321

Abstract An advanced Boundary Element method (BEM) accelerated via Adaptive Cross Approximation (ACA) and Hierarchical Matrices (HM) techniques is presented for the solution of large-scale elastostatic problems with multi-connected domains like in fiber reinforced composite materials. Although the proposed ACA/ BEM is demonstrated for two-dimensional (2D) problems, it is quite general and it can be used for 3D problems. Different forms of ACA technique are employed for exploring their efficiency when they combined with a BEM code. More precisely, the fully and partially pivoted ACA with and without recompression are utilized, while the solution of the More >

J. Useche¹, H. Alvarez¹

CMES-Computer Modeling in Engineering & Sciences, Vol.107, No.4, pp. 277-296, 2015, DOI:10.3970/cmes.2015.107.277

Abstract Dynamic stress analysis of laminated composites plates represents a relevant task in designing of aerospace, shipbuilding and automotive components where impulsive loads can lead to sudden structural failure. The mechanical complexity inherent to these kind of components makes the numerical modeling an essential engineering analysis tool. This work deals with dynamic analysis of stresses and deformations in laminated composites thick plates using a new Boundary Element Method formulation. Composite laminated plates were modeled using the Reissner’s plate theory. We propose a direct time-domain formulation based on elastostatic fundamental solution for symmetrical laminated thick plates. Formulation More >

Chengxi Li¹

CMES-Computer Modeling in Engineering & Sciences, Vol.107, No.2, pp. 125-153, 2015, DOI:10.3970/cmes.2015.107.125

Abstract Applying an efficient numerical wave absorption method is very important for realization of an open sea condition especially for long time numerical fluid structure interaction simulation. This paper proposed a hybrid numerical wave absorption method for the fully nonlinear fluid structure simulation. A numerical code “QBEM” which is based on the quadratic boundary element method is applied to evaluate the efficiency of the wave absorption methods and damping schemes by checking the energy conservation and wave elevation in the computational domain. Specifically, we conduct a 3D numerical wave tank experiment to find the bestperforming damping More >

Z. Karakaya¹, B. Baranoğlu², B. Çetin³, A. Yazici⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.104, No.3, pp. 227-249, 2015, DOI:10.3970/cmes.2015.104.227

Abstract A new formulation for tracking multiple particles in slow viscous flow for microfluidic applications is presented. The method employs the manipulation of the boundary element matrices so that finally a system of equations is obtained relating the rigid body velocities of the particle to the forces applied on the particle. The formulation is specially designed for particle trajectory tracking and involves successive matrix multiplications for which SMP (Symmetric multiprocessing) parallelisation is applied. It is observed that present formulation offers an efficient numerical model to be used for particle tracking and can easily be extended for More >

Xiaomin An¹, Min Xu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.98, No.6, pp. 601-629, 2014, DOI:10.3970/cmes.2014.098.601

Abstract Nonlinear aeroelasticity, caused by the interaction between nonlinear fluid and geometrically nonlinear structure, is studied by an improved CFD and CSD coupled program. An AUSMpw+ flux splitting scheme, combined with an implicit time marching technology and geometric conservation law, is utilized to solve unsteady aerodynamic pressure; The finite element co-rotational theory is applied to model geometrically nonlinear two-dimensional and three-dimensional panels, and a predictor-corrector program with an approximately energy conservation is developed to obtain nonlinear structure response. The two solvers are connected by Farhat’s second order loosely coupled method and the aerodynamic loads and structural More >

A.E. Kampitsis¹, E.J. Sapountzakis²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.6, pp. 367-409, 2014, DOI:10.3970/cmes.2014.103.367

Abstract In this paper a Boundary Element Method (BEM) is developed for the geometrically nonlinear inelastic analysis of Timoshenko beams of arbitrary doubly symmetric simply or multiply connected constant cross-section, resting on inelastic tensionless Winkler foundation. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading, while its edges are subjected to the most general boundary conditions. To account for shear deformations, the concept of shear deformation coefficients is used. A displacement based formulation is developed and inelastic redistribution More >