露天矿矿岩量网格化计算算法的改进研究

基于露天矿网格化地质模型矿岩量计算方法的思想,在分析了各种可能的采掘范围线与算量网格相互关系的基础上,确立了一种逻辑算法,并通过编程方法实现。该算法可以严格确定出单个网块内构成计算矿岩量多边形的点,从而提高了算量精度和可靠性。

A Parallel Algorithm for Adaptive Local Refinement of Tetrahedral Meshes Using Bisection

Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp：//lsec. cc. ac. cn/phg/）, a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.

Location of anemometer along Lanzhou-Xinjiang railway

Using structured mesh to discretize the calculation region, the wind velocity and pressure distribution in front of the wind barrier under different embankment heights are investigated based on the Detached Eddy Simulation(DES) with standard SpalartAllmaras(SA) model. The Reynolds number is 4.0×105 in this calculation. The region is three-dimensional. Since the wind barrier and trains are almost invariable cross-sections, only 25 m along the track is modeled. The height of embankment ranges from 1 m to 5 m and the wind barrier is 3 m high. The results show that the wind speed changes obviously before the wind barrier on the horizontal plane, which is 4.5 m high above the track. The speed of wind reduces gradually while approaching the wind barrier. It reaches the minimum value at a distance about 5 m before the wind barrier, and increases dramatically afterwards. The speed of wind at this location is linear with the speed of far field. The train aerodynamic coefficients decrease sharply with the increment of the embankment height. And they take up the monotonicity. Meanwhile, when the height increases from 3 m to 5 m, they just change slightly. It is concluded that the optimum anemometer location is nearly 5 m in front of the wind barrier.

Unstructured Grid Immersed Boundary Method for Numerical Simulation of Fluid Structure Interaction

This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.

Coupled flood and sediment transport modelling with adaptive mesh refinement

Coupled flood and sediment transport modelling in large-scale domains has for long been hindered by the high computational cost.Adaptive mesh refinement is one of the viable ways to solving this problem without degrading the accuracy.This goal can be accomplished through mesh adaptation,e.g.,mesh coarsening and refining based on the dynamic regime of the flow and sediment transport along with bed evolution.However,previous studies in this regard have been limited to cases either without involving sediment transport or featuring flow-sediment-bed decoupling and the assumption of sediment transport capacity,which are not generally justified.Here,a coupled hydrodynamic and non-capacity sediment transport model is developed on adaptive non-uniform rectangular mesh.The proposed model is validated against experimental tests and numerical results based on the fixed meshes.It is demonstrated that the proposed model can properly capture shock waves,resolve the wetting/drying transition and reproduce morphological evolution.Compared with models based on the fixed meshes,the proposed model features great advantage in computational efficiency and holds promise for wide applications.

ARAP++:an extension of the local/global approach to mesh parameterization

Mesh parameterization is one of the fundamental operations in computer graphics(CG) and computeraided design(CAD). In this paper, we propose a novel local/global parameterization approach, ARAP++, for singleand multi-boundary triangular meshes. It is an extension of the as-rigid-as-possible(ARAP) approach, which stitches together 1-ring patches instead of individual triangles. To optimize the spring energy, we introduce a linear iterative scheme which employs convex combination weights and a fitting Jacobian matrix corresponding to a prescribed family of transformations. Our algorithm is simple, efficient, and robust. The geometric properties(angle and area)of the original model can also be preserved by appropriately prescribing the singular values of the fitting matrix. To reduce the area and stretch distortions for high-curvature models, a stretch operator is introduced. Numerical results demonstrate that ARAP++ outperforms several state-of-the-art methods in terms of controlling the distortions of angle, area, and stretch. Furthermore, it achieves a better visualization performance for several applications, such as texture mapping and surface remeshing.

Numerical method of static aeroelastic correction and jig-shape design for large airliners

In this paper,a coupled CFD-CSD method based on N-S equations is described for static aeroelastic correction and jig-shape design of large airliners.The wing structural flexibility matrix is analyzed by a finite element method with a double-beam model.The viscous multi-block structured grid is used in aerodynamic calculations.Flexibility matrix interpolation is fulfilled by use of a surface spline method.The load distributions on wing surface are evaluated by solving N-S equations with a parallel algorithm.A flexibility approach is employed to calculate the structural deformations.By successive iterations between steady aerodynamic forces and structural deformations,a coupled CFD-CSD method is achieved for the static aeroelastic correction and jig-shape design of a large airliner.The present method is applied to the static aeroelastic analysis and jig-shape design for a typical large airliner with engine nacelle and winglet.The numerical results indicate that calculations of static aeroelastic correction should employ tightly coupled CFD-CSD iterations,and that on a given cruise shape only one round of iterative design is needed to obtain the jig-shape meeting design requirements.