Search results

Recently, the original Wavenumber approach was introduced for dynamic analysis of dam-reservoir systems in frequency domain in the context of pure finite element programming. But its main disadvantages are that it cannot be implemented in time domain. The purpose of this paper is to propose an approximation to the original approach which enables one to carry out this effective method in time domain as well as in frequency domain. Based on the present investigation, it is proven that the Approximate Wavenumber approach has inherent characteristics, which allows it to be envisaged as an effective technique for calculating the response of concrete gravity dam-reservoir systems in time domain.

The method is described initially. Subsequently, the response of an idealized triangular dam-reservoir system is obtained by the proposed approach as well as by applying two other well-known absorbing conditions which are widely utilized in practice. The results are also controlled against the corresponding exact responses. It should be emphasized that all results presented herein are obtained by the FE-FE method under different absorbing conditions applied on the truncation boundary. These include two well-known absorbing conditions referred to as Sommerfeld and Sharan as well as the proposed approach of the present study (i.e. Approximate Wavenumber condition).

It is concluded that the maximum error for the Approximate Wavenumber approach is in the range of 10 percent at the major peaks of the response. This occurs mainly for the very low reservoir lengths under full reflective reservoir base condition and vertical excitation. This is a remarkable result for any kind of robust truncation boundary simulation that one may expect. The fundamental frequency of the system is captured correctly for the Approximate Wavenumber approach, even in cases of low reservoir length.

Based on this investigation, it is proven that the Approximate Wavenumber approach has inherent characteristics, which allows it to be envisaged as an effective technique for calculating the response of concrete gravity dam-reservoir systems in time domain. It is concluded that the maximum error for the Approximate Wavenumber approach is in the range of 10 percent at the major peaks of the response. This occurs mainly for the very low reservoir lengths under full reflective reservoir base condition and vertical excitation. This is a remarkable result for any kind of robust truncation boundary simulation that one may expect.

This paper aims to propose numerical-based and experiment-based identification processes, accounting for uncertainties to identify structural parameters, in a wave propagation framework.

A variant of the inhomogeneous wave correlation (IWC) method is proposed. It consists on identifying the propagation parameters, such as the wavenumber and the wave attenuation, from the frequency response functions. The latters can be computed numerically or experimentally. The identification process is thus called numerical-based or experiment-based, respectively. The proposed variant of the IWC method is then combined with the Latin hypercube sampling method for uncertainty propagation. Stochastic processes are consequently proposed allowing more realistic identification.

The proposed variant of the IWC method permits to identify accurately the propagation parameters of isotropic and composite beams, whatever the type of the identification process in which it is included: numerical-based or experiment-based. Its efficiency is proved with respect to an analytical model and the Mc Daniel method, considered as reference. The application of the stochastic identification processes shows good agreement between simulation and experiment-based results and that all identified parameters are affected by uncertainties, except damping.

The proposed variant of the IWC method is an accurate alternative for structural identification on wide frequency ranges. Numerical-based identification process can reduce experiments’ cost without significant loss of accuracy. Statistical investigations of the randomness of identified parameters illustrate the robustness of identification against uncertainties.

The purpose of this paper is to introduce a fast and simple method to calculate an estimation of parameters of interest of microstrip antennas, such as the resonance frequencies for example.

The Trefftz collocation method will be used to solve the governing differential equations of the problem. This method uses trial functions that satisfy, in a certain region the governing differential equations. Complete sets of solutions of such equations are required so that completeness and convergence can be guaranteed. The values of the wavenumbers for which the solution of the governing equation is unbounded, are those correspondent to the resonance frequencies of the antenna. After finding the wavenumbers, with the help of empirical correction formulas (because of the effect of the fringing field), the actual resonance frequencies are determined.

The Trefftz collocation method was found to be a very simple, fast and accurate method for the computation of the electric field under the patch of a microstrip antenna. Results obtained from this method showed excellent accuracy with less computational effort than other methods previously used.

Although the resonance wavenumbers may be accurate for any shape of antenna (because of the method convergence), the resonance frequencies might not be so accurate for irregular shapes since the parameters of the empirical formulas are approximated. Also the resonant cavity model is only valid for antennas made of thin substrates.

This formulation of the Trefftz method was for the first time applied to this problem, showing promising results.

This paper aims to present the development of a numerical continuum-discrete approach for computing the sensitivity of the waves propagating in periodic composite structures. The work can be directly used for evaluating the sensitivity of the structural dynamic performance with respect to geometric and layering structural modifications.

A structure of arbitrary layering and geometric complexity is modelled using solid finite element (FE). A generic expression for computing the variation of the mass and the stiffness matrices of the structure with respect to the material and geometric characteristics is hereby given. The sensitivity of the structural wave properties can thus be numerically determined by computing the variability of the corresponding eigenvalues for the resulting eigenproblem. The exhibited approach is validated against the finite difference method as well as analytical results.

An intense wavenumber dependence is observed for the sensitivity results of a sandwich structure. This exhibits the importance and potential of the presented tool with regard to the optimization of layered structures for specific applications. The model can also be used for computing the effect of the inclusion of smart layers such as auxetics and piezoelectrics.

The paper presents the first continuum-discrete approach specifically developed for accurately and efficiently computing the sensitivity of the wave propagation data for periodic composite structures irrespective of their size. The considered structure can be of arbitrary layering and material characteristics as FE modelling is used.

The purpose of this paper is to extend statistical energy analysis (SEA)‐like modeling to fluid‐structure coupled systems.

An equivalent approach of aerodynamic loads is applied to a SEA‐like modeling of a panel‐cavity coupled system with rain‐on‐the‐roof excitation. Two aerodynamic excitations are presented: turbulent boundary layer (TBL) and diffuse field excitation. The energetic description of the coupled system is studied with both aerodynamic excitations, taking in account the coincidence effects. In order to extent the approach to more general systems, some parameters of the coupled system are also modified and the accuracy of the coupled system modeling is investigated.

The boundary conditions of the panel and the coupling strength between the panel and the cavity have been modified. As it was expected, the accuracy of equivalent approach is shown to be independent of such modifications. The interest of such calculation is thus highlighted: modelings of systems and aerodynamic excitations are independent, and can be treated separately.

This result is interesting in the space industry, for launch vehicles are excited by different types of random excitations. Those excitations can be modeled by SEA‐like with low calculation time and memory and applied to a unique system modeling.

The purpose of this paper is to investigate the effect of initial stress and heterogeneity on the propagation of torsional waves in dissipative medium. The problem consists of dry sand poroelastic half-space embedded between heterogeneous self-reinforced half-space and poroelastic medium. The frequency equation is derived in the framework of Biot's theory with some variants.

Torsional wave propagation in dry sand poroelastic half-space embedded between self-reinforced half-space and poroelastic medium. All the constituents here are assumed to be dissipative, heterogeneous and initial stressed.

Phase velocity and attenuation are computed against wavenumber for various values of self-reinforcement parameter, inhomogeneity parameter and initial stress. Particular cases are discussed in absence of dissipation. The numerical results are presented graphically.

Initial stress and heterogeneity effects on torsional waves in dry sand half-space between reinforced half-space and poroelastic medium are investigated. The frequency equation is derived, and which intern gives the phase velocity and attenuation coefficient for various values of initial stress, self-reinforcement parameter and heterogeneity parameter. From the numerical results, it is clear that as wavenumber varies phase velocity and attenuation are periodic in nature for all the cases. Particular cases are discussed in absence of dissipation. This kind of analysis can be extended to any elastic solid by taking magnetic, thermo and piezoelectric effects into account.

Microstereolithography is capable of producing millimeter-scale polymer parts having micron-scale features. Material properties of the cured polymers can vary depending on build parameters such as exposure. Current techniques for determining the material properties of these polymers are limited to static measurements via micro/nanoindentation, leaving the dynamic response undetermined. The purpose of this paper is to demonstrate a method to measure the dynamic response of additively manufactured parts to infer the dynamic modulus of the material in the ultrasonic range.

Frequency-dependent material parameters, such as the complex Young’s modulus, have been determined for other relaxing materials by measuring the wave speed and attenuation of an ultrasonic pulse traveling through the materials. This work uses laser Doppler velocimetry to measure propagating ultrasonic waves in a solid cylindrical waveguide produced using microstereolithography to determine the frequency-dependent material parameters of the polymer. Because the ultrasonic wavelength is comparable with the part size, a model that accounts for both geometric and viscoelastic dispersive effects is used to determine the material properties using experimental data.

The dynamic modulus in the ultrasonic range of 0.4-1.3 MHz was determined for a microstereolithography part. Results were corroborated by using the same experimental method for an acrylic part with known properties and by evaluating the natural frequency and storage modulus of the same microstereolithography part with a shaker table experiment.

The paper demonstrates a method for determining the dynamic modulus of additively manufactured parts, including relatively small parts fabricated with microstereolithography.

The purpose of this study is to investigate the accuracy and applicability of the Flowfield Dependent Variation (FDV) method for large-eddy simulations (LES) of decaying isotropic turbulence.

In an earlier paper, the FDV method was successfully demonstrated for simulations of laminar flows with speeds varying from low subsonic to high supersonic Mach numbers. In the current study, the FDV method, implemented in a finite element framework, is used to perform LESs of decaying isotropic turbulence. The FDV method is fundamentally derived from the Lax–Wendroff Scheme (LWS) by replacing the explicit time derivatives in LWS with a weighted combination of explicit and implicit time derivatives. The increased implicitness and the inherent numerical dissipation of FDV contribute to the scheme’s numerical stability and monotonicity. Understanding the role of numerical dissipation that is inherent to the FDV method is essential for the maturation of FDV into a robust scheme for LES of turbulent flows. Accordingly, three types of LES of decaying isotropic turbulence were performed. The first two types of LES utilized explicit subgrid scale (SGS) models, namely, the constant-coefficient Smagorinsky and dynamic Smagorinsky models. In the third, no explicit SGS model was employed; instead, the numerical dissipation inherent to FDV was used to emulate the role played by explicit SGS models. Such an approach is commonly known as Implicit LES (ILES). A new formulation was also developed for quantifying the FDV numerical viscosity that principally arises from the convective terms of the filtered Navier–Stokes equations.

The temporal variation of the turbulent kinetic energy and enstrophy and the energy spectra are presented and analyzed. At all grid resolutions, the temporal profiles of kinetic energy showed good agreement with t^(−1.43) theoretical scaling in the fully developed turbulent flow regime, where t represents time. The energy spectra also showed reasonable agreement with the Kolmogorov’s k^(−5/3) power law in the inertial subrange, with the spectra moving closer to the Kolmogorov scaling at higher-grid resolutions. The intrinsic numerical viscosity and the dissipation rate of the FDV scheme are quantified, both in physical and spectral spaces, and compared with those of the two SGS LES runs. Furthermore, at a finite number of flow realizations, the numerical viscosities of FDV and of the Streamline Upwind/Petrov–Galerkin (SUPG) finite element method are compared. In the initial stages of turbulence development, all three LES cases have similar viscosities. But, once the turbulence is fully developed, implicit LES is less dissipative compared to the two SGS LES runs. It was also observed that the SUPG method is significantly more dissipative than the three LES approaches.

Just as any computational method, the limitations are based on the available computational resources.

Solving problems involving turbulent flows is by far the biggest challenge facing engineers and scientists in the twenty-first century, this is the road that the authors have embarked upon in this paper and the road ahead of is very long.

Understanding turbulence is a very lofty goal and a challenging one as well; however, if the authors succeed, the rewards are limitless.

The derivation of an explicit expression for the numerical viscosity tensor of FDV is an important contribution of this study, and is a crucial step forward in elucidating the fundamental properties of the FDV method. The comparison of viscosities for the three LES cases and the SUPG method has important implications for the application of ILES approach for turbulent flow simulations.

The purpose of this paper is to investigate the thermal convection inside a spatially modulated domain.

The governing equations are mapped onto an infinite strip, allowing Fourier expansion of the flow and temperature in the streamwise direction.

Similar to Rayleigh‐Benard convection, conduction is lost to convection at a critical Rayleigh number, which depends strongly on both the modulation amplitude and the wavenumber. The effect of modulation is found to be destabilizing (stabilizing) for conduction for relatively large (small) modulation wavelength. Oscillatory convection sets in as the Rayleigh number is increased.

This paper presents novel results.

The onset of Darcy–Benard convection in a fluid-saturated porous layer within channeled walls can be through either the convectively amplified perturbations or absolutely unstable modes. The convective type route to formation of cellular forms has received much more interest in the literature than the absolute mode. This paper aims to explore the absolute instability mechanism on triggering the Benard convection in the presence of a uniform magnetic field acting perpendicular to the channel walls.

Pursuing a completely theoretical linear stability approach, the locus of wavenumbers and critical Rayleigh numbers with respect to varying Peclet numbers leading to zero complex group velocity, and hence the absolute instability onset, is formulated in closed form. The formulae also incorporate the influence of fluid movement through the porous layer.

Wanenumbers are found to be increasing in magnitude under the effect of the magnetic field. Although the magnetic field expectedly behaves in favor of stabilizing the convection through both convective and absolute instabilities by pushing the threshold value of Rayleigh number to higher values, a peculiarity exists in such a manner that the stabilizing effect of magnetic field can no longer compete against the absolute instability mechanism within the magnetized field after some critical location possessing negative part of the wavenumber.

The significant outcome is attributed to the exchange of instabilities as far as the absolute instability mechanism is concerned. Magnetic field is always stabilizing and no such trend is observed regarding the convective type instability.