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.

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 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 the paper is to propose some modifications to the SIMPLE (semi-implicit method for pressure-linked equations) algorithm. These modifications can ensure the numerical robustness and optimize computational efficiency. They remarkably promote the ability of the SIMPLE algorithm for incompressible DNS (direct numerical simulation) of multiscale problems, such as transitional flows and turbulent flows, by improving the properties of dispersion and dissipation.

The MDCD (minimized dispersion and controllable dissipation) scheme and MMIM (modified momentum interpolation method) are introduced. Six typical test cases are used to validate the modified algorithm, including the linear convective flow, lid-driven cavity flow, laminar boundary layer, Taylor vortex and DHIT (decaying homogenous isotropic turbulence). Particularly, a highly unsteady DNS of separated-flow transition in turbomachinery is precisely predicted by the modified algorithm.

The numerical examples show the distinct superiority of the modified algorithm in both internal flows and external flows. The advantages of the MDCD scheme and MMIM make the SIMPLE algorithm a promising method for DNS.

Some effective modifications to the SIMPLE algorithm are addressed. It is the first attempt to introduce the MDCD approach into the SIMPLE-type algorithms. The new algorithm is especially suitable for the incompressible DNS of convection-dominated flows.

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 centrifugal instability mechanism of boundary layers over concave surfaces is responsible for the development of quasi-periodic, counter-rotating vortices aligned in a streamwise direction known as Görtler vortices. By distorting the boundary layer structure in both the spanwise and the wall-normal directions, Görtler vortices may modify heat transfer rates. The purpose of this study is to conduct spatial numerical simulation experiments based on a vorticity–velocity formulation of the incompressible Navier–Stokes system of equations to quantify the role of the transition in the heat transfer process.

Experiments are conducted using an in-house, parallel, message-passing code. Compact finite difference approximations and a spectral method are used to approximate spatial derivatives. A fourth-order Runge–Kutta method is adopted for time integration. The Poisson equation is solved using a geometric multigrid method.

Results show that the numerical method can capture the physics of transitional flows over concave geometries. They also show that the heat transfer rates in the late stages of the transition may be greater than those for either laminar or turbulent ones.

The numerical method can be considered as a robust alternative to investigate heat transfer properties in transitional boundary layer flows over concave surfaces.

The computational accuracy and performance of finite-difference time-domain (FDTD) methods are affected by the implementation of approximating derivative formulae in diverse ways. This study aims to focus on FDTD models featuring material dispersion with negligible losses and investigates two specific aspects that, until today, are usually examined in the context of non-dispersive media only. These aspects pertain to certain abnormal characteristics of coarsely resolved electromagnetic waves and the selection of the proper time-step size, in the case of a high-order discretization scheme.

Considering a Lorentz medium with negligible losses, the propagation characteristics of coarsely resolved waves is examined first, by investigating thoroughly the numerical dispersion relation of a typical discretization scheme. The second part of the study is related to the unbalanced space-time errors in FDTD schemes with dissimilar space-time approximation orders. The authors propose a remedy via the suitable choice of the time-step size, based on the single-frequency minimization of an error expression extracted, again, from the scheme’s numerical dispersion formula.

Unlike wave propagation in free space, there exist two parts of the frequency spectrum where waves in a Lorentz medium experience non-physical attenuation and display non-changing propagation constants, due to coarse discretization. The authors also show that an optimum time-step size can be determined, in the case of the (2,4) FDTD scheme, which minimizes the selected error formula at a specific frequency point, promoting more efficient implementations.

Unique characteristics displayed by discretized waves, which have been known for non-dispersive media, are examined and verified for the first time in the case of dispersive materials, thus completing the comprehension of the space-time discretization impact on simulated quantities. In addition, the closed-form formula of the optimum time-step enables the efficient implementation of the (2,4) FDTD method, minimizing the detrimental influence of the low-order temporal integration.

The purpose of this study is to develop an efficient model to solve the electromagnetic forward problem using a novel semi-analytical approach to compute the electromagnetic fields because of the presence of a scatterer.

The proposed model involves a novel formulation of a complete orthonormal set of radiating/nonradiating polarization currents. Furthermore, an integral equation-based representation is derived, and the appropriate boundary conditions are imposed to get the scattered electromagnetic field. An error term is introduced to evaluate the obtained results.

The proposed model was tested using several examples at different frequencies. The results of this study show that the novel representation exhibits fast convergence behavior and achieves highly accurate results, when compared to the results provided by the transmission line method.

The derived formulations presented in this study are significant in the electromagnetic forward modelling field because of the meaningful physical representation they provide. This is an important aspect that leads to precise calculation of electromagnetic fields for various applications.

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 present the results of experimental analysis of the elastic-guided wave mode conversion phenomenon in glass fiber-reinforced polymers. The results of this research presented in this paper are strictly focused on S0/A0’ mode conversion phenomenon caused by discontinuities in the form of circular Teflon inserts (artificial delaminations) and impact damage. Results of this research could be useful in problems of damage detection and localization.

In the research, guided waves are excited using a piezoelectric transducer and sensed in a non-contact manner using a scanning laser Doppler vibrometer. Full wavefield measurements are analyzed. Analysis of the influence of investigated discontinuities on S0/A0’ mode conversion is based on the elastic wave mode filtration in frequency-wavenumber domain. Mode filtration process allows us to remove the effects of the propagation of unwanted type of mode in forward or backward direction. Effects of S0/A0’ mode conversion are characterized by a mode conversion indicator (MCI) based on the amplitude of new mode A0’ and the amplitude of incident S0 mode.

It was noticed that the magnitude of MCI depends on the depth at which the Teflon inserts were located for all analyzed excitation frequencies and diameters of inserts (10 and 20 mm). The magnitude of MCI also increases with increasing impact energies. The S0/A0’ mode conversion phenomenon could be utilized for the detection of surface and internal located discontinuities.

This paper presents the original results of this research related to the influence of discontinuity location with respect to the sample thickness and severity of discontinuity on S0/A0’ mode conversion.