The frequency response function of the creep compliance

Motivated from the need to convert time-dependent rheometry data into complex frequency response functions, this paper studies the frequency response function of the creep compliance that is coined the complex creep function. While for any physically realizable viscoelastic model the Fourier transform of the creep compliance diverges in the classical sense, the paper shows that the complex creep function, in spite of exhibiting strong singularities, it can be constructed with the calculus of generalized functions. The mathematical expressions of the real and imaginary parts of the Fourier transform of the creep compliance of simple rheological networks derived in this paper are shown to be Hilbert pairs; therefore, returning back in the time domain a causal creep compliance. The paper proceeds by showing how a measured creep compliance of any solid-like or fluid-like viscoelastic material can be decomposed into elementary functions with parameters that can be identified from best fit of experimental data. The proposed technique allows for a direct determination of the sufficient parameters needed to approximate an experimentally measured creep compliance and the presented mathematical formulae offers dependable expressions of the corresponding complex-frequency response functions.

     

Non-linear wave propagation solutions by Fourier transform perturbation

A Fourier transform perturbation method is developed and used to obtain uniformly valid asymptotic approximations of the solution of a class of one-dimensional second order wave equations with small non-linearities. Multiple time scales are used and the initial-value problem on the infinite line is solved by Fourier transforming the wave equation and expanding the Fourier transform in powers of the small parameter. The non-linearity involves only the first partial derivatives of the dependent variable and the determination of the leading approximation is reduced to the solution of a pair of coupled non-linear ordinary differential equations in Fourier space. Examples are given involving a convolution non-linearity and a Van-der-Pol non-linearity.     

On the Fourier representation of elastic immersed boundaries

This paper presents an efficient treatment of fluid/elastic–structure interactions that takes advantage of the Fourier representation of immersed boundaries. We assume that the fluid is incompressible with uniform density and viscosity and that the immersed boundaries have fixed topologies. These elastic bodies can have large deformations and evolve anywhere within the fluid domain. They may be thick and are assumed to be piecewise smooth. We process the fluid–structure coupling with the immersed boundary method of C. S. Peskin. We can take advantage of the Fourier representation of the immersed bodies in many ways. First, the use of Fourier expansions allows us to filter out the high frequencies of the spatial oscillations along the boundary vectors. Second, we can work with a smaller number of boundary points to represent the interface, while preserving the same level of accuracy as long as enough points are used in the force spreading process. Finally, the harmonic information gathered by the Fourier coefficients is useful to control some global properties of the immersed boundaries. For example, we introduce a technique that corrects the volume conservation issue of closed immersed boundaries by performing constrained optimization in the Fourier space. We illustrate our method with two applications: one is a suspension flow with a large number of elastic ‘bubbles’, the other is an interesting case of artificial motion based on inertia rather than on flapping fins or flagella. Copyright © 2009 John Wiley & Sons, Ltd.     

Prediction of Gear Dynamics Using Fast Fourier Transform of Static Transmission Error*

ABSTRACT

An analytical computer simulation procedure for dynamic modeling of low-contact-ratio spur gear systems is presented. The procedure computes the gear static transmission error and uses a Fast Fourier Transform (FFT) to generate its frequency spectrum at various tooth profile modifications. The dynamic loading response of an unmodified (perfect involut) gear pair is compared with that of gears with profile modifications. Correlations are found between several profile modifications and the resulting dynamic loads. An effective error, obtained from frequency domain anal     

The Fractional Fourier Transform and Harmonic Oscillation

The ath-order fractional Fourier transform is a generalization ofthe ordinary Fourier transform such that the zeroth-order fractionalFourier transform operation is equal to the identity operation and thefirst-order fractional Fourier transform is equal to the ordinaryFourier transform. This paper discusses the relationship of thefractional Fourier transform to harmonic oscillation; both correspondto rotation in phase space. Various important properties of thetransform are discussed along with examples of commontransforms. Some of the applications of the transform are brieflyreviewed.     

Fast Fourier homogenization for elastic wave propagation in complex media

In the context of acoustic or elastic wave propagation, the non-periodic asymptotic homogenization method allows one to determine a smooth effective medium and equations associated with the wave propagation in a given complex elastic or acoustic medium down to a given minimum wavelength. By smoothing all discontinuities and fine scales of the original medium, the homogenization technique considerably reduces meshing difficulties as well as the numerical cost associated with the wave equation solver, while producing the same waveform as for the original medium (up to the desired accuracy). Nevertheless, finding the effective medium requires one to solve the so-called “cell problem”, which corresponds to an elasto-static equation with a finite set of distinct loadings. For general elastic or acoustic media, the cell problem is a large problem that has to be solved on the whole domain and its resolution implies the use of a finite element solver and a mesh of the fine scale medium. Even if solving the cell problem is simpler than solving the wave equation in the original medium (because it is time and source independent, based on simple tetrahedral meshes and embarrassingly parallel) it is still a challenge. In this work, we present an alternative method to the finite element approach for solving the cell problem. It is based on a well-known method designed by H. Moulinec and P. Suquet in 1998 in structural mechanics. This iterative technique relies on Green functions of a simple reference medium and extensively uses Fast Fourier Transforms. It is easy to implement, very efficient and relies on a simple regular gridding of the medium. Through examples we show that the method gives excellent results, even, under some conditions, for discontinuous media.     

From single-scale turbulence models to multiple-scale and subgrid-scale models by Fourier transform

A theoretical method based on mathematical physics formalism that allows transposition of turbulence modeling methods from URANS (unsteady Reynolds averaged Navier–Stokes) models, to multiple-scale models and large eddy simulations (LES) is presented. The method is based on the spectral Fourier transform of the dynamic equation of the two-point fluctuating velocity correlations with an extension to the case of non-homogenous turbulence. The resulting equation describes the evolution of the spectral velocity correlation tensor in wave vector space. Then, we show that the full wave number integration of the spectral equation allows one to recover usual one-point statistical closure whereas the partial integration based on spectrum splitting gives rise to partial integrated transport models (PITM). This latter approach, depending on the type of spectral partitioning used, can yield either a statistical multiple-scale model or subfilter transport models used in LES or hybrid methods, providing some appropriate approximations are made. Closure hypotheses underlying these models are then discussed by reference to physical considerations with emphasis on identification of tensorial fluxes that represent turbulent energy transfer or dissipation. Some experiments such as the homogeneous axisymmetric contraction, the decay of isotropic turbulence, the pulsed turbulent channel flow and a wall injection induced flow are then considered as typical possible applications for illustrating the potentials of these models.      

Sorting-free Hill-based stability analysis of periodic solutions through Koopman analysis

In this paper, we aim to study nonlinear time-periodic systems using the Koopman operator, which provides a way to approximate the dynamics of a nonlinear system by a linear time-invariant system of higher order. We propose for the considered system class a specific choice of Koopman basis functions combining the Taylor and Fourier bases. This basis allows to recover all equations necessary to perform the harmonic balance method as well as the Hill analysis directly from the linear lifted dynamics. The key idea of this paper is using this lifted dynamics to formulate a new method to obtain stability information from the Hill matrix. The error-prone and computationally intense task known by sorting, which means identifying the best subset of approximate Floquet exponents from all available candidates, is circumvented in the proposed method. The Mathieu equation and an n-DOF generalization are used to exemplify these findings.

     

Fourier Law,Phase Transitions and the Stationary Stefan Problem

We study the one-dimensional stationary solutions of the integro-differential equation which, as proved in Giacomin and Lebowitz (J Stat Phys 87:37–61, 1997; SIAM J Appl Math 58:1707–1729, 1998), describes the limit behavior of the Kawasaki dynamics in Ising systems with Kac potentials. We construct stationary solutions with non-zero current and prove the validity of the Fourier law in the thermodynamic limit showing that below the critical temperature the limit equilibrium profile has a discontinuity (which defines the position of the interface) and satisfies a stationary free boundary Stefan problem. Under-cooling and over-heating effects are also studied: we show that if metastable values are imposed at the boundaries then the mesoscopic stationary profile is no longer monotone and therefore the Fourier law is not satisfied. It regains its validity however in the thermodynamic limit where the limit profile is again monotone away from the interface.     

Hard-material Adhesion: Which Scales of Roughness Matter?

Background

Surface topography strongly modifies adhesion of hard-material contacts, yet roughness of real surfaces typically exists over many length scales, and it is not clear which of these scales has the strongest effect. Objective: This investigation aims to determine which scales of topography have the strongest effect on macroscopic adhesion.

Methods

Adhesion measurements were performed on technology-relevant diamond coatings of varying roughness using spherical ruby probes that are large enough (0.5-mm-diameter) to sample all length scales of topography. For each material, more than 2000 measurements of pull-off force were performed in order to investigate the magnitude and statistical distribution of adhesion. Using sphere-contact models, the roughness-dependent effective values of work of adhesion were measured, ranging from 0.08 to 7.15 mJ/m² across the four surfaces. The data was more accurately fit using numerical analysis, where an interaction potential was integrated over the AFM-measured topography of all contacting surfaces.

Results

These calculations revealed that consideration of nanometer-scale plasticity in the materials was crucial for a good quantitative fit of the measurements, and the presence of such plasticity was confirmed with AFM measurements of the probe after testing. This analysis enabled the extraction of geometry-independent material parameters; the intrinsic work of adhesion between ruby and diamond was determined to be 46.3 mJ/m². The range of adhesion was 5.6 nm, which is longer than is typically assumed for atomic interactions, but is in agreement with other recent investigations. Finally, the numerical analysis was repeated for the same surfaces but this time with different length-scales of roughness included or filtered out.

Conclusions

The results demonstrate a critical band of length-scales—between 43 nm and 1.8 µm in lateral size—that has the strongest effect on the total adhesive force for these hard, rough contacts.

     

Summation of trigonometric series by Fourier transforms

This paper gives the theorems concerning the summation of trigonometric series with the help of Fourier transforms. By means of the known results of Fourier transforms, many difficult and complex problems of summation of trigonometric series can be solved. This method is a comparatively unusual way to find the summation of trigonometric series, and has been used to establish the comprehensive table of summation of trigonometric series. In this table 10 thousand series are given and most of them are new.     

大体积混凝土拆模后表层温度的Fourier震荡

将大体积混凝土拆模后很快予以保护的过程简化为大气温度的线性脉冲。研究了脉冲时间对Fourier震荡的影响。成果可以解释混凝土温度在边界层的不合理现象,指出了修正Fourier震荡的途径。     

Application of the Fourier transform in electronic speckle photography

In speckle photograph technology, to determine the displacement of the points on the surface of the measured body, the conventional method is to put the film which has recorded the speckle patterns before and after displacement into a system of optical Fourier transforms. After filtering on the spectrum plane, the experimentalist can obtain the displacement information from the interference pattern on the image plane. Instead of setting up a complex optical Fourier transform system, we consider the speckle field as a light intensity function of 2 dimensions, which will change with different positioning of the points. After working on the function's discrete Fourier transform (DFT), according to one of the properties of Fourier transformation, the displacement of the measured point is involved in the phase of its spectrum. Having extracted the displacement information from the phase, we obtain the distribution of the displacement field. In this paper, we deduce the expression for the displacement field by using Fourier transformation under conditions of both equal and unequal displacement and show their applications.     

A wavelet transform analysis applied to unsteady aspects of supersonic jet screech resonance

 The multiple acoustic modes and shear layer instability waves which characterize a supersonic underexpanded rectangular jet are investigated experimentally via the Morlet wavelet transform. Because of its quasi-locality in both physical-space and Fourier space, the wavelet transformation allows one to track the time evolution of the various scales in both acoustic and velocity fluctuation signals. Using this technique it is demonstrated that multiple acoustic modes produced by the jet coexist and are not due to a mode switching mechanism. Unsteady screech tone modulation is observed and a mechanism for its occurrence is proposed. Received: 9 February 1996 / Accepted: 17 June 1996     

Numerical Investigation of Turbulent Kinetic Energy Dynamics in Chemically-Reacting Homogeneous Turbulence

In this work, we investigate numerically the temporal evolution of Turbulent Kinetic Energy (TKE) of a chemically-reacting n-heptane and air mixture in statistically Homogeneous Isotropic Turbulence (HIT). Our specific focus is on the concurrent view of TKE evolution in both physical and scale (Fourier) spaces to identify the impact of reaction-induced heat release on turbulence. The simulation parameters are selected to represent the combustion characteristics of heavy hydrocarbon fuels under engine conditions. Results indicate that pressure dilatation work dominates the TKE evolution during the period of strong heat release and its dominance is attributed to the strong volumetric dilatation associated with the presence of reaction fronts in physical space. Viscous dissipation and viscous dilatation terms become much stronger with increasing heat release, primarily due to the increase in strain-rate and dilatation at the vicinity of the reaction fronts, but their magnitudes are still small compared to that of pressure dilatation work. In addition, the analysis in Fourier space shows that pressure dilatation work dominates the evolution of TKE not only in the mean, but also over a wide range of scales. The spectrum of pressure dilatation shows a power-law behavior, which is a direct consequence of the localized sheet-like reaction fronts in physical space. It is also shown that viscous dissipation spectrum initially removes kinetic energy at small scales when heat release is weak, but starts to remove kinetic energy at intermediate and later at large scales due to the presence of localized reaction fronts during the strong heat release period. More interestingly, it is observed that the inter-scale kinetic energy transfer spectrum moves energy from less dissipative scales (small scales) to scales where kinetic energy is more effectively removed by viscous dissipation work (large scales) during the period of strong heat release, which indicates possible up-scale kinetic energy transfer in Fourier space.     

Fractal Pattern for Multiscale Digital Image Correlation

Background:

Digital Image Correlation (DIC) is based on the matching, between reference and deformed state images, of features contained in patterns that are deposited on test sample surfaces. These features are often suitable for a single scale, and there is a current lack of multiscale patterns capable of providing reliable displacement measurements over a wide range of scales.

Objective:

Here, we aim to demonstrate that a pattern based on a fractal (self-affine) surface would make a suitable pattern for multiscale DIC.

Methods:

A method to numerically generate patterns directly from a desired auto-correlation function is introduced. It is then enhanced by a Mean Intensity Gradient (MIG) improvement process based on grey level redistribution. Numerical experiments at multiple scales are performed for two different imposed displacement fields and results for one of the patterns generated are compared with those obtained for a random pattern and a Perlin noise one.

Results:

The proposed pattern is shown to lead to DIC errors comparable to those found with the two others for the first scales, but has much greater robustness. More importantly, the pattern generated here exhibits stable errors and robustness with respect to the scale whereas these two outputs become significantly degraded for the other two patterns as the scale increases.

Conclusions:

As a result, scale invariance properties of the pattern based on fractal surfaces correspond to scale invariance in DIC errors as well. This is of great interest regarding the use of such patterns in multiscale DIC.

     

The Fourier pseudo-spectral method for the SRLW equation

In this paper we present a Fourier pseudo-spectral method with a restraint operator for the SRLW equation. We prove the stability of the schemes and give optimum error estimates.     

Invariant Analysis of Turbulent Pipe Flow

The anisotropy analysis of Lumley provides a useful tool to quantify the degree of anisotropy in turbulent flows. Also included in the analysis are relations which may be used to check if the flow is axisymmetric or two-dimensional. However, the method does not provide any scale information about the structures. The analysis has therefore been extended here to Fourier space, which allows scale information to be derived. The method was applied to fully developed pipe flow and it was shown that the large-scale motion is everywhere close to axisymmetric. The intermediate scales are strongly influenced by the restrictions posed by the pipe walls. At the centre line, the flow structure appears axisymmetric at all scales, but the measurement sindicate that true axisymmetry is lost very quickly away from the centre line. The structure of the smallest scales could not be determined reliably due to a singularity in the analysis which develops as the scales go to zero. This revised version was published online in July 2006 with corrections to the Cover Date.     

Global measurement of water waves by Fourier transform profilometry

In this paper, we present an optical profilometric technique that allows for single-shot global measurement of free-surface deformations. This system consists of a high-resolution system composed of a videoprojector and a digital camera. A fringe pattern of known characteristics is projected onto the free surface and its image is registered by the camera. The deformed fringe pattern arising from the surface deformations is later compared to the undeformed (reference) one, leading to a phase map from which the free surface can be reconstructed. Particularly, we are able to project wavelength-controlled sinusoidal fringe patterns, which considerably increase the overall performance of the technique and the quality of the reconstruction compared to that obtained with a Ronchi grating. In comparison to other profilometric techniques, it allows for single-shot non-intrusive measurement of surface deformations over large areas. In particular, our measurement system and analysis technique is able to measure free surface deformations with sharp slopes up to 10 with a 0.2 mm vertical resolution over an interrogation window of size 450 × 300 mm² sampled on approximately 6.1 × 10⁶ measurement points. Some illustrative examples of the application of this measuring system to fluid dynamics problems are presented.

---