On the use of some numerical damping methods of spurious oscillations in the case of elastic wave propagation

The use of finite element and finite difference methods of spatial and temporal discretization for solving structural dynamics problems gives rise to purely numerical errors. Among the many numerical methods used to damp out the spurious oscillations occurring in the high frequency domain, it is proposed here to analyse and compare the well-known Bulk Viscosity method, which modifies the stresses calculations and a method recently presented by Tchamwa and Wielgosz, which is based on a modification of an explicit time integration algorithm. The efficiency of both methods is evaluated in a 2-D axisymmetric compressive test.     

Dynamics and wave propagation in dilatant granular materials

The equations of motion for dilatant granular material are obtained from a Hamiltonian variational principle of local type in the conservative case. The propagation of nonlinear waves in a region with uniform state is studied by means of an asymptotic approach that has already appeared useful in an investigation on wave propagation in bubbly liquids and in fluid mixtures. When the grains are assumed to be incompressible, it is shown that the material behaves as a continuum with latent microstructure.

---

Sommario Si ricavano le equazioni di moto per i materiali granulari dilatanti da un principio variazionale Hamiltoniano di tipo locale nel caso conservativo. Si studia la propagazione delle onde non lineari in una regione di stato costante per mezzo di un approccio asintotico già rivelatosi utile nello studio della propagazione di onde nei liquidi con bolle e nelle miscele di fluidi. Quando si supponga che i granuli siano incomprimibili, si dimostra che il materiale si comporta come un continuo con microstruttura latente.

     

Numerical study of wave processes in a pressure-wave refrigerator

The paper provides some results concerning the numerical study of the strongly transient gasdynamic processes in a pressure-wave refrigerator (PWR). A hierarchical set of numerical models from the simplest one-dimensional to a fully three-dimensional formulation is introduced. The computations show that one-dimensional solutions give a reasonable foundation for the understanding of PWR operational principles but cannot satisfactorily predict the refrigeration efficiency. Good agreement with experiment is achieved by considering two- and three-dimensional effects of gas mixing by overlap of rotating nozzles and expansion tubes. Received Received 14 February 1996 / Accepted 1 June 1996     

The study of wave motion in the Talbot interferometer with a lens

We present the study of the wave motion in the Talbot interferometer with an additional element such as a lens for all related audiences. Our solutions are in the analytic form. A general principle of the Talbot effect, which is the optical near-field effect, is the Fresnel diffraction. The Fresnel integral is rather complicated. We therefore introduce an alternative method which is based on the wave propagation through the transmission functions of the grating and the lens. Our method has been proved by a simple experimental setup.     

Timoshenko beam theory: A perspective based on the wave-mechanics approach

This paper reviews Timoshenko beam theory from the point of view of wave mechanics. Vibration of beam structures can be studied in terms of either normal modes or propagating waves. The latter wave approach has two distinct features: first, it gives rise to clear physical understanding of beam vibration; second, it leads to exact methods for vibration analysis of beam structures, especially in the mid-frequency range. In this paper, the work on wave solutions of an infinite Timoshenko beam is first discussed. The work on the splitting effect of spinning on wave solutions is also reviewed. The wave is treated as constitutive components of standing waves (i.e. normal modes), and a discussion on how the wave components formulate various standing waves is presented. Finally, several numerical examples are presented to illustrate the pros and cons of using different wave approaches to tackle vibration analysis of finite-length Timoshenko beams.     

Steady-state countercurrent flow in one dimension

Two phase countercurrent steady-state flow through permeable media in one dimension is discussed. For steady-state countercurrent flow in water wet porous media, a saturation profile is predicted with the water saturation decreasing in the direction that the water phase is flowing. The de la Cruz and Spanos equations predict that the Muskat relative permeability curves for countercurrent flow will be less than the Muskat relative permeability curves for steady-state cocurrent flow. This result has immediate implications regarding the use of external drive techniques to determine relative permeabilities based on the Buckley-Leverett theory and Muskat's equations. These equations and current experimental evidence involving countercurrent flow indicate that Muskat's equations do not adequately describe the multiphase flow of immiscible fluids.     

On the use of the reach-back characteristics method for calculation of dispersion

The Holly-Preissmann two-point finite difference scheme (HP method) has been popularly used for solving the advection equation. The key idea of this scheme is to solve the dependent variable (i.e. the concentration for the pollutant transport problem) by the method of characteristics with the use of cubic interpolation on the spatial axis. The interpolating polynomials of higher order are constructed by use of the dependent variable and its derivatives at two adjacent grid points. In this paper a new interpolating technique is introduced for incorporation with the Holly-Preissmann two-point method. The new method is denoted herein as the Holly-Preissmann reach-back method (HPRB) and allows the characteristics to project back several time steps beyond the present time level. Through stability analyses it has been observed that the increase of the reach-back time step numbers for the characteristics indeed reduces the numerical damping and dispersive phenomena. A schematic model has been constructed to demonstrate the merits of this new technique for the calculation of the pure advection and dispersion equations. Numerical experiments and comparisons with analytical solutions which support and demonstrate this new technique are presented.     

CFD modelling of the cross-flow through normal triangular tube arrays with one tube undergoing forced vibrations or fluidelastic instability

A CFD methodology involving structure motion and dynamic re-meshing has been optimized and applied to simulate the unsteady flow through normal triangular cylinder arrays with one single tube undergoing either forced oscillations or self-excited oscillations due to damping-controlled fluidelastic instability. The procedure is based on 2D URANS computations with a commercial CFD code, complemented with user defined functions to incorporate the motion of the vibrating tube. The simulation procedure was applied to several configurations with experimental data available in the literature in order to contrast predictions at different calculation levels. This included static conditions (pressure distribution), forced vibrations (lift delay relative to tube motion) and self-excited vibrations (critical velocity for fluidelastic instability). Besides, the simulation methodology was used to analyze the propagation of perturbations along the cross-flow and, finally, to explore the effect on the critical velocity of the Reynolds number, the pitch-to-diameter ratio and the degrees of freedom of the vibrating cylinder.     

Dispersive waves in microstructured solids

The wave motion in micromorphic microstructured solids is studied. The mathematical model is based on ideas of Mindlin and governing equations are derived by making use of the Euler–Lagrange formalism. The same result is obtained by means of the internal variables approach. Actually such a model describes internal fields in microstructured solids under external loading and the interaction of these fields results in various physical effects. The emphasis of the paper is on dispersion analysis and wave profiles generated by initial or boundary conditions in a one-dimensional case.     

Analysis of micro-structural effects on phononic waves in layered elastic media with periodic nonsmooth coordinates

Propagation of elastic phononic waves in layered composite materials is analyzed by introducing nonsmooth periodic coordinates associated with structural specifics of the materials. Spatial scales of the original (smooth) coordinates are estimated by the wave lengths. In terms of the new coordinates, the homogenization procedure occurs naturally from the continuity conditions imposed on elastic displacements and forces at layer interfaces. As a result, higher-order asymptotic approximations describing spatiotemporal ‘macro’- and ‘micro’-effects of wave propagation are obtained in closed form. Such solutions provide visualizations for the wave shapes illustrating their structure induced local details. In particular, beat-wise mode shapes and effective anisotropy of acoustic wave propagation are revealed. The subharmonic beating in wave modes occur when wave lengths orthogonal to layers is about to ‘resonate’ with layer’ thickness. If the wave speed has a non-zero projection along the layers, then phase shifts between the beats are observed in different cross sections perpendicular to the layers.     

Wave propagation and dispersion in elasto-plastic microstructured materials

A Mindlin continuum model that incorporates both a dependence upon the microstructure and inelastic (nonlinear) behavior is used to study dispersive effects in elasto-plastic microstructured materials. A one-dimensional equation of motion of such material systems is derived based on a combination of the Mindlin microcontinuum model and a hardening model both at the macroscopic and microscopic level. The dispersion relation of propagating waves is established and compared to the classical linear elastic and gradient-dependent solutions. It is shown that the observed wave dispersion is the result of introducing microstructural effects and material inelasticity. The introduction of an internal characteristic length scale regularizes the ill-posedness of the set of partial differential equations governing the wave propagation. The phase speed does not necessarily become imaginary at the onset of plastic softening, as it is the case in classical continuum models and the dispersive character of such models constrains strain softening regions to localize.     

Wave propagation in carbon nanotubes embedded in an elastic matrix

This paper reports the results of an investigation into the characteristics of wave propagation in carbon nanotubes embedded in an elastic matrix, based on an exact shell model. Each of the concentric tubes of multi-walled carbon nanotubes is considered as an individual elastic shell and coupled together through the van der Waals forces between two adjacent tubes. The matrix surrounding carbon nanotubes is described as a spring element defined by the Winkler model. The effects of rotatory inertia and elastic matrix on the wave velocity, the critical frequency, and the amplitude ratio between two adjacent tubes are described and discussed through numerical examples. The results obtained show that wave propagation in carbon nanotubes may appear in a critical frequency at which the wave velocity changes suddenly; the elastic matrix surrounding carbon nanotubes debases the critical frequency and the wave velocity, and changes the vibration modes between two adjacent tubes; the rotatory inertia based on an exact shell model obviously influences the wave velocity at some wave modes. Finally, a comparison of dispersion solutions from different shell models is given. The present work may serve as a useful reference for the application and the design of nano-electronic and nano-drive devices, nano-oscillators, and nano-sensors, in which carbon nanotubes act as basic elements.     

FLEXURAL WAVE PROPAGATION IN NARROW MINDLIN＇S PLATE

Appling Mindlin＇s theory of thick plates and Hamilton system to propagation of elastic waves under free boundary condition, a solution of the problem was given. Dispersion equations of propagation mode of strip plates were deduced from eigenfunction expansion method. It was compared with the dispersion relation that was gained through solution of thick plate theory proposed by Mindlin. Based on the two kinds of theories, the dispersion curves show great difference in the region of short waves, and the cutoff frequencies are higher in Hamiltonian systems. However, the dispersion curves are almost the same in the region of long waves.     

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle flows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained. Received 14 June 1996 / Accepted 19 October 1996     

Wave propagation in double helical rods

The propagation of waves in helical rods has been studied extensively. However, studying the wave propagation in double helical rods have received less attention although this can be useful in multiple fields of science and engineering. Obtaining an analytical model for a double helical rod is challenging since the curvature and tortuosity are not constant. Thus, resolving the wave behaviour analytically is nearly impossible. In this paper, wave propagation in a double helical rod will be studied using the wave and finite element method which is a technique that can be used to model homogeneous and periodic one and two dimensional structures based on the periodic structure theory. For modelling a double helical rod, the finite element model of a single turn is processed using Bloch waves. The dispersion curves and wavemodes are obtained and the similarities and differences of waves in helical and double helical rods are highlighted.     

Seismic gyrotropy as a result of rock-microstructure dissymmetry

Seismic investigations show that rocks exhibit not only anisotropic properties but also gyrotropic ones. Recently, the phenomenological theory has been developed to describe the phenomenon ‘seismic gyrotropy'. The paper gives a possible explanation of the gyrotropic properties of rocks on a microlevel. A gyrotropic model of grainy rock with dissymetry in its microstructure is presented. The model is built on a principle ‘azimuthal turn plus translation'. The mathematical and physical formulations of the problem are discussed. An algorithm is constructed to find the displacements, their first and second derivatives with respect to spatial coordinates as well as stresses inside a grain, an element of the model. Mathematical modelling of elastic-wave propagation in the model considered enabled to demonstrate that the model possesses gyrotropic properties, i.e., it rotates on a polarization vector. Gyration constants are determined.     

Primary wave transmission in systems of elastic rods with granular interfaces

We study analytically and numerically primary pulse transmission in one dimensional systems of identical linearly elastic non-dispersive rods separated by identical homogeneous granular layers composed of n beads. The beads interact elastically through a strongly (essentially) nonlinear Hertzian contact law. The main challenge in studying pulse transmission in such strongly nonlinear media is to analyze the ‘basic problem’, namely, the dynamical response of a single intermediate granular layer, confined from both ends by barely touching linear elastic rods subject to impulsive excitation of the left rod. The analysis of the basic problem is carried out under two basic assumptions; namely, of sufficiently small duration of the shock excitation applied to the first layer of the system, and of sufficiently small mass of each bead in the granular interface compared to the mass of each rod. In fact, the smallness of the mass of the bead defines the small parameter in the asymptotic analysis of this problem. Both assumptions are reasonable from the point of view of practical applications. In the analysis we focus only in primary pulse propagation, by neglecting secondary pulse reflections caused by wave scattering at each granular interface and considering only the transmission of the main (primary) pulse across the interface to the neighboring elastic rod. Two types of shock excitations are considered. The first corresponds to fixed time duration (but still much smaller compared to the characteristic time of pulse propagation through the length of each rod), whereas the second type corresponds to a pulse duration that depends on the small parameter of the problem. The influence of the number of beads of the granular interface on the primary wave transmission is studied, and it is shown that at granular interfaces with a relatively low number of beads fast time scale oscillations are excited with increasing amplitudes with increasing number of beads. For a larger number of beads, primary pulse transmission is by means of solitary wave trains resulting from the dispersion of the original shock pulse; in that case fast oscillations result due to interference phenomena caused by the scattering of the main pulse at the boundary of the interface. Considering a periodic system of rods we demonstrate significant reduction of the primary pulse when transmitted through a sequence of granular interfaces. This result highlights the efficacy of applying granular interfaces for passive shock mitigation in layered elastic media.     

On the determination of the Young’s modulus of thin films using indentation tests

The main difficulty with the characterization of thin coatings using depth-sensing indentation tests is related to the determination of the contributions of the substrate and the film to the measured properties. In this study, three-dimensional numerical simulations of the Vickers hardness test are used in order to examine the influence of the elastic and plastic properties of the substrate and the film on the composite’s Young’s modulus results. The hardness of the film is equal to or higher than the substrate hardness. A study of the stress distributions and the indentation geometry of composites, film/substrate, was performed, taking into account the relative mechanical properties of the film and substrate. In addition, stress evolution during indentation was studied, in order to quantify the critical indentation depth under which the substrate is not elastically deformed. The accurate evaluation of the Young’s modulus of the films using weight functions is also examined: some of these have previously been proposed and one was introduced for this study. Two different fitting procedures were used to compare the results obtained from eight fictive film/substrate combinations using six weight functions. The first procedure, commonly used, considers the substrate’s modulus as a known parameter in the fitting process. In the second, the film and the substrate’s modulus are considered as unknown variables that are calculated simultaneously during the fitting process. The validity of the conclusions obtained using the fictive materials was checked by applying the weight functions to four real composites.