A flexible multiscale algorithm based on an improved smoothed particle hydrodynamics method for complex viscoelastic flows

Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.

NUMERICAL SIMULATION OF MORPHOLOGY OF POLYMER CHAIN COILS IN COMPLEX FLOWS

A new coupled finite element formulation is proposed to calculate a conformation tensor model in two complex flows： a planar contraction flow and a planar flow around a symmetrically placed cylinder. The components of conformation tensor are first computed together with the velocity and pressure to describe the change of morphology of polymer chain coils in flow fields. Macroscopic quantities of viscoelastic flow are then calculated based on the conformation tensor. Comparisons between the numerical simulations and experiments for stress patterns and velocity profiles are carried out to prove the validity of the method.

Experimental Study of Sediment Incipience Under Complex Flows

Sediment incipience under flows passing a backward-facing step was studied. A series of experiments were conducted to measure scouring depth, probability of sediment incipience, and instantaneous flow velocity field downstream of a backward-facing step. Instantaneous flow velocity fields were measured by using Particle Image Velocimetry (PIV), and an image processing method for determining probability of sediment incipience was employed to analyze the experimental data. The experimental results showed that the probability of sediment incipience was the highest near the reattachment point, even though the near-wall instantaneous flow velocity and the Reynolds stress were both much higher further downstream of the backward-facing step. The possible me- chanisms are discussed for the sediment incipience near the reattachment point.

Hybrid Methods for Computing the Streamfunction and Velocity Potential for Complex Flow Fields over Mesoscale Domains

Three types of previously used numerical methods are revisited for computing the streamfunctionψand velocity potentialχfrom the horizontal velocity v in limited domains.The first type,called the SOR-based method,uses a classical successive over-relaxation(SOR)scheme to computeψ(orχ)first with an arbitrary boundary condition(BC)and thenχ(orψ)with the BC derived from v.The second type,called the spectral method,uses spectral formulations to construct the inner part of(ψ,χ)-the inversion of(vorticity,divergence)with a homogeneous BC,and then the remaining harmonic part of(ψ,χ)with BCs from v.The third type,called the integral method,uses integral formulas to compute the internally induced(ψ,χ)-the inversion of domain-internal(vorticity,divergence)using the free-space Greenꞌs function without BCs and then the remaining harmonicψ(orχ)with BCs from v minus the internally-induced part.Although these methods have previously been successfully applied to flows in large-scale and synoptic-scale domains,their accuracy is compromised when applied to complex flows over mesoscale domains,as shown in this paper.To resolve this problem,two hybrid approaches,the integral-SOR method and the integral-spectral method,are developed by combining the first step of the integral method with the second step adopted from the SOR-based and spectral methods,respectively.Upon testing these methods on real-case complex flows,the integral-SOR method is significantly more accurate than the integral-spectral method,noting that the latter is still generally more accurate than the three previously-used methods.The integral-SOR method is recommended for future applications and diagnostic studies of complex flows.

Advances in LES of Two-phase Combustion(Ⅱ) LES of Complex Gas-Particle Flows and Coal Combustion

Large-eddy simulation(LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering.Large-eddy simulation of two-phase flows and combustion is particularly important for engineering applications.Some investigators,including the present authors,give their review on LES of spray combustion in gas-turbine combustors and internal combustion engines.However,up to now only a few papers are related to the state-of-the-art on LES of gas-particle flows and combustion.In this paper a review of the advances in LES of complex gas-particle flows and coal combustion is presented.Different sub-grid scale(SGS) stress models and combustion models are described,some of the main results are summarized,and some research needs are discussed.

Application Progress of Computational Fluid Dynamic Techniques for Complex Viscous Flows in Ship and Ocean Engineering

Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.

AN EXPLORATION TO THE MODELLING OF TWO-DIMENSIONAL COMPLEX TURBULENT SHEAR FLOWS

The regions with shear stress and mean velocity gradient of opposite sign often exist in complex turbulent shear flows.In these cases,the eddy viscosity hypothesis breaks down.Hinze regards the,departure from eddy viscosity hypothesis as a result from transportation of mean momentum over distance by the large structures and arrives at a shear stress expression including the second order derivatives of the mean velocity.However,his expression greatly overestimates the shear stress.This implies that the flow particles are unlikely to have enough memory of the mean momentum over distance.By assuming the departure from eddy viscosity hypothesis as a result from transportation of the shear stress contained in smaller eddies over distance by the large structures,the present author has arrived at a new shear stress expression.The shear stress estimated so far is in good agreement with the experiments.

Investigation on Liquid Holdup in Vertical Zero Net-Liquid Flow

Zero net-liquid flow (ZNLF) is a special case of upward gas-liquid two-phase flow. It is a phenomenon observed as a gas-liquid mixture flows in a conduit but the net liquid flow rate is zero. Investigation on the liquid holdup of ZNLF is conducted in a vertical ten-meter tube with diameter of 76 mm, both for Newtonian and nonNewtonian fluids. The gas phase is air. The Newtonian fluid is water and the non-Newtonian fluids are water-based guar gel solutions. The correlations developed for predicting liquid holdup on the basis of Lockhart-Martinelli parameter are not suitable to ZNLF. A constitutive correlation for the liquid holdup of vertical ZNLF was put forward by using the mass balance. It is found that the liquid holdup in ZNLF is dependent on both the gas flow rate and the flow distribution coefficient.

DIGITAL ANALYSIS TECHNOLOGY FOR MORPHOLOGY OF POLYMER CHAIN COILS IN FLOW FIELDS

Polymer chain coils with entanglement is a crucial scale of structures in polymer materials since their relaxationtimes are matching practical processing times.Based on the phenomenological model of polymer chain coils and a new finiteelement approach,we have designed a computer software including solver,pre-and post-processing modules,and developeda digital analysis technology for the morphology of polymer chain coils in flow fields(DAMPC).Using this technology wemay simulate the morphology development of chain coils in various flow fields,such as simple shear flow,elongational flow,and any complex flow at transient or steady state.The applications made up to now show that the software predictions arecomparable with experimental results.

Laser Doppler Measurements of Twin Impinging Jets Aligned with a Crossflow

This paper presents a detailed analysis of the complex flow beneath two impinging jets aligned with a low-velocity crossflow which is relevant for the future F-35 VSTOL configuration, and provides a quantitative picture of the main features of interest for impingement type of flows. The experiments were carried out for a Reynolds number based on the jet exit conditions of Rej = 4.3 × 10^4, an impingement height of 20.1 jet diameters and for a velocity ratio between the jet exit and the crossflow VR = V/Uo of 22.5. The rear jet is located at S = 6 D downstream of the first jet. The results show a large penetration of the first （upstream）jet that is deflected by the crossflow and impinges on the ground, giving rise to a ground vortex due to the collision of the radial wall and the crossflow that wraps around the impinging point like a scarf. The rear jet （located downstream） it is not so affected by the crossflow in terms of deflection, but due to the downstream wall jet that flows radially from the impinging point of the first jet it does not reach the ground. The results indicate a new flow pattern not yet reported so far, that for a VSTOL aircraft operating in ground vicinity with front wind or small forward movement may result in enhanced under pressures in the aft part of the aircraft causing a suction down force and a change of the pitching moment towards the ground.

Turbulent Energy Budgets of a Ground Vortex Flow

Turbulent kinetic energy budgets are presented for a highly curved flow generated by the collision of plane wall turbulent jet with a low-velocity boundary layer. The different terms are obtained in the vertical plane of symmetry by quadratic interpolation of the LDV （Laser Doppler Velocimetry） measurements, for a wall jet-to-boundary layer velocity ratio of 2. The results, which have relevance to flows encountered in powered-lift aircraft operating in ground effect, quantify the structure of the complex ground vortex flow. The analysis of turbulent energy equation terms using the measured data revealed that production by normal and shear stresses are both very important to the turbulent structure of the impact zone of the ground vortex. This is an indication that the modeling of turbulence of a ground vortex requires a good representation of the production by normal stresses which is most important in the collision zone.

终端区交通复杂性评估及演变分析

终端区内飞行密集、环境多变,已成为国家空域系统运行能力提升的重要瓶颈。精准量化终端区内交通复杂性并把握其演化规律是实现空中交通精细化管理的基础,对于提升终端区运行能力有重要意义。终端区空中交通复杂性描述了某一时间点或时间段内空中交通流特征的复杂程度。基于动态密度思想,针对终端区运行特征,构建了多维度的空中交通复杂性评价指标,并运用熵权法确定指标权重,反映评估指标的无序化程度,评定信息量的大小。构建复杂度计算模型,从客观上利用多指标反映航空器实时复杂性演变。采用多智能体仿真建模工具NetLogo搭建了终端区仿真环境,以广州终端区的实际进场数据搭建了基础场景,模拟了终端区交通复杂性的宏观演化过程。通过调节进场移交点交通流需求比例、设计不同天气情景,分析了不同场景下交通复杂性演化特征。结果表明,该评估体系使用较少的指标和计算时间,能够较为准确、客观地评估时间片内的复杂性,符合实际情况。NetLogo建模工具能够仿真终端区交通流,直观地观察交通流运行行为。当进场移交点交通流流量分布更为均匀时,复杂性较低。终端区受恶劣天气影响时,其受影响程度取决于恶劣天气覆盖航段的航班流量分布情况、影响的范围以及受影响航段到跑道距离。部分交通流存在迟滞现象,反映了复杂性的蔓延时空特征。该复杂性评估方法能够有效评估航段、航线、终端区复杂性及演变规律,为开展终端区精细化空中交通管理提供参考。

面向任务的重叠联盟结构生成计算复杂性

传统的重叠联盟形成问题大都聚焦智能体,鲜有从任务视角出发.为此,本文首先构建了一种面向任务的重叠联盟结构生成模型,并分析了其解空间和相关决策问题的计算复杂性.此外,基于流网络分别设计了相应的孤立联盟、重叠联盟、重叠联盟结构成功性判别算法和最优重叠联盟结构生成算法.分析结果表明,判别孤立联盟、重叠联盟、重叠联盟结构的成功性的时间复杂度均与智能体数和任务数呈多项式关系,而搜索最优重叠联盟结构的时间复杂度与智能体数和任务数呈指数关系.最后,通过仿真实验验证了上述结果.

Complex granular flows of sticky-wet material on flip-flow screens:Calibration of discrete element simulations

Discrete element method(DEM)is an effective approach for studying the screening process of flip-flow screens.However,there have been few studies focusing on the thick layer of sticky-wet particles on flip-flow screens.To achieve accurate simulations of the thick layer of sticky-wet particles on a flip-flow screen,firstly,the movement law of particle flow was studied,and a multi-regime combination cali-bration method based on characteristics of particle flow regimes was proposed.Based on the Plackett-Burman experiment,the curse of dimensionality caused by multi-state and multi-contact parameters was overcome.Subsequently,the lifting cylinder,rotating drum,and trampoline tests were carried out to obtain macroscopic reference values under various granular flow regimes.The calibration results were then determined using the response surface method and climbing algorithm.Finally,the calibration results were tested at both macroscopic and mesoscopic scales and compared with a commonly used calibration method.The results demonstrate that the calibration method,which considers the multi-state characteristics,improves simulation accuracy by 2%-10%and reduces the simulation error to less than 10%,thus meeting the requirements for engineering optimization of flip-flow screens.

基于潮流和复杂网络理论的综合能源系统关键环节识别

综合能源系统(integrated energy system, IES)的关键环节识别对于系统的脆弱性研究和改善有重要意义。针对传统还原论方法无法解释系统层面且复杂网络分析中未考虑能流带来影响的问题,提出一种基于潮流和复杂网络结构的IES关键环节识别方法。首先,建立基于潮流的综合能源复杂网络模型,将IES的数学模型参数和潮流结果引入复杂网络的建模以使模型更加贴合现实的运行工况;其次,基于潮流计算结果和复杂网络参数,提出节点能量度和能量边介数2个评价指标,提高了关键环节识别精度;最后,根据所建立的模型对关键环节进行破坏并观察网络效率的变化,与其他攻击模式进行对比分析,验证了所提方法的可行性。

ENERGY FLOWS IN COMPLEX ECOLOGICAL SYSTEMS:A REVIEW

Energy flow drives the complex systems to evolve.The allometric scaling as the universalenergy flow pattern has been found in different scales of ecological systems.It reflects the generalpower law relationship between flow and store.The underlying mechanisms of energy flow patterns areexplained as the branching transportation networks which can be regarded as the result of systematicoptimization of a biological target under constraints.Energy flows in the ecological system may bemodelled by the food web model and population dynamics on the network.This paper reviews thelatest progress on the energy flow patterns,explanatory models for the allometric scaling and modellingapproach of flow and network evolution dynamics in ecology.Furthermore,the possibility of generalizingthese flow patterns,modelling approaches to other complex systems is discussed.

高流速复杂水域下的海底输水管道施工技术

以舟山市大陆引水三期工程(金塘岛引水工程)为例,对浅水铺管船增配的锚泊系统定位能力进行了理论分析和计算,同时介绍了高流速复杂水域环境下的长距离海底输水管道铺设总体施工工艺和关键施工技术,以期为类似工程提供借鉴。

基于复杂系统理论的天津市环境与经济关系分析

基于区域环境与经济系统开放复杂性特点,利用最大流原理建立区域环境与经济复杂系统模型分析二者间的关系.以天津市为研究对象,研究了其环境与经济系统在2001—2010年的发展过程,借助Matlab R2010a求ξ(系统稳定性系数),通过ξ的变化规律分析各组元间的相互作用和系统的整体发展状况.结果表明:2001—2010年天津市环境与经济系统中各组元间的相互作用逐年增强,ξ逐年上升并在2010年达到最大值(18.615 4),表明该系统中环境与经济各组元间的相互作用逐渐协调,天津市的环境与经济发展呈相互促进的关系.

汽轮机低压排汽系统气动性能分析

本文描述了来流条件和几何形状对大功率汽轮机低压排汽系统内部复杂三维流动的影响,研究了均匀和畸变两种来流条件下某排汽系统的气动性能;根据实际运行环境下高性能透平排汽系统的气动设计原则设计了新的排汽系统,用数值模拟的方法对新设计进行了气动性能分析并与原设计进行了对比。

载机军舰回收舰载机的运行模型

对良好气象条件下载机军舰回收舰载机的模型进行研究,以便合理安排在回收阶段发生舰载机复飞、逃逸时空中交通,保证载机军舰回收舰载机的流量与等待航线、着舰航线、甲板阻拦区及甲板停机区等各阶段的容量相适应。对回收舰载机流程进行研究,基于系统动力学理论,建立影响空中交通流量管理的存量流量图及运行模型,并以美国载机军舰回收舰载机为例对该运行模型进行了验证。运行结果表明,本运行模型在舰载机复飞、逃逸的条件下能进行载机军舰空中交通流量管理的趋势预测,可为空中交通流量管理方案的制定提供技术基础。