Vadim V. Romanuke¹

CMES-Computer Modeling in Engineering & Sciences, Vol.108, No.2, pp. 113-134, 2015, DOI:10.3970/cmes.2015.108.113

Abstract The problem of solving infinite noncooperative games approximately is considered. The game may either have solution or have no solution. The existing solution may be unknown as well. Therefore, an approach of obtaining the approximate solution of the infinite noncooperative game on the unit hypercube is suggested. The unit-hypercubic game isomorphism to compact infinite noncooperative games allows to disseminate the approximation approach on a pretty wide class of noncooperative games. The approximation intention is in converting the infinite game into a finite one, whose solution methods are easier rather than solving infinite games. The conversion starts with sampling the players’… More >

Marcin Maździarz¹, Marcin Gajewski²

CMES-Computer Modeling in Engineering & Sciences, Vol.105, No.2, pp. 123-150, 2015, DOI:10.3970/cmes.2015.105.123

Abstract In this paper, the choice and parametrisation of finite deformation polyconvex isotropic hyperelastic models to describe the behaviour of a class of defect-free monocrystalline metal materials at the molecular level is examined. The article discusses some physical, mathematical and numerical demands which in our opinion should be fulfilled by elasticity models to be useful. A set of molecular numerical tests for aluminium and tungsten providing data for the fitting of a hyperelastic model was performed, and an algorithm for parametrisation is discussed. The proposed models with optimised parameters are superior to those used in non-linear mechanics of crystals. More >

K. Kakuda¹, Y. Hayashi¹, J. Toyotani¹

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.4, pp. 229-249, 2014, DOI:10.3970/cmes.2014.103.229

Abstract In this paper, we present the application of the particle-based simulations to complicated fluid flow problem with free surfaces. The particle approach is based on the MPS (Moving Particle Simulation) method using hyperbolic-type weighting function to stabilize the spurious oscillatory solutions for solving the Poisson equation with respect to the pressure fields. The hyperbolic-type weighting function is constructed by differentiating the characteristic function based on neural network framework. The weighting function proposed herein is collaterally applied to the kernel function in the SPH-framework. Numerical results demonstrate the workability and validity of the present MPS approach through the dambreaking flow problem. More >

H. S. Shukla¹, Mohammad Tamsir¹, Vineet K. Srivastava², Jai Kumar³

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.1, pp. 1-17, 2014, DOI:10.3970/cmes.2014.103.001

Abstract The objective of this article is to carry out an approximate analytical solution of the time fractional order Cauchy-reaction diffusion equation by using a semi analytical method referred as the fractional-order reduced differential transform method (FRDTM). The fractional derivative is illustrated in the Caputo sense. The FRDTM is very efficient and effective powerful mathematical tool for solving wide range of real world physical problems by providing an exact or a closed approximate solution of any differential equation arising in engineering and allied sciences. Four test numerical examples are provided to validate and illustrate the efficiency of FRDTM. More >

George A. Gravvanis¹, Christos K. Filelis-Papadopoulos¹, Paschalis I.Matskanidis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.4, pp. 323-345, 2014, DOI:10.3970/cmes.2014.100.323

Abstract Since the introduction of the Algebraic MultiGrid algorithm (AMG) over twenty years ago, significant progress has been made in improving the coarsening and the convergence behavior of the method. In this paper, an AMG method is introduced that utilizes a new generic approximate inverse algorithm as a smoother in conjunction with common coarsening techniques, such as classical Ruge-Stüben coarsening, CLJP and PMIS coarsening. The proposed approximate inverse scheme, namely Generic Approximate Banded Inverse (GenAbI), is a banded approximate inverse based on Incomplete LU factorization with zero fill–in (ILU(0)). The new class of Generic Approximate Banded Inverse can be computed for… More >

Xiaokui Yue¹, Yong Xu², Jianping Yuan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.4, pp. 277-291, 2013, DOI:10.3970/cmes.2013.093.277

Abstract In this paper, a stochastic averaging method is derived for a class of non-linear stochastic systems under parametrical Poisson white noise excitation, which may be used to model the beam-beam interaction models in particle accelerators. The averaged Generalized Fokker-Planck equation is derived and the approximate stationary solution of the averaged Generalized Fokker-Planck equation is solved by using perturbation method. The present method applied in this paper can reduce the dimensions of stochastic ODE from 2n to n, which simplify the complex stochastic ODE, and then the analytical stationary solutions can be obtained. An example is employed to demonstrate the procedure… More >

P. K. Sahu¹, S. Saha Ray^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.2, pp. 91-112, 2013, DOI:10.3970/cmes.2013.093.091

Abstract In this paper, nonlinear integral equations have been solved numerically by using B-spline wavelet method and Variational Iteration Method (VIM). Compactly supported semi-orthogonal linear B-spline scaling and wavelet functions together with their dual functions are applied to approximate the solutions of nonlinear Fredholm integral equations of second kind. Comparisons are made between the variational Iteration Method (VIM) and linear B-spline wavelet method. Several examples are presented to compare the accuracy of linear B-spline wavelet method and Variational Iteration Method (VIM) with their exact solutions. More >

N Pimprikar¹, J Teresa², D Roy^1,3, R M Vasu⁴, K Rajan⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.92, No.1, pp. 33-61, 2013, DOI:10.3970/cmes.2013.092.033

Abstract Nearly pollution-free solutions of the Helmholtz equation for k-values corresponding to visible light are demonstrated and verified through experimentally measured forward scattered intensity from an optical fiber. Numerically accurate solutions are, in particular, obtained through a novel reformulation of the H1 optimal Petrov-Galerkin weak form of the Helmholtz equation. Specifically, within a globally smooth polynomial reproducing framework, the compact and smooth test functions are so designed that their normal derivatives are zero everywhere on the local boundaries of their compact supports. This circumvents the need for a priori knowledge of the true solution on the support boundary and relieves the… More >

Christos K. Filelis-Papadopoulos¹, George A. Gravvanis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.90, No.3, pp. 233-253, 2013, DOI:10.3970/cmes.2013.090.232

Abstract During the last decades, multigrid methods have been extensively used in order to solve large scale linear systems derived from the discretization of partial differential equations using the finite difference method. Approximate Inverses in conjunction with Richardon’s iterative method could be used as smoothers in the multigrid method. Thus, a new class of smoothers based on approximate inverses could be derived. Effectiveness of explicit approximate inverses relies in the fact that they are close approximants to the inverse of the coefficient matrix and are fast to compute in parallel. Furthermore, the class of finite difference approximate inverses proposed in conjunction… More >

G.A. Gravvanis¹, K.M. Giannoutakis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.2, pp. 69-82, 2006, DOI:10.3970/cmes.2006.016.069

Abstract A new class of normalized explicit optimized approximate inverse finite element matrix techniques, based on normalized finite element approximate factorization procedures, for solving sparse linear systems resulting from the finite element discretization of partial differential equations in three space variables are introduced. A new parallel normalized explicit preconditioned conjugate gradient square method in conjunction with normalized approximate inverse finite element matrix techniques for solving efficiently sparse finite element linear systems on distributed memory systems is also presented along with theoretical estimates on speedups and efficiency. The performance on a distributed memory machine, using Message Passing Interface (MPI) communication library, is… More >