Multi‐stage high order semi‐Lagrangian schemes for incompressible flows in Cartesian geometries

Efficient transport algorithms are essential to the numerical resolution of incompressible fluid‐flow problems. Semi‐Lagrangian methods are widely used in grid based methods to achieve this aim. The accuracy of the interpolation strategy then determines the properties of the scheme. We introduce a simple multi‐stage procedure, which can easily be used to increase the order of accuracy of a code based on multilinear interpolations. This approach is an extension of a corrective algorithm introduced by Dupont & Liu (2003, 2007). This multi‐stage procedure can be easily implemented in existing parallel codes using a domain decomposition strategy, as the communication pattern is identical to that of the multilinear scheme. We show how a combination of a forward and backward error correction can provide a third‐order accurate scheme, thus significantly reducing diffusive effects while retaining a non‐dispersive leading error term. Copyright © 2016 John Wiley & Sons, Ltd.     

The efect of the couple stress parameter and Prandtl number on the transient natural convection flow over a vertical cylinder

An analysis is performed to study transient free convective boundary layer flow of a couple stress fluid over a vertical cylinder, in the absence of body couples. The solution of the time-dependent non-linear and coupled governing equations is carried out with the aid of an unconditionally stable Crank–Nicolson type of numerical scheme. Numerical results for the steady-state velocity, temperature as well as the time histories of the skin-friction coefcient and Nusselt number are presented graphically and discussed. It is seen that for all flow variables as the couple stress control parameter, Co, is amplified, the time required for reaching the temporal maximum increases but the steady-state decreases.     

Viscoelastic response and interlaminar delamination resistance of epoxy/glass fiber/functionalized graphene oxide multi-scale composites

Graphene oxide (GO) was functionalized using three different diamines, namely ethylenediamine (EDA), 4,4′-diaminodiphenyl sulfone (DDS) and p-phenylenediamine (PPD) to reinforce an epoxy/glass fiber (EP/GF) composite laminate, with the aim of improving the overall composite mechanical performance. Different mechanical characterization techniques were used to determine the mechanical performance, including: tensile stress strain, double cantilever beam (DCB) mode-I fracture toughness and dynamic mechanical thermal analysis (DMTA). Scanning electron microscopy (SEM) was used to support the results and conclusions. The results demonstrated remarkable enhancements in the mechanical performance of EP/GF composite laminates by incorporation of functionalized graphene oxide (FGO) nanofiller, whilst the mechanical performance of the GO reinforced composite only improved marginally. Finally, the mechanical performance of the EP/GF/FGO multi-scale composites was found to be dependent on the type of FGO functional groups; of which EDA exhibited the highest performance. These investigations confirmed that the EDA-FGO-reinforced EP/GF composites possess excellent potential to be used as multifunctional engineering materials in industrial applications.     

Experimental and numerical investigation of aramid fibre reinforced laminates subjected to low velocity impact

The low velocity impact performance of domestic aramid fibre reinforced laminates is investigated experimentally and numerically. Laminates with different thicknesses are impacted by drop-weight test machine under different impact energies. The time histories of impact force are recorded and ultrasonic C-scan technology is used to inspect the internal damage of the laminates. Numerical simulation is conducted using finite element method (FEM), taking into account both intralaminar and interlaminar damage. The intralaminar damage model is based on the continuum damage mechanics (CDM) approach, which consists of the strain-based Hashin failure criteria and the exponential damage evolution law, and considers the nonlinear shear behaviour of the material. The interlaminar damage is simulated by interface elements with cohesive zone model. The numerical results show good agreements with the experiments, thus verifying the validity of the presented numerical model.     

FULLY SYMBOLIC GENERATION OF COMPLEX MULTIBODY MODELS*

Abstract

In this paper, a quasistatic finite element model of a spur gear pair is developed. A node-to-target contact formulation is given, where calculations of initial gaps are based upon the actual geometry of the gear flanks rather than upon a contact node and a facet or a line segment. By using a special contact search algorithm, profile modifications and mounting errors are easily incorporated in the analysis. The problem, which also includes friction, is solved by using a nonsmooth Newton method. The static transmission error can be calculated with accuracy with a relatively small number of nodes along the gear flanks. Several examples are given in order to demonstrate the model.     

A new formulation of regularized meshless method applied to interior and exterior anisotropic potential problems

This paper proposes a new formulation of regularized meshless method (RMM), which differs from the traditional RMM in that the traditional formulation generates the diagonal elements of influence matrix via null-field integral equations, while our new one directly employs the boundary integral equations at the domain point to evaluate the diagonal elements. We test the present RMM formulation to two-dimensional anisotropic potential problems in finite and infinite domains in comparison with the traditional RMM. Numerical results show that the present RMM sharply outperforms the traditional RMM in the solution of interior problems, while the latter is clearly superior for exterior problems. A rigorous theoretical analysis of circular domain case also corroborates such numerical experiment observations and is provided in the appendix of this paper.     

Mixed mode fracture of a piezoelectric–piezomagnetic bi-layer structure with two un-coaxial cracks parallel to the interface and each in a layer

The purpose of the present work is to study the mixed mode fracture of a piezoelectric–piezomagnetic composite with two un-coaxial cracks parallel to the interface and each in a layer. Methods of generalized dislocation simulation, Green’s function, Cauchy singular integral equation and Lobatto–Chebyshev collocation are combined together to get the numerical results of mechanical strain energy release rate (MSERR). Three kinds of effects are revealed by parametric studies, i.e., the free-surface effect, the shielding effect and the interference effect, and they are used to interpret the characteristics of COD and MSERR curves. In addition, the effects of shear loading, magnetic loading and electric loading on MSERR are also disclosed, respectively, by varying the corresponding loading factor.     

Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements

Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials' auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.     

Probabilistic approaches to compute uncertainty intervals and sensitivity factors of ultrasonic simulations of a weld inspection

For comprehension purpose, numerical computations are more and more used to simulate the propagation phenomena observed during experimental inspections. However, the good agreement between experimental and simulated data necessitates the use of accurate input data and thus a good characterization of the inspected material. Generally the input data are provided by experimental measurements and are consequently tainted with uncertainties. Thus, it becomes necessary to evaluate the impact of these uncertainties on the outputs of the numerical model. The aim of this study is to perform a probabilistic analysis of an ultrasonic inspection of an austenitic weld containing a manufactured defect based on advanced techniques such as polynomial chaos expansions and computation of sensitivity factors (Sobol, DGSM). The simulation of this configuration with the finite element code ATHENA2D was performed 6000 times with variations of the input parameters (the columnar grain orientation and the elastic constants of the material). The 6000 sets of input parameters were obtained from adapted statistical laws. The output parameters (the amplitude and the position of the defect echo) distributions were then analyzed and the 95% confidence intervals were determined.