移动火焰锋面(MFF)模型火焰锋面确定条件的进一步改进

考虑碳粒表面还原反应的移动火焰锋面模型(MFF-Ⅰ)在推导火焰锋面位置时,仅考虑了碳粒表面氧化反应的影响.本文综合考虑碳粒表面的氧化反应、还原反应对火焰锋面位置的影响,给出更严格的推导,即MFF-Ⅱ模型.在通常的燃烧工况下,改进模型MFF-Ⅱ与原始模型MFF-Ⅰ相比,预报结论基本一致.但对于较低氧浓度的工况,比如再燃、气化问题,改进模型MFF-Ⅱ同样可以给出可信的预报结果,而这是MFF-Ⅰ模型所不及的.     

Resolution of computational aeroacoustics problems on unstructured grids with a higher-order finite volume scheme

Computational fluid dynamics (CFD) has become increasingly used in the industry for the simulation of flows. Nevertheless, the complex configurations of real engineering problems make the application of very accurate methods that only work on structured grids difficult. From this point of view, the development of higher-order methods for unstructured grids is desirable. The finite volume method can be used with unstructured grids, but unfortunately it is difficult to achieve an order of accuracy higher than two, and the common approach is a simple extension of the one-dimensional case. The increase of the order of accuracy in finite volume methods on general unstructured grids has been limited due to the difficulty in the evaluation of field derivatives. This problem is overcome with the application of the Moving Least Squares (MLS) technique on a finite volume framework. In this work we present the application of this method (FV-MLS) to the solution of aeroacoustic problems.     

THE COMPACTION OF TIME-DEPENDENT VISCOUS GRANULAR MATERIALS CONSIDERING INERTIAL FORCES

In the present paper, compactions of time-dependent viscous granular materials are simulated step by step using the automatic adaptive mesh generation schemes. Inertial forces of the viscous incompressible aggregates are taken into account. The corresponding conservation equations, the weighted-integral formulations, and penalty finite element model are investigated. The fully discrete finite element equations for the simulation are derived. Polygonal particles of aggregates are simplified as mixed three-node and four-node elements. The automatic adaptive mesh generation schemes include contact detection algorithms, and mesh upgrade schemes. Solutions of the numerical simulation are in good agreement with some results from literatures. With minor modification, the proposed numerical model can be applied in several industries, including the pharmaceutical, ceramic, food, and household product manufacturing.     

A conservative unstructured scheme for rapidly varied flows

This article introduces a new semi‐implicit, staggered finite volume scheme on unstructured meshes for the modelling of rapidly varied shallow water flows. Rapidly varied flows occur in the inundation of dry land during flooding situations. They typically involve bores and hydraulic jumps after obstacles such as road banks. Near such sudden flow transitions, the grid resolution is often low compared with the gradients of the bathymetry. Locally the hydrostatic pressure assumption may become invalid. In these situations, it is crucial to apply the correct conservation properties to obtain accurate results. An important feature of this scheme is therefore its ability to conserve momentum locally or, by choice, preserve constant energy head along a streamline. This is achieved using a special interpolation method and control volumes for momentum. The efficiency of inundation calculations with locally very high velocities, and in the case of unstructured meshes locally very small grid distances, is severely hampered by the Courant condition. This article provides a solution in the form of a locally implicit time integration for the advective terms that allows for an explicit calculation in most of the domain, while maintaining unconditional stability by implicit calculations only where necessary. The complex geometry of flooded urban areas asks for the flexibility of unstructured meshes. The efficient calculation of the pressure gradient in this, and other semi‐implicit staggered schemes, requires, however, an orthogonality condition to be put on the grid. In this article a simple method is introduced to generate unstructured hybrid meshes that fulfil this requirement. Copyright © 2008 John Wiley & Sons, Ltd.     

Investigation of critical influential speed for moving mass problems on beams

A traveling mass due to its mass inertia has significant effects on the dynamic response of the structures. According to recent developments in structural materials and constructional technologies, the structures are likely to be affected by sudden changes of masses and substructure elements, in which the inertia effect of a moving mass is not negligible. The transverse inertia effects have been a topic of interest in bridge dynamics, design of railway tracks, guide way systems and other engineering applications such as modern high-speed precision machinery process. In this study an analytical–numerical method is presented which can be used to determine the dynamic response of beams carrying a moving mass, with various boundary conditions. It has been shown that the Coriolis acceleration, associated with the moving mass as it traverses along the vibrating beam shall be considered as well. Influences regarding the speed of the moving mass on the dynamic response of beams with various boundary conditions were also investigated. Results illustrated that the speed of a moving mass has direct influence on the entire structural dynamic response, depending on its boundary conditions. Critical influential speeds in the moving mass problems were introduced and obtained in numerical examples for various BC’s.     

拟矩形剖分上的样条空间的维数

矩形剖分~(记为$\Delta_{QR}$)~是指在矩形剖分~(记为$\Delta_{R}$)的基础上进行局部修改后得到的剖分,通常包括T-剖分~(记为$\Delta_{T}$)~和L-剖分~(记为$\Delta_{L}$).本文利用光滑余因子协调方法讨论了该剖分上的二元样条空间$S^\mu_k(\Delta_{QR})$的维数.在满足一定约束条件下, 得到了仅依赖于样条空间的次数,光滑度和剖分拓扑结构的显式维数公式.     

Evaluation of adaptive mesh refinement and coarsening for the computation of compressible flows on unstructured meshes

The compressible gas flows of interest to aerospace applications often involve situations where shock and expansion waves are present. Decreasing the characteristic dimension of the computational cells in the vicinity of shock waves improves the quality of the computed flows. This reduction in size may be accomplished by the use of mesh adaption procedures. In this paper an analysis is presented of an adaptive mesh scheme developed for an unstructured mesh finite volume upwind computer code. This scheme is tailored to refine or coarsen the computational mesh where gradients of the flow properties are respectively high or low. The refinement and coarsening procedures are applied to the classical gas dynamic problems of the stabilization of shock waves by solid bodies. In particular, situations where oblique shock waves interact with an expansion fan and where bow shocks arise around solid bodies are considered. The effectiveness of the scheme in reducing the computational time, while increasing the solution accuracy, is assessed. It is shown that the refinement procedure alone leads to a number of computational cells which is 20% larger than when alternate passes of refinement and coarsening are used. Accordingly, a reduction of computational time of the same order of magnitude is obtained. Copyright © 2005 John Wiley & Sons, Ltd.     

Sliding mesh algorithm for CFD analysis of helicopter rotor–fuselage aerodynamics

The study of rotor–fuselage interactional aerodynamics is central to the design and performance analysis of helicopters. However, regardless of its significance, rotor–fuselage aerodynamics has so far been addressed by very few authors. This is mainly due to the difficulties associated with both experimental and computational techniques when such complex configurations, rich in flow physics, are considered. In view of the above, the objective of this study is to develop computational tools suitable for rotor–fuselage engineering analysis based on computational fluid dynamics (CFD). To account for the relative motion between the fuselage and the rotor blades, the concept of sliding meshes is introduced. A sliding surface forms a boundary between a CFD mesh around the fuselage and a rotor‐fixed CFD mesh which rotates to account for the movement of the rotor. The sliding surface allows communication between meshes. Meshes adjacent to the sliding surface do not necessarily have matching nodes or even the same number of cell faces. This poses a problem of interpolation, which should not introduce numerical artefacts in the solution and should have minimal effects on the overall solution quality. As an additional objective, the employed sliding mesh algorithms should have small CPU overhead. The sliding mesh methods developed for this work are demonstrated for both simple and complex cases with emphasis placed on the presentation of the inner workings of the developed algorithms. Copyright © 2008 John Wiley & Sons, Ltd.     

Fourier series and integral equation method for the exterior Stokes problem

The exterior Stokes problem between two parallel planes that are separated by a prismatic cylinder is extended to the interior of the prism by requiring the continuity of the velocity across the lateral faces. The well‐posedness of the exterior–interior problem is proved in suitable weighted Sobolev spaces. The solution is represented by Fourier series in the z‐variable. The Fourier coefficients, solutions of auxiliary two‐dimensional exterior–interior problems, are analyzed by viewing them as boundary integral equations of potential theory and global regularity of the densities, is established in weighted Sobolev spaces of traces. A boundary element method, with suitably refined mesh size, is implemented for the numerical treatment of the Fourier coefficients. This provides optimal convergent semi‐ and fully‐discrete spectral methods of Fourier–Galerkin type. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Comments on ‘Simplex space–time meshes in finite element simulations’

Some comments are provided on the citations offered in a recent paper (M. Behr, Int. J. Numer. Meth. Fluids 2008; 57 :1421–1434) that describes space–time finite element computations of advection of ‘Gaussian hills’, including computations with mesh refinement in the time direction. Copyright © 2008 John Wiley & Sons, Ltd.