Darin Koblick^1,2,3, Shujing Xu⁴, Joshua Fogel⁵, Praveen Shankar¹

CMES-Computer Modeling in Engineering & Sciences, Vol.111, No.1, pp. 1-27, 2016, DOI:10.3970/cmes.2016.111.001

Abstract The Modified Picard-Chebyshev Method (MPCM) is implemented as an orbit propagation solver for a numerical optimization method that determines minimum time orbit transfer trajectory of a satellite using a series of multiple impulses at intermediate waypoints. The waypoints correspond to instantaneous impulses that are determined using a nonlinear constrained optimization routine, SNOPT with numerical force models for both Two-Body and J2 perturbations. It is found that using the MPCM increases run-time performance of the discretized lowthrust optimization method when compared to other sequential numerical solvers, such as Adams-Bashforth-Moulton and Gauss-Jackson 8th order methods. More >

Weiwei Zhang¹, Yufeng Nie¹, Li Cai¹, Nan Qi²

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.1, pp. 55-82, 2014, DOI:10.3970/cmes.2014.102.055

Abstract In this paper, we derive an adaptive mesh generation method for discretizing the incompressible flow using node-based local grids. The flow problem is described by the Stokes equations which are solved by a stabilized low-order P1-P1 (linear velocity, linear pressure) mixed finite element method. The proposed node-based adaptive mesh generation method consists of four components: mesh size modification, a node placement procedure, a node-based local mesh generation strategy and an error estimation technique, which are combined so as to guarantee obtaining a conforming refined/coarsened mesh. The nodes are considered as particles with interaction forces, which are generated by dynamic simulation… More >

G Devaraj¹, Shashi Narayan¹, Debasish Roy^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.1, pp. 1-54, 2014, DOI:10.3970/cmes.2014.102.001

Abstract This work sets forth a 'hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (C^p-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background… More >

Sanderson L. Gonzaga de Oliveira¹, Mauricio Kischinhevsky², João Manuel R. S. Tavares³

CMES-Computer Modeling in Engineering & Sciences, Vol.95, No.2, pp. 119-141, 2013, DOI:10.3970/cmes.2013.095.119

Abstract A novel graph-based adaptive mesh refinement technique for triangular finite-volume discretizations in order to solve second-order partial differential equations is described. Adaptive refined meshes are built in order to solve timedependent problems aiming low computational costs. In the approach proposed, flexibility to link and traverse nodes among neighbors in different levels of refinement is admitted; and volumes are refined using an approach that allows straightforward and strictly local update of the data structure. In addition, linear equation system solvers based on the minimization of functionals can be easily used; specifically, the Conjugate Gradient Method. Numerical and analytical tests were carried… More >

Jianxin Zhu¹, Zheqi Shen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.4, pp. 347-362, 2011, DOI:10.3970/cmes.2011.071.347

Abstract It is known that the perfectly matched layer (PML) is a powerful tool to truncate the unbounded domain. Recently, the PML technique has been introduced in the computation of nonlinear Schrödinger equations (NSE), in which the nonlinearity is separated by some efficient time-splitting methods. A major task in the study of PML is that the original equation is modified by a factor c which varies fast inside the layer. And a large number of grid points are needed to capture the profile of c in the discretization. In this paper, the possibility is discussed for using some nonuniform finite difference… More >

Delfim Soares Jr.¹

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.3, pp. 227-248, 2010, DOI:10.3970/cmes.2010.066.227

Abstract In this work, a meshless method based on the local Petrov-Galerkin approach is proposed for the solution of pore-dynamic problems considering elastic and elastoplastic materials. Formulations adopting the Heaviside step function as the test functions in the local weak form are considered. The moving least-square method is used for the approximation of physical quantities in the local integral equations. After spatial discretization is carried out, a nonlinear system of time-domain ordinary differential equations is obtained. This system is solved by Newmark/Newton-Raphson techniques. The present work is based on the u-p formulation and the incognita fields of the coupled analysis in… More >

Delfim Soares Jr.¹

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.2, pp. 177-200, 2010, DOI:10.3970/cmes.2010.061.177

Abstract In this work, meshless methods based on the local Petrov-Galerkin approach are employed for the time-domain dynamic analysis of porous media. For the spatial discretization of the pore-dynamic model, MLPG formulations adopting Gaussian weight functions as test functions are considered, as well as the moving least square method is used to approximate the incognita fields. For time discretization, the generalized Newmark method is adopted. The present work is based on the u-p formulation and the incognita fields of the coupled analysis in focus are the solid skeleton displacements and the interstitial fluid pore pressures. Independent spatial discretization is considered for… More >

R. Vodička¹, V. Mantič², F. París²

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.3, pp. 173-204, 2007, DOI:10.3970/cmes.2007.017.173

Abstract An original approach to solve domain decomposition problems by the symmetric Galerkin boundary element method is developed. The approach, based on a new variational principle for such problems, yields a fully symmetric system of equations. A natural property of the proposed approach is its capability to deal with nonconforming discretizations along straight and curved interfaces, allowing in this way an independent meshing of non-overlapping subdomains to be performed. Weak coupling conditions of equilibrium and compatibility at an interface are obtained from the critical point conditions of the energy functional. Equilibrium is imposed through local traction (Neumann) boundary conditions prescribed on… More >

Tarun Kant¹, Sandeep S. Pendhari², Yogesh M. Desai³

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.2, pp. 135-162, 2007, DOI:10.3970/cmes.2007.017.135

Abstract A two-point boundary value problem (BVP) is formed in the present work governed by a set of first-order coupled ordinary differential equations (ODEs) in terms of displacements and the transverse stresses through the thickness of laminate (in domain -h/2 < z < h/2) by introducing partial discretization methodology only in the plan area of the three dimensional (3D) laminate. The primary dependent variables in the ODEs are those which occur naturally on a plane z=a constant. An effective numerical integration (NI) technique is utilized for tackling the two-point BVP in an efficient manner. Numerical studies on cross-ply and angle-ply composite… More >

C.-W. Chang¹, C.-S. Liu², J.-R. Chang^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 13-38, 2005, DOI:10.3970/cmes.2005.010.013

Abstract In this paper, the inverse heat conduction problem governed by sideways heat equation is investigated numerically. The problem is ill-posed because the solution, if it exists, does not depend continuously on the data. To begin with, this ill-posed problem is analyzed by considering the stability of the semi-discretization numerical schemes. Then the resulting ordinary differential equations at the discretized times are numerically integrated towards the spatial direction by the group preserving scheme, and the stable range of the index r = 1/2ν Δt is investigated. When the numerical results are compared with exact solutions, it is found that they are… More >