Xiao Zhao¹, Haitian Yang^1,2, Qiang Gao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.2, pp. 139-160, 2011, DOI:10.3970/cmes.2011.074.139

Abstract A piecewised adaptive algorithm in the time domain is presented to solve the transient convection-diffusion heat transfer problem. By expanding all variables at a time interval, an initial and boundary value problem is decoupled into a series of recursive boundary value problems which can be solved by FEM or other well developed numerical schemes to deal with boundary value problems. A steady computing accuracy can be adaptively maintained via the power increase of the expansion, particularly when the step size varies in the whole computing process. Additionally for the nonlinear cases, there is no requirement of iteration and additional assumption… More >

A. Gakwaya¹, H. Sharifi², M. Guillot¹, M. Souli³, F. Erchiqui⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.3, pp. 209-266, 2011, DOI:10.3970/cmes.2011.073.209

Abstract A mechanical computer aided design and engineering system can be used to reduce the design-to-manufacture cycle time in metal forming process. Such a system could be built upon a solid modeling geometry engine and an efficient finite element (FE) solver. The maintenance of a high-quality mesh throughout the analysis is an essential feature of an efficient finite element simulation of large strain metal forming problems. In this paper, a mesh adaptation technique employing the Arbitrary Lagrangian-Eulerian formulation (ALE) is applied to some industrial metal forming problems. An ACIS boundary representation of the solid model is employed. This type of representation… More >

David B. Carrington¹, Xiuling Wang², Darrell W. Pepper³

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.1, pp. 23-44, 2010, DOI:10.3970/cmes.2010.061.023

Abstract A two-equation turbulence closure model (k-ω) using an h-adaptive grid technique and finite element method (FEM) has been developed to simulate low Mach flow and heat transfer. These flows are applicable to many flows in engineering and environmental sciences. Of particular interest in the engineering modeling areas are: combustion, solidification, and heat exchanger design. Flows for indoor air quality modeling and atmospheric pollution transport are typical types of environmental flows modeled with this method. The numerical method is based on a hybrid finite element model using an equal-order projection process. The model includes thermal and species transport, localized mesh refinement… More >

J. Trevelyan¹and G. Coates¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 147-170, 2010, DOI:10.3970/cmes.2010.055.147

Abstract The terminology "wave boundary elements" relates to boundary elements enriched in the Partition of Unity sense by a multiple plane wave basis for the analysis of the propagation of short wavelength waves. This paper presents a variant of this approach in which the plane wave basis is selected adaptively according to an error indicator. The error indicator is residual based, and exhibits useful local and global properties. Model improvement in each adaptive iteration is carried out by the addition of new plane waves with no h-refinement. The convergence properties of the scheme are demonstrated. More >

Nicolas Ali Libre^1,2, Arezoo Emdadi², Edward J. Kansa^3,4, Mohammad Shekarchi², Mohammad Rahimian²

CMES-Computer Modeling in Engineering & Sciences, Vol.50, No.2, pp. 161-190, 2009, DOI:10.3970/cmes.2009.050.161

Abstract We present a wavelet based adaptive scheme and investigate the efficiency of this scheme for solving nearly singular potential PDEs over irregularly shaped domains. For a problem defined over Ω∈ℜ^d, the boundary of an irregularly shaped domain, Γ, is defined as a boundary curve that is a product of a Heaviside function along the normal direction and a piecewise continuous tangential curve. The link between the original wavelet based adaptive method presented in Libre, Emdadi, Kansa, Shekarchi, and Rahimian (2008, 2009) or LEKSR method and the generalized one is given through the use of simple Heaviside masking procedure. In addition… More >

Alexandre de Castro^1,2, Carlos Frederico Fronza², Domingos Alves^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.45, No.1, pp. 83-96, 2009, DOI:10.3970/cmes.2009.045.083

Abstract Immunological systems have been an abundant inspiration to contemporary computer scientists. Problem solving strategies, stemming from known immune system phenomena, have been successfully applied to chall enging problems of modern computing. Simulation systems and mathematical modeling are also beginning use to answer more complex immunological questions as immune memory process and duration of vaccines, where the regulation mechanisms are not still known sufficiently (Lundegaard, Lund, Kesmir, Brunak, Nielsen, 2007). In this article we studiedin machinaa approach to simulate the process of antigenic mutation and its implications for the process of memory. Our results have suggested that the durability of the… More >

A. Brancati¹, M. H. Aliabadi¹, I. Benedetti^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.2, pp. 149-172, 2009, DOI:10.3970/cmes.2009.043.149

Abstract In this paper a new Rapid Acoustic Boundary Element Method (RABEM) is presented using a Hierarchical GMRES solver for 3D acoustic problems. The Adaptive Cross Approximation is used to generate both the system matrix and the right hand side vector. The ACA is also used to evaluate the potential and the particle velocity values at selected internal points. Two different GMRES solution strategies (without preconditioner and with a block diagonal preconditioner) are developed and tested for low and high frequency problems. Implementation of different boundary conditions (i.e. Dirichlet, Neumann and mixed Robin) is also described. The applications presented include the… More >

G. C. Bourantas¹, E. D. Skouras^2,3,4, G. C. Nikiforidis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.1, pp. 1-26, 2009, DOI:10.3970/cmes.2009.043.001

Abstract The extent of application of meshfree methods based on point collocation (PC) techniques with adaptive support domain for strong form Partial Differential Equations (PDE) is investigated. The basis functions are constructed using the Moving Least Square (MLS) approximation. The weak-form description of PDEs is used in most MLS methods to circumvent problems related to the increased level of resolution necessary near natural (Neumann) boundary conditions (BCs), dislocations, or regions of steep gradients. Alternatively, one can adopt Radial Basis Function (RBF) approximation on the strong-form of PDEs using meshless PC methods, due to the delta function behavior (exact solution on nodes).… More >

Shang Ma¹, Xiong Zhang^1,2, Yanping Lian¹, Xu Zhou³

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.2, pp. 101-124, 2009, DOI:10.3970/cmes.2009.039.101

Abstract Numerical simulation of high explosive explosion problems is a big challenge to traditional numerical methods because explosion usually involves extremely large deformation and multi-material interaction of different phases. Recently developed meshfree methods show much advantages over mesh-based method for problems associated with very large deformation. Some of them have been successfully applied to impact and explosion problems, such as smoothed particle hydrodynamics (SPH). Similar to SPH, material point method (MPM) is an efficient meshfree particle method solving continuum problems. With combination of the advantages of Eulerian and Lagrangian methods, MPM is a promising numerical tool for solving large deformation problems,… More >

Nicolas Ali Libre^1,2, Arezoo Emdadi², Edward J. Kansa^3,4, Mohammad Shekarchi², Mohammad Rahimian²

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.3, pp. 263-284, 2008, DOI:10.3970/cmes.2008.038.263

Abstract We present a wavelet based adaptive scheme and investigate the efficiency of this scheme for solving nearly singular potential PDEs. Multiresolution wavelet analysis (MRWA) provides a firm mathematical foundation by projecting the solution of PDE onto a nested sequence of approximation spaces. The wavelet coefficients then were used as an estimation of the sensible regions for node adaptation. The proposed adaptation scheme requires negligible calculation time due to the existence of the fast Discrete Wavelet Transform (DWT). Certain aspects of the proposed adaptive scheme are discussed through numerical examples. It has been shown that the proposed adaptive scheme can detect… More >