Integrated numerical model for vegetated surface and saturated subsurface flow interaction

The construction of an integrated numerical model is presented in this paper to deal with the interactions between vegetated surface and saturated subsurface flows. A numerical model is built by integrating the previously developed quasi-three-dimensional （Q3D） vegetated surface flow model with a two-dimensional （2D） saturated groundwater flow model. The vegetated surface flow model is constructed by coupling the explicit finite volume solution of 2D shallow water equations （SWEs） with the implicit finite difference solution of Navier-Stokes equations （NSEs） for vertical velocity distribution. The subsurface model is based on the explicit finite volume solution of 2D saturated groundwater flow equations （SGFEs）. The ground and vegetated surface water interaction is achieved by introducing source-sink terms into the continuity equations. Two solutions are tightly coupled in a single code. The integrated model is applied to four test cases, and the results are satisfactory.

NUMERICAL SIMULATION OF BUBBLE FLOW INTERACTIONS

Bubble flow interaction can be important in many practical engineering applications. For instance, cavitation is a problem of interaction between nuclei and local pressure field variations including turbulent oscillations and large scale pressure variations. Various types of behaviours fundamentally depend on the relative sizes of the nuclei and the length scales of the pressure variations as well as the relative importance of bubble natural periods of oscillation and the characteristic time of the field pressure variations. Similarly, bubbles can significantly affect the performance of lifting devices or propulsors. We present here some fundamental numerical studies of bubble dynamics and deformation, then a practical method using a multi-bubble Surface Averaged Pressure （DF-Multi-SAP ） to simulate cavitation inception and scaling, and connect this with more precise 3-D simulations. This same method is then extended to the study of two-way coupling between a viscous compressible flow and a bubble population in the flow field.

Quasi-stationary planetary wave-mean flow interactions in the Northern Hemisphere stratosphere and their responses to ENSO events

Based on the ERA-40 reanalysis data from the European Centre for Medium-Range Weather Forecasts and the output of ECHAM5/MPI-OM, this study investigated the interactions between the quasi-stationary planetary wave （SPW） and mean flow, and their responses to E1 Nifio-Southern Oscillation （ENSO） events in the northern hemispheric stratosphere. Results show that the activity of SPW is the strongest in winter, when the SPW propagates along the polar waveguide into the stratosphere and along the low-latitude waveguide to the subtropical tropopause. The analysis of three dimensional SPW structure indicates that the main sources of SPW activity are located over the Eurasian continent and the North Pacific north of 45°N. On the one hand, the two waveguides of the SPW reflect the influence of mean flow on the propagation of the SPW. On the other hand, the upward propagating SPW can interact with the stratospheric mean flow, leading to deceleration of the zonal mean westerly. Furthermore, the SPW exhibits clear responses to ENSO events. During E1 Nifio winters, the SPW in the strat- osphere tends to propagate more upward and poleward. Its interactions with mean flow can induce a dipole pattern in zonal mean zonal winds, with accelerated westerly winds at low-middle latitudes and decelerated westerly winds at high latitudes. The ECHAM5/MPI-OM model reproduces the climatology of the SPW well. Although the simulated SPW is slightly weaker than the observations in the stratosphere, the model＇s performance has significant improvements compared with other GCMs used in previous studies. However, there are still some problems in the responses of the SPW to ENSO in the model. Although the model reproduces the responses of both the amplitude and the SPW-mean flow interactions to ENSO well in the troposphere, the stratospheric responses are quite weak. Therefore, further studies are needed to improve the simulation of the stratospheric atmospheric circulation and related dynamical processes.

The theory of interaction between wave and basic flow

This paper investigates the interaction between transient wave and non-stationary and non-conservative basic flow. An interaction equation is derived from the zonally symmetric and non-hydrostatic primitive equations in Cartesian coordinates by using the Momentum-Casimir method. In the derivation, it is assumed that the transient disturbances satisfy the linear perturbation equations and the basic states are non-conservative and slowly vary in time and space. The diabatic heating composed of basic-state heating and perturbation heating is also introduced. Since the theory of wave-flow interaction is constructed in non-hydrostatic and ageostrophic dynamical framework, it is applicable to diagnosing the interaction between the meso-scale convective system in front and the background flow. It follows from the local interaction equation that the local tendency of pseudomomentum wave-activity density depends on the combination of the perturbation flux divergence second-order in disturbance amplitude, the local change of basic-state pseudomomentum density, the basic-state flux divergence and the forcing effect of diabatic heating. Furthermore, the tendency of pseudomomentum wave-activity density is opposite to that of basic-state pseudomomentum density. The globally integrated basic-state pseudomomentum equation and wave-activity equation reveal that the global development of basic-state pseudomomentum is only dominated by the basic-state diabatic heating while it is the forcing effect of total diabatic heating from which the global evolution of pseudomomentum wave activity results. Therefore, the interaction between the transient wave and the non-stationary and non-conservative basic flow is realized in virtue of the basic-state diabatic heating.

Interannual variations in energy conversion and interaction between the mesoscale eddy field and mean flow in the Kuroshio south of Japan

Using 19-year satellite altimetric data, variations in the eddy kinetic energy, energy exchanges and interaction between the eddy fi eld and mean fl ow are discussed for the Kuroshio south of Japan. In the seasonal cycle, the eddy kinetic energy level is a minimum in December/January and a maximum in April/May. In addition to seasonal variations, the eddy kinetic energy undergoes interannual changes. The energy transfers mainly from the mean fl ow to the eddy fi eld in the Kuroshio south of Japan, and dominant energy exchanges mainly occur along the Kuroshio path south of Japan in each year from 1993 to 2011. In addition, there is often barotropic instability south of Honshu. Regarding interactions between the eddy fi eld and mean fl ow, cyclonic and anticyclonic accelerations are also found along the Kuroshio path and they fl ank each other. There is cyclonic acceleration always imposed on southeast of Kyushu, and anticyclonic acceleration dominates south of Honshu from 2001 to mid-2005. Reynolds stress is used to explain the dynamic process of energy exchange. Furthermore, lag-correlation and linear regression analysis show that variability of the energy conversion rate and Reynolds stress involve responses to eddy acceleration at two time scales. The enhanced eddy acceleration induces large Reynolds stress, and enhanced Reynolds stress or barotropic instability further enforces energy transfer from the mean fl ow to the eddy fi eld.

Determination of Slip Length in Couette Flow Based on an Analytical Simulation Incorporating Surface Interaction

An analytical simulation based on a new model incorporating surface interaction is conducted to study the slip phenomenon in the Couette flow at different scales. The velocity profile is calculated by taking account of the micro-force between molecules and macro-force from the viscous shearing effect, as they contribute to the achieve- ment of the slip length. The calculated results are compared with those obtained from the molecular dynamics simulation, showing an excellent agreement. Further, the effect of the shear rate on the slip is investigated. The results can well predict the fluid flow behaviors on a solid substrate, but has to be proved by experiment.

DIAMOND PORT JET INTERACTION WITH SUPERSONIC FLOW

Interaction flow field of the sonic air jet through diamond shaped orifices at different incidence angles （10 degrees, 27.5 degrees, 45 degrees and 90 degrees） and total pressures （0.10 MPa and 0. 46 MPa） with a Mach 5.0 freestream was studied experimentally. A 90 degrees circular injector was examined for comparison. Crosssection Mach number contours were acquired by a Pitot-cone five-hole pressure probe. The results indicate that the low Mach semicircular region close to the wall is the wake region. The boundary layer thinning is in the areas adjacent to the wake. For the detached case, the interaction shock extends further into the freestream, and the shock shape has more curvature, also the low-Mach upwash region is larger. The vortices of the plume and the height of the jet interaction shock increase with increasing incidence angle and jet pressure. 90 degrees diamond and circular injector have stronger plume vorticity, and for the circular injector low-Mach region is smaller than that for the diamond injector. Tapered ramp increases the plume vorticity, and the double ramp reduces the level of vorticity. The three-dimensional interaction shock shape was modeled from the surface shock shape, the center plane shock shape, and crosssectional shock shape. The shock total pressure was estimated with the normal component of the Mach number using normal shock theory. The shock induced total pressure losses decrease with decreasing jet incidence angle and injection pressure, where the largest losses are incurred by the 90 degrees, circular injector.

Air interaction around outdoor air-cooled condensers

In order to increase cooling or heating efficiency,a porous computational fluid dynamics（CFD）model is employed to predict the thermo-fluid status and optimize the placement of outdoor units.A full scale model is established to validate the accuracy of CFD simulation in terms of velocity and temperature distributions.The comparison between the measurement and the simulation shows a good agreement.By evaluating the condensers＇ sucked air temperature with CFD for three units installed in a row,it is found that the minimum separation distance among neighboring units is 0.2 m;a vertical wall should be apart from the unit line by at least 0.8 m;and large different operating pressures among units do not impact the flow rate and the heat transfer of the other units meaningfully.

Aerodynamic interaction between forewing and hindwing of a hovering dragonfly

The phase change between the forewing and hindwing is a distinct feature that sets dragonfly apart from other insects.In this paper,we investigated the aerodynamic effects of varying forewing-hindwing phase di ff erence with a60 inclined stroke plane during hovering flight.Force measurements on a pair of mechanical wing models showed that in-phase flight enhanced the forewing lift by 17%and the hindwing lift was reduced at most phase differences.The total lift of both wings was also reduced at most phase di ff erences and only increased at a phase range around in-phase.The results may explain the commonly observed behavior of the dragonfly where 0 is employed in acceleration.We further investigated the wing-wing interaction mechanism using the digital particle image velocimetry（PIV）system,and found that the forewing generated a downwash flow which is responsible for the lift reduction on the hindwing.On the other hand,an upwash flow resulted from the leading edge vortex of the hindwing helps to enhance lift on the forewing.The results suggest that the dragonflies alter the phase di ff erences to control timing of the occurrence of flow interactions to achieve certain aerodynamic effects.

An application of interacting shear flows theory: exact solution for unsteady oblique stagnation point flow

An analytical solution of the governing equations of the interacting shear flows for unsteady oblique stagnation point flow is obtained. It has the same form as that of the exact solution obtained from the complete NS equations and physical analysis and relevant discussions are then presented.

Effects of bubbly flow on bending moment acting on the shaft of a gas sparged vessel stirred by a Rushton turbine

The bending moment acting on the overhung shaft of a gas-sparged vessel stirred by a Rushton turbine,as one of the results of fluid and structure interactions in stirred vessels,was measured using a moment sensor equipped with digital telemetry.An analysis of the shaft bending moment amplitude shows that the amplitude distribution of the bending moment,which indicates the elasticity nature of shaft material against bending deformation,follows the Weibull distribution.The trends of amplitude mean,standard deviation and peak deviation characteristics manifest an "S" shape versus gas flow.The "S" trend of the relative mean bending moment over gas flow rate,depending on the flow regime in gas-liquid stirred vessels,resulted from the competition among the nonuniformity of bubbly flow around the impeller,the formation of gas cavities behind the blades,and the gas direct impact on the impeller when gas is introduced.A further analysis of the bending moment power spectral density shows that the rather low frequency and speed frequency are evident.The low-frequency contribution to bending moment fluctuation peaks in the complete dispersion regime.

Gao's interacting shear flows( ISF) theory and its inferences and their applications

Gao＇s viscous/in-viscid interacting shear flows （ISF） theory, proposed by professor Gao Zhi in Institute of Mechanics, China Academy of Science, and its inferences and their applications in computational fluid dynamics （CFD） are reviewed and some subjects worthy to be studied are pro- posed in this paper. The flow-field and motion law of ISF, mathematics definition of strong viscous shear layer flow in ISF, ISF equations, wall-surface compatibility criteria （Gao＇s criteria ）, space scale variety law of strong viscous shear layer reveals flow mechanism and local space small scale triggered by strong interaction that cause some abnormal severe local pneumatic heating phenomenon in hypersonic flow. Gao＇s ISF theory was used in near wall flow, free ISF flow simulation and design of computing grids, Gao＇s wall-surface criteria were used to verify calculation reliability and accuracy of near wall flows, ISF theory approximate analytical result of shock waves-boundary layer interac- tion and ISF equations were used to obtain the numerical exact solution of local area flow （ such as stationary point flow）. Some new subjects, such as, improving near-wall turbulent models according to the turbulent flow simulation satisfying the wall-criteria and illustrating relation between grid-con- vergence based on the wall criteria and other convergence tactics, are suggested. The necessity of applying Gao＇s ISF theory and wall criteria is revealed. Difficulties and importance of hypersonic vis- cous/in-viscid interaction phenomenon were also emphasized.

NON-LINEAR WATER WAVES ON SHEARING FLOWS

This article is devoted to the study of the propagations of the non- linear water waves on the shear flows. Assuming μ = kh is small and ε/μ~2 ～ 0 (1), and the base flow is uniformly sheared, the modified Boussinesq equation is obtained. We calculate propagations of the single sohtary wave with vorticity Γ = 0, >0 and <0. The influences of the vorticity are manifested. At the end examples of the interactions of two solitary waves, moving in opposite and the same directions, are given. Besides the phase shift, there also occur second wavelets after head-on collision.

An Interdecadal Change in the Influence of the NAO on Atlantic-Induced Arctic Daily Warming around the Mid-1980s

After approaching 0℃owing to an Atlantic storm at the end of 2015,the Arctic temperature approached freezing again in 2022,indicating that Arctic daily warming events remain a concern.The NCEP/NCAR Reanalysis dataset was used to investigate the influence of the NAO on the Arctic winter daily warming events induced by Atlantic storms,known as the Atlantic pattern-Arctic Rapid Tropospheric Daily Warming(Atlantic-RTDW)event.Atlantic-RTDW events are triggered by Atlantic storms that transport warm and humid air masses moving into the Arctic.Furthermore,an interdecadal change in the influence of NAO on Atlantic-RTDW-event frequency was observed around the mid-1980s.Specifically,before the mid-1980s(pre-transition period),500-hPa southerly(northerly)wind anomalies occupied the North Atlantic(NA)in the positive(negative)phase of NAO,which increased(decreased)the Atlantic-RTDW events occurrence by steering Atlantic storms into(away from)the Arctic;thus,the NAO could potentially influence the Atlantic-RTDW-event frequency.However,the relationship between the NAO and the Atlantic-RTDW-event frequency has weakened since the mid-1980s(post-transition period).In the post-transition period,such 500-hPa southerly(northerly)wind anomalies over the NA hardly existed in the positive(negative)phase of NAO,which was attributed to a stronger Atlantic Storm Track(AST)activity intensity than that in the pre-transition period.During this period,the strong AST induced an enhanced NAOrelated cyclone via transient eddy-mean flow interactions,resulting in the disappearance of southerly and northerly wind anomalies over the NA.

The Propagation and Transport Effect of Planetary Waves in the Northern Hemisphere Winter

Based on the transformed Eulerian-mean equations, the dynamics of planetary waves are discussed. Both observations and simulations indicate that in the Northern Hemisphere winter there are two waveguides for the meridional propagation of quasi-stationary planetary waves. One is the high latitude waveguide, and the other is the low latitude waveguide. These results are in good agreement with theoretical analysis. Moreover, the convergence of EP flux indicates that the stratospheric sudden warming is the result of anomalous planetary wave propagation along the high latitude waveguide and its interaction with mean flows. The tropical quasi-biennial oscillation (QBO) winds, which represent one significant variation of zonal flow in the lower stratosphere at low latitudes, can influence the low latitude waveguide of planetary wave propagation. Our results of the wave-mean flow coupled model show that these tropical winds can also modulate the high latitude waveguide significantly in the case of wave-mean flow interaction.The transport effect of planetary waves on ozone is also analyzed. The residual mean circulation forced by planetary waves indicates that there is strong transport circulation for the dissipative planetary waves. Under the forcing of northward eddy heat transport, a positive transport circulation can result which rises at low latitudes and sinks at high latitudes. At the same time, the modification of planetary wave propagation by the equatorial QBO winds is shown to have an important impact on the transport circulation. The model results indicate that the meridional transport is amplified during the easterly phase of the QBO. This mechanism may explain the interannual variability of ozone in the stratosphere at high latitudes.

Experimental and Analytical Investigation on Aerodynamic Characteristics of Helicopter Scissors Tail Rotor

An experimental and analytical investigation on the aerodynamic characteristics of a helicopter scissors tail rotor has been made. The experiments are carried out on a model rotor rig and the variations of the rotor thrust and torque with scissors angles are investigated. The effects of two different configurations, i.e., Configuration L in which the lower blade is in front of the upper blade and Configuration U in which the upper blade is in front of the lower blade, on rotor performance are compared. The experimental results have shown that the thrust for Configuration L is greater than that for Configuration U, and the strong blade-vortex interaction may occur at some scissors angles. A free-wake analytical model is then developed for predicting the aerodynamic characteristics of a scissors tail rotor and validated by numerical examples. By using the model, the blade induced velocity, lift distribution and tip vortex displacement of the scissors rotor are calculated. Based upon the calculations, the variations of the thrust with scissors angles in the experiments are analyzed and the strong blade-vortex interaction occurring in the experiments is explained. Finally, several conclusions are presented.

The Propagation of Wave Packets and Its Relationship with the Subtropical Jet over Southern China in January 2008

The propagation of wave packets and its relationship with the subtropical jet was investigated for the period 26 29 January 2008 over southern China using ECMWF Interim re-analysis data. Wave packets propagated from the north to the south side of an upper front with eastward development along the upper front during this period. Due to the eastward development of propagation, the acceleration of geostrophic westerly winds shifted eastward along the front. There were two primary sources of the propagation of wave packets at around 30°N. The first was the temperature inversion layer below 500 hPa, and the second was baroclinic zones located along the polarward flank of the subtropical jet in the middle and upper troposphere. Most wave packets propagated horizontally from the baroclinic zones and then converged on the zero meridional gradients of zonal winds.

ANALYSIS OF STABILITY CHARACTERISTICS OF NONLINEAR FLUID STRUCTURE INTERACTION SYSTEM IN INTERNAL FLOWS

A set of nonlinear partial differential equations was used to describe the motions of an internal flowing system, with the consideration of transient fluid structure interaction (FSI). The stability characteristics of the system were analyzed on the basis of the nonlinear FSI-equation model. The results show that the stability of the nonlinear FSI system become very complicated because of the influence of strongly coupled interaction between flowing liquid and vibrating pipe. The flowing velocity, pressure, the piping stiffness, and constrained conditions, etc. play a very important role in developing procedure of the system stabilization. In an analysis example given in this paper, for example, when the flowing pressure is comparatively small, the system exhibits flutter instability, and when the flowing pressure is increases to a specified value, a subcritical divergence and a supercritical flutter bifurcation take successively place in the same velocity range.

Analysis of Pseudomomentum Wave-Activity Density in a Heavy Rainfall Event in East China

Taking into account moisture in virtue of general potential temperature,the author derive a three-dimensional(3D) pseudomomentum wave-activity relation for the moist atmosphere from the primitive equations in Cartesian coordinates using the Momentum-Casimir method.Since the wave-activity relation is constructed in an ageostrophic and non-hydrostatic dynamical framework,it may be applicable to diagnosing the evolution and propagation of mesoscale systems leading to heavy rainfall.The theoretical analysis shows that,besides the local change of wave-activity flux divergence and source or sink,the wave-activity relation includes two additional forcing terms.The first is the zonal gradient of difference between general potential temperature and potential temperature perturbations,and the second is the covariance of the solenoid and gradient of water vapor,denoting the direct influence of moisture on wave-activity density.The wave-activity density was applied to a heavy precipitation event occurring in the Jianghuai region of China.The calculation showed that the wave-activity density was consistent with 6-h accumulated precipitation observations,in terms of both spatial distribution and temporal tendency.This suggested that the disturbance represented by wave-activity density was closely related to the heavy precipitation.Although the wave-activity flux divergence and the covariance of the solenoid and gradient of water vapor made the primary contribution to the local change of wave-activity density,the covariance was more remarkable.The zonal gradient of difference between general potential temperature and potential temperature perturbations made a weaker contribution to the waveactivity density.

Performance analysis and flow mechanism of channel wing considering propeller slipstream

The development of the Urban Air Mobility concept has proposed stringent requirements for the fixed-wing vehicle’s Short/Vertical Take-Off and Landing(S/VTOL)performance.With its lift-enhancing impact,the channel wing,can improve aircraft low-speed performance and Short Take-Off and Landing(STOL)capability.It is critical for the performance analysis and flow mechanism study of a channel wing that considers the influence of propeller slipstream to guide the design of S/VTOL aircraft.The law and mechanism for the effect of the enclosing angle of the arc wing,the phase angle of the propeller,the tip clearance,the rotational speed and the chordwise position of the propeller on the channel wing are explored utilizing the quasi-steady multi reference frame method.A channel wing with a larger enclosing angle has a better ability to enhance lift and reduce drag using propeller slipstream.The effect of propeller phase angle on channel wing aerodynamic forces is periodic and weak.Increasing propeller rotational speed is helpful to enhance lift and resist flow separation for channel wing.It can reduce drag for tractor configuration but increase drag for pusher configuration.However,the nose-down pitching moment of a channel wing will grow dramatically,making longitudinal trimming in aircraft layout design challenging.