• 한국분무공학회지
• /
• 제29권1호
• /
• pp.23-31
• /
• 2024

This study conducted a 3D thermo-acoustic analysis based on the helmholtz solver to analyze the major resonance modes causing combustion instability in a single-can combustor. The experimental investigations were carried out on a test rig designed by the Korea Institute of Machinery & Materials (KIMM) under various conditions of hydrogen co-firing and fuel staging. Through these experiments, two primary unstable frequencies were identified. To determine the resonance modes of these frequencies, a 3D thermo-acoustic analysis was conducted using temperature information from the test rig. The results confirmed that the unstable frequencies observed in the experiments were all longitudinal modes. Additionally, the mode shapes identified in the analysis facilitated a simplification of the exit geometry for the low-order network model, confirming that this did not significantly affect the fundamental resonance modes.

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

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2017년도 제48회 춘계학술대회논문집
• /
• pp.68-77
• /
• 2017

고주파 연소불안정은 비정상 화염의 열방출율 섭동과 연소실 내부에서 공진되는 음향파의 상호 결합으로 발생하는 열음향 문제로, 다양한 해석적 접근방법이 존재한다. 본 논문에서는 주파수 영역에서 선형음향 가정과 시간지연 이론을 이용한 3차원 FEM Helmholtz solver의 개발 사례를 소개하였으며, 가변길이 희박 예혼합 연소기의 자발 연소불안정 예측과 수동제어기구(배플, 음향공진기)의 설계분석 결과를 제시하였다. 또한 시간 영역에서 시간지연 이론을 이용한 압축성 유동 해석코드를 통해, 고진폭 압력섭동에 의해 야기되는 비선형 음향 특성과 한계사이클 현상을 분석하였다.

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

• 한국추진공학회지
• /
• 제19권4호
• /
• pp.15-23
• /
• 2015

본 연구는 유한요소해석 기반의 Helmholtz solver인 ASCI3D를 사용하여 예혼합 가스터빈 연소기에서 연소불안정 현상과 관련된 주요 결과들을 예측하였다. 해석결과는 실험결과와 비교할 때 일반적인 불안정 특성 예측에는 성공하였으나, 불안정 구간이 나타나는 영역을 다소 과대 예측하는 경향이 발견되었다. 이를 개선하기 위한 노력으로 반사계수가 불안정 특성에 미치는 영향이 분석되었다. 결과로부터 반사계수는 불안정 특성에 큰 영향을 미치고, 정확한 예측 결과를 얻기 위해서는 정확한 음향 경계 조건의 정의가 필수적인 것으로 나타났다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제2권1호
• /
• pp.27-39
• /
• 1998

We describe a fast numerical method for non-separable elliptic equations in self-adjoin form on irregular adaptive domains. One of the most successful results in numerical PDE is developing rapid elliptic solvers for separable EPDEs, for example, Fourier transformation methods for Poisson problem on a square, however, it is known that there is no rapid elliptic solvers capable of solving a general nonseparable problems. It is the purpose of this paper to present an iterative solver for linear EPDEs in self-adjoint form. The scheme discussed in this paper solves a given non-separable equation using a sequence of solutions of Poisson equations, therefore, the most important key for such a method is having a good Poison solver. High performance is achieved by using a fast high-order adaptive Poisson solver which requires only about 500 floating point operations per gridpoint in order to obtain machine precision for both the computed solution and its partial derivatives. A few numerical examples have been presented.

• 한국항공우주학회지
• /
• 제38권5호
• /
• pp.445-455
• /
• 2010

본 연구에서는 실제 로켓엔진 및 가스터빈용 연소기 내부의 열음향 불안정을 효과적으로 예측하기 위하여, 헬름홀츠 방정식과 시간지연모델을 이용한 3차원 유한요소법 해석코드를 개발하였다. 연소응답항에 의해 수치적으로 야기되는 비선형성은 반복법으로 선형화 하였으며, Arnoldi 방법을 사용하여 대용량 고유치 문제를 해석하였다. 해석결과인 복소각주파수와 음향 압력장을 통해 각 음향모드의 공진주파수, 진폭의 증폭/감쇠 여부 그리고 모드 형태를 예측할 수 있다. 이론해가 존재하는 두 가지 문제를 통해 출구 임피던스와 예혼합 화염이 종 방향 음향장에 미치는 영향에 대한 예측 정확도를 평가하였으며, 배플 유무에 따른 횡 방향 음향 모드의 주파수 변이를 상온 음향시험 결과와 비교/검증하였다.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제18권4호
• /
• pp.305-316
• /
• 2014

The Hodge-Helmholtz decomposition splits a vector field into the unique sum of a divergence-free vector field (solenoidal part) and a gradient field (irrotational part). In a bounded domain, a boundary condition needs to be supplied to the decomposition. The decomposition with the non-penetration boundary condition is equivalent to solving the Poisson equation with the Neumann boundary condition. The Gibou-Min method is an application of the Poisson solver by Purvis and Burkhalter to the decomposition. Using the $L^2$-orthogonality between the error vector and the consistency, the convergence for approximating the divergence-free vector field was recently proved to be $O(h^{1.5})$ with step size h. In this work, we analyze the convergence of the irrotattional in the decomposition. To the end, we introduce a discrete version of the Poincare inequality, which leads to a proof of the O(h) convergence for the scalar variable of the gradient field in a domain with general intersection property.

• 한국분무공학회지
• /
• 제25권3호
• /
• pp.119-125
• /
• 2020

In this paper, we used Helmholtz solver based on 3D finite element method to quantitatively analyze the effects of change of gain, time delay and time delay spread, which are the main variables of flame transfer function, on combustion instability in gas turbine combustor. The effects of the variable of flame transfer function on the frequency and growth rate, which are the main results of combustion instability, were analyzed by applying the conventional heat release fluctuation model and modified one considering the time spread. The analysis results showed that the change of gain and time delay in the same resonance mode affected the frequency of the given resonance modes as well as growth rate of the feedback instability, however, the effect of time delay spread was not relatively remarkable, compared with the dominant effect of time delay.