• Journal of ILASS-Korea
• /
• v.29 no.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
• /
• 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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.68-77
• /
• 2017

High-frequency combustion instability results from a feedback coupling between the unsteady heat release rate and the acoustic waves formed resonantly in the combustion chamber. It can be modeled as thermoacoustic problems with various degrees of the assumptions and simplifications. This paper presents numerical analysis of self-excited combustion instabilities in a variable-length lean-premixed combustor and designs of passive control devices such as baffle and acoustic resonators in a framework of 3-D FEM Helmholtz solver. Nonlinear behaviors such as steep-fronted shock waves and a finite amplitude limit cycle are also investigated with a compressible flow simulation technique.

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

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.4
• /
• pp.15-23
• /
• 2015

This study predicts the basic characteristics of combustion instabilities in a gas turbine lean premixed combustor using ASCI3D code which is a FEM(Finite Element Method)-based Helmholtz solver. The prediction results show the good agreement with the measured data in modeling the overall combustion instability features, however, the code is found to overpredict the unstable conditions. As one of the efforts to improve the model accuracy, the effects of acoustic boundary conditions on the instability growth rate are analyzed. As a result, it is shown that the acoustic reflection coefficient has a great impact on the instability and the prediction accuracy can be enhanced by defining the precise acoustic conditions.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.5
• /
• pp.445-455
• /
• 2010

In order to effectively predict thermo-acoustic instability within real combustors of rocket engines and gas turbines, in the present study, the Helmholtz equation in conjunction with the time lag hypothesis is discretized by the finite element method on three-dimensional hybrid unstructured mesh. Numerical nonlinearity caused by the combustion response term is linearized by an iterative method, and the large-scale eigenvalue problem is solved by the Arnoldi method available in the ARPACK. As a consequence, the final solution of complex valued eigenfrequency and acoustic pressure field can be interpreted as resonant frequency, growth rate, and modal shape for acoustic modes of interest. The predictive capabilities of the present method have been validated against two academic problems with complex impedance boundary and premixed flame, as well as an ambient acoustic test for liquid rocket combustion chamber with/without baffle.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.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.

• Journal of ILASS-Korea
• /
• v.25 no.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.