Chein-Shan Liu¹, Chih-Wen Chang²

CMC-Computers, Materials & Continua, Vol.25, No.2, pp. 107-134, 2011, DOI:10.3970/cmc.2011.025.107

Abstract We consider two inverse problems for estimating radiative coefficients α(x) and α(x, y), respectively, in T_t(x, t) = T_xx(x, t)-α(x)T(x, t), and T_t(x, y, t) = T_xx(x, y, t) + T_yy(x, y, t)-α(x, y)T(x, y, t), where a are assumed to be continuous functions of space variables. A Lie-group adaptive method is developed, which can be used to find a at the spatially discretized points, where we only utilize the initial condition and boundary conditions, such as those for a typical direct problem. This point is quite different from other methods, which need the overspecified final time data. Three-fold advantages can be gained More >

S. Gerace¹, K. Erhart¹, E. Divo^1,2, A. Kassab¹

CMES-Computer Modeling in Engineering & Sciences, Vol.81, No.1, pp. 35-68, 2011, DOI:10.3970/cmes.2011.081.035

Abstract The focus of this work is to demonstrate a novel approach to true CFD automation based on an adaptive Cartesian point distribution process coupled with a Meshless flow solution algorithm. As Meshless method solutions require only an underlying nodal distribution, this approach works well even for complex flow geometries with non-aligned domain boundaries. Through the addition of a so-called shadow layer of body-fitted nodes, application of boundary conditions is simplified considerably, eliminating the stair-casing issues of typical Cartesian-based techniques. This paper describes the approach taken to automatically generate the Meshless nodal distribution, along with the More >

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 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… More >

Cheng-Te Chi¹, Ming-Hsiao Lee², Wen-Hwa Chen^1,3

CMC-Computers, Materials & Continua, Vol.24, No.3, pp. 239-256, 2011, DOI:10.3970/cmc.2011.024.239

Abstract A novel three-dimensional adaptive element-free Galerkin method (EFGM) based on a uniform background grid is proposed to cope with the problems with extremely large deformation. On the basis of this uniform background grid, an interior adaptive strategy through an error estimation within the analysis domain is developed. By this interior adaptive scheme, additional adaptive nodes are inserted in those regions where the solution accuracy needs to be improved. As opposed to the fixed uniform background grid, these inserted nodes can move along with deformation to describe the particular local deformation of the structure. In addition,… More >

M. Al-Ajmi¹, A. Benjeddou²

CMC-Computers, Materials & Continua, Vol.23, No.3, pp. 265-286, 2011, DOI:10.3970/cmc.2011.023.265

Abstract A new discrete layer finite element (DLFE) is presented for electro-mechanically coupled analyses of moderately thick piezoelectric adaptive composite plates. The retained kinematics is based on layer-wise first-order shear deformation theory, and considers the plies top and bottom surfaces in-plane displacements and the plate transverse deflection as mechanical unknowns. The former are assumed in-plane Lagrange linear, while the latter is assumed in-plane full (Lagrange) quadratic; this results in a nine nodes quadrangular (Q9) DLFE. The latter is validated in free-vibrations, first numerically against ANSYS three-dimensional piezoelectric finite elements for a cantilever moderately thick aluminum plate… More >

Yihua Xiao¹, Xu Han^1,2, Dean Hu¹

CMC-Computers, Materials & Continua, Vol.23, No.1, pp. 9-34, 2011, DOI:10.3970/cmc.2011.023.009

Abstract For impact computations, it is efficient to model small and large deformation regions by Finite Element Method (FEM) and Smoothed Particle Hydrodynamics (SPH), respectively. However, it requires an effective algorithm to couple FEM and SPH calculations. To fulfill this requirement, an alternative coupling algorithm is presented in this paper. In the algorithm, the coupling between element and particle regions are achieved by treating elements as imaginary particles and applying equivalent tractions to element sides on coupling interfaces. In addition, an adaptive coupling technique is proposed based on the algorithm to improve the computational efficiency of More >

Wenqing Wang¹, Thomas Schnicke², Olaf Kolditz³

CMC-Computers, Materials & Continua, Vol.21, No.3, pp. 217-236, 2011, DOI:10.3970/cmc.2011.021.217

Abstract This work focuses on parallel finite element simulation of thermal hydraulic and mechanical (THM) coupled processes in porous media, which is a common phenomenon in geological applications such as nuclear waste repository and CO2 storage facilities. The Galerkin finite element method is applied to solve the derived partial differential equations. To deal with the coupling terms among the equations, the momentum equation is solved individually in a monolithic manner, and moreover their solving processes are incorporated into the solving processes of nonisothermal hydraulic equation and heat transport equation in a staggered manner. The computation task… More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.21, No.1, pp. 17-40, 2011, DOI:10.3970/cmc.2011.021.017

Abstract Only the left-boundary data of temperature and heat flux are used to estimate an unknown parameter function α(x) in T_t(x,t) = ∂(α(x)T_x)/∂x + h(x,t), as well as to recover the right-boundary data. When α(x) is given the above problem is a well-known inverse heat conduction problem (IHCP). This paper solves a mixed-type inverse problem as a combination of the IHCP and the problem of parameter identification, without needing to assume a function form of α(x) a priori, and without measuring extra data as those used by other methods. We use the one-step Lie-Group Adaptive Method (LGAM) More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.64, No.3, pp. 299-326, 2010, DOI:10.3970/cmes.2010.064.299

Abstract We solve the Calderón inverse conductivity problem [Calderón (1980, 2006)], for an elliptic type equation in a rectangular plane domain, to recover an unknown conductivity function inside the domain, from the over-specified Cauchy data on the bottom of the rectangle. The Calderón inverse problem exhibitsthree-fold simultaneous difficulties: ill-posedness of the inverse Cauchy problem, ill-posedness of the parameter identification, and no information inside the domain being available on the impedance function. In order to solve this problem, we discretize the whole domain into many sub-domains of finite strips, each with a small height. Thus the Calderón… More >