• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.35-47
• /
• 1996

A surface panel method treating a boundary-value problem of the Dirichlet type is presented to design a three dimensional body with free surface corresponding to a prescribed pressure distribution. An integral equation is derived from Green's theorem, giving a relation between total potential of known strength and the unknown local flux. Upon discretization, a system of linear simultaneous equations is formed including free surface boundary condition and is solved for an assumed geometry. The pseudo local flux, present due to the incorrect positioning of the assumed geometry, plays a role f the geometry corrector, with which the new geometry is computed for the next iteration. Sample designs for submerged spheroids and Wigley hull and carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. For the calculation of the wave resistance, normal dipoles and Rankine sources are distributed on the body surface and Rankine sources on the free surface. The free surface boundary condition is linearized with respect to the oncoming flow. Four-points upwind finite difference scheme is used to compute the free surface boundary condition. A hyperboloidal panel is adopted to represent the hull surface, which can compensate the defects of the low-order panel method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but by experiment.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.137-145
• /
• 2014

We present a design sensitivity analysis(DSA) method for multiscale problems based on bridging scale decomposition. In this paper, we utilize a bridging scale method for the coupled system analysis. Since the analysis of full MD systems requires huge amount of computational costs, a coupled system of MD-level and continuum-level simulation is usually preferred. The information exchange between the MD and continuum levels is taken place at the MD-continuum boundary. In the bridging scale method, a generalized Langevin equation(GLE) is introduced for the reduced MD system and the GLE force using a time history kernel is applied at the boundary atoms in the MD system. Therefore, we can separately analyze the MD and continuum level simulations, which can accelerate the computing process. Once the simulation of coupled problems is successful, the need for the DSA is naturally arising for the optimization of macro-scale design, where the macro scale performance of the system is maximized considering the micro scale effects. The finite difference sensitivity is impractical for the gradient based optimization of large scale problems due to the restriction of computing costs but the analytical sensitivity for the coupled system is always accurate. In this study, we derive the analytical design sensitivity to verify the accuracy and applicability to the design optimization of the coupled system.

• Journal of the Korean Chemical Society
• /
• v.54 no.5
• /
• pp.594-602
• /
• 2010

Since the requirement of the high density integration and thin package technique of semiconductor have been increasing, the main package type of semiconductor will be a chip scale package (CSP). The changes of diffusion coefficient and moisture content ratio of epoxy resin systems according to the change of liquid type epoxy resin and fillers for CSP applications were investigated. The epoxy resins used in this study are RE-304S, RE310S, and HP-4032D, and Kayahard MCD as hardener and 2-methylimidazole as catalyst were used in these epoxy resin systems. The micro-sized and nano-sized spherical type fused silica as filler were used in order to study the moisture absorption properties of these epoxy molding compound (EMC) according to the change of filler size. The temperature of glass transition (Tg) of these EMC was measured using Dynamic Scanning Calorimeter (DSC), and the moisture absorption properties of these EMC according to the change of time were observed at $85^{\circ}C$ and 85% relative humidity condition using a thermo-hygrostat. The diffusion coefficients in these EMC were calculated in terms of modified Crank equation based on Ficks' law. An increase of diffusion coefficient and maximum moisture absorption ratio with Tg in these systems without filler can be observed, which are attributed to the increase of free volume with Tg. In the EMC with filler, the changes of Tg and maximum moisture absorption ratio with the filler content can be hardly observed, however, the diffusion coefficients of these systems with filler content show the outstanding changes according to the filler size. The diffusion via free volume is dominant in the EMC with micro-sized filler; however, the diffusion with the interaction of absorption according the increase of the filler surface area is dominant in the EMC with nano-sized filler.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.529-538
• /
• 2016

The simplified plate theory is presented for static and free vibration analysis of power-law(P) and sigmoid(S) Functionally Graded Materials(FGM) plates. This theory considers the parabolic distribution of the transverse shear stress, and satisfies the condition that requires the transverse shear stress to be zero on the upper and lower surfaces of the plate, without the shear correction factor. The simplified plate theory uses only four unknown variables and shares strong similarities with classical plate theory(CPT) in many aspects such as stress-resultant expressions, equation of motion and boundary conditions. The material properties of the plate are assumed to vary according to the power-law and sigmoid distributions of the volume fractions of the constituents. The Hamilton's principle is used to derive the equations of motion and Winkler-Pasternak elastic foundation model is employed. The results of static and dynamic responses for a simply supported FGM plate are calculated and a comparative analysis is carried out. The results of the comparative analysis with the solutions of references show relevant and accurate results for static and free vibration problems of FGM plates. Analytical solutions for the static and free vibration problems are presented so as to reveal the effects of the power law index, elastic foundation parameter, and side-to-thickness ratio.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.6
• /
• pp.64-72
• /
• 2017

The longitudinal dynamic instability which can occur in the fueling process of a space launch vehicle is called pogo. It is caused by coupling between the fuselage and propulsion system and they would be formed as a closed-loop system. so that the amplitude of the response may increase or decrease. In this paper, a mathematical model which is applicable to the systematic pogo analysis of a general launch vehicle is developed for an example of space shuttle. The formulations are composed of the linearized second-order differential equation for the propulsion system, and of the pressure, weight displacement, and generalized displacement. Those are important parameters for pogo analysis, are derived through eigenvalue analysis. By the formulation suggested in this paper, it is expected that mathematical modeling method of the pogo system can be obtained and systematic pogo stability analysis for any launch vehicle will be enabled.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.159-173
• /
• 1991

This paper describes a potential-based panel method formulated for the analysis of a super-cavitating two-dimensional hydrofoil. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring thats the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lilting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged. Characteristics of iteration and discretization of the present numerical method are much faster and more stable than the existing nonlinear theories. The theory shows good correlations with the existing theories and experimental results for the super-cavitating flow. In the region of small angles of attack, the present prediction shows and excellent comparison with the Geurst's linear theory. For the long cavity, the method recovers the trends of the Wu's nonlinear theory. In the intermediate regions of the short super-cavitation, the method compares very well with the experimental results of Parkin and also those of Silberman.

• Polymer(Korea)
• /
• v.36 no.2
• /
• pp.245-250
• /
• 2012

Branched polypropylenes (PP) with long chain branch were prepared by solid state reaction with three different branching agent of 0.3 wt% content. The chemical structures, non-isothermal crystallization behavior and complex viscosity of the branched PP were investigated by FTIR, DSC, optical microscope, and dynamic rheological measurement. The chemical structure of the branched PP was confirmed by the existence of =C-H stretching peak of the branching agent at 3100 $cm^{-1}$. There was no distinct change in melting temperature in case of PP-D-0-3 and PP-F-0-3, but PP-H-0-3 indicated a decrease in melting temperature. The decrease in melting temperature was interpreted by the fact that the degradation reaction of PP was more dominant than branched reaction, and confirmed by a decrease in complex viscosity. The non-isothermal crystallization behavior of the branched PP was analyzed using by Avrami equation. The Avrami exponent of PP was 3, and the values of the branched PP with DVB and FS were below 3. The activation energy of PP calculated by Kissinger method was 25 kJ/mol, and there were no big difference in activation energies of the branched PPs compared to PP.

• The Journal of Engineering Geology
• /
• v.18 no.4
• /
• pp.381-391
• /
• 2008

The eigenvalue technique is introduced to overcome the problem of truncation errors caused by temporal discretization of numerical analysis. The eigenvalue technique is different from simulation in that only the space is discretized. The spatially discretized equation is diagonized and the linear dynamic system is then decoupled. The time integration can be done independently and continuously for any nodal point at any time. The results of eigenvalue technique are compared with the exact solution and FEM numerical solution. The eigenvalue technique is more efficient than the FEM to the computation time and the computer storage in the same conditions. This technique is applied to the solute transport analysis in nonuniform flow fields around underground storage caverns. This method can be very useful for time consuming simulations. So, a sensitivity analysis is carried out by using this method to analyze the safety of caverns from nearly located contaminant sources. According to the simulations, the reaching time from source to the nearest cavern may takes 50 years with longitudinal dispersivity of 50 m and transversal dispersivity of 5 m, respectively.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.2
• /
• pp.70-80
• /
• 2000

Princeton Ocean Model (POM) is modified to construct a three-dimensional, semi-implicit hydro¬dynamic model with a wetting-and-drying scheme. The model employs semi-implicit treatment of the barotropic pressure gradient terms and the vertical mixing terms in the momentum equations, and the velocity divergence term in the vertically-integrated continuity equation. Such treatment removes the external mode and thus the mode splitting scheme in POM, allowing the semi-implicit model to use a larger time step. Applied to hypothetical systems, both the semi-implicit model and POM give nearly the same results. The semi-implicit model, however, runs approximately 4.4 times faster than POM showing its improved computational efficiency. Applied to a hypothetical system with intertidal flats, POM employing the mode splitting scheme produces noises at the intertidal flats, that propagate into the main channel resulting in unstable current velocities. Despite its larger time step, the semi-implicit model gives stable current velocities both at the intertidal flats and main channel. The semi-implicit model when applied to Kyeonggi Bay gives a good reproduction of the observed tides and tidal currents throughout the modeling domain, demonstrating its prototype applicability.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.1
• /
• pp.75-90
• /
• 1995

A mathematical model for the dynamic analysis of a riser system with the inclusion of internal flow and nonlinear effects due to large structural displacements is developed to investigate the effect of internal flow on marine riser dynamics. The riser system accounts fir the nonlinear boundary conditions and includes a steady flow inside the pipe which is modeled as an extensible or inextensible. tubular beam subject to nonlinear three dimensional hydrodynamic loads such as current or wave excitation. Galerkin's finite element approximation and time incremental operator are implemented to derive the matrix equation of equilibrium for the finite element system and the extensibility or inextensibility condition is used to reduce degree of freedom of the system and the required computational time in the case of a nonlinear model. The algorithm is implemented to develop computer programs used in several numerical applications. The investigations of the effect of infernal flow on riser vibration due to current or wave loading are performed according to the change of various parameters such as top tension, internal flow velocity, current velocity, wave period, and so on. It is found that the effect of internal flow can be controlled by the increase of top tension. However, careful consideration has to be given in the design point particularly for the long riser under the harmonic loading such as waves. And it is also found that the consideration of nonlinear effects due to large structural displacements increases the effect of internal flow on riser dynamics.