• Journal of the Korean Society of Safety
• /
• v.19 no.3 s.67
• /
• pp.124-129
• /
• 2004

To evaluate the numerical method in simulation of diffusion flames and to see the effects of strain rate and fuel concentration on the flame radius and thickness, the nonpremixed methane-air counterflow flames in microgravity were simulated axisymmetrically by using the MST Fire Dynamics Simulator (FDS). The $1000^{\circ}C$ based flame radius and thickness were investigated for the mole fraction of methane in the fuel stream, $X_m=20,\;50,\;and\;80\%$ and the global strain rates $a_g=20,\;60,\;and\;90s^{-1}$ for each mole fraction. The flame radius increased with the global strain rate while the flame thickness decreased linearly as the global strain rate increased. The flame radius decreased as the mole fraction increased, but it was not so sensitive to the mole fraction compared with the global strain rate. Since there was good agreement in the nondimensional flame thickness obtained with OPPDIF and FDS respectively, it was confirmed that FDS is capable of predicting well the counterflow flames in a wide range of strain rate and fuel concentration.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.4
• /
• pp.297-309
• /
• 2013

Design optimization is a mathematical process to find an optimal solution through the use of formal optimization algorithms. Design plays a vital role in the engineering field; therefore, using design tools in education and research is becoming more and more important. Recently, numerical design optimization in fluid mechanics, which uses computational fluid dynamics (CFD), has numerous applications in the engineering field, because of the rapid development of high-performance computing resources. However, it is difficult to find design optimization software and contents for educational purposes in aerospace engineering. In the present study, we have developed an aerodynamic design framework specifically for an airfoil, based on the EDucation-research Integration through Simulation On the Net (EDISON) portal. The airfoil design framework is composed of three subparts: a geometry kernel, CFD flow analysis, and an optimization algorithm. Through a seamless interface among the subparts, an iterative design process is conducted. In addition, the CFD flow analysis and the design framework are provided through a web-based portal system, while the computation is taken care of by a supercomputing facility. In addition to the software development, educational contents are developed for lectures associated with design optimization in aerospace and mechanical engineering education programs. The software and content developed in this study is expected to be used as a tool for e-learning material, for education and research in universities.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.4
• /
• pp.277-284
• /
• 2020

The remaining range between missile and target can be used to calculate the guidance command as well as to determine the explosion time of the warhead. Since the range, however, is not directly measured by on-board sensors of the missile, it is usually estimated by filter-based techniques using angle-only measurements. Conventional filter-based techniques are complex and require huge computation. In this paper, we propose a range estimation algorithm based on the geometrical triangulation principle for two points of missiles and a target. The proposed algorithm has a simple structure but the accuracy is largely dependent on the measurement errors. To improve the accuracy of estimation, Digital Fading Memory Filter (DFMF) is applied. The performance of the proposed algorithm is analyzed through numerical simulations.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.21 no.2
• /
• pp.495-504
• /
• 1996

In this paper, series-fed mirostrip array anteenas with coupling-slots are prposed and their operating characteristics are analyzed based on analyzed based on analytical and experimental results. An accurate analysis method for the slot-coupled feed structure is based on using both circuit coupling between anteenas and full-wave analysis which use travling wave mode and non-traveling wave mode on feed line. The basis functions that used for the numerical analysis bas been determined depending upon the accuracy, convergence properties of the solution, and the computation time:The patch uses 3 EB mode, the slot uses 1PWS mode, and feeders use 5 PWS mode. Series-fed array antennas have been designed, built, and tested in a standing-wave configuration. Using the results of the full-wave analysis, the chebyshev array antennas consisting of 8 elements are designed and fabricated changing the amount excited to each array element by adjusting slot length and by the slot position relative to the feeder. Experiment results show that the series-fed array antenna designed by adjusting the slot position relative to the feeder is superior to that designed by slot length.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.14 no.8
• /
• pp.871-877
• /
• 2003

In this paper, the design parameters for the magnetic field source at extremely low frequency are proposed. This facility can be used fur in vivo experiments with small animals to investigate biological response to the driving magnetic fields. In case that the exposed animals are motionless, the animals may be affected by the directivity of driving field. To avoid this effect, a 2-axis ELF magnetic field driving apparatus was designed. The optimum location and number of turns of each coil were obtained by numerical analysis. Applying these data to the MATLAB code(for computation), the magnetic field distribution was obtained. The calculation result fur a well-designed facility showed that the space in which the amplitude of the magnetic field lies within the 95 % of the magnetic field distribution was more than 60 % of each axis length.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.4
• /
• pp.1-6
• /
• 2005

The ability to predict radar return from flying vehicles becomes a critical technology issue in the development of stealth configurations. Toward developing a CEM code based on Maxwell's equations for analysis of RCS reduction schemes, an explicit upwind scheme suitable for multidisciplinary design is presented. The DFFT algorithm is utilized to convert the time-domain field values to the frequency-domain. A Green's function based on near field-to-far field transformation is also employed to calculate the bistatic RCS. To verify the numerical calculation the two-dimensional field around a perfectly conducting cylinder is considered. Finally results are obtained for the scattering electromagnetic field around an airfoil in order to illustrate the feasibility of applying CFD based methods to CEM.

• Structural Engineering and Mechanics
• /
• v.4 no.5
• /
• pp.553-568
• /
• 1996

Finite element stiffness matrix methods are presented for finding natural frequencies (or buckling loads) and modes of repetitive structures. The usual approximate finite element formulations are included, but more relevantly they also permit the use of 'exact finite elements', which account for distributed mass exactly by solving appropriate differential equations. A transcendental eigenvalue problem results, for which all the natural frequencies are found with certainty. The calculations are performed for a single repeating portion of a rotationally or linearly (in one, two or three directions) repetitive structure. The emphasis is on rotational periodicity, for which principal advantages include: any repeating portions can be connected together, not just adjacent ones; nodes can lie on, and members along, the axis of rotational periodicity; complex arithmetic is used for brevity of presentation and speed of computation; two types of rotationally periodic substructures can be used in a multi-level manner; multi-level non-periodic substructuring is permitted within the repeating portions of parent rotationally periodic structures or substructures and; all the substructuring is exact, i.e., the same answers are obtained whether or not substructuring is used. Numerical results are given for a rotationally periodic structure by using exact finite elements and two levels of rotationally periodic substructures. The solution time is about 500 times faster than if none of the rotational periodicity had been used. The solution time would have been about ten times faster still if the software used had included all the substructuring features presented.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.99-104
• /
• 2009

We propose a waveform inversion method for SH-wave data obtained in a shallow seismic refraction survey, to determine a 2D inhomogeneous S-wave profile of shallow soils. In this method, a 2.5D equation is used to simulate SH-wave propagation in 2D media. The equation is solved with the staggered grid finite-difference approximation to the 4th-order in space and 2nd-order in time, to compute a synthetic wave. The misfit, defined using differences between calculated and observed waveforms, is minimised with a hybrid heuristic search method. We parameterise a 2D subsurface structural model with blocks with different depth boundaries, and S-wave velocities in each block. Numerical experiments were conducted using synthetic SH-wave data with white noise for a model having a blind layer and irregular interfaces. We could reconstruct a structure including a blind layer with reasonable computation time from surface seismic refraction data.

• Korean Journal of Optics and Photonics
• /
• v.19 no.1
• /
• pp.1-8
• /
• 2008

Shack-Hartmann sensors are widely employed as a wavefront measuring device in various applications. Adaptive optics is one of the major applications. Since an adaptive optics system should be operated in real-time, high-speed wavefront sensing is essential. In high-speed operation, integration time of an image detector is very short. In this case, noises such as readout noise and photon noise greatly influence the accuracy of wavefront sensing. Therefore a fast and noise-insensitive centroid finding algorithm is required for the real-time wavefront sensing. In this paper, the multi-resolution correlation method is proposed. By employing multi-resolution images, this method greatly reduces the computation time when compared to the fast Fourier transform (FFT) correlation method. The verification is performed through the computational simulation. In this paper, the center of mass method, correlation method and multi-resolution correlation method are employed to compare the measurement accuracy of the centroid finding algorithms. The accuracy of a Shack-Hartmann wavefront sensor using the proposed algorithm is proved to be comparable to that of the conventional correlation method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.391-399
• /
• 2008

In this paper, we further generalize the work of Lin and Abel to the case of the first and the second order derivatives of energy release rates for two-dimensional, multiply cracked systems. The direct integral expressions are presented for the energy release rates and their first and second order derivatives. The salient feature of this numerical method is that the energy release rates and their first and second order derivatives can be computed in a single analysis. It is demonstrated through a set of examples that the proposed method gives expectedly decreasing, but acceptably accurate results for the energy release rates and their first and second order derivatives. The computed errors were approximately 0.5% for the energy release rates, $3\sim5%$ for their first order derivatives and $10\sim20%$ for their second order derivatives for the mesh densities used in the examples. Potential applications of the present method include a universal size effect model and a probabilistic fracture analysis of cracked structures.