Geometric Aspects of Degenerate Modulation Equations

Geometric aspects of degenerate modulation equations associated with spatially reversible systems are considered. Our primary observation is that stationary solutions of such equations always have a Poisson structure that is reminiscent of the equations governing the rigid body in mechanics. The Poisson structure is used to study the geometry of “spatial” phase space: A nontrivial Casimir of the Poisson structure provides a foliation of the phase space, spatially periodic states are given by critical points on level sets of the Casimir and stability type is given by the rate of change of the Casimir function. The bifurcation of spatially periodic states is then studied using singularity theory. The case where branches intersect transversely is treated in detail.     

A parameter estimation scheme for a class of sequential hybrid systems

It may happen that the equations governing the response of dynamical systems have some parameters whose values may not be known a priori and have to be obtained using parameter estimation schemes. In this article, we present a parameter estimation scheme for a class of sequential hybrid systems. By hybrid systems, we refer to those systems whose response is described by different governing equations corresponding to various regimes/modes of operation along with some criteria to switch between the same. In a sequential hybrid system, the different modes are arranged in a specific sequence and the system can switch from a given mode to either the previous mode or the following mode in this sequence. Here, we consider those systems whose governing equations consist of ordinary differential equations and algebraic equations. The conditions for switching between the various modes (referred to as transition conditions) are in the form of linear inequalities involving the system output. We shall first consider the case where the transition conditions are known completely. We present a parameter update scheme along with sufficient conditions that will guarantee bounded parameter estimation errors. Then, we shall consider the case where the transition conditions are not known in the sense that some parameters in these conditions are not known. We present a parameter estimation scheme for this case. We illustrate the performance of the parameter estimation scheme in both cases with some examples.     

Plane Sub-Critical Flow Past a Lifting Aerofoil

In this paper a numerical technique is used in the solutionof the equations governing the two-dimensional flow past a liftingaerofoil. The method as originated by Sells is modified by theuse of an extra conformal mapping such that more mesh pointsoccur in the region where the flow varies most. Thus more accurateresults are obtained without having to use substantially morecomputer storage or time. Results illustrating the method are presented for the 10% RAE101 section, at zero and non-zero incidence.     

Computer-based proof of the existence of superconvergence points in the finite element method; superconvergence of the derivatives in finite element solutions of Laplace's,Poisson's,and the elasticity equations

In this article we address the problem of the existence of superconvergence points for finite element solutions of systems of linear elliptic equations. Our approach is quite different from all other studies of superconvergence. We prove that the existence of superconvergence points can be guaranteed by a numerical algorithm, which employs a finite number of operations (provided that there is no roundoff-error). By employing this approach, we can reproduce all known results on superconvergence of finite element solutions for linear elliptic problems and we can obtain many new results. Here, in particular, we address the problem of the superconvergence points for the gradient of finite element solutions of Laplace's and Poisson's equations and we show that the sets of superconvergence points are very different for these two cases. We also study the superconvergence of the components of the gradient of the displacement, the strain and stress for finite element solutions of the equations of elasticity. For Laplace's and Poisson's equations (resp. the equations of elasticity), we consider meshes of triangular as well as square elements of degree p, 1 ? p ? 7 (resp. 1 ? p ? 4). For the meshes of triangular elements we investigate the effect of the geometry of the mesh by considering four mesh patterns that typically occur in practical meshes, while in the case of square elements, we study the effect of the element type (tensor-product, serendipity, or other). © 1996 John Wiley & Sons, Inc.     

Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems

The governing dynamics of fluid flow is stated as a system of partial differential equations referred to as the Navier-Stokes system. In industrial and scientific applications, fluid flow control becomes an optimization problem where the governing partial differential equations of the fluid flow are stated as constraints. When discretized, the optimal control of the Navier-Stokes equations leads to large sparse saddle point systems in two levels. In this paper, we consider distributed optimal control for the Stokes system and test the particular case when the arising linear system can be compressed after eliminating the control function. In that case, a system arises in a form which enables the application of an efficient block matrix preconditioner that previously has been applied to solve complex-valued systems in real arithmetic. Under certain conditions, the condition number of the so preconditioned matrix is bounded by 2. The numerical and computational efficiency of the method in terms of number of iterations and execution time is favorably compared with other published methods.     

The Design of axisymmetric ducts for incompressible flow with vorticity.

In this paper a numerical algorithm is described for solving the boundary value problem associated with axisymmetric, inviscid, incompressible, rotational (and irrotational) flow in order to obtain duct wall shapes from prescribed wall velocity distributions. The governing equations are formulated in terms of the stream function and the function as independent variables where for irrotational flow can be recognized as the velocity potential function, for rotational flow ceases being the velocity potential function but does remain orthogonal to the stream lines. A numerical method based on finite differences on a uniform mesh is employed. The technique described is capable of tackling the so–called inverse problem where the velocity wall distributions are prescribed from which the duct wall shape is calculated, as well as the direct problem where the velocity distribution on the duct walls are calculated from prescribed duct wall shapes. The two different cases as outlined in this paper are in fact boundary value problems with Neumann and Dirichlet boundary conditions respectively. Even though both approaches are discussed, only numerical results for the case of the Dirichlet boundary conditions are given. A downstream condition is prescribed such that cylindrical flow, that is flow which is independent of the axial coordinate, exists. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Cubic autocatalysis in a reaction–diffusion annulus: semi-analytical solutions

Semi-analytical solutions for cubic autocatalytic reactions are considered in a circularly symmetric reaction–diffusion annulus. The Galerkin method is used to approximate the spatial structure of the reactant and autocatalyst concentrations. Ordinary differential equations are then obtained as an approximation to the governing partial differential equations and analyzed to obtain semi-analytical results for this novel geometry. Singularity theory is used to determine the regions of parameter space in which the different types of steady-state diagram occur. The region of parameter space, in which Hopf bifurcations can occur, is found using a degenerate Hopf bifurcation analysis. A novel feature of this geometry is the effect, of varying the width of the annulus, on the static and dynamic multiplicity. The results show that for a thicker annulus, Hopf bifurcations and multiple steady-state solutions occur in a larger portion of parameter space. The usefulness and accuracy of the semi-analytical results are confirmed by comparison with numerical solutions of the governing partial differential equations.     

杆件轴向受迫振动的Galerkin有限元EEP法自适应求解

基于单元能量投影(element energy projection,EEP)法自适应分析在杆件静力问题以及离散系统运动方程组中所取得的成果,以直杆轴向受迫振动为例,研究并建立了一种在时间域和一维空间域同时实现自适应分析的方法.该方法在时间和空间两个维度都采用连续的Galerkin有限元法(finite element method,FEM)进行求解,根据半离散的思想,由空间有限元离散将模型问题的偏微分控制方程转化为离散系统运动方程组,对该方程组进行时域有限元自适应求解;然后再基于空间域超收敛计算的EEP解对空间域进行自适应,直至最终的时空网格下动位移解答的精度逐点均满足给定误差限要求.文中对其基本思想、关键技术和实施策略进行了阐述,并给出了包括地震波输入下的典型算例以展示该法有效可靠.     

离散的Burgers-Ginzburg-Landau方程组的吸引子

文章通过对空间变量的有限差分方法离散了具有周期边值的Ｂｕｒｇｅｒｓ Ｇｉｎｚｂｕｒｇ Ｌａｎｄａｕ方程组．研究了这个离散方程组初值问题解的适定性．证明了当差分网格足够大时离散方程组存在吸引子，并得到了吸引子的Ｈａｕｓｄｏｒｆｆ维数和分形维数的上界估计．这个上界不会随着网格的加细而无限增大，因此数值分析离散的有限维系统的吸引子可以近似探讨原无限维系统的吸引子．     

Alternative T-complete systems of shape functions applied in analytical Trefftz finite elements

The authors compare the behavior of hybrid Trefftz p-elements with two different types of shape functions identically fulfilling governing differential equations. Numerical examples include several boundary problems for Laplace, Poisson, and plane elasticity equations. Accuracy of the solutions, convergence properties, numerical stability and sensitivity for mesh distortion are investigated. It is shown that both systems of the functions can be efficiently applied, although they have different properties. © 1995 John Wiley & Sons, Inc.     

Dynamic motion of whirling rods with Coriolis effect

Longitudinal vibrations coupled with transverse vibrations of whirling rods are investigated. It is known that longitudinal and transverse vibrations are governed by second and fourth order differential equations, respectively. Due to the Coriolis effect, a system of equations that governs the longitudinal and transverse displacements will be constructed by coupling these two equations together. Solutions of the equations assume small oscillations of vibration being superimposed on the steady state of the whirling rod. Exact and approximate solutions are obtained from the proposed governing equations, where the approximate solutions on displacements and natural frequencies are acquired by neglecting the Coriolis effect. A proposed numerical scheme known as complete function collocation method is implemented to solve the governing equations coupled with longitudinal and transverse displacements. The approximate results on both longitudinal and transverse natural frequencies show that natural frequencies are decreasing while the angular velocity of the rod is increasing. Exact and numerical results on both longitudinal and transverse natural frequencies show that there are no predictable trends whether natural frequencies are increasing or decreasing while the angular velocity of the rod is increasing.     

二维流体力学变分通用公式

推导得到了二维流体力学变分通用公式,该公式适用于任何二维守恒型流体力学方程,得到的泛函受约于所谓的参数约束方程(控制方程中各参数间的相互关系式)。消除参数约束,我们可以十分方便地从通用公式导得广义变分原理。几个实例证明这种方法是有效的、简单的,并具有普遍的意义。     

Analysis of a bilinear finite element for shallow shells I: Approximation of inextensional deformations

We consider a bilinear reduced-strain finite element formulation for a shallow shell model of Reissner-Naghdi type. The formulation is closely related to the facet models used in engineering practice. We estimate the error of this scheme when approximating an inextensional displacement field. We make the strong assumptions that the domain and the finite element mesh are rectangular and that the boundary conditions are periodic and the mesh uniform in one of the coordinate directions. We prove then that for sufficiently smooth fields, the convergence rate in the energy norm is of optimal order uniformly with respect to the shell thickness. In case of elliptic shell geometry the error bound is furthermore quasioptimal, whereas in parabolic and hyperbolic geometries slightly enhanced smoothness is required, except for the degenerate cases where the characteristic lines are parallel with the mesh lines. The error bound is shown to be sharp.

     

Numerical simulations for the contraction flow using grid generation

We study the incomprssible Navier Stokes equations for the flow inside contraction geometry. The governing equations are expressed in the vorticity-stream function formulations. A rectangular computational domain is arised by elliptic grid generation technique. The numerical solution is based on a technique of automatic numerical generation of acurvilinear coordinate system by transforming the governing equation into computational plane. The transformed equations are approximated using central differences and solved simultaneously by successive over relaxation iteration. The time dependent of the vorticity equation solved by using explicit marching procedure. We will apply the technique on several irregularshapes.     

On the whole spectrum of Timoshenko beams. Part II: further applications

The problem of free vibrations of the Timoshenko beam model has been addressed in the first part of this paper. A careful analysis of the governing equations has shown that the vibration spectrum consists of two parts, separated by a transition frequency, which, depending on the applied boundary conditions, might be itself part of the spectrum. Here, as an extension, the case of a doubly clamped beam is considered. For both parts of the spectrum, the values of natural frequencies are computed and the expressions of eigenmodes are provided: this allows to acknowledge that the nature of vibration modes changes when moving across the transition frequency. This case is a meaningful example of more general ones, where the wave-numbers equation cannot be written in a factorized form and hence must be solved by general root-finding methods for nonlinear transcendental equations. These theoretical results can be used as further benchmarks for assessing the correctness of the numerical values provided by several numerical techniques, e.g. finite element models.     

General method for producing a boundary-fitted orthogonal curvilinear mesh

A general method is described for computing an orthogonal mesh fitted to a two-dimensional physical domain with arbitrary closed boundary. The method allows optimum control of mesh spacing through the introduction of arbitrary (with weak constraints) ‘packing’ functions into the elliptic governing equations. Two particular aspects are addressed: first, the condition on a scaling factor which normalizes the mesh aspect ratio; second, the condition for avoiding run-out of the mesh beyond the boundaries of the physical domain.Conversion of the equations to finite difference form and appropriate iterative techniques are discussed. Finally applications of the method in the context of flow across a bundle of rods are presented.     

A Geometric Theory of Nonlinear Morphoelastic Shells

Many thin three-dimensional elastic bodies can be reduced to elastic shells: two-dimensional elastic bodies whose reference shape is not necessarily flat. More generally, morphoelastic shells are elastic shells that can remodel and grow in time. These idealized objects are suitable models for many physical, engineering, and biological systems. Here, we formulate a general geometric theory of nonlinear morphoelastic shells that describes both the evolution of the body shape, viewed as an orientable surface, as well as its intrinsic material properties such as its reference curvatures. In this geometric theory, bulk growth is modeled using an evolving referential configuration for the shell, the so-called material manifold. Geometric quantities attached to the surface, such as the first and second fundamental forms, are obtained from the metric of the three-dimensional body and its evolution. The governing dynamical equations for the body are obtained from variational consideration by assuming that both fundamental forms on the material manifold are dynamical variables in a Lagrangian field theory. In the case where growth can be modeled by a Rayleigh potential, we also obtain the governing equations for growth in the form of kinetic equations coupling the evolution of the first and the second fundamental forms with the state of stress of the shell. We apply these ideas to obtain stress-free growth fields of a planar sheet, the time evolution of a morphoelastic circular cylindrical shell subject to time-dependent internal pressure, and the residual stress of a morphoelastic planar circular shell.     

Nonlinear Evolution Equations and the Braiding of Weakly Transporting Flows Over Gravel Beds

The flow of a river that transports sediment in the form of gravel as bedload is investigated for the case when the transport is small. The linear stability of such flows is discussed and used to formulate some strongly nonlinear investigations describing the interaction of bar instabilities that are known to occur. The key spatial scales in the asymptotic limit of small transport are identified, and highly nonlinear evolution equations derived for each case. A generalized KDV equation is found to govern the nonlinear evolution at small wavenumbers, while at O(1) wavenumbers an infinite set of "triad-like" amplitude equations describes the flow. The interactions demonstrate the natural tendency of rivers of width significantly higher than the critical width at which instability first occurs to form complex patterns that may be associated with braided rivers. The weak transport limit used in our anaysis makes our work directly relevant to rivers experiencing flood conditions where the onset of a flood causes transport to begin. The results shown suggest that in the highly nonlinear stages, bars take the form of slabs tilted in the flow direction with steep edges. In addition, it is found that there is no equilibrium state. These findings are consistent with observations.     

Finite Element Simulations with Adaptively Moving Mesh for the Reaction Diffusion System

A moving mesh method is proposed for solving reaction-diffusion equations. The finite element method is used to solving the partial different equation system, and an efficient scheme is applied to implement mesh moving. In the practical calculations, the moving mesh step and the problem equation solver are performed alternatively. Serveral numerical examples are presented, including the Gray-Scott, the Activator-Inhibitor and a case with a growing domain. It is illustrated numerically that the moving mesh method costs much lower, compared with the numerical schemes on a fixed mesh. Even in the case of complex pattern dynamics described by the reaction-diffusion systems, the adapted meshes can capture the details successfully.     

Application of a Krylov subspace iterative method in a multi-level adaptive technique to solve the mild-slope equation in nearshore regions

A multi-level adaptive numerical technique is applied to a nonlinear formulation of the mild-slope equation, to obtain the nearshore wave field, where the dominant processes of wave transformation are shoaling, refraction and diffraction. The advantage of this formulation over the traditional elliptic, parabolic and hyperbolic formulations is to require a lower minimum number of grid nodes per wavelength, thus, its capacity to predict the wave field for larger coastal areas. The efficiency of the interactions between the grid mesh levels, where two robust Krylov subspace iterative methods, the Bi-CGSTAB and the GMRES, are applied to solve the governing equation, is tested, for several hierarchies of grid mesh levels. The results show that the multi-level adaptive technique is efficient only if the GMRES iterative method is applied, and that for six grid mesh levels good results can be achieved for a residual as low as 10⁻³ for the finest grid.