• Nuclear Engineering and Technology
• /
• v.45 no.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.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.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
• /
• v.5 no.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.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.107-115
• /
• 1996

Displacement response is one of the important parameter to determine vibration characteristics of bridge structure. Reliable estimate of displacement response is obtained economically from integration of measured acceleration data in frequency domain. Proper sampling rate of frequency in discretization process of measured acceleration is proposed. Comparison of integrated and directly measured displacement response from laboratory experiment for a cantilever beam shows good agreement each other. Mode shape obtained from estimated displacement response also closely match with analytical result, thus the developed method is proved to be effective in practical use.

• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.125-131
• /
• 1990

A finite element technique for the time dependent large structural problems is presented. It is based on the error estimation for the bases of solution spaces. An a-posteriori energy norm of residual error serves as the error indicator. Mode shapes which are calculated by scaling the Ritz vectors are applied to discretize the continuous spatial domain. Finally, the performance of the proposed methods is demonstrated by solving simple examples.

• Ocean Systems Engineering
• /
• v.7 no.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.

• Bulletin of the Korean Mathematical Society
• /
• v.51 no.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.05a
• /
• 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.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.25 no.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
• /
• v.36 no.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.