Computational experience with general equilibrium problems

We report on computational experience with an implementation of three algorithms for the general economic equilibrium problem. As a result we get that the projection algorithm for variational inequalities increases the size of solvable models by a factor of 5–10 in comparison with the classical homotopy method. As a third approach we implemented a simulated annealing heuristic which might be suitable to estimate equilibria for very large models.Supported by the German Research Association (Deutsche Forschungsgemeinschaft, SFB 303).     

A 3D-BPP approach for optimising stowage plans and terminal productivity

This paper addresses the problem of determining stowage plans for containers in a ship, that is the so-called master bay plan problem (MBPP).     

多级轴流压气机全工况特性计算

本文使用三维粘性流动计算软件Fine／Numeca,对某十五级轴流压气机进行了内流流场和全工况特性的数值计算尝试。分析了该压气机在设计工况和非设计工况的性能,同时把整机计算结果和前七级叶片的计算结果进行了比较。计算结果表明,当计算的级数较少时,目前的软件和硬件平台可以比较合理地预测压气机的全工况特性;而当计算的级数较多时,准确的数值模拟仍需要更为准确的多级模型和数值方法。     

Meshless simulation of equilibrium swelling/deswelling of pH‐sensitive hydrogels

Hydrogels have been widely used in microelectromechanical systems (MEMS) and Bio‐MEMS devices. In this article, the equilibrium swelling/deswelling of the pH‐stimulus cylindrical hydrogel in the microchannel is studied and simulated by the meshless method. The multi‐field coupling model, called multi‐effect‐coupling pH‐stimulus (MECpH) model, is presented and used to describe the chemical field, electric field, and the mechanical field involved in the problem. The partial differential equations (PDEs) describing these three fields are either nonlinear or coupled together. This multi‐field coupling and high nonlinear characteristics produce difficulties for the conventional numerical methods (e.g., the finite element method or the finite difference method), so an alternative—meshless method is developed to discretize the PDEs, and the efficient iteration technique is adopted to solve the nonlinear problem. The computational results for the swelling/deswelling diameter of the hydrogel under the different pH values are firstly compared with experimental results, and they have a good agreement. The influences of other parameters on the mechanical properties of the hydrogel are also investigated in detail. It is shown that the multi‐field coupling model and the developed meshless method are efficient, stable, and accurate for simulation of the properties of the stimuli‐sensitive hydrogel. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 326–337, 2006     

The nanogranular nature of C-S-H

Despite its ubiquitous presence as binding phase in all cementitious materials, the mechanical behavior of calcium-silicate-hydrates (C-S-H) is still an enigma that has deceived many decoding attempts from experimental and theoretical sides. In this paper, we propose and validate a new technique and experimental protocol to rationally assess the nanomechanical behavior of C-S-H based on a statistical analysis of hundreds of nanoindentation tests. By means of this grid indentation technique we identify in situ two structurally distinct but compositionally similar C-S-H phases heretofore hypothesized to exist as low density (LD) C-S-H and high density (HD) C-S-H, or outer and inner products. The main finding of this paper is that both phases exhibit a unique nanogranular behavior which is driven by particle-to-particle contact forces rather than by mineral properties. We argue that this nanomechanical blueprint of material invariant behavior of C-S-H is a consequence of the hydration reactions during which precipitating C-S-H nanoparticles percolate generating contact surfaces. As hydration proceeds, these nanoparticles pack closer to center on-average around two characteristic limit packing densities, the random packing limit (η=64%) and the ordered face-centered cubic (fcc) or hexagonal close-packed (hcp) packing limit (η=74%), forming a characteristic LD C-S-H and HD C-S-H phase.     

Constant Q-curvature metrics in arbitrary dimension

Working in a given conformal class, we prove existence of constant Q-curvature metrics on compact manifolds of arbitrary dimension under generic assumptions. The problem is equivalent to solving a nth-order non-linear elliptic differential (or integral) equation with variational structure, where n is the dimension of the manifold. Since the corresponding Euler functional is in general unbounded from above and below, we use critical point theory, jointly with a compactness result for the above equation.     

On the evolution of two-dimensional patterns in an inviscid liquid sheet

We perform an analysis of the pattern formation for a moving sheet of inviscid fluid. The sheet, which is assumed to have an infinite horizontal extent, moves at some prescribed velocity into a passive surrounding gas. The sheet’s thickness is assumed much smaller than the horizontal scale of the fluid motion. By considering a system that is symmetric with respect to the horizontal planes, long scale asymptotics are used to reduce the full governing equations in three dimensions to a set of three coupled nonlinear partial differential equations for the horizontal components of the velocity field and the height of the interface profile. The interfacial conditions consisting of the kinematic and normal stress balance are incorporated into these evolution equations. Investigations are carried out as function of the sole dimensionless parameter, namely the Weber number. A small amplitude stability analysis around the planar gas–liquid interface reveals that wave patterns in the form of traveling plane waves occur subcritically, and are therefore unstable. The reduced evolution equations are solved numerically for fixed values of the Weber number. Since the reduced system of equations is homogeneous, the wave motion is generated by initial conditions. Five initial conditions have been imposed: one-dimensional rolls, two-dimensional squares, two-dimensional hexagons, two-dimensional ridges, and smooth peaks. The ensuing evolution of the liquid sheet’s shape and corresponding flow fields are described by illustrations of the changes in the sheet’s morphology with time.     

两种群竞争模型的奇摄动群

本文应用多重尺度法构造出非线性微分方程组的解的渐近展开式。并用微分不等式的技巧，证明原问题的解的存在性，且给出解的一致有效渐近估计．     

An evolutionary method for optimization of plate buckling resistance

Optimization of plate buckling resistance is very complicated, because the in-plane stress resultants in the prebuckled state of a plate are functions of thickness distribution. This paper discusses the problem of finding the optimum thickness distribution of isotropic plate structures, with a given volume and layout, that would maximise the buckling load. A simple numerical method using the finite-element analysis is presented to obtain the optimum thickness distribution. Optimum designs of compression-loaded rectangular plates with different boundary conditions and plate aspect ratios are obtained by using the proposed method. Optimum designs from earlier studies and the methods of buckling analysis used to attain these results are discussed and compared with the designs from the proposed method. This paper also examines the reliability of the optimality criterion generally used for plate buckling optimization, which is based on the uniform strain energy density.     

An analysis of reduced Hessian methods for constrained optimization

We study the convergence properties of reduced Hessian successive quadratic programming for equality constrained optimization. The method uses a backtracking line search, and updates an approximation to the reduced Hessian of the Lagrangian by means of the BFGS formula. Two merit functions are considered for the line search: the ₁ function and the Fletcher exact penalty function. We give conditions under which local and superlinear convergence is obtained, and also prove a global convergence result. The analysis allows the initial reduced Hessian approximation to be any positive definite matrix, and does not assume that the iterates converge, or that the matrices are bounded. The effects of a second order correction step, a watchdog procedure and of the choice of null space basis are considered. This work can be seen as an extension to reduced Hessian methods of the well known results of Powell (1976) for unconstrained optimization.This author was supported, in part, by National Science Foundation grant CCR-8702403, Air Force Office of Scientific Research grant AFOSR-85-0251, and Army Research Office contract DAAL03-88-K-0086.This author was supported by the Applied Mathematical Sciences subprogram of the Office of Energy Research, U.S. Department of Energy, under contracts W-31-109-Eng-38 and DE FG02-87ER25047, and by National Science Foundation Grant No. DCR-86-02071.