H. Nguyen-Van¹, N.Mai-Duy¹, T. Tran-Cong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.1, pp. 65-96, 2008, DOI:10.3970/cmes.2008.036.065

Abstract A novel node-based smoothing element for triangular and quadrilateral meshes is presented for static analysis of planar piezoelectric structures. In contrast to the smoothed finite element formulation that was based on sub-cells within an original quadrilateral element, this new method transforms a general original finite element mesh into a mesh of new smoothing cells individually associated with a single node which is termed as node-based elements. The displacement fields of the element are approximated by the linear interpolation functions of the original mesh while the approximations of mechanical strains and electric potential fields are normalized using the stabilized conforming nodal… More >

A. P. S. Selvadurai¹, Wenjun Dong

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 41-54, 2006, DOI:10.3970/cmes.2006.012.041

Abstract A Fourier analysis conducted on both the spatial and the temporal discretizations of the governing partial differential equation shows that the Courant number as well as the time marching scheme have significant influences on the numerical behaviour of a Modified Least Squares (MLS) method for the solution of the advection equation. The variations of the amplification factor and the relative phase velocity with the Courant number and the dimensionless wave number indicate that when Courant number is equal to unity, the MLS method with the specified time-weighting and upwind function gives accurate results. This conclusion is confirmed by the numerical… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.3, pp. 171-198, 2012, DOI:10.3970/cmes.2012.086.171

Abstract In the optimal control theory, the Hamiltonian formalism is a famous one to find an optimal solution. However, when the performance index is complicated or for a degenerate case with a non-convexity of the Hamiltonian function with respect to the control force the Hamiltonian method does not work to find the solution. In this paper we will address this important issue via a quite different approach, which uses the optimal control problem of nonlinear Duffing oscillator as a demonstrative example. The optimally controlled vibration problem of nonlinear oscillator is recast into a nonlinear inverse problem by identifying the unknown heat… More >

M.G. Armentano¹, C. Padra², R. Rodríguez³, M. Scheble²

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.4, pp. 359-382, 2012, DOI:10.3970/cmes.2012.084.359

Abstract In this paper we propose an hp finite element method to solve a two-dimensional fluid-structure vibration problem. This problem arises from the computation of the vibration modes of a bundle of parallel tubes immersed in an incompressible fluid. We use a residual-type a posteriori error indicator to guide an hp adaptive algorithm. Since the tubes are allowed to be different, the weak formulation is a non-standard generalized eigenvalue problem. This feature is inherited by the algebraic system obtained by the discretization process. We introduce an algebraic technique to solve this particular spectral problem. We report several numerical tests which allow… More >

M.R. Myers¹, A.B. Jorge², D.E. Yuhas³, D.G. Walker¹

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.3, pp. 221-248, 2012, DOI:10.3970/cmes.2012.083.221

Abstract An adaptive extended Kalman filter is developed and investigated for a transient heat transfer problem in which a high heat flux spot source is applied on one side of a thin plate and ultrasonic pulse time of flight is measured between spatially separated transducers on the opposite side of the plate. The novel approach is based on the uncertainty in the state model covariance and leverages trends in the extended Kalman filter covariance to drive changes to the state model covariance during convergence. This work is an integral part of an effort to develop a system capable of locating the… 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 details of an automatic local… 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 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 >