• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.3
• /
• pp.279-285
• /
• 2011

We perform a numerical study to determine the time of onset of natural convection in a transient hot wire (THW) device for measuring the thermal conductivity of nanofluids. The samples used in this simulation are water-based $Al_2O_3$ nanofluids with volume fractions of 1%, 4%, and 10%, and the properties are calculated by theoretical models and experimental correlations. The THW apparatus using coated wire is modeled by the control-volume-based finite difference method, and the start of natural convection is determined by observing the temperature rise of the wire under a gravity field. The onset time is 11.5 s for water and 41.6 s for water-based $Al_2O_3$ nanofluids predicted by Maxwell thermal conductivity model with a 10% volume fraction. We confirm that the onset time of natural convection of nanofluids in the cylinder increases with the nanoparticle volume fraction. We suggest a correlation for predicting the onset time on the basis of the numerical results. Finally, it is shown that the measurement error due to natural convection is negligible if the measurement using the transient hot wire method is completed before the onset of natural convection in the base fluid.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.18 no.6
• /
• pp.569-579
• /
• 2016

In order to ensure the stability of tunnel portal slope, reinforcement method such as anchors, soil nails and rock bolts have been used in Korea. When selecting slope reinforcement methods in tunnel portal area such as reinforcement arrangement and length, trial and error method can be very time-consuming and it was also not easy to verify the selection of an optimum condition. In this study, using the FISH language embedded in the finite difference code FLAC3D program, the optimization technique was developed with the Differential Evolution Algorithm (DEA). After building a database on the soil nailing method in tunnel portal area, this system can be selected to an optimum arrangement plan based on the factor of safety through the FLAC3D analysis. Through the results of numerical analysis, it was confirmed that the number of analysis was decreased by about 8 times when DEA based optimization technique was used compared to the full combination (FC). In case of the design of slope reinforcement in tunnel portal area, if this built-system is used, it is expected that the selection of an optimum arrangement plan can be relatively easier.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.15 no.5
• /
• pp.473-479
• /
• 2004

In this paper, power density in near field is calculated about analytic object which has comparatively large volume in considering used wavelength such as cellular base station antenna. Panel sector antenna which is used widespreadly in domestic cellular wireless communication system is modeled and electromagnetic field distribution in reactive near field region is calculated by FDTD (Finite Difference Time Domain) method. After that, antenna gain in far field region is obtain by near to far transformation. Power spectral density in radiated near field is calculated in applying to gain-based model with antenna gain in far field. Finally, compliance distance is obtained in considering the result from radiated near field calculation and basic restrictions on occupational and general public exposure limits in ICNIRP guideline. In the center of main radiating position, the result from gain-based model is -14.55 ㏈m and the result from surface scanning method is -15.75 ㏈m. When the losses from cables and connectors used in measurement are considered, the results from gain-based model and surface scanning method are nearly coincident.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.2
• /
• pp.63-70
• /
• 2014

Design sensitivity analysis and topology design optimization for a piezoelectric crystal resonator are developed. The piezoelectric crystal resonator is deformed mechanically when subjected to electric charge on the electrodes, or vice versa. The Mindlin plate theory with higher-order interpolations along thickness direction is employed for analyzing the thickness-shear vibrations of the crystal resonator. Thin electrode plates are masked on the top and bottom layers of the crystal plate in order to enforce to vibrate it or detect electric signals. Although the electrode is very thin, its weight and shape could change the performance of the resonators. Thus, the design variables are the bulk material densities corresponding to the mass of masking electrode plates. An optimization problem is formulated to find the optimal topology of electrodes, maximizing the thickness-shear contribution of strain energy at the desired motion and restricting the allowable volume and area of masking plates. The necessary design gradients for the thickness-shear frequency(eigenvalue) and the corresponding mode shape(eigenvector) are computed very efficiently and accurately using the analytical design sensitivity analysis method using the eigenvector expansion concept. Through some demonstrative numerical examples, the design sensitivity analysis method is verified to be very efficient and accurate by comparing with the finite difference method. It is also observed that the optimal electrode design yields an improved mode shape and thickness-shear energy.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.441-449
• /
• 2015

In this paper we perform a linearized buckling analysis using the Kirchhoff plate theory and the von Karman nonlinear strain-displacement relation. Design sensitivity analysis(DSA) expressions for plane elasticity and buckling problems are derived with respect to Young's modulus and thickness. Using the design sensitivity, we can formulate the topology optimization method for minimizing the compliance and maximizing eigenvalues. We develop a topology optimization method applicable to plate buckling problems using the prestress for buckling analysis. Since the prestress is needed to assemble the stress matrix for buckling problem using the von Karman nonlinear strain, we introduced out-of-plane motion. The design variables are parameterized into normalized bulk material densities. The objective functions are the minimum compliance and the maximum eigenvalues and the constraint is the allowable volume. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with the finite difference ones and the topology optimization yields physically meaningful results.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.5
• /
• pp.612-620
• /
• 2019

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has results in considerable effort being directed toward developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have traditionally been obtained by the standard finite difference. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method(DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as excessive use of storage due to conditions of complex geometry and loading. In this study, the in-plane extensional vibration for asymmetric curved beams with linearly varying cross-section is analyzed using the DQM. Fundamental frequency parameters are calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results are compared with the result by other methods for cases in which they are available. According to the analysis of the solutions, the DQM, used only a limited number of grid points, gives results which agree very well with the exact ones.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.5
• /
• pp.129-137
• /
• 2008

An active application of concrete track is being expected far the future constructions of Korean railroad. In comparison with the existing ballasted tract, a concrete track is very susceptible for the settlement, since its rehabilitation requires much time and cost. When a concrete track is constructed on fine-grained subgrade soil, excessive cumulative plastic settlements due to repetitive train road may occur. In this case, the settlement of the concrete track may be effectively reduced by installing a small number of small-diameter piles beneath the track. This paper presents the effect of pile installation on the reduction of cumulative plastic settlement of concrete track. A method combining experiential equation and numerical method is proposed. Using an existing experiential equation and the estimated earth pressure distribution, the cumulative plastic strain was calculated. From the results, it is verified that the effects of the pile installation is significant to effectively reduce the cumulative plastic settlement of concrete track. The reduction effects of the cumulative plastic settlement according to the pile number and pile arrangement are presented.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.1
• /
• pp.53-66
• /
• 1980

The research work is concerned with the analytical and experimental studies on the heat transfer phenomenon around the underground concrete digester used for biogas production Systems. A mathematical and computational method was developed to estimate heat losses from underground cylindrical concrete digester used for biogas production systems. To test its feasibility and to evaluate thermal parameters of materials related, the method was applied to six physical model digesters. The cylindrical concrete digester was taken as a physical model, to which the model,atical model of heat balance can be applied. The mathematical model was transformed by means of finite element method and used to analyze temperature distribution with respect to several boundary conditions and design parameters. The design parameters of experimental digesters were selected as; three different sizes 40cm by 80cm, 80cm by 160cm and l00cm by 200cm in diameter and height; two different levels of insulation materials-plain concrete and vermiculite mixing in concrete; and two different types of installation-underground and half-exposed. In order to carry out a particular aim of this study, the liquid within the digester was substituted by water, and its temperature was controlled in five levels-35。 C, 30。 C, 25。 C, 20。C and 15。C; and the ambient air temperature and ground temperature were checked out of the system under natural winter climate conditions. The following results were drawn from the study. 1.The analytical method, by which the estimated values of temperature distribution around a cylindrical digester were obtained, was able to be generally accepted from the comparison of the estimated values with the measured. However, the difference between the estimated and measured temperature had a trend to be considerably increased when the ambient temperature was relatively low. This was mainly related variations of input parameters including the thermal conductivity of soil, applied to the numerical analysis. Consequently, the improvement of these input data for the simulated operation of the numerical analysis is expected as an approach to obtain better refined estimation. 2.The difference between estimated and measured heat losses was shown to have the similar trend to that of temperature distribution discussed above. 3.It was found that a map of isothermal lines drawn from the estimated temperature distribution was very useful for a general observation of the direction and rate of heat transfer within the boundary. From this analysis, it was interpreted that most of heat losses is passed through the triangular section bounded within 45 degrees toward the wall at the bottom edge of the digesten Therefore, any effective insulation should be considered within this region. 4.It was verified by experiment that heat loss per unit volume of liquid was reduced as the size of the digester became larger For instance, at the liquid temperature of 35˚ C, the heat loss per unit volume from the 0. 1m$^3$ digester was 1, 050 Kcal/hr m$^3$, while at for 1. 57m$^3$ digester was 150 Kcal/hr m$^3$. 5.In the light of insulation, the vermiculite concrete was consistently shown to be superior to the plain concrete. At the liquid temperature ranging from 15。 C to 350 C, the reduction of heat loss was ranged from 5% to 25% for the half-exposed digester, while from 10% to 28% for the fully underground digester. 6.In the comparison of heat loss between the half-exposed and underground digesters, the heat loss from the former was fr6m 1,6 to 2, 6 times as much as that from the latter. This leads to the evidence that the underground digester takes advantage of heat conservation during winter.

• Nuclear Engineering and Technology
• /
• v.22 no.2
• /
• pp.116-127
• /
• 1990

The fission gas release model used In the SPEAR-BETA fuel performance code was modified by use of effective thermal conductivity for cracked fuel and by laking Into account axial fission-gas mixing between the fuel-clad gap and the plenum. With use of this modified model the fission gas release was analyzed under various power ramping conditions of P$_{max}$ and $\Delta$.fP. Effective fuel thermal conductivity that accounts for the effect of fuel tracking was used in calculation of the fuel temperature distribution and the Internal gas pressure under power ramping conditions. Mixing and dilution effects due to axial gas flow were also considered in computing the width and the thermal conductivity of the gap. The effect of axial gas flow w3s solved by the Crank-Nicholson method. The finite difference method was used to save running time in the calculation. The present modified fission-gas release model was validated by comparing its predicted results with experimental data from various lamping tests In the literature and calculated results with use of the models used In the SPEAR-BETA and FEMAXI-IV codes. Results obtained with use of the present modified model showed better agreement with experimental data reported in the literature than those results with use of the latter codes. The fuel centerline temperature calculated with introduction of effective thermal conductivity for centerline temperature calculated with Introduction of effective thermal conductivity for cracked fuel was 200 higher fission gas release predicted with use of the modified model was nearly 6% larger on the average than that calculated by use of the unmodified model used in the SPEAR-BETA code.e SPEAR-BETA code.e.