Search results

The spectral boundary element method for solving two‐dimensional transient heat conduction problems is developed. This method is combined with the fast Fourier transform (FFT) to convert the solution between the time and frequency domains. The fundamental solutions in the frequency domain, required for the present method, are discussed. The resulting line integrations in the frequency domain are discretized using constant boundary elements and used in a Fortran boundary element program. Three examples are used to illustrate the accuracy and effectiveness of the method in both the frequency and time domains. First, the frequency domain solution procedure is verified using the steady‐state example of a semi‐infinite half space with a heat flux applied to a patch of the surface. This spectral boundary element method is then applied to the problem of a circular hole in an infinite solid subjected to a time‐varying heat flux, and solutions in both the frequency and time domains are presented. Finally, the method is used to solve the circular hole problem with a convection boundary condition. The accurary of these results leads to the conclusion that the spectral boundary element method in conjunction with the FFT is a viable option for transient problems. In addition, this spectral approach naturally produces frequence domain information which is itself of interest.

This paper aims to present a framework for deriving analytical and semi‐numerical models for coupled conductive‐convective heat transfer processes in a borehole heat exchanger subjected to general initial and boundary conditions.

The discrete Fourier transform and the spectral element method have been utilized for deriving two spectral models for a single U‐tube borehole heat exchanger in contact with a soil mass.

Verification and numerical examples have shown that the developed models are accurate and computationally very efficient. It is illustrated that one spectral element is capable of producing results which are more accurate than those produced by 200 finite elements.

The gained computational efficiency and accuracy will boost considerably the possibilities for more insight into geothermal analysis, which will improve the procedure for designing competitive energy extraction systems.

The models are capable of calculating exactly the temperature distribution in all involved borehole heat exchanger components and their thermal interactions.

Deals with the formulation and application of temporal and spatial spectral element approximations for the solution of convection‐diffusion problems. Proposes a new high‐order splitting space‐time spectral element method which exploits space‐time discontinuous Galerkin for the first hyperbolic substep and space continuous‐time discontinuous Galerkin for the second parabolic substep. Analyses this method and presents its characteristics in terms of accuracy and stability. Also investigates a subcycling technique, in which several hyperbolic substeps are taken for each parabolic substep; a technique which enables fast, cost‐effective time integration with little loss of accuracy. Demonstrates, by a numerical comparison with other coupled and splitting space‐time spectral element methods, that the proposed method exhibits significant improvements in accuracy, stability and computational efficiency, which suggests that this method is a potential alternative to existing schemes. Describes several areas for future research.

The purpose of the study is to present a frequency domain spectral finite element model (SFEM) based on fast Fourier transform (FFT) for wave propagation analysis of smart laminated composite beams with embedded delamination. For generating and sensing high-frequency elastic waves in composite beams, piezoelectric materials such as lead zirconate titanate (PZT) are used because they can act as both actuators and sensors. The present model is used to investigate the effects of parametric variation of delamination configuration on the propagation of fundamental anti-symmetric wave mode in piezoelectric composite beams.

The spectral element is derived from the exact solution of the governing equation of motion in frequency domain, obtained through fast Fourier transformation of the time domain equation. The beam is divided into two sublaminates (delamination region) and two base laminates (integral regions). The delamination region is modeled by assuming constant and continuous cross-sectional rotation at the interfaces between the base laminate and sublaminates. The governing differential equation of motion for delaminated composite beam with piezoelectric lamina is obtained using Hamilton’s principle by introducing an electrical potential function.

A detailed study of the wave response at the sensor shows that the A0 mode can be used for delamination detection in a wide region and is more suitable for detecting small delamination. It is observed that the amplitude and time of arrival of the reflected A0 wave from a delamination are strongly dependent on the size, position of the delamination and the stacking sequence. The degraded material properties because of the loss of stiffness and density in damaged area differently alter the S0 and A0 wave response and the group speed. The present method provides a potential technique for researchers to accurately model delaminations in piezoelectric composite beam structures. The delamination position can be identified if the time of flight of a reflected wave from delamination and the wave propagation speed of A0 (or S0) mode is known.

Spectral finite element modeling of delaminated composite beams with piezoelectric layers has not been reported in the literature yet. The spectral element developed is validated by comparing the present results with those available in the literature. The spectral element developed is then used to investigate the wave propagation characteristics and interaction with delamination in the piezoelectric composite beam.

The purpose of this paper is to develop a new formulation using spectral approach, which can predict the wave behavior to uncertain parameters in mid and high frequencies.

The work presented is based on a hybridization of a spectral method called the “wave finite element (WFE)” method and a non-intrusive probabilistic approach called the “polynomial chaos expansion (PCE).” The WFE formulation for coupled structures is detailed in this paper. The direct connection with the conventional finite element method allows to identify the diffusion relation for a straight waveguide containing a mechanical or geometric discontinuity. Knowing that the uncertainties play a fundamental role in mid and high frequencies, the PCE is applied to identify uncertainty propagation in periodic structures with periodic uncertain parameters. The approach proposed allows the evaluation of the dispersion of kinematic and energetic parameters.

The authors have found that even though this approach was originally designed to deal with uncertainty propagation in structures it can be competitive with its low time consumption. The Latin Hypercube Sampling (LHS) is also employed to minimize CPU time.

The approach proposed is quite new and very simple to apply to any periodic structures containing variabilities in its mechanical parameters. The Stochastic Wave Finite Element can predict the dynamic behavior from wave sensitivity of any uncertain media. The approach presented is validated for two different cases: coupled waveguides with and without section modes. The presented results are verified vs Monte Carlo simulations.

A splitting technique for solutions of the Navier—Stokes and the energy equations, in Boussinesq approximately, is presented. The equations are first integrated in time using a splitting procedure and then discretized spatially by means of a high‐order spectral element method. The whole technique is validated on the flow in a differentially‐heated cavity at intermediate and transitional Rayliegh numbers. The results are in a very good agreement with other available numerical solutions.

The purpose of this paper is to develop a numerical framework based on the accurate spectral element method (SEM) to simulate the mixed convective heat transfer within a porous enclosure with three adiabatic thin baffles of different lengths in nine cases and analyze the effects of several parameters.

The authors develop an improved time-splitting method to solve the Darcy–Brinkman–Forchheimer model. No extra assumptions are introduced for the intermediate velocity, and the final velocity field satisfies the incompressible constraint strictly compared with the classical method. The governing equations are split into a pure convection problem, a Stokes problem and a thermal diffusion problem. The least-squares variation is adopted for the Stokes problem, and the Galerkin variation is used for the other two problems, such that the pressure and velocity can be discretized with the same interpolation order, which benefits the numerical accuracy and program design.

Regarding the method, the excellent spectral accuracy, the capability of discretizing complex computational regions and the improved time-splitting methods make SEM an effective tool to accurately predict the non-Darcy convective heat transfer; as for the numerical tests, it is observed that weakened convection and heat transfer are induced by the increasing length of the baffles. The flow and heat transfer in channel 1 is only related to the length of baffle 1 because of the downward-driven right sidewall, and it is more difficult for baffle 3 to form the secondary flow on its tip.

A novel numerical framework for Darcy–Brinkman–Forchheimer model is developed, expanding the application of SEM for simulating non-Darcy convective heat transfer to improve the numerical accuracy. Numerical results and analysis for flow and heat fields could help designers understand the control of heat transfer using adiabatic baffles better.

This paper aims to present a high-order algorithm implemented with the modal spectral element method and simulations of three-dimensional thermal convective flows by using the full viscous dissipation function in the energy equation. Three benchmark problems were solved to validate the algorithm with exact or theoretical solutions. The heated rotating sphere at different temperatures inside a cold planar Poiseuille flow was simulated parametrically at varied angular velocities with positive and negative rotations.

The fourth-order stiffly stable schemes were implemented and tested for time integration. To provide the hp-refinement and spatial resolution enhancement, a modal spectral element method using hierarchical basis functions was used to solve governing equations in a three-dimensional space.

It was found that the direction of rotation of the heated sphere has totally different effects on drag, lateral force and torque evaluated on surfaces of the sphere and walls. It was further concluded that the angular velocity of the heated sphere has more influence on the wall normal velocity gradient than on the wall normal temperature gradients and therefore, more influence on the viscous dissipation than on the thermal dissipation.

This paper concerns incompressible fluid flow at constant properties with up to medium temperature variations in the absence of thermal radiation and ignoring the pressure work.

This paper contributes a viable high-order algorithm in time and space for modeling convective heat transfer involving an internal heated rotating sphere with the effect of viscous heating.

Results of this paper could provide reference for related topics such as enhanced heat transfer forced convection involving rotating spheres and viscous thermal effect.

The merits include resolving viscous dissipation and thermal diffusion in stationary and rotating boundary layers with both h- and p-type refinements, visualizing the viscous heating effect with the full viscous dissipation function in the energy equation and modeling the forced advection around a rotating sphere with varied positive and negative angular velocities subject to a shear flow.

The research is directed at development of an efficient and accurate technique for modelling incompressible free surface flows in which viscous effects may not be neglected. The paper describes the methodology, and gives illustrative results for simple geometries.

The pseudospectral matrix element method of discretisation is selected as the basis for the CFD technique adopted, because of its high spectral accuracy. It is implemented as a means of solving the Navier‐Stokes equations coupled with the modified compressibility method.

The viscous solver has been validated for the benchmark cases of uniform flow past a cylinder at a Reynolds number of 40, and 2D cavity flows. Results for sloshing of a viscous fluid in a tank have been successfully compared with those from a linearised analytical solution. Application of the method is illustrated by the results for the interaction of an impulsive wave with a surface piercing circular cylinder in a cylindrical tank.

The paper demonstrates the viability of the approach adopted. The limitation of small amplitude waves should be tackled in future work.

The results will have particular significance in the context of validating computations from more complex schemes applicable to arbitrary geometries.

The new methodology and results are of interest to the community of those developing numerical models of flow past marine structures.

The purpose of this paper is to develop a spectral element modeling to predict electromechanical admittance in the surface-bonded piezoelectric wafer and beam structure considering temperature effects.

For modeling the beam, the axial and transverse vibrations of the beam have been considered, and temperature-dependent mechanical and electromechanical properties of piezoelectric wafer active sensor and aluminum have been analyzed. The influences of temperature effects on electromechanical admittance are investigated.

The results show that a frequency left shift and a decrease in amplitude of admittance in any natural frequencies with increasing temperature have been observed. The mechanism of these changes is discussed.

The numerical results may be considered helpful for structural health monitoring using electromechanical impedance technique.