On streamline diffusion arising in Galerkin FEM with predictor/multi‐corrector time integration

In the present paper, the author shows that the predictor/multi‐corrector (PMC) time integration for the advection–diffusion equations induces numerical diffusivity acting only in the streamline direction, even though the equations are spatially discretized by the conventional Galerkin finite element method (GFEM). The transient 2‐D and 3‐D advection problems are solved with the PMC scheme using both the GFEM and the streamline upwind/Petrov Galerkin (SUPG) as the spatial discretization methods for comparison. The solutions of the SUPG‐PMC turned out to be overly diffusive due to the additional PMC streamline diffusion, while the solutions of the GFEM‐PMC were comparatively accurate without significant damping and phase error. A similar tendency was seen also in the quasi‐steady solutions to the incompressible viscous flow problems: 2‐D driven cavity flow and natural convection in a square cavity. Copyright © 2002 John Wiley & Sons, Ltd.     

Cohomological aspects of gauge theories: superfield formalism

Some of the key cohomological features of the two (1 + 1)-dimensional (2D) free Abelian- and self-interacting non-Abelian gauge theories (having no interaction with matter fields) are briefly discussed first in the language of symmetry properties of the Lagrangian densities and the same issues are subsequently addressed in the framework of superfield formulation on the four (2 + 2)-dimensional supermanifold. Special emphasis is laid on the on-shell- and off-shell nilpotent (co-)BRST symmetries that emerge after the application of (dual) horizontality conditions on the supermanifold. The (anti-)chiral superfields play a very decisive role in the derivation of the on-shell nilpotent symmetries. The study of the present superfield formulation leads to the derivation of some new symmetries for the Lagrangian density and the symmetric energy-momentum tensor. The topological nature of the above theories is captured in the framework of superfield formulation and the geometrical interpretations are provided for some of the topologically interesting quantities.     

On stress rate and plasticity constitutive equations referred to a body-fixed coordinate system

Investigators have viewed the stress rate in two different ways: the material (body-fixed) point of view and the Eulerian point of view. We discuss the Zaremba–Jaumann rate and Oldroyd’s rate from the material viewpoint and apply them to the material formulation of a theory of plasticity for materials undergoing anisotropic plastic deformation. Significant advantages of the material formulation are that the derivation of equations is straight forward, the distortion of yield surface can be easily accounted for, and the issue of self-consistent elastic equation does not arise.     

Error analysis of the L2 least‐squares finite element method for incompressible inviscid rotational flows

In this article we analyze the L² least‐squares finite element approximations to the incompressible inviscid rotational flow problem, which is recast into the velocity‐vorticity‐pressure formulation. The least‐squares functional is defined in terms of the sum of the squared L² norms of the residual equations over a suitable product function space. We first derive a coercivity type a priori estimate for the first‐order system problem that will play the crucial role in the error analysis. We then show that the method exhibits an optimal rate of convergence in the H¹ norm for velocity and pressure and a suboptimal rate of convergence in the L² norm for vorticity. A numerical example in two dimensions is presented, which confirms the theoretical error estimates. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Navier-Stokes equations with imposed pressure and velocity fluxes

Boundary value problems for Stokes and Navier-Stokes equations with non-standard boundary conditions are studied. Included is the case where the pressure or its normal derivative is given on some part of the boundary or the pressure is given up to a constant but given velocity flux. First, a variational formulation is introduced which is shown to be equivalent to the Stokes equations with the non-standard boundary conditions under consideration. The existence and uniqueness of the solution of the variational problem are studied. Secondly, most of the results obtained for the Stokes equations are extended to the case of the Navier-Stokes equations. The final section is devoted to numerical experiments, flows in pipes and physiological flows.     

Solvability of the fractional hyperbolic Keller–Segel system

We study a new nonlocal approach to the mathematical modelling of the chemotaxis problem, which describes the random motion of a certain population due to a substance concentration. Considering the initial–boundary value problem for the fractional hyperbolic Keller–Segel model, we prove the solvability of the problem. The solvability result relies mostly on fractional calculus and kinetic formulation of scalar conservation laws.     

On Non-Symmetric Deformations of Neo-Hookean Solids

Deformations possible (i.e., those satisfying the governing three-dimensional equations of equilibrium and the incompressibility constraint) within a class of non-symmetric deformations for a neo-Hookean nonlinearly elastic body were determined in [1], where it was found that only three special cases of the class of deformation fields considered could be solutions. One of these is the trivial solution, one the solution describing radially symmetric deformation, and the other a (non-symmetric, non-homogeneous) deformation contained within a family of universal deformations. In this paper, the results reported in [1] are shown to hold for a substantially broadened deformation field. This revised version was published online in July 2006 with corrections to the Cover Date.     

Global structural behaviour of thin and moderately thick monoclinic spherical shells using a mixed shear deformation model

Summary The static and dynamic responses of anisotropic spherical shells under a uniformly distributed transverse load are investigated. Analytical solutions using the mixed variational formulation are presented for spherical shells subjected to various boundary conditions. Numerical results of a refined mixed first-order shear deformation theory for natural frequencies, critical buckling, center deflections and stresses are compared with those obtained using the classical shell theory. A variety of simply-supported and clamped boundary conditions are considered and comparisons with the existing literature are made. The sample numerical results presented herein for global structural behaviour of monoclinic spherical shells should serve as references for future comparisons.     

Numerical solution of Navier–Stokes equations using multiquadric radial basis function networks

A numerical method based on radial basis function networks (RBFNs) for solving steady incompressible viscous flow problems (including Boussinesq materials) is presented in this paper. The method uses a ‘universal approximator’ based on neural network methodology to represent the solutions. The method is easy to implement and does not require any kind of ‘finite element‐type’ discretization of the domain and its boundary. Instead, two sets of random points distributed throughout the domain and on the boundary are required. The first set defines the centres of the RBFNs and the second defines the collocation points. The two sets of points can be different; however, experience shows that if the two sets are the same better results are obtained. In this work the two sets are identical and hence commonly referred to as the set of centres. Planar Poiseuille, driven cavity and natural convection flows are simulated to verify the method. The numerical solutions obtained using only relatively low densities of centres are in good agreement with analytical and benchmark solutions available in the literature. With uniformly distributed centres, the method achieves Reynolds number Re = 100 000 for the Poiseuille flow (assuming that laminar flow can be maintained) using the density of , Re = 400 for the driven cavity flow with a density of and Rayleigh number Ra = 1 000 000 for the natural convection flow with a density of . Copyright © 2001 John Wiley & Sons, Ltd.     

Interlaminar hygrothermal stresses in laminated plates

Within the elasticity formulation the most general displacement field for hygrothermal problems of long laminated composite plates is presented. The equivalent single-layer theories are then employed to determine the global deformation parameters appearing in the displacement fields of general cross-ply, symmetric, and antisymmetric angle-ply laminates under thermal and hygroscopic loadings. Reddy’s layerwise theory is subsequently used to determine the local deformation parameters of various displacement fields. An elasticity solution is also developed in order to validate the efficiency and accuracy of the layerwise theory in predicting the interlaminar normal and shear stress distributions. Finally, various numerical results are presented for edge-effect problems of several cross-ply, symmetric, and antisymmetric angle-ply laminates subjected to uniform hygrothermal loads. All results indicate high stress gradients of interlaminar normal and shear stresses near the edges of laminates.