• Journal of KSNVE
• /
• v.10 no.4
• /
• pp.610-614
• /
• 2000

Cellular Automata(CA)s are used as a simple mathematical model to investigate self-organization in statistical mechanics, which are originally introduced by von Neumann and S. Ulam at the end of the 1940s. CAs provide a framework for a large class of discrete models with homogeneous interactions, which are characterized by the following fundamental properties: 1) CAs are dynamical systems in which space and time are discrete. 2) The systems consist of a regular grid of cells. 3) Each cell is characterized by a state taken from a finite set of states and updated synchronously in discrete time steps according to a local, identical interaction rule. 4) The state of a cell is determined by the previous states of a surrounding neighborhood of cells. A cellular automaton has been attracted wide interest in modeling physical phenomena, which are described generally, partial differential equations such as diffusion and wave propagation. This paper describes one and two-dimensional analysis of wave propagation phenomena modeled by CA, where the local interaction rules were derived referring to the Lattice Gas Model reported by Chen et al., and also including finite difference scheme. Modeling processes by using CA are discussed and the simulation results of wave propagation with one wave source are compared with that by finite difference method.

• Aerospace Engineering and Technology
• /
• v.3 no.1
• /
• pp.1-8
• /
• 2004

In this study, the method of computing and extracting the generalized unsteady aerodynamic matrices using MSC/NASTRAN and MSC/NASTRAN DMAP ALTER has been suggested for the analysis and control of aeroelastic phenomena such as flutter and gust response analysis. In addition to that, the method of approximating the generalized unsteady aerodynamic matrices using minimum state approximation method has been proposed in order to cast the aeroelastic equations of motion in state space form for aeroelastic analysis and control application. Simplified aircraft wing box model has been used for the validation of the methods suggested in this study.

• Nuclear Engineering and Technology
• /
• v.43 no.3
• /
• pp.243-256
• /
• 2011

This paper considers the concept of analog computing based on a cellular neural network (CNN) paradigm to simulate nuclear reactor dynamics using a time-dependent second order form of the neutron transport equation. Instead of solving nuclear reactor dynamic equations numerically, which is time-consuming and suffers from such weaknesses as vulnerability to transient phenomena, accumulation of round-off errors and floating-point overflows, use is made of a new method based on a cellular neural network. The state-of-the-art shows the CNN as being an alternative solution to the conventional numerical computation method. Indeed CNN is an analog computing paradigm that performs ultra-fast calculations and provides accurate results. In this study use is made of the CNN model to simulate the space-time response of scalar flux distribution in steady state and transient conditions. The CNN model also is used to simulate step perturbation in the core. The accuracy and capability of the CNN model are examined in 2D Cartesian geometry for two fixed source problems, a mini-BWR assembly, and a TWIGL Seed/Blanket problem. We also use the CNN model concurrently for a typical small PWR assembly to simulate the effect of temperature feedback, poisons, and control rods on the scalar flux distribution.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1318-1325
• /
• 2017

Direct Numerical Simulation (DNS) serves as an irreplaceable tool to probe the complexities of multiphase flow and identify turbulent mechanisms that elude conventional experimental measurement techniques. The insights unlocked via its careful analysis can be used to guide the formulation and development of turbulence models used in multiphase computational fluid dynamics simulations of nuclear reactor applications. Here, we perform statistical analyses of DNS bubbly flow data generated by Bolotnov ($Re_{\tau}=400$) and LueTryggvason ($Re_{\tau}=150$), examining single-point statistics of mean and turbulent liquid properties, turbulent kinetic energy budgets, and two-point correlations in space and time. Deformability of the bubble interface is shown to have a dramatic impact on the liquid turbulent stresses and energy budgets. A reduction in temporal and spatial correlations for the streamwise turbulent stress (uu) is also observed at wall-normal distances of $y^+=15$, $y/{\delta}=0.5$, and $y/{\delta}=1.0$. These observations motivate the need for adaptation of length and time scales for bubble-induced turbulence models and serve as guidelines for future analyses of DNS bubbly flow data.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.53-59
• /
• 2011

In this study, steady and unsteady aerodynamic analyses of a huge rocket configuration have been conducted in a transonic flow region. The launch vehicle structural response are coupled with the transonic flow state transitions at the nose of the payload fairing. Before performing the coupled fluid-structure transonic aeroealstic simulations transonic aerodynamic characteristics are investigated for the pitching motions of the rocket at finite angle-of-attack. An unsteady CFD analysis method with a moving grid technique based on the Reynolds-averaged Navier-Stokes equations with the k-w SST transition turbulence model is applied to accurately predict the transonic loads of the rocket at pitching motion. It is shown that the fluctuating amplitude of the lateral aerodynamic loads imposed on the rocket due to the pitching motion can be significantly increased in the transonic flow region.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.8
• /
• pp.18-25
• /
• 2005

This paper describes a supersonic separation of air-launching rocket from the mother plane. Three dimensional Euler equations were numerically solved to analyze steady/unsteady state fluid flows. The results of simulation clearly demonstrate effect of shock-expansion wave interaction between the rocket and the mother plane. Moreover, important influential factors at separating stage of the rocket were extracted with a comprehensive analysis. Finally, from the consideration of supersonic-separation, a guideline to safety-separation is given to the design of supersonic air-launching rocket.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.1
• /
• pp.13-26
• /
• 2006

When the stagnation temperature of a perfect gas increases, the specific heat for constant pressure and ratio of the specefic heats do not remain constant any more and start to vary with this temperature. The gas remains perfect: its state equation remains always valid, with exception that it will be named by calorically imperfect gas. The aim of this research is to develop the relations of the necessary thermodynamics and geometrical ratios. and to study the supersonic flow at high temperature. lower than the threshold of dissociation. The results are found by the resolution of nonlinear algebraic equations and integration of complex analytical functions where the exact calculation is impossible. The dichotomy method is used to solve the nonlinear equation. and the Simpson algorithm for the numerical integration of the found integrals. A condensation of the nodes is used. Since. the functions to be integrated have a high gradient at the extremity of the interval of integration. The comparison is made with the calorifcally perfect gas to determine the error made by this last. The application is made for the air in a supersonic nozzle.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.1
• /
• pp.42-48
• /
• 2003

To get a periodic trim solution from Level II helicopter flight dynamic equations, DAE based PPT A (partial Periodic Trimming Algorithm) has been proposed. Iterative update of state variables from PPT A can cause a numerical instability in DAE solver which needs compatible initial conditions. By simply adjusting the order of DAE solver a periodic trim can be obtained with good accuracy. Application for CBM (Common Baseline Model) helicopter showed the same trim result as harmonic balance method and the effective elimination of unrealistic initial responses at the start of flight simulation.

• Advances in aircraft and spacecraft science
• /
• v.9 no.3
• /
• pp.263-278
• /
• 2022

The time-varying structural reliability of an aeroelastic launch vehicle subjected to stochastic parameters is investigated. The launch vehicle structure is under the combined action of several stochastic loads that include aerodynamics, thrust as well as internal combustion pressure. The launch vehicle's main body structural flexibility is modeled via the normal mode shapes of a free-free Euler beam, where the aerodynamic loadings on the vehicle are due to force on each incremental section of the vehicle. The rigid and elastic coupled nonlinear equations of motion are derived following the Lagrangian approach that results in a complete aeroelastic simulation for the prediction of the instantaneous launch vehicle rigid-body motion as well as the body elastic deformations. Reliability analysis has been performed based on two distinct limit state functions, defined as the maximum launch vehicle tip elastic deformation and also the maximum allowable stress occurring along the launch vehicle total length. In this fashion, the time-dependent reliability problem can be converted into an equivalent time-invariant reliability problem. Subsequently, the first-order reliability method, as well as the Monte Carlo simulation schemes, are employed to determine and verify the aeroelastic launch vehicle dynamic failure probability for a given flight time.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.11a
• /
• pp.519-522
• /
• 2006

Accurate modeling of a dynamic system from experimental data is the bases for the model updating or heath monitoring of the system. Modal analysis or modal test is a routine process to get the modal parameters of a dynamic system. The modal parameters include the natural frequencies, damping ratios and mode shapes. This paper presents a new method that can derive the equations of motion for a dynamic system from the modal parameters obtained by the modal analysis or modal test. The present method based on the relation between the eigenvalues and eigenvectors of the state space equation derives the mass, damping and stiffness matrices of the system. The modeling of a cantilevered beam from modal parameters is an example to prove the efficiency and accuracy of the present method. Using the lateral displacements only, not the rotations, gives limited information for the system. The numerical verification up to now gives reasonable results and the verification with the test data is scheduled.