A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

定常的热传导-对流问题的Galerkin/Petrov最小二乘混合元方法

1.引言 热传导-对流问题是大气动力学中的一个重要的方程,这个方程组也称为强迫耗散的非线性系统方程组,其较Navier-Stokes方程多了一个未知函数温度场,且温度与速度和压力之间存在着复杂的非线性关系.从热动力学可知,任何运动都会产生热量即有温度,而且温度与速度和压力之间必定互相转化,因此对该非线性系统的研究更具有实际意义.[1]先对     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

Fully adaptive multiresolution finite volume schemes for conservation laws

The use of multiresolution decompositions in the context of finite volume schemes for conservation laws was first proposed by A. Harten for the purpose of accelerating the evaluation of numerical fluxes through an adaptive computation. In this approach the solution is still represented at each time step on the finest grid, resulting in an inherent limitation of the potential gain in memory space and computational time. The present paper is concerned with the development and the numerical analysis of fully adaptive multiresolution schemes, in which the solution is represented and computed in a dynamically evolved adaptive grid. A crucial problem is then the accurate computation of the flux without the full knowledge of fine grid cell averages. Several solutions to this problem are proposed, analyzed, and compared in terms of accuracy and complexity.

     

多周期公用工程系统运行的模型,优化方法与应用

针对多周期公用工程系统的运行优化问题，考虑了设备的启停费用的情况下。建立了混合整数非线性规划模型并证明了最优解的存在性。针对该运行优化问题本将其分解成若干子问题，然后利用改进的Hooke-Jeeves优化算法求解每个子问题。应用于具体实例，其数值结果与其它方法得到的相比。运行时间短，且更适合多周期公用工程问题的求解。     

Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.     

A Hybrid Approach to Scheduling with Earliness and Tardiness Costs

A hybrid technique using constraint programming and linear programming is applied to the problem of scheduling with earliness and tardiness costs. The linear model maintains a set of relaxed optimal start times which are used to guide the constraint programming search heuristic. In addition, the constraint programming problem model employs the strong constraint propagation techniques responsible for many of the advances in constraint programming for scheduling in the past few years. Empirical results validate our approach and show, in particular, that creating and solving a subproblem containing only the activities with direct impact on the cost function and then using this solution in the main search, significantly increases the number of problems that can be solved to optimality while significantly decreasing the search time.     

Banach空间中拟-似变分包含解的存在与逼近问题

研究了Banach空间中拟-似变分包含解的存在与逼近问题.给出了一种寻求解的新的迭代算法,建立了具混合误差的Ishikawa型迭代序列强收敛到解的充要条件.所得结果推广了一些相关的结果.     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.