LiviuMarin ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.37, No.3, pp. 203-242, 2008, DOI:10.3970/cmes.2008.037.203

Abstract In this paper, a meshless method for the stable solution of direct and inverse problems associated with the two-dimensional Laplace equation in the presence of boundary singularities and noisy boundary data is proposed. The governing equation and boundary conditions are discretized by the method of fundamental solutions (MFS), whilst the existence of the boundary singularity is taken into account by subtracting from the original MFS solution the corresponding singular solutions, as given by the asymptotic expansion of the solution near the singular point. However, even in the case when the boundary singularity is accounted for, the numerical solutions obtained by… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.3, pp. 261-286, 2008, DOI:10.3970/cmes.2008.036.261

Abstract The inverse Sturm-Liouville problem finds its applications in the identification of mechanical properties and/or geometrical configurations of a vibrating continuous medium; however, this problem is hard to solve, either theoretically or numerically. Previously, Liu (2008a) has constructed a Lie-group shooting method to determine the eigenvalues, and the corresponding eigenfunctions, for the direct Sturm-Liouville problem. In this study, we are concerned with solving the inverse Sturm-Liouville problem, by developing a Lie-group of SL(2,R) to construct nonlinear algebraic equations (NAEs), when discrete eigenvalues are specified. Our purpose here is to use these NAEs to solve the unknown function in the Sturm-Liouville operator.… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.91, No.1, pp. 17-42, 2013, DOI:10.3970/cmes.2013.091.017

Abstract In the present paper, we propose a novel two-side equilibration method to properly reduce the condition number of a given non-singular matrix only through a few operations. Then, two different conditioners together with the conjugate gradient method (CGM) are developed, which can overcome the defect of CGM, being not vulnerable to noisy disturbance exerted on an ill-posed linear system. The twoside CGM (TSCGM) and the pre-conditioning CGM (PrCGM) are convergent fast and accurate in solving linear inverse problems and the linear Hilbert problem under a large random noise. More >

W. Zhang¹, X. Han^1,2, J. Liu¹, R. Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.5, pp. 385-408, 2012, DOI:10.3970/cmes.2012.086.385

Abstract A novel inverse technique is presented for quantifying the uncertainty of the identified the results in an encounter condition identification problem. In this technique, the polynomial response surface method based on the structure-selection technique is first adopted to construct the approximation model of the projectile/target system, so as to reduce the computational cost. The Markov Chain Monte Carlo method is then used to identify the encounter condition parameters and their confidence intervals based on this cheap approximation model with Bayesian perspective. The results are demonstrated through the simulated test cases, which show the utility and efficiency of the proposed technique.… More >

Gusein Sh. Guseinov¹

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.4, pp. 301-320, 2012, DOI:10.3970/cmes.2012.086.301

Abstract This paper deals with the inverse spectral problem for two spectra of finite order complex Jacobi matrices (tri-diagonal symmetric matrices with complex entries). The problem is to reconstruct the matrix using two sets of eigenvalues, one for the original Jacobi matrix and one for the matrix obtained by replacing the first diagonal element of the Jacobi matrix by some another number. The uniqueness and existence results for solution of the inverse problem are established and an explicit algorithm of reconstruction of the matrix from the two spectra is given. 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 >

Chein-Shan Liu

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.5, pp. 499-526, 2012, DOI:10.3970/cmes.2012.083.499

Abstract The iterative algorithms based on the concept of best vector are proposed to solve an ill-conditioned linear system: Bx-b=0, which might be a discretization of linear inverse problem. In terms of r:=Bx-b and a monotonically increasing positive function Q(t) of a time-like variable t, we define a future cone in the Minkowski space, wherein the discrete dynamics of the proposed algorithm is evolved. We propose two methods to approximate the best vector B-1r, and obtain three iterative algorithms for solving x, which we label them as the steepest-descent and optimal vectors iterative algorithm (SOVIA), the mixed optimal iterative algorithm (MOIA),… More >

Gongsheng Li¹, De Yao², Hengyi Jiang³, Xianzheng Jia¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.2, pp. 83-108, 2011, DOI:10.3970/cmes.2011.074.083

Abstract This paper deals with an inverse problem of determining a time-depen -dent reaction coefficient arising from a disturbed soil-column infiltrating experiment based on measured breakthrough data. A purpose of doing such experiment is to simulate and study transport behaviors of contaminants when they vertically penetrating through the soils. Data compatibility of the inverse problem is discussed showing a sufficient condition to the solution's monotonicity and positivity with the help of an adjoint problem. Furthermore, an optimal perturbation regularization algorithm is applied to solve the inverse problem, and two typical numerical examples are presented to support the inversion algorithm. Finally, transport… More >

V.C. Mariani¹, V. J. Neckel², L. S. Coelho³

CMES-Computer Modeling in Engineering & Sciences, Vol.68, No.2, pp. 167-184, 2010, DOI:10.3970/cmes.2010.068.167

Abstract In this study an inverse heat conduction problem using two optimization methods to estimate apparent thermal diffusivity at different drying temperatures is solved. Temperature and moisture versus time were obtained numerically using heat and mass transfer equations with drying temperatures in the range between 20°C to 70°C. The solution of the partial differential equation is made with a finite difference method coupled to optimization techniques of Differential Evolution (DE) and Particle Swarm Optimization (PSO) used in inverse problem. Statistical analysis shows no significant differences between reported and estimated curves, and no remarkable differences between results obtained using DE and PSO… More >

Arif Amirov¹, Fikret Gölgeleyen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.2, pp. 179-192, 2010, DOI:10.3970/cmes.2010.065.179

Abstract In this work, we derive the solvability conditions for an inverse problem for the kinetic equation and develop a new symbolic algorithm to obtain the approximate solution of the problem. The computational experiments show that proposed method provides highly accurate numerical solutions even subjecting to a large noise in the given data. More >