A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

A novel method for automated grid generation of ice shapes for local‐flow analysis

Modelling a complex geometry, such as ice roughness, plays a key role for the computational flow analysis over rough surfaces. This paper presents two enhancement ideas in modelling roughness geometry for local flow analysis over an aerodynamic surface. The first enhancement is use of the leading‐edge region of an airfoil as a perturbation to the parabola surface. The reasons for using a parabola as the base geometry are: it resembles the airfoil leading edge in the vicinity of its apex and it allows the use of a lower apparent Reynolds number. The second enhancement makes use of the Fourier analysis for modelling complex ice roughness on the leading edge of airfoils. This method of modelling provides an analytical expression, which describes the roughness geometry and the corresponding derivatives. The factors affecting the performance of the Fourier analysis were also investigated. It was shown that the number of sine–cosine terms and the number of control points are of importance. Finally, these enhancements are incorporated into an automated grid generation method over the airfoil ice accretion surface. The validations for both enhancements demonstrate that they can improve the current capability of grid generation and computational flow field analysis around airfoils with ice roughness. Copyright © 2004 John Wiley & Sons, Ltd.     

An improved progressive preconditioning method for steady non‐cavitating and sheet‐cavitating flows

An improved progressive preconditioning method for analyzing steady inviscid and laminar flows around fully wetted and sheet‐cavitating hydrofoils is presented. The preconditioning matrix is adapted automatically from the pressure and/or velocity flow‐field by a power‐law relation. The cavitating calculations are based on a single fluid approach. In this approach, the liquid/vapour mixture is treated as a homogeneous fluid whose density is controlled by a barotropic state law. This physical model is integrated with a numerical resolution derived from the cell‐centered Jameson's finite volume algorithm. The stabilization is achieved via the second‐and fourth‐order artificial dissipation scheme. Explicit four‐step Runge–Kutta time integration is applied to achieve the steady‐state condition. Results presented in the paper focus on the pressure distribution on hydrofoils wall, velocity profiles, lift and drag forces, length of sheet cavitation, and effect of the power‐law preconditioning method on convergence speed. The results show satisfactory agreement with numerical and experimental works of others. The scheme has a progressive effect on the convergence speed. The results indicate that using the power‐law preconditioner improves the convergence rate, significantly. Copyright © 2010 John Wiley & Sons, Ltd.     

基于多重网格分层四边形单元有限元空间的预处理方法

用有限元方法解二阶问题时,常选用三角单元和四边形单元.对于基于三角形单元的有限元方程,采用多重网格分层函数的预处理方法,Yserentant,H.已作了详细的讨论.本文对基于四边形单元的有限元方程,采用类似于Yserentant,H.的多重网格分层函数的预处理方法,使得系数矩阵的条件数由原来的O((1/h)~2)缩小为O((log(1/h)~3).     

Anomalous electronic transport in quasicrystals and related complex metallic alloys

We analyse the transport properties in approximants of quasicrystals α-AlMnSi, 1/1-AlCuFe and for the complex metallic phase λ-AlMn. These phases present strong analogies in their local atomic structures and are related to existing quasicrystalline phases. Experimentally, they present unusual transport properties with low conductivities and a mix of metallic-like and insulating-like characteristics. We compute the band structure and the quantum diffusion in the perfect structure without disorder and introduce simple approximations that allow us to treat the effect of disorder. Our results demonstrate that the standard Bloch–Boltzmann theory is not applicable to these intermetallic phases. Indeed their dispersion relations are flat, indicating small band velocities, and corrections to quantum diffusion, which are not taken into account in the semi-classical Bloch–Boltzmann scheme, become dominant. We call this regime the small velocity regime. A simple relaxation time approximation to treat the effect of disorder allows us to reproduce the main experimental facts on conductivity qualitatively and even quantitatively.     

Rapidly convergent series for high-accuracy calculation of the Voigt function

A rapidly convergent series, based on Fourier expansion of the exponential multiplier, is presented for highly accurate approximation of the Voigt function (VF). The corresponding algorithm enables the rapid calculation, required for its implementation as a subprogram in an interpolation approach. The numerical analysis of this VF approximation suggests that it may be more accurate than 10⁻⁹ in the Humlí?ek regions 3 and 4.     

Synthesis methods of multiple phase-cycled SSFP images to reduce the band artifact and noise more reliably

The band artifact in steady-state free precession can be reduced by synthesizing the multiple images obtained through different phase increments of successive radiofrequency pulses. Even though the complex summation method was reported to be effective in reducing the band artifact, it has the pitfalls of intensity abnormality and sensitivity to the phase abnormality. Two new methods have been developed for more reliable reduction of the band artifact than the complex summation method. One method is to sum the complex images partially and to take the maximum intensity of the partially summed images. The other method is to sum the free induction decay (FID) and primary echo components of the Fourier series that are obtained through Fourier analysis of the complex base images. Both proposed methods were compared with other magnitude (maximum intensity projection, spectrally decomposed synthesis, sum-of-squares, nonlinear averaging) and complex-based (complex summation, magnitude-weighted complex summation) methods experimentally at 3 T for the phantom and volunteer's head imaging. Both proposed methods were confirmed to maintain the advantage of the complex summation in reducing both the dark and bright band artifacts while reducing the intensity abnormality and sensitivity to the phase abnormality from that of the complex summation method over a wide range of flip angles and relaxation times.     

复杂混合溶液成分高光谱分析的可行性

采用高光谱技术对复杂混合溶液进行检测分析,同时利用被测物质的吸光度和散射特性信息以提高光谱的信噪比。实验设计了高光谱采集装置,采集生物组织模拟液(Intralipid-10%)的漫反射高光谱图像,并用Monte Carlo方法和漫射近似理论对其进行了正向和反向推导,获得了632nm波长下,Intralipid-10%吸收系数为0.002 0cm-1,与标准参数相对误差为11.1%;约化散射系数为63.35cm-1,与标准参数相对误差为6.49%,基本符合标准参数的误差范围,验证了该高光谱检测系统的准确性。还利用该高光谱系统对不同厂家出品的牛奶、果汁等样本进行了高光谱采集,得到不同样本间差异较传统二维光谱更为明显的结果,充分证明了高光谱方法在复杂混合溶液成分分析中具有很强的可行性。     

The equilibrium limit of a constitutive model for two-phase granular mixtures and its numerical approximation

In this paper we analyze the equilibrium limit of the constitutive model for two-phase granular mixtures introduced in Papalexandris (2004) [13], and develop an algorithm for its numerical approximation. At, equilibrium, the constitutive model reduces to a strongly coupled, overdetermined system of quasilinear elliptic partial differential equations with respect to the pressure and the volume fraction of the solid granular phase. First we carry a perturbation analysis based on standard hydrostatic-type scaling arguments which reduces the complexity of the coupling of the equations. The perturbed system is then supplemented by an appropriate compatibility condition which arises from the properties of the gradient operator. Further, based on the Helmholtz decomposition and Ladyzhenskaya’s decomposition theorem, we develop a projection-type, Successive-Over-Relaxation numerical method. This method is general enough and can be applied to a variety of continuum models of complex mixtures and mixtures with micro-structure. We also prove that this method is both stable and consistent hence, under standard assumptions, convergent. The paper concludes with the presentation of representative numerical results.     

基于正弦光栅调制的相位恢复

提出一种基于正弦光栅调制的相位恢复方法.首先,在空间光调制器上加载两个不同方向的正弦光栅对光场复振幅进行调制;然后,推导强度信息和相位偏导数之间的关系式;最后由强度信息计算相位偏导数,进而恢复相位.实验分别测试了该方法对简单图像、复杂图像以及含噪声图像的恢复能力.实验结果表明,该方法可以通过方便地控制图像移动量来有效恢复相位信息,适用于不同尺度的相位物体的计算.