• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.1014-1017
• /
• 2007

Sound focusing technologies has been studied for various purposes from early 1990s. As a result, these technologies make us possible to apply in many uses. For example, we can treat tumors using focused ultrasonic waves without surgical knife and communicate in the ocean using time reversal array. Also applications for personal audio system become issues. Recently, as technologies are developing, in some applications, needs for regional focusing become increasing because previously suggested focusing methods, such as phase conjugation, time reversal and inverse filtering, were all about a point focusing. Therefore, studies on regional focusing method are essentially needed. Regional focusing method was firstly mentioned by Choi and Kim in 2002: acoustic contrast control. However, in regional focusing, physical interpretations between control variables and results are still not easy because of its complexity. In this regard, we tried to understand the relations between control variables and results in wavenumber domain and suggested a solution method for regional focusing: wavenumber spectrum matching. We also showed how to make spatially focused sound shape using the suggested method from the simplest case: line focusing.

• Journal of the Korean Society for Railway
• /
• v.18 no.4
• /
• pp.289-300
• /
• 2015

An important source of noise from railways is rolling noise caused by wheel and rail vibrations induced by acoustic roughness at the wheel-rail contact. In conventional approaches to predicting rail noise, the rail is regarded as placed in a free space so that the reflection from the ground is not included. However, in order to predict rail noise close to the rail, the effect of the ground should be contained in the analysis. In this study the rail noise reflected from the ground is investigated using the wavenumber domain finite element and boundary element methods. First, two rail models, one using rail attached to the rigid ground and one using rail located above rigid ground, are considered and examined to determine the rigid ground effect in terms of the radiation efficiency. From this analysis, it was found that the two models give considerably different results, so that the distance between the rail and the ground is an important factor. Second, an impedance condition was set for the ground and the effect of the ground impedance on the rail noise was evaluated for the two rail models.

• Journal of the Korean earth science society
• /
• v.38 no.3
• /
• pp.173-181
• /
• 2017

The stability of the steady Rossby-Haurwitz wave (R-H wave) in the nondivergent barotropic model (NBM) on the sphere was investigated with the normal mode method. The linearized NBM equation with respect to the R-H wave was formulated into the eigenvalue-eigenvector problem consisting of the huge sparse matrix by expanding the variables with the spherical harmonic functions. It was shown that the definite threshold R-H wave amplitude for instability could be obtained by the normal mode method. It was revealed that some unstable modes were stationary, which tend to amplify without the time change of the spatial structure. The maximum growth rate of the most unstable mode turned out to be in almost linear proportion to the R-H wave amplitude. As a whole, the growth rate of the unstable mode was found to increase with the zonal- and total-wavenumber. The most unstable mode turned out to consist of more-than-one zonal wavenumber, and in some cases, the mode exhibited a discontinuity over the local domain of weak or vanishing flow. The normal mode method developed here could be readily extended to the basic state comprised of multiple zonalwavenumber components as far as the same total wavenumber is given.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.12 s.105
• /
• pp.1378-1388
• /
• 2005

Spatial control of sound is essential to deliver better sound to the listener's position in space. As it can be experienced in many listening environments. the quality of sound can not be manifested over every Position in a hall. This motivates us to control sound in a region we select. The primary focus of the developed method has to do with the brightness and contrast of acoustic image in space. In particular, the acoustic brightness control seeks a way to increase loudness of sound over a chosen area, and the contrast control aims to enhance loudness difference between two neighboring regions. This enables us to make two different kinds of zone - the zone of quiet and the zone of loud sound - at the same time. The other perspective of this study is on the direction of sound. It is shown that we can control the direction of perceived sound source by focusing acoustic energy in wavenumber domain. To begin with, the proposed approaches are formulated for pure-tone case. Then the control methods are extended to a more general case, where the excitation signal has broadband spectrum. In order to control the broadband signal in time domain, an inverse filter design problem is defined and solved in frequency domain. Numerical and experimental results obtained in various conditions certainly validate that the acoustic brightness, acoustic contrast, direction of wave front can be manipulated for some finite region in space and time.

• Structural Engineering and Mechanics
• /
• v.75 no.4
• /
• pp.425-434
• /
• 2020

Practical non-synoptic fluctuating wind often exhibits nonstationary features and should be modeled as nonstationary random processes. Generally, the coherence function of the fluctuating wind field has time-varying characteristics. Some studies have shown that there is a big difference between the fluctuating wind field of the coherent function model with and without time variability. Therefore, it is of significance to simulate nonstationary fluctuating wind field with time-varying coherent function. However, current studies on the numerical simulation of nonstationary fluctuating wind field with time-varying coherence are very limited, and the proposed approaches are usually based on the traditional spectral representation method with low simulation efficiency. Especially, for the simulation of multi-variable wind field of large span structures such as transmission tower-line, not only the simulation is inefficient but also the matrix decomposition may have singularity problem. In this paper, it is proposed to conduct the numerical simulation of nonstationary fluctuating wind field in one-spatial dimension with time-varying coherence based on the wavenumber-frequency spectrum. The simulated multivariable nonstationary wind field with time-varying coherence is transformed into one-dimensional nonstationary random waves in the simulated spatial domain, and the simulation by wavenumber frequency spectrum is derived. So, the proposed simulation method can avoid the complicated Cholesky decomposition. Then, the proper orthogonal decomposition is employed to decompose the time-space dependent evolutionary power spectral density and the Fourier transform of time-varying coherent function, simultaneously, so that the two-dimensional Fast Fourier transform can be applied to further improve the simulation efficiency. Finally, the proposed method is applied to simulate the longitudinal nonstationary fluctuating wind velocity field along the transmission line to illustrate its performances.

• Journal of the Optical Society of Korea
• /
• v.19 no.1
• /
• pp.55-62
• /
• 2015

In this study we demonstrate ultrahigh-resolution spectral domain optical coherence tomography (UHR SD-OCT) with a linear-wavenumber (k) spectrometer, to accelerate signal processing and to display two-dimensional (2-D) images in real time. First, we performed a numerical simulation to find the optimal parameters for the linear-k spectrometer to achieve ultrahigh axial resolution, such as the number of grooves in a grating, the material for a dispersive prism, and the rotational angle between the grating and the dispersive prism. We found that a grating with 1200 grooves and an F2 equilateral prism at a rotational angle of $26.07^{\circ}$, in combination with a lens of focal length 85.1 mm, are suitable for UHR SD-OCT with the imaging depth range (limited by spectrometer resolution) set at 2.0 mm. As guided by the simulation results, we constructed the linear-k spectrometer needed to implement a UHR SD-OCT. The actual imaging depth range was measured to be approximately 2.1 mm, and axial resolution of $3.8{\mu}m$ in air was achieved, corresponding to $2.8{\mu}m$ in tissue (n = 1.35). The sensitivity was -91 dB with -10 dB roll-off at 1.5 mm depth. We demonstrated a 128.2 fps acquisition rate for OCT images with 800 lines/frame, by taking advantage of NVIDIA's compute unified device architecture (CUDA) technology, which allowed for real-time signal processing compatible with the speed of the spectrometer's data acquisition.

• Geophysics and Geophysical Exploration
• /
• v.21 no.1
• /
• pp.1-7
• /
• 2018

In the two-dimensional (2D) modeling of electrical method, the potential in the space domain is reconstructed with the calculated potentials in the wavenumber domain using inverse Fourier transform. The inverse Fourier integral is numerically evaluated using the transformed potential at different wavenumbers. In order to improve the precision of the integration, either the logarithmic or exponential approximation has been used depending on the size of wavenumber. Two numerical methods have been generally used to evaluate the integral; interval integration and Gaussian quadrature. However, both methods do not consider the distance from the current source. Thus the resulting potential in the space domain shows some error. Especially when the distance from the current source is very small or large, the error increases abruptly and the evaluated potential becomes extremely unstable. In this study, we developed a new method to calculate the integral accurately by introducing the distance from the current source to the rescaled Gauss abscissa and weight. The numerical tests for homogeneous half-space model show that the developed method can yield the error level lower than 0.4 percent over the various distances from the current source.

• Geophysics and Geophysical Exploration
• /
• v.7 no.1
• /
• pp.19-24
• /
• 2004

Frequency-space domain full-wave tomography is a promising technique for delineating detailed subsurface structure with high resolution. However, this method requires criteria for the selection of a set of optimal temporal frequency components, to achieve stability in the sequence of inversion processes together with computational efficiency. We propose a method of selecting optimal temporal frequencies, based on wavenumber continuity. The proposed method is tested numerically and is shown to be able to select an optimal set of frequency components that are sufficient to image the anomalies.

• The Korean Journal of Petroleum Geology
• /
• v.1 no.1 s.1
• /
• pp.47-52
• /
• 1993

A method and results of computations are presented for the 2-D seismic migration process in the frequency-wavenumber domain for the laterally and vertically inhomogeneous medium. In order to take the intrinsic attenuation effect into account in the migration process the complex-valued wave velocity is used in the wavefield extrapolation operator, improving the generalized frequency-wavenumber migration technique. The imaginary part of the complex-valued wave velocity includes the seismic quality factor Q value. In derivation of the solution of the wave equation for the medium of inhomogeneous wave velocity and anelasticity, the inhomogeneous medium is mathematically converted to an equivalent system which consists of a homogeneous medium of averaged slowness and an inhomogeneous distribution of hypothetical wave source. The strength of the hypothetical wave source depends on the deviation of squared slowness from the averaged value of the medium. Results of numerical computation using the technique show more distinct geologic images than those using the convensional generalized frequency-wavenumber migration. Especially, the obscured images due to the wave attenuation by anelasticity are restored to show sharp boundaries of structures. The method will be useful in the imaging of the reflection data obtained in the regions of possible petroleum or natural gas reservoir and of fractured zone.

• 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.