The Transport Properties of Hot Electrons in a Xenon/Methane Mixture

The results are given of numerical solution of the Boltzmann equation in a binomial approximation in view of elastic and inelastic electron collisions in Xe+CH₄ mixtures and in pure methane. The electron energy distribution functions obtained are used to calculate the electron transport coefficients for E/N values of up to several Townsends, i.e., the drift velocity, mobility, average and characteristic energies, diffusion coefficient. The results of calculation for pure methane fit the available experimental data. A similarity rule is found for the electron transport coefficients in a Xe+CH₄ mixture with different concentrations of methane molecules, which enables one to determine the values of transport coefficients in a mixture with a minor (less than 30 percents) methane content.     

The Difference in Turbulent Diffusion between the Active and Passive Scalars in a Thermally Stably Stratified Medium

The difference in the turbulent diffusion between the active (heat) and passive (mass) scalars in a thermally stably stratified medium is investigated. The axisymmetric problem is treated on the formation of a turbulent circulation flow above a heated disk and on the turbulent diffusion of a passive scalar (impurity) from a continuous surface source in a stably stratified medium. The results indicate that the thermal stratification causes appreciable differences in the coefficients of turbulent transfer between the active (heat) and passive (mass) scalars. This means that the assumption of the identity of the coefficient of turbulent diffusion of heat and mass, employed in conventional models of turbulence, produces significant errors in estimating the heat and mass transfer in a thermally stably stratified medium.     

Bifurcation points and bifurcated branches by an asymptotic numerical method and Padé approximants

The aim of this work is to develop a reliable and fast algorithm to compute bifurcation points and bifurcated branches. It is based upon the asymptotic numerical method (ANM) and Padé approximants. The bifurcation point is detected by analysing the poles of Padé approximants or by evaluating, along the computed solution branch, a bifurcation indicator well adapted to ANM. Several examples are presented to assess the effectiveness of the proposed method, that emanate from buckling problems of thin elastic shells. Especially problems involving large rotations are discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

A three‐dimensional model describing stress‐temperature induced solid phase transformations: solution algorithm and boundary value problems

An always increasing knowledge on material properties as well as a progressively more sophisticated production technology make shape memory alloys (SMA) extremely interesting for the industrial world. At the same time, SMA devices are typically characterized by complex multi‐axial stress states as well as non‐homogeneous and non‐isothermal conditions both in space and time. This aspect suggests the finite element method as a useful tool to help and improve application design and realization. With this aim, we focus on a three‐dimensional macroscopic thermo‐mechanical model able to reproduce the most significant SMA features (Int. J. Numer. Methods Eng. 2002; 55 : 1255–1264), proposing a simple modification of such a model. However, the suggested modification allows the development of a time‐discrete solution algorithm, which is more effective and robust than the one previously discussed in the literature. We verify the computational tool ability to simulate realistic mechanical boundary value problems with prescribed temperature dependence, studying three SMA applications: a spring actuator, a self‐expanding stent, a coupling device for vacuum tightness. The effectiveness of the model to solve thermo‐mechanical coupled problems will be discussed in a forthcoming work. Copyright © 2004 John Wiley & Sons, Ltd.     

A direct numerical simulation of transition phenomena in a mixed convection channel flow

This study intends to provide an increased understanding of the laminar-turbulent transition phenomena for the buoyancy-assisted heated vertical channel flow during the early transient stage. The spectral method with weak formulation is applied in the direct numerical simulation. Initial disturbances consist of the finite-amplitude two-dimensional TS wave and a pair of three-dimensional oblique waves for the K-type disturbances. The results from the harmonic energy competitions of different wave modes show that for the buoyancy-assisted heated flow, the (k_x=1, k_z=1) or (1,1) and (1,0) modes would gain energy immediately and start to rise at almost the same rate. This phenomenon is different from that of the buoyancy-opposed flow, where the (1,1) mode decays slowly in the beginning until other modes gain enough energy and then it begins to grow quickly and overtakes the (1,0) mode after a short time period. These different transition patterns match with the experimental results that the flow transition is supercritical and subcritical for the buoyancy-assisted and -opposed flows, respectively. Buoyancy-assisted heated flow transition follows the general trend of an isothermal flow in the beginning, but the thermal-buoyant force is crucial in accelerating the instability and also causing notable differences during the subsequent transition process. All of the results for the vortex structures development, kinetic energy budget of the disturbances, flow visualization by tagged fluid particles, and the local temperature fluctuations are consistent in pointing to a clear pattern for the buoyancy-assisted heated flow transition.     

引水发电系统无压隧洞中弃水堰位置的数值模拟研究

水电站引水发电系统无压隧洞中的水力过渡过程是近年来电站设计中出现的一个新课题,特别是流域梯级电站开发中经常遇到。结合工程实例应用扩散法,给出了无压隧洞水力过渡过程数值模拟的基本方程以及上下游及弃水堰等边界附加方程,模拟了弃水堰布置在不同位置下,无压隧洞内的水位变化过程。研究成果和讨论对实际工程应用具有参考意义。     

Extension of a combined analytical/numerical initial value problem solver for turbulent mixing with combustion

Here, we describe the development of a reacting flow multi-species/combustion methodology, implemented as an extension to the differential reduced ejector analysis (DREA) computer program [Mathematical and computer modeling, vol. 31, 2000, p. 21; Appl. Math. Model. 25 (2001) 427; Comput. Math. Appl. 43(10–11); NASA Contractor Report, 1998]. Use of the single fluid IVP solver framework that was developed for the original DREA model has been directly coupled into the combustion formulation. With these modifications, the analysis has an elementary single step reaction Fuel+Oxidizer→Product combustion capability. Though approximate in nature, the simplicity and efficiency of the DREA formulation make it suitable for its original niche, namely design and preliminary design environments where more complex and expensive models may be inappropriate.     

Active set truncated-Newton algorithm for simultaneous optimization of distillation column

An active set truncated-Newton algorithm (ASTNA) is proposed to solve the large-scale bound constrained sub-problems. The global convergence of the algorithm is obtained and two groups of numerical experiments are made for the various large-scale problems of varying size. The comparison results between ASTNA and the subspace limited memory quasi-Newton algorithm and between the modified augmented Lagrange multiplier methods combined with ASTNA and the modified barrier function method show the stability and effectiveness of ASTNA for simultaneous optimization of distillation column.     

EXPERIMENTAL STUDIES AND NUMERICAL CALCULATION ON TURBULENT BOUNDARY LAYER OF WATER FLOW ALONG CURVED SURFACE

This paper will introduce experimental studies and numerical calculation onturbulent boundary layer of water flow along curved surface in our country in recent years.Onthe basis of the experimental studies,the effects of curvature and roughness on velocitydistribution and pressure distribution and the change of turbulent flow boundary layer onoverflow bucket concave surface is discussed.We proposed the empirical formula of velocity,pressure and the change of turbulent flow boundary layer on outlet bucket concave.According tothe momentum principle,we have deduced the momentum integral equation full water depthboundary layer and using the element as control unit inside the boundary layer on concavesurface of bucket.Combining with continuity equation,we have computed the boundary layerdevelopment on the bucket of a spillway.Compared with the field experimental data,thecalculation results are available.Under polar coordinates,a mathematical model for simulatingtime-averaged flow characteristics in concave surface of bucket is established.The turbulent flowfield on concave surface of bucket is calculated by SIMPLE method and this mathematicalmodel.The flow velocity field,pressure field,distribution of turbulent kinetic energy,distribution of turbulent energy dissipating rate and distribution of shear stress are available.Thecalculation value is consistent with measured test data.     

饱和—非饱和土体非稳定渗流数值分析

本文对一般的非饱和渗流有限元计算方法加以改进，有效地消除了非饱和渗流数值计算中存在的数值弥散现象。同时还提出了一种简便有效的逸出面处理新方法，并给出了非饱和非稳定渗流计算的实例。