A modified integral global optimization method and its asymptotic convergence

In this paper, we propose a new integral global optimization algorithm for finding the solution of continuous minimization problem, and prove the asymptotic convergence of this algorithm. In our modified method we use variable measure integral, importance sampling and main idea of the cross-entropy method to ensure its convergence and efficiency. Numerical results show that the new method is very efficient in some challenging continuous global optimization problems.     

POLLUTION, GOVERNMENT EXPENDITURE, TAXES AND STOCHASTIC GROWTH

This paper studies a stochastic endogenous growth model with pollution. It introduces government expenditure and exogenous pollution abatement technology to eliminate pollution and proves that under appropriate equilibrium conditions the main economic indexes （including economic growth rate, the optimal government expenditure rate） in the centrally planned economy and decentralized economy can be expressed by the parameters of the model uniquely. The optimal tax policy is analyzed ,and the optimal pollution is derived.     

AN EFFECTIVE CONTINUOUS ALGORITHM FOR APPROXIMATE SOLUTIONS OF LARGE SCALE MAX-CUT PROBLEMS

An effective continuous algorithm is proposed to find approximate solutions of NP-hardmax-cut problems.The algorithm relaxes the max-cut problem into a continuous nonlinearprogramming problem by replacing n discrete constraints in the original problem with onesingle continuous constraint.A feasible direction method is designed to solve the resultingnonlinear programming problem.The method employs only the gradient evaluations ofthe objective function,and no any matrix calculations and no line searches are required.This greatly reduces the calculation cost of the method,and is suitable for the solutionof large size max-cut problems.The convergence properties of the proposed method toKKT points of the nonlinear programming are analyzed.If the solution obtained by theproposed method is a global solution of the nonlinear programming problem,the solutionwill provide an upper bound on the max-cut value.Then an approximate solution to themax-cut problem is generated from the solution of the nonlinear programming and providesa lower bound on the max-cut value.Numerical experiments and comparisons on somemax-cut test problems(small and large size)show that the proposed algorithm is efficientto get the exact solutions for all small test problems and well satisfied solutions for mostof the large size test problems with less calculation costs.     

SEQUENTIAL CONVEX PROGRAMMING METHODS FOR SOLVING LARGE TOPOLOGY OPTIMIZATION PROBLEMS： IMPLEMENTATION AND COMPUTATIONAL RESULTS

In this paper, we describe a method to solve large-scale structural optimization problems by sequential convex programming （SCP）. A predictor-corrector interior point method is applied to solve the strictly convex subproblems. The SCP algorithm and the topology optimization approach are introduced. Especially, different strategies to solve certain linear systems of equations are analyzed. Numerical results are presented to show the efficiency of the proposed method for solving topology optimization problems and to compare different variants.     

Distance Between Two Vertices of Maximum Matching Graphs

The maximum matching graph of a graph has a vertex for each maximum matching and an edge for each pair of maximum matchings which differ by exactly one edge. In this paper, we obtain a lower bound of distance between two vertices of maximum matching graph, and give a necessary and sufficient condition that the bound can be reached.     

A GENERALIZED PROJECTION-SUCCESSIVE LINEAREQUATIONS ALGORITHM FOR NONLINEARLY EQUALITy AND INEQUALITY CONSTRAINED OPTIMIZATION AND ITS RATE OF CONVERGENCE

In this paper, a new superlinearly convergent algorithm is presented for optimization problems with general nonlineer equality and inequality Constraints, Comparing with other methods for these problems, the algorithm has two main advantages. First, it doesn‘t solve anyquadratic programming (QP), and its search directions are determined by the generalized projection technique and the solutions of two systems of linear equations. Second, the sequential points generated by the algoritbh satisfy all inequity constraints and its step-length is computed by the straight line search,The algorithm is proved to possesa global and auperlinear convergence.     

带常数项的最优加权组合预测

We propose a model based on the optimal weighted combinational forecasting with constant terms,give formulae of the weights and the average errors as well as a rela- tion of the model and the corresponding model without constant terms,and compare these models.Finally an example was given,which showed that the fitting precision has been enhanced.     

THE OPTIMAL ORDER ERROR ESTIMATES FOR FINITEELEMENT APPROXIMATIONS TO HYPERBOLIC PROBLEMS

In this paper, the linear finite element approximation to the positive and symmetric,linear hyperbolic systems is analyzed and an O(h^2) order error estimate is established under the conditions of strongly regular triangulation and the H^3-regularity for the exact solutions. The convergence analysis is based on some superclose estimates derived in this paper. Our method and result here are also applicable to general hyperbolic problems.Finally, we discuss the linearized shallow water system of equations.     

Extensions and Refinements ofAdamovic＇s Inequality

In this article, by means of tile theory of majorization, Adamovic‘s inequality isextended to the cases of the general elementary symmetric functions and its duals, and therefined and reversed forms are also given. As applications, some new inequalities for simplexare established.     

A Filter-Based Pattern Search Method for Unconstrained Optimization

We discuss a filter-based pattern search method for unconstrained optimization in this paper. For the purpose to broaden the search range we use both filter technique and frames, which are fragments of grids, to provide a new criterion of iterate acceptance. The convergence can be ensured under some conditions. The numerical result shows that this method is practical and efficient.     

RECENT PROGRESS ON SPHERICAL HARMONIC APPROXIMATION MADE BY BNU RESEARCH GROUP ——In memory of Professor Sun Yongsheng

As early as in 1990, Professor Sun Yongsheng, suggested his students at Beijing Normal University to consider research problems on the unit sphere. Under his guidance and encouragement his students started the research on spherical harmonic analysis and approximation. In this paper, we incompletely introduce the main achievements in this area obtained by our group and relative researchers during recent 5 years （2001-2005）. The main topics are： convergence of Cesaro summability, a.e. and strong summability of Fourier-Laplace series; smoothness and K-functionals; Kolmogorov and linear widths.     

Boundary-Value Problems of the Theory of Thin and Nonthin Orthotropic Shells with Account of Nonlinearly Elastic Properties and Low Shear Rigidity of Composite Materials

The basic geometric and physical relations and resolving equations of the theory of thin and nonthin orthotropic composite shells with account of nonlinear properties and low shear rigidity of their materials are presented. They are derived based on two theories, namely the theory of anisotropic shells employing the Timoshenko or Kirchhoff-Love hypothesis and the nonlinear theory of elasticity and plasticity of anisotropic media in combination with the Lagrange variational principle. The procedure and algorithm for the numerical solution of nonlinear (linear) problems are based on the method of successive approximations, the difference-variational method, and the Lagrange multiplier method. Calculations of the stress-strain state for a spherical shell with a circular opening loaded with internal pressure are presented. The effect of transverse shear strains and physical nonlinearity of the material on the distribution of maximum deflections and circumferential stresses in the shell, obtained according to two variants of the shell theories, is studied. A comparison of the results of the problem solution in linear and nonlinear statements with and without account of the shell shear strains is given. The numerical data obtained for thin and nonthin (medium thick) composite shells are analyzed.     

Active control of secondary resonances piezoelectric sandwich beams

Secondary resonances of piezoelectric/elastic/piezoelectric sandwich beams submitted to active control are studied in this paper. The proportional and derivative nonlinear potential feedback controls via piezoelectric sensor and actuator layers are used. The dynamics of the beam is modelled by a highly nonlinear ordinary-differential equation. The method of multiple scales is applied and approximate solutions are obtained for hard excitations. Analytical frequency and phase-amplitude relationships as well as the time response are explicitly given for various super- and subharmonic resonances. Static and dynamic stability criteria are elaborated and critical displacement and excitation amplitudes associated to the resulting unstable zones are analytically given. The feedback parameters effects on the subharmonic and superharmonic resonances and on their stability are investigated.     

Nonlinear and Geometric Programming – Current Status

Great strides have been made in nonlinear programming (NLP) in the last 5 years. In smooth NLP, there are now several reliable and efficient codes capable of solving large problems. Most of these implement GRG or SQP methods, and new software using interior point algorithms is under development. NLP software is now much easier to use, as it is interfaced with many modeling systems, including MSC/NASTRAN, and ANSYS for structural problems, GAMS and AMPL for general optimization, Matlab and Mathcad for general mathematical problems, and the widely used Microsoft Excel spreadsheet. For mixed integer problems, branch and bound and outer approximation codes are now available and are coupled to some of the above modeling systems, while search methods like Tabu Search and Genetic algorithms permit combinatorial, nonsmooth, and nonconvex problems to be attacked.     

A 3D modeling of static and forward smoldering combustion in a packed bed of materials

A three-dimensional (3D) model based on the first principles of mass, momentum and energy was developed that numerically simulates the processes of static and forward smoldering in a porous packed bed of plant materials. The packed bed contains cellulose material or tobacco (cigarette) wrapped in a porous paper and surrounded by an ambient air. Other major characteristics of the model are including the effects of buoyancy forces in the flow field, separate treatment of solid and gas in a thermally non-equilibrium environment, and use of multi-precursor kinetic models for the pyrolysis of staring material and oxidation of char. The changes in porosity due to pyrolysis and char oxidation and the effect of porosity on the bed permeability and gas diffusivity are included. The mass, momentum, energy, and species transport equations are solved in a discretized computational domain using a commercially available computational fluid dynamics (CFD) code. The simulation results show that the model reasonably reproduces the major features of a burning cigarette during smoldering and puffing and are in a good agreement with the existing experimental results for cigarettes. Results include the velocity profiles, gas and solid temperatures, coal shape, burn rates, profile and transport of gas and vapor species throughout the packed bed, dilution through the wrapper paper and ventilation in the filter section, and the mass fraction of some pyrolysis and oxidation products in the mainstream and sidestream flows.     

A system-theoretic model for cooperation, interaction and allocation

A system-theoretic approach to cooperation, interaction and allocation is presented that simplifies, unifies and extends the results on classical cooperative games and their generalizations. In particular, a general Weber theory of linear values is obtained and a new theory for local cooperation and general interaction indices is established. The model is dynamic and based on the notion of states of cooperation that change under actions of agents. Careful distinction between “local” states of cooperation and general “system” states leads to a notion of entropy for arbitrary non-negative and efficient allocations and thus to a new information-theoretic criterion for fairness of allocation mechanisms. Shapley allocations, for instance, are exhibited as arising from random walks with maximal entropy. For a large class of cooperation systems, a characterization of game symmetries in terms of λ-values is given. A concept for cores and Weber sets is proposed and it is shown that a Weber set of a game with selection structure always contains the core.     

拓扑空间上不动点和具有上下界的变分不等式

ωω根据广义凸空间上的KKM型定理和Fan-Browder型不动点定理, 得到了没有凸和线性结构且没有紧致框架的拓扑空间上的Φ -映射和弱Φ -映射的若干个新的不动点定理. 作为应用, 在非紧致的拓扑空间上讨论了具有上下界的变分不等式解的存在性问题.     

CFD modeling of hydrodynamic characteristics of a gas–liquid two-phase stirred tank

A three-dimensional CFD model was developed in this work to simulate hydrodynamic characteristics of a gas–liquid two-phase stirred tank with two six-bladed turbines and four baffles, coupling of the Multiple Size Group model to determine bubble size distribution. Important hydrodynamic parameters of the multi-phase system such as volume-averaged overall and time-averaged local gas holdups and axial liquid velocities along time and transversal courses were simulated and analyzed in detail, under varied operating conditions (inlet air flow rate and impeller rotation speed). Model predictions of local transient gas holdup and liquid velocity distributions on vertical and horizontal sections of the tank were also carried out. The overall flow patterns were discussed in detail to assess the mixing. Bubble size distributions were further predicted to reveal the unique properties of gas phase. Experimental measurements of overall gas holdups and local axial liquid velocities were used to validate the developed model.     

Financial Risk Modelling and Econometric Inference

In finance, the explicit modelling of uncertainty takes on a particularly important role. The values of financial derivatives increase in the return volatility of the underlying security. This notion requires a concept of volatility and hence uncertainty. In addition, the choice between modelling in discrete and continuous time is not arbitrary, since it corresponds to a distinction between incomplete and complete markets, respectively, and this distinction matters for asset pricing, financial risk modelling, and inference. Risk and volatility are closely connected, and implied volatility, volatility forecasting, volatility in term structure models, stochastic volatility, and portfolio analysis are considered and related to a more general interplay between cross-sectional and dynamic aspects in finance. Stocks, bonds, and options are considered and placed in the context of efficiency and separation in inference.     

Effects of permeability and viscosity in linear polymeric gels

We propose and analyze a mathematical model of the mechanics of gels, consisting of the laws of balance of mass and linear momentum of the polymer and liquid components of the gel. We consider a gel to be an immiscible and incompressible mixture of a nonlinearly elastic polymer and a fluid. The problems that we study are motivated by predictions of the life cycle of body‐implantable medical devices. Scaling arguments suggest neglecting inertia terms, and therefore, we consider the quasi‐static approximation to the dynamics. We focus on the linearized system about stress‐free states, uniform expansions, and compressions and derive sufficient conditions for the solvability of the time‐dependent problems. These turn out to be conditions that guarantee local stability of the equilibrium solutions. We also consider non‐stress free equilibria and states with residual stress and derive an energy law for the corresponding time‐dependent system. The conditions that guarantee stability of solutions provide a selection criteria of the material parameters of devices. The boundary conditions that we consider are of two types, displacement‐traction and permeability of the gel surface to the fluid. We address the cases of viscous and inviscid solvent, assume Newtonian dissipation for the polymer component, and establish existence of weak solutions for the different boundary permeability conditions and viscosity assumptions. We present two‐dimensional, finite element numerical simulations to study stress concentration on edges, this being the precursor to debonding of the gel from its substrate. Copyright © 2015 John Wiley & Sons, Ltd.