Optimal control of attitude for coupled-rigid-body spacecraft via Chebyshev-Gauss pseudospectral method

The attitude optimal control problem(OCP) of a two-rigid-body spacecraft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a discrete nonlinear programming(NLP) problem using the Chebyshev-Gauss pseudospectral method(CGPM). Solutions of the NLP problem can be obtained using the sequential quadratic programming(SQP) algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss(CG) points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform(FFT). To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and control variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.     

Meshfree weak–strong (MWS) form method and its application to incompressible flow problems

A meshfree weak–strong (MWS) form method has been proposed by the authors' group for linear solid mechanics problems based on a combined weak and strong form of governing equations. This paper formulates the MWS method for the incompressible Navier–Stokes equations that is non‐linear in nature. In this method, the meshfree collocation method based on strong form equations is applied to the interior nodes and the nodes on the essential boundaries; the local Petrov–Galerkin weak form is applied only to the nodes on the natural boundaries of the problem domain. The MWS method is then applied to simulate the steady problem of natural convection in an enclosed domain and the unsteady problem of viscous flow around a circular cylinder using both regular and irregular nodal distributions. The simulation results are validated by comparing with those of other numerical methods as well as experimental data. It is demonstrated that the MWS method has very good efficiency and accuracy for fluid flow problems. It works perfectly well for irregular nodes using only local quadrature cells for nodes on the natural boundary, which can be generated without any difficulty. Copyright © 2004 John Wiley & Sons, Ltd.     

微分求积法的特性及其改进

微分求积法已在科学和工程计算中得到了广泛应用。然而,有关时域微分求积法的数值稳定性、计算精度即阶数等基本特性,仍缺乏系统性的分析结论。依据微分求积法的基本原理,推导证明了微分求积法的权系数矩阵满足V-变换这一重要特性;利用微分求积法和隐式Runge-Kutta法的等值性,证明了时域微分求积法是A-稳定、s级s阶的数值方法。在此基础上,为进一步提高传统微分求积法的计算精度,利用待定系数法和Padé逼近,推导出了一类新的s级2s阶的微分求积法。数值计算对比结果验证了所提出的新微分求积法比传统的微分求积法具有更高的计算精度。     

极坐标系下Fourier-Legendre谱元方法与有限差分法数值扩散的比较

提出一种Fourier-Legendre谱元方法用于求解极坐标系下的Navier-Stokes方程,其中极点所在单元的径向采用Gauss-Radau积分点,避免了r=0处的1/r坐标奇异性。时间离散采用时间分裂法,引入数值同位素模型跟踪同位素的输运过程验证数值模拟的精度,分别利用谱元法和有限差分法的迎风差分格式求解匀速和加速坩埚旋转流动中的同位素方程。计算结果表明,有限差分法中的一阶迎风差分格式存在严重的数值假扩散,二阶迎风差分格式的数值结果较精确,增加节点可以有效地缓解数值扩散。然而,谱元法具有以较少节点得到高精度解的优势。     

Multiscale diffusion method for simulations of long-time defect evolution with application to dislocation climb

In many problems of interest to materials scientists and engineers, the evolution of crystalline extended defects (dislocations, cracks, grain boundaries, interfaces, voids, precipitates) is controlled by the flow of point defects (interstitial/substitutional atoms and/or vacancies) through the crystal into the extended defect. Precise modeling of this behavior requires fully atomistic methods in and around the extended defect, but the flow of point defects entering the defect region can be treated by coarse-grained methods. Here, a multiscale algorithm is presented to provide this coupling. Specifically, direct accelerated molecular dynamics (AMD) of extended defect evolution is coupled to a diffusing point defect concentration field that captures the long spatial and temporal scales of point defect motion in the presence of the internal stress fields generated by the evolving defect. The algorithm is applied to study vacancy absorption into an edge dislocation in aluminum where vacancy accumulation in the core leads to nucleation of a double-jog that then operates as a sink for additional vacancies; this corresponds to the initial stages of dislocation climb modeled with explicit atomistic resolution. The method is general and so can be applied to many other problems associated with nucleation, growth, and reaction due to accumulation of point defects in crystalline materials.     

Quadrature formulas of maximum algebraic accuracy for cauchy type integrals with complex exponents of the Jacobi weight function

The present paper presents a Gauss type quadrature formula for a Cauchy type integral whose density is the product of a Hölder function by the weight function (1 ? x) ^α (1 + x) ^β (Re α, Reβ > ?1) of orthogonal Jacobi polynomials. It is shown that at the roots of the function of the second kind corresponding to the Jacobi polynomial P _n ^(α,β) (x), the quadrature formula with n nodes gives the exact value of a Cauchy type integral for an arbitrary polynomial of order k ≤ 2n. This formula was tested when solving several contact and mixed problems of the theory of elasticity.     

In-plane vibration analysis of plates in curvilinear domains by a differential quadrature hierarchical finite element method

Free in-plane vibration analysis of plates is carried out by a differential quadrature hierarchical finite element method (DQHFEM). The NURBS (Non-Uniform Rational B-Splines) patches of geometries were first transformed into differential quadrature hierarchical (DQH) patches, and then the elastic field was discretized by the same DQH basis. The DQHFEM solved the compatibility problem caused by different parametrization of neighbouring patches of isogeometric analysis using NURBS. And mesh refinement in DQHFEM does not propagate from patch to patch. The DQHFEM matrices also have the embedding property as the hierarchical finite element method (HFEM). In-plane vibration analyses of plates of several planforms showed that the DQHFEM is similar as the fixed interface mode synthesis method that can analyse a structure using a few nodes on the boundary of substructure elements and only several clamped modes inside each substructure element, but the DQHFEM does not need modal analysis and is of high accuracy. The accuracy and convergence of the DQHFEM were validated through comparison with exact and approximate results in literatures and computed by the authors.     

压电-弹性层合梁自由振动分析的新方法

将状态空间法和微分求积法相结合,分析了压电-弹性层合梁的自由振动. 通过微分求积把状态方程在每一个节点处离散,进而求得解答. 选用不同的节点数目, 分析了方法的收敛性. 计算结果与相关文献的结果能较好地符合. 该方法 对于分析压电-弹性层合梁的工程振动问题非常方便.     

新型带虚点的径向基函数微分求积法及其在薄板弯曲中的应用

本文提出了新型带虚点的径向基函数微分求积法,并将其应用于模拟薄板弯曲问题。带虚点的径向基函数微分求积法是一种基于传统径向基函数微分求积法的新型无网格方法,传统方法只将中心点放在计算域内,而该方法扩展了中心点的区域,使其既位于计算域内又位于计算域外,在不增加计算量和存储量的基础上,显著提高计算精度。本文首次尝试将此方法应用于求解薄板弯曲问题,并与解析解和传统方法进行对比,验证了此方法的优越性     

Free surface seepage analysis based on the element-free method

Element-free method (EFM) for seepage analysis with a free surface presented in the paper is based on the moving least square method which needs only the information at nodes. It avoids troublesome modifications of the mesh as in the finite element method. Being irrelative to the nodes, the mesh for quadrature is fixed throughout the iterations in determining the free surface. And the nodes can be easily added, moved or deleted in the iterations. Steady seepage and transient seepage in a uniform earth dam were analyzed in the paper. The examples show that the proposed method gives satisfactory results.     

Quadrature method of moments for modeling multi-component spray vaporization

This article puts forward the quadrature method of moments (QMoM) for modeling droplet composition during the spray vaporization process. This method is implemented for solving the Continuous Thermodynamic Model (CTM) of multi-component droplet vaporization, an advantageous alternative to the classical Discrete Component Model (DCM) when the droplet is formed of a great number of components. The CTM approach consists in modeling the droplet’s composition using a probability density function (PDF). This method was first tried out for vaporizing droplets by Hallett, who assumed a Gamma-function for the PDF. However, Harstadt et al. underlined some problems in the case of vapor condensation on the droplet surface, since the Gamma-PDF model presumes the PDF’s mathematical form. The QMoM which does not require this hypothesis is studied in this article, according to Lage’s research dealing with QMoM application to phase equilibria. The numerical features of QMoM are investigated in detail, and then the method is implemented for the difficult test case of vapor condensation. The results are analyzed to illustrate the application of QMoM to multi-component droplet vaporization modeling and to provide a better understanding of the QMoM main advantages and limitations.     

微分求积单元法在结构工程中的应用

微分求积法（Differential Quadrature Method）是求鳃偏微分方程和积分-微分方程的一种数值方法，该法具有计算简便、精度较高和易于实现等优点。微分求积单元法（Differential Quadrature Element Method）是在微分求积法的基础上结合区域分割和集成规则而形成的一种新的数值计算方法，能通过自适应地选取微分求积网点数目正确模拟构件的刚度和荷载性质，其精度可通过细分单元或增加离散点数目加以提高。微分求积单元法是一种可供选择的、性能优越的数值计算方法。本文将详细论述这一数值方法的基本原理，并通过数值算例说明该方法的应用过程及其优越性，为这一方法在结构工程中的推广应用提供参考。     

Direct Quadrature Conditional Moment Closure for Modelling of Turbulent Combustion

We present a method of direct quadrature conditional moment closure (DQCMC) for the treatment of realistic turbulence-chemistry interaction in computational fluid dynamics (CFD) software. The method which is based on the direct quadrature method of moments (DQMOM) coupled with the conditional moment closure (CMC) equations is in simplified form and easily implementable in existing CMC formulation for CFD. The observed fluctuations of scalar dissipation around the conditional mean values are captured by the treatment of a set of mixing environments, each with its pre-defined weight. Unlike the early versions of the DQCMC method the resulting equations are similar to that of the first-order CMC, and the ??diffusion?? term is strictly positive and no correction factors are used. We present results for two mixing environments where the resulting matrices of the DQCMC can be inverted analytically. We have performed this analysis for a simple hydrogen flame using a multi species chemical scheme containing nine species. The effects of the fluctuations around the conditional means are captured accurately and the predicted results are in very good agreement with observed trends from direct numerical simulations. Furthermore, the differences between the first order CMC and DQCMC are discussed.     

Stationary Boundary Points for a Laplacian Growth Problem in Higher Dimensions

It is known that corners of interior angle less than π/2 in the boundary of a plane domain are initially stationary for Hele–Shaw flow arising from an arbitrary injection point inside the domain. This paper establishes the corresponding result for Laplacian growth of domains in higher dimensions. The problem is treated in terms of evolving families of quadrature domains for subharmonic functions.     

无网格径向点插值法在拓扑优化中的应用

论文基于无网格径向点插值法(RPIM)对连续体结构进行拓扑优化设计.以高斯点的相对密度为设计变量,以结构的柔度最小化为目标函数,利用带惩罚的各向同性固体材料模型(SIMP)和优化准则法,建立了设计变量的迭代格式.利用灵敏度过滤技术有效地消除了点状棋盘格现象.给出了相应的计算流程,并用Fortran程序语言实现其算法.算例证明,应用无网格径向点插值法能够有效地对连续体结构进行拓扑优化设计.     

极坐标系下的Legendre谱元方法求解Poisson-型方程

r=0处的坐标奇异性是求解极坐标下Poisson-型方程的关键。本文提出一种极坐标系下基于Galerkin变分的Legendre谱元方法用于求解圆形区域内的Poisson-型方程,物理区域的径向和周向划分若干单元,计算单元均采用Legendre多项式展开;圆心所在单元的径向使用LGR(Legendre Gauss Radau)积分点,其他单元径向使用LGL(Legendre Gauss Lobatto)积分点,从而避免了极点处1/r坐标奇异性,周向单元均采用LGL积分点。利用区域分解技术,可以避免节点在极点附近聚集;最后求解了多个Dirichlet或Neumann边界条件下的Poisson-型方程算例。数值结果表明,谱元方法具有很高的精度。     

Numerical simulations of viscous flows using a meshless method

This paper uses the element‐free Galerkin (EFG) method to simulate 2D, viscous, incompressible flows. The control equations are discretized with the standard Galerkin method in space and a fractional step finite element scheme in time. Regular background cells are used for the quadrature. Several classical fluid mechanics problems were analyzed including flow in a pipe, flow past a step and flow in a driven cavity. The flow field computed with the EFG method compared well with those calculated using the finite element method (FEM) and finite difference method. The simulations show that although EFG is more expensive computationally than FEM, it is capable of dealing with cases where the nodes are poorly distributed or even overlap with each other; hence, it may be used to resolve remeshing problems in direct numerical simulations. Flows around a cylinder for different Reynolds numbers are also simulated to study the flow patterns for various conditions and the drag and lift forces exerted by the fluid on the cylinder. These forces are calculated by integrating the pressure and shear forces over the cylinder surface. The results show how the drag and lift forces oscillate for high Reynolds numbers. The calculated Strouhal number agrees well with previous results. Copyright © 2008 John Wiley & Sons, Ltd.     

轴向均布载荷下压杆稳定问题的DQ解

叙述了微分求积法(differential quadrature method)的一般方法，研究用微分求积法求解在均布轴向载荷下细长杆的稳定问题．通过Newton-Raphson法求解非线性方程组，以及对问题进行线性假设后求解广义特征值方程，得到了精度很高的后屈曲挠度数值和临界载荷数值．与解析解和其他近似解相比，微分求积法具有较高的精度和简便性．     

Online monitoring of gas–solid two-phase flow using projected CG method in ECT image reconstruction

Electrical capacitance tomography (ECT) is a promising technique for multi-phase flow measurement due to its high speed, low cost and non-intrusive sensing. Image reconstruction for ECT is an inverse problem of finding the permittivity distribution of an object by measuring the electrical capacitances between sets of electrodes placed around its periphery. The conjugate gradient (CG) method is a popular image reconstruction method for ECT, in spite of its low convergence rate. In this paper, an advanced version of the CG method, the projected CG method, is used for image reconstruction of an ECT system. The solution space is projected into the Krylov subspace and the inverse problem is solved by the CG method in a low-dimensional specific subspace. Both static and dynamic experiments were carried out for gas–solid two-phase flows. The flow regimes are identified using the reconstructed images obtained with the projected CG method. The results obtained indicate that the projected CG method improves the quality of reconstructed images and dramatically reduces computation time, as compared to the traditional sensitivity, Landweber, and CG methods. Furthermore, the projected CG method was also used to estimate the important parameters of the pneumatic conveying process, such as the volume concentration, flow velocity and mass flow rate of the solid phase. Therefore, the projected CG method is considered suitable for online gas–solid two-phase flow measurement.