• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.234-239
• /
• 2001

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isometric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.161-164
• /
• 2004

A finite volume numerical model is developed for simulating non-equilibrium electroosmotic flow in micro- and nanochannels. The Guoy-Chapman model is adopted to compute the flow and electric potential. The Nernst-Planck equation is employed to trace unsteady transports of ionic species, i.e., time-dependent net charge density. A new set of boundary conditions based on surface charge density are designed rather than using the conventionally-employed zeta potential. A few issues for an efficient computation of electroosmotic flows are discussed. Representative computational examples are given to illustrate the robustness of the numerical model.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.47-50
• /
• 2006

In the aeroacoustic application of computational fluid dynamics, the physical phenomena like the crackle in the unsteady compressible jets should be based on very time-accurate numerical solution. The accuracy of the present numerical scheme is extended to the fifth order, using the WENO filter to the sixth-order central difference computation. However, the computational capacity is very restricted by the environment of computational power, so therefore the quadrilateral adaptive grids technique is introduced for this high-order accuracy scheme. The first problem is the multi-dimensional interpolation between fine and coarse grids. Some general benchmark problems are solved to show the effectiveness of this method.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.53-59
• /
• 2010

Preliminary design of a ground plate, a device installed close to the aircraft model for wind tunnel test to simulate the ground effect, was performed by a numerical simulation. A two-dimensional numerical study was performed initially to decide the optimal leading edge and flap configurations. Then, three-dimensional studies were conducted to decide the optimal flap deflection angle for pressure distribution reduction since the plate and the plate supporting system generate static pressure difference between the upper and lower flow regions. Three-dimensional simulation additionally studied the effect of the clearance between the plate and the wind tunnel side wall. For the efficiency of computation, half model was simulated and a symmetric boundary condition was applied on the center plane. Based on the preliminary design, a ground plate was designed, manufactured and tested at the Korea Aerospace Research Institute(KARI) wind tunnel. The measured pressure differences versus flap deflection angle agreed well with the predicted results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.05a
• /
• pp.457-464
• /
• 2001

The aeroacoustic phenomena in the simple rectangular open cavity are well published by many researchers. But the geometry shapes of aircraft landing gear wells, weapon bays, etc. are more complicate than that of the simple retangular cavity. They are more similar to the cavity having cover-plates at adges, or Helmholtz resonator. Therefore, the effects of cover-plates existing on edges of rectangular open cavity are numerically investigated in this paper. The compressible Navier-Stokes equations are solved for two-dimensional cavities with laminar boundary layers upstream. The high-order and high-resolution numerical schemes are used for the evaluation of spatial derivatives and the time integration. Physically correct numerical boundary conditions and buffer zone techniques are implemented to produce time-accurate solutions in the whole computation domain. The computational domain is large enough to directly resolve a portion of the radiated acoustic field. Results show that the cover-plates existing on edges of cavity reduce the noise convected from cavity, make the frequency of noise become higher, and change the directivity pattern. So these results can be used in the design of a low noise cavity.

• 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 for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.104-109
• /
• 2000

The modes of oscillation and radiated acoustic fields of compressible flows over open cavities are investigated computationally. The compressible Navier-Stokes equations are solved for two-dimensional cavities with laminar boundary layers upstream. The high-order and high-resolution numerical schemes are used for the evaluation of spatial derivatives and the time integration. Physically correct numerical boundary conditions are implemented to produce time-accurate solutions in the whole computation domain. The computational domain is large enough to directly resolve a portion of the radiated acoutic field. The results show a transition from a shear layer mode, for shorter cavities and lower Mach numbers, to a wake mode for longer cavities and higher Mach numbers. The shear layer mode is well characterized by Rossiter modes and these oscillations lead to intense upstream acoustic radiation dominated by a single frequency. The wake mode is characterized instead by a large-scale vortex shedding. Acoustic radiation is more intense, with multiple frequencies present.

• Journal of the Korean Society for Precision Engineering
• /
• v.7 no.1
• /
• pp.28-36
• /
• 1990

In this paper, the flow field with steady, axisymmetric flow characteristics in cylinder has been numerically investigated. The prediction of flow in cylinder was compared to experimental result. The turbulence closure being applied was K-${\epsilon}$ model and numerical scheme was Hybrid and Power-law scheme. The results of numerical computation showed some deviation from exper- imental data in the in cial region of cylinder, where streamline curvature is significant. However, the computational results agree qualitatively well with Mores's experimental results and the difference between Hybrid and Power-law schows similar results

• Communications for Statistical Applications and Methods
• /
• v.27 no.3
• /
• pp.377-383
• /
• 2020

The likelihood-based inference in a nonlinear generalized mixed model often requires computing moments of truncated multivariate normal random variables. Many methods have been proposed for the computation using a recurrence relation or the moment generating function; however, these methods rely on high dimensional numerical integrations. The numerical method is known to be inefficient for high dimensional integral in accuracy. Besides the accuracy, the methods demand too much computing time to use them in practical analyses. In this note, a moment calculation method is proposed under an assumption of a certain covariance structure that occurred mostly in generalized mixed models. The method needs only low dimensional numerical integrations.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.36-41
• /
• 1999

Shock wave propagating in the particle suspension has important applications. Examples are shock waves occurring in the solid rocket plume and detonation of dusty particles by shock waves. Experimental and numerical investigations on this subject have drawn much attention. More recently, Sivier et al. numerically simulated the experiment of Sommerfeld using the unstructured adaptive grid. They used the Eulerian-Eulerian approach based on the continuum assumption for both gas and particles. In the present paper, a new numerical method using the Lagrangian particle tracing technique and unstructured particle-adaptive grid for the polydisperse system is presented. It is explained why the existing numerical calculation has showed discrepancy with the experimental results by Sommerfeld.