Numerical simulation of viscous oscillatory flow past four cylinders in square arrangement

The results of a numerical study of the viscous oscillating flow past four circular cylinders, for a constant frequency parameter equal to 50 and KC ranging between 0.2 and 10, are presented. The cylinders were placed on the vertices of a square, two sides of which were perpendicular and two parallel to the oncoming flow, for pitch ratios, P/D, ranging between 2 and 5. The finite-element method was employed for the solution of the Navier-Stokes equations, in the formulation where the stream function and the vorticity are the field variables, whereas the pressure distribution throughout the computational domain was obtained from the solution of Poisson’s equation. When the Keulegan-Carpenter number is lower than 4, the flow remains symmetrical with respect to the horizontal axis of symmetry of the solution domain and periodic at consecutive cycles. As KC increases to 4, the flow becomes aperiodic in different cycles, although symmetry with respect to the horizontal central line of the domain is preserved. For KC equal to 5, asymmetries appear intermittently in the flow, which are eventually amplified as KC increases still further. These asymmetries, in association with the aperiodicity of flow in different cycles, lead to an almost chaotic configuration as KC grows larger. For characteristic cases the flow pattern and the time histories of the in-line and transverse forces exerted on the cylinders are presented. The mean transverse forces acting on the cylinders, the r.m.s. values of the in-line and transverse forces, together with the drag and inertia coefficients of the in-line force, were evaluated for each pitch ratio in the range of Keulegan-Carpenter numbers examined and are presented in diagrams.     

Analyse thermomécanique des problèmes de fissure fixe sous chargement dynamique

The purpose of this work is the study of thermomechanical coupling in the mechanisms of dynamic rupture in a fixed crack under dynamic loading.     

Numerical analysis of the stress concentration near holes originating in previously loaded viscoelastic bodies at finite strains

The state of stress and strain of previously loaded viscoelastic bodies with holes originating in them, successively or simultaneously, is analyzed under finite plane deformations. The problem statement and solution are based on the theory of repeatedly superimposed large deformations. The material mechanical properties are described using integral relations of the convolution type over time with a weakly singular kernel. The problem solving is based on the finite-element method. To calculate the integral of the convolution type, a recurrence formula is used that can be obtained by approximating the initial kernel with a linear combination of exponential functions (the truncated Prony’s series). The nonlinear effects and the effect of the interaction between holes on the stress concentration are analyzed. For the dynamic problems, the results for incompressible and weakly compressible materials are compared.     

Effects of thermal relaxations on thermoelastic interactions in an infinite orthotropic elastic medium with a cylindrical cavity

Thermoelastic interactions in an infinite orthotropic elastic medium with a cylindrical cavity are studied. The cavity surface is subjected to ramp-type heating of its internal boundary, which is assumed to be traction free. Lord–Shulman and Green–Lindsay models for the generalized thermoelasticity theories are selected since they allow for second-sound effects and reduce to the classical model for an appropriate choice of the parameters. The temperature, radial displacement, radial stress, and hoop stress distributions are computed numerically using the finite-element method (FEM). The results are presented graphically for different values of the thermal relaxation times using the three different theories of generalized thermoelasticity. Excellent agreement is found between the finite-element analysis and analytical and classical solutions.     

Modelling a recirculating density-driven turbulent flow

A mathematical model of turbulent density-driven flows is presented and is solved numerically. A form of the k–? turbulence model is used to characterize the turbulent transport, and both this non-linear model and a sediment transport equation are coupled with the mean-flow fluid motion equations. A partitioned, Newton–Raphson-based solution scheme is used to effect a solution. The model is applied to the study of flow through a circular secondary sedimentation basin.     

应用于液压传感的光子晶体光纤特性

为实现结构紧凑、高灵敏度的光纤压力(液压)传感器,提出了一种应用于液压传感的边孔结构光子晶体光纤.基于全矢量有限元方法,研究了传统光子晶体光纤和边孔结构光子晶体光纤的有效折射、模式等特性以及在液压情况下的应力和应力特性.根据光弹效应给出了传统光子晶体光纤和边孔结构光子晶体光纤在液压情况下的折射率变化特性.模拟结果表明边孔结构光子晶体光纤可以获得更大的液压传感灵敏度,增大边孔半径可以提高液压传感灵敏度,因此结构优化的边孔结构光子晶体光纤可以实现高灵敏度的光纤压力(液压)压力传感器.     

A simple four-node quadrilateral finite element for plates

A simple and accurate four-node quadrilateral finite element based on the Mindlin plate theory and Kirchhoff constraints is presented for general thin plate bending applications. The derivation of the element stiffness properties is straightforward, starting with a specified eight-node interpolation; usual discrete Kirchhoff (DK) constraints are employed to constrain out the four midside nodes of the element. The present resulting DK element passes patch tests with elements of arbitrary and even highly distorted mesh types. Numerical studies of the element convergence behaviours are undertaken for various plate bending problems so far investigated. It is indicated from comparative examples that fairly good convergence characteristics have been achieved.     

On computational proofs of the existence of solutions to nonlinear parabolic problems

This paper is an extension of the preceding study (Nakao, this journal, 1991) in which we described a numerical verification method of the solution for one-space dimensional parabolic problems, to the several-space dimensional case. Here, numerical verification means the automatic proof of the existence of solutions to the problems by some numerical techniques on a computer. We reformulate the verification condition for nonlinear parabolic initial boundary value problems using the fixed-point problem of a compact operator on certain function spaces. As in the preceding study based upon a simple C⁰ finite-element approximation and its constructive a priori error estimates, a numerical verification procedure is presented with some numerical examples.     

Numerical analysis of a moving boundary problem in coastal hydrodynamics

A finite element model to tackle the moving boundary problem of wave run-up on moderately steep slopes is developed. The special aspects considered in this study are (1) the modification of shallow water equations to accommodate the effect of vertical accelerations and (2) the use of Lagrangian acceleration coupled with an element that adapts itself to the moving boundary closely. The pressure term in the one-dimensional momentum equation is derived using the Eulerian equation in the vertical direction. This takes care of the vertical accelerations which are significant during the motion of a wave on moderately steep slopes. The element near the boundary is allowed to change its dimension so that the fluid boundary is closely followed. Such a flexible element precludes the need for approximation of the variables with regard to the indefinite position of the boundary. This element is split into two when its dimension becomes unduly large compared to the unchanging elements. The need for such a splitting is shown by an examination of the entries in the global matrix. Results of water profile as a wave runs up a structure are given. A brief history of the work on similar problems is outlined.     

基于小圆孔结构纤芯的高双折射光子晶体光纤

为了实现高双折射光子晶体光纤,提出了一种在纤芯中引入微小圆孔的方法.利用全矢量有限元方法和完美匹配层条件研究了基于圆孔微细结构纤芯的光子晶体光纤的双折射特性.讨论了纤芯圆孔数量、孔径、间隔距离对光纤双折射特性的影响;设计了一种双折射达到10－2量级的光子晶体光纤.模拟结果表明采用三个以上圆孔可以获得较大的双折射,增大外包层数目可以有效减小约束损耗.