Calculation of impulsively started incompressible viscous flows

The paper's focus is the calculation of unsteady incompressible 2D flows past airfoils. In the framework of the primitive variable Navier–Stokes equations, the initial and boundary conditions must be assigned so as to be compatible, to assure the correct prediction of the flow evolution. This requirement, typical of all incompressible flows, viscous or inviscid, is often violated when modelling the flow past immersed bodies impulsively started from rest. Its fulfillment can however be restored by means of a procedure enforcing compatibility, consisting in a pre‐processing of the initial velocity field, here described in detail. Numerical solutions for an impulsively started multiple airfoil have been obtained using a finite element incremental projection method. The spatial discretization chosen for the velocity and pressure are of different order to satisfy the inf–sup condition and obtain a smooth pressure field. Results are provided to illustrate the effect of employing or not the compatibility procedure, and are found in good agreement with those obtained with a non‐primitive variable solver. In addition, we introduce a post‐processing procedure to evaluate an alternative pressure field which is found to be more accurate than the one resulting from the projection method. This is achieved by considering an appropriate ‘unsplit’ version of the momentum equation, where the velocity solution of the projection method is substituted. Copyright © 2004 John Wiley & Sons, Ltd.     

Damage and fracture evaluation of granular composite materials by digital image correlation method

This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steelfiber reinforced concrete, sandstone and crystal-polymer composite. The deformation fields of the composite materials resulted from stress localization were obtained by the correlation computation of the surface images with loading steps and thus the related damage prediction and fracture parameters were evaluated. The correlation searching could be performed either directly based on the gray levels of the digital images or from the wavelet transform (WT) coefficients of the transform spectrum. The latter was developed by the authors and showed higher resolution and sensitivity to the singularity detection. Because the displacement components came from the rough surfaces of the composite materials without any coats of gratings or fringes of optical interferometry, both surface profiles and the deformation fields of the composites were visualized which was helpful to compare each other to analyze the damage of those heterogeneous materials. The project supported by the National Natural Science Foundation of China (10125211 and 10072002), the Scientific Committee of Yunnan Province for the Program of Steel Fiber Reinforced Concrete, and the Institute of Chemical Materials, CAEP at Mianyang     

低速二维流动的粒子仿真研究

本文运用信息保存法对低速二维的流动现象进行模拟，考察了低速条件下的有限平板绕流以及微槽道气体流动问题。研究表明：在对低速流动的模拟过程中，运用IP法在能够获得较好的结果的同时，具有比DSMC方法更高的计算效率。     

Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.     

The fifth dimension: Theodor Kaluza's ground‐breaking idea

Theodor Kaluza (1885–1954) attracted the attention of the physical community since 1921 with his unified field theory of gravitation and electromagnetism in five dimensions. Despite Einstein's great interest in Kaluza's theory, 50 years elapsed before it contributed toward a paradigm shift in modern theoretical physics. The biography of this still unknown scientist is briefly presented along with an outline of his work in physics. A short history of the theories of unification and the dimensionality of space‐time is followed by a discussion of the significance of Kaluza's five‐dimensional unified theory in modern physics from the point of view of superstring and M‐theory.     

Problem of the free oscillations of a capillary bridge

The free oscillations of a capillary bridge whose equilibrium shape is determined by the surface tension forces and the static gravity field are investigated. The values of 25 “lower” levels of the spectrum of natural oscillations of the capillary bridge are found for various control parameters in accordance with the experimental conditions.     

有限元节点编号优化

从分析最简单最优的矩形单元节点编号着手，分析其特点，得出各个节点的相邻单元节点编号总和除以该节点相邻单元数所得的节点商以及相邻单元最小最大节点编号和都为升序排列。总结出一种新型的有限元节点重编号方法，并举例说明该算法能减少有限元的带宽。     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

Maximal σ‐polynomials of connected 3‐chromatic graphs

In the set of graphs of order n and chromatic number k the following partial order relation is defined. One says that a graph G is less than a graph H if c_i(G) ≤ c_i(H) holds for every i, k ≤ i ≤ n and at least one inequality is strict, where c_i(G) denotes the number of i‐color partitions of G. In this paper the first ? n/2 ? levels of the diagram of the partially ordered set of connected 3‐chromatic graphs of order n are described. © 2003 Wiley Periodicals, Inc. J Graph Theory 43: 210–222, 2003     

Dispersion analysis of the least‐squares finite‐element shallow‐water system

The frequency or dispersion relation for the least‐squares mixed formulation of the shallow‐water equations is analysed. We consider the use of different approximation spaces corresponding to co‐located and staggered meshes, respectively. The study includes the effect of Coriolis, and the dispersion properties are compared analytically and graphically with those of the mixed Galerkin formulation. Numerical solutions of a test problem to simulate slow Rossby modes illustrate the theoretical results. Copyright © 2003 John Wiley & Sons, Ltd.