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.     

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

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

粉尘薄膜爆轰

与普通爆轰不同，粉尘薄膜爆轰波前氧化剂与燃料处于分离状态．由于激波作用下壁面燃料的注入才于壁面附近形成燃烧的粉尘云，燃烧释放的热使激波自持形成粉尘薄膜爆轰。本文在长2、8m，横截面20mm×20mm的方型水平激波管中对此类爆轰进行了实验研究．并根据作者建立的模型对其二维结构进行了数值模拟。     

热分层流动计算中若干差分格式的比较

本文应用常规的守恒方程及湍流模型，对二维池内热分层流动建立了数学模型，并在Ｐａｔａｎｋａｒ的压力校正法基础上对该数学模型进行了数值计算，计算中的差分格式为乘方格式、一阶迎风、二阶迎风和ＱＵＩＣＫ。通过与实验的比较分析，建议计算应采用ＱＵＩＣＫ格式，并在进行风格独立性试验的基础上使用加密网络。