• Journal of the Korean earth science society
• /
• v.34 no.5
• /
• pp.393-401
• /
• 2013

Spherical harmonics power spectrum of the geopotential field of Gaussian-bell type on the sphere was investigated using integral formula that is associated with Legendre polynomials. The geopotential field of Gaussian-bell type is defined as a function of sine of angular distance from the bell's center in order to guarantee the continuity on the global domain. Since the integral-formula associated with the Legendre polynomials was represented with infinite series of polynomial, an estimation method was developed to make the procedure computationally efficient while preserving the accuracy. The spherical harmonics power spectrum was shown to vary significantly depending on the scale parameter of the Gaussian bell. Due to the accurate procedure of the new method, the power (degree variance) spanning over orders that were far higher than machine roundoff was well explored. When the scale parameter (or width) of the Gaussian bell is large, the spectrum drops sharply with the total wavenumber. On the other hand, in case of small scale parameter the spectrum tends to be flat, showing very slow decaying with the total wavenumber. The accuracy of the new method was compared with theoretical values for various scale parameters. The new method was found advantageous over discrete numerical methods, such as Gaussian quadrature and Fourier method, in that it can produce the power spectrum with accuracy and computational efficiency for all range of total wavenumber. The results of present study help to determine the allowable maximum scale parameter of the geopotential field when a Gaussian-bell type is adopted as a localized function.

• Journal of the Korean earth science society
• /
• v.41 no.4
• /
• pp.344-355
• /
• 2020

Tropical cyclone scale vortices and associated Rossby waves were investigated numerically using high-resolution barotropic models on the global domain. The equations of the barotropic model were discretized using the spectral transform method with the spherical harmonics function as orthogonal basis. The initial condition of the vortex was specified as an axisymmetric flow in the gradient wind balance, and four types of basic zonal states were employed. Vortex tracks showed similar patterns as those on the beta-plane but exhibited more eastward displacement as they moved northward. The zonal-mean flow appeared to control not only the west-east translation but also the meridional translation of the vortex. Such a meridional influence was revealed to be associated with the beta gyre and the Rossby wave, which are formed around the vortex due to the beta effect. In the case of the basic zonal state of climatological mean, the meridional translation speed reached the maximum value when the vortex underwent recurving.

• Atmosphere
• /
• v.33 no.4
• /
• pp.387-398
• /
• 2023

The double Fourier series (DFS) spectral dynamical core is evaluated for the two idealized test cases in comparison with the spherical harmonics (SPH) spectral dynamical core. A new approach in calculating the meridional expansion coefficients of DFS, which was recently developed to alleviate a computational error but only applied to the 2D spherical shallow water equation, is also tested. In the 3D deformational tracer transport test, the difference is not conspicuous between SPH and DFS simulations, with a slight outperformance of the new DFS approach in terms of undershooting problem. In the baroclinic wave development test, the DFS-simulated wave pattern is quantitatively similar to the SPH-simulated one at high resolutions, but with a substantially lower computational cost. The new DFS approach does not offer a salient advantage compared to the original DFS while computation cost slightly increases. This result suggests that the current DFS spectral method can be a practical and alternative dynamical core for high-resolution global modeling.

• Journal of the Korean earth science society
• /
• v.37 no.2
• /
• pp.107-116
• /
• 2016

Stability of the barotropic Rossby-Haurwitz wave is investigated using the numerical models on the global domain. The Rossby-Haurwitz wave under investigation is composed of the basic zonal flow of super-rotation and a finite amplitude spherical harmonic wave. The Rossby-Haurwitz wave is given as either steady or unsteady wave by adjusting the strength of the super-rotating zonal flow. Stability as well as the growth rate of the wave in the numerical simulation is determined by comparing the perturbation amplitude at two different time stages. Unstable modes of the Rossby-Haurwitz wave exhibited a horizontal structure composing of various zonal-wavenumber components. The vorticity perturbation for some modes showed a discontinuity around the area of weak flow, which was found robust regardless of the horizontal resolution of the model. Fourier finite element model was shown to generate the unstable mode in earlier stage of the time integration due to less accuracy compared to the spherical harmonic spectral model. Taking the overall accuracy of the models into consideration, the time by which the unstable mode begin to dominate over the spherical harmonic wave was estimated.

• Journal of the Korean earth science society
• /
• v.28 no.7
• /
• pp.936-946
• /
• 2007

Analytical geopotential field in balance with the sectoral mode (the first symmetric mode with respect to the equator) of the Rossby-Haurwitz wave on the inclined rotation axis was derived in presence of superrotation background flow. The balanced field was obtained by inverting the divergence equation with the time derivative being zero. The inversion consists of two steps, i.e., the evaluation of nonlinear forcing terms and the finding of analytical solutions based on the Poisson's equation. In the second step, the forcing terms in the from of Legendre function were readily inverted due to the fact that Legendre function is the eigenfunction of the spherical Laplacian operator, while other terms were solved either by introducing a trial function or by integrating the Legendre equation. The balanced field was found to be expressed with six zonal wavenumber components, and shown to be of asymmetric structure about the equator. In association with asymmetricity, the advantageous point of the balanced field as a validation method for the numerical model was addressed. In special cases where the strength of the background flow is a half of or exactly the same as the rotation rate of the Earth it was revealed that one of the zonal wavenumber components vanishes. The analytical balanced field was compared with the geopotential field which was obtained using a spherical harmonics spectral model. It was found that the normalized difference lied in the order of machine rounding, indicating the reliability of the analytical results. The stability of the sectoral mode of Rossby-Haurwitz wave and the associated balanced field was discussed, comparing with the flrst antisymmetric mode.

• Investigative Magnetic Resonance Imaging
• /
• v.10 no.2
• /
• pp.108-116
• /
• 2006

Purpose : High-resolution spiral-scan imaging is performed at 3 Tesla MRI system. Since the gradient waveforms for the spiral-scan imaging have lower slopes than those for the Echo Planar Imaging (EPI), they can be implemented with the gradient systems having lower slew rates. The spiral-scan imaging also involves less eddy currents due to the smooth gradient waveforms. The spiral-scan imaging method does not suffer from high specific absorption rate (SAR), which is one of the main obstacles in high field imaging for rf echo-based fast imaging methods such as fast spin echo techniques. Thus, the spiral-scan imaging has a great potential for the high-speed imaging in high magnetic fields. In this paper, we presented various high-resolution images obtained by the spiral-scan methods at 3T MRI system for various applications. Materials and Methods : High-resolution spiral-scan imaging technique is implemented at 3T whole body MRI system. An efficient and fast higher-order shimming technique is developed to reduce the inhomogeneity, and the single-shot and interleaved spiral-scan imaging methods are developed. Spin-echo and gradient-echo based spiral-scan imaging methods are implemented, and image contrast and signal-tonoise ratio are controlled by the echo time, repetition time, and the rf flip angles. Results : Spiral-scan images having various resolutions are obtained at 3T MRI system. Since the absolute magnitude of the inhomogeneity is increasing in higher magnetic fields, higher order shimming to reduce the inhomogeneity becomes more important. A fast shimming technique in which axial, sagittal, and coronal sectional inhomogeneity maps are obtained in one scan is developed, and the shimming method based on the analysis of spherical harmonics of the inhomogeneity map is applied. For phantom and invivo head imaging, image matrix size of about $100{\times}100$ is obtained by a single-shot spiral-scan imaging, and a matrix size of $256{\times}256$ is obtained by the interleaved spiral-scan imaging with the number of interleaves of from 6 to 12. Conclusion : High field imaging becomes increasingly important due to the improved signal-to-noise ratio, larger spectral separation, and the higher BOLD-based contrast. The increasing SAR is, however, a limiting factor in high field imaging. Since the spiral-scan imaging has a very low SAR, and lower hardware requirements for the implementation of the technique compared to EPI, it is suitable for a rapid imaging in high fields. In this paper, the spiral-scan imaging with various resolutions from $100{\times}100$ to $256{\times}256$ by controlling the number of interleaves are developed for the high-speed imaging in high magnetic fields.