Moving Least-Squares Differential Quadrature Method for Free Vibration of Antisymmetric Laminates

In this paper, the moving least-squares differential quadrature (MLSDQ) method is employed for free vibration of thick antisymmetric laminates based on the first-order shear deformation theory. The generalized displacements of the laminates are independently approximated with the centered moving least-squares (MLS) technique within each domain of influence. The MLS nodal shape functions and their partial derivatives are computed quickly through back-substitutions after only one LU decomposition. Subsequently, the weighting coefficients in the MLSDQ discretization are determined with the nodal partial derivatives of the MLS shape functions. The MLSDQ method combines the merits of both the differential quadrature and meshless methods which can be conveniently applied to complex domains and irregular discretizations without loss of implementation efficiency and numerical accuracy. The natural frequencies of the laminates with various edge conditions, ply angles, and shapes are calculated and compared with the existing solutions to study the numerical accuracy and stability of the MLSDQ method. Effects of support size, order of completeness of basis functions, and node irregularity on the numerical accuracy are investigated in detail.     

带参数的一阶脉冲微分方程边值问题

利用上下解和单调迭代法。讨论了带参数的一阶脉冲微分方程。获得了单边Lipschitz条件下这类方程边值问题的解的存在性。     

Estimation and Optimization of the Principal H-Maser Parameters

The scale along the coordinate axes in the equations of motion hitherto used for an H-maser (in relative units) is shown to depend on Q as the quality factor for an H-maser. This hinders both analysis of the operation of the H-maser and optimization of its parameters. The modified equations of motion proposed here do not have that shortcoming. We determine the atomic flux into the bulb that gives the H-maser the best metrological characteristics.     

明渠水流泥沙运动水深平均及力矩方程的理论和验证

首次提出了水流泥沙运动水深平均及其力矩方程的概念 ,并从理论上推导出这套方程组。提出水流泥沙运动力矩方程的目的有两个 ,一是在不增加空间维数的情况下可以有效地恢复水流泥沙因子沿水深分布信息 ,克服目前普遍使用的水深平均方程过多地忽略沿水深分布信息的缺点 ;二是增强模型的模拟能力。针对这两个目的是否能真正得到体现 ,分别从理论和实验室资料两个方面对新推导的力矩方程进行了全面的检验。理论验证结果表明 ,一维的力矩方程除了能模拟水深平均的流速和含沙量 ,还能够较好地模拟水流泥沙因子沿水深分布信息 ,这充分说明了一维力矩方程的二维特性。为了检验力矩方程的模拟能力 ,特别选取了一些极端的水流泥沙运动实验 ,如清水冲刷、点源加沙和深槽回淤等。模拟结果表明 ,力矩方程对于这些极端的水沙运动具有较好的模拟能力     

Interaction of free-surface waves with floating flexible strips

The method developed by the author to derive a set of algebaic equations to solve the interaction of free-surface waves with a single floating rigid or flexible two-dimensional platform with small draft is extended to the case that the platform consists of strips with different constant flexural rigidity and mass. The method is based on the application of Green's theorem, with a specific choice of the Green function to arrive at a differential-integral equation along the platform. This equation can be solved exactly by means of superposition of exponential functions, a standard method to solve a set of linear differential equations. After integration with respect to the space coordinate the residue theorem leads to both the dispersion relation along each individual strip and an algebraic equation for the coefficients. Due to very fast convergence with respect to the number of coefficients taken into account the series are truncated. Depending on the water-depth, in each series three to ten terms are taken into account. Results are shown for a structure consisting of several strips that are tightly connected and for disjoint strips. In the latter case the computation of the water level between the strips is also computed. The water level and the reflection and transmission coefficients are not unknowns in the algebraic equation, but are computed afterwards by means of Green's theorem.     

On solving systems of ordinary differential equations on MIMD-computers

The error of solution of Cauchy problems for systems of ordinary differential equations is estimated in the case where the input data are approximate. It is shown how to prepare a program for computing the right-hand sides of the system automatically and simultaneously. Diagrams are presented to illustrate the efficiency of parallelization. __________ Translated from Kibernetika i Sistemnyi Analiz, No. 2, pp. 175–182, March–April 2007.     

Coupling of free surface and groundwater flows

In this paper we present some preliminary results about the coupling of shallow water equations for free surface flows and Darcy equation for groundwater flows. A suitable set of interface conditions is discussed: the Beavers and Joseph formula for the bottom stress is used. An iterative algorithm to solve the coupled problem is proposed and some numerical results are presented.     

Stabilization Methods for Spectral Element Computations of Incompressible Flows

A stabilization method for the spectral element computation of incompressible flow problems is investigated. It is based on a filtering procedure which consists in filtering the velocity field by a spectral vanishing Helmholtz-type operator at each time step. Relationship between this filtering procedure and SVV-stabilization method, introduced recently in [JCP, 2004, 196(2), p680], is established. A number of numerical examples are presented to show the accuracy and stabilization capability of the method.     

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.