Integrable cases of a rigid body dynamics and integrable systems on the ellipsoids

Infinite-dimensional sets of integrable cases are found for the equations of a rigid body rotation around a fixed point in an axially symmetric potential field and also in more complicated fields in the presence of some symmetry of the rigid body inertia tensor.     

Influence of viscosity on the diffraction of sound by a circular aperture in a plane screen

The linearized equations of viscous fluid flow are used to analyze the diffraction of a time-harmonic acoustic plane wave by a circular aperture in a rigid plane screen. Arbitrary aperture size and arbitrary angle of incidence are considered. Sets of dual integral equations are derived for the diffracted velocity and pressure fields, and are solved by analytic reduction to sets of linear algebraic equations. In the case of normal incidence, numerical results are presented for the fluid velocity in the aperture and the power absorption due to viscous dissipation. The theoretical results for power absorption are compared to previously obtained results from high amplitude acoustic experiments in air. The conditions under which the dissipation predicted by linear theory becomes significantare quantified in terms of the fluid viscosity and sound speed, the acoustic frequency, and the aperture radius.     

Free vibration of a single-walled carbon nanotube containing a fluid flow using the Timoshenko beam model

The flexural vibration of the fluid-conveying single-walled carbon nanotube (SWCNT) is derived by the Timoshenko beam model, including rotary inertia and transverse shear deformation. The effects of the flow velocity and the aspect ratio of length to diameter on the vibration frequency and mode shape of the SWCNT are analyzed. Results show that the effects of rotary inertia and transverse shear deformation result in a reduction of the vibration frequencies, especially for higher modes of vibration and short nanotubes. The frequency is also compared with the previous study based on Euler beam model. In addition, if the ratio of length to diameter increased to 60, the influence of the shear deformation and rotary inertia on the mode shape and the resonant frequencies can be neglected. However, the influence is very obvious when the ratio decreased to 20. As the flow velocity of the fluid increases in the vicinity of 2π, the SWCNT reveals the divergence instability. It regains stability when the flow velocity reaches about 9. As the velocity increases further, the SWCNT undergoes a coupled-mode flutter and results in a larger amplitude.     

PIV method for research of the structure of pulsating flow in a smooth duct

The system for Particle Image Velocimetry (PIV) measurements in the flow with superimposed pulsations of the fluid (air) has been developed. Measurements of velocity and vorticity fields in a smooth duct in certain phases of superimposed pulsations have been performed. Statistics of a turbulent pulsating flow: velocity profiles, turbulent pulsations, and Reynolds stresses has been obtained.     

New Equation for Predicting Pipe Friction Coefficients Using the Statistical Based Entropy Concepts

In general, this new equation is significant for designing and operating a pipeline to predict flow discharge. In order to predict the flow discharge, accurate determination of the flow loss due to pipe friction is very important. However, existing pipe friction coefficient equations have difficulties in obtaining key variables or those only applicable to pipes with specific conditions. Thus, this study develops a new equation for predicting pipe friction coefficients using statistically based entropy concepts, which are currently being used in various fields. The parameters in the proposed equation can be easily obtained and are easy to estimate. Existing formulas for calculating pipe friction coefficient requires the friction head loss and Reynolds number. Unlike existing formulas, the proposed equation only requires pipe specifications, entropy value and average velocity. The developed equation can predict the friction coefficient by using the well-known entropy, the mean velocity and the pipe specifications. The comparison results with the Nikuradse’s experimental data show that the R2 and RMSE values were 0.998 and 0.000366 in smooth pipe, and 0.979 to 0.994 or 0.000399 to 0.000436 in rough pipe, and the discrepancy ratio analysis results show that the accuracy of both results in smooth and rough pipes is very close to zero. The proposed equation will enable the easier estimation of flow rates.     

Independent particle Schrödinger fluid: Moments of inertia

The philosophy of the single particle Schrödinger fluid, especially as regards the velocity fields which find such a natural role therein, is applied to the study of the moments of inertia of independent fermion system. It is shown that three simplified systems exhibit the rigid-body rotational velocity field in the limit of large A, and that the leading deviations, both on the average and fluctuating, from this large-A limit can be described analytically, and verified numerically. For a single particle in a Hill-Wheeler box the moments are studied numerically, and their large fluctuations identified with the specific energy level degeneracies of its parallelepiped shape. The full assemblage of these new and old results is addressed to the question of the necessary and sufficient condition that the moment has the rigid value. Counter examples are utilized to reject some conditions, and the conjecture is argued that unconstrained shape equilibrium might be the necessary and sufficient condition. The spheroidal square-well problem is identified as a promising test case.     

Measurement of temperature and velocity fields in a heater unit by liquid crystal thermometry and particle image velocimetry

Temperature and velocity fields in a heating unit for automobiles are measured through a model experiment in water tunnel using flow visualizations and image analysis to investigate the mixing mechanism of the flow that has passed through the heater and not passed through it inside the unit. The temperature fields are measured by the liquid crystal visualization technique combined with the field calibration using full color spline fitting technique, and the velocity fields are evaluated by a particle imaging technique with a cross-correlation algorithm. These results indicate an enhanced flow mixing at larger mix-door angles, which results in a shorter mixing distance of temperature and velocity downstream of the mix door. The enhanced flow mixing is caused by the high velocity fluctuations created by the flow separation over the mix door.     

Momentengleichungen für Plasmen mit anisotropem Neutralgasdruck und anisotropem Ionendruck

Strongly non-Maxwellian and non-isotropic velocity distributions of the neutral atoms and of the ions occur in collisionless plasmas at high degrees of ionization, especially in gas discharges at low pressures and high current densities and in high temperature plasmas. The velocity distributions and the related velocity moments for the neutral gas and the ion gas are calculated. The influence of the magnetic fields on the ions is neglected. Especially, the pressure tensors and the heat flow tensors are investigated. The differential equations are given for the velocity moments of the velocity distribution. Additional terms occur in the equation of motion, if the pressure is non-isotropic and non-Cartesian coordinates are used. A heat flow tensor is evaluated that closes the system of differential equations for the neutral gas consistently and allows to rederive typical formulas of the molecular neutral gas flow. The heat flow tensor essentially determines the type of the differential equation system for the velocity moments. It is shown, that the neutral gas temperature is not constant across the plasma. Different statements deal with the heat flow tensor in the ion gas. In particular, non-vanishing ion temperature on the axis and a system of differential equations for the positive column under free-fall conditions are investigated. The inertia terms for the ion gas and the neutral gas must be taken into account in the pressure balance of the plasma.     

A numerical method for fully resolved simulation (FRS) of rigid particle–flow interactions in complex flows

A fictitious-domain based formulation for fully resolved simulations of arbitrary shaped, freely moving rigid particles in unsteady flows is presented. The entire fluid–particle domain is assumed to be an incompressible, but variable density, fluid. The numerical method is based on a finite-volume approach on a co-located, Cartesian grid together with a fractional step method for variable density, low-Mach number flows. The flow inside the fluid region is constrained to be divergence-free for an incompressible fluid, whereas the flow inside the particle domain is constrained to undergo rigid body motion. In this approach, the rigid body motion constraint is imposed by avoiding the explicit calculation of distributed Lagrange multipliers and is based upon the formulation developed by Patankar [N. Patankar, A formulation for fast computations of rigid particulate flows, Center for Turbulence Research Annual Research Briefs 2001 (2001) 185–196]. The rigidity constraint is imposed and the rigid body motion (translation and rotational velocity fields) is obtained directly in the context of a two-stage fractional step scheme. The numerical approach is applied to both imposed particle motion and fluid–particle interaction problems involving freely moving particles. Grid and time-step convergence studies are performed to evaluate the accuracy of the approach. Finally, simulation of rigid particles in a decaying isotropic turbulent flow is performed to study the feasibility of simulations of particle-laden turbulent flows.     

On acoustic radiation by a rigid object in a fluid flow

A problem of sound radiation by an absolutely rigid object, moving with respect to the surrounding fluid, is considered on the basis of the Lighthill's equation for aerodynamic sound. An integral representation of the radiated acoustic field is utilized, where the field is characterized as the sum of three fields, generated by a volume distribution of monopoles and by distributions of monopoles and dipoles on the surface of the rigid object. It is shown that, due to a discontinuity of Lighthill's stress tensor on the rigid boundary, a layer of surface divergence of hydrodynamic stresses on the boundary must be taken into account when evaluating the volume integral over Lighthill's quadrupole sources. When the contribution of the surface divergence is included in the solution of Lighthill's equation, amplitudes of the monopole and dipole sound radiated by the rigid object are shown to depend on the potential components of the normal velocity and the pressure on the rigid surface. The obtained solution is compared with Curle's solution for this problem, which establishes that the sound radiation by a rigid object is determined by the force exerted by the object upon the fluid. Both solutions are applied to two known problems of sound scattering and radiation by a rigid sphere in variable pressure and velocity fields. It is shown that predictions based on the obtained solution are equivalent to the results known from literature, whereas Curle's solution gives predictions contradicting the known results. It is also shown that the Ffowcs Williams and Hawkings equation, which coincides with Curle's equation for an immoveable rigid object, does not lead to the correct predictions as well.     

Stereoscopic PIV measurements of flow in and around axial flow fan

To analyze the complex three-dimensional flow structure of an axial flow fan and determine the validity of its application, PIV is used to provide detailed space and time resolved experimental data for understanding and control of flow field. The high resolution stereoscopic PIV system was successfully employed in this study for the investigation of flow structure around the axial flow fan. Using the once-per-revolution signal from the rotor, image fields were captured at a fixed position of the blades and hence provides the ability to do phase-averaging. The three-dimensional instantaneous velocity fields, phase-averaged velocity fields, instantaneous and mean vorticity distributions of the stereoscopic PIV measurement results were represented at typical planes of the flow field. Phaseaveraged velocity fields were calculated based on 200 frames of the instantaneous stereoscopic PIV measurement results. From the velocity distribution, the vorticity and turbulent intensity distribution, which are known to be major factors of fan noise, were calculated and its diffusion was discussed as they travel downstream. From the reconstructed three-dimensional velocity iso-surface at 8 cross planes of the outlet flow fields, the three-dimensional features can be seen clearly.     

A boundary condition capturing immersed interface method for 3D rigid objects in a flow

To simulate the flow around an object, we can replace the object with the fluid enclosed by a singular force. We can then simulate the flow on a fixed domain with a fluid–fluid interface supporting the singular force. In this paper, we present a boundary condition capturing approach to determine the singular force for a 3D rigid object. We apply a discontinuous body force to enforce the rigid motion of the fluid replacing the object and compute the singular force based on the kinematics of the object. Due to the singular force and the body force, the flow is not smooth across the interface. We solve the flow using the immersed interface method. Our boundary condition capturing immersed interface method is very efficient and stable, and its accuracy based on the infinity norm is near second order for the velocity and above first order for the pressure.     

Quantitative analysis of PC MRI velocity maps: pulsatile flow in cylindrical vessels.

The accuracy of MR phase contrast (PC) velocity mapping, and the subsequent derivation of wall shear stress (WSS) values, has been quantitatively assessed. Using a retrospectively gated PC gradient-echo technique, the temporal-spatial velocity fields were measured for pulsatile flow in a rigid cylindrical vessel. The experimental data were compared with values derived from the Womersley solution of the Navier-Stokes equations. For a sinusoidal waveform, the overall root-mean-square (rms) difference between the measured and analytical velocities corresponded to 13% of the peak fluid velocity. The WSS derived from the data displayed a 14% rms difference with the analytical model. As an example of a more complicated flow, a triangular saw-tooth waveform was deconstructed into its Fourier components. Velocity maps and the WSS were calculated by the superposition of the individual solutions, weighted by the Fourier series coefficient, for each harmonic. The velocity and experimentally derived WSS agreed with the analytical results (4% and 12% rms difference, respectively). Evaluation of the analytical models allowed an estimate of the inherent accuracy in the measurement of velocity maps and WSS values.     

气液两相流中的声速研究

气液两相流在工业各领域中广泛存在,而声速是描述其声学性质的一个重要参数。本文从流体的体积弹性模量的定义出发,推导了气液两相流中的声速随含气率的变化关系式,即混合流体的Wood声速公式,将其声速的部分计算结果和其他作者的实验数据进行了比较,吻合良好。并通过COMSOL有限元模拟软件得到不同气体分布下圆管谐振腔最低阶模式的共振频率,间接数值模拟研究了含气率对声速的影响。模拟结果与理论计算结果一致,当气液两相流中含气率较低时,声速随含气率的增大急剧减小。本研究结果为确定声速与气液两相流中的含气率间的关系提供了参考依据。     

平板热管蒸发段多孔介质内流动和传热分析

通过假定平板热管多孔芯内流体的压力分布,对流动和传热进行了分析,得出结论:边界对于速度分布和流量计算的影响比惯性影响大,对于液层的分布二者影响都很小;多孔介质底部温度分布均匀,热流密度方向基本垂直于边界.     

The method of images for regularized Stokeslets

The image system for the method of regularized Stokeslets is developed and implemented. The method uses smooth localized functions to approximate a delta distribution in the derivation of the fluid flow due to a concentrated force. In order to satisfy zero-flow boundary conditions at a plane wall, the method of images derived for a standard (singular) Stokeslet is extended to give exact cancellation of the regularized flow at the wall. As the regularization parameter vanishes, the expressions reduce to the known images for singular Stokeslets. The advantage of the regularized method is that it gives bounded velocity fields even for isolated forces or for distributions of forces along curves. These are useful in the simulation of ciliary beats, flagellar motion, and particle suspensions. The expression relating force and velocity can be inverted to find the forces that generate a given velocity boundary condition. The latter is exemplified by modeling a cilium as a filament moving in a three-dimensional flow. The cilium velocity at various times is constructed from known data and used to determine the force field along the filament. Those forces can then reproduce the flow everywhere. The validity of the method is evaluated by computing the drag on a sphere moving near a wall. Comparisons with known expressions for the drag show that the method gives accurate results for spheres even within a distance from the wall equal to the surface discretization size.     

传感器和Origin7.0软件在电脑控制韦伯福斯摆的实验研究中应用

研究用电脑控制韦伯福斯摆的实验中,用光电传感器、位移传感器和力传感器测量摆的振动周期和摆幅、角速度等及弹簧的劲度系数,由此来学习一种测量刚体转动惯量的方法。     

Coherent structures over riblets in turbulent boundary layer studied by combining time-resolved particle image velocimetry(TRPIV),proper orthogonal decomposition(POD),and finite-time Lyapunov exponent(FTLE)

Time-resolved particle image velocimetry(TRPIV) experiments are performed to investigate the coherent structure's performance of riblets in a turbulent boundary layer(TBL) at a friction Reynolds number of 185. To visualize the energetic large-scale coherent structures(CSs) over a smooth surface and riblets, the proper orthogonal decomposition(POD) and finite-time Lyapunov exponent(FTLE) are used to identify the CSs in the TBL. Spatial-temporal correlation is implemented to obtain the characters and transport properties of typical CSs in the FTLE fields. The results demonstrate that the generic flow structures, such as hairpin-like vortices, are also observed in the boundary layer flow over the riblets, consistent with its smooth counterpart. Low-order POD modes are more sensitive to the riblets in comparison with the high-order ones,and the wall-normal movement of the most energy-containing structures are suppressed over riblets. The spatial correlation analysis of the FTLE fields indicates that the evolution process of the hairpin vortex over riblets are inhibited. An apparent decrease of the convection velocity over riblets is noted, which is believed to reduce the ejection/sweep motions associated with high shear stress from the viscous sublayer. These reductions exhibit inhibition of momentum transfer among the structures near the wall in the TBL flows.     

High-frequency diffraction corrections to backscattering cross sections from a rough three-dimensional surface

Diffraction corrections to scalar wave fields at perfectly free and rigid rough surfaces were derived by two iterations of the corresponding integral equations. These diffraction corrections to the pressure or normal velocity (which, in the geometrical optics limit, are doubled at perfectly rigid and free surfaces, respectively) were obtained with an accuracy of approximately 1k(2), where k is the wave number of incidence radiation. Based on these corrections to the surface fields, the backscattering cross sections at normal incidence from the statistically rough Gaussian surfaces were derived. It was found that for the gentle roughness, diffraction results in effective "smoothing" of roughness for rigid and free surfaces and increasing of the backscattering cross sections, but for a rigid surface with steep roughness, the "fictitious" surface can be more rough than the real one, and the diffraction corrections become negative.     

刚体转动惯量测定方法改进

对测定刚体转动惯量的方法进行改进,避免匀减速运动时初始角速度的随意性,降低实验时阻力矩随角速度变化对实验结果的影响,减小实验结果误差。