Numerical Methods in Engineering with Python

"Zienkiewicz, Prof. Olgierd Cecil, (18 May 1921–2 Jan. 2009), Professor and Head of Civil Engineering Department, 1961–88, and Director, Institute for Numerical Methods in Engineering, 1976–88, University of Wales at Swansea, then Professor Emeritus; UNESCO Professor of Numerical Methods in Engineering, Universidad Politécnica de Cataluña, Barcelona; J. Walter Professor of Engineering, University of Texas, Austin, 1988–98" published on by Oxford University Press.

It is known that problems arise over existence and uniqueness of solution in quasi–static contact problems involving large coefficients of Coulomb friction. This paper considers the behaviour of an elastically supported mass with three translational degrees of freedom that can make contact with a rigid Coulomb friction support. A critical coefficient of friction is identified, above which the quasi–static solution can be non–unique. A numerical solution of the corresponding dynamic problem shows that the state realized then depends on the initial conditions. Even below the critical coefficient of friction, the dynamic solution can deviate from the quasi–static even for arbitrarily small loading rates. Typical dynamic responses include limit–cycle oscillations in velocity, oscillations involving a brief period of zero velocity (stick) and ‘stick–slip’ motion in which the mass spends significant periods in a state of stick, interspersed with short rapid–slip periods. All these non–steady motions involve non–rectilinear motion, even in cases where the time derivative of the applied load is constant in direction. A perturbation analysis is performed on the quasi–static frictional slip solution and predicts instability for certain slip directions, depending on functions of the off–diagonal stiffness components and the coefficient of friction. These results are presented in dimensionless form in stability diagrams. It is also shown that quasi–static slip that is stable when there is no damping can be destabilized if the damping matrix has sufficiently large off–diagonal components.

Free Access Other titles from ISTE in Numerical Methods in Engineering Alain Dervieux, Alain DervieuxSearch for more papers by this authorFrédéric Alauzet, Frédéric AlauzetSearch for more papers by this authorAdrien Loseille, Adrien LoseilleSearch for more papers by this authorBruno Koobus, Bruno KoobusSearch for more papers by this author Book Author(s):Alain Dervieux, Alain DervieuxSearch for more papers by this authorFrédéric Alauzet, Frédéric AlauzetSearch for more papers by this authorAdrien Loseille, Adrien LoseilleSearch for more papers by this authorBruno Koobus, Bruno KoobusSearch for more papers by this author First published: 15 August 2022 https://doi.org/10.1002/9781394164028.scard AboutPDFPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShareShare a linkShare onFacebookTwitterLinked InRedditWechat Mesh Adaptation for Computational Fluid Dynamics 2: Unsteady and Goal-oriented Adaptation RelatedInformation

Communications in Numerical Methods in Engineering, founded by Roland W. Lewis in 1985, will change its title to the ‘International Journal for Numerical Methods in Biomedical Engineering’ and have a revised Aims and Scope. Copyright © 2009 John Wiley & Sons, Ltd.

One of the fundamental problems in simulating the motion of sharp interfaces between immiscible fluids is a description of the transition that occurs when the interfaces merge and reconnect. It is well known that classical methods involving sharp interfaces fail to describe this type of phenomena. Following some previous work in this area, we suggest a physically motivated regularization of the Euler equations which allows topological transitions to occur smoothly. In this model, the sharp interface is replaced by a narrow transition layer across which the fluids may mix. The model describes a flow of a binary mixture, and the internal structure of the interface is determined by both diffusion and motion. An advantage of our regularization is that it automatically yields a continuous description of surface tension, which can play an important role in topological transitions. An additional scalar field is introduced to describe the concentration of one of the fluid components and the resulting system of equations couples the Euler (or Navier–Stokes) and the Cahn–Hilliard equations. The model takes into account weak non–locality (dispersion) associated with an internal length scale and localized dissipation due to mixing. The non–locality introduces a dimensional surface energy; dissipation is added to handle the loss of regularity of solutions to the sharp interface equations and to provide a mechanism for topological changes. In particular, we study a non–trivial limit when both components are incompressible, the pressure is kinematic but the velocity field is non–solenoidal (quasi–incompressibility). To demonstrate the effects of quasi–incompressibility, we analyse the linear stage of spinodal decomposition in one dimension. We show that when the densities of the fluids are not perfectly matched, the evolution of the concentration field causes fluid motion even if the fluids are inviscid. In the limit of infinitely thin and well–separated interfacial layers, an appropriately scaled quasi–incompressible Euler–Cahn–Hilliard system converges to the classical sharp interface model. In order to investigate the behaviour of the model outside the range of parameters where the sharp interface approximation is sufficient, we consider a simple example of a change of topology and show that the model permits the transition to occur without an associated singularity.

Two higher-order compact finite difference approaches, the Hermitian and the Lax-Wendroff are examined by applying them to the viscous Burgers' equation. The difference equations obtained by the two methods were integrated in time through the third-order Strong-Stability-Preserving Runge-Kutta scheme. Absolute errors are computed by using an exact solution. The results are also compared with a second order central difference scheme. The Hermitian approach is far easier to implement. On uniform grids the Lax-Wendroff scheme produces smaller errors during the initial stages, but both methods are equally good for larger durations of integration. The convergence rate of the Hermitian scheme is slightly higher than the Lax-Wendroff scheme although both are of fourth order. Both schemes are unstable beyond a certain step size in time and space. When numerical boundary conditions are imposed, second-order conditions produce the best results whereas linear extrapolation proved to be the worst. It was also observed that large domains were required to implement the numerical boundary conditions properly. There was no detrimental effect on the accuracy of the results obtained through either of the two schemes when the size of the domain was greatly increased. Both schemes showed remarkable improvement in accuracy when clustered grids were employed. However much smaller time steps are required for stable solutions.

The International Journal for Numerical Methods in Engineering keeps its first place in impact factor ranking for 2008 with the historic high impact factor of 2.229, as announced in the 2008 Journal Citation Reports® (Thomson Reuters, 2009). Copyright © 2009 John Wiley & Sons, Ltd.

This work presents a new technique for joining dissimilar meshes in non‐linear finite element analysis, with the novel result that the adjoining surfaces do not need to be spatially coincident for patch tests to be passed. The work builds in some sense on results previously presented by Dohrmann et al. (International Journal for Numerical Methods in Engineering 1991; 40:1205), but considerably generalizes upon the results presented there. Most importantly, it is shown that consistency of the tied finite element formulation may be achieved such that implementation of element formulations is not affected, through the use of a three‐field variational principle applied over the contacting region. A key conclusion of the work is that although the approach advocated bears some resemblance to mortar element schemes for contact analysis (see, for example, McDevitt and Laursen, International Journal for Numerical Methods in Engineering 2000; 48:1525), a straightforward application of a mortaring scheme is not sufficient to guarantee patch test passage when dissimilar meshing of curved surfaces leads to local open volumes (or overlaps) between tied surfaces in their reference configurations. The proposed algorithm is demonstrated through a number of numerical examples performed within a matrix‐free non‐linear equation solving framework. Published in 2003 by John Wiley & Sons, Ltd.

On the use of a reduced enhanced solid-shell (RESS) element for sheet forming simulations

We present a new approach to simulating unsteady fluid flows, with only very few degrees of freedom, by employing directly eigenmodes extracted from digital particle image velocimetry experimental data. In particular, we formulate standard Galerkin and nonlinear Galerkin approximations of the incompressible Navier–Stokes equations using hierarchical empirical eigenfunctions extracted from an ensemble of flow snapshots. We demonstrate that standard Galerkin approaches produce simulations capable of capturing the short–term dynamics of the flow, but nonlinear Galerkin projections are more effective in capturing both the short– and long–term dynamics, leading to bounded solutions. These findings are documented by applying these approaches to flow past a stationary circular cylinder at Reynolds number 610.

The third edition of this successful text describes and evaluates a range of widely used numerical methods, with an emphasis on problem solving. Every method is discussed thoroughly and illustrated with problems involving both hand computation and programming. MATLAB® M-files accompany each method and are available on the book's web page. Code is made simple and easy to understand by avoiding complex book-keeping schemes, while maintaining the essential features of the method. The third edition features a new chapter on Euler's method, a number of new and improved examples and exercises, and programs which appear as function M-files. Numerical Methods in Engineering with MATLAB®, 3rd edition is a useful resource for both graduate students and practicing engineers.

The third edition of this successful text describes and evaluates a range of widely used numerical methods, with an emphasis on problem solving. Every method is discussed thoroughly and illustrated with problems involving both hand computation and programming. MATLAB® M-files accompany each method and are available on the book's web page. Code is made simple and easy to understand by avoiding complex book-keeping schemes, while maintaining the essential features of the method. The third edition features a new chapter on Euler's method, a number of new and improved examples and exercises, and programs which appear as function M-files. Numerical Methods in Engineering with MATLAB®, 3rd edition is a useful resource for both graduate students and practicing engineers.

including an appendix by F

International Conference on Steel and Concrete StructuresBuilt Heritage Management and PresentationCongress on Numerical Methods in Engineering

Numerical methods in engineering and applied science : B.M. Irons and N.G. Shrive: Ellis Horwood, Chichester, United Kingdom, 1987, 248 pages, UK £10.95 (paper), UK £29.95/U.S. $54.20 (Cloth), ISBN 0-7458-0101-3.