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.     

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

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

Emission of trace toxic metals during pulverized fuel combustion of Czech coals

A study of the trace elements emission (As, Se, Cd, Co, Cr, Cu, Zn, Hg, Tl, Pb, Ni, Sn, Sb, V, Mn and Fe) from pulverized coal combustion has been made at six heating and power stations situated in the Czech Republic. The amount of chlorine in coal has considerable influence on volatilization of some elements such as Zn, Cu, Pb, Hg and Tl, which is explained by the formation of thermodynamically stable compounds of these elements with chlorine. Generally, the affinities for Cl follows the order Tl > Cu > Zn > Pb > Co > Mn > Sn > Hg. The experimental data indicates enrichment of some of the trace toxic elements in the emissions (Cu, Zn, As, Se, Cd, Sn, Sb, Hg and Pb) and good agreement was obtained by thermodynamic equilibrium calculations with a few exceptions. In the case of Fe, Mn, Co, Cr and Sn calculated values are overestimated in the bottom ash and there are zero predicted amounts of these elements in the fly ash. In comparison, the results from experiments show up to 80% of these elements retained in fly ash. This implies that there exist additional steps leading to the enrichment by Fe, Mn, Co, Cr and Sn of small particles. Such mechanisms could include the ejection during devolatilization of small inorganic particles from the coal of bottom ash particles, or disintegration of the char containing these metals to small particles of fly ash. On the other hand, there are slightly overestimated or similar values of relative enrichment factors for As, V, Cu, Cd, Sb, Tl and Pb in the fly ashes and zero predicted values for bottom ashes. Our experimental results show about 5% or less of these elements are retained in bottom ashes, so they probably remain in the bottom ash inside unburned parts of coal. Copyright © 2003 John Wiley & Sons, Ltd.     

油气集输系统新增产能建设拓扑优化方法研究

在油气集输新增产能建设中,为充分利用已建管网系统的处理能力,达到降低投资,提高经济效益的目的,进行了油气集输新增产能建设布局优化方法的研究。针对新建集输管网为树状结构的情况,建立了油气集输系统新增产能建设拓扑优化的数学模型,讨论了计算复杂性问题。根据模型的特点,将问题分解为拓扑级优化和几何级优化两个子问题,并通过它们之间的迭代进行求解。实例计算表明,所建立的数学模型正确,求解方法可行,能用于工程实际设计。     

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.     

一种弧长参数曲线及其计算机实现

本文提出了一种弧长参数曲线的数值算法，这种算法具有收敛迅速，迭代过程自动确定等优点。