On the plane problem of the flow around a submerged beam

We consider the problem of the steady flow of an ideal heavy fluid around a submerged beam. The problem is obtained from the free-boundary problem of the flow past a submerged obstacle in the limit of bodies of vanishing thickness. We introduce a special Sobolev space formulation of the problem in term of a perturbed stream function and prove its unique solvability for every value of the unperturbed flow velocity, with the possible exception of a discrete set depending on the geometry of the domain. The asymptotic properties of the solutions are discussed.     

A multigrid procedure for Cartesian ghost‐cell methods

This paper proposes a multigrid technique for Cartesian grid flow solvers. A recently developed ghost body‐cell method for inviscid flows is combined with a nested‐level local refinement procedure, which employs multigrid to accelerate convergence to steady state. Different from standard multigrid applications for body‐fitted grids, a fictitious residual needs to be defined in the ghost cells to perform a correct residual collection and thus to avoid possible stalling of the multigrid procedure. The efficiency of the proposed local refinement multigrid Cartesian method is demonstrated for the case of the inviscid subsonic flow past a circular body. Copyright © 2008 John Wiley & Sons, Ltd.     

Subspace Barzilai-Borwein Gradient Method for Large-Scale Bound Constrained Optimization

An active set subspace Barzilai-Borwein gradient algorithm for large-scale bound constrained optimization is proposed. The active sets are estimated by an identification technique. The search direction consists of two parts: some of the components are simply defined; the other components are determined by the Barzilai-Borwein gradient method. In this work, a nonmonotone line search strategy that guarantees global convergence is used. Preliminary numerical results show that the proposed method is promising, and competitive with the well-known method SPG on a subset of bound constrained problems from CUTEr collection. This work was supported by the 973 project granted 2004CB719402 and the NSF project of China granted 10471036.     

Sign-changing solutions on a kind of fourth-order Neumann boundary value problem

In this paper, we consider the fourth-order Neumann boundary value problem u(4)(t)−2u^″(t)+u(t)=f(t,u(t)) for all t∈[0,1] and subject to u^′(0)=u^′(1)=u^?(0)=u^?(1)=0. Using the fixed point index and the critical group, we establish the existence theorem of solutions that guarantees the problem has at least one positive solution and two sign-changing solutions under certain conditions.     

Numerical analysis of turbulent flow separation in a rectangular duct with a sharp 180‐degree turn by algebraic Reynolds stress model

Turbulent flow in a rectangular duct with a sharp 180‐degree turn is difficult to predict numerically because the flow behavior is influenced by several types of forces, including centrifugal force, pressure‐driven force, and shear stress generated by anisotropic turbulence. In particular, this type of flow is characterized by a large‐scale separated flow, and it is difficult to predict the reattachment point of a separated flow. Numerical analysis has been performed for a turbulent flow in a rectangular duct with a sharp 180‐degree turn using the algebraic Reynolds stress model. A boundary‐fitted coordinate system is introduced as a method for coordinate transformation to set the boundary conditions next to complicated shapes. The calculated results are compared with the experimental data, as measured by a laser‐Doppler anemometer, in order to examine the validity of the proposed numerical method and turbulent model. In addition, the possibility of improving the wall function method in the separated flow region is examined by replacing the log‐law velocity profile for a smooth wall with that for a rough wall. The analysis results indicated that the proposed algebraic Reynolds stress model can be used to reasonably predict the turbulent flow in a rectangular duct with a sharp 180‐degree turn. In particular, the calculated reattachment point of a separated flow, which is difficult to predict in a turbulent flow, agrees well with the experimental results. In addition, the calculation results suggest that the wall function method using the log‐law velocity profile for a rough wall over a separated flow region has some potential for improving the prediction accuracy. Copyright © 2007 John Wiley & Sons, Ltd.     

CartGen: Robust,efficient and easy to implement uniform/octree/embedded boundary Cartesian grid generator

An efficient and easy to implement method to generate Cartesian grids is presented. The presented method generates various kinds of Cartesian grids such as uniform, octree and embedded boundary grids. It supports the variation of grid size along each spatial direction as well as anisotropic and non‐graded refinements. The efficiency and ease of implementation are the main benefits of the presented method in contrast to the alternative methods. Regarding octree grid generation, applying a simple and efficient data compression method permits to store all grid levels without considerable memory overhead. The presented method generates octree grids up to a 13‐level refinement (8192³ grids on the finest level) from a complicated geometry in a few minutes on the traditional desktop computers. The FORTRAN 90 implementation of the presented method is freely available under the terms of the GNU general public license. Copyright © 2007 John Wiley & Sons, Ltd.     

Comparison of Lagrangian bounds for one class of generalized assignment problems

Classical and modified Lagrangian bounds for the optimal value of optimization problems with a double decomposable structure are examined. For the class of generalized assignment problems, this property of constraints is used to design a Benders algorithm for solving the modified dual problem. Numerical results are presented that compare the quality of classical and modified bounds.     

Evolution of 3-D crystals from supersaturated vapor with modified Stefan condition: Galerkin method approach

We examine the evolution of crystals in three dimensions. We assume that the Wulff shape is a prism with a hexagonal base. We include the Gibbs-Thomson law on the crystal surface and the so-called Stefan condition. We show local in time existence of solutions assuming that the initial crystal has admissible shape.     

An effective ant colony optimization algorithm (ACO) for multi-objective resource allocation problem (MORAP)

The multi-objective resource allocation problem (MORAP) addresses the important issue which seeks to find the expected objectives by allocating the limited amount of resource to various activates. Resources may be manpower, assets, raw material or anything else in limited supply which can be used to accomplish the goals. The goals may be objectives (i.e., minimizing costs, or maximizing efficiency) usually driven by specific future needs. In this paper, in order to obtain a set of Pareto solution efficiently, we proposed a modified version of ant colony optimization (ACO), in this algorithm we try to increase the efficiency of algorithm by increasing the learning of ants. Effectiveness and efficiency of proposed algorithm was validated by comparing the result of ACO with hybrid genetic algorithm (hGA) which was applied to MORAP later.     

Energy Estimate for the Type-Ⅱ Superconducting Film

We study the minimizers of the Ginzburg-Landau model for variable thickness, superconducting, thin films with high k, placed in an applied magnetic field hex, when hex is of the order of the ＂first critical field＂, i.e. of the order of ｜lnε｜. We obtain the asymptotic estimates of minimal energy and describe the possible locations of the vortices.