中子输运方程DSN方法的稳定性、收敛性和误差估计

中子输运方程的离散s_N方法(以下简称DSN方法)有显式递推求解方便,程序逻辑简单,存储量节省而又能保证一定精确度等优点,应用甚多,是解一维和多维中子输运方程的一种有效方法. 197l年,P.Lascaux和P.A.Raviart在[1]中,用能量法证明了标准方程     

Unsteady behaviour of a heterogeneous elastic beam floating on shallow water

The unsteady behaviour of a thin elastic Euler beam with heterogeneous structural properties, floating freely on the surface of an ideal incompressible liquid is investigated using the linear theory. The unsteady behaviour of the beam is due to the incidence of a localized wave on its surface or initial deformation. Two methods of solving the problem are proposed in which the sagging of the beam is sought in the form of an expansion in eigenfunctions of the oscillations of a heterogeneous beam (the first method) or of a homogeneous beam (the second method) in the void. In both methods the problem is reduced to solving an infinite system of ordinary differential equations for the unknown amplitudes. The effect of different actions on a beam having a piecewise-constant distribution of the cylindrical stiffness and the specific mass is investigated. The eigenvalues of the systems of differential equations are determined.     

A New Exponential Compact Scheme for the Two-Dimensional Unsteady Nonlinear Burgers' and Navier-Stokes Equations in Polar Cylindrical Coordinates

In this article, a new compact difference scheme is proposed in exponential form to solve two-dimensional unsteady nonlinear Burgers' and Navier-Stokes equations of motion in polar cylindrical coordinates by using half-step discretization. At each time level by using only nine grid points in space, the proposed scheme gives accuracy of order four in space and two in time. The method is directly applicable to the equations having singularities at boundary points. Stability analysis is explained in detail and many benchmark problems like Burgers', Navier-Stokes and Taylor-vortex problems in polar cylindrical coordinates are solved to verify the accuracy and efficiency of the scheme.     

The motion of a cylindrical body in a stratified fluid under the action of a radiation force

The free motion of a thin cylindrical body is investigated based on a previously derived expression for the radiation force acting on moving point sources in a stratified fluid. The fundamental equations of motion are derived, the limits of applicability of the approximation used are indicated and the results of calculations of typical trajectories of a body which begins to move with a specified velocity from a position of neutral buoyancy at an angle to the horizon are presented. Calculations of the trajectory of motion of a thin cylindrical body in a stratified fluid when the total radiation force is taken into account show that the effect of the lateral component of this force is considerable and leads not only to quantitative corrections but also to qualitative effects (for example, to an increase in the oscillations of the body and a change in its direction of motion). The results obtained pertain both to the motion of solids in fluids and to the translational motion of vortex dipoles in weakly stratified media.     

三重介质轴对称二维不定常渗流的精确解

本文求得了三重介质轴对称二维不定常渗流的线源解和有限封闭地层的面源解.它们不仅概括和发展了已有的均质和双重介质的主要结果,而且给出多重介质弹性渗流的基本特征.     

Space-charge effects on current transport in media having cylindrical and spherical symmetries

Space-charge effects on transport in solids under a given appliedpotential are examined for cylindrical and spherical symmetries.It has been proven for these symmetries that a predominant spacecharge of the same sign as that of the mobile species producesa retardation of the current, whereas a dominant space chargeof opposite sign to that of the mobile species enhances thecurrent. Analytical expressions are developed which can be usedto evaluate quantitatively the effects of the retardation orthe enhancement for these geometries with arbitrary space-chargedistributions.     

A solution method for the sub-surface stresses and local deflection of a semi-infinite inhomogeneous elastic medium

This paper proposes analytical Fourier series solutions (based on the Airy stress function) for the local deflection and subsurface stress field of a two-dimensional graded elastic solid loaded by a pre-determined pressure distribution. We present a selection of numerical results for a simple sinusoidal pressure which indicates how the inhomogeneity of the solid affects its behaviour. The model is then adapted and used to derive an iterative algorithm which may be used to solve for the contact half width and pressure induced from contact with a rigid punch. Finally, the contact of a rigid cylindrical stamp is studied and our results compared to those predicted by Hertzian theory. It is found that solids with a slowly varying elastic modulus produce results in good agreement with those of Hertz whilst more quickly varying elastic moduli which correspond to solids that become stiffer below the surface give rise to larger maximum pressures and stresses whilst the contact pressure is found to act over a smaller area.     

Harmonic decomposition analysis of contact mechanics of bonded layered elastic solids

The paper presents an iterative method for obtaining footprint, pressure distribution, local deformation and sub-surface stress field for the contact between a rigid cylindrical indenter and an elastic flat substrate. The methodology is applicable for semi-infinite, as well as for thin or thick bonded elastic layered solids with high or low elastic moduli. All findings are in accord with the observed behaviour of hard wear resistant and soft solid lubricating coatings. It is shown that the decomposed contact pressure distribution into a series of harmonic waves induces sub-surface stress fields that decay into the depth of the solid according to their wavelengths. Consequently, conditions vis-à-vis fatigue spalling and adhesion performance may be predicted for given thickness of layered bonded elastic solids.     

Exact solutions of problems of the three-dimensional flow of ideal viscous incompressible fluids near cylindrical surfaces

Three-dimensional flows of an incompressible fluid, the parameters of which depend on two coordinates and time, are considered. The stream surfaces of such flows are cylindrical. The equations of continuity and the Navier-Stokes equations can be transformed to relations, one of which is the equation for the stream function the other is the integral of the equations relating the pressure and the stream function, and the third is a linear equation for the projection of the velocity vector onto the axis parallel to the generatrix of the cylindrical surfaces. The problems of modelling the flows are considered on the basis of the exact solutions of the Navier-Stokes equations and Euler's equations using examples. Relations for the distribution of the flow parameters in the channel created by hyperbolical cylinders are derived for the case of unsteady inviscid flow. The streamlines of these flows are situated on the side surfaces of the hyperbolical cylinders and intercept the generatrices of the cylinders at certain indirect angles. The flow around a circular cylinder and the flow of fluid inside an elliptic cylinder are considered in the case of steady inviscid flow. The streamlines on the circular cylinder are arranged transverse to the cylinder (the projection of the velocity vector onto the coordinate axis, parallel to the generatrix of the cylinder, is equal to zero). Far from the cylinder the streamlines are also situated on a cylindrical surfaces, but not transverse to the cylinder, making certain indirect angles with the generatrix. Viscous three-dimensional flows, possessing a certain symmetry, are considered. In the case of radial symmetry the streamlines are helical lines. The non-planar Couette flow between parallel moving planes is characterized by the fact that the velocity vectors, being situated in the same plane, are collinear, while the velocity vectors in parallel planes are not collinear. Relations for viscous steady three-dimensional flows, using well-known relations, obtained for the stream function of two-dimensional flows, are given.     

The plane contact problem of steady thermoelasticity taking heat generation into account

The plane steady contact problem of thermoelasticity when there is heat generation from friction, which arises when an infinite cylindrical punch moves over the surface of an elastic half-space along its generatrix, is considered. It is assumed that heat exchange between the free boundary of the half-space and the surrounding medium obeys Newton's law, while the condition for ideal thermal contact exists in the region in which the solids interact. The problem is reduced to a system of three integral equations in the heat fluxes and temperature. The effect of the thermal and mechanical properties of the cylinder and the half-space on the main contact characteristics is investigated numerically.     

Transient flow in a rotating porous annulus due to suction at the walls

A highly porous material occupies the annular region between two coaxial infinitely long cylinders. A viscous incompressible fluid fills this porous medium and is initially in a state of rigid rotation together with the medium. The flow has been disturbed by imposing suction/injection at the outer/inner cylindrical boundaries respectively. The Brinkman's law has been used to represent the fluid motion. The exact solution for the resulting unsteady flow is obtained by Laplace transformation technique. The transient evolution of the boundary layers and the response of steady boundary layers to the resistance of the medium are discussed.     

Unlimited cumulation under one-dimensional unsteady compression of an ideal gas

Various methods of unlimited cumulation (UC) of an ideal (inviscid and non-heat-conducting) gas subject to one-dimensional unsteady compression by a plane, cylindrical, or spherical piston are considered. The most perfect method, namely, UC with isentropic compression from rest to rest, which is referred to as “ideal” (IUC), is compared with three other methods of UC, which correspond to well-known self-similar solutions of one-dimensional gas compression. The most effective of these is UC with a reflected shock wave, behind which the compressed gas is homogeneous and at rest, as in IUC. The efficiency of various methods of UC is estimated by the ratio of the work done during compression to the work in the case of IUC, the ratio of the internal energy to the total energy of the compressed gas, and the degree of gas homogeneity with respect to the Lagrangian variable. Computations of these characteristics are carried out for a perfect gas with various adiabatic exponents.     

Elastodynamics in parabolic cylinders

Explicit solutions are constructed of an elasto-dynamical boundary value problem for cylindrical solids.

---

Zusammenfassung Es werden explizite Lösungen eines elastodynamischen Randwertproblems für zylindrische Festkörper konstruiert.

     

Viscous flow and temperature distribution in rotating annulus with suction and injection

Part A.—Steady viscous incompressible flow in a rotating coaxial cylindrical annulus with suction and injection is studied. The unsteady flow is also considered. Part B.—An exact solution for temperature distribution at different constant wall temperatures is obtained. It is assumed that the rate of injection at the inner wall equals the rate of suction at the outer wall.     

Dynamic strength around a coated nanowire with surfaces/interfaces in a half solid

In this paper, the multiple scattering of anti-plane shear waves around a coated nanowire with surfaces/interfaces embedded in a half solid is studied, and the dynamic stress at the two surfaces/interfaces is presented. The boundary condition at the edge of the half solid is satisfied by the image method. The analytical solutions of displacements in the two half solids, in the coating layers, and inside the nanowire are expressed by wave function expansion method. The expanded mode coefficients are determined by satisfying the boundary conditions at the two surfaces/interfaces of the coated nanowire and the straight edge of the structure. The addition theorem for cylindrical wave function is employed to accomplish the superposition of displacement fields in the two half solids. Analyses show that the properties of the outer and inner interfaces show different effect on the dynamic stress around the nanowire. The dynamic stress distribution around the nanowire is also significantly related to the interfacial properties at the edge and the position of the nanowire.     

On a Control Problem for a Linear System with Measurements of a Part of Phase Coordinates

We consider a three-dimensional unsteady flow with a rotating detonation wave arising in an annular gap of an axially symmetric engine between two parallel planes perpendicular to its symmetry axis. The corresponding problem is formulated and studied. It is assumed that there is a reservoir with quiescent homogeneous propane–air combustible mixture with given stagnation parameters; the mixture flows from the reservoir into the annular gap through its external cylindrical surface toward the symmetry axis, and the parameters of the mixture are determined by the pressure in the reservoir and the static pressure in the gap. The detonation products flow out from the gap into a space bounded on one side by an impermeable wall that is an extension of a side of the gap. Through a hole on the other side of the gap and through a conical output section with a half-opening angle of 45°, the gas flows out from the engine into the external space. We formulate a model of detonation initiation by energy supply in which the direction of rotation of the detonation wave is defined by the position of the energy-release zone of the initiator with respect to the solid wall situated in a plane passing through the symmetry axis. After a while, this solid wall disappears (burns out). We obtain and analyze unsteady shock-wave structures that arise during the formation of a steady rotating detonation. The analysis is carried out within single-stage combustion kinetics by the numerical method based on the Godunov scheme with the use of an original software system developed for multiparameter calculations and visualization of flows. The calculations were carried out on the Lomonosov supercomputer at Moscow State University.     

绕振动机翼非定常气动力迟滞特性的模拟研究

利用数值模拟方法开展了机翼在不同条件下,由于俯仰振动引起的非定常气动力迟滞特性的模拟计算研究.根据有限体积方法对非定常欧拉方程进行数值求解,以确定相应问题的流场及气动力特性;同时,分别以具有NACA-0012翼型的矩形机翼和带有65°后掠角的三角翼为例,研究了机翼绕不同转轴或以不同频率振动的非定常气动力迟滞特性.     

Two-dimensional ultrasound Doppler velocimeter for velocity field measurements of liquid metal flows

We present a novel fluid flow measurement system based on the pulsed-wave ultrasound Doppler velocimetry being able to determine two-dimensional velocity fields. It applies for the measurement of unsteady liquid metal flows driven by electromagnetic forces concerning the research field of magnetohydrodynamics. The application of advanced processing techniques enable high data acquisition rates and concurrently a high spatial resolution facilitating to resolve transient liquid metal flow structures which could not been acquired so far. An experimental setup utilizing liquid metal in a cubic vessel exposed to a stationary rotating magnetic field was used to validate the reliability of the measurement system. The swirling fluid motion in its horizontal section could be resolved into a velocity field grid of 24 × 24 vectors while achieving frame rates of about 30 fps. Results from a further study driving liquid metal in a cylindrical vessel by a pulsed rotating magnetic field are presented. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling of the deformation processes and the localization of plastic deformations in the torsion–tension of solids of revolution

Generalized two-dimensional problems of the torsion of elastoplastic solids of revolution of arbitrary shape for large deformations under non-uniform stress-strain conditions are formulated and a method for their numerical solution is proposed. The use of this method to construct strain diagrams of materials based on experiments on the torsion of axisysmmetric samples of variable thickness until fracture occurs is described. Experimental and numerical investigations of processes of elastoplastic deformation, loss of stability and supercritical behaviour of solid cylindrical steel samples of variable thickness under conditions of monotonic kinematic loading with a torque, a tension and a combined load are presented. The mutual influence of torsion and tension on the deformation process and the limit states is estimated, and the universality (the independence of the form of the stress-strain state) of the “stress intensity – Odqvist parameter” diagram for steel for large deformations is proved.     

关于绕任意机翼非定常流动的一种无条件稳定的欧拉方程解

本文给出了绕二维与三维刚性或弹性振动机翼非定常无粘流动的欧拉方程解。首先利用Jameson的有限体积方法建立了求解欧拉方程的Runge-Kutta方法。为了提高受Runge-Kutta法稳定性限制的时间步长,文中采用了变系数的残值光顺方法。该方法避免了常系数残值光顺引起局部流场损失较大的问题。同时可在保证原计算格式精度要求下,大幅度提高计算时间步长,从而提高了计算效率。文中以二维与三维矩形机翼为例,分别对其在跨音速流场中作则性或弹性振动的非定常气动力进行了计算,研究了不同振动频率对流动产生的影响。部分计算结果与相应实验结果进行了比较。结果证明本方法是可靠的,可以用于求解绕任意运动机翼非定常流动问题。