人体上呼吸道内气流运动特性的数值模拟分析

运用计算流体动力学(CFD)方法对人体上呼吸道内的气流运动特性进行了数值模拟,通过PIV实验对数值模拟结果进行了验证。研究结果表明:气流在咽部外壁、气管外壁发生分离现象,气流在气管内壁形成局部高速区,支气管内的气流在分叉处发生分离,靠近支气管内壁速度较高,并且在支气管边界层的外缘速度达到最大值。气管和支气管内的二次涡流运动和轴向速度的分布使得气管支气管内壁受到的剪应力较大,内壁粘膜更容易受到损伤。     

On the recovery of the reference trajectory of rotation

The boundary value problem of determining the time dependence of the spacecraft angular velocity vector from measurements performed by uniaxial transducers rigidly fixed in the body-fixed frame is stated and solved. The proposed algorithm ensures that the mean and variance of the variations in the estimates of the angular velocity vector are smaller than the corresponding accuracy characteristics of uniaxial transducers. The implementation of the proposed technique permits increasing the accuracy of inertial measurement devices.     

Wing/body kinematics measurement and force and moment analyses of the takeoff flight of fruitflies

In the paper, we present a detailed analysis of the takeoff mechanics of fruitflies which perform voluntary takeoff flights. Wing and body kinematics of the insects during takeoff were measured using high-speed video techniques.Based on the measured data, inertia force acting on the insect was computed and aerodynamic force and moment of the wings were calculated by the method of computational fluid dynamics. Subtracting the aerodynamic force and the weight from the inertia force gave the leg force. The following has been shown. In its voluntary takeoff, a fruitfly jumps during the first wingbeat and becomes airborne at the end of the first wingbeat. When it is in the air, the fly has a relatively large "initial" pitch-up rotational velocity(more than5 000°/s) resulting from the jumping, but in about 5 wingbeats, the pitch-up rotation is stopped and the fly goes into a quasi-hovering flight. The fly mainly uses the force of jumping legs to lift itself into the air(the force from the flapping wings during the jumping is only about 5%–10% of the leg force). The main role played by the flapping wings in the takeoff is to produce a pitch-down moment to nullify the large "initial" pitch-up rotational velocity(otherwise, the fly would have kept pitching-up and quickly fallen down).     

Announcements

Recent numerical investigations on pressure surges during pump trip in pumping installations showed that by including an air entrainment variable wave speed model, reasonable predictions of transient pressure surges with proper phasing and attenuation of pressure peaks can be obtained. These calculated results are consistent with similar field measurements made with the pumps operating at low pump cut-out levels, when air entrainment due to an attached surface vortex was observed. However, in the numerical calculation procedures it is assumed that the inertia of the moving elements of the check valve is small and that the check valve closes at zero reverse flow velocity. In practice, check valves seldom close precisely at zero reverse flow velocity. With the check valves not closing at zero reverse velocity, the present numerical computations show that the air content in a fluid system can adversely affect the check valve performance. With the fluid system operating within a critical range of air entrainment values, the present analysis showed that there is a possibility of ‘check valve slamming’ when the check valves are selected based only on the analysis of an air-free system. This phenomenon is confirmed through field observations.     

Generalization of the mathematical model of lungs for describing the intensity of the tracheal sounds during forced expiration

The possibility to relate the nature of the forced expiration tracheal sounds with the sound radiation by a separated flow that arises in the region of dynamic trachea constriction during forced expiration is investigated. A mathematical model of forced expiration is used for estimating. The calculated form of the time dependence of the sound intensity during forced expiration qualitatively corresponds to the experimental dependence obtained experimentally for normal subjects. The results should be taken into account in the physical explanation of tracheal sound generation mechanisms and in the justification of using the tracheal sound characteristics in the diagnostics of human lung pathologies.     

Effect of near-wall turbulence enhancement on the mechanisms of particle deposition

The modification of deposition mechanisms of small particles in wall turbulence due to enhanced near-wall fluctuations is presented. The direct numerical simulation database of turbulent air flow over a water surface populated by gravity-capillary waves of small wave slope was used to mimic the enhancement in fluctuation intensity. Lagrangian tracking of particles is performed under the assumption of one-way coupling between the particles and the flow. Two sets of particles have been considered with inertial response times of 5 and 15, respectively, normalized using the friction velocity at the air–water interface and the kinematic viscosity of air. Compared to wall-bounded flow, the particle deposition rates on the interface were found to be considerably higher; specifically for the low-inertia particles, an eightfold increase was observed. The deposition rate for particles of higher inertia increased by only 60%. The correlation characterizing particle deposition rates for wall-bounded flows, where the deposition rate is proportional to the square of the particle response time, was found to be invalid for the flow with enhanced near-wall turbulence. Comparison with experimental results on particle deposition onto rough walls showed better correlation. Depositing particles were divided into free-flight and diffusional deposition populations. Since the primary effect of the interfacial waves is to increase the turbulence intensity in the near-interface region with high particle concentration, a remarkable increase in diffusional deposition is observed. As in wall-bounded flows, diffusional deposition is seen to be the dominant mechanism of deposition. The free-flight mechanism, where particles acquire velocities high enough to travel directly to the interface, remains unaffected by enhanced near-wall velocity fluctuations.     

Leading-Order Dynamics of Gravity-Driven Flows in?Free-Standing Soap Films

We derive the leading-order equations that govern the dynamics of the flow in a falling, free-standing soap film. Starting with the incompressible Navier?CStokes equations, we carry out an asymptotic analysis using parameters that correspond to a common experimental setup. We account for the effects of inertia, surface elasticity, pressure, viscous stresses, gravity, and air drag. We find that the dynamics of the flow is dominated by the effects of inertia, surface elasticity, gravity, and air drag. We solve the leading-order equations to compute the steady-state profiles of velocity, thickness, and pressure in an experiment in which the film is in the Marangoni elasticity regime. The computational results, which include a Marangoni shock, are in good accord with the experimental measurements.     

Dynamic behaviour of an axially moving plate undergoing small cylindrical deformation submerged in axially flowing ideal fluid

The out-of-plane dynamic response of a moving plate, travelling between two rollers at a constant velocity, is studied, taking into account the mutual interaction between the vibrating plate and the surrounding, axially flowing ideal fluid. Transverse displacement of the plate (assumed cylindrical) is described by an integro-differential equation that includes a local inertia term, Coriolis and centrifugal forces, the aerodynamic reaction of the external medium, the vertical projection of membrane tension, the bending resistance, and external perturbation forces. In the two-dimensional model thus set up, the aerodynamic reaction is found analytically as a functional of the cylindrical displacement, using the techniques of complex analysis. The resulting integro-differential problem is discretized in space with the Fourier-Galerkin method, and integrated in time with the diagonalization method. Examples are computed with physical parameters corresponding to air and some paper materials. The effects of the surrounding fluid on the critical velocity and first natural frequency are investigated, for stationary air, for an air mass moving with the plate, and for some arbitrary axial fluid velocities. The obtained results are applicable for both an ideal membrane and a plate with nonzero bending rigidity.     

Regression Modeling of the Acousto-Biomechanical Characteristics of Wheezes during Human Forced Expiration

Within the framework of a specified acousto-biomechanical model, a possibility of explaining wheezes during forced expiration on the basis of a vortex separation mechanism is tested and the localization of the zones in the bronchial tree where forced expiratory wheezes are generated is refined. As an experimental model, a group of 18 healthy volunteers from 18 to 44 years old (median is equal 19) was used. On the basis of linear regression modeling the relationship between the principal spectral frequency of medium-frequency forced expiratory wheezes (400–600 Hz) recorded on the trachea and the standard volume flow-rates of expired air measured in computer spirometry is analyzed. The data obtained show that vortex separation in the air flow at the bronchial tree bifurcations (where a stepwise increase in the cross-section area takes place) is a probable mechanism of medium-frequency forced expiratory wheeze generation and that during the forced expiratory maneuver the vortex separation zone tends to be displaced deeper into the bronchial tree.     

周期性结构热动力时间-空间多尺度分析

研究一种时间-空间多尺度渐近均匀化分析方法，模拟不同的极端热和动力载荷下微尺度多 相周期性结构中热动力响应问题，并建立一个广义的波动函数场控制方程描述热动力响应. 通过引入一个放大空间尺度和两个缩小时间尺度，在不同时间尺度上获得由空间非均匀性引 起的波动效应和非局部效应. 根据高阶均匀化理论在空间和时间上进行均匀化，获得高阶非 局部函数场波动方程. 并进一步用C0连续修正了高阶非局部函数场波动方程的有限元近 似解，使问题的求解避免了对有限元离散的C1连续性要求. 并与经典的空间均匀化方法 相比较，指出了经典的空间均匀化方法的局限性，进一步以一维非傅立叶热传导和热动力问 题为例，讨论了各种情况下方法的正确性与有效性.     

Dynamics of PIV seeding particles in turbulent premixed flames

Particle image velocimetry (PIV) estimates the fluid velocity field measuring the displacement of small dispersed particles between two successive instants separated by a small time interval. The accuracy of the measurements depends on the ability of the particles to accommodate their velocity to the fluid fluctuations. When the fluid is subjected to extreme accelerations, the small but finite inertia prevents the particles from following the fluid, originating a substantial relative velocity. This effect is shown to be crucial for applications of PIV to turbulent premixed combustion, particularly in the product region at locations just behind the instantaneous flame front. The issuing inaccuracy may easily spoil the estimate of certain statistical observables which are of crucial importance in the theory of turbulent premixed combustion. By exploiting the direct numerical simulation of a model air/methane flame, a suitable criterion for proper particle seeding is validated and compared with the corresponding experiments with a combined PIV/OH-LIF (laser-induced fluorescence) system. The proposed parameter, the flamelet Stokes number, depends on particle properties and thermochemical conditions of the flame and substantially restricts the particle dimensions required for a reliable estimate of the relevant flow statistics.     

Basic properties of natural vibrations of an extended segment of a pipeline

Natural transverse vibrations of an extended segment of a pipeline containing a uniformly moving fluid are considered. The mechanical model under study takes into account the inertial forces of the pipe and environment and the moment of Coriolis and centrifugal forces arising because of the medium motion. It is proved that all natural frequencies of the pipeline rigidly clamped at both ends are real (and hence no flutter can arise in this model). For the first three modes, the dependence of the eigenvalues on the fluid flow velocity (varying from zero to the buckling velocity) are constructed, and their properties depending on the inertia parameter are studied. Families of vibration mode shapes of the pipeline are constructed and investigated.     

Transient free‐surface flow of a viscoelastic fluid in a narrow channel

The interplay between inertia and elasticity is examined for transient free‐surface flow inside a narrow channel. The lubrication theory is extended for the flow of viscoelastic fluids of the Oldroyd‐B type (consisting of a Newtonian solvent and a polymeric solute). While the general formulation accounts for non‐linearities stemming from inertia effects in the momentum conservation equation, and the upper‐convected terms in the constitutive equation, only the front movement contributes to non‐linear coupling for a flow inside a straight channel. In this case, it is possible to implement a spectral representation in the depthwise direction for the velocity and stress. The evolution of the flow field is obtained locally, but the front movement is captured only in the mean sense. The influence of inertia, elasticity and viscosity ratio is examined for pressure‐induced flow. The front appears to progress monotonically with time. However, the velocity and stress exhibit typically a strong overshoot upon inception, accompanied by a plug‐flow behaviour in the channel core. The flow intensity eventually diminishes with time, tending asymptotically to Poiseuille conditions. For highly elastic liquids the front movement becomes oscillatory, experiencing strong deceleration periodically. A multiple‐scale solution is obtained for fluids with no inertia and small elasticity. Comparison with the exact (numerical) solution indicates a wide range of validity for the analytical result. Copyright © 2004 John Wiley & Sons, Ltd.     

A second strain gradient elasticity theory with second velocity gradient inertia – Part II: Dynamic behavior

This paper is the sequel of a companion Part I paper devoted to the constitutive equations and to the quasi-static behavior of a second strain gradient material model with second velocity gradient inertia. In the present Part II paper, a multi-cell homogenization procedure (developed in the Part I paper) is applied to a nonhomogeneous body modelled as a simple material cell system, in conjunction with the principle of virtual work (PVW) for inertial actions (i.e. momenta and inertia forces), which at the macro-scale level takes on the typical format as for a second velocity gradient inertia material model. The latter (macro-scale) PVW is used to determine the equilibrium equations relating the (ordinary, double and triple) generalized momenta to the inertia forces. As a consequence of the surface effects, the latter inertia forces include (ordinary) inertia body forces within the bulk material, as well as (ordinary and double) inertia surface tractions on the boundary layer and (ordinary) inertia line tractions on the edge line rod; they all depend on the acceleration in a nonstandard way, but the classical laws are recovered in the case of no higher order inertia. The classical linear and angular momentum theorems are extended to the present context of second velocity gradient inertia, showing that the extended theorems—used in conjunction with the Cauchy traction theorem—lead to the local force and moment (stress symmetry) motion equations, just like for a classical continuum. A gradient elasticity theory is proposed, whereby the dynamic evolution problem for assigned initial and boundary conditions is shown to admit a Hamilton-type variational principle; the uniqueness of the solution is also discussed. A few simple applications to wave propagation and dispersion problems are presented. The paper indicates the correct way to describe the inertia forces in the presence of higher order inertia; it extends and improves previous findings by the author [Polizzotto, C., 2012. A gradient elasticity theory for second-grade materials and higher order inertia. Int. J. Solids Struct. 49, 2121–2137]. Overall conclusions are drawn at the end of the paper.     

Oscillating flows in capillaries filled with immiscible liquids

Hydrodynamic flows generated by mechanical vibrations of a capillary filled with immiscible liquids are investigated. Air bubbles are contained at the hermetically sealed ends of the capillary. Equations for the change in the volumes of the air bubbles as functions of time and velocity distribution in the liquids are obtained for the case when the radius of the capillary is much less than the lengths of the liquid columns. Results of numerical calculations are given for a capillary filled with two liquids: water and mercury. Amplitude-frequency dependences of the change in volumes of the air bubbles are constructed which have a resonance nature. Graphs of the dependence of the velocity of the water and the mercury on the radial coordinate at different times are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 13–18, September–October, 1987.     

Measurement of shock waves velocity in the air plasma of capacitively coupled RF discharge

The velocity of shock wave propagation in the air plasma of stationary capacitively coupled RF discharge at different gas pressure and charged particles concentration has been measured. It is shown, that the velocity of the shock wave increases at the increase of the concentration. Measurement results are brought to the universal dependence. Received 17 August 1998 / Accepted 10 December 1999     

Flow rate through a long channel with a rotating damper at the end

The experimental data on the flow rate (velocity) of a fluctuating air flow are presented on a wide fluctuation frequency range at a constant pressure difference between the channel entry and exit. The superimposed flow fluctuations were produced by periodic cut-off of the exit section by a rotating damper. A considerable dependence of the mean flow rate (velocity) on the wave structure of the flow is established. A flow rate minimum corresponds to resonance flow modes with a maximum relative amplitude of the flow velocity fluctuations.     

Modeling of thermally activated dislocation glide and plastic flow through local obstacles

A unified phenomenological model is developed to study the dislocation glide through weak obstacles during the first stage of plastic deformation in metals. This model takes into account both the dynamical responses of dislocations during the flight process and thermal activations while dislocations are bound by obstacle arrays. The average thermal activation rate is estimated using an analytical model based on the generalized Friedel relations. Then, the average flight velocity after an activation event is obtained numerically by discrete dislocation dynamics (DD). To simulate the dynamical dislocation behavior, the inertia term is implemented into the equation of dislocation motion within the DD code. The results from the DD simulations, coupled with the analytical model, determine the total dislocation velocity as a function of the stress and temperatures. By choosing parameters typical of the face centered cubic metals, the model reproduces both obstacle control and drag control motion in low and high velocity regimes, respectively. As expected by other string models, dislocation overshoots of obstacles caused by the dislocation inertia at the collisions are enhanced as temperature goes down.     

Dynamics of a prestressed incompressible layered half-space under moving load

The paper addresses a plane problem: a concentrated force acts on a plate resting on an elastic half-space with homogeneous prestrain. The equations of motion of the plate incorporate shear and rotary inertia. The half-space is assumed to be incompressible and isotropic in the natural state. The elastic potential is given in general form and is only specified for numerical purposes. The dependence of the critical velocity of the load and the stress-strain state on the prestresses is analyzed for different ratios between the stiffnesses of the layer and half-space and different contact conditions. The calculations are carried out for a half-space with Bartenev-Khazanovich potential __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 3, pp. 36–54, March 2008.