• Magazine of the Korean Society of Agricultural Engineers
• /
• v.40 no.5
• /
• pp.91-99
• /
• 1998

The widespread use of thin shell structures has created a need for a systematic method of analysis which can adequately account for arbitrary geometric form and boundary conditions as well as arbitrary general type of loading. Therefore, the stress and analysis of thin shell has been one of the more challenging areas of structural mechanics. A wide variety of numerical methods have been applied to the governing differential equations for spherical and cylindrical structures with a few results applicable to practice. The analysis of axisymmetric spherical shell is almost an every day occurrence in many industrial applications. A reliable and accurate finite element analysis procedure for such structures was needed. Dynamic loading of structures often causes excursions of stresses well into the inelastic range and the influence of geometry changes on the response is also significant in many cases. Therefore both material and geometric nonlinear effects should be considered. In general, the shell structures designed according to quasi-static analysis may fail under conditions of dynamic loading. For a more realistic prediction on the load carrying capacity of these shell, in addition to the dynamic effect, consideration should also include other factors such as nonlinearities in both material and geometry since these factors, in different manner, may also affect the magnitude of this capacity. The objective of this paper is to demonstrate the dynamic characteristics of spherical shell. For these purposes, the spherical shell subjected to uniformly distributed step load was analyzed for its large displacements elasto-viscoplastic static and dynamic response. Geometrically nonlinear behaviour is taken into account using a Total Lagrangian formulation and the material behaviour is assumed to elasto-viscoplastic model highly corresponding to the real behaviour of the material. The results for the dynamic characteristics of spherical shell in the cases under various conditions of base-radius/central height(a/H) and thickness/shell radius(t/R) were summarized as follows : The dynamic characteristics with a/H. 1) AS the a/H increases, the amplitude of displacement in creased. 2) The values of displacement dynamic magnification factor (DMF) were ranges from 2.9 to 6.3 in the crown of shell and the values of factor in the mid-point of shell were ranged from 1.8 to 2.6. 3) As the a/H increases, the values of DMF in the crown of shell is decreased rapidly but the values of DMF in mid-point shell is increased gradually. 4) The values of DMF of hoop-stresses were range from 3.6 to 6.8 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.3 to 2.6, and the values of DMF of stress were larger than that of displacement. The dynamic characteristics with t/R. 5) With the thickness of shell decreases, the amplitude of the displacement and the period increased. 6) The values of DMF of the displacement were ranged from 2.8 to 3.6 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.1 to 2.2.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.2
• /
• pp.79-93
• /
• 2015

Three-dimensional structures of large ocean rings in the Gulf Stream region are investigated using the HYbrid Coordinate Ocean Model (HYCOM). Numerically simulated flow structures around four selected cyclonic and anticyclonic rings are compared with two different horizontal resolutions: $1/12^{\circ}$ and $1/48^{\circ}$. The vertical distributions of Lagrangian Coherent Structures (LCSs) are analyzed using Finite Size Lyapunov Exponent (FSLE) and Okubo-Weiss parameters (OW). Curtain-shaped FSLE ridges are found in all four rings with extensions of surface ridges throughout the water columns, indicating that horizontal stirring is dominant over vertical motions. Near the high-resolution rings, many small-scale flow structures with size O(1~10) km are observed while these features are rarely found near the low-resolution rings. These small-scale structures affect the flow pattern around the rings as flow particles move more randomly in the high-resolution models. The dispersion rates are also affected by these small-scale structures as the relative horizontal dispersion coefficients are larger for the high-resolution models. The absolute vertical dispersion rates are, however, lower for the high-resolution models, because the particles tend to move along inclined eddy orbits when the resolution is low and this increases the magnitude of absolute vertical dispersion. Since relative vertical dispersion can reduce this effect from the orbital trajectories of particles, it gives a more reasonable magnitude range than absolute dispersion, and so is recommended in estimating vertical dispersion rates.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.10
• /
• pp.5-16
• /
• 2005

Geosynthetic-reinforced and pile-supported embankments have been increasingly used and researched around the world. The inclusion of one or multiple geosynthetic reinforcements over the pile is intended to enhance the efficiency of load transfer from soft ground to piles, to reduce total and differential settlement and increase global or local stability. In this paper, the reinforcement effectiveness and arching effect of the geogrid-reinforced and pile-supported embankments have been studied in terms of field model tests and numerical analysis with varying the space between piles and reinforcement. 2-dimensional numerical analysis has been conducted using the FLAC (Fast Lagrangian Analysis of Continua) program. And load transfer mechanisms between soil-piles-geogrid were investigated. The mechanisms of load transfer can be considered as a combination of embankment soil arching, tension geogrid, and stress concentration due to the stiffness difference between pile and soft ground. Based on the field model test and numerical analysis results, it was found that the geosynthetic reinforcement slightly interferes with soil arching, and helps reduce differential settlement of the soft ground. Also. at the D/b=3 (D: spacing of pile cap, b: diameter of pile), the total settlement is reduced by about $40\%$ compared to that without reinforcement. For $D/b{\ge}6$, the effectiveness of geogrid reinforcement in reducing settlement is negligible.

• The Korean Journal of Food And Nutrition
• /
• v.17 no.1
• /
• pp.41-46
• /
• 2004

This study was designed to investigate the effective ratio of healthy bread with jasmin tea powder. The weight, volume, and color difference of bread with various jasmin tea powder was measured and the picture of bread was taken and the panel test was taken. The weight of bread was increased and the volume was decreased significantly according to increase of jasmin tea powder. To avoid the decrease of volume, the active gluten was added, but there was not shown the increse of bread volume by active gluten. In the result of panel test and volume measurement, addition of 3% jasmin tea powder was appropriate. Therefore, using the jasmin tea powder, the diversity of bread by functionality, nutritional composition, taste, color, flavor of jasmin tea powder can be developed. This result will support the background the development of jasmin tea bread of good quality.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.20 no.4
• /
• pp.335-344
• /
• 2014

To understand the drift-diffusion of HAB(Harmful Algal Bloom) in this paper, we used three-dimensional hydrodynamic model POM(Pringceton Ocean Model) and Lagrangian particle track module. First, the results of residual flow that considered tide, wind, temperature, salinity, and TWC(Tsushima Warm Current) effect was tend to northeast in the coastal area and the flow in the offshore region showed results similar to TWC. To understand of HAB's movement, released each area that southern Kamak bay(Case 1), Mijo coast(Case 2), and southern Mireukdo coast(Case 3) assumption that red tide occurred. The areas where the HAB occurs frequently. As a result of HAB occurred in southern Kamak Bay(Case 1), mainly drifts to Narodo coast and Yeoja bay that located on the west side. Case 2 was mainly drifts to Yokjido coast and Saryangdo coast Especially, HAB occurred in Mireukdo coast(Case 3) relatively many particles drift to eastward as the influence of the TWC.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.9-16
• /
• 2014

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The "Hamilton-Jacobi(H-J)" equation and computationally robust numerical technique of "up-wind scheme" lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H-J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes are not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• The KIPS Transactions:PartA
• /
• v.14A no.4
• /
• pp.191-196
• /
• 2007

In the field of computer graphics, Navier-Stokes equations would be used for realistic simulations of smokes and currents. However, implementations derived from these equations are hard to achieve for real-time simulations, mainly due to its massive and complex calculations. Thus, there have been various attempts to approximate these equations for real-time simulation of smokes and others. When the advection terms of the equations are approximated by the Semi-Lagrange methods, the fluid density can be rapidly reduced and small-scale vorticity phenomena are easy to be missed, mainly due to the numerical losses over time. In this paper, we propose an improved numerical method to approximately calculate the advection terms, and thus eliminate these problems. To calculate the advection terms, our method starts to set critical regions around the target grid points. Then, among the grid points in a specific critical region, we search for a grid point which will be advected to the target grid point, and use the velocity of this grid point as its advection vector. This method would reduce the numerical losses in the calculation of densities and vorticity phenomena, and finally can implement more realistic smoke simulations. We also improve the overall efficiency of vector calculations and related operations through GPU-based implementation techniques, and thus finally achieve the real-time simulation.

• Nuclear Engineering and Technology
• /
• v.50 no.5
• /
• pp.665-672
• /
• 2018

A venturi scrubber is an important element of Filtered Containment Venting System (FCVS) for the removal of aerosols in contaminated air. The present work involves computational fluid dynamics (CFD) study of dust particle removal efficiency of a venturi scrubber operating in self-priming mode using ANSYS CFX. Titanium oxide ($TiO_2$) particles having sizes of 1 micron have been taken as dust particles. CFD methodology to simulate the venturi scrubber has been first developed. The cascade atomization and breakup (CAB) model has been used to predict deformation of water droplets, whereas the Eulerian-Lagrangian approach has been used to handle multiphase flow involving air, dust, and water. The developed methodology has been applied to simulate venturi scrubber geometry taken from the literature. Dust particle removal efficiency has been calculated for forced feed operation of venturi scrubber and found to be in good agreement with the results available in the literature. In the second part, venturi scrubber along with a tank has been modeled in CFX, and transient simulations have been performed to study self-priming phenomenon. Self-priming has been observed by plotting the velocity vector fields of water. Suction of water in the venturi scrubber occurred due to the difference between static pressure in the venturi scrubber and the hydrostatic pressure of water inside the tank. Dust particle removal efficiency has been calculated for inlet air velocities of 1 m/s and 3 m/s. It has been observed that removal efficiency is higher in case of higher inlet air velocity.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.199-206
• /
• 2003

Characteristics of the static shift are discussed by comparing the three-dimensional MT inversion with/without static shift parameterization. The galvanic distortion by small-scale shallow feature often leads severe distortion in inverted resistivity structures. The new inversion algorithm is applied to four numerical data sets contaminated by different amount of static shift. In real field data interpretations, we generally do not have any a-priori information about how much the data contains the static shift. In this study, we developed an algorithm for finding both Lagrangian multiplier for smoothness and the trade-off parameter for static shift, simultaneously in 3-D MT inversion. Applications of this inversion routine for the numerical data sets showed quite reasonable estimation of static shift parameters without any a-priori information. The inversion scheme is successfully applied to all the four data sets, even when the static shift does not obey the Gaussian distribution. Allowing the static shift parameters have non-zero degree of freedom to the inversion, we could get more accurate block resistivities as well as static shifts in the data. When inversion does not consider the static shift as inversion parameters (conventional MT inversion), the block resistivities on the surface are modified considerably to match possible static shift. The inhomogeneous blocks on the surface can generate the static shift at low frequencies. By those mechanisms, the conventional 3-D MT inversion can reconstruct the resistivity structures to some extent in the deeper parts even when moderate static shifts are in the data. As frequency increased, however, the galvanic distortion is not frequency independent any more, and thus the conventional inversion failed to fit the apparent resistivity and phase, especially when strong static shift is added. Even in such case, however, reasonable estimation of block resistivity as well as static shift parameters were obtained by 3-D MT inversion with static shift parameterization.