• Nuclear Engineering and Technology
• /
• 제45권6호
• /
• pp.707-720
• /
• 2013

We analyze piecewise homogenization with flux-weighted cross sections and preservation of averaged currents at the boundary of the homogenized domain. Introduction of a set of flux discontinuity ratios (FDR) that preserve reference interface currents leads to preservation of averaged region reaction rates and fluxes. We consider the class of numerical discretizations with one degree of freedom per volume and per surface and prove that when the homogenization and computing meshes are equal there is a unique solution for the FDRs which exactly preserve interface currents. For diffusion submeshing we introduce a Jacobian-Free Newton-Krylov method and for all cases considered obtain an 'exact' numerical solution (eight digits for the interface currents). The homogenization is completed by extending the familiar full assembly homogenization via flux discontinuity factors to the sides of regions laying on the boundary of the piecewise homogenized domain. Finally, for the familiar nodal discretization we numerically find that the FDRs obtained with no submesh (nearly at no cost) can be effectively used for whole-core diffusion calculations with submesh. This is not the case, however, for cell-centered finite differences.

• 대한조선학회논문집
• /
• 제49권1호
• /
• pp.60-67
• /
• 2012

This study considers the numerical analysis on parametric roll for container ships. As a method of numerical simulation, an impulse-response-function approach is applied in time domain. A systematic study is carried out for the parametric roll of two container ships, particularly observing the sensitivity of computational results to some parameters which can affect the analysis of parametric roll. The parameters to be considered are metacentric height (GM), simulation time window, and the discretization of wave spectrum. Based on the result of parametric roll simulation, numerical sensitivity and uncertainty in computational analysis are discussed.

• Coupled systems mechanics
• /
• 제5권3호
• /
• pp.235-254
• /
• 2016

A numerical model for fluid-structure interactions (abbr. FSI) is presented in the context of sloshing effects in movable, partially filled tanks to improve understanding of interactions between the fluid and the dynamics of a tank flexibly attached to a vehicle. The purpose of this model is to counteract the penalizing impact of the added mass effect on classical partitioned FSI coupling scheme: the proposed investigation is based on an added mass corrected version of the classical strongly coupled partitioned scheme presented in (Song et al. 2013). Results show that this corrected version systematically allows convergence to the coupled solution. In the rare cases where convergence is already obtained, the corrected version significantly reduces the number of iterations required. Finally, it is shown that the convergence limit imposed by added mass effect for the non-corrected coupling scheme, is directly dependent on the aspect ratio of the fluid domain and highly related to the precision order of the temporal discretization scheme.

• 전산구조공학
• /
• 제9권4호
• /
• pp.107-115
• /
• 1996

변위응답은 교량구조의 진동특성을 결정하는데 중요한 인자 중의 하나이다. 계측된 가속도 데이타를 진동수 영역에서 적분하여 변위응답을 경제적이고 합리적으로 구할 수 있다. 이를 위해 계측된 가속도 데이타를 이산화하기 위해 적절한 표본추출 진동수가 제시되었다. 캔틸레버 보를 이용한 실내시험에서 직접 계측된 변위와 적분된 변위는 서로 잘 일치하였다. 평가된 변위응답으로부터 구한 모우드 형상도 해석치와 근접하므로 개발된 방법은 실제로 효율적으로 사용될 수 있음을 입증하였다.

• 전산구조공학
• /
• 제3권3호
• /
• pp.125-131
• /
• 1990

본 논문에서는 시간 종속 하중을 받는 대형 전산구조응력 해석 문제를 위한 유한요소 모델링 기법이 소개된다. 유한요소모델링의 분할기준은 문제에 대한 해석 결과에 대한 오차선정에 근거를 둔다. 이 오차산정은 해석결과치에 의한 잔유 에너지의 크기를 유한요소별로 산정한다. 이의 시간존장 구조물에 대한 응용은, 구조연속체의 Ritz 고유진동 모드를 계산하고 이들 진동 모드 중에서 저주파에 상응하는 진동모드에 대해 잔유 에너지의 크기가 구조체 전체영역에서 평형을 유지하도록 유한요소 모델링을 수행한다. 마지막으로, 여기서 제안된 알고리즘이 몇 예제들을 통해서 검증된다.

• Ocean Systems Engineering
• /
• 제7권3호
• /
• pp.195-223
• /
• 2017

This paper presents 3D numerical simulations of a Free Standing Hybrid Riser under Vortex Induced Vibration, with prescribed motion on the top to replace the motion of the buoyancy can. The model is calculated using a fully implicit discretization scheme. The flow field around the riser is computed by solving the Navier-Stokes equations numerically. The fluid domain is discretized using the overset grid approach. Grid points in near-wall regions of riser are of high resolution, while far field flow is in relatively coarse grid. Fluid-structure interaction is accomplished by communication between fluid solver and riser motion solver. Simulation is based on previous experimental data. Two cases are studied with different current speeds, where the motion of the buoyancy can is approximated to a 'banana' shape. A fully three-dimensional CFD approach for VIV simulation for a top side moving Riser has been presented. This paper also presents a simulation of a riser connected to a platform under harmonic regular waves.

• 대한수학회보
• /
• 제51권3호
• /
• pp.717-737
• /
• 2014

We study the effect of numerical quadrature in space on semidiscrete and fully discrete piecewise linear finite element methods for parabolic interface problems. Optimal $L^2(L^2)$ and $L^2(H^1)$ error estimates are shown to hold for semidiscrete problem under suitable regularity of the true solution in whole domain. Further, fully discrete scheme based on backward Euler method has also analyzed and optimal $L^2(L^2)$ norm error estimate is established. The error estimates are obtained for fitted finite element discretization based on straight interface triangles.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.5-9
• /
• 2010

A numerical wind tunnel simulation is performed in order to predict wind loads acting on a building. The aim of the present study is to suggest a guideline for the numerical wind tunnel analysis, which could provide more detail wind load distributions compared to the wind code and expensive wind tunnel experiments. To validate the present numerical simulation, wind-induced loads on a 6 m cube model is predicted. Atmospheric boundary layer is used as a inlet boundary condition. Various effect of numerical methods are investigated such as size of computational domain, grid density, turbulence model and discretization scheme. The appropriate procedure for the numerical wind tunnel analysis is suggested through the present study.

• 한국소음진동공학회논문집
• /
• 제25권11호
• /
• pp.779-786
• /
• 2015

The Kansa method, which is used for various free vibration problems of arbitrarily shaped plates including membranes, discretizes the domain of a plate using only nodes without elements unlike FEM. The method requires a small amount of computation relative to FEM thanks to this discretization scheme but has limit in the accuracy of its solution. This paper reveals the reason of the limit and, to overcome the limit, proposes the practical method of calculating the singularity of a system matrix and extracting accurate natural frequencies. Case studies for a rectangular plate and an arbitrarily shaped plate validate the proposed method.

• Structural Engineering and Mechanics
• /
• 제36권3호
• /
• pp.343-361
• /
• 2010

The Discrete Element Method adopting particles for the domain discretization has recently been adopted in fracture studies of non-homogeneous continuous media such as concrete and rock. A model is proposed in which the reinforcement is modelled by 1D rigid-spring discrete elements. The rigid bars interact with the rigid circular particles that simulate the concrete through contact interfaces. The DEM enhanced model with reinforcement capabilities is evaluated using three point bending and four point bending tests on reinforced concrete beams without stirrups. Under three point bending, the model is shown to reproduce the expected final crack pattern, the crack propagation and the load displacement diagram. Under four point bending, the model is shown to match the experimental ultimate load, the size effect and the crack propagation and localization.