定常的热传导-对流问题的Galerkin/Petrov最小二乘混合元方法

1.引言 热传导-对流问题是大气动力学中的一个重要的方程,这个方程组也称为强迫耗散的非线性系统方程组,其较Navier-Stokes方程多了一个未知函数温度场,且温度与速度和压力之间存在着复杂的非线性关系.从热动力学可知,任何运动都会产生热量即有温度,而且温度与速度和压力之间必定互相转化,因此对该非线性系统的研究更具有实际意义.[1]先对     

Stability analysis of abnormal multiplication of plankton using parameter identification technique

This paper presents the mathematical approach for the abnormal multiplication of plankton. An abnormal multiplication can be expressed as an unstable problem and the stability of the system is investigated by introducing eigenvalues of a mathematical equation. The stability of the system can be judged by an eigenvalue based on the Lyapunov's stability theory. In this paper, the Arnoldi‐QR method is used to obtain eigenvalues and eigenvectors of the system. The mode superposition method is employed to create spatial distribution needed to analyse the stability. To obtain the objective eigenvalue, the parameter identification technique is employed. The finite element method is used for the discretization in space. Lake Kasumigaura, which is located in Ibaraki Prefecture in Japan, is selected and actual data in 1975, 1976, 1991 and 2000 are used in order to investigate the stability of the specified lake in Japan. Copyright © 2004 John Wiley & Sons, Ltd.     

Indirect boundary element method for unsteady linearized flow over prolate and oblate spheroids and hemispheroidal protuberances

The indirect boundary element method was used to study the hydrodynamics of oscillatory viscous flow over prolate and oblate spheroids, and over hemispheroidal bodies hinged to a plate. Analytic techniques, such as spheroidal coordinates, method of images, and series representations, were used to make the numerical methods more efficient. A novel method for computing the hydrodynamic torque was used, since for oscillatory flow the torque cannot be computed directly from the weightings. Instead, a Green's function for torque was derived to compute the torque indirectly from the weightings. For full spheroids, the method was checked by comparing the results to exact solutions at low and high frequencies, and to results computed using the singularity method. For hemispheroids hinged to a plate, the method for low frequencies was checked by comparing the results to previous results, and to exact solutions at high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

A Constitutive Model in Finite Viscoelasticity of Particle-reinforced Rubbers

Two series of tensile relaxation tests are performed on natural rubber filled with high abrasion furnace black. To fit observations, constitutive equations are derived for the nonlinear viscoelastic behavior of a particle-reinforced elastomer. A filled rubber is modeled as a composite medium, where inclusions with low concentrations of junctions are randomly distributed in the host matrix. The inclusions are treated as equivalent networks of macromolecules, where strands can separate from temporary junctions as they are thermally agitated. The bulk medium is thought of as a permanent network of chains. Unlike conventional concepts of transient networks, the concentration of strands in inclusions is assumed to be affected by mechanical factors: under active loading, inter-chain interactions weaken and some strands that were prevented from detachment from their junctions in a stress-free compound become free to separate from the junctions in a deformed medium. Unloading strengthens interactions between macromolecules, which results in an increase in the number of permanent strands. By using the laws of thermodynamics, stress–strain relations for a particle-reinforced rubber are developed. Adjustable parameters in the constitutive equations are found by fitting the experimental data. It is demonstrated that mechanical pre-loading and annealing of specimens at an elevated temperature noticeably affect concentrations of inclusions with various activation energies for rearrangement of strands.     

Non-elliptic deformation field near the tip of a mixed-mode crack in a compressible hyperelastic material

This paper gives an asymptotic analysis of the deformation field near the tip of an arbitrary mixed-mode crack in a compressible hyperelastic harmonic material which loses ellipticity at sufficiently large deformations. It is found that the near-tip deformation field is characterized by a localized non-elliptic deformation band issuing from the crack-tip and bounded by two curves of discontinuous deformation gradient. Explicit expression for the near-tip deformation field is obtained both inside and outside the localized deformation band. In particular, a simple relation is derived that determines the orientation of the deformation band in terms of two complex governing parameters of the near-tip fields inside and outside the deformation band, respectively.     

Analysis of a finite-difference scheme for a singularly perturbed problem with two small parameters

We study a model linear convection-diffusion-reaction problem where both the diffusion term and the convection term are multiplied by small parameters ε_d and ε_c, respectively. Depending on the size of the parameters the solution of the problem may exhibit exponential layers at both end points of the domain. Sharp bounds for the derivatives of the solution are derived using a barrier-function technique. These bounds are applied in the analysis of a simple upwind-difference scheme on Shishkin meshes. This method is established to be almost first-order convergent, independently of the parameters ε_d and ε_c.     

Performance of noise barriers with various edge shapes and acoustical conditions

The performances of barriers having different shapes and surface conditions were tested using the boundary element method in a well-controlled environment. The heights and widths of the barriers were standardized and the insertion losses for six receiver positions were averaged and compared. Figures displaying the results allow for straightforward barrier performance estimation. It was shown that absorbing and soft edges significantly improve the efficiency of the barrier, but configuration modifications provide only a slight improvement. The soft T-shaped barrier produces the highest performance. A 3 m high T-shaped barrier provides the same performance as a 10 m high plain barrier. The spectral efficiency was also investigated. The insertion loss spectra for the absorbing and the soft barriers exhibit a similar shape, but the rigid barrier differs from these two.     

O～3岁儿童微量元素检测及纠正150例观察与体会

对西安地区150例0～3岁儿童进行了六种微量元素：钙、锌、铜、铁、镁、硒的检测，发现微量元素缺乏可导致急慢性感染、免疫力低下、小儿营养不良、贫血、发育迟缓等症状。并对小儿营养不良、喂养不当进行矫正。其中铁、钙、锌同时缺乏占50％。     

A study on the applications of the acoustic design sensitivity analysis of vibrating bodies

The determination of the sensitivity of the acoustical characteristics of vibrating systems with respect to the variation of the design parameters predicting these characteristics is a necessary and important step of the acoustic design and optimization process. Acoustic design sensitivity analysis includes the computation and evaluation of the sensitivity information required for this procedure. In this study, a boundary element code performing the sensitivity analysis of the acoustic pressure by using the matrix sensitivities with respect to different design variables has been developed. The effect of the precision of boundary element discretization on the acoustic pressure sensitivity is examined via this code. The formulation is applied to a multi-source system and the dimension sensitivity analysis of near field pressures of two-dilating-spherical source is performed. The last application is devoted to a real sound source: a washing machine sitting on the floor. Sensitivity of the field pressures to the machine’s dimensions (size), surface velocity and frequency is examined on the bases of the boundary element model of the machine and half-space condition. The impacts of these variables are compared; and a limiting speed for the machine responding both the acoustical and operational requirements is determined.