Analysis of plane and axisymmetric flows of incompressible fluids with the stream tube method: Numerical simulation by trust-region optimization algorithm

New concepts for the study of incompressible plane or axisymmetric flows are analysed by the stream tube method. Flows without eddies and pure vortex flows are considered in a transformed domain where the mapped streamlines are rectilinear or circular. The transformation between the physical domain and the computational domain is an unknown of the problem. In order to solve the non-linear set of relevant equations, we present a new algorithm based on a trust region technique which is effective for non-convex optimization problems. Experimental results show that the new algorithm is more robust compared to the Newton-Raphson method.     

A geometric method in nonlinear programming

A differential geometric approach to the constrained function maximization problem is presented. The continuous analogue of the Newton-Raphson method due to Branin for solving a system of nonlinear equations is extended to the case where the system is under-determined. The method is combined with the continuous analogue of the gradient-projection method to obtain a constrained maximization method with enforced constraint restoration. Detailed analysis of the global behavior of both methods is provided. It is shown that the conjugate-gradient algorithm can take advantage of the sparse structure of the problem in the computation of a vector field, which constitutes the main computational task in the methods.This is part of a paper issued as Stanford University, Computer Science Department Report No. STAN-CS-77-643 (Ref. 45), which was presented at the Gatlinburg VII Conference, Asilomar, California, 1977. This work was supported in part by NSF Grant No. NAT BUR OF ECON RES/PO No. 4369 and by Department of Energy Contract No. EY-76-C-02-0016.The main part of this work was presented at the Japan-France Seminar on Functional Analysis and Numerical Analysis, Tokyo, Japan, 1976. The paper was prepared in part while the author was a visitor at the Department of Mathematics, North Carolina State University, Raleigh, North Carolina, 1976–77, and was completed while he was a visitor at the Computer Science Department, Stanford University, Stanford, California, 1977. He acknowledges the hospitality and stimulating environment provided by Professor G. H. Golub, Stanford University, and Professors N. J. Rose and C. D. Meyer, North Carolina State University.     

Initial guess by improved population-based intelligent algorithms for large inter-frame deformation measurement using digital image correlation

Digital image correlation (DIC) has received a widespread research and application in experimental mechanics. In DIC, the performance of subpixel registration algorithm (e.g., Newton-Raphson method, quasi-Newton method) relies heavily on the initial guess of deformation. In the case of small inter-frame deformation, the initial guess could be found by simple search scheme, the coarse-fine search for instance. While for large inter-frame deformation, it is difficult for simple search scheme to robustly estimate displacement parameters and deformation parameters simultaneously with low computational cost. In this paper, we proposed three improving strategies, i.e. Q-stage evolutionary strategy (T), parameter control strategy (C) and space expanding strategy (E), and then combined them into three population-based intelligent algorithms (PIAs), i.e. genetic algorithm (GA), differential evolution (DE) and particle swarm optimization (PSO), and finally derived eighteen different algorithms to calculate the initial guess for qN. The eighteen algorithms were compared in three sets of experiments including large rigid body translation, finite uniaxial strain and large rigid body rotation, and the results showed the effectiveness of proposed improving strategies. Among all compared algorithms, DE-TCE is the best which is robust, convenient and efficient for large inter-frame deformation measurement.     

Mathematical modeling on the effect of equivalence ratio in emission characteristics of compression ignition engine with hydrogen substitution

The investigation presented in this paper concerns on the computational simulation of emissions characteristics in compression ignition engine with hydrogen substitution. Combustion process has been modeled based on Equilibrium Constants Method (ECM) with MATLAB program to calculate the mole fractions of 18 combustion products when hydrogen is burnt along with diesel fuel at variable equivalence ratios. It can be observed that hydrogen substitution causes significant increase in NH₃, H₂, atom H emissions during rich combustion and OH, NO₂, HNO₃ emissions during lean combustion. As the equivalence ratio increases during rich combustion, mole fractions of HCN, CH₄, CO and atom C decreases with increment of hydrogen substitution. N₂, atom N and CO₂ emissions decrease whereas no significant changes in O₂, NO, O₃ and atom O emissions throughout all equivalence ratios as hydrogen is added to the combustion.     

On the numerical determination of optimal inputs

The design of optimal inputs for linear and nonlinear system identification involves the maximization of a quadratic performance index subject to an input energy constraint. In the classical approach, a Lagrange multiplier is introduced whose value is an unknown constant. In recent papers, the Lagrange multiplier has been determined by plotting a curve of the Lagrange multiplier as a function of the critical interval length or a curve of input energy versus the interval length. A new approach is presented in this paper in which the Lagrange multiplier is introduced as a state variable and evaluated simultaneously with the optimal input. Numerical results are given for both a linear and a nonlinear dynamic system.     

Penalized Regressions: The Bridge versus the Lasso

Abstract

Bridge regression, a special family of penalized regressions of a penalty function Σ|β_j|^γ with γ ≤ 1, considered. A general approach to solve for the bridge estimator is developed. A new algorithm for the lasso (γ = 1) is obtained by studying the structure of the bridge estimators. The shrinkage parameter γ and the tuning parameter λ are selected via generalized cross-validation (GCV). Comparison between the bridge model (γ ≤ 1) and several other shrinkage models, namely the ordinary least squares regression (λ = 0), the lasso (γ = 1) and ridge regression (γ = 2), is made through a simulation study. It is shown that the bridge regression performs well compared to the lasso and ridge regression. These methods are demonstrated through an analysis of a prostate cancer data. Some computational advantages and limitations are discussed.     

Photoionization cross sections of H2(+) and H2 with complex Gaussian-type basis functions optimized for the frequency-dependent polarizabilities

Within the framework of the complex basis function method, the photoionization cross sections of H(2)(+) and H(2) were calculated based on the variational principle for the frequency-dependent polarizabilities. In these calculations, complex orbital exponents of Gaussian-type basis functions for the final state continuum wavefunctions were fully optimized for each photon energy with the numerical Newton-Raphson method. In most cases, the use of only one or two complex Gaussian-type basis functions was enough to obtain excellent agreement with previous high precision calculations and available experimental results. However, there were a few cases, in which the use of complex basis functions having various angular momentum quantum numbers was crucial to obtain the accurate results. The behavior of the complex orbital exponents as a function of photon energy was discussed in relation to the scaling relation and the effective charge for photoelectron. The success of this method implies the effectiveness of the optimization of orbital exponents to reduce the number of basis functions and shows the possibility to calculate photoionization cross sections of general molecules using only Gaussian-type basis functions.     

缺失数据下含几何分布的对数线性模型的EM算法

本文研究缺失数据下对数线性模型参数的极大似然估计问题.通过Monte-Carlo EM算法去拟合所提出的模型.其中,在期望步中利用Metropolis-Hastings算法产生一个缺失数据的样本,在最大化步中利用Newton-Raphson迭代使似然函数最大化.最后,利用观测数据的Fisher信息得到参数极大似然估计的渐近方差和标准误差.     

The Tikhonov regularization method in elastoplasticity

The numerical simulation of the mechanical behavior of industrial materials is widely used for viability verification, improvement and optimization of designs. Elastoplastic models have been used to forecast the mechanical behavior of different materials. The numerical solution of most elastoplastic models comes across problems of ill-condition matrices. A complete representation of the nonlinear behavior of such structures involves the nonlinear equilibrium path of the body and handling of singular (limit) points and/or bifurcation points. Several techniques to solve numerical problems associated to these points have been disposed in the specialized literature. Two examples are the load-controlled Newton–Raphson method and displacement controlled techniques. However, most of these methods fail due to convergence problems (ill-conditioning) in the neighborhood of limit points, specially when the structure presents snap-through or snap-back equilibrium paths. This study presents the main ideas and formalities of the Tikhonov regularization method and shows how this method can be used in the analysis of dynamic elastoplasticity problems. The study presents a rigorous mathematical demonstration of existence and uniqueness of the solution of well-posed dynamic elastoplasticity problems. The numerical solution of dynamic elastoplasticity problems using Tikhonov regularization is presented in this paper. The Galerkin method is used in this formulation. Effectiveness of Tikhonov’s approach in the regularization of the solution of elastoplasticity problems is demonstrated by means of some simple numerical examples.     

Three-Dimensional Nonlinear Analysis of Drill String Structurein Annulus

ＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＮＯＮＬＩＮＥＡＲＡＮＡＬＹＳＩＳＯＦＤＲＩＬＬＳＴＲＩＮＧＳＴＲＵＣＴＵＲＥＩＮＡＮＮＵＩＵＳＬｉｕＹａｎ－ｑｉａｎｇ（刘延强）ＬｕＹｉｎｇ－ｍｉｎ（吕英民）（ＤｅｐａｒｔｍｅｎｔｏｆＥｎｇｉｎｅｅｒｉｎｇｍｅｃｈ...