A fast multi-resolution approach to tomographic PIV

Tomographic particle image velocimetry (Tomo-PIV) is a recently developed three-component, three-dimensional anemometric non-intrusive measurement technique, based on an optical tomographic reconstruction applied to simultaneously recorded images of the distribution of light intensity scattered by seeding particles immersed into the flow. Nowadays, the reconstruction process is carried out mainly by iterative algebraic reconstruction techniques, well suited to handle the problem of limited number of views, but computationally intensive and memory demanding. The adoption of the multiplicative algebraic reconstruction technique (MART) has become more and more accepted. In the present work, a novel multi-resolution approach is proposed, relying on the adoption of a coarser grid in the first step of the reconstruction to obtain a fast estimation of a reliable and accurate first guess. A performance assessment, carried out on three-dimensional computer-generated distributions of particles, shows a substantial acceleration of the reconstruction process for all the tested seeding densities with respect to the standard method based on 5 MART iterations; a relevant reduction in the memory storage is also achieved. Furthermore, a slight accuracy improvement is noticed. A modified version, improved by a multiplicative line of sight estimation of the first guess on the compressed configuration, is also tested, exhibiting a further remarkable decrease in both memory storage and computational effort, mostly at the lowest tested seeding densities, while retaining the same performances in terms of accuracy.     

基于TDLAS的层析成像技术TDLAT

可调谐二极管吸收光谱技术（tunable diode laser absorption spectroscopy,TDLAS）已成为超声速燃烧的重要测量手段之一. 为提高其空间分辨能力,需发展基于TDLAS,结合层析成像术的二维断层成像术（tunablediode laser absorption tomography, TDLAT）. 设计了一套基于6 平行光束-旋转测量的新型TDLAT 系统,吸收波长为7 185.6 cm-1 和7 444.3 cm-1 双线,采用分时-直接吸收探测策略. 重建中,使用代数重建算法,先分别反演计算两吸收线的吸收率和吸收比分布,再获得温度和浓度分布. 利用该系统,在CH₄/Air 预混平面燃烧炉上开展初步验证试验. 结果表明,TDLAT 系统可以反演出温度和浓度分布特征,反演的温度分布结果与热电偶测量值吻合较好. 进一步改进该系统,可用于超燃直连台中,测量燃烧室出口气流的温度和组分浓度分布.      

Three-dimensional tomographic reconstruction of unstable ejection phenomena of droplets for electrohydrodynamic jet

The electrohydrodynamic (EHD) jet is becoming increasingly popular within industrial printing areas based on phenomena induced by electrical potentials. Regardless of the physical observations of unstable ejection phenomena in regions possessing high electric potential, quantitative visualization is still necessary; no report exists exemplifying quantitative visualization. Thus, the size, shape and position of EHD droplets were reconstructed in this study using developed three-dimensional tomography methods. Two computer-synthesized phantoms for the liquid meniscus containing small satellite droplets were generated according to actual images captured by two high-speed cameras. These droplets were made in order to numerically reconstruct droplet behavior. Four three-dimensional tomography methods, such as the algebraic reconstruction technique (ART), the adaptive algebraic reconstruction technique (AART), the simultaneous iterative reconstruction technique (SIRT) and the multiplicative algebraic reconstruction technique (MART), were developed to accurately mimic droplet movement using multiple image views. Four basis functions including the cubic B-spline, cosine, o-Moms and Keys basis functions were adopted in order to improve the performance of the tomographic reconstructions. After completing a comparison of the four tomography results, the MART method in association with the cubic cosine basis function was selected as the means to significantly improve reconstruction accuracy. Additionally, it was the applied method for the reconstruction of the droplet behavior from experimental projections by two cameras.     

基于TDLAS的层析成像技术TDLAT简

可调谐二极管吸收光谱技术（tunable diode laser absorption spectroscopy，TDLAS）已成为超声速燃烧的重要测量手段之一. 为提高其空间分辨能力，需发展基于TDLAS，结合层析成像术的二维断层成像术（tunablediode laser absorption tomography, TDLAT）. 设计了一套基于6 平行光束-旋转测量的新型TDLAT 系统，吸收波长为7 185.6 cm^-1 和7 444.3 cm^-1 双线，采用分时-直接吸收探测策略. 重建中，使用代数重建算法，先分别反演计算两吸收线的吸收率和吸收比分布，再获得温度和浓度分布. 利用该系统，在CH₄/Air 预混平面燃烧炉上开展初步验证试验. 结果表明，TDLAT 系统可以反演出温度和浓度分布特征，反演的温度分布结果与热电偶测量值吻合较好. 进一步改进该系统，可用于超燃直连台中，测量燃烧室出口气流的温度和组分浓度分布.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法,用以求解强吸气作用下旋转圆筒表面边界层流动.首先,理论分析了同心圆筒间的旋转流体运动,外筒静止、内筒旋转且为多孔吸气条件.强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动,基于这一事实得到了周向速度分布的解析表达式.其次,通过引入新参数扩展Cebeci-Smith代数湍流模型,使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素.针对这些因素的综合影响,采用解析修正和经验参数对模型进行调整.同时,基于Reynolds应力湍流模型的仿真结果,校准代数湍流模型中的经验参数.最后,给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法,该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况.计算了不同旋转速度和吸气强度组合工况下的边界层流动,其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近.并且表明,当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时,该方法也再现了相同初始条件下的层流边界层.     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

Predicting Buoyant Shear Flows Using Anisotropic Dissipation Rate Models

This paper examines the modeling of two-dimensional homogeneous stratified turbulent shear flows using the Reynolds-stress and Reynolds-heat-flux equations. Several closure models have been investigated; the emphasis is placed on assessing the effect of modeling the dissipation rate tensor in the Reynolds-stress equation. Three different approaches are considered; one is an isotropic approach while the other two are anisotropic approaches. The isotropic approach is based on Kolmogorov's hypothesis and a dissipation rate equation modified to account for vortex stretching. One of the anisotropic approaches is based on an algebraic representation of the dissipation rate tensor, while another relies on solving a modeled transport equation for this tensor. In addition, within the former anisotropic approach, two different algebraic representations are examined; one is a function of the Reynolds-stress anisotropy tensor, and the other is a function of the mean velocity gradients. The performance of these closure models is evaluated against experimental and direct numerical simulation data of pure shear flows, pure buoyant flows and buoyant shear flows. Calculations have been carried out over a range of Richardson numbers (Ri) and two different Prandtl numbers (Pr); thus the effect of Pr on the development of counter-gradient heat flux in a stratified shear flow can be assessed. At low Ri, the isotropic model performs well in the predictions of stratified shear flows; however, its performance deteriorates as Ri increases. At high Ri, the transport equation model for the dissipation rate tensor gives the best result. Furthermore, the results also lend credence to the algebraic dissipation rate model based on the Reynolds stress anisotropy tensor. Finally, it is found that Pr has an effect on the development of counter-gradient heat flux. The calculations show that, under the action of shear, counter-gradient heat flux does not occur even at Ri = 1 in an air flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

An extended algebraic reconstruction technique (ART) for density-gradient projections: laser speckle photographic tomography

An extended algebraic reconstruction technique (ART) is presented for tomographic image reconstruction from the density-gradient projections, such as laser speckle photography. The essence of the extended ART is that the density-gradient projection data of speckle photography (Eq. (1)) are first numerically integrated to the algebraic representation of interferometric fringe number data (Eq. (12)), which ART can readily reconstruct into the cross-sectional field. The extended ART is numerically examined by using two computer synthesized phantom fields, and experimentally by using asymmetric single and double helium jets in air. The experimentally reconstructed images were also compared with the direct measurements of helium concentration using an oxygen analyzing probe. The extended ART method shows an improved accuracy and is proposed to use to tomographically reconstruct the density-gradient projections over the previous Fourier convolution (FC) method (Liu et al. 1989). Received: 26 June 1998/Accepted: 18 March 1999     

Influence of contact resistance and natural convection effects on non-dimensional effective thermal conductivity estimation of two-phase materials

In this article, the development of geometry dependent resistance model by considering contact resistance and natural convection effects are used to estimate the effective thermal conductivity of two-phase materials based on the unit cell approach. The algebraic equations have been derived based on isotherm approach for 2-D and 3-D spatially periodic medium. Comparison study has been made between developed models and experimental data. The result agrees well with experimental values.     

Problème spectral pour la propagation conique des ondes élastiques dans un réseau de galeriesA spectral problem for conically propagating elastic waves through an array of cylindrical channels

This Note is devoted to the analysis of elastic waves conically propagating through a doubly periodic array of cylindrical channels. A new method, based on a multiple scattering approach, has been proposed to reduce the problem to an algebraic system of the Rayleigh type. We obtain an eigenvalue problem formulation that enables us to construct the high-order dispersion curves and to study phononic band gap structures in oblique propagation. We note an effect of singular perturbation associated with a small angle of conical propagation. To cite this article: S. Guenneau et al., C. R. Mecanique 330 (2002) 491–497.     

An efficient and accurate algebraic interface capturing method for unstructured grids in 2 and 3 dimensions: The THINC method with quadratic surface representation

We present in this paper an efficient and accurate volume of fluid (VOF) type scheme to compute moving interfaces on unstructured grids with arbitrary quadrilateral mesh elements in 2D and hexahedral elements in 3D. Being an extension of the multi‐dimensional tangent of hyperbola interface capturing (THINC) reconstruction proposed by the authors in Cartesian grid, an algebraic VOF scheme is devised for arbitrary quadrilateral and hexahedral elements. The interface is cell‐wisely approximated by a quadratic surface, which substantially improves the numerical accuracy. The same as the other THINC type schemes, the present method does not require the explicit geometric representation of the interface when computing numerical fluxes and thus is very computationally efficient and straightforward in implementation. The proposed scheme has been verified by benchmark tests, which reveal that this scheme is able to produce high‐quality numerical solutions of moving interfaces in unstructured grids and thus a practical method for interfacial multi‐phase flow simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

Portable implementation model for CFD simulations. Application to hybrid CPU/GPU supercomputers

Nowadays, high performance computing (HPC) systems experience a disruptive moment with a variety of novel architectures and frameworks, without any clarity of which one is going to prevail. In this context, the portability of codes across different architectures is of major importance. This paper presents a portable implementation model based on an algebraic operational approach for direct numerical simulation (DNS) and large eddy simulation (LES) of incompressible turbulent flows using unstructured hybrid meshes. The strategy proposed consists in representing the whole time-integration algorithm using only three basic algebraic operations: sparse matrix–vector product, a linear combination of vectors and dot product. The main idea is based on decomposing the nonlinear operators into a concatenation of two SpMV operations. This provides high modularity and portability. An exhaustive analysis of the proposed implementation for hybrid CPU/GPU supercomputers has been conducted with tests using up to 128 GPUs. The main objective consists in understanding the challenges of implementing CFD codes on new architectures.     

Tomographic particle tracking velocimetry using telecentric imaging

The goal of this article is to discuss 3D Particle Tracking Velocimetry (PTV) in a tomographic reconstructed voxel space with at least doubling the spatial resolution compared to classical 3D PTV. For this purpose, a new tomographic reconstruction technique based on telecentric imaging in combination with the epipolar geometry is presented. The method overcomes the need for memory intensive weighting matrices or cost intensive iterations, which are necessary in iterative algebraic reconstruction techniques. A characteristic of tomographic reconstruction is the reconstruction of ghost particles. As the aim of PTV is the reconstruction of true particle paths, this article focuses on the removal of ghost particles and ghost trajectories. The method is validated via a synthetic turbulent flow field and via the benchmark experiment of a vortex ring.     

A PDF propagation approach to model turbulent dispersion in swirling flows

A computationally efficient approach that solves for the spatial covariance matrix along the dense particle ensemble-averaged trajectory has been successfully applied to describe turbulent dispersion in swirling flows. The procedure to solve for the spatial covariance matrix is based on turbulence isotropy assumption, and it is analogous to Taylor's approach for turbulent dispersion. Unlike stochastic dispersion models, this approach does not involve computing a large number of individual particle trajectories in order to adequately represent the particle phase; a few representative particle ensembles are sufficient to describe turbulent dispersion. The particle Lagrangian properties required in this method are based on a previous study (Shirolkar and McQuay, 1998). The fluid phase information available from practical turbulence models is sufficient to estimate the time and length scales in the model. In this study, two different turbulence models are used to solve for the fluid phase – the standard k–ε model, and a multiple-time-scale (MTS) model. The models developed here are evaluated with the experiments of Sommerfeld and Qiu (1991). A direct comparison between the dispersion model developed in this study and a stochastic dispersion model based on the eddy lifetime concept is also provided. Estimates for the Reynolds stresses required in the stochastic model are obtained from a set of second-order algebraic relations. The results presented in the study demonstrate the computational efficiency of the present dispersion modeling approach. The results also show that the MTS model provides improved single-phase results in comparison to the k–ε model. The particle statistics, which are computed based on the fundamentals of the present approach, compare favorably with the experimental data. Furthermore, these statistics closely compare to those obtained using a stochastic dispersion model. Finally, the results indicate that the particle predictions are relatively unaffected by whether the Reynolds stresses are based on algebraic relations or on the turbulence isotropy assumption.     

SR-CT技术乘型迭代算法的研究

算法为SR-CT(synchrotron radiation computed tomography)技术中的重要组成部分,其中乘型迭代算法为该技术的一种有效算法。本文对乘型迭代算法进行了研究,并对该算法中的一些重要参数,如迭代步长、迭代次数及初始解进行了优化分析,给出了这些参数对重建图像及运算时间的影响关系以及这些参数间的相互影响关系,同时给出了获取较优初始解的方法以及迭代步长和迭代次数的最佳取值范围。     

An all-at-once algebraic multigrid method for finite element discretizations of Stokes problem

We consider numerical solution of finite element discretizations of the Stokes problem. We focus on the transform-then-solve approach, which amounts to first apply a specific algebraic transformation to the linear system of equations arising from the discretization, and then solve the transformed system with an algebraic multigrid method. The approach has recently been applied to finite difference discretizations of the Stokes problem with constant viscosity, and has recommended itself as a robust and competitive solution method. In this work, we examine the extension of the approach to standard finite element discretizations of the Stokes problem, including problems with variable viscosity. The extension relies, on one hand, on the use of the successive over-relaxation method as a multigrid smoother for some finite element schemes. On the other hand, we present strategies that allow us to limit the complexity increase induced by the transformation. Numerical experiments show that for stationary problems our method is competitive compared to a reference solver based on a block diagonal preconditioner and MINRES, and suggest that the transform-then-solve approach is also more robust. In particular, for problems with variable viscosity, the transform-then-solve approach demonstrates significant speed-up with respect to the block diagonal preconditioner. The method is also particularly robust for time-dependent problems whatever the time step size.     

Tomographic shadowgraphy for three-dimensional reconstruction of instantaneous spray distributions

Tomographic shadowgraphy is an image-based optical technique capable of reconstructing the three dimensional instantaneous spray distributions within a given volume. The method is based on a multiple view imaging setup with inline illumination provided by current-pulsed LEDs, which results in droplet shadows being projected onto multiple sensor planes. Each camera records image pairs with short inter-framing times that allow the trajectories of the individual droplets to be estimated using conventional three-dimensional particle tracking approaches. The observed volume is calibrated with a traversed micro-target. A comparison is made between several photogrammetric and polynomial least-square camera calibration techniques regarding their accuracy in deep volume calibration at magnifications close to unity. A calibration method based on volume calibration from multiple planar homographies at equally spaced z-planes was found to produce the most reliable calibration. The combination of back-projected images at each voxel plane efficiently reproduces the droplet positions in three-dimensional space by line-of-sight (LOS) intensity reconstruction. Further improvement of the reconstruction can be achieved by iterative tomographic reconstruction, namely simultaneous multiplicative algebraic reconstruction technique (SMART). The quality of spray reconstruction is investigated using experimental data from multiple view shadowgraphs of hollow cone and flat fan water sprays. The investigations are further substantiated with simulations using synthetic data.     

Lower-Dimensional Manifold (Algebraic) Representation of Reynolds Stress Closure Equations

For complex turbulent flows, Reynolds stress closure modeling (RSCM) is the lowest level at which models can be developed with some fidelity to the governing Navier–Stokes equations. Citing computational burden, researchers have long sought to reduce the seven-equation RSCM to the so-called algebraic Reynolds stress model which involves solving only two evolution equations for turbulent kinetic energy and dissipation. In the past, reduction has been accomplished successfully in the weak-equilibrium limit of turbulence. In non-equilibrium turbulence, attempts at reduction have lacked mathematical rigor and have been based on ad hoc hypotheses resulting in less than adequate models.?In this work we undertake a formal (numerical) examination of the dynamical system of equations that constitute the Reynolds stress closure model to investigate the following questions. (i) When does the RSCM equation system formally permit reduced representation? (ii) What is the dimensionality (number of independent variables) of the permitted reduced system? (iii) How can one derive the reduced system (algebraic Reynolds stress model) from the full RSCM equations? Our analysis reveals that a lower-dimensional representation of the RSCM equations is possible not only in the equilibrium limit, but also in the slow-manifold stage of non-equilibrium turbulence. The degree of reduction depends on the type of mean-flow deformation and state of turbulence. We further develop two novel methods for deriving algebraic Reynolds stress models from RSCM equations in non-equilibrium turbulence. The present work is expected to play an important role in bringing much of the sophistication of the RSCM into the realm of two-equation algebraic Reynolds stress models. Another objective of this work is to place the other algebraic stress modeling efforts in the lower-dimensional modeling context. Received 19 November 1999 and accepted 3 August 2000     

A Subgrid-Scale Stabilized Finite Element Method for Multicomponent Reactive Transport through Porous Media

Standard Galerkin finite element methods (GFEM) lack stability in solving advection-dominated solute transport in porous media. They usually require prohibitively fine grids and extremely small time steps to solve for advection-dominated problems. The algebraic subgrid-scale stabilized (ASGS) finite element method has been proved to overcome such problems for single-species reactive transport. Its potential for dealing with multicomponent reactive transport has not yet been explored. Here we present a numerical formulation of ASGS for steady and transient multicomponent reactive transport. Subgrid-scale transport equations are solved first by using an ASGS approximation and their solutions are substituted back into the grid-scale equations. A sequential iteration approach (SIA) is used to solve for coupled transport and chemical equations. Coupling of ASGS and SIA, ASGS+SIA, has been implemented in a reactive transport code, CORE^2D V4, and verified for conservative solute transport. ASGS+SIA has been tested for a wide range of 1-D transient multicomponent reactive transport problems involving different types of chemical reactions such as: (1) Kinetically controlled aqueous species degradation, (2) Kinetic mineral dissolution, (3) Serial-parallel decay networks, and (4) Cation exchange and pyrite oxidation at local equilibrium. ASGS+SIA always provides accurate solutions and therefore offers an efficient option to solve for advection-dominated multicomponent reactive transport problems.