• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.281-288
• /
• 2009

Using analytical method, intrinsic absorption and scattering attenuation was lately separated from total quality factor $(Q^{-1}_t)$ on the seismic data of Korean Peninsula. However, we should use numerical method rather than the analytical method to consider depth dependent structure of scattering. The direct simulation Monte Carlo (DSMC) method, as a kind of the numerical method, is good option due to its extended availability from 1 to 3-dimensional model; but there is few study to use it. In this paper, we introduced the analytical method and the DSMC method, and compared the results of the two analysis applied to the isotropic scattering model. While the scattering attenuation coefficients $(\eta_s)$ are identical, the intrinsic absorption coefficients $(\eta_i)$ for the analytical method are larger than those for DSMC method and have large errors. In addition, the $(Q^{-1}_t)$ by the previous studies show closer to DSMC method than analytical method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.10
• /
• pp.917-923
• /
• 2010

Direct numerical simulation (DNS) of a spatially developing turbulent boundary layer (TBL) with regularly arrayed cubical roughness elements was performed to investigate the effects of three-dimensional (3D) surface elements. The staggered cubes downstream were periodically arranged in the streamwise and spanwise directions with pitches of $p_x$/k=8 and $p_z$/k=2, where $p_x$ and $p_z$ are the streamwise and spanwise spacings of the cubes; the roughness height (k) was k=$1.5{\theta}_{in}$, where ${\theta}_{in}$ is the momentum thickness at the inlet. Spatially developing characteristics over the 3D cubical roughness were compared with the data obtained from the DNS over the two-dimensional (2D) rod roughened wall and smooth wall. Introduction of the cubical roughness on the TBL affected the turbulent Reynolds stresses not only in the roughness sublayer but also in the outer layer; and these effects are consistent with those observed over the 2D rough wall.

• Composites Research
• /
• v.19 no.1
• /
• pp.1-8
• /
• 2006

Prediction of damage caused by low-velocity impact in laminated composite plate is an important problem faced by designers using composites. Not only the inplane stresses but also the interlaminar normal and shear stresses playa role in estimating the damage caused. But it is well known that the conventional approach based on the homogenization has the limit in description of damage. The work reported here is an effort in getting better predictions of dynamic behavior and damage in composite plate using DNS approach. In the DNS model, we discretize the composite plates through separate modeling of fiber and matrix for the local microscopic analysis. In the view of microscopic mechanics with DNS model, interlaminar stress behaviors in the inside of composite materials are investigated and compared with the results of the homogenized model which has been used in the conventional approach to impact analysis. Also the multiscale model based on DNS concept is developed in order to enhance the effectiveness of impact analysis, and we present the results of multiscale analysis considering micro and macro structures simultaneously.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.5
• /
• pp.297-309
• /
• 2018

After an accident involving mooring link failures in an offloading buoy, verification of the fatigue safety in terms of the out-of-plane bending (OPB) and in-plane bending (IPB) moments has become a key engineering item in the design of various floating offshore units. The mooring links for an 8 MW floating offshore wind turbine were selected for this study. To identify the OPB stiffness (OPB moment versus interlink angle), a numerical simulation model, called the 3-link model, is usually composed of three successive chain links closest to the fairlead or chain hawse. This paper introduces two numerical simulation techniques for the 3-link analyses. The conventional and advanced approaches are both based on the prescribed rotation approach (PRA) and direct tension approach (DTA). Comparisons of the nominal stress distributions, OPB stiffnesses, hotspot stress curves, and stress concentration curves are presented. The multiple link analyses used to identify the tension angle versus interlink angle require the OPB stiffness data from the 3-link analyses. A convergence study was conducted to determine the minimum number of links for a multi-link analysis. It was proven that 10 links were sufficient for the multi-link analysis. The tension angle versus interlink angle relations are presented based on multi-link analyses with 10 links. It was found that the subsequent results varied significantly according to the 3-link analysis techniques.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.7
• /
• pp.56-66
• /
• 1999

Thermoacoustic oscillation is a significant problem in cylindrical-type combustors such as common internal combustion engines, industrial furnaces, gas turbine, etc. This kind of low frequency oscillation can lead to serious consequences such as destruction of the combustor and production of strong noise. The accurate numerical simulation of thermoacoustic phenomena is a complex and challenging problem, especially when considering the chemical reaction of mixtures. As with other simulations of aerodynamics and aeroacoustics, the direct computation of thermoacoustic phenomena requires that Navier-Stokes equations be solved using accurate numerical differentiation and time-marching schemes, with non-reflecting boundary conditions. The numerical approach used here aims at qualitative analysis and efficient prediction of those problems, not at the development of an accurate scheme. The numerical prediction developed in this work is shown to be reasonably matched with experimental result.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.3
• /
• pp.282-291
• /
• 2017

The main aim of this study is to investigate the effect of berm breakwater on wave run-up. A total of 200 numerical analysis tests have been carried out in this paper to investigate the effect of berm width, wave height, and wave period on the wave run-up, using an integrating technique of Computer-Aided Design (CAD) and Computational Fluid Dynamics (CFD). Direct application of Navier Stokes equations within the berm width has been used to provide a more reliable approach for studying the wave run-up over berm breakwaters. A well tested Reynolds-averaged Navier-Stokes (RANS) code with the Volume of Fluid (VOF) scheme was adopted for numerical computations. The computational results were compared with theoretical data to validate the model outputs. Numerical results showed that the simulation method can provide accurate estimations for wave run-up over berm breakwaters. It was found that the wave run-up may be decreased by increasing the berm width up to about 36 percent. Furthermore, the wave run-up may increase by increasing the wave height and wave period up to about 53 and 36 percent, respectively. These results may convince the engineers to use this model for design of berm breakwater in actual scale by calculating the Reynolds numbers.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.05a
• /
• pp.843-846
• /
• 2004

Cavity tone is generated due to the feedback between flow and acoustic wave. It is recognized that the period is determined by the time required for the flow convection in one direction, the time required for the acoustic propagation in the other direction and the time for phase shift depending on the flows and mode. Most of the phenomena have been investigated by experiments and a simple but fundamental theory. But the cause of the phase shift and the correctness of the theory have not been clearly explained so far. In this paper, the phenomena are calculated numerically to obtain detail information of flow and acoustic wave to explain the mechanism including the phase. High order high resolution scheme of optimized high order compact is used to resolve the small acoustic quantities and large flow quantities at the same time. The data are reduced using cross correlation function in space and time and cross spectral density function which has phase information. Abrupt change in pressure near corner in cavity is observed and is relate to phase variation. The time required for the feedback between the flow and acoustic wave is calculated after the numerical simulation f3r various modes. The periods based on the time calculated using the above method and direct observation from the acoustic waves generated and propagated in the numerical simulation are compared. It is found that no phase shift is required if we examine the time required carefully. Rossiter's formula for the cavity tone used for quick estimation needs to be modified far some modes.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.25 no.4
• /
• pp.1-11
• /
• 2021

To develop technologies for the stable operation of electric propulsion systems, the exhaust plume behavior of electric thrusters was studied using PIC-DSMC(particle-in-cell and direct simulation Monte Carlo). For the numerical analysis, the Simple Electron Fluid Model using Boltzmann relation was employed, and the charge and momentum exchanges due to atom-ion collisions were considered. The results of this study agreed with the plasma potentials measured experimentally. Near the thruster exit, active collisions among particles and charge exchanges created slow ions and fast atoms, which were expected to significantly affect the trajectory and velocity of the thruster exhaust plume.

• 한국연소학회:학술대회논문집
• /
• 2006.10a
• /
• pp.177-186
• /
• 2006

Flame surface area is a critical parameter determining turbulent flame speed. Three-dimensionaldirect numerical simulations (DNS) were conducted to figure out the evolution process of flame surface area. Fully compressible Navier-Stokes equations are solved to reproduce premixed flame embedded in isotropic decaying turbulent flow. The tangential straining and curvature of propagating surface affect development of flame area. In this study, four different turbulent intensity flows and three different Le number flames are investigated to force changes in straining and curvature effects. Consistent results are obtained for the probability density functions (PDF) of strain and curvature with previous researches. It is revealed that displacement speed, which is a speed of flame surface relative to unburnt flow, controls the balance between sink and source of flame surface area.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.1
• /
• pp.1-8
• /
• 2000

A modified low-Reynolds-number Reynolds stress model is developed for the calculation of drag-reducing turbulent flows induced by polymer injection. The results without polymer injection are compared with the results of direct numerical simulation to ensure the validity of the basic model. In case of drag reduction, profiles of mean velocity and Reynolds stress components, in two-dimensional channel flow, obtained with a proper value of viscosity ratio are presented and discussed. Computed mean velocity profile is in very good agreement with experimental data. And, the qualitative behavior of Reynolds stress components with the viscosity ratio is also reasonable.