1.Introduction

Rotorcraft plays an important role in aviation for its unique ability in hovering and vertical descending flight.However,an accurate prediction of the rotor flow field remains a challenging problem.The three-dimensional rotor wake is unsteady and complex.In a hovering condition,highly energetic tip vortices shed from the blade tip region,swirl downward,and then undergo vortex breakdown in the far wake.In descending flight,tip vortices constantly persist near the rotor disk and interact with the blades,which may cause a high level of fuselage vibration and remarkable induced power consumption.Moreover,when a rotorcraft is operated over certain ranges of descent rate,convection of the wake would be inhibited.It results in a doughnut-shaped ring around the rotor disk,which is known as the Vortex Ring State1(VRS,see Fig.1).The VRS is regarded as the roughest flight regime of a rotorcraft.According to statistical data,2at least 32 helicopter accidents were attributed to this dangerous regime between 1982 and 1997.

Over the past few decades,studies about the unsteady characteristics of rotor wake in hovering and descending flight have been focused on experiments3–11and qualitative analysis.12–15Meanwhile,relatively sparse numerical simulations have been paralleled.The aerodynamics near the blade surface,like compressible and viscous effects,are predicted well by conventional Eulerian-based CFD methods, but they are computationally expensive,and inherent numerical dissipation makes the rotor wake diffuse too early.16To date,conventional CFD methods are insufficient in far wake capturing of the rotor.

Lagrangian-based models can address the problem of nonphysical wake diffusion,and they are more computationally efficient than CFD.In recent years,coupled Eulerian/Lagrangian simulation methods have shown promise in rotor wake simulations and received much attention.17–20However,simple Lagrangian models(e.g.,prescribed wake models21and free wake models22)rely heavily on empirical parameters such as vortex core size and decay factor.Furthermore,they cannot provide detailed information of the wake structure.Due to the progress in wake modeling techniques,more advanced Lagrangian models have been developed.Those models are referred to as Vorticity Transport Models(VTMs).23–25They can explicitly conserve wake vorticity without any artificial dissipation and cancel the restriction of empirical parameters.

The main work of the present study is to couple a novel Lagrangian-based VTM model with an Eulerian-based CFD solver.The hybrid solver proposed in this paper combines the merits of CFD and the VTM.A CFD module is used to resolve the compressible blade near body aerodynamics and a VTM module is used to predict the complex wake convection.Details of the CFD/VTM hybrid solver are described in Section 2.Numerical simulations of rotorcraft in hovering and vertical descending flight are performed in Section 3.Results show good predictions of both the blade near body aerodynamics and the wake structure in hover.Induced in flow and time history of rotor thrust in descending flight are also investigated.Main conclusions are summarized in Section 4.

2012年，在巴西的永续发展大会上，各国代表轮番发言，探讨地球环境和人类社会的永续发展。轮到乌拉圭总统穆希卡时，这位捐出大半收入、住在郊区农场、被人民亲切地称为“佩佩”的总统，用质朴的语言，给出了一场直抵人心的演讲。“我们追求更加方便和富裕的生活，可是我们幸福吗？”这是一册让我们深刻反思幸福内涵的绘本。

2.CFD/VTM model description

Fig.1 Schematic of the air flow for a rotorcraft in the vortex ring state.

Fig.2 Schematic of the computational zone.

The computational zone is decomposed into two domains(Fig.2):the blade body- fitted Eulerian domain,which covers a relatively small region near the blade and follows a C-O topology grid(Fig.3),and the background Lagrangian domain,which employs a set of particles to model the wake vorticity.These two domains are solved by CFD and a VTM,respectively.Since the VTM only solves incompressible flow,the C-O grid of the Eulerian domain extends far enough(over two chord lengths)to ensure that outside the grid,air compressibility could be neglected.

2.1.CFD solution of blade body- fitted Eulerian domain

The simulation of the rotor blade near flow field follows the Eulerian description,which relates to the grid-based solutions of compressible Reynold Averaged Navier-Stokes(RANS)equations.RANS equations are solved in terms of the conservation forms of mass,momentum,and energy,and can be written in tensor form as26

Fig.3 Rotor blade near body C-O topology grid.

where ρ,uj,p, τij,E,and qjare the fluid density,the velocity,the pressure,the stress tensor,the total energy per unit mass,and the heat flux,respectively.

In general,the blade rotates at a constant speed,so Eqs.(1)–(3)can be expressed as

2.3.2.Transmission of boundary conditions from VTM to CFD

where u，v，w are the velocity components.

A structured,cell-centered finite volume scheme is employed for CFD calculation and to model the formation of blade tip vortices.The spatial discretization is based on Roe’s flux-difference splitting scheme.27Convective fluxes are evaluated at the face of the control volume by solution reconstruction from the cell center,e.g.,

主产区小麦价格持续走强，尤其华北局部价格已经逼近政策“天花板”。考虑到政策调整及制粉企业的接受能力，后期麦价持续上涨的几率已经不大。

where ARoeis the Roe matrix,and the Spalart-Allmaras oneequation model is chosen for turbulence modeling.

The dual time stepping technique is applied for unsteady simulation,and the implicit Lower-Upper Symmetric Gauss-Seidel(LU-SGS)scheme is adopted for pseudo time discretization.This discretization scheme is defined by

where D is the diagonal matrix,L the strictly lower triangular matrix,U the strictly upper triangular matrix,and R the residual.The implementation of the LU-SGS scheme consists of the following two steps:

利用河流行洪能力复核计算，将导流建筑、堤坝施工进行设计。施工需要在围堰维护下进行，首先对岸边进行低水围堰。采用袋装沙土进行围堰叠筑，叠筑要整齐、密实。利用基础底设计的坡度进行基坑排水，在下游集中抽排。排水系统必须整段贯穿，并有一定的备用量，保证排水的有效性。

Forward sweep:

Backward sweep:

以孝与廉应为做人根本、基础与底线开展道德教育。如果一个人不慈不孝、不廉不洁，可以说其还没有成为真正的人，更不用说是一个好人、一个健康的人，这种伦理就是一种人之为人的底线伦理。西方伦理学所谓的底线伦理或者本底伦理，就是所有人之所以为人都要遵循的东西。以孝廉文化作为底线伦理开展道德教育也体现了社会主义核心价值观的要求。

where γ is the ratio of specific heat coefficient,and Ma is the Mach number.

2.2.VTM solution of background Lagrangian domain

There are two categories of VTMs:one is grid-based and uses an adaptive Cartesian grid to model the rotor wake,23,24while in the other category,wake vorticity is represented by a set of viscous vortex particles.25In the current study,the grid-free Lagrangian description that utilizes vortex particles is adopted.The governing vorticity dynamics equation is written as follows:

The above equation is derived from incompressible Navier-Stokes equations and is written in a vorticity-velocity form.The first term on the right-hand side ω ·∇u is the so-called vortex stretching term,which denotes the dilation and rotation effects of the vortex filaments within the flow.The second term ν∇2ω is the viscous diffusion term,which expresses the effects of air viscosity.In Eq.(10),ω=∇×u is the vorticity,denotes the material derivative,ν is the air kinematic viscosity coefficient,and the source S is the vorticity shed and trailed from each rotor blade segment as

where Γbis the blade bound circulation vector,and ubis the local velocity of the blade segment relative to the air.

Assume that at time t and position x,the vorticity of the ith particle is defined by

whereis the cutoff function,ε the smoothing parameter,ξ(|x|)the distribution function,and a Gaussian distribution is adopted here.For the non-dimensional distance parameter x between the target point i and the source point j,the Gaussian distribution is as follows:

The velocity vector is computed by the application of the Poisson equation as

Substituting Green’s function G(x′)into Eq.(14)yields the generalized Biot-Savart equation28as

Green’s function is given by

where the error function is written as

7.Xin H,Gao ZH.A prediction of the vortex-ring state boundary based on model tests.Trans Nanjing Univ Aeronaut Astronaut 1994;11(2):159–64.

where subscript ‘st” denotes the vortex ‘stretching term”.Eq.(18)represents a dense algebraic system.If the Lagrangian domain consists of N vorticity-containing particles,25the task of calculating all the pairwise interactions by invoking Eq.(15)on each particle is equivalent to the classical N-body problem.The CPU requirements will be of the order O(N2)if no accelerating technique is utilized.Such a direct summation approach would quickly become prohibitive as the number of particles N increases.To address this problem,the Fast Multipole Method(FMM)29based on kernel evaluations is implemented in the current study.It does not require multipole expansions and works well for highly non-uniform particle distributions.Furthermore,it reduces the complexity from O(N2)to asymptotic O(N),and remarkably accelerates the computation.

Fig.4 Flowchart of the CFD/VTM coupling method.

The Particle Strength Exchange(PSE)scheme30is selected to solve the viscous diffusion term.One significant advantage of PSE lies in that this scheme provides enough flexibility in dealing with viscous effects and does not require any viscous splitting.In addition,vorticity conservation is naturally satisfied when two particles ‘exchange” their strengths.This conservation property is crucial to simulation accuracy because a non-physical variation of vorticity can seriously influence the predicted wake structure.The rationale of the PSE scheme is to solve an approximate integral equation for the Laplacian operator,which can be expressed as follows:

where ζεis the PSE kernel,which is identical to the Gaussian distribution in Eq.(13).The approximate manipulation of the Laplacian operator avoids any numerical differentiation,and the integral operator can be discretized by midpoint quadrature.

In the present study,a predictor-corrector time marching scheme is used to determine the new position xn+1and vorticity of the particles ωn+1,and can be written as follows:

3.4 术后应用中医护理技术的重要性 中医护理技术在手术患者的护理中有着不可忽视的重要性，对患者术后的康复起到重要的作用。对于手术患者而言，有效预防术后并发症是治疗成功的关键，可减轻病痛，缩短住院时间，降低住院费用。本研究表明，围手术期的中医护理技术干预，在泌尿外科微创手术患者中的应用发挥了较好的作用，得到患者及家属的接受，弘扬了祖国医学，发挥了中医的特色护理。

(1)Predictor step

(2)Corrector step

where subscript ‘vdt”denotes the ‘viscous diffusion term”.

2.3.Flowchart and coupling methodology

Fig.5 Time marching scheme for CFD/VTM coupling.

The flowchart of CFD/VTM is shown in Fig.4.After flow field initialization,the iterative procedure is repeated until converged results are obtained.An interpolation algorithm allowing data communication on the interface between CFD and the VTM is implemented.The key idea of this data exchange process is as follows:the VTM module provides outer velocity boundary conditions vinfor the CFD module,and receives the updated vortex sources from the CFD solution.The time marching scheme is given in Fig.5.The Lagrangian description enables the VTM to adopt a larger time-step size than that of CFD(i.e.,ΔtVTM＞ ΔtCFD).When conducting unsteady simulations,CFD corresponds to an azimuthal time step of 1°and the VTM has a 4°time step.

2.3.1.Transmission of vortex sources from CFD to VTM

In the current study,distributed airloads along the blade span are treated as lifting line vorticity source generators.This technique is referred to as ‘integrated vorticity source method”by He and Zhao.25Integrating the blade surface pressure calculated from CFD,we obtain the blade local lift vector Lb.Then the blade bound circulation can be determined by the Kutta-Joukowski theorem as

New viscous vortex particles are released into the flow field to account for blade bound circulation variations at each time step,and the vorticity strength of the new vortex particles is computed by Eq.(11).

References

Table 1 Caradonna-Tung rotor parameters.

Parameter Metric Airfoil NACA0012 Number of blade 2 Rotor radius(m) 1.143 Blade chord(m) 0.1905 Blade twist(°) 0

Fig.6 Half of the overset grid system used for full CFD simulation.

Fig.7 Four representative span wise locations of the Caradonna-Tung rotor blade.

Fig.8 Blade chord wise pressure distributions.

where Q is the conservative variable vector,Fc,Gc,and Hcare the convective fluxes,and Fv,Gv,and Hvare the viscous fluxes as follows:

The VTM specifies flow variables on the outer boundary of the CFD gird.31The velocity components u,v,w of the wake can be directly sent to CFD.One the other hand,since the VTM resolution is based on incompressible Navier-Stokes equations,fluid density and pressure cannot be specified like the velocity components.Considering that the flow is approximately is entropic for low-Mach number flow simulation,the density and pressure of the fluid are calculated through the following is entropic equations:

《金属焊接》涉及到焊接电源、焊接方法和设备、焊接材料、焊接冶金、焊接工艺、焊接应力与变形、焊接质量检验等内容，集电学、材料学、冶金学、力学等学科于一体。由于课程内容需要浓缩整个焊接领域中的基本知识点和各知识结构的精华部分，因此，课程内容主要突出体现了以下几个特色：

Finally,the conservative variables are updated by

本文从不同角度对广西与东盟国家跨境贸易发展及跨境人民币结算现状进行深入分析。同时，选取当年累计出口额、当年累计进口额及跨境人民币结算量等具有代表性的指标，运用协整检验、构建VAR模型和脉冲响应函数等分析方法，对广西跨境人民币结算ln（crmb）与广西—东盟进出口贸易ln（ix）、ln（ex）三个指标变量的关系进行实证分析。

3.Numerical results and discussions

3.1.Code validation and wake kinematics

Fig.9 RMS values of blade chordwise pressure distributions.

In the following subsections,the proposed CFD/VTM solver is validated against hover cases of the Caradonna-Tung model rotor,an OH-13E full-scale helicopter rotor,and a scaled V-22 tiltrotor.Simulation results are compared with available test data to show the capability of the presented hybrid solver.

3.1.1.Caradonna-Tung model rotor

The main parameters of this two-bladed model rotor are listed in Table 1.This rotor has no blade twist and the chord length c is constant.The rotor is operated at a collective pitch setting θ0of 8°with two different blade tip Mach numbers Matip,0.612 and 0.877.The near-body C-O type grid for the Eulerian domain consists of 225×43×74 grid points,where 225 points are used in the wrap around direction,43 points in the normal direction,and 74 points in the spanwise direction.The spacing between the outer boundary of the C-O grid and the blade surface is 3c.The number of vortex particles in the Lagrangian domain is about 52000.

In order to compare with the CFD/VTM coupling method,a full CFD simulation employing an overset grid system is also performed.Half of the overset grid system is shown in Fig.6.Blade chordwise pressure distributions at four representative spanwise locations(see Fig.7,r/R=0.68,0.80,0.89,and 0.96)are illustrated in Fig.8.The dots are from Caradonna and Tung’s experiment.32It can be seen that the results of both CFD/VTM(solid lines)and full CFD(dashed lines)methods are in good agreement with the experimental data,only with minor discrepancies at the blade leading edge.A Root Mean Square(RMS)analysis is carried out to evaluate the accuracy of these two methods,and the RMS is computed by

(3)In respect of vertical descending flight with a fixed collective pitch,the rotor thrust time history shows that the thrust starts to fluctuate at a moderate descent rate.A sudden loss of the thrust manifests the emergence of a deep VRS.Non-dimensional thrust behaviors of different rotor con figurations are similar.

Fig.10 presents the wake structure predicted by CFD/VTM.As shown in the 3D wake visualization(Fig.10(a)),concentrated tip vortices are preserved for a relatively long time.Both the near wake and the far wake could be clearly identified.In Fig.10(b)and(c),the rotor wake immediately contracts below the rotor disk plane,and demonstrates an ordered helicoidal structure in the first few revolutions.As the wake convects downstream further,the tip vortex filaments begin to bundle with each other.This mutual interaction phenomenon,or the so-called ‘vortex paring instability”,causes a breakdown of the vortex structure in the far wake.In contrast,only near-wake characteristics could be observed in the wake structure simulated by the full CFD method because of the well-known numerical dissipation (see Fig.11).The Lagrangian-based VTM model demonstrates its advantage on vortex tracking in the far wake.

Fig.10 Wake structure for the Caradonna-Tung rotor in hover(the CFD/VTM method).

Fig.11 Wake structure for the Caradonna-Tung rotor in hover(the full CFD method).

Fig.12 shows the tip vortex evolution in the blade near body domain at three selected planes downstream from the trailing edge(spacing=0.2c,0.5c,1.0c).On Slice 1,both the CFD/VTM and full CFD methods capture the tip vortex very well,while on Slices 2 and 3,CFD/VTM preserves the vorticity better than the full CFD method.Since the VTM captures the far wake structure better than CFD,it provides a more accurate feedback to the blade near body domain,and hence results in less diffused tip vortices.

明成祖永乐四年(1406)七月，明廷于交趾(即今越南)发生了历时一年的战争，朝廷派朱能、沐晟、张辅统率官军分道平叛。刘清从征，三次立功。因“从征交趾功，升陕西都指挥使”。［2］

The simulation time of these two methods are documented in Table 2.Computations are carried out on Intel eight-core processors with a 3.4 GHz CPU speed.CFD/VTM and full CFD simulations are performed with the same blade bodyfitted C-O grid.The number of background grid nodes is about 5.14 million in full CFD simulation.As shown,it requires 28.7 h while CFD/VTM only needs 5.6 h.A noticeable reduction of the computation time is achieved by the CFD/VTM solver,which is mainly because that the wake structure is simulated by particles.A dense background gird33is no longer neededin this hybrid solver.In addition,a larger step size that the VTM allows is also conducive to improve the computational efficiency.

Table 2 Simulation time of the CFD/VTM coupling method and the full CFD method.

Method CFD/VTM Full CFD Blade grid node 2×0.72 million 2×0.72 million Background flow field 52000 particles 5.14 million grid nodes Simulation time(h) 5.6 28.7

3.1.2.OH-13E helicopter rotor

To further validate the CFD/VTM hybrid solver,the downwash of a full-scale OH-13E rotor is calculated.The main characteristics of this rotor are given in Table 3.The rotor is operated at a tip Mach number of 0.403 with two different collective pitches,6.25°and 10.75°.Fig.13 shows the radial distributions of the down wash at four sampling stations below the rotor disk(Z=-0.1R,-0.3R,-0.5R,-0.7R),and videnotes the velocity of the down wash.As shown,the calculated data agrees well with the measured data.34The wake contracted near the radial position at r/R=0.8,and the downwash distributions are characterized by sharp velocity peaks at this location.Comparisons between the axial and radial components of tip vortex positions are given in Fig.14.As shown,correlations between the predicted and measured data are favorable.

Table 3 OH-13E rotor parameters.

Parameter Metric Airfoil NACA0015 Number of blade 2 Rotor radius(m) 5.353 Blade root chord(m) 0.3557 Blade tip chord(m) 0.2576 Blade twist(°) -4.25

Fig.12 Vorticity contours at selected planes downstream from the trailing edge of the blade tip.

Fig.13 Radial distributions of the rotor downwash at different sampling stations.

3.1.3.Scaled V22 ‘Osprey”tiltrotor

Fig.14 Tip vortex position for the OH-13E rotor in hover.

Table 4 0.658-scale V22 rotor parameters.

Parameter Metric Airfoil XN series Number of blade 3 Rotor radius(m) 3.81 Blade root chord(m) 0.615 Blade tip chord(m) 0.4 Blade twist(°) -47.5

Fig.15 Figure of merit for the scaled V22 tiltrotor in hover.

Anothervalidation caseisa 0.658-scaleisolated V-22‘Osprey” tiltrotor.The rotor parameters are tabulated in Table 4.This scaled tiltrotor has a very high blade twist angle of-47.5°and with a precone angle of 2.5°.More information about this rotor con figuration can be seen in Ref.35Fig.15 compares the figure of merit when the tiltrotor is operated at various hover conditions,where CTand σ represent the rotor thrust and solidity,respectively.It can be clearly seen that CFD/VTM-predicted results match well with the test data.36A typical wake structure for this three-bladed model tiltrotor is shown in Fig.16.

3.2.Descending flight simulation

Numerical simulations of the Caradonna-Tung rotor and the scaled V-22 tiltrotor in descending flight are performed using the developed CFD/VTM hybrid solver.In the following simulations,the vortex ring state is approached from a hovering state with a fixed collective pitch.The simulation is initialized at the hovering state,which isfollowed by computations at gradually increased vertical descent rates(μz=0.2，0.4，···，1.0，1.2).Initial flow field data at each descent rate is obtained from the previous one.The descent rate is normalized by the uniform induced velocity in hover,which is calculated by the formulaIn present cases,the value of vi0is 9.925 m/s for the Caradonna-Tung rotor and 14.62 m/s for the scaled V-22 tiltrotor.

Fig.16 Wake structure for the scaled V22 tiltrotor in hover.

Fig.17 Rotor thrust time histories at different descent rates with a fixed collective pitch.

Fig.18 Mean rotor induced velocity at different descent rates.

Fig.19 Difference between hover and VRS.

From Fig.17,it is seen that the non-dimensional thrust time histories of the two rotor con figurations are similar.However,the tiltrotor con figuration is less susceptible to the VRS with a higher vi0.At a low descent rate(μz＜ 0.4),the rotor thrust is not significantly in fluenced,and it is slightly higher than that in hover.As the descent rate increases,the rotor thrust starts to fluctuate,but a nearly converged solution can still be obtained(an incipient VRS).When the descent rate is higher than 1.0,a remarkable thrust loss occurs,which indicates the emergence of a deep VRS.

Fig.18 shows the normalized mean in flow of the rotor disk.As shown,the classical momentum theory is no more valid in descending flight.Predicted data matches well with the flight test of Taghizad et al.,9the wind tunnel test of Betzina,5and the empirical relation given by Castles and Gray.12Rotorinduced flow increases as the descent rate increases,and it grows faster than the opposing descent rate.Fig.19 shows the difference between hover and the VRS at a specific blade section.Compared with the hovering state,the high induced flow in the VRS dramatically reduces the effective Angle of Attack(AOA),and causes sudden loss of the rotor thrust.

魏金枝：《回忆“五四”时期的学生生活》，新民晚报副刊部著《夜光杯文粹(1946～1966)》，上海远东出版社1999年8月第1版，第466页。

其中：l1+l5=450 mm，且选定 l1=300 mm，l5=150 mm，根据整个装置的长度和一般家庭使用的座便器的尺寸选定l2一个合适的值为700 mm.

Fig.20 shows the snapshots of particle distributions at different descent rates around the rotor disk.As shown,the vortex particles tend to move toward the disk plane,and hence the convection of vortices is inhibited.The accumulated vortices interact with rotor blades,which leads to a fluctuation of the rotor thrust.Furthermore,they can dramatically increase the induced velocity(as shown in Fig.18).At a descent rate of μz=1.2,most particles are transported above the rotor disk plane.Fig.21 shows the flow regime at this fully developed VRS.The wake structure is characterized by a toroidal‘vortex ring”.This‘vortex ring” lingers within the vicinity of the rotor disk(Fig.21(a)),and induces a large-scale recirculating flow(Fig.21(b)),which corresponds to the schematic of air flow shown in Fig.1.

Fig.20 Snapshots of particle distributions at different descent rates(the Caradonna-Tung rotor).

Fig.21 Flow regime in a deep vortex ring state,μz=1.2(the Caradonna-Tung rotor).

4.Conclusions

An Eulerian-Lagrangian hybrid solver that combines a vorticity transport model with a CFD model has been developed in this paper.From numerical cases of hovering and vertical descending flight,the following conclusions can be drawn:

(1)In the hovering state,the blade surface pressure distribution,tip vortex position,induced velocity,and figure of merit predicted by the presented hybrid solver are in close agreement with experimental data,which indicates that it provides a high- fidelity analytical tool for rotorcraft application.

(2)The CFD/VTM coupling method demonstrates its exclusive capability in far wake modeling and vortex tracking.Rotor wake instabilities(i.e.,bundling of vortex filaments and wake breakdown in the far field)in hover that are hard to be captured by conventional CFD methods,are successfully simulated in the current study.In addition,since the far wake is modeled by viscous vortex particles,the computational cost of the hybrid solver is not expensive.The present coupling method will be added to the Chinese Laboratory of Rotorcraft Navier-Stokes(CLORNS)code in the near future.37

where ΔCpis the difference between the predicted data and the measured data in the experiment at each pressure transducer position.Results of the RMS analysis are shown in Fig.9.It is demonstrated that the predicted Cpdistributions of CFD/VTM are slightly better than those of the full CFD method.

转身帮泥巴扣好了带子。左小龙开着摩托车载着她走。当时是春天。是春天的中旬，是一个独立的气候。阳光洒满，云朵从云朵里穿透过来，空中的风就像是裙子撩动的气流，左小龙默默的载着泥巴到了一个垃圾站前。他把泥巴放下车，摘掉自己的头盔，再取下泥巴的头盔，问：“你是不是言情小说看多了？”

(4)When the descent rate increases,the wake vorticity tends to move upward to the rotor disk plane.In contrast with the helicoidal wake structure in hover,the deep VRS is characterized by a toroidal vortex ring.This ring induces a largescale recirculating flow around the rotor disk,which dramatically reduces the blade effective angle of attack,and eventually causes the sudden loss of the rotor thrust.

Acknowledgements

This work was co-supported by the Funding of Jiangsu Innovation Program for Graduate Education of China(No.KYLX16_0389)and the Fundamental Research Funds for the Central Universities of China.

建筑设施的屋顶设计是设计的重点内容，屋顶是建筑设施的第五立面，屋顶的形式具有多元化的特性，在不同的居住设施中，屋顶形式种类繁多，但是大部分是平屋顶的建筑形式。屋顶在实际的使用中，会出现漏水的问题，通过对坡屋顶和平屋顶渗漏率进行统计，坡屋顶的渗透率约为10%，平屋顶的渗透率约为70%，所以在改造旧建筑设施的过程中，为了处理好屋顶渗漏问题，选择将平屋顶改为坡屋顶不啻是一个好的解决方案。

由图10可知，截割头静止情况下，视觉测量的抬升角跳动在0.3°以内，回转角跳动在0.25°以内，满足可视化辅助截割系统设计要求。

1.Brand A,Dreier M,Kisor R,Wood T.The nature of vortex ring state.J Am Helicopter Soc 2011;56(2):1–14.

2.Varnes DJ,Duren RW,Wood ER.An onboard warning system to prevent hazardous ‘vortex ring state” encounters.Proceedings of the 26th European rotorcraft forum;2000 Sep 26–29;The Hague,The Netherlands.2000.

3.Landgrebe AJ.The wake geometry of a hovering helicopter rotor and its in fluence on rotor performance.J Am Helicopter Soc 1972;17(4):3–15.

4.Brotherhood P.Flow through a helicopter rotor in vertical descent.London:Her Majesty’s Stationery Office;1949.Report No.:N2735.

5.Betzina MD.Tiltrotor descent aerodynamics:a small scale experimental investigation of vortex ring state.Proceedings of the 57th American Helicopter Society annual forum;2001 May 9–11;Washington,D.C.,USA.Alexandria:The AHS International,Inc.;2001.

6.Washizu K,Azuma A,Koo J,Oka T.Experiments on a model helicopter rotor operating in the vortex ring state.J Aircr 1966;3(3):225–30.

The stretching term can be computed by substituting Eqs.(12)and(15)into ω ·∇u,which gives rise to

8.Green RB,Gillies EA,Brown RE.The flow field around a rotor in axial descent.J Fluid Mech 2005;534:237–61.

9.Taghizad A,Jimenez J,Binet L,HeuzéD.Experimental and theoretical investigations to develop a model of rotor aerodynamics adapted to steep descents.Proceedings of the 58th American Helicopter Society annual forum;2002 Jun 11–13;Montréal,Canada.Alexandria:The AHS International,Inc.；2002.

10.Kisor R,Blyth R,Brand A,MacDonald T.V–22 low-speed/high rate of descent(HROD)test results.Proceedings of the 60th American Helicopter Society annual forum;2004 Jun 7-10;Baltimore,USA.Alexandria:The AHS International,Inc.；2004.

11.Brand A,Kisor R,Blyth R,Mason D,Host D.V–22 high rate of descent(HROD)test procedures and long record analysis.Proceedings of the 60th American Helicopter Society annual forum;2004 Jun 7–10;Baltimore,USA.Alexandria:The AHS International,Inc.;2004.

12.Castles W,Gray R.Empirical relation between induced velocity,thrust,and rate of descent of a helicopter rotor as determined by wind tunnel test on four model rotors.Washington,D.C.:NACA;1951.Report No.:NACA TN 2474.

13.Wang SC.Analytical approach to the induced flow of a helicopter rotor in vertical descent.J Am Helicopter Soc 1990;35(1):92–8.

14.Basset PM,Prasad JVR.Study of the vortex ring state using bifurcation theory.Proceedings of the 58th American Helicopter Society annual forum;2002 Jun 11–13;Montréal,Canada.Alexandria:The AHS International,Inc.;2002.

15.Basset P,Chen C,Prasad JVR,Kolb S.Prediction of vortex ring state boundary of a helicopter in descending flight by simulation.J Am Helicopter Soc 2008;53(2):139–51.

16.Hoinville E,Renaud T.CFD simulation of helicopter rotor in the vortex ring state regime.Proceedings of the 63th American Helicopter Society annual forum;2007 May 1–3;Virginia Beach,US.Alexandria:The AHS International,Inc.;2007.

17.Cao YH,Yu ZQ,Su Y,Kang K.Combined free wake/CFD methodology for prediction transonic rotor flow in hover.Chin J Aeronaut 2002;15(2):65–71.

18.Shi YJ,Zhao QJ,Fan F,Xu GH.A new single-blade based hybrid CFD method for hovering and forward- flight rotor computation.Chin J Aeronaut 2011;24(2):127–35.

19.Tan JF,Wang HW.Panel/full-span free-wake coupled method for unsteady aerodynamics of helicopter rotor blade.Chin J Aeronaut 2013;26(3):535–43.

20.Stock M,Gharakhani A,Stone CP.Modeling rotor wakes with a hybrid OVERFLOW-vortex method on a GPU cluster.Proceedings of the 63th AIAA applied aerodynamics conference;2010 Jun 28–Jul 1.Chicago,USA.Reston:AIAA;2010.

21.Sugiura M,Tanabe Y,Sugawara H.Development of a hybrid method of CFD and prescribed wake model for helicopter BVI noise prediction.Proceedings of the 69th American Helicopter Society annual forum;2013 May 21–23;Phoenix,USA.Alexandria:The AHS International,Inc.;2013.

22.Leishman JG,Bhagwat MJ,Ananthan S.Free-vortex wake predictions of the vortex ring state for single-rotor and multirotor con figurations.Proceedings of the 58th American Helicopter Society annual forum;2002 Jun 11–13;Montréal,Canada.Alexandria:The AHS International,Inc.;2002.

23.Brown RE,Line AJ.efficient high-resolution wake modeling using the vorticity transport equation.AIAA J 2005;43(7):1434–43.

24.White house GR.Investigation of hybrid grid–based computational fluid dynamics methods for rotorcraft flow analysis.J Am Helicopter Soc 2011;56(3):1–10.

25.He CJ,Zhao JG.Modeling rotor wake dynamics with viscous vortex particle method.AIAA J 2009;47(4):902–15.

26.Conlisk AT.Modern helicopter rotor aerodynamics.Prog Aerosp Sci 2001;37(5):419–76.

27.Roe PL.Approximate Riemann solvers,parameter vectors,and difference schemes.J Comput Phys 1981;43(2):357–72.

28.Yokota R,Narumi T,Sakamaki R,Yasuoka K.Fast multipole methods on a cluster of GPUs for the meshless simulation of turbulence.Comput Phys Commun 2009;180(11):2066–78.

29.Ying LX,Biros G,Zorin D.A kernel-independent adaptive fast multipole algorithm in two and three dimensions.J Comput Phys 2004;196(2):591–626.

30.Eldredge JD,Leonard A,Colonius T.A general deterministic treatment of derivatives in particle methods.J Comput Phys 2002;180(2):686–709.

31.Shi YJ,Xu Y,Xu GH,Wei P.A coupling VWM/CFD/CSD method for rotor airload prediction.Chin J Aeronaut 2017;30(1):204–15.

32.Caradonna FX,Tung C.Experimental and analytical studies of a model helicopter rotor in hover.Vertica 1980;5(2):149–61.

33.Ye L,Zhang Y,Yang S,Zhu XL,Dong J.Numerical simulation of aerodynamic interaction for a tilt rotor aircraft in helicopter mode.Chin J Aeronaut 2016;29(4):843–54.

34.Boat wright DW.Measurements of velocity components in the wake of a full-scale helicopter rotor in hover.Virginia:U.S.Army Air Mobility Research and Development Laboratory;1972.Report No.:TR 72–33.

35.Acree CW.JVX proprotor performance calculations and comparisons with hover and airplane-mode test data.Washington,D.C.:NASA;2009.Report No.:NASA/TM–2009-215380.

36.Felker FF,Signor DB,Young LA,Betzina MD.Performance and loads data from a hover test of a 0.658-scale V-22 rotor and wing.Washington,D.C.:NASA;1987.Report No.:NASA-TM-89419.

37.Zhao QJ,Zhao GQ,Wang B,Wang Q,Shi YJ,Xu GH.Robust Navier-Stokes method for predicting unsteady flow field and aerodynamic characteristics of helicopter rotor.Chin J of Aeronaut 2018;31(2):214–24.