波动方程数值模拟的隐式差分法

在有限差分波动方程数值模拟中,通常采用高阶差分方法来提高空间导数的数值逼近精度,以实现降低数值频散,提高数值模拟精度的目的。首先对差分频散进行了理论分析;然后讨论了估计一阶空间导数的隐式差分格式,并与通常采用的高阶精度显式差分格式进行了对比分析,结果表明,隐式差分格式能够在更宽的波数范围使差分频散控制在可接受的水平,如8阶精度的显式差分格式所适应的波数带宽约为O．55kmax,而隐式差分格式所适应的波数带宽约为0．7kmax;最后通过模型试算,对隐式差分格式的有效性进行了验证。模拟结果表明,用隐式差分格式在一定程度上降低了差分频散,提高了模拟精度。     

A multigrid accelerated hybrid unstructured mesh method for 3D compressible turbulent flow

 A cell vertex finite volume method for the solution of steady compressible turbulent flow problems on unstructured hybrid meshes of tetrahedra, prisms, pyramids and hexahedra is described. These hybrid meshes are constructed by firstly discretising the computational domain using tetrahedral elements and then by merging certain tetrahedra. A one equation turbulence model is employed and the solution of the steady flow equations is obtained by explicit relaxation. The solution process is accelerated by the addition of a multigrid method, in which the coarse meshes are generated by agglomeration, and by parallelisation. The approach is shown to be effective for the simulation of a number of 3D flows of current practical interest. Sponsored by The Research Council of Norway, project number 125676/410 Dedicated to the memory of Prof. Mike Crisfield, a respected colleague     

Prediction of brittle-to-ductile transitions in polystyrene

In this study it is attempted to predict brittle-to-ductile transitions (BDTs) in polystyrene blends, induced either by an increase in temperature or by a decrease in inter-particle distance. A representative, two-dimensional volume element (RVE) of a polystyrene matrix with 20% circular voids, is deformed in tension. During deformation a hydrostatic-stress based craze-nucleation criterion [1] is evaluated. The simulations demonstrate that crazes initiate at low temperatures while a transition from crazing to shear yielding (BDT) is found around 75 °C. The numerical results correlate well with tensile tests on similar heterogeneous polystyrene. The presence of an absolute length, as experimentally found, is more difficult to explain. Near a free surface a T_g-depression is measured for polystyrene and also the resistance to indentation in polystyrene is lower than expected from bulk properties. Both observations are rationalised by an enhanced segmental mobility of chains near a free surface. As a consequence of these findings, an absolute length-scale could be incorporated in the numerical simulations. For simplicity, the length-scale is modelled by taking a temperature gradient over a thin layer near the internal free surfaces of the RVE. Deformation of the RVE with different absolute length-scales shows that indeed also the experimentally found brittle-to-ductile transition can be predicted if the ligament thickness between the inclusions (‘voids’) in polystyrene is below a critical value of ca. 15 nm.     

Field Static Load Test on Kao-Ping-Hsi Cable-Stayed Bridge

Field load testing is an effective method for understanding the behavior and fundamental characteristics of a cable-stayed bridge. This paper presents the results of field static load tests on the Kao-Ping-Hsi cable-stayed bridge, the longest cable-stayed bridge in Taiwan, before it was open to traffic. A total of 40 loading cases, including the unit and distributed bending and torsion loading effects, were conducted to investigate the bridge behavior. The atmospheric temperature effect on the variations of the main girder deflections was also monitored. The results of static load testing include the main girder deflections, the flexural strains of the prestressed concrete girder, and the variations of the cable forces. A three-dimensional finite-element model was developed. The results show that the bridge under the planned load test conditions has linear superposition characteristics and the analytical model shows a very good agreement with the bridge responses. Further discussion of deflection and cable forces of the design specifications for a cable-stayed bridge is also presented.     

Finite element piezothermoelasticity analysis and the active control of FGM plates with integrated piezoelectric sensors and actuators

 An efficient finite element model is presented for the static and dynamic piezothermoelastic analysis and control of FGM plates under temperature gradient environments using integrated piezoelectric sensor/actuator layers. The properties of an FGM plate are functionally graded in the thickness direction according to a volume fraction power law distribution. A constant displacement-cum-velocity feedback control algorithm that couples the direct and inverse piezoelectric effects is applied to provide active feedback control of the integrated FGM plate in a closed loop system. Numerical results for the static and dynamic control are presented for the FGM plate, which consists of zirconia and aluminum. The effects of the constituent volume fractions and the influence of feedback control gain on the static and dynamic responses of the FGM plates are examined. Received: 13 March 2002 / Accepted: 5 March 2003 The work described in this paper was supported by a grant awarded by the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. CityU 1024/01E).     

Numerical Study of Finite Element Method Based Solutions for Propagation of Wetting Fronts in Unsaturated Soil

The accurate prediction of the propagation of a wetting front in an unsaturated soil subjected to surficial infiltration is of practical importance to many geotechnical and geoenvironmental problems. The finite element method is the most common solution technique as the hydraulic soil properties are highly nonlinear. Two important issues are often found to create difficulties in such analyses. First, numerical oscillations are usually observed in the calculated pore pressures at the wetting front. Second, when a reasonable mesh size and time step are used, the elevation of the wetting front may be seriously overpredicted. This paper is focused on the second issue. The under-relaxation (UR) technique used in the iterative process within each time step is found to have a serious impact on rate of convergence with refinement in mesh size and time step. Two different techniques are typically used; the first evaluates the hydraulic conductivity using an average of heads calculated from the preceding time node and the most recent iteration of the current time node (UR1), and the second evaluates the hydraulic conductivity using the average of heads calculated from the two most recent iterations of the current time nodes (UR2). The study shows that UR1, which is adopted in programs such as SEEP/W, ensures that the solution converges rapidly to a stable solution within a time step, but may converge to the wrong wetting front at a given elapsed time unless a sufficiently refined mesh is used. UR2 converges much more slowly within a time step, but the error in the wetting front is smaller than that generated by UR1.     

Nonreactive atom/molecule-surface scattering within the finite basis wave packet method

A quantum wave packet code for studying nonreactive scattering of closed-shell atoms or diatomic molecules from a rigid surface is described. The time evolution relies on the Chebychev propagator. Up to 5 collider degrees of freedom, 3 in translation and 2 in rotation, are treated in a pseudospectral way with the momentum or finite basis representation as the primary space. Potential matrix elements are efficiently evaluated by means of sequential 1D transformations between momentum and coordinate spaces. Fast Fourier transforms are performed for the translational and azimuthal coordinates whereas a Gauss-Legendre transform is used for the polar coordinate. This pseudospectral strategy minimizes memory requirements because no off-diagonal Hamiltonian matrix elements need to be stored. In addition, a wide variety of physical systems can be studied since no particular functional form is imposed for the interaction potential.     

联产供冷与电力供冷能耗比较分析

文中将两种供冷能耗比较分解成用电煤耗差、用热煤耗量、增加产汽用电煤耗量三部分进行分析，按照热量法得出在目前一般情况下联产供冷比电力供冷通常费能和定性定量结论，为全面研究比较热电冷联产与分产的能耗奠定一个基础。