• Composites Research
• /
• v.18 no.5
• /
• pp.15-20
• /
• 2005

In this paper, damage induced acoustic emission in the composite plate in numerically simulated by using the three dimensional finite element method and explicit time integration. Acoustic source is modeled by equivalent volume source. To verify the proposed method, dynamic displacements due to the elastic wave are compared with the experiment when the fiber is broken in the single fiber embedded isotropic plate. For the laminated composite plates, the results are compared between homogenized model and DNS approach which models fibers and matrix separately. To capture high frequencies in the elastic wave, small time step size and a large number of meshes are required. The parallel computing technology is introduced to solve a large scale problem efficiently.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1634-1638
• /
• 2008

하천의 2차원 흐름 해석, 유사이동 해석, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있다. 유체에 대한 수치해석 기법으로 전통적으로 가장 많이 사용되고 있는 방법은 유한차분법이지만, 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 다양한 형태의 하천 해석에는 더욱 적합할 것이다. 본 연구에서는 비구조적 요소망을 advanced front method를 이용하여 생성해 보았다. Advanced front method는 해석하고자 하는 영역에 적절한 절점들을 생성한 후 삼각 요소망을 구성하는 grid based advanced front method와 절점들을 생성하지 않고 해석 영역에 삼각 요소를 바로 구성하는 element based advanced front method로 나눌 수 있다. Grid based advanced front method에서 해석 영역에 적절한 절점을 생성하는 방법으로는 일반적인 격자 구조의 절점 생성 방법을 적용하였으며 경계와의 거리가 가까운 절점은 생성되지 않으며, 삼각 요소를 구성할 때에는 두 개의 인접 절점을 비교하여 최적의 삼각 요소를 구성하게 된다. 단 두 개의 인접 절점만을 비교함으로서 비교적 빠른 시간 안에 최적의 삼각 요소망을 구성할 수 있다. 삼각 요소망을 생성한 후에는 Laplacian smoothing을 이용하여 삼각 요소망의 형질을 개선하였다. Element based advanced front method는 외부 경계에서부터 시작된 Front가 내부 영역으로 확대되어지며 각 Front에서 적절한 절점을 직접 생성하여 바로 삼각 요소를 구성하게 된다. Front에서 생성된 절점은 인접 절점들이 있는지 검색하여 인접 절점이 있다면 생성된 절점은 삭제되어지며 인접 절점이 삼각 요소를 위한 나머지 한 점으로 채택되어진다. Front는 외부 경계와 교차되어지지 않아야 하며 또한 연속된 Front를 효율적으로 관리하기 위해 list 자료 구조를 활용하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.337-337
• /
• 2015

홍수 저감, 생태계 복원, 위락 등 다양한 목적의 충족을 위해 강변에 저류지, 즉 다목적 유수지(detention basin)를 조성하는 사례가 나타나고 있다. 하천에서 홍수의 발생으로 수위가 어떤 기준보다 높아지면, 흐름의 일부를 돌려 저류지로 보냄으로써 본류의 부담을 덜 수 있다. 이때, 흐름의 분기를 위해 설치되는 하천구조물 중 하나가 측면 위어(side weir) 또는 횡월류 위어(side discharge/overflow weir)이다. 하천의 계획과 설계에서 위어가 적용될 때, 위어에 대한 수위-유량 관계, 즉 그 형식과 제원에 적합한 유량계수(discharge coefficient)의 결정이 관건이 된다. 일반적인 위어와 달리 흐름 양상이 복잡한 측면 위어의 경우, 이론과 실제의 괴리가 아직까지 해소되지 않아 실물 또는 3차원 수치 모형을 이용한 시험으로 수위-유량 관계를 수립할 필요가 있다. 이렇게 결정된 수위-유량 관계는 1차원 또는 수심적분 2차원 모형의 내부 또는 외부 경계로 사용되며, 본류의 수위 증감에 따른 측면 위어의 횡월류량을 통해 저류지의 홍수 조절 능력을 평가할 수 있다. 이 연구에서는, 측면 위어의 수위-유량 관계가 알려지지 않더라도, 저류지에 의한 홍수 조절 효과를 평가할 수 있는 2차원 수치모의에 대해 검토하였다. 수치해법으로서 2차원 천수방정식에 대해 유한체적법을 적용하고, 흐름률(flux)의 정확한 계산을 위해 근사 Riemann 해법을 도입하였다. 먼저, 측면 위어가 없는 실험 조건에 대해 수로 내 한 측선에서 측정된 수위와 유량을 모의 결과와 비교하여 모형을 검증하였다. 이때, 경계조건으로 상류 끝에 측정 유량을, 하류 끝에 측정 수위를 부여하였으며, Manning의 조도계수를 0.014로 설정하였다. 또한, 측면 위어가 설치된 수로에 대해 계산 영역을 340개의 삼각형 격자로 분할하고 측면 위어가 없는 경우와 동일한 조건을 두어 모의하였다. 측면 위어의 하류에 위치한 측선에서 측정치에 대한 평균 제곱근(root mean square) 오차가 수위에 대해 1.9 mm, 유량에 대해 $2.2{\ell}/s$로서 그림과 같이 모의 결과는 실험의 그것과 잘 일치하였다. 이로써, 측면 위어에 대한 수위-유량 관계의 수립을 위한 실물 모형 시험 없이 수심적분 2차원 수치모의를 통해 저류지의 홍수 조절 효과를 평가할 수 있음이 확인되었다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.4
• /
• pp.297-303
• /
• 2014

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that $C^1$-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5B
• /
• pp.429-439
• /
• 2009

When Catastrophic extreme flood occurs due to dam break, the response time for flood warning is much shorter than for natural floods. Numerical models can be powerful tools to predict behaviors in flood wave propagation and to provide the information about the flooded area, wave front arrival time and water depth and so on. But flood wave propagation due to dam break can be a process of difficult mathematical characterization since the flood wave includes discontinuous flow and dry bed propagation. Nevertheless, a lot of numerical models using finite volume method have been recently developed to simulate flood inundation due to dam break. As Finite volume methods are based on the integral form of the conservation equations, finite volume model can easily capture discontinuous flows and shock wave. In this study the numerical model using Riemann approximate solvers and finite volume method applied to the conservative form for two-dimensional shallow water equation was developed. The MUSCL scheme with surface gradient method for reconstruction of conservation variables in continuity and momentum equations is used in the predictor-corrector procedure and the scheme is second order accurate both in space and time. The developed finite volume model is applied to 2D partial dam break flows and dam break flows with triangular bump and validated by comparing numerical solution with laboratory measurements data and other researcher's data.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.2 no.1
• /
• pp.30-37
• /
• 1999

Transmissivity and the signal loss in soil are computed. An electric field expression for the inhomogeneous complex media modelled by two layers is shown as an integral form. Volume scattering occurs in inhomogeneous media, and iterative Born approximation is used to analyze this scattering effect. The degree of randomness is controlled by specifying the variance and correlation length. Expression for the transmissivity and the signal loss is presented as the parameter of soil moisture contents, soil particle radius, temperature and frequency. The analysis shows that big deviation in signal loss depends on the temperature variation remarkably and the physical reason of unusual level is explained.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.5
• /
• pp.405-412
• /
• 2012

A flow control on airfoil installed a vortex cell for high efficiency wind turbine blade in stationary and dynamic stall conditions have been numerically investigated by solving the compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with Roe's flux-difference splitting and an implicit time-integration method coupled with dual time step sub-iteration. The computed result for the airfoil in the stationary showed that lift-drag ratio increases due to low pressure by the vortex cell. The oscillating airfoil with the vortex cell showed that the magnitude of hysteresis loop is reduced due to the enhanced vortex in the cell.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.18 no.7
• /
• pp.819-827
• /
• 2007

A numerical algorithm for estimating precise backscattering coefficients of rice fields is proposed and its accuracy is verified in this paper. After a bunch of rice plants above water surface is modeled with a bunch of randomly oriented lossy dielectric bodies above an impedance surface and the equivalent volume currents of the lossy dielectrics are computed using the moment method. Then, the scattered fields of a rice field with many bunches are computed with a Monte Carlo method, and consequently the backscattering coefficient of the rice field is computed for various incidence angles and polarizations. Finally, the backscattering coefficient of a rice field is measured at 1.85 GHz using an R-band scatterometer system, and these experimental data are used to verify the numerical algorithm proposed in this paper. It is found that the numerical computation results agree well with the measurement data.

• The Journal of the Acoustical Society of Korea
• /
• v.42 no.3
• /
• pp.250-261
• /
• 2023

In this study, the noise contributions by the duct, stator and rotor, which are the propulsor components, are evaluated to identify the flow noise source in cavitation and non-cavitation conditions on pumpjet propulsion and the noise levels in both conditions are compared. The unsteady incompressible Reynolds averaged Navier-Stokes (RANS) equation based on the homogeneous mixture assumption is applied on the suboff submarine hull and pumpjet propeller in the cavitation tunnel, and the Volume of Fluid (VOF) method and Schnerr-Sauer cavitation model are used to describe the two-phase flow. Based on the flow simulation results, the acoustic analogy formulated by Ffowcs Williams and Hawkings (FW-H) equation is applied to predict the underwater radiated noise. The noise contributions are evaluated by using the three types of impermeable integral surface on the duct, stator and rotor, and the two types of permeable integral surface surrounding the propulsor. As a result of noise prediction, the contribution by the stator is insignificant, but it affects the generation of flow noise source due to flow separation in the duct and rotor, and the noise is predominantly radiated into the upward and right where the flow separations are. Also, the noise is radiated into the thrust direction due to pressure fluctuation between suction and pressure sides on the rotor blades, and the it can be seen that the cavitation effect into the noise can be considered through the permeable integral surface.