STATISTICAL ENERGY ANALYSIS OF COUPLED PLATE SYSTEMS WITH LOW MODAL DENSITY AND LOW MODAL OVERLAP

Finite element methods, experimental statistical energy analysis (ESEA) and Monte Carlo methods have been used to determine coupling loss factors for use in statistical energy analysis (SEA). The aim was to use the concept of an ESEA ensemble to facilitate the use of SEA with plate subsystems that have low modal density and low modal overlap. An advantage of the ESEA ensemble approach was that when the matrix inversion failed for a single deterministic analysis, the majority of ensemble members did not encounter problems. Failure of the matrix inversion for a single deterministic analysis may incorrectly lead to the conclusion that SEA is not appropriate. However, when the majority of the ESEA ensemble members have positive coupling loss factors, this provides sufficient motivation to attempt an SEA model. The ensembles were created using the normal distribution to introduce variation into the plate dimensions. For plate systems with low modal density and low modal overlap, it was found that the resulting probability distribution function for the linear coupling loss factor could be considered as lognormal. This allowed statistical confidence limits to be determined for the coupling loss factor. The SEA permutation method was then used to calculate the expected range of the response using these confidence limits in the SEA matrix solution. For plate systems with low modal density and low modal overlap, relatively small variation/uncertainty in the physical properties caused large differences in the coupling parameters. For this reason, a single deterministic analysis is of minimal use. Therefore, the ability to determine both the ensemble average and the expected range with SEA is crucial in allowing a robust assessment of vibration transmission between plate systems with low modal density and low modal overlap.     

Extension of SEA model to subsystems with non-uniform modal energy distribution

In order to widen the application of statistical energy analysis (SEA), a reformulation is proposed. Contrary to classical SEA, the model described here, statistical modal energy distribution analysis (SmEdA), does not assume equipartition of modal energies.Theoretical derivations are based on dual modal formulation described in Maxit and Guyader (Journal of Sound and Vibration 239 (2001) 907) and Maxit (Ph.D. Thesis, Institut National des Sciences Appliquées de Lyon, France 2000) for the general case of coupled continuous elastic systems. Basic SEA relations describing the power flow exchanged between two oscillators are used to obtain modal energy equations. They permit modal energies of coupled subsystems to be determined from the knowledge of modes of uncoupled subsystems. The link between SEA and SmEdA is established and make it possible to mix the two approaches: SmEdA for subsystems where equipartition is not verified and SEA for other subsystems.Three typical configurations of structural couplings are described for which SmEdA improves energy prediction compared to SEA: (a) coupling of subsystems with low modal overlap, (b) coupling of heterogeneous subsystems, and (c) case of localized excitations.The application of the proposed method is not limited to theoretical structures, but could easily be applied to complex structures by using a finite element method (FEM). In this case, FEM are used to calculate the modes of each uncoupled subsystems; these data are then used in a second step to determine the modal coupling factors necessary for SmEdA to model the coupling.     

A novel hybrid ES-FE-SEA for mid-frequency prediction of Transmission losses in complex acoustic systems

The widely-used numerical modeling approaches such as the finite element method (FEM) and statistical energy analysis (SEA) often have limited applicability to the transmission loss prediction in mid-frequency range. In this paper, a novel hybrid edge-based smoothed FEM coupled with statistical energy analysis (ES-FE-SEA) method is proposed to further improve the accuracy of “mid-frequency” transmission loss predictions. The application of ES-FEM will “soften” the well-known ‘‘overly-stiff’’ behavior in the standard FEM solution and reduce the inherent numerical dispersion error. While the SEA approach deals with the physical uncertainty in the relatively higher frequency range. The plate of interest is appropriately described by an ES-FEM model, due to its relative robustness to perturbations. Its adjacent reverberation cavities are modeled by employing the SEA approach, because of their high model density. The coupling and interaction between SEA subsystems and the FE subsystem is governed by the “reciprocity relationship” theorem. A standard numerical example for benchmarking is examined and excellent agreement was achieved between the prediction and reference results. The proposed ES-FE-SEA is also verified by various numerical examples. The method is finally applied to the modeling a complicated engineering problem–acoustic fields on both sides of the front windshield in a passenger car.     

Analysis of the modal energy distribution of an excited vibrating panel coupled with a heavy fluid cavity by a dual modal formulation

This paper describes the modal interaction between a panel and a heavy fluid cavity when the panel is excited by a broad band force in a given frequency band. The dual modal formulation (DMF) allows describing the fluid–structure coupling using the modes of each uncoupled subsystem. After having studied the convergence of the modal series on a test case, we estimate the modal energies and the total energy of each subsystem. An analysis of modal energy distribution is performed. It allows us to study the validity of SEA assumptions for this case. Added mass and added stiffness effects of the fluid are observed. These effects are related to the non-resonant acoustic modes below and above the frequency band of excitation. Moreover, the role of the spatial coupling of the resonant cavity modes with the non-resonant structure modes is also highlighted. As a result, the energy transmitted between the structure and the heavy fluid cavity generally cannot be deduced from the SEA relation established for a light fluid cavity.     

Estimating ensemble average power delivered by a piezoelectric patch actuator to a non-deterministic subsystem

Engineering systems such as aircraft, ships and automotive are considered built-up structures. Dynamically they are taught of as being fabricated from many components that are classified as ‘deterministic subsystems’ (DS) and ‘non-deterministic subsystems’ (Non-DS). Structures' response of the DS is deterministic in nature and analysed using deterministic modelling methods such as finite element (FE) method. The response of Non-DS is statistical in nature and estimated using statistical modelling technique such as statistical energy analysis (SEA). SEA method uses power balance equation, in which any external input to the subsystem must be represented in terms of power. Often, input force is taken as point force and ensemble average power delivered by point force is already well-established. However, the external input can also be applied in the form of moments exerted by a piezoelectric (PZT) patch actuator. In order to be able to apply SEA method for input moments, a mathematical representation for moment generated by PZT patch in the form of average power is needed, which is attempted in this paper. A simply-supported plate with attached PZT patch is taken as a benchmark model. Analytical solution to estimate average power is derived using mobility approach. Ensemble average of power given by the PZT patch actuator to the benchmark model when subjected to structural uncertainties is also simulated using Lagrangian method and FEA software. The analytical estimation is compared with the Lagrangian model and FE method for validation. The effects of size and location of the PZT actuators on the power delivered to the plate are later investigated.     

Implications of heterogeneous inputs and connectivity on the synchronization in excitable networks

We study the combined implications of connectivity and heterogeneous inputs on the synchronization features of a one-dimensional chain of diffusively coupled FitzHugh Nagumo (FHN) systems. The uncoupled systems are triggered into a regime of chaotic firing by periodic parametric forces modeling external stimuli. Due to the parameter dispersion involved in randomly distributed amplitudes and/or phases of the forces the units are nonidentical and the firing events on the chain of uncoupled units will be asynchronous leading to a distribution of the spiking times. Interest is focused on mutually synchronized spikings arising through the coupling where the connectivity of the network may range from nearest-neighbor interaction to global interactions. From our studies we conclude that increasing the interaction radius does not necessarily entail better spike synchrony and the coupling strength plays a more important role than connectivity. It is found that for driving with random amplitudes together with random phases a critical interaction radius exists beyond which firing becomes suppressed if the coupling between the units is too strong. In such cases of ‘firing death’ the units perform only small-amplitude oscillations which are mutually synchronous. The optimal coupling for spike synchrony is of intermediate strength and altering the connectivity does not really matter for the degree of spike synchrony. Distinct to that, when all the phases are equal and only the amplitudes of the forces are randomly distributed enhanced spike synchrony is achieved for sufficiently strong coupling regardless of the interaction radius.     

Extension of the statistical modal energy distribution analysis for estimating energy density in coupled subsystems

The present article deals with an extension of the Statistical modal Energy distribution Analysis (SmEdA) method to estimate kinetic and potential energy density in coupled subsystems. The SmEdA method uses the modal bases of uncoupled subsystems and focuses on the modal energies rather than the global energies of subsystems such as SEA (Statistical Energy Analysis). This method permits extending SEA to subsystems with low modal overlap or to localized excitations as it does not assume the existence of modal energy equipartition. We demonstrate that by using the modal energies of subsystems computed by SmEdA, it is possible to estimate energy distribution in subsystems. This approach has the same advantages of standard SEA, as it uses very short calculations to analyze damping effects. The estimation of energy distribution from SmEdA is applied to an academic case and an industrial example.     

A theory of cosmic rays

We present a theory of non-solar cosmic rays (CRs) in which the bulk of their observed flux is due to a single type of CR source at all energies. The total luminosity of the Galaxy, the broken power-law spectra with their observed slopes, the position of the ‘knee(s)’ and ‘ankle’, and the CR composition and its variation with energy are all predicted in terms of very simple and completely ‘standard’ physics. The source of CRs is extremely ‘economical’: it has only one parameter to be fitted to the ensemble of all of the mentioned data. All other inputs are ‘priors’, that is, theoretical or observational items of information independent of the properties of the source of CRs, and chosen to lie in their pre-established ranges. The theory is part of a ‘unified view of high-energy astrophysics’ — based on the ‘Cannonball’ model of the relativistic ejecta of accreting black holes and neutron stars. The model has been extremely successful in predicting all the novel properties of Gamma Ray Bursts recently observed with the help of the Swift satellite. If correct, this model is only lacking a satisfactory theoretical understanding of the ‘cannon’ that emits the cannonballs in catastrophic processes of accretion.     

Validity diagrams of statistical energy analysis

This paper is concerned with the validity domain of statistical energy analysis (SEA) which is defined in terms of four criteria. The mode count N and the modal overlap M must be high, the normalized attenuation factor and the coupling strength γ must be small. The application of dimensional analysis on the governing equations of plates gives the space of dimensionless parameters in which the validity domain of SEA must be delimited. This domain is discussed on the basis of geometry of the surfaces delimiting it. The diagrams of validity of SEA are introduced and discussed. A numerical simulation on a couple of rectangular plates coupled along one edge illustrates the theoretical approach.     

Gravitational instabilities of isothermal spheres in the presence of a cosmological constant

Gravitational instabilities of isothermal spheres are studied in the presence of a positive or negative cosmological constant, in the Newtonian limit. In gravity, the statistical ensembles are not equivalent. We perform the analysis both in the microcanonical and the canonical ensembles, for which the corresponding instabilities are known as ‘gravothermal catastrophe’ and ‘isothermal collapse’, respectively. In the microcanonical ensemble, no equilibria can be found for radii larger than a critical value, which is increasing with increasing cosmological constant. In contrast, in the canonical ensemble, no equilibria can be found for radii smaller than a critical value, which is decreasing with increasing cosmological constant. For a positive cosmological constant, characteristic reentrant behavior is observed.     

Correlation of finite element models of multi-physics systems

The modal assurance criterion (MAC) and normalized cross-orthogonality (NCO) check are widely used to assess the correlation between the experimentally determined modes and the finite element model (FEM) predictions of mechanical systems. Here, their effectiveness in the correlation of FEM of two types of multi-physics systems, namely, viscoelastic damped systems and a shunted piezoelectric system are investigated using the dynamic characteristics obtained from a nominal FEM, that are considered as the ‘true’ or experimental characteristics and those obtained from the inaccurate FEMs. The usefulness of the MAC and NCO check in the prediction of the overall loss factor of the viscoelastic damped system, which is an important design tool for such systems, is assessed and it is observed that these correlation methods fail to properly predict the damping characteristics, along with the responses under base excitation. Hence, base force assurance criterion (BFAC) is applied by comparing the ‘true’ dynamic force at the base and inaccurate FEM predicted force such that the criterion can indicate the possible error in the acceleration and loss factor. The effect of temperature as an uncertainty on the MAC and NCO check is also studied using two viscoelastic systems. The usefulness of MAC for the correlation of a second multi-physics FEM that consists of a shunted piezoelectric damped system is also analyzed under harmonic excitation. It has been observed that MAC has limited use in the correlation and hence, a new correlation method – current assurance criterion – based on the electric current is introduced and it is demonstrated that this criterion correlates the dynamic characteristics of the piezoelectric system better than the MAC.     

Vibration characteristics of a cantilever plate with attached spring-mass system

Coupled free vibration analysis has been performed on a cantilever thin plate carrying a spring-mass system attached on an arbitrary point by using Rayleigh-Ritz method. Influence of an attached ‘spring-mass’ system, i.e., attached position, relative values of mass and spring constant, on the coupled vibration characteristics of the system has been clarified comparing with those of uncoupled ones. Optimal attached position to maximize coupled plate natural frequency is also investigated and shown in contour diagrams. The influence of an attached mass has also been investigated, as the limiting case whereby the spring stiffness of the ‘spring-mass’ system approaches infinity.     

The Rayleigh-Plesset equation in terms of volume with explicit shear losses

The most common nonlinear equation of motion for the damped pulsation of a spherical gas bubble in an infinite body of liquid is the Rayleigh-Plesset equation, expressed in terms of the dependency of the bubble radius on the conditions pertaining in the gas and liquid (the so-called ‘radius frame’). However over the past few decades several important analyses have been based on a heuristically derived small-amplitude expansion of the Rayleigh-Plesset equation which considers the bubble volume, instead of the radius, as the parameter of interest, and for which the dissipation term is not derived from first principles. So common is the use of this equation in some fields that the inherent differences between it and the ‘radius frame’ Rayleigh-Plesset equation are not emphasised, and it is important in comparing the results of the two equations to understand that they differ both in terms of damping, and in the extent to which they neglect higher order terms. This paper highlights these differences. Furthermore, it derives a ‘volume frame’ version of the Rayleigh-Plesset equation which contains exactly the same basic physics for dissipation, and retains terms to the same high order, as does the ‘radius frame’ Rayleigh-Plesset equation. Use of this equation will allow like-with-like comparisons between predictions in the two frames.     

Sound transmission loss of foam-filled honeycomb sandwich panels using statistical energy analysis and theoretical and measured dynamic properties

Comparisons between the experimental and predicted sound transmission loss values obtained from statistical energy analysis are presented for two foam-filled honeycomb sandwich panels. Statistical energy analysis (SEA) is a modeling procedure which uses energy flow relationships for the theoretical estimation of the sound transmission through structures in resonant motion. The accuracy of the prediction of the sound transmission loss using SEA greatly depends on accurate estimates of: (1) the modal density, (2) the internal loss factor, and (3) the coupling loss factor parameters of the structures. A theoretical expression for the modal density of sandwich panels is developed from a sixth-order governing equation. Measured modal density estimates of the two foam-filled honeycomb sandwich panels are obtained by using a three-channel spectral method with a spectral mass correction to allow for the mass loading of the impedance head. The effect of mass loading of the accelerometer is corrected in the estimations of both the total loss factor and radiation loss factor of the sandwich panels.     

Non-Abelian hydrodynamics and the flow of spin in spin-orbit coupled substances

Motivated by the heavy ion collision experiments there is much activity in studying the hydrodynamical properties of non-Abelian (quark-gluon) plasmas. A major question is how to deal with color currents. Although not widely appreciated, quite similar issues arise in condensed matter physics in the context of the transport of spins in the presence of spin-orbit coupling. The key insight is that the Pauli Hamiltonian governing the leading relativistic corrections in condensed matter systems can be rewritten in a language of SU(2) covariant derivatives where the role of the non-Abelian gauge fields is taken by the physical electromagnetic fields: the Pauli system can be viewed as Yang-Mills quantum-mechanics in a ‘fixed frame’, and it can be viewed as an ‘analogous system’ for non-Abelian transport in the same spirit as Volovik’s identification of the He superfluids as analogies for quantum fields in curved space time. We take a similar perspective as Jackiw and coworkers in their recent study of non-Abelian hydrodynamics, twisting the interpretation into the ‘fixed frame’ context, to find out what this means for spin transport in condensed matter systems. We present an extension of Jackiw’s scheme: non-Abelian hydrodynamical currents can be factored in a ‘non-coherent’ classical part, and a coherent part requiring macroscopic non-Abelian quantum entanglement. Hereby it becomes particularly manifest that non-Abelian fluid flow is a much richer affair than familiar hydrodynamics, and this permits us to classify the various spin transport phenomena in condensed matter physics in an unifying framework. The “particle based hydrodynamics” of Jackiw et al. is recognized as the high temperature spin transport associated with semiconductor spintronics. In this context the absence of faithful hydrodynamics is well known, but in our formulation it is directly associated with the fact that the covariant conservation of non-Abelian currents turns into a disastrous non-conservation of the incoherent spin currents of the high temperature limit. We analyze the quantum-mechanical single particle currents of relevance to mesoscopic transport with as highlight the Ahronov-Casher effect, where we demonstrate that the intricacies of the non-Abelian transport render this effect to be much more fragile than its abelian analog, the Ahronov-Bohm effect. We subsequently focus on spin flows protected by order parameters. At present there is much interest in multiferroics where non-collinear magnetic order triggers macroscopic electric polarization via the spin-orbit coupling. We identify this to be a peculiarity of coherent non-Abelian hydrodynamics: although there is no net particle transport, the spin entanglement is transported in these magnets and the coherent spin ‘super’ current in turn translates into electric fields with the bonus that due to the requirement of single valuedness of the magnetic order parameter a true hydrodynamics is restored. Finally, ‘fixed-frame’ coherent non-Abelian transport comes to its full glory in spin-orbit coupled ‘spin superfluids’, and we demonstrate a new effect: the trapping of electrical line charge being a fixed frame, non-Abelian analog of the familiar magnetic flux trapping by normal superconductors. The only known physical examples of such spin superfluids are the ³He A- and B-phase where unfortunately the spin-orbit coupling is so weak that it appears impossible to observe these effects.     

S-duality in 3D gravity with torsion

The deformation of the connection in three spacetime dimensions by the kinematically equivalent coframe is shown to induce a duality between the (Lorentz-) rotational and translational field momenta, for which the coupling to the deformation parameter is inverted. This new kind of strong/weak duality, pertinent to 3D, is instrumental for studying exact solutions of the 3D Poincaré gauge field equations and the particle content of propagating modes on a background of constant curvature. For a topological Chern-Simons model of gravity, the propagating modes ‘living’ on an Anti-de Sitter (AdS) background correspond to real massive particles. Yang-Mills type generalizations and new cubic Lagrangians are found and completely classified in 3D. AdS or black hole type solutions with constant axial torsion emerge, also for these higher-order Lagrangians with new ‘exotic’ torsion-curvature couplings. Their pattern complies with our S-duality, with new repercussions for the field redefinition of the metric in 3D quantum gravity and the cosmological constant problem.     

Spin scissors mode and the fine structure of M1 states in nuclei

The coupled dynamics of low lying modes, including the scissors mode, and various giant quadrupole resonances are studied with the help of the Wigner Function Moments method generalized to take into account spin degrees of freedom. Equations of motion for collective variables are derived on the basis of Time Dependent Hartree-Fock equations in the harmonic oscillator model including spin-orbit potential plus quadrupole-quadrupole residual interaction. Introducing spin allows one to consider new types of nuclear collective motion where the nucleons with spin ‘up’ oscillate against nucleons with spin ‘down’.     

Lumped parameter models representing impedance functions at the interface of a rod on a viscoelastic medium

In this study, a lumped parameter model that properly simulates the impedance characteristics at the extremity of a uniform, isotropic, homogeneous rod on a viscoelastic medium is proposed. The lumped parameter model consists of springs, dashpots, and so-called “gyro-mass elements”. The gyro-mass elements generate a reaction force proportional to the relative acceleration of the nodes between which they are placed. This model consists of units arranged in series, each unit consisting of a spring, a damper, and a gyro-mass element arranged in parallel. A formula is proposed for determining the properties of the elements in the units calculated from a closed-form solution based on a modal expansion. The impedance function simulated by the proposed model shows good agreement with the rigorous impedance function derived from the differential equation of motion of the rod. The results obtained by employing this model in some example applications show that the accuracy of the model is appreciably high when compared with conventional finite element models. A great advantage of this model is that a significant reduction of the number of degrees of freedom can be achieved for solving recent vibration problems with high-frequency excitations, such as ultrasonic vibrations.     

Quantum correlation cost of the weak measurement

Quantum correlation cost (QCC) characterizing how much quantum correlation is used in a weak-measurement process is presented based on the trace norm. It is shown that the QCC is related to the trace-norm-based quantum discord (TQD) by only a factor that is determined by the strength of the weak measurement, so it only catches partial quantumness of a quantum system compared with the TQD. We also find that the residual quantumness can be ‘extracted’ not only by the further von Neumann measurement, but also by a sequence of infinitesimal weak measurements. As an example, we demonstrate our outcomes by the Bell-diagonal state.