• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.3
• /
• pp.174-180
• /
• 2010

The inhomogeneous Helmholtz equation is introduced for variable water depth and potential function and separation of variables are introduced for the derivation. Only harmonic wave motions are considered. The governing equation composed of the potential function for irrotational flow is directly applied to the still water level, and the inhomogeneous Helmholtz equation for variable water depth is obtained. By introducing the wave amplitude and wave phase gradient the governing equation with complex potential function is transformed into two equations of real variables. The transformed equations are the first and second-order ordinary differential equations, respectively, and can be solved in a forward marching manner when proper boundary values are supplied, i.e. the wave amplitude, the wave amplitude gradient, and the wave phase gradient at a side boundary. Simple spatially-centered finite difference numerical schemes are adopted to solve the present set of equations. The equation set is applied to two test cases, Booij’ inclined plane slope profile, and Bragg’ wavy bed profile. The present equations set is satisfactorily verified against other theories including the full linear equation, Massel's modified mild-slope equation, and Berkhoff's mild-slope equation etc.

• Journal of the Korean Society of Combustion
• /
• v.22 no.3
• /
• pp.41-46
• /
• 2017

This study has numerically investigated the flame-acoustics interactions in the turbulent partially premixed flame field. In the present approach, in order to analyze the combustion instability, the present approach has employed the LES-based combustion model as well as the Helmholtz solver. Computations are made for the validation case of the partially premixed LIMOUSINE burner. In terms of the FFT data, numerical results are compared with experimental data. Moreover, Helmholtz equation in frequency domain is solved by combining CFD field data including the flight time from a nozzle to the flame zone. Based on numerical results, the detailed discussions are made for the essential features of the combustion instability encountered in the partially premixed burner.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.1-4
• /
• 1999

The Helmholtz equation is very important in physics and engineering. However, solution of the Helmholtz equation is in general known as a very difficult phenomenon. For if the ${\omega}$ is negative, the FDM discretized linear system becomes indefinite, whose solution by iterative method requires a very clever preconditioner. In this paper we assume that ${\omega}$ is nonnegative, and determine the optimal ${\rho}$ parameter for the three dimensional ADI iteration for the Helmholtz equation. The ADI(Alternating Direction Implicit) method is also getting new attentions due to the fact that it is very suitable to the vector/parallel computers, for example, as a preconditioner to the Krylov subspace methods. However, classical ADI was developed for two dimensions, and for three dimensions it is known that its convergence behaviour is quite different from that in two dimensions. So far, in three dimensions the so-called Douglas-Rachford form of ADI was developed. It is known to converge for a relatively wide range of ${\rho}$ values but its convergence is very slow. In this paper we determine the necessary conditions of the ${\rho}$ parameter for the convergence and optimal ${\rho}$ for the three dimensional ADI iteration of the Peaceman-Rachford form for the real Helmholtz equation with nonnegative ${\omega}$. Also, we conducted some experiments which is in close agreement with our theory. This straightforward extension of Peaceman-rachford ADI into three dimensions will be useful as an iterative solver itself or as a preconditioner to the the Krylov subspace methods, such as CG(Conjugate Gradient) method or GMRES(m).

• Communications of the Korean Mathematical Society
• /
• v.16 no.3
• /
• pp.427-436
• /
• 2001

It is the purpose of this paper to study the reflection of solutions of Helmholtz equation with Neumann boundary data. In detail let u be a solution of Helmholtz equation in the exterior of a ball in R$^3$ with exterior Neumann data ∂(sub)νu = 0 on the boundary of the ball. We prove that u can be extended to R$^3$ except the center of the ball. As a corollary, we prove that a sound hard ball can be identified by the scattering amplitude corresponding to a single incident direction and as single frequency.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.35D no.6
• /
• pp.8-14
• /
• 1998

The derivation of the dual integral equation in the spectral domain, which has total fields of the interfaces as unknowns, is investigated. It is analytically shown that the derivation of the dual integral equation is equivalent to deriving the Helmholtz-Kirchhoff integral equation from the Wiener-Hopf integral equation.

• Journal of KSNVE
• /
• v.7 no.6
• /
• pp.975-984
• /
• 1997

A direct boundary element method (DBEM) is developed for thin bodies whose surfaces are rigid or compliant. The Helmholtz integral equation and its normal derivative integral equation are adoped simultaneously to calculate the pressure on both sides of the thin body, instead of the jump values across it, to account for the different surface conditions of each side. Unlike the usual assumption, the normal velocity is assumed to be discontinuous across the thin body. In this approach, only the neutral surface of the thin body has to be discretized. The method is validated by comparison with analytic and/or numerical results for acoustic scattering and radiation from several surface conditions of the thin body; the surfaces are rigid when stationary or vibrating, and part of the interior surface is lined with a sound-absoring material.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.2
• /
• pp.101-113
• /
• 2018

The dual singular function method(DSFM) is a numerical algorithm to get optimal solution including corner singularities for Poisson and Helmholtz equations. In this paper, we apply DSFM to solve heat equation which is a time dependent problem. Since the DSFM for heat equation is based on DSFM for Helmholtz equation, it also need to use Sherman-Morrison formula. This formula requires linear solver n + 1 times for elliptic problems on a domain including n reentrant corners. However, the DSFM for heat equation needs to pay only linear solver once per each time iteration to standard numerical method and perform optimal numerical accuracy for corner singularity problems. Because the Sherman-Morrison formula is rather complicated to apply computation, we introduce a simplified formula by reanalyzing the Sherman-Morrison method.

• Honam Mathematical Journal
• /
• v.41 no.2
• /
• pp.403-420
• /
• 2019

The Helmholtz equation represents acoustic or electromagnetic scattering phenomena. The Method of Lines are known to have many advantages in simulation of forward and inverse scattering problems due to the usage of angle rays and Bessel functions. However, the method does not account for the jump phenomena on obstacle boundary and the approximation includes many high order Bessel functions. The high order Bessel functions have extreme blow-up or die-out features in resonance region obstacle boundary. Therefore, in particular, when we consider shape reconstruction problems, the method is suffered from severe instabilities due to the logical confliction and the severe singularities of high order Bessel functions. In this paper, two approximation formulas for the Helmholtz equation are introduced. The formulas are new and powerful. The derivation is based on Method of Lines, Huygen's principle, boundary jump relations, Addition Formula, and the orthogonality of the trigonometric functions. The formulas reduce the approximation dimension significantly so that only lower order Bessel functions are required. They overcome the severe instability near the obstacle boundary and reduce the computational time significantly. The convergence is exponential. The formulas adopt the scattering jump phenomena on the boundary, and separate the boundary information from the measured scattered fields. Thus, the sensitivities of the scattered fields caused by the boundary changes can be analyzed easily. Several numerical experiments are performed. The results show the superiority of the proposed formulas in accuracy, efficiency, and stability.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.8
• /
• pp.411-417
• /
• 2008

An alternative formulation of the Helmholtz integral equation derived to express the pressure field explicitly in terms of the velocity vector of a radiating surface is used to solve acoustic radiation and fluid/structure interaction problems. This formulation, derived for arbitrary sources, is similar in form to the Rayleigh's formula for planar sources. Because the surface pressure field is expressed explicitly as a surface integral of the surface velocity, which can be implemented numerically using standard Gaussian quadratures, there is no need to use BEM to solve a set of simultaneous equations for the surface pressure at the discretized nodes. Furthermore the non-uniqueness problem inherent in methods based on Helmholtz integral equation is avoided. Validation of this formulation is demonstrated for some simple geometries.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.41 no.1
• /
• pp.89-94
• /
• 1992

In this paper, the static limit of Helmholtz equation is discussed for the analysis of wave scattering in a wave scattering in a waveguide. Boundary integral equation method is used to formulato the scattering process in the exterior of the scatterer and finite element method in the interior of the scatterer. And hybrid ray-mode method is used to provide the Green's function in the waveguide. The proposed algorithm is applied algarithm is applied to a sample problem with arbitrary scatterer in a waveguide. The results are compared with those of static analysis.