• Journal of Mechanical Science and Technology
• /
• v.14 no.1
• /
• pp.103-112
• /
• 2000

Numerical simulations of unsteady opposed-flow flames are performed using an adaptive time integration method designed for differential-algebraic systems. The compressibility effect is considered in deriving the system of equations, such that the numerical difficulties associated with a high-index system are alleviated. The numerical method is implemented for systems with detailed chemical mechanisms and transport properties by utilizing the Chemkin software. Two test simulations are performeds hydrogen/air diffusion flames with an oscillatory strain rate and transient ignition of methane against heated air. Both results show that the rapid transient behavior is successfully captured by the numerical method.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.03a
• /
• pp.852-855
• /
• 2010

Numerical simulations for large scale triaxial tests with large diameter rockfill materials are conducted using distinct element method. For generation of compacted assembly with specific grain size distribution and initial material porosity, the clump logic method and expansion of generated particles are adapted. With micro parameters which are chosen by calibration process, discrete particle modelling of triaxial test in case of other confining stress and cyclic loading condition were conducted. Also numerical simulations of fluid injection into particulate materials were conducted to observe cavity initiation and propagation using distinct element method. The fluid scheme solves the continuity and Navior-Stokes equations numerically, then derives pressure and velocity vectors for fixed grid by considering the existence of particles within the fluid cell.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.428-433
• /
• 2011

In this paper, a prediction method for dynamic stability derivatives is studied using steady state simulations in rotational coordinates. The simulations require the extension of a standard CFD formulations based on inertial coordinate. A new CFD code based on the method are developed. Flows induced by steady circular motions of airfoils with a constant pitch rate are simulated with the code. From the numerical simulations, the pitch rate derivatives are obtained at various Mach numbers, and the results are compared with other numerical results. The numerical simulations show that the new code are capable of predicting dynamic stability derivatives.

• Advances in aircraft and spacecraft science
• /
• v.6 no.4
• /
• pp.333-348
• /
• 2019

Low fuel regression rate is the main drawback of hybrid rocket which should be overcome. One of the improvement techniques to this problem is usage of a solid fuel grain with a complicated geometry port, which has been promoted owing to the recent development of additive manufacturing technologies. In the design of a hybrid rocket fuel grain with a complicated geometry port, the understanding of fuel regression behavior is very important. Numerical investigations of fuel regression behavior requires a capturing method of solid fuel surface, i.e. gas-solid interface. In this study, level set method is employed as such a method and the preliminary numerical tool for capturing a hybrid rocket solid fuel surface is developed. At first, to test the adequacy of the numerical modeling, the simulation results for circular port are compared to the experimental results in open literature. The regression rates and oxidizer to fuel ratios show good agreements between the simulations and the experiments, after passing enough time. However, during the early period of combustion, there are the discrepancies between the simulations and the experiments, owing to transient phenomena. Second, the simulations of complicated geometry ports are demonstrated. In this preliminary step, a star shape is employed as complicated geometry of port. The slot number effect in star port is investigated. The regression rate decreases with increasing the slot number, except for the star port with many slots (8 slots) in the latter half of combustion. The oxidizer to fuel ratio increases with increasing the slot number.

• Geomechanics and Engineering
• /
• v.22 no.1
• /
• pp.97-108
• /
• 2020

The time-consuming and less objectivity are the main problems of conventional micromechanical parameters calibration method of Particle Flow Code simulations. Thus this study aims to address these two limitation of the conventional "trial-and-error" method. A new calibration method for the linear parallel bond model (CM-LPBM) is proposed. First, numerical simulations are conducted based on the results of the uniaxial compression tests on limestone. The macroscopic response of the numerical model agrees well with the results of the uniaxial compression tests. To reduce the number of the independent micromechanical parameters, numerical simulations are then carried out. Based on the results of the orthogonal experiments and the multi-factor variance analysis, main micromechanical parameters affecting the macro parameters of rocks are proposed. The macro-micro parameter functions are ultimately established using multiple linear regression, and the iteration correction formulas of the micromechanical parameters are obtained. To further verify the validity of the proposed method, a case study is carried out. The error between the macro mechanical response and the numerical results is less than 5%. Hence the calibration method, i.e., the CM-LPBM, is reliable for obtaining the micromechanical parameters quickly and accurately, providing reference for the calibration of micromechanical parameters.

• Communications of the Korean Mathematical Society
• /
• v.22 no.3
• /
• pp.453-464
• /
• 2007

A computationally efficient numerical scheme is presented for the phase-field model of two-phase systems for anisotropic interfacial energy. The scheme is solved by using a nonlinear multigrid method. When the coefficient for the anisotropic interfacial energy is sufficiently high, the interface of the system shows corners or missing crystallographic orientations. Numerical simulations with high and low anisotropic coefficients show excellent agreement with exact equilibrium shapes. We also present spinodal decomposition, which shows the robustness of the pro-posed scheme.

• Journal of Ship and Ocean Technology
• /
• v.5 no.4
• /
• pp.29-43
• /
• 2001

In this paper, free surface flows around an advancing two-dimensional rectangular cylinder piercing the free surface are studied using numerical and experimental methods. Especially, wave breaking phenomenon around the cylinder is treated in detail. A series of numerical simulations and experiments were performed for the purpose of comparison. For the numerical simulations, a finite difference method was adopted with a rectangular grid system, and the variation of the free surface was computed by the marker density method. The computational results are compared with the experiments. It is confirmed that the present numerical method is useful for the numerical simulation of nonlinear free surface waves around a piercing body.

• Journal of The Korean Astronomical Society
• /
• v.48 no.1
• /
• pp.67-73
• /
• 2015

We present a novel method to implement time-delayed propagation of radiation fields in cosmological radiative transfer simulations. Time-delayed propagation of radiation fields requires construction of retarded-time fields by tracking the location and lifetime of radiation sources along the corresponding light-cones. Cosmological radiative transfer simulations have, until now, ignored this "light-cone effect" or implemented ray-tracing methods that are computationally demanding. We show that radiative transfer calculation of the time-delayed fields can be easily achieved in numerical simulations when periodic boundary conditions are used, by calculating the time-discretized retarded-time Green's function using the Fast Fourier Transform (FFT) method and convolving it with the source distribution. We also present a direct application of this method to the long-range radiation field of Lyman-Werner band photons, which is important in the high-redshift astrophysics with first stars.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.215-224
• /
• 2001

We present the whole basis of numerical method and useful formulae for general relativistic magnetohydrodynamic simulations in Kerr space-time.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.206-218
• /
• 2014

The wave attenuation by floating breakwaters in high amplitude waves, which can lead to wave overtopping and breaking, is examined by numerical simulations. The governing equations, the Navier-Stokes equations and the continuity equation, are calculated in a fixed Cartesian grid system. The body boundaries are defined by the line segment connecting the points where the grid line and body surface meet. No-slip and divergence free conditions are satisfied at the body boundary cell. The nonlinear waves near the moving body is defined using the modified marker-density method. To verify the present numerical method, vortex induced vibration on an elastically mounted cylinder and free roll decay are numerically simulated and the results are compared with those reported in the literature. Using the present numerical method, the wave attenuations by three kinds of floating breakwaters are simulated numerically in a regular wave to compare the performance.