CHARACTERISTICS OF SUBDIFFERENTIALS OF FUNCTIONS

ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ（郭兴明）ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ￥ＧｕｏＸｉｎｇｍｉｎｇ（ＲｅｃｅｉｖｅｄＪｕｎｅ１６，１９９５...     

CHARACTERIZATION OF POWDER FLOWABILITY USING MEASUREMENT OF ANGLE OF REPOSE

In response to the need in industry for a quick and reproducible method of measuring the flowability of powders in processes involving transport and storage, the apparently simple idea of measuring angle of repose has been revived. The simple tester first used has evolved over a number of years into the present version which is shown to be capable of handling quite small samples of powders that are even slightly cohesive. Experimental data are presented and results shown to compare consistently with the better known Hausner ratio.     

TORSIONAL RIGIDITY OF SHELLS OF REVOLUTION

In this paper,the general equations of equilibrium for axisymmetrical deformationincluding the torsional deformation of revolutional shells are derived.It is shown that theshearing stress distribution due to torsion is independent of other stress componentsincluding those of membrane stress and bending stress.In this paper,the torsional deformation is considered to be represented by membraneaction only,and also by the combined action of bending membrane deformation.It is shownthat the main contribution of torsional rigidity is that related to membrane action.     

THE THEORY OF SIMILARITY AND THE PROFILE OF VERTICAL DISTRIBUTION OF CONCENTRATION OF SAND PARTICLES

In this article,we treat the problem of two-dimensionaluniform steady channel flow of turbid water with theory ofsimilarity.Under the condition of similarity of turbulentfluctuation velocity and fluctuation of concentration ofsand particles,we obtain the profile of the vertical dis-tribution of concentration of sand particles.This profileof vertical distribution of concentration of sand particlesis slightly different from that obtained by diffusion theory,and departs from that obtained by gravitational theory.     

ANALYSIS OF THE MOTION OF A GYRO-THEODOLITE

With the method of analytical mechanics,this paper studies the motions of a gyro-theodolite under the action of (1) the torque of gravity only,(2)the torque applied by theband suspension.(3)the torque of the band suspension with air damping considered,theequations of motion are then established and their solutions are found.Furthermore,analysis of the law of motion and the behaviour of gyro-theodolite during the orientation ismade.     

THEORY OF ELASTICITY OF AN ANISOTROPIC BODY FOR THE BENDING OF BEAMS

A theory of elasticity for the bending of orthogonal anisotropic beams has been developed by analogy with the special case, which can be obtained by applying the theory of elasticity for bending of transversely isotropic plates to the problems of two deminsions. In this paper, we present a method to solve the problems of bending of orthogonal anisotropic beams and a new theory of the deep-beam whose ratio of depth to length is larger. It is pointed out that Reissner's theory to account for the effect of transverse shear deformation is not very approximate in the components of stress,     

APPLICATION 0F THE THEORY OF BRANCHING OF SOLUTIONS OF NONLINEAR EQUATIONS IN THE ESTIMATION OF POLE PERTURBANCE OF LINEAR CLOSED-LOOP SYSTEMS

When the state and input matrices of a multivariable lineartime-invariant system are perturbed.the problem of the es-timation of pole perturbance of closed-loop system is con-sidered by making use of the theory of branching of solutionsfor nonlinear equations.     

BLOW-UP OF SOLUTIONS OF NONLINEAR PSEUDO-HYPERBOLIC EQUATIONS OF GENERALIZED NERVE CONDUCTION TYPE

This paper deals with the two type of mixed problems with respect to Neumannboundary and Dirichlet boundary for nonlinear pseudo-hyperbolic equations of generalizednerve conduction typeu_1-△u_1=F(x,t,u,(?)u,u_t,(?)i_t)when the nonlinear part F(x.t,u,(?)u,u_t,(?)u_t)and the initial values satisfy someconditions,the blow-up properties of the solutions are obtained.     

COMPUTER SIMULATION OF THE MOTION OF THE BULLET BELT OF AIRPLANE GUN

ＣＯＭＰＵＴＥＲＳＩＭＵＬＡＴＩＯＮＯＦＴＨＥＭＯＴＩＯＮＯＦＴＨＥＢＵＬＬＥＴＢＥＬＴＯＦＡＩＲＰＬＡＮＥＧＵＮＺｈａｎｇＤｉｎｇｇｕｏ（章定国）（ＲｅｃｅｉｖｅｄＭａｒｃｈ１７，１９９５；ＣｏｍｍｕｎｉｃａｔｅｄｂｙＷｕＲｕｉｆｅｎｇ）Ａｂｓｔｒ...     

EFFECT OF POLYETHYLENEIMINE ON THE PROPERTIES OF CONCENTRATED SUSPENSIONS OF TITANIUM DIOXIDE MICROBEADS

The effect of polyethyleneimine (PEI) concentration on the properties of titanium dioxide (TiO2) suspensions is studied with or without the addition of an electrolyte (barium acetate). Measurements of the apparent viscosity and the stability of TiO2 suspensions showed that PEI is an effective dispersant for TiO2 particles in suspension in the absence of an electrolyte, not only reducing the viscosity of the suspension but also increasing its stability. In the presence of an electrolyte, however, small quantities of polyethyleneimine could neither disperse the TiO2 particles nor de-crease the viscosity of the TiO2 suspensions; only PEI concentrations beyond saturation adsorption could perceptively improve the stabilitv of TiO2 suspensions.     

Solution of the connection problems between a finite holed plate and a stiffener by using the partitioning concept of the generalized variational method

In this paper a new finite element method is presented, in which complex functions are chosen to be the finite element model and the partitioning concept of the generalized variational method is utilized. The stress concentration factors for a finite holed plate welded by a stiffener are calculated and the analytical solutions in series form are obtained. From some computer trials it is demonstrated that the problem of displacement compatibility and continuity of tractions between the holed plate and the stiffener is successfully analysed by using this method. Since only three elements need to be formulated, relatively less storage is required than the usual finite element methods. Furthermore, the accuracy of solutions is improved and the computer time requirements are considerably reduced. Numerical results of stress concentration factors and stresses along the welded-line which may be referential to engineers are shown in tables.     

An exact element method for plane problem

In this paper, based on the step reduction method^[1] and exact analytic method^[2] anew method-exact element method for constructing finite element, is presented. Since the new method doesn ’t need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a quadrilateral noncompatible element with 8 degrees of freedom is derived for the solution of plane problem. Since Jacobi ’s transformation is not applied, the present element may degenerate into a triangle element. It is convenient to use the element in engineering. In this paper, the convergence is proved. Numerical examples are given at the endof this paper, which indicate satisfactory results of stress and displacements can be obtained and have higher numerical precision in nodes.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.     

A Bubble Element for the Compressible Euler Equations

In this paper, we present a new Galerkin finite element method with bubble function for the compressible Euler equations. This method is derived from the scaled bubble element for the advection-diffusion problems developed by Simo and his colleagues, which is based on the equivalence between the Galerkin method employing piecewise linear interpolation with bubble functions and the Streamline-Upwind/Petrov Galerkin (SUPG) finite element method using P₁ approximation in the steady advection-diffusion problem. Simo and this author have applied this approach to transient advection-diffusion problems by using a special scaled bubble function called P-scaled bubble, which is designed to work in the transient advection-diffusion problems for any Peclet number from 0 to ∞. The method presented in this paper is an application of this p-scaled bubble element to a pure hyperbolic system.     

The dynamic stiffness matrix of the finite annular plate element

The dynamic deformation of harmonic vibration is used as the shape functions of the finite annular plate element, and sonic integration difficulties related to the Bessel’s functions are solved in this paper. Then the dynamic stiffness matrix of the finite annular plate element is established in closed form and checked by the direct stiffness method. The paper has given wide convcrage for decomposing the dynamic matrix into the power series of frequency square. By utilizing the axial symmetry of annular elements, the modes with different numbers of nodal diameters at s separately treated. Thus some terse and complete results are obtained as the foundation of structural characteristic analysis and dynamic response compulation.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

A finite element method for the two‐dimensional extended Boussinesq equations

A new numerical method for Nwogu's (ASCE Journal of Waterway, Port, Coastal and Ocean Engineering 1993; 119 :618)two‐dimensional extended Boussinesq equations is presented using a linear triangular finite element spatial discretization coupled with a sophisticated adaptive time integration package. The authors have previously presented a finite element method for the one‐dimensional form of these equations (M. Walkley and M. Berzins (International Journal for Numerical Methods in Fluids 1999; 29 (2):143)) and this paper describes the extension of these ideas to the two‐dimensional equations and the application of the method to complex geometries using unstructured triangular grids. Computational results are presented for two standard test problems and a realistic harbour model. Copyright © 2002 John Wiley & Sons, Ltd.     

A new stable space–time formulation for two‐dimensional and three‐dimensional incompressible viscous flow

A space–time finite element method for the incompressible Navier–Stokes equations in a bounded domain in ?^d (with d=2 or 3) is presented. The method is based on the time‐discontinuous Galerkin method with the use of simplex‐type meshes together with the requirement that the space–time finite element discretization for the velocity and the pressure satisfy the inf–sup stability condition of Brezzi and Babu?ka. The finite element discretization for the pressure consists of piecewise linear functions, while piecewise linear functions enriched with a bubble function are used for the velocity. The stability proof and numerical results for some two‐dimensional problems are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Two-level stabilized finite element method for Stokes eigenvalue problem

A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered.This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h = O(H 2),which can still maintain the asymptotically optimal accuracy.It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution,which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h.Hence,the two-level stabilized finite element method can save a large amount of computational time.Moreover,numerical tests confirm the theoretical results of the present method.     

Numerical simulation of transient planar flow of a viscoelastic material with two moving free surfaces

In this study, we examine the numerical simulation of transient viscoelastic flows with two moving free surfaces. A modified Galerkin finite element method is implemented to the two-dimensional non-steady motion of the fluid of the Oldroyd-B type. The fluid is initially placed between two parallel plates and bounded by two straight free boundaries. In this Lagrangian finite element method, the spatial mesh deforms in time along with the moving free boundaries. The unknown shape of the free surfaces is determined with the flow field u, v, τ, p by the deformable finite element method, combined with a predictor-corrector scheme in an uncoupled fashion. The moving free surfaces and fluid motion of both Newtonian and non-Newtonian flows are investigated. The results include the influence of surface tension, fluid inertia and elasticity.