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.     

Markov analysis of qualitative dynamics1

Common sense sometimes predicts events to be likely or unlikely rather than merely possible. We extend methods of qualitative reasoning to predict the relative likelihoods of possible qualitative behaviors by viewing the dynamics of a system as a Markov chain over its transition graph. This involves adding qualitative or quantitative estimates of transition probabilities to each of the transitions and applying the standard theory of Markov chains to distinguish persistent states from transient states and to calculate recurrence times, settling times, and probabilities for ending up in each state. Much of the analysis depends solely on qualitative estimates of transition probabilities, which follow directly from theoretical considerations and which lead to qualitative predictions about entire classes of systems. Quantitative estimates for specific systems are derived empirically and lead to qualitative and quantitative conclusions, most of which are insensitive to small perturbations in the estimated transition probabilities. The algorithms are straightforward and efficient.     

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     

On a three-dimensional axisymmetric boundary-value problem of nonlinear elastic deformation: Asymptotic solution and exponentially small error

In this paper, we study a three-dimensional axisymmetric boundary-value problem of a slender cylinder composed of a nonlinearly elastic material subjected to an axial force. Starting from the field equations, after a transformation and proper scalings, we identify a small variable and two small parameters, which characterize the present problem. Then, by an approach involving compound series-asymptotic expansions, a nonlinear ODE is derived, which governs the axial strain (the first-term in the series expansion). By imposing the zero radial displacement conditions at two ends, we manage to get the analytical solution of the axial strain, from which all other physical quantities can be deduced and thus the three-dimensional displacement field can be determined. Graphical results are presented, which show that there are two boundary layers near the two ends while the middle part is in a state of almost uniform extension. The asymptotic structure of the analytical solution is derived, which offers clear explanations to the structure of the deformed configuration and shows that the thickness of both boundary layers is of the order of the radius. We also point out the relevance of the present results to the St. Venant’s problem. In particular, we obtain the explicit uniformly-valid exponentially small error term, when the obtained deformed configuration is compared to the configuration of a uniform extension.     

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.     

Viscoelastic State for the Solutions of Spherically Symmetric Viscoelasticity Problems

The solutions of the spherically symmetric, linear, isothermal, and transient viscoelasticity problems via reciprocity theorem have been investigated for a specific material. The integral form of stress–strain relations has been used. The Laplace transform of a viscoelastic state, which is necessary for the integral equation arising as a result of reciprocity theorem, has been derived. This integral equation has been solved by Laplace transform. A sample problem has been solved to test the presented formulation. A numerical application of the analytic solution of this problem has been given.     

基于Delphi的报童模型的实现与动态仿真

已往库存论中的决策模型,它是把现实过程用数学的方法进行分析,得到决策值,但不能看到决策的过程,不能反映出现实生活中的不确定性、动态性。运用Delphi 7.0和SQL server 2000对销售过程进行了动态仿真,利用Newsboy模型的推广(S,s)订购策略来监控库存并做出决策,这样就体现出此过程的不确定性、动态性;且举出了一个仿真商品销售的实例,对其中的一些关键技术,如何仿真、算法的实现等进行了阐述。     

一种复合型数值方法的改进与算法实现

文中对原有的复合方法进行了改进,提出了一种新的复合型数值方法,这种方法解决了原方法求解含源项的问题时存在数值震荡的缺点。为了验证算法的有效性,对计算流体力学中一阶浅水波问题的两个典型的算例进行了数值模拟,在解的光滑性、锐利性等几个方面做了分析,结果表明改进后的方法既消除了数值震荡现象,又能准确地逼近真解,求出锐利的间断,是一种行之有效的数值方法。     

PMMA／PBA乳胶IPN阻尼材料的研究

用乳液聚合方法合成了一系列的互穿聚合物网络，其中包括：ＰＭＭＡ／ＰＢＡ和Ｐ（ＭＭＡ－ｃｏ－ＰＭＡ）／ＰＢＡ乳胶互穿聚合物网络（ＬＩＰＮ），ＰＭＭＡ／ＰＢＡ乳胶双向互穿聚合物网络（ＬＢＩＰＮ）以及ＰＭＭＡ／ＰＢＡ界面交联互穿聚合物网络（ＩＣＬＩＰＮ）。动态力学谱和拉伸试验结果表明：（１）ＰＭＭＡ／ＰＢＡ互穿聚合物网络的两组分相容性和阻尼性能可通过ＭＭＡ和ＰＭＡ共聚来改善；（２）双向互穿技术是一种行之