• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.4
• /
• pp.301-311
• /
• 2016

This paper describes the whole procedures to determine ice-induced global loads on the ship using measured full-scale data in accordance with the method proposed by the Canadian Hydraulics Centre of the National Research Council of Canada. Ship motions of 6 degrees of freedom (dof) are found by processing the commercial sensor signals named Motion Pak II under the assumption of rigid body motion. Linear accelerations as well as angular rates were measured by Motion Pak II data. To eliminate the noise of the measured data and the staircase signals due to the resolution of the sensor, a band pass filter that passes frequencies between 0.001 and 0.6 Hz and cubic spline interpolation resampling had been applied. 6 dof motions were computed by the integrating and/or differentiating the filtered signals. Added mass and damping force of the ship had been computed by the 3-dimensional panel method under the assumption of zero frequency. Once the coefficients of hydrodynamic and hydrostatic data as well as all the 6 dof motion data had been obtained, global ice loads can be computed by solving the fully coupled 6 dof equations of motion. Full-scale data were acquired while the ARAON rammed old ice floes in the high Arctic. Estimated ice impact forces for two representative events showed 7e15 MN when ship operated in heavy ice conditions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.9
• /
• pp.997-1002
• /
• 2012

With the increasing demand for lightweight materials to reduce fuel consumption, especially in the transportation industry, magnesium alloys are being widely studied. However, there are several limitations to the large-scale application of magnesium alloys in a structure because of their low formability and strong anisotropy. In order to take into account both the strong anisotropy and tension-compression asymmetry of AZ31 sheet alloy, the Cazacu-Plunkett-Barlat yield criterion (Cazacu, 2006) was adopted in material modeling. The variation of the anisotropic coefficients that describe the yield surface evolution of AZ31 is optimized using an interpolation function based on specific calibration results. It generates continuous yield surfaces, which makes it possible to describe different hardening rates in tension and compression as well as the tension-compression asymmetry of magnesium alloys. The performance of the CPB06 yield criterion for simulating an axial crushing test was tested and compared with that of the Hill (1948) yield criterion.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3732-3753
• /
• 2023

The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, developed based on the Multi-stage Response-function Transport (MRT) methodology, enables real-time simulation of nuclear systems such as reactor cores, spent nuclear fuel pools and casks, and sub-critical facilities. This paper presents the application of a novel fission matrix-based burnup methodology to the well-characterized JSI TRIGA Mark II research reactor. This methodology allows for calculation of nuclear fuel depletion by combination and interpolation of RAPID's burnup dependent fission matrix (FM) coefficients to take into account core changes due to burnup. The methodology is compared to experimentally validated Serpent-2 Monte Carlo depletion calculations. The results show that the burnup methodology for RAPID (bRAPID) implemented into RAPID is capable of accurately calculating the k_eff burnup changes of the reactor core as the average discrepancies throughout the whole burnup interval are 37 pcm. Furthermore, capability of accurately describing 3D fission source distribution changes with burnup is demonstrated by having less than 1% relative discrepancies compared to Serpent-2. Good agreement is observed for axially and pin-wise dependent fuel burnup and nuclear fuel nuclide composition as a function of burnup. It is demonstrated that bRAPID accurately describes burnup in areas with high gradients of neutron flux (e.g. vicinity of control rods). Observed discrepancies for some isotopes are explained by analyzing the neutron spectrum. This paper presents a powerful depletion calculation tool that is capable of characterization of spent nuclear fuel on the fly while the reactor is in operation.

• Journal of Environmental Science International
• /
• v.6 no.5
• /
• pp.541-550
• /
• 1997

A least-square regression analysis is applied for the estimation of velocity streamfunction field based on discretely sampled current meter data. The coefficients of a streamfuunction that is expanded in terms of trigonometric basis function are obtained by enforcing the horizontal non-divergence of two-dimensional flow field. This method avoids Interpolation and gives a root-mean-square (rms) residual of fit which Includes the divergent part and noisiness of oceanic data. The Implementation of the method Is done by employing a boundary-fitted, curvilinear orthogonal coordinate which facilitates the specification of boundary conditions. An application is successfully made to the Texas-Louisiana shelf using the 32 months current meter data (31 moorings) observed as a part of the Texas-Louisiana Shelf and Transport Processes Study (LATEX). The rms residual of the fitting is relatively small for the shelf, which indicates the field Is Ivell represented by the streamnfunction.

• Structural Engineering and Mechanics
• /
• v.68 no.4
• /
• pp.385-398
• /
• 2018

In the present paper, we offer a new flat shell finite element. It is the result of the combination of a membrane element and a bending element, both based on the strain-based formulation. It is known that $C^{\circ}$ plane membrane elements provide poor deflection and stress for problems where bending is dominant. In addition, they encounter continuity and compliance problems when they connect to C1 class plate elements. The reach of the present work is to surmount these problems when a membrane element is coupled with a thin plate element in order to construct a shell element. The membrane element used is a triangular element with four nodes, three nodes at the vertices of the triangle and the fourth one at its barycenter. Each node has three degrees of freedom, two translations and one rotation around the normal. The coefficients related to the degrees of freedom at the internal node are subsequently removed from the element stiffness matrix by using the static condensation technique. The interpolation functions of strain, displacements and stresses fields are developed from equilibrium conditions. The plate element used for the construction of the present shell element is a triangular four-node thin plate element based on Kirchhoff plate theory, the strain approach, the four fictitious node, the static condensation and the analytic integration. The shell element result of this combination is robust, competitive and efficient.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.577-582
• /
• 2007

A variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions in analysis domain arc generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Thus. the solution space can be represented in terms of the same functions to represent the geometry. The coefficients of basis functions or the control variables play the role of degrees-of-freedom. Furthermore, due to h-. p-, and k-refinement schemes, the high order geometric features can be described exactly and easily without tedious re-meshing process. The isogeometric sensitivity analysis method enables us to analyze arbitrarily shaped structures without re-meshing. Also, it provides a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling. To obtain precise shape sensitivity, the normal and curvature of boundary should be taken into account in the shape sensitivity expressions. However, in conventional finite element methods, the normal information is inaccurate and the curvature is generally missing due to the use of linear interpolation functions. A continuum-based adjoint sensitivity analysis method using the isogeometric approach is derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of boundary. In isogeometric analysis, however, the geometric properties arc already embedded in the B-spline shape functions and control points. The perturbation of control points in isogeometric analysis automatically results in shape changes. Using the conventional finite clement method, the inter-element continuity of the design space is not guaranteed so that the normal vector and curvature arc not accurate enough. On tile other hand, in isogeometric analysis, these values arc continuous over the whole design space so that accurate shape sensitivity can be obtained. Through numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Journal of the Korea Institute of Information Security & Cryptology
• /
• v.15 no.5
• /
• pp.93-103
• /
• 2005

In this paper, we propose a RST resistant robust logo embedding technique for multimedia copyright protection Geometric manipulations are challenging attacks in that they do not introduce the quality degradation very much but make the detection process very complex and difficult. Watermark embedding in the normalized image directly suffers from smoothing effect due to the interpolation during the image normalization. This can be avoided by estimating the transform parameters using an image normalization technique, instead of embedding in the normalized image. Conventional RST resistant schemes that use full frame transform suffer from the absence of effective perceptual masking methods. Thus, we adopt $8\times8$ block DCT and calculate masking using a spatio-frequency localization of the $8\times8$ block DCT coefficients. Simulation results show that the proposed algorithm is robust against various signal processing techniques, compression and geometrical manipulations.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.8 no.2
• /
• pp.126-136
• /
• 1997

When we apply the spectrally efficient quadrature amplitude modulation(QAM) to mobile communications, it is necessary to estimate and compensate the channel charac- teristics. In this paper, a new type fading estimation method for PSAM using sinc function is presented. Gaussian interpolation method has a drawback that the performance degrades rapidly if pilot symbol period increases even though pilot sysbol period is less than Nyquist sampling rate. The Wiener filter method does not degrade until pilot symbol period is equal to the Nyquist sampling rate. It is difficult for Wiener filter method to be applied to real system because autocorrelation function of channel gain, Doppler frequency and SNR(signal to noise ratio) must be known to optimize the filter coefficients. But proposed method has a similar performance to the Wiener filter method, and does not need to know the autocorrelation function of channel gain, the doppler frequency and SNR. Therefore the proposed method cna be applied to real system easily.

• Proceedings of the Korea Society of Information Technology Applications Conference
• /
• 2006.06a
• /
• pp.730-736
• /
• 2006

In this paper, we propose a rotation resistant robust logo embedding watermarking technique. Geometric manipulations make the detection process very complex and difficult. Watermark embedding ill the normalized image directly suffers from smoothing effect due to the interpolation during the image normalization. This can be avoided by estimating the transform parameters using image normalization technique, instead of embedding in the normalized image. Conventional rotation resistant schemes that use full frame transform. In this paper, we adopt $8{\times}8$ block DCT and calculate masking using a spatio-frequency localization of the $8{\times}8$ block DCT coefficients. Experimental results show that the proposed algorithm is robust against rotation process.

• Journal of Korean Society of Disaster and Security
• /
• v.11 no.2
• /
• pp.37-43
• /
• 2018

The storm surge is caused by an typhoons and it is not easy to predict the location, strength, route of the storm. Therefore, research using a scenario for storms occurrence has been conducted. In Korea, hazard maps for various scenarios were produced using the storm surge numerical simulation. Such a method has a disadvantage in that it is difficult to predict when other scenario occurs, and it is difficult to cope with in real time because the simulation time is long. In order to compensate for this, we developed a method to predict the storm surge damage by using research database. The risk grade prediction for the storm surge was performed predominantly in the study area of the East coast. In order to estimate the equation, COMSOL developed by COMSOL AB Corporation was utilized. Using some assumptions and limitations, the form of the basic equation was derived. the constants and coefficients in the equation were estimated by the trial and error method. Compared with the results, the spatial distribution of risk grade was similar except for the upper part of the map. In the case of the upper part of the map, it was shown that the resistance coefficient, k was calculated due to absence of elevation data. The SIND model is a method for real-time disaster prediction model and it is expected that it will be able to respond quickly to disasters caused by abnormal weather.