Compact finite difference - Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

A Novel High Order Space-Time Spectral Method for the Time Fractional Fokker--Planck Equation

The fractional Fokker-Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the non-local property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker-Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker-Planck equation attains the same approximation order as the time fractional diffusion equation developed in [23] by using the present method. That indicates an exponential decay may be achieved if the exact solution is sufficiently smooth. Finally, some numerical results are given to demonstrate the high order accuracy and efficiency of the new numerical scheme. The results show that the errors of the numerical solutions obtained by the space-time spectral method decay exponentially.

Analysis of β-amyloid (Aβ) deposition in the temporal lobe in Alzheimer's disease using Fourier (spectral) analysis

To determine the spatial pattern of ß-amyloid (Aß) deposition throughout the temporal lobe in Alzheimer's disease (AD). Methods: Sections of the complete temporal lobe from six cases of sporadic AD were immunolabelled with antibody against Aß. Fourier (spectral) analysis was used to identify sinusoidal patterns in the fluctuation of Aß deposition in a direction parallel to the pia mater or alveus. Results: Significant sinusoidal fluctuations in density were evident in 81/99 (82%) analyses. In 64% of analyses, two frequency components were present with density peaks of Aß deposits repeating every 500–1000 µm and at distances greater than 1000 µm. In 25% of analyses, three or more frequency components were present. The estimated period or wavelength (number of sample units to complete one full cycle) of the first and second frequency components did not vary significantly between gyri of the temporal lobe, but there was evidence that the fluctuations of the classic deposits had longer periods than the diffuse and primitive deposits. Conclusions: (i) Aß deposits exhibit complex sinusoidal fluctuations in density in the temporal lobe in AD; (ii) fluctuations in Aß deposition may reflect the formation of Aß deposits in relation to the modular and vascular structure of the cortex; and (iii) Fourier analysis may be a useful statistical method for studying the patterns of Aß deposition both in AD and in transgenic models of disease.

A Crank--Nicolson ADI spectral method for a two-dimensional riesz space fractional nonlinear reaction-diffusion equation

In this paper, a new alternating direction implicit Galerkin--Legendre spectral method for the two-dimensional Riesz space fractional nonlinear reaction-diffusion equation is developed. The temporal component is discretized by the Crank--Nicolson method. The detailed implementation of the method is presented. The stability and convergence analysis is strictly proven, which shows that the derived method is stable and convergent of order $2$ in time. An optimal error estimate in space is also obtained by introducing a new orthogonal projector. The present method is extended to solve the fractional FitzHugh--Nagumo model. Numerical results are provided to verify the theoretical analysis.

Fourier spectral methods for fractional-in-space reaction-diffusion equations

Fractional differential equations are becoming increasingly used as a powerful modelling approach for understanding the many aspects of nonlocality and spatial heterogeneity. However, the numerical approximation of these models is demanding and imposes a number of computational constraints. In this paper, we introduce Fourier spectral methods as an attractive and easy-to-code alternative for the integration of fractional-in-space reaction-diffusion equations described by the fractional Laplacian in bounded rectangular domains ofRn. The main advantages of the proposed schemes is that they yield a fully diagonal representation of the fractional operator, with increased accuracy and efficiency when compared to low-order counterparts, and a completely straightforward extension to two and three spatial dimensions. Our approach is illustrated by solving several problems of practical interest, including the fractional Allen–Cahn, FitzHugh–Nagumo and Gray–Scott models, together with an analysis of the properties of these systems in terms of the fractional power of the underlying Laplacian operator.

A novel variable resolution global spectral method on the sphere

A variable resolution global spectral method is created on the sphere using High resolution Tropical Belt Transformation (HTBT). HTBT belongs to a class of map called reparametrisation maps. HTBT parametrisation of the sphere generates a clustering of points in the entire tropical belt; the density of the grid point distribution decreases smoothly in the domain outside the tropics. This variable resolution method creates finer resolution in the tropics and coarser resolution at the poles. The use of FFT procedure and Gaussian quadrature for the spectral computations retains the numerical efficiency available with the standard global spectral method. Accuracy of the method for meteorological computations are demonstrated by solving Helmholtz equation and non-divergent barotropic vorticity equation on the sphere. (C) 2011 Elsevier Inc. All rights reserved.

Fourier Transform Method for Measuring the Micro Amplitu of Water Surface Waves

采用一种非接触的光学方法傅立叶变换莫尔法(Fourier transform method),结合数字图像处理技术,对微幅振荡的水表面波的振幅进行测量.它是对全场中每一个像素点进行测量,比接触测量法具有更高的灵敏度.它为微幅水表面波振幅的测量提供了一种手段.通过将计算机生成的周期性光栅图像经投影机直接投影到被测物体的参考平面,经CCD摄像头、图像板捕捉存储形成数字化的光栅图像,利用傅立叶变换莫尔法处理光栅图像,从而获得包含有水表面波的振幅的相位信息,再经适当的几何变换获得振幅信息.我们在垂直振荡装置上进行了不同激励频率和不同振幅的表面波的振幅测量.

Measurement of gain spectrum for Fabry-Perot semiconductor lasers by the Fourier transform method with a deconvolution process

To improve the accuracy of measured gain spectra, which is usually limited by the resolution of the optical spectrum analyzer (OSA), a deconvolution process based on the measured spectrum of a narrow linewidth semiconductor laser is applied in the Fourier transform method. The numerical simulation shows that practical gain spectra can be resumed by the Fourier transform method with the deconvolution process. Taking the OSA resolution to be 0.06, 0.1, and 0.2 nm, the gain-reflectivity product spectra with the difference of about 2% are obtained for a 1550-nm semiconductor laser with the cavity length of 720 pm. The spectra obtained by the Fourier transform method without the deconvolution process and the Hakki-Paoli method are presented and compared. The simulation also shows that the Fourier transform method has less sensitivity to noise than the Hakki-Paoli method.

Study on tapered crossed subwavelength gratings by Fourier modal method

Fourier modal method incorporating staircase approximation is used to study tapered crossed subwavelength gratings in this paper. Three intuitive formulations of eigenvalue functions originating from the prototype are presented, and their convergences are compared through numerical calculation. One of them is found to be suitable in modeling the diffraction efficiency of the circular tapered crossed subwavelength gratings without high absorption, and staircase approximation is further proven valid for non-highly-absorption tapered gratings. This approach is used to simulate the "moth-eye" antireflection surface on silicon, and the numerical result agrees well with the experimental one.

Fast and accurate SCFT calculations for periodic block-copolymer morphologies using the spectral method with Anderson mixing

We study the numerical efficiency of solving the self-consistent field theory (SCFT) for periodic block-copolymer morphologies by combining the spectral method with Anderson mixing. Using AB diblock-copolymer melts as an example, we demonstrate that this approach can be orders of magnitude faster than competing methods, permitting precise calculations with relatively little computational cost. Moreover, our results raise significant doubts that the gyroid (G) phase extends to infinite $\chi N$. With the increased precision, we are also able to resolve subtle free-energy differences, allowing us to investigate the layer stacking in the perforated-lamellar (PL) phase and the lattice arrangement of the close-packed spherical (S$_{cp}$) phase. Furthermore, our study sheds light on the existence of the newly discovered Fddd (O$^{70}$) morphology, showing that conformational asymmetry has a significant effect on its stability.

Linear readout of dynamic phase modulation index in an optical voltage sensor using the new J(1) ... J(3) spectral method

The J(1)...J(3) is a recent optical method for linear readout of dynamic phase modulation index in homodyne interferometers. In this work, the J(1)... J(3) method is applied to measure voltage in an optical voltage sensor. Based on the classical J(1)...J(4) method, the J(1)... J(3) technique shows to be more stable to phase drift and simpler for implementation than the original one. The sensor dynamic range is enhanced. The agreement between theoretical and experimental results, based on 1/f noise, is demonstrated.

Time and Frequency Domain Finite Element Models for Axial Wave Analysis in Hyperelastic Rods

The propagation of axial waves in hyperelastic rods is studied using both time and frequency domain finite element models. The nonlinearity is introduced using the Murnaghan strain energy function and the equations governing the dynamics of the rod are derived assuming linear kinematics. In the time domain, the standard Galerkin finite element method, spectral element method, and Taylor-Galerkin finite element method are considered. A frequency domain formulation based on the Fourier spectral method is also developed. It is found that the time domain spectral element method provides the most efficient numerical tool for the problem considered.

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.