Numerical Study of a Modified Time-Stepping θ-Scheme for Incompressible Flow Simulations

In [Turek (1996). Int. J. Numer. Meth. Fluids 22, 987–1011], we had performed numerical comparisons for different time stepping schemes for the incompressible Navier–Stokes equations. In this paper, we present the numerical analysis in the context of the Navier–Stokes equations for a modified time-stepping θ-scheme which has been recently proposed by Glowinski [Glowinski (2003). In: Ciarlet, P. G., and Lions, J. L. (eds.), Handbook of Numerical Analysis, Vol. IX, North-Holland, Amsterdam, pp. 3–1176]. Like the well-known classical Fractional-Step-θ-scheme which had been introduced by Glowinski [Glowinski (1985). In Murman, E. M. and Abarbanel, S. S. (eds.), Progress and Supercomputing in Computational Fluid Dynamics, Birkh?user, Boston MA; Bristeau et al. (1987). Comput. Phys. Rep. 6, 73–187], too, and which is still one of the most popular time stepping schemes, with or without operator splitting techniques, this new scheme consists of 3 substeps with nonequidistant substepping to build one macro time step. However, in contrast to the Fractional-Step-θ-scheme, the second substep can be formulated as an extrapolation step for previously computed data only, and the two remaining substeps look like a Backward Euler step so that no expensive operator evaluations for the right hand side vector with older solutions, as for instance in the Crank–Nicolson scheme, have to be performed. This modified scheme is implicit, strongly A-stable and second order accurate, too, which promises some advantageous behavior, particularly in implicit CFD simulations for the nonstationary Navier–Stokes equations. Representative numerical results, based on the software package FEATFLOW [Turek (2000). FEATFLOW Finite element software for the incompressible Navier–Stokes equations: User Manual, Release 1.2, University of Dortmund] are obtained for typical flow problems with benchmark character which provide a fair rating of the solution schemes, particularly in long time simulations.Dedicated to David Gottlieb on the occasion of his 60th anniversary     

Implicit–Explicit Schemes for BGK Kinetic Equations

In this work a new class of numerical methods for the BGK model of kinetic equations is presented. In principle, schemes of any order of accuracy in both space and time can be constructed with this technique. The methods proposed are based on an explicit–implicit time discretization. In particular the convective terms are treated explicitly, while the source terms are implicit. In this fashion even problems with infinite stiffness can be integrated with relatively large time steps. The conservation properties of the schemes are investigated. Numerical results are shown for schemes of order 1, 2 and 5 in space, and up to third-order accurate in time.     

Viscoelastic State for the Solutions of Spherically Symmetric Viscoelasticity Problems

The solutions of the spherically symmetric, linear, isothermal, and transient viscoelasticity problems via reciprocity theorem have been investigated for a specific material. The integral form of stress–strain relations has been used. The Laplace transform of a viscoelastic state, which is necessary for the integral equation arising as a result of reciprocity theorem, has been derived. This integral equation has been solved by Laplace transform. A sample problem has been solved to test the presented formulation. A numerical application of the analytic solution of this problem has been given.     

石蜡微胶囊中脲醛树脂壁材渗透性的研究

利用差示扫描量热法研究了石蜡微胶囊中脲醛树脂壁材对水及甲苯的渗透性 ,结果表明 ,水不能透过脲醛树脂壁材 ,而甲苯则易透过脲醛树脂壁材。因此 ,该石蜡微胶囊不能用于芳香族溶剂存在的场合。     

一种复合型数值方法的改进与算法实现

文中对原有的复合方法进行了改进,提出了一种新的复合型数值方法,这种方法解决了原方法求解含源项的问题时存在数值震荡的缺点。为了验证算法的有效性,对计算流体力学中一阶浅水波问题的两个典型的算例进行了数值模拟,在解的光滑性、锐利性等几个方面做了分析,结果表明改进后的方法既消除了数值震荡现象,又能准确地逼近真解,求出锐利的间断,是一种行之有效的数值方法。     

试论双点后方交会简便公式

根据双点后方交会的图形特点和已知条件，列出了方程组，推导出了双点后方交会直接计算座标的简便公式，同时求出了精度计算公式，以供测量人员参考     

Nonlinear electrical transport in K-TCNQ at low temperatures

A remarkable nonlinearity has been observed for electric transport in K-TCNQ (potassium tetracyanoquinodimethane) quasi-one-dimensional conductor. The negative differential resistance region appears afterV _m (turn over voltage), beyond which voltage goes down with increase of current. A possible mechanism is discussed in terms of dynamics of charged solitons and domain walls in one-dimensional molecular stacks of these types of crystals.     

AHMED experiments on hygroscopic and inert aerosol behaviour in LWR containment conditions: experimental results

Hygroscopic NaOH, CsI, CsOH and inert Ag aerosol behaviour at different temperatures and relative humidities (RH) has been studied in a well instrumented and controlled vessel of 1.81 m³ total free volume. Homogeneous thermal-hydraulic conditions for aerosol measurement in the vessel were achieved. The aerosol number and mass concentration were measured continuously during the experiments using a Condensation Nucleus Counter and a Tapered Element Oscillating Microbalance. The particle size distribution and chemical composition in the test conditions were measured by Berner low pressure impactors. In the case of NaOH the half life of the aerosol mass concentration was more than four times longer at low RH (22%) as compared to high RH (96%). The half lives of the CsOH and CsI aerosols were only twice as long at low RH as compared to high RH. Thus at high RH (96–97%) the half lives of CsOH and CsI were twice as long as the half life for the NaOH aerosol. The faster decay of the NaOH aerosol is due to the smaller density decrease of NaOH during water condensation. CsOH particles grew rapidly to their equilibrium size at all humidities. The measured equilibrium size for CsOH aerosol agree well with the calculated particle size at different RHs. Experimental results were also compared with calculations obtained by severe accident computer codes. These calculated results will be presented in a later paper.     

Effect of tap spacing on observed differential pressure fluctuations in a pneumatic transport line

The standard deviation of differential pressure fluctuations between pressure taps in a 31.7 mm i.d. vertical downer was found to be proportional to the square root of the distance between the taps when measured under otherwise constant conditions. This finding confirms the prediction of a theoretical model based on the Central Limit Theorem of sampling statistics. Although the paper only presents experimental confirmation for downward flow, the theoretical model should be valid for flows in all orientations. Except at very low solids flows, the pressure gradient was found to be positive in the flow direction, indicating that rise in pressure due to static pressure more than compensated for frictional pressure losses.     

泥沙输移矢量式及其应用

本文基于Parker等人泥沙输移矢量式的分析和研究，增加了上举力的作用，并考虑了泥沙颗粒受力非平衡条件，采用牛顿力学的矢量分析法，推导了泥沙输移的通用矢量式。经对公式简化，且用Ikeda和笔者试验资料的验证，证明计算值与试验值吻合较好。