Steady-state distribution of charged particles in the diffusion approximation

The steady-state distribution of charged particles in a weakly ionized plasma is examined for the case in which volume ionization, recombination, and diffusion in a space-charge field occur. A joint solution is obtained for the equation for charged-particle balance and the Poisson equations for the case of planar and cylindrical plasma configurations satisfying the Schottky condition at the boundaries of the region. A solution is also found for the case in which the ionization is localized in a spherically symmetric volume and in which the Schottky condition is satisfied at infinity. The condition for the existence of a steady-state solution is given and analyzed.     

Further Development of a Dynamic Intermittency Model For Wake-Induced Transition

The transition model studied in this paper uses an equation for free-stream intermittency and one for near-wall intermittency, combined with the SST turbulence model. The model was already assessed in earlier work for wake-induced transition in boundary layers in separated state and in attached state. For separated state transition, the model predicts the start and growth of transition very well, except for the impossibility to describe the breakdown of the roll-up vortices in a 2D RANS simulation. In attached state transition, the model has the tendency to generate a delayed start of transition and a too slow growth rate in the initial phase of transition. We demonstrate that a delayed start of transition is inherent for a RANS simulation. We propose a repair for the too slow growth rate in the initial stage of the transition. Additionally, both for wake-induced transition in separated state and in attached state, the original model has the tendency to predict a state too near to laminar in between wakes at the trailing edge. We demonstrate that this is caused by the too rapid destruction of the near-wall intermittency during relaminarization and in the lack of activation of near-wall intermittency for low turbulence level in the free stream. We propose a repair for both deficiencies. The fourth improvement is that we take into account the influence of the free stream turbulence length scale in the criteria for onset of transition. We demonstrate the high quality of the improved model for wake-induced transition on a steam turbine stator blade for outlet Reynolds number 600,000 and two levels of turbulence in the background flow: 3% and 0.4%.     

The re-examination of the weakly nonlinear theory of hydrodynamic stability

The weakly nonlinear theory has been widely applied in the problem of hydrodynamicstability and also in other fields.However,although its application has been successful forsome problems,yet,for other problems,the results obtained are not satisfactory,especiallyfor problems like transition or the evolution of the vortex in the free shear flow,for whichthe goal of the theoretical investigation is not seeking for a steady state,but predicting anevolutional process.In this paper,we shall examine the reason for the unsuccessfulness andsuggest ways for its amendment.     

2017年度力学科学处基金项目受理情况

介绍了2017年度国家自然科学基金委员会数理科学部力学科学处受理的面上项目、青年科学基金项目、地区科学基金项目、重点项目、优秀青年科学基金、国家杰出青年科学基金和海外及港澳学者合作研究基金项目的申请情况及与2016年度的基金项目申请情况的比较.     

The singular behaviour in the vicinity of intersection between the body and free surface

The singular behaviour in the vicinity of intersection between the body and free surface is presented.It is shown that in the linear regime the singularity of velocity potential for transient problem is in d~2|nd.The singular behaviour for harmonic problem is the same as the result for the transient problem.In particular,the singularity for the harmonic problem with infinite frequency is in d~2 lnd for velocity potential(d is the distance between field point and intersection).     

A new evaluation method for Nusselt numbers in naphthalene sublimation experiments in rotating-disk systems

A comparative analysis of theoretical and empirical relations, as well as experimental data for mass transfer of a rotating disk in laminar, transitional and turbulent flows for naphthalene sublimation in air was done. New correlations between local and average Sherwood numbers for the entire disk were offered. A new evaluation approach for Nusselt numbers based on the experimental data for naphthalene sublimation in laminar, transitional and turbulent flows was developed.     

2015年度力学科学处基金项目受理情况介绍

介绍了2015 年度国家自然科学基金委员会数理科学部力学科学处受理的面上项目、青年科学基金项目、地区科学基金项目、重点项目、优秀青年科学基金、国家杰出青年科学基金和海外及港澳学者合作研究基金项目的申请情况.      

The use of displacement potentials in second order elasticity

This paper is concerned with the use of a representation in terms of displacement potentials in second order elasticity for equilibrium problems of homogeneous and isotropic materials. After justifying the adoption of an existing representation for linear elasticity for the purpose at hand, appropriate representations for solutions of second order elasticity problems in terms of displacement potentials (for both compressible and incompressible materials) are discussed. The use of the representations in obtaining complete solutions for equilibrium boundary-value problems is then illustrated by application to two examples of plane strain problems of compressible materials.     

Weakly-singular,weak-form integral equations for cracks in three-dimensional anisotropic media

Singularity-reduced integral relations are developed for displacement discontinuities in three-dimensional, anisotropic linearly elastic media. An isolated displacement discontinuity is considered first, and a systematic procedure is followed to develop relations for the displacement and stress fields induced by the discontinuity. The singularity-reduced relation for the stress is particularly important since it is in a form which allows a weakly-singular, weak-form traction integral equation to be readily established. The integral relations obtained for a general displacement discontinuity are then specialized to an isolated crack and to dislocations; the relations for dislocations are introduced to emphasize their direct connection to corresponding results for cracks and to allow earlier independent findings for these two types of discontinuities to be put into proper context. Next, the singularity-reduced integral equations obtained for an isolated crack are extended to allow treatment of cracks in a finite domain, and a pair of weakly-singular, weak-form displacement and traction integral equations is established. These integral equations can be combined to obtain a final formulation which is in a symmetric form, and in this way they serve as the basis for a weakly-singular, symmetric Galerkin boundary element method suitable for analysis of cracks in anisotropic media.     

Numerical investigation of particle lateral migration in straight channel flows using a direct-forcing immersed boundary method

Inertia-induced cross-stream migration has been recently exploited for precise position of particles in confined channel flows. In this work, a three-dimensional finite volume based immersed boundary method has been developed to study the lateral migration and hydrodynamic self-assembly of neutrally-buoyant particles in pressure-driven flows. Simulation results show that, in 2D channel flows, the equilibrium position for a circular particle is closer to the centreline for larger particle Reynolds number due to the increasing flow rate, while in 3D square duct flow, the equilibrium position for a spherical particle is near a face centre and is closer to the wall for larger particle Reynolds number. Self-assembly of a pair of particles is observed in 3D square duct flows but not in 2D channel flows. Mechanisms for the self-assembly are discussed.     

Fluid flow in an arbitrary cascade on an axisymmetric stream surface in a variable thickness layer

A solution is given for the problem of flow past a cascade on an axisymmetric stream surface in a layer of variable thickness, which is a component part of the approximate solution of the three-dimensional problem for a three-dimensional cascade. Generalized analytic functions are used to obtain the integral equation for the potential function, which is solved via iteration method by reduction to a system of linear algebraic equations. An algorithm and a program for the Minsk-2 computer are formulated. The precision of the algorithm is evaluated and results are presented of the calculation of an example cascade.In the formulation of [1, 3] the problem of flow past a three-dimensional turbomachine cascade is reduced approximately to the joint solution of two-dimensional problems of the averaged axisymmetric flow and the flow on an axisymmetric stream surface in an elementary layer of variable thickness.In the following we solve the second problem for an arbitrary cascade with finite thickness rotating with constant angular velocity in ideal fluid flow: the solution is carried out on a Minsk-2 computer.Many studies have been devoted to this problem. A method for solving the direct problem for a cascade of flat plates in a hyperbolic layer was presented in [2]. Methods were developed in [1, 3] for constructing the flow for the case of a channel with variable thickness; these methods are approximately applicable for dense cascades but yield considerable error for small-load turbomachine cascades. The solution developed in [4], somewhat reminiscent of that of [2], is applicable for thin, slightly curved profiles in a layer with monotonically varying thickness. A solution has been given for a circular cascade for layers varying logarithmically [5] and linearly [6]. Approximate methods for slightly curved profiles in a monotonically varying layer with account for layer variability only in the discharge component were examined in [7–9]. A solution is given in [10] for an arbitrary layer by means of the relaxation method, which yields a roughly approximate flow pattern. The general solution of the problem by means of potential theory and the method of singularities presented in [11] is in error because of neglect of the crossflow through the skeletal line. The computer solution of [12] contains an unassessed error for the calculations in an arbitrary layer. The finite difference method is used in [13] to solve the differential equation of flow, which is illustrated by numerical examples for monotonie layers of axial turbomachines. The numerical solution of [13] is very complex.The solution presented below is found in the general formulation with respect to the geometric parameters of the cascade and the axisymmetric surface and also in terms of the layer thickness variation law.The numerical solution requires about 15 minutes of machine time on the Minsk-2 computer.     

The influence of constraints on the properties of acceleration waves in isotropic thermoelastic media

The properties of acceleration waves are investigated for situations in which the waves propagate in isotropic heat-conducting elastic media subject to arbitrary sets of constraints. Conditions under which waves may exist in the presence of constraints are investigated for classes of constraints broad enough to encompass all those encountered in practice. Attention is focussed on principal waves, and results are presented for the growth of the amplitudes of such waves first for fronts of arbitrary curvature, and subsequently by specialisation for plane, cylindrical and spherical waves travelling in material which has undergone one-dimensional plane deformation, cylindrically symmetric and spherically symmetric deformation, respectively.     

Stress-softening Effects in the Transverse Vibration of a Non-Gaussian Rubber String

The Mullins effect in the small amplitude transverse vibration of a rubber cord is investigated. The fundamental frequency is determined for a specific class of stress-softening materials. Analytical relations for the cord vibration frequency are illustrated graphically for three phenomenological models. These results demonstrate the role of the material parameters and exhibit response characteristic of those reported in experiments by others and subsequently described here in new experiments. Frequency versus stretch results for two kinds of non-Gaussian molecular network models for rubber elasticity are compared with experimental data for four varieties of rubber cords, for each of which only three experimentally determined material constants are needed. It is shown that the theoretical predictions stand in excellent agreement with test data.     

Stochastic equations for continuum and determination of hydraulic drag coefficients for smooth flat plate and smooth round tube with taking into account intensity and scale of turbulent flow

The stochastic equations of continuum are used for determining the hydraulic drag coefficients. As a result, the formulas for the hydraulic drag coefficients dependent on the turbulence intensity and scale instead of only on the Reynolds number are proposed for the classic flows of an incompressible fluid along a smooth flat plate and a round smooth tube. It is shown that the new expressions for the classical drag coefficients, which depend only on the Reynolds number, should be obtained from these new general formulas if to use the well-known experimental data for the initial turbulence. It is found that the limitations of classical empirical and semiempirical formulas for the hydraulic drag coefficients and their deviation from the experimental data depend on different parameters of initial fluctuations in the flow for different experiments in a wide range of Reynolds numbers. On the basis of these new dependencies, it is possible to explain that the differences between the experimental results for the fixed Reynolds number are caused by the difference in the values of flow fluctuations for each experiment instead of only due to the systematic error in the processing of experiments. Accordingly, the obtained general dependencies for the smooth flat plate and the smooth round tube can serve as the basis for clarifying the results of experiments and the experimental formulas, which used for continuum flows in different devices.     

2014年度力学科学处基金项目受理情况介绍

介绍了2014年度国家自然科学基金委员会数理科学部力学科学处受理的面上项目、青年科学基金项目、地区科学基金项目、重点项目、优秀青年科学基金、国家杰出青年科学基金和海外及港澳学者合作研究基金项目的申请情况.      

Higher-order and length-scale statistics of velocity and temperature fluctuations in turbulent boundary layer along a heated vertical flat plate

Time-developing direct numerical simulation (DNS) was performed to clarify the higher-order turbulent behaviors in the thermally-driven boundary layers both in air and water along a heated vertical flat plate. The predicted statistics of the heat transfer rates and the higher-order turbulent behaviors such as skewness factors, flatness factors and spatial correlation coefficients of the velocity and temperature fluctuations in the natural-convection boundary layer correspond well with those obtained from experiments for space-developing flows. The numerical results reveal that the turbulent structures of the buoyancy-driven boundary layers are mainly controlled by the fluid motions in the outer region of the boundary layer, and these large-scale structures are strongly connected with the generation of turbulence in the thermally-driven boundary layers, in accordance with the actual observations for space-developing flows. Moreover, to specify the turbulence structures of the boundary layers, the cross-correlation coefficients and the characteristic length scales are examined for the velocity and thermal fields. Consequently, it is found that with a slight increase in freestream velocity, the cross-correlation coefficient for the Reynolds shear stress and turbulent heat flux increases for opposing flow and decreases for aiding flow, and the integral scales for the velocity and temperature fields become larger for opposing flow and smaller for aiding flow compared with those for the pure natural-convection boundary layer.     

Boiling heat transfer in a vertical tube with freon 114

Test results for boiling heat transfer coefficients for R114 in vertical stainless tubes are reported both for upflow and downflow direction. Results are compared with different formulas given in the literature. A Recommendation of the numerical value ofC _sf in Rohsenow's formula for fully developed nucleation boiling for R114/stainless steel combination is given.     

Rate form of the Eshelby and Hill tensors

Expressions are derived for the rates of change of the S and P tensors for transformed homogeneous inclusions in an anisotropic comparison medium undergoing prescribed changes of its elastic moduli. General results are obtained for ellipsoids and then reduced to yield explicit expressions in terms of the Stroh eigenvalues for cylindrical and disk-shaped inclusions in anisotropic solids and for spherical inclusions in isotropic solids. Applications are illustrated by solving the rate problem for an inhomogeneity in a large volume of a comparison medium, which is shown to be readily adaptable to standard averaging techniques for predictions of rates of change of overall moduli of composite materials experiencing evolution of phase moduli.