scholarly journals A Dynamical Core with Double Fourier Series: Comparison with the Spherical Harmonics Method

2006 ◽  
Vol 134 (4) ◽  
pp. 1299-1315 ◽  
Author(s):  
Hyeong-Bin Cheong

Abstract A dynamical core of a general circulation model with the spectral method using double Fourier series (DFS) as basis functions is presented. The model uses the hydrostatic balance approximation and sigma coordinate system in the vertical direction and includes no topography. The model atmosphere is divided into 25 layers with equal sigma depths. Prognostic equations for the vorticity, divergence, temperature, and logarithmic surface pressure are solved by the DFS spectral-transform method with the Fourier filtering at middle and high latitudes. A semi-implicit time-stepping procedure, which deals with the eigendecomposition and inversion of the 3D Helmholtz equation associated with the gravity wave terms, is incorporated for the gravity wave–related terms. The DFS model is tested in terms of the solution of the 3D Helmholtz equation, balanced initial state, developing baroclinic waves, and short- and long-term Held–Suarez–Williamson simulations for T42, T62, T84, and T106 resolutions. It is found that the DFS model is stable and accurate and produces almost the same results as the spherical harmonics method (SHM). The normalized difference (i.e., L2 norm error) measured from the results of highest-resolution SHM-T106 showed a desirable convergence of the DFS solution with the resolution. The convergence property, however, varies with the test case and prognostic variables. The total mass (or global integrated surface pressure) is conserved to a good approximation in the long-term simulations. Computing on the high-performance computer NEC SX-5 (parallel-vector architecture) indicated that DFS is more efficient than the SHM and the efficiency increases with the resolution, for example, by factors of 2.09 and 7.68 for T212 and T1022, respectively.

2013 ◽  
Vol 141 (9) ◽  
pp. 3052-3061 ◽  
Author(s):  
Hoon Park ◽  
Song-You Hong ◽  
Hyeong-Bin Cheong ◽  
Myung-Seo Koo

Abstract This study describes an application of the double Fourier series (DFS) spectral method developed by Cheong as an alternative dynamical option in a model system that was ported into the Global/Regional Integrated Model System (GRIMs). A message passing interface (MPI) for a massive parallel-processor cluster computer devised for the DFS dynamical core is also presented. The new dynamical core with full physics was evaluated against a conventional spherical harmonics (SPH) dynamical core in terms of short-range forecast capability for a heavy rainfall event and seasonal simulation framework. Comparison of the two dynamical cores demonstrates that the new DFS dynamical core exhibits performance comparable to the SPH in terms of simulated climatology accuracy and the forecast of a heavy rainfall event. Most importantly, the DFS algorithm guarantees improved computational efficiency in the cluster computer as the model resolution increases, which is consistent with theoretical values computed from a dry primitive equation model framework. The current study shows that, at higher resolutions, the DFS approach can be a competitive dynamical core because the DFS algorithm provides the advantages of both the spectral method for high numerical accuracy and the gridpoint method for high performance computing in the aspect of computational cost.


2011 ◽  
Vol 60 (2) ◽  
pp. 137-148
Author(s):  
Igor Korotyeyev ◽  
Beata Zięba

Steady-state modelling method for matrix-reactance frequency converter with boost topologyThis paper presents a method intended for calculation of steady-state processes in AC/AC three-phase converters that are described by nonstationary periodical differential equations. The method is based on the extension of nonstationary differential equations and the use of Galerkin's method. The results of calculations are presented in the form of a double Fourier series. As an example, a three-phase matrix-reactance frequency converter (MRFC) with boost topology is considered and the results of computation are compared with a numerical method.


2018 ◽  
pp. 1
Author(s):  
Amira Ali Een Faied

A function of two variables may be expanded in a double Fourier series, as a function of one variable is expanded in an ordinary Fourier series. Purpose that the function f ( x, y ) possesses a double Lebesgue integral over the square (– π < π ; – π < y < π ). Then the general term of the double Fourier series of this function is given by cos = є mn { a mn cos mx cos ny + b mn sin mx sin ny + c mn cos mx sin ny + d mn sin mx cos ny } There є 00 = ¼, є m0 = ½ ( m > 0), є 0n = ½ ( n > 0), є ms = 1 ( m > 0, n >0). the coefficients are given by the formulæ a mn = 1/ π 2 ∫ π -π ∫ π -π f ( x, y ) cos mx cos ny dx dy , obtained by term-by-term integration, as in an ordinary Fourier series. Ti sum of a finite number of terms of the series may also be found as in the ordinary theory. Thus ∫ ms = Σ m μ = 0 Σ n v = 0 A μ v = 1/π 2 ∫ π -π ∫ π -π f (s, t) sin( m +½) ( s - x ) sin ( n + ½) ( t - y )/2 sin ½ ( s - x ) 2 sin ½ ( t - y ) if f ( s , t ) is defined outside the original square by double periodicity, we have sub S ms = 1/π 2 ∫ π 0 ∫ π 0 f ( x + s , y + t ) + f ( x + s , y - t ) + f ( x - s , y + t ) + f ( x - s , y - t ) sin ( m + ½) s / 2 sin ½ s sin ( n + ½) t / 2 sin ½ t ds dt .


2018 ◽  
Vol 2018 ◽  
pp. 1-8
Author(s):  
Salem Neily ◽  
Sami Dhouibi ◽  
Roland Bonnet

Inclined threading dislocations (TDs) piercing the oriented free surface of a crystal are currently observed after growth of oriented thin films on substrates. Up to date the unique way to treat their anisotropic elastic properties nearby the free surface region is to use the integral formalism, which assumes no dislocation core size and needs numerical double integrations. In a first stage of the work, a new and alternative approach to the integral formalism is developed using double Fourier series and the concept of a finite core size, which is often observed in high-resolution transmission electron microscopy. In a second stage, the integral formalism and the Fourier series approaches are applied to the important case of a TD piercing the basal free surface of a hexagonal crystal. For this particular geometry, easy-to-use expressions are derived and compared to a third approach previously known for a plate-like crystal. Finally, the numerical interest and the convergence of these approaches are tested using the basal free surface of the GaN compound, in particular for TDs with Burgers vectors c and (a + c).


1977 ◽  
Vol 4 (4) ◽  
pp. 405-411
Author(s):  
A. Farah ◽  
I. M. Ibrahim ◽  
R. Green

Formulations are presented for enhancing the serviceability of one-way floor systems subjected to dynamic loading through the use of constrained viscoelastic layers. The constrained layers are combined with the floors and the resulting systems are analyzed as sandwich structures using a double Fourier series approach. Results indicate that the damping of the resulting sandwich beams is governed by factors related to the elastic and geometric properties of the constrained layer and the constraining system (i.e. cover plates and beams) and the loss modulus of the viscoelastic material, and is highly influenced by the location of the viscoelastic layer in the sandwich beams. Optimum designs of the sandwich beams are obtained using the box algorithm optimization technique.


Sign in / Sign up

Export Citation Format

Share Document