• 충청수학회지
• /
• 제23권2호
• /
• pp.271-277
• /
• 2010

We prove the global (in time) vorticity existence for the 2-D Euler equations of a perfect incompressible fluid in $B^0_{{\infty},1}({\mathbb{R}}^2){\cap}L^p({\mathbb{R}}^2)$ with 1 < p < 2. Moreover, we prove that the particle trajectory map X(x, t) satisfies the following estimate: for some positive constant C $${\parallel}X^{\pm1}(\cdot,\;t)-id(\cdot){\parallel}_{B^1_{\infty,1}}{\leq}Ce^{e^{Ct}}$$, where id represents the identity map on ${\mathbb{R}}^2$.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1595-1600
• /
• 2004

When a shock wave arrives at a duct, an impulsive wave is discharged from the duct exit and causes serious noise and vibration problems. In the current study, the characteristics of the impulsive wave discharged from a partial closed duct exit is numerically investigated using a CFD method. The Yee-Roe- Davis's total variation diminishing(TVD) scheme is used to solve the axisymmetric, unsteady, compressible Euler equations. With several partial closed duct exits, the Mach number of the incident shock wave $M_s$ and the distance L/D between the duct exit and a flat plate are varied in the range of $M_s$ = 1.01 ${\sim}$ 1.50 and L/D = 1.0 ${\sim}$ 4.0, respectively. The results obtained show that the magnitude of the impulsive wave impinging upon the flat plate strongly depends upon $M_s$, L/D and the partial closure of duct exit. The impulsive wave on the flat plate can be considerably alleviated by the partial closure of duct exit and, thus, the present method can be a passive control for the impulsive wave.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• pp.462-467
• /
• 2001

The present study addresses a computational work of the weak shock wave propagatings inside a silencer system of automobile exhaust pipe. Four different types of the silencer systems and the initial shock wave Mach number $M_s$ of $1.01\sim1.30$ are applied to investigate their effects on the noise reduction and the flow field in a silencer system. The results obtained from the present computational work are compared with the experimental results. The second order total variation diminishing (TVD) scheme is employed to solve the two dimensional, compressible, unsteady, Euler equations. The present computational results predict the experimental results with a quite good accuracy. Of the four silencer systems applied, the most desirable silencer system to reduce the peak pressure at the exit of the exhaust pipe is discussed.

• 한국소음진동공학회논문집
• /
• 제12권1호
• /
• pp.29-35
• /
• 2002

A simulation program for vehicle dynamic analysis was developed. The Cartesisn coordinate system was used for translational motion and the Euler angle system was used for rotational motion. A three dimensional multi-wheeled vehicle model and equations of motion were derived. Also static equilibrium analysis was added for initial vehicle condition setting. The program user can describe the exact characteristics of suspension spring force and damping force in the user subroutine. A wheel-ground contact model which represents geometrical effect was developed. Two cases of simulation for 16 D.O.F. vehicle model were conducted to validate the developed program by comparing the simulation results with the experimental data.

• 한국소음진동공학회논문집
• /
• 제12권12호
• /
• pp.943-949
• /
• 2002

The propagation of the impulse wave discharged from the Inclined exit of a pipe is investigated through shock tube experiment and numerical computations. The pressure histories and directivities of the impulse wave propagating outside from the exit of pipe with several different configurations are analyzed for the range of the incident shock wave Mach number between 1.1 and 1.4. In the shock tube experiments, the impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. Computations using the two-dimensional. unsteady, compressible, Euler equations are carried out to represent the experimented impulse waves. Computed Schlieren images predict the experimented impulse waves with a good accuracy. The results obtained show that for the radial direction the peak pressure of the impulse wave discharged depends upon the Inclined angle of the exit of the pipe. but for the axial direction it is almost constant regardless of the inclined angle of the pipe exit.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.151-154
• /
• 2002

The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to understand detailed flow physics associated with the twin impulse wave, compared with those in a single impulse wave. In the current study, the merging phenomena and propagation characteristics of the impulse waves are investigated using a shock tube experiment and by numerical computations. The Harten-Yee's total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The Mach number $M_{s}$, of incident shock wave is changed below 1.5 and the distance between two-parallel tubes, L/d, is changed from 1.2 to 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. The results obtained show that on the symmetric axis between two parallel tubes, the peak pressure produced by the twin-impulse waves and its location strongly depend upon the distance between two parallel tubes, L/d and the incident shock Mach number, $M_{s}$. The predicted Schlieren images represent the measured twin-impulse wave with a good accuracy.

• 한국전산구조공학회논문집
• /
• 제16권3호
• /
• pp.251-260
• /
• 2003

원형단면의 깊은 테이퍼봉과 보의 진동수와 모드형상을 결정하는 3차원 해석방법이 제시되었다. 수학적으로 1차원인 전통적인 봉과 보이론과는 달리, 본 연구에서는 3차원 동탄성방정식을 근간으로 하였다. 반경방향(r), 원주방향(θ), 축방향(z)으로의 변위성분인 u/sup r/, u/sub θ/, u/sub z/를 시간에 대해서는 정현적으로, θ에 대해서는 주기적으로, r과 z방향으로는 다수다항식의 형태로 표현하였다. 봉과 보의 위치(변형률)에너지와 운동에너지를 정식화하고, 고유치문제를 해결하기 위해 Ritz법을 사용하였으며, 진동수의 최소화과정을 통해 엄밀해의 상위경계치의 진동수를 구하였다. 이때 다항식의 차수를 증가시키면 진동수는 엄밀해에 수렴하게 된다. 봉과 보의 하위 5개의 진동수에 대해서 유효숫자 4자리까지의 수렴성 연구가 이루어졌다. 축방향으로 1차 직선적, 2차 및 3차 곡선으로 테이퍼된 9가지 형상의 봉과 보의 수치결과를 3차원 이론을 이용하여 최초로 계산하였다. 또한 선형 테이퍼 보의 예를 통해 3차원 Ritz법과 고전적인 1차원 Euler-Bernoulli 보이론과의 비교가 이루어졌다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2008년도 추계학술대회논문집
• /
• pp.71-81
• /
• 2008

This paper proposes a spectral element which can represent dynamic responses in high frequency domain such as Lamb waves on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by piezoelectric layer (PZT layer) bonded on a base plate. In the two layer beam model, a PZT layer is assumed to be rigidly bonded on a base beam. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with electro mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are formulated through equations of motions converted into frequency domain. A detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through comparison results with the conventional 2-D FEM and the previously developed spectral elements.

• Smart Structures and Systems
• /
• 제18권2호
• /
• pp.355-374
• /
• 2016

This paper presents and compares a one-dimensional (1D) bending theory for piezoelectric thin beam-type structures with resistive-inductive electrodes to ANSYS$^{(R)}$ three-dimensional (3D) finite element (FE) analysis. In particular, the lateral deflections and vibrations of slender piezoelectric beams are considered. The peculiarity of the piezoelectric beam model is the modeling of electrodes in such a manner that is does not fulfill the equipotential area condition. The case of ideal, perfectly conductive electrodes is a special case of our 1D model. Two-coupled partial differential equations are obtained for the lateral deflection and for the voltage distribution along the electrodes: the first one is an extended Bernoulli-Euler beam equation (second-order in time, forth order in space) and the second one the so-called Telegrapher's equation (second-order in time and space). Analytical results of our theory are validated by 3D electromechanically coupled FE simulations with ANSYS$^{(R)}$. A clamped-hinged beam is considered with various types of electrodes for the piezoelectric layers, which can be either resistive and/or inductive. A natural frequency analysis as well as quasi-static and dynamic simulations are performed. A good agreement between the extended beam theory and the FE results is found. Finally, the practical relevance of this type of electrodes is shown. It is found that the damping capability of properly tuned resistive or resistive-inductive electrodes exceeds the damping performance of beams, where the electrodes are simply linked to an optimized impedance.

• Nuclear Engineering and Technology
• /
• 제55권6호
• /
• pp.2230-2245
• /
• 2023

KANT (KAIST Advanced Nuclear Tachygraphy) is a PWR core simulator recently developed at Korea Advance Institute of Science and Technology, which solves three-dimensional steady-state and transient multigroup neutron diffusion equations under Cartesian geometries alongside the incorporation of thermal-hydraulics feedback effect for multi-physics calculation. It utilizes the standard Nodal Expansion Method (NEM) accelerated with various Coarse Mesh Finite Difference (CMFD) methods for neutronics calculation. For thermal-hydraulics (TH) calculation, a single-phase flow model and a one-dimensional cylindrical fuel rod heat conduction model are employed. The time-dependent neutronics and TH calculations are numerically solved through an implicit Euler scheme, where a detailed coupling strategy is presented in this paper alongside a description of nodal equivalence, macroscopic depletion, and pin power reconstruction. For validation of the steady, transient, and depletion calculation with pin power reconstruction capacity of KANT, solutions for various benchmark problems are presented. The IAEA 3-D PWR and 4-group KOEBERG problems were considered for the steady-state reactor benchmark problem. For transient calculations, LMW (Lagenbuch, Maurer and Werner) LWR and NEACRP 3-D PWR benchmarks were solved, where the latter problem includes thermal-hydraulics feedback. For macroscopic depletion with pin power reconstruction, a small PWR problem modified with KAIST benchmark model was solved. For validation of the multi-physics analysis capability of KANT concerning large-sized PWRs, the BEAVRS Cycle1 benchmark has been considered. It was found that KANT solutions are accurate and consistent compared to other published works.