V.G.Yakhno¹, T.M. Yakhno²

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 233-250, 2011, DOI:10.3970/cmes.2011.080.233

Abstract The time-dependent Maxwell's equations in non-dispersive homogeneous bi-anisotropic materials are considered in the paper. These equations are written as a symmetric hyperbolic system. A new method of the computation of the electric and magnetic fields arising from electric current is suggested in the paper. This method consists of the following. The Maxwell's equations are written in terms of the Fourier transform with respect to the space variables. The Fourier image of the obtained system is a system of ordinary differential equations whose coefficients depend on the 3D Fourier parameter. The formula for the solution of More >

M.M. Kaminski

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.2, pp. 113-140, 2011, DOI:10.3970/cmes.2011.080.113

Abstract The main aim of this work is to demonstrate a solution to the transient problems for the statistically homogeneous media with random physical parameters. This is done with the use of the stochastic perturbation technique based on the general order Taylor series expansions and the additionally modified implementation of the Finite Element Method. Now, both the Direct Differentiation Method as well as the Response Function Method are employed to form and solve up to the nth order state equations. Computational implementation of both approaches is illustrated using two examples - by determination of the probabilistic More >

A. Khosravifard¹, M.R. Hematiyan^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.3&4, pp. 185-208, 2011, DOI:10.3970/cmes.2011.078.185

Abstract An inverse method for optimal control of the freezing front motion in the solidification of pure materials is presented. The inverse technique utilizes the idea of a pseudo heat source to account for the latent heat effects. The numerical formulation of this inverse method is based on a formerly introduced meshless technique. In this method, the flux and the velocity of the liquid-solid interface are treated as secondary variables and the liquid and solid domains are modeled simultaneously. Some numerical examples are provided to demonstrate the efficiency of the presented method. The effects of regularization More >

Chein-Shan Liu¹, Chung-Lun Kuo²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.1, pp. 57-80, 2011, DOI:10.3970/cmes.2011.077.057

Abstract In this paper, the solutions of inverse Cauchy problems for quasi-linear elliptic equations are resorted to an unusual mixed group-preserving scheme (MGPS). The bottom of a finite rectangle is imposed by overspecified boundary data, and we seek unknown data on the top side. The spring-damping regularization method (SDRM) is introduced by converting the governing equation into a new one, which includes a spring term and a damping term. The SDRM can further stabilize the inverse Cauchy problems, such that we can apply a direct numerical integration method to solve them by using the MGPS. Several More >

A. Rezaei Mojdehi^1,2, A. Darvizeh³, A. Basti²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.1, pp. 1-32, 2011, DOI:10.3970/cmes.2011.077.001

Abstract In this paper, a meshless method based on the local petrov-galerkin approach is proposed for the three dimensional (3D) elasto-plastic problems. Galerkin weak-form formulation is applied to derive the discrete governing equations. A weak formulation for the set of governing equations is transformed into local integral equations on local sub-domains by using a unit test function. Nodal points are distributed in the 3D analyzed domain and each node is surrounded by a cubic sub-domain to which a local integral equation is applied. Three dimensional Moving Least-Square (MLS) approximation is used as shape function to approximate 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 >

K. Jayabal¹, A. Menzel^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.73, No.2, pp. 183-208, 2011, DOI:10.3970/cmes.2011.073.183

Abstract Naturally evolving Voronoi discretisation of two-dimensional plane domains renders representative microstructures that turn out to be useful for the modelling and simulation of polycrystalline materials. Hybrid finite element approaches are employed on such polygonal discretisations to solve, for instance, mechanical and electromechanical problems within a finite element context. In view of solving mechanical problems, varying order of polynomial functions are suggested in the literature to sufficiently approximate stresses within the polygonal finite elements. These are, in addition to the order of the approximation functions for the displacements, characterised by the number of edges in the… More >

D.-A. An-Vo¹, N. Mai-Duy¹, T. Tran-Cong¹

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.4, pp. 299-336, 2011, DOI:10.3970/cmes.2011.072.299

Abstract This paper presents a new control-volume discretisation method, based on Cartesian grids and integrated-radial-basis-function elements (IRBFEs), for the solution of second-order elliptic problems in one and two dimensions. The governing equation is discretised by means of the control-volume formulation and the division of the problem domain into non-overlapping control volumes is based on a Cartesian grid. Salient features of the present method include (i) an element is defined by two adjacent nodes on a grid line, (ii) the IRBF approximations on each element are constructed using only two RBF centres (a smallest RBF set) associated… More >

Veturia Chiroiu¹, Ligia Munteanu², Ioan Ursu³

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.3, pp. 229-246, 2011, DOI:10.3970/cmes.2011.072.229

Abstract Chaotic behavior of uncertain nonlinear systems offers a rich variety of orbits, which can be controlled by bounding the signals involved in closed-loop systems. In this paper, systems with nonlinear uncertainties with no prior knowledge of their bounds, unmodeled dynamic law and rapidly varying disturbances are analyzed in order to propose a stabilization controller of the chaotic behavior via the fuzzy logic systems. More >

H. T. Wang^1,2, Z. H. Yao³

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.2, pp. 115-148, 2011, DOI:10.3970/cmes.2011.072.115

Abstract A fast boundary element solver for the analysis of three-dimensional general crack problems is presented. In order to effectively model the embedded or edge cracked structures a dual boundary integral equation (BIE) formulation is used. By implementing the fast multipole method (FMM) to the discretized BIE, structures containing a large number of three-dimensional cracks can be readily simulated on one personal computer. In the FMM framework, a multipole expansion formulation is derived for the hyper-singular integral in order that the multipole moments of the dual BIEs containing the weakly-, strongly- and hyper-singular kernels are collected More >