A vortex axis and vortex core border grid adaptation algorithm

In the design process of hydrodynamical and aerodynamical technical applications, the numerical simulation of massively separated vortical flow is crucial for predicting, for example, lift or drag. To obtain reliable numerical results, it is mandatory to accurately predict the physical behavior of vortices. Thus, the dominant vortical flow structures have to be resolved in detail, which requires a local grid refinement and certain adaptation techniques. In this paper, a vortex flow structure adaptation algorithm is presented, which is particularly designed for local grid refinement at vortex axes positions and associated vortex core border locations. To this end, a fast and efficient vortex axis detection scheme is introduced and the algorithm for the vortex core border determination is explained. As the interaction between vortices makes the assignment of grid points to a certain vortex axis difficult, a helicity‐based vortex distinction approach in combination with a geometrical rotational sensor is developed. After describing the combined different techniques in detail, the vortex feature adaptation algorithm is applied to analytical and more realistic examples, which show that the described grid adaptation algorithm is able to enhance the grid cell resolution locally such that all significant vortical flow phenomena are resolved. Copyright © 2008 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.     

Adaptive Runge–Kutta discontinuous Galerkin method for complex geometry problems on Cartesian grid

A Cartesian grid method using immersed boundary technique to simulate the impact of body in fluid has become an important research topic in computational fluid dynamics because of its simplification, automation of grid generation, and accuracy of results. In the frame of Cartesian grid, one often uses finite volume method with second order accuracy or finite difference method. In this paper, an h‐adaptive Runge–Kutta discontinuous Galerkin (RKDG) method on Cartesian grid with ghost cell immersed boundary method for arbitrarily complex geometries is developed. A ghost cell immersed boundary treatment with the modification of normal velocity is presented. The method is validated versus well documented test problems involving both steady and unsteady compressible flows through complex bodies over a wide range of Mach numbers. The numerical results show that the present boundary treatment to some extent reduces the error of entropy and demonstrate the efficiency, robustness, and versatility of the proposed approach. Copyright © 2013 John Wiley & Sons, Ltd.     

一种特殊设计约束面的有限元网格优化问题

网格设计质量的优劣是有限元分析的关键，对一些特殊的设计约束面，网格的优化设计是建模时需要首先研究解决的问题，本文对两种材料在几何尺寸差异很大情况下的有限元网格进行了优化设计并对其进行了讨论。     

The structure of water-air bubble grid turbulence in a square duct

Local measurements of void fraction and continuous phase velocity field in water-air bubble, grid turbulence were conducted in a channel of vertical, square test section. The measured statistics indicate that, due mainly to the interaction of mean shear with the dispersed phase, the turbulence structure of the flow is modified. The observed change is characterized by a strong spatial dependence of void fraction and liquid flow properties, and the emergence of two spatial regions controlled by different physical processes. Intensity measurements indicate significant departure from isotropy in the flow. Two distinct regimes corresponding to low and high values of void fraction have been also identified. The autocorrelation and spectra measurements indicate that for low void fraction the scales of turbulence decrease while for higher values of void fraction increase again and inverse cascade is observed.     

A Poisson–Boltzmann solver on irregular domains with Neumann or Robin boundary conditions on non-graded adaptive grid

We introduce a second-order solver for the Poisson–Boltzmann equation in arbitrary geometry in two and three spatial dimensions. The method differs from existing methods solving the Poisson–Boltzmann equation in the two following ways: first, non-graded Quadtree (in two spatial dimensions) and Octree (in three spatial dimensions) grid structures are used; Second, Neumann or Robin boundary conditions are enforced at the irregular domain’s boundary. The irregular domain is described implicitly and the grid needs not to conform to the domain’s boundary, which makes grid generation straightforward and robust. The linear system is symmetric, positive definite in the case where the grid is uniform, nonsymmetric otherwise. In this case, the resulting matrix is an M-matrix, thus the linear system is invertible. Convergence examples are given in both two and three spatial dimensions and demonstrate that the solution is second-order accurate and that Quadtree/Octree grid structures save a significant amount of computational power at no sacrifice in accuracy.     

Numerical modeling of transient two-dimensional viscoelastic waves

This paper deals with the numerical modeling of transient mechanical waves in linear viscoelastic solids. Dissipation mechanisms are described using the generalized Zener model. No time convolutions are required thanks to the introduction of memory variables that satisfy local-in-time differential equations. By appropriately choosing the relaxation parameters, it is possible to accurately describe a large range of materials, such as solids with constant quality factors. The evolution equations satisfied by the velocity, the stress, and the memory variables are written in the form of a first-order system of PDEs with a source term. This system is solved by splitting it into two parts: the propagative part is discretized explicitly, using a fourth-order ADER scheme on a Cartesian grid, and the diffusive part is then solved exactly. Jump conditions along the interfaces are discretized by applying an immersed interface method. Numerical experiments of wave propagation in viscoelastic and fluid media show the efficiency of this numerical modeling for dealing with challenging problems, such as multiple scattering configurations.     

三维目标与粗糙面复合散射的广义稀疏矩阵平面迭代及规范网格算法

提出了快速计算二维导体粗糙面与面上金属目标复合散射的广义稀疏矩阵平面迭代及规范网格法(G-SMFSIA/CAG).推导了二维导体粗糙面与面上目标相互作用的耦合积分方程,用稀疏矩阵平面迭代及规范网格法(SMFSIA/CAG)求解粗糙面部分的表面积分方程,而用基于RWG基函数的矩量法(MOM)计算目标部分的表面积分方程,并通过更新方程的激励项迭代求解目标与粗糙面的相互耦合作用.结合Monte-Carlo方法产生具有PM(Pierson-Moskowitz)海浪谱的随机海洋粗糙面,数值分析了海面上不同形状导体目 关键词： 复合散射 广义稀疏矩阵平面迭代及规范网格法 随机海洋粗糙面 双站散射系数