Unified semi‐analytical wall boundary conditions for inviscid,laminar or turbulent flows in the meshless SPH method

Wall boundary conditions in smoothed particle hydrodynamics (SPH) is a key issue to perform accurate simulations. We propose here a new approach based on a renormalising factor for writing all boundary terms. This factor depends on the local shape of a wall and on the position of a particle relative to the wall, which is described by segments (in two‐dimensions), instead of the cumbersome fictitious or ghost particles used in most existing SPH models. By solving a dynamic equation for the renormalising factor, we significantly improve traditional wall treatment in SPH, for pressure forces, wall friction and turbulent conditions. The new model is demonstrated for cases including hydrostatic conditions for still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. The latter case is also compared with a finite volume and volume‐of‐fluid scheme. The performance of the model for a two‐dimensional laminar flow in a channel is demonstrated where the profiles of velocity are in agreement with the theoretical ones, demonstrating that the derived wall shear stress balances the pressure gradient. Finally, the performance of the model is demonstrated for flow in a schematic fish pass where both the velocity field and turbulent viscosity fields are satisfactorily reproduced compared with mesh‐based codes. Copyright © 2012 John Wiley & Sons, Ltd.     

An immersed boundary method wall model for high‐Reynolds‐number channel flow over complex topography

High‐Reynolds‐number channel flows regularly encounter topographies composed of multiple length scales and that protrude into the boundary layer. Physically, the presence of immersed obstacles leads to increased velocity gradients, turbulence production, and manifestation of wakes. Considerable challenges are associated with numerically describing the presence of obstacles in channel flows. Common approaches include generation of a computational mesh that is uniquely designed for the flow and obstacle, the immersed boundary method, and terrain‐following coordinates. There are challenges and limitations associated with each of these techniques. Specification of boundary conditions representing the perimeter of solid obstacles is a primary challenge of the immersed boundary method. In this document, a simplistic canopy stress‐like wall model is used to impose boundary conditions. The model isolates aerodynamically relevant local frontal areas through evaluation of the gradient of the topographic height field. The gradient of the height field describes both the surface‐normal direction and the frontal area, making it ideal for detecting areas on which the flow impinges. The model is tested in numerical simulations of turbulent half‐channel flow over topographies with different obstacles affixed–right prisms, rectangular prisms, ellipsoidal mounds, and sinusoids. In all cases, the performance is strong relative to datasets presented in the literature. Results are finally presented for numerical simulation of flow over complex synthetic fractal‐like topography and a synthetic city. These results show interesting trends in how the turbulent multiscale flow field responds to multiscale topography. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison between computational fluid dynamics,fluid–structure interaction and computational structural dynamics predictions of flow-induced wall mechanics in an anatomically realistic cerebral aneurysm model

Haemodynamically induced stress plays an important role in the progression and rupture of cerebral aneurysms. The current work describes computational fluid dynamics (CFD), fluid–structure interaction (FSI) and computational structural dynamics (CSD) simulations in an anatomically realistic model of a carotid artery with two saccular cerebral aneurysms in the ophthalmic region. The model was obtained from three-dimensional (3D) rotational angiographic imaging data. CFD and FSI were studied under a physiologically representative waveform of inflow. The arterial wall was assumed elastic or hyperelastic, as a 3D solid or as a shell depending on the type of modelling used. The flow was assumed to be laminar, non-Newtonian and incompressible. The CFD, FSI and CSD models were solved with the finite elements package ADINA. Predictions of velocity field and wall shear stress (WSS) on the aneurysms made using CFD and FSI were compared. The CSD model of the aneurysms using complete geometry was compared with isolated aneurysm models. Additionally, the effects of hypertensive pressure on CSD aneurysm models are also reported. The vortex structure, WSS, effective stress, strain and displacement of the aneurysm walls showed differences, depending on the type of modelling used.     

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

LES of turbulent flow in a concentric annulus with rotating outer wall

Fully-developed turbulent flow in a concentric annulus, r₁/r₂ = 0.5, Re_h = 12,500, with the outer wall rotating at a range of rotation rates N = U_θ,wall/U_b from 0.5 up to 4 is studied by large-eddy simulations. The focus is on the effects of moderate to very high rotation rates on the mean flow, turbulence statistics and eddy structure. For N up to ∼2, an increase in the rotation rate dampens progressively the turbulence near the rotating outer wall, while affecting only mildly the inner-wall region. At higher rotation rates this trend is reversed: for N = 2.8 close to the inner wall turbulence is dramatically reduced while the outer wall region remains turbulent with discernible helical vortices as the dominant turbulent structure. The turbulence parameters and eddy structures differ significantly for N = 2 and 2.8. This switch is attributed to the centrifuged turbulence (generated near the inner wall) prevailing over the axial inertial force as well as over the counteracting laminarizing effects of the rotating outer wall. At still higher rotation, N = 4, the flow gets laminarized but with distinct spiralling vortices akin to the Taylor–Couette rolls found between the two counter-rotating cylinders without axial flow, which is the limiting case when N approaches to infinity. The ratio of the centrifugal to axial inertial forces, Ta/Re² ∝ N² (where Ta is the Taylor number) is considered as a possible criterion for defining the conditions for the above regime change.     

Chaotic fractional‐order model for muscular blood vessel and its control via fractional control scheme

This article studies the chaotic and complex behavior in a fractional‐order biomathematical model of a muscular blood vessel (MBV). It is shown that the fractional‐order MBV (FOMBV) model exhibits very complex and rich dynamics such as chaos. We show that the corresponding maximal Lyapunov exponent of the FOMBV system is positive which implies the existence of chaos. Strange attractors of the FOMBV model are depicted to validate the chaotic behavior of the system. We change the fractional order of the model and investigate the dynamics of the system. To suppress the chaotic behavior of the model, we propose a single input fractional finite‐time controller and prove its stability using the fractional Lyapunov theory. In addition, the effects of the model uncertainties and external disturbances are taken into account and a robust fractional finite‐time controller is constructed. The upper bound of the chaos suppression time is also given. Some computer simulations are presented to illustrate the findings of this article. © 2014 Wiley Periodicals, Inc. Complexity 20: 37–46, 2014     

基于有限分析法的平面铅垂壁羽流的相似解

应用计及浮力对湍动能及其耗散率的影响的k-epsilon湍流模式,结合有限分析法对密度差引起的平面铅垂紊动壁羽流进行了数值分析．在均匀环境条件下控制该类流动的连续性方程、流动方向的动量方程、浓度扩散方程、湍动能及其耗散率方程存在相似解．考虑到浮力通量守恒条件,应用有限分析法给出了铅直壁羽流的速度、相对密度差、湍动能及耗散率的分布,进而给出了各物理量最大值沿主流方向变化的关系式．湍流Schmidt数为1.0时的计算结果与实验资料吻合较好,表明应用有限分析法分析铅垂平面羽流是有效的,即在分析壁羽流时浮力对湍动能及其耗散率的影响应该予以考虑．     

A combined analytical approach applied to Medieval wall paintings from Puglia (Italy): The study of painting techniques and its conservation state

A combined analytical approach has been applied to the wall paintings, dated from 10th to 14th centuries, of the Santi Stefani crypt at Vaste (Lecce, Southern Italy). These paintings are a precious testimony of Medieval art in Southern Italy. However, the church shows problems of damp as well as clear evidences of flora, fungi and mold presence, and there is little knowledge of the pictorial methodologies used. Raman spectroscopy allowed to determine the palette and to reconstruct the worksite and the chronological sequence of the various paint layers. Kaolinite, calcite, carbon black, hematite, massicot, goethite, indigo and azurite were identified as pigments along with synthetic pigments, like phthalocyanine blue and chrome yellow. Attenuated total reflectance‐Fourier transform infrared spectroscopy suggested the presence of egg as a binder in some pictorial layers. The conservation state of the crypt is poor, and detachments of pigmented layers are frequent because of the presence of subflorescence and efflorescence: nitrate, sulfate and chloride salts have been identified spectroscopically and quantified by ion chromatography. The extensive use of kaolinite in Santi Stefani, actually not uncommon in Medieval art, is observed for the first time in a crypt of Puglia: its use to stabilize some pigments and to improve their adhesion on substrate is proposed. Copyright © 2015 John Wiley & Sons, Ltd.     

贴壁风对脱硫反应器内液滴粘壁的影响

反应器内液滴粘壁对半干法烟气脱硫效率和运行稳定性具有重要影响.本文采用数值方法研究了不同贴壁风风速和风口高度对反应器内液滴粘壁现象的影响规律.结果表明,贴壁风对减轻液滴粘壁具有一定作用.贴壁风风速越高,液滴贴壁数量越少.贴壁风衰减较快,其风口高度存在最佳值,贴壁风风口与液滴喷射口的高度距离是按照液滴直线运动粘壁时液滴高度行程的1.2倍.     

基于超声红外技术对金属管内壁缺陷的检测

为提高管道的运送效率,及时检测排除管道内壁的缺陷非常重要。提出利用超声红外无损检测方法对管道内部进行检测。超声主动热激励试件,高频红外热像仪记录试件表面温度变化,结合了超声摩擦生热和红外热成像的优点。对壁厚约为3.3mm检测难度较大的金属管内壁缺陷进行了检测。通过对采集数据和热图的处理分析,准确快速地确定缺陷所在的位置。实物对比分析表明：超声热激励红外无损检测技术可以对金属管道内壁缺陷进行准确检测定位。