Transient Stochastic Analysis of Biodegradable Contaminant Transport: First-Order Decay

A stochastic methodology was used to analyze the field-scale transport of solutes in heterogeneous aquifers with first-order biodegradation. Spectral methods and perturbation techniques were utilized to develop expressions for the field-scale effective parameters in the mean transport model. Expressions were obtained for the longitudinal and transverse macrodispersivity coefficients, and effective velocity and an effective decay parameter for statistically anisotropic, and isotropic, heterogeneous porous medium, respectively. The behavior of these parameters was described as function of time and log K correlation scale. The expressions for asymptotic values of the dispersivity coefficients and effective velocity and decay parameters were also derived for the isotropic case.     

Compact finite difference-Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported. The project supported by the National Natural Science Foundation of China     

透射边界条件在波动谱元模拟中的实现:二维波动

将邢浩洁和李鸿晶提出的多次透射公式(multi-transmitting formula,MTF)的谱元格式应用于均匀介质中线弹性SH波动问题的谱元模拟.假定紧邻人工边界的一层谱单元为具有直线边界的四边形单元,以保证每个人工边界节点都唯一对应一条指向内域的离散网格线.人工边界节点在某时刻的位移由该离散网格线上的节点在前若干时刻的位移确定,按照MTF谱元格式进行计算.通过平面波以一定角度传播的外源问题算例和点源脉冲自由扩散的内源问题算例,验证了方法的可行性以及对实际复杂波动问题的适用性.通过不同类型初值问题算例,在时域内分析了插值多项式阶次、人工波速和透射阶次三个参数对反射误差的影响.结果表明:插值多项式阶次较高的格式会表现出更好的精度,但总体上对反射误差的影响较小;人工波速对反射误差具有显著影响,当人工波速小于介质物理波速时反射误差较大,而当人工波速等于或稍大于介质物理波速时反射误差处于较低水平;透射阶次对反射误差具有决定性影响,表现在不失稳的情形下提高透射阶次能够迅速降低反射误差,但内源问题从三阶MTF开始出现飘移失稳,外源问题从二阶MTF开始出现轻微的飘移失稳.     

基于谱元法的空间刚架动力学特性分析

以空间刚架结构为研究对象,应用谱元法研究其动力学特性,依次推导出杆、轴和铁木辛柯梁的动力学刚度阵,进而得到空间刚架的动力学刚度阵和频域响应的理论方程,拓展了谱元法应用范围。编程计算获得空间刚架的固有频率和点的谐响应曲线,与有限元分析相比较,进一步得到载荷下的瞬态响应。研究表明,谱元法在刚架结构的动力学特性分析中,比有限元方法高效高精度,更适宜解决中高频振动问题。     

NEW METHOD FOR LOW ORDER SPECTRAL MODEL AND ITS APPLICATION

In order to overcome the deficiency in classical method of low order spectral model, a new method for low order spectral model was advanced. Through calculating the multiple correlation coefficients between combinations of different functions and the recorded data under the least square criterion, the truncated functions which can mostly reflect the studied physical phenomenon were objectively distilled from these data. The new method overcomes the deficiency of artificially selecting the truncated functions in the classical low order spectral model. The new method being applied to study the inter-annual variation of summer atmospheric circulation over Northern Hemisphere, the truncated functions were obtained with the atmospheric circulation data of June 1994 and June 1998. The mechanisms for the two-summer atmospheric circulation variations over Northern Hemisphere were obtained with two-layer quasi-geostrophic baroclinic equation.     

GENERALIZED FINITE SPECTRAL METHOD FOR 1D BURGERS AND KDV EQUATIONS

A generalized finite spectral method is proposed. The method is of high-order accuracy. To attain high accuracy in time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Her-mite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation (nonlinear convection- diffusion problem) and KdV equation (single solitary and 2-solitary wave problems), where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.     

Homogeneous turbulence subjected to rotation-dominated plane distortion

A new configuration of distorted homogeneous turbulence is investigated in the domain of rotation-dominated elliptical flows. The experimental results are compared with exact numerical solutions obtained in Fourier space for the linear part of the problem. The inherently periodic character of the flow results in typical oscillatory variations of the anisotropy and the pressure-strain correlations. A new two-point closure approach of the EDQNM type is proposed for the nonlinear problem; it is based on a representation of the spectral tensor by scalar functions expanded in terms of spherical harmonics.     

Analysis of periodic-quasiperiodic nonlinear systems via Lyapunov-Floquet transformation and normal forms

In this paper a general technique for the analysis of nonlinear dynamical systems with periodic-quasiperiodic coefficients is developed. For such systems the coefficients of the linear terms are periodic with frequency ω while the coefficients of the nonlinear terms contain frequencies that are incommensurate with ω. No restrictions are placed on the size of the periodic terms appearing in the linear part of system equation. Application of Lyapunov-Floquet transformation produces a dynamically equivalent system in which the linear part is time-invariant and the time varying coefficients of the nonlinear terms are quasiperiodic. Then a series of quasiperiodic near-identity transformations are applied to reduce the system equation to a normal form. In the process a quasiperiodic homological equation and the corresponding ‘solvability condition’ are obtained. Various resonance conditions are discussed and examples are included to show practical significance of the method. Results obtained from the quasiperiodic time-dependent normal form theory are compared with the numerical solutions. A close agreement is found.     

Inertial manifolds for nonautomous infinite dimensional dynamical systems

In this paper, the long time behavior of nonautonomous infinite dimensional dynamical systems is discussed. Under the spectral gap condition, It is proved that there exist inertial manifolds for a class of nonautonomous evolution equations. Project supported by the National Natural Science Foundation of China     

Lamb波在复合材料板中传播的谱有限元模拟

众所周知Lamb波在复合材料中的传播呈各向异性的特点,经典有限元法模拟这类问题效率不高,所以,本文采用谱有限元法进行研究。先建立了一种新的谱有限板单元,该单元以Gauss-Lobatto-Legendre点作为节点,使质量矩阵是对角矩阵;另外,该单元采用了扩展的位移场,能够较好地模拟板结构的三维特性。然后,对复合材料板结构中Lamb波在对称模式与反对称模式下的传播速度进行了求解,将计算结果与Mindlin板谱单元的结果以及三维弹性理论的结果进行了比较,并讨论了Lamb波在反对称层合板中的传播特点。最后,模拟了Lamb波在含和不含损伤复合材料层合板中的传播,数值结果表明所建立的谱有限板单元可以较好地模拟出Lamb波在复合材料板结构中的传播特性。