• 한국해양환경ㆍ에너지학회지
• /
• 제13권3호
• /
• pp.174-180
• /
• 2010

변동 수심에서의 파랑변형을 비균질 Helmholtz 방정식을 이용하여 계산하였다. 포텐셜 함수가 존재한다고 가정하였으며, 변수분리를 적용하였다. 본 논문에서는 조화파만을 고려하였다. 포텐셜 함수로 구성된 지배방정식을 정수면에 직접 적용하였고, 변동 수심에 대한 비균질 Helmholtz 방정식을 얻었다. 파랑의 진폭과 위상차로 얻어진 복합 포텐셜 함수의 지배방정식을 실수형 변수로 된 두 방정식으로 분리하였다. 분리된 방정식들은 각각 1차와 2차 상미분 방정식이며, 이 방정식들을 단순한 형태의 중앙차분 수치기법을 이용하여 차분식으로 변형하였다. 측면 경계조건에서의 파랑의 진폭, 진폭경사, 그리고 위상경사를 경계면에 적용하여 전방진행방법으로 전 영역에서 해를 구하였다 Booij의 경사면 있는 저면의 경우와 Bragg의 물결모양이 있는 저면의 경우에 적용하였다. 본 연구로 도출된 비균질 Helmholtz 방정식은 완전 선형방정식 계산 결과, Massel의 수정 완경사 방정식, 그리고 Berkhoff의 완경사 방정식의 적용 결과와 비교하였으며, 만족스러운 결과를 얻었다.

• 한국연소학회지
• /
• 제22권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
• /
• 제3권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).

• 대한수학회논문집
• /
• 제16권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.

• 전자공학회논문지D
• /
• 제35D권6호
• /
• pp.8-14
• /
• 1998

Wiener-Hopf 적분방정식으로부터 경계면 위의 전체파를 미지수로 하는 파수영역에서의 쌍적분 방정식을 얻는 기존의 유도과정을 검토하였다. 이러한 기존의 유도 과정은 결국 Wiener-Hopf 적분방정식으로부터 Helmholtz-Kirchhoff 적분방정식을 유도하는 과정임을 해석적으로 보였다.

• 소음진동
• /
• 제7권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
• /
• 제22권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.

• 호남수학학술지
• /
• 제41권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.

• 한국음향학회지
• /
• 제27권8호
• /
• pp.411-417
• /
• 2008

본 논문에서는 헬름홀쯔 적분 방정식에서 유도된 식을 이용하여 구조물의 표면 압력을 구조진동 성분에 대한 단순한 적분형태로 표현하여 음향방사 및 구조/음향 연성 문제를 수치적으로 푸는 방법에 대하여 다룬다. 이 식은 임의의 형상에 대하여 유도된 식으로 Rayleigh 식과 유사한 형태를 갖는다. 이 식을 이용하면 표면 압력을 구조물의 속도에 대한 단순 적분 형태로 나타낼 수 있기 때문에 경계요소법과 같이 연립방정식에 대한 행렬식을 풀 필요가 없다. 또한 헬름홀쯔 적분 방정식에 기반을 둔 다른 방법 들이 가지는 해의 유일성 문제도 갖지 않는 장점이 있다. 본 논문에서는 구형 셀에 대하여 수치해와 정해를 비교하여 제안한 방법의 타당성을 검증하였다.

• 대한전기학회논문지
• /
• 제41권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.