• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.3
• /
• pp.139-148
• /
• 2000

In this paper, the interaction of oblique incident waves with a bottom-mounted and fluid-filled flexible membrane structure is investigated in the frame of linear hydro-elastic theory. The static shape of a membrane structure containing the fluid of a specific density is initially unknown and must be calculated before the hydrodynamic analysis. To solve hydrodynamic problem, the fluid domain is divided into the inner and outer region. The inner solution based on discrete membrane dynamic model and simple-source distribution over the entire fluid boundaries is matched to the outer solution ba~ed on an eigenfunction expansion method. The numerical results were compared to a series of Ohyama's experimental results. The measured reflection and tran¬smission coefficients reasonably follow the trend of predicted values. Using the computer program developed, the performance of a bottom-mounted and fluid-filled flexible membrane strocture is tested with various system parameters (membrane shape, internal pressure, density ratio) and wave characteristics (wave frequencies, incident wave angle). It is found that a bottom-mounted and fluid-filled flexible membrane structure can be an effel;tive wave barrier if properly designed.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.4
• /
• pp.121-128
• /
• 1998

In this study, when tension and bending stress act on plate simultaneously, stress intensity factor is analyzed at crack tip with using BEM(Boundary Element Method). In this analysis, stress intensity factors(S.I.F) are defined for variable ligament, aspect and stress ratio($\sigma$T/$\sigma$B). Consequently, predicted that crack grow to depth direction at low aspect and ligament ratio in tension stress and to surface direction in bending stress. Tension and bending stress act on plate same time, effect of tension stress in the first stage and effect of bending stress in the after stage was to observed. The outbreak of secondary crack in backside is under the control of stress amplitude and predict that the point of outbreak is mear backside.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1408-1416
• /
• 1990

A powerful and efficient method for finding approximate solutions to initial-boundary-value problems in the transient coupled thermoelasticity is formulated in time domain using the finite element technique with time-marching strategy. The final system equations can be derived by the Guritin's variational principle using the definition of convolution integral. But, the finite element formulation for the equations of motion is modified by differentiating in time. Numerical results to some test problems are compared with analytical and other sophisticated approximate solutions. Stable responces are observed in all the given examples irrespective of incremental time steps and mesh shapes. In addition, it is shown that good numerical results are obtained even in coarser mesh or larger time step comparing to other numerical methods.

• Tunnel and Underground Space
• /
• v.18 no.2
• /
• pp.134-148
• /
• 2008

The state of in situ stress is an important factor in considering the suitability of a site as a geological repository for nuclear waste. In this study, three-dimensional distinct numerical analysis was conducted to investigate the effect of deformation zones on the state of stress in the Oskarshamn area, which is one of two candidate sites in Sweden. A discontinuum numerical model was constructed by explicitly representing the numerous deformation zones identified from site investigation and far-field tectonic stress was applied in the constructed model. The numerical model successfully captured the variation of measured stress often observed in the rock mass containing large-scale fractures, which shows that numerical analysis can be an effective tool in improving the understanding of the state of stresses. Discrepancies between measured and modelled stress are attributed to the inconsistent quality of measured stress, uncertainty in geological geometry. and input data for fractures.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.3
• /
• pp.1-13
• /
• 1998

In the present paper, the hydrodynamic characteristics of three dimensional hydrofoils moving with a constant speed below the free surface using a higher-order boundary element method based on 9-node Lagrangian curvilinear elements are investigated. A bi-quadratic spline scheme is employed to improve the numerical results on the free surface. To validate the present scheme, the calculated results are compared with the analytic solutions for a submerged sphere and a spheroid showing a good agreement. For the validation of the hydrofoil study, the computed lift and drag of a hydrofoil having $NACA64_{1}A412$ section with aspect ratio(A.R.) of 4 are compared with the experimental data by Wadlin et al.[28]. The comparison covers a number of variations of angle of attack and submergence depth. Then, using an A.R. hydrofoil with NACA0012 section, the free surface on the lift and drag are investigated and these are compared with the previous results. The wave elevations and patterns created by the aforementioned submerged bodies are also investigated with Froude numbers and submergences.

• Journal of the Korea Concrete Institute
• /
• v.16 no.2 s.80
• /
• pp.155-162
• /
• 2004

An analytical method to predict the time dependent flexural behavior of composite girder is presented based on sectional analysis. The time dependent constitutive relation accounting for the early-age concrete properties including maturing of elastic modulus, creep and shrinkage is derived in an incremental format by the first order Taylor series expansion. The sectional analysis calculates the axial and curvature strains based on the force and moment equilibriums. The deflection curve of the girder approximated by the quadratic polynomial function is calculated by applying to the proper boundary conditions in the consecutive segments. Numerical applications are made for the 3-span double composite steel box girder which is a composite bridge girder filled with concrete at the bottom of the steel box in the negative moment region. The calculated results are compared with those by finite element analysis results. Close agreement is observed between the two approaches.

• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.227-235
• /
• 1998

A method of calculating the natural frequency corresponding to the first mode of vibration of beams and tower structures, with irregular cross-sections and with arbitrary boundary conditions was developed and reported by D. H. Kim in 1974. This method has been developed for two-dimensional problems including the laminated composite plates and was proved to be very effective for the plates with arbitrary boundary conditions and irregular sections. In this paper, the result of application of this method to the building slabs with passive and active control devices is presented. Finite difference method is used to obtain the deflection influence surfaces needed for this vibration analysis in this paper. The influence of the modulus of the foundation on the natural frequency is thoroughly studied.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.18 no.2
• /
• pp.101-108
• /
• 1982

In this paper the result of heat transfer analysis around the square heat source of interior solid by using the F. E. M is reported. Calculation for temperature distribution and each element was used by F. E. M. the solid is sub-divided into system of equal size triangular shape. These values of temperature distribution will valuable for design of jet engine and steam generator and the results gained are as follow; 1. Calculation by F. E. M is identified with the experiment. 2. Temperature distribution on the horizontal surface is $\theta$=0.698 in model 4 and the other hand $\theta$=0.401 in model 6 for X=16cm. In intermediates surface between heat source and bottom surface, the influence of L is more greater than that of height in the temperature difference. 3. Temperature distribution on the vertical surface for model 2 is resulted strong influence by K. In the case of Y=4cm is identified with $\theta$=0.0790 for K=7 and also $\theta$=0.0036 for K=0.3. In the difference of temperature distribution.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.33 no.3
• /
• pp.110-115
• /
• 2023

During the Hg₂Br₂ physical vapor transport process, with increasing the partial pressure of component B, P_B from 40 Torr to 200 Torr, a unicellular convective flow structures move from the crystal growth region to the center region in the vapor phase. The boundary layer flow is dominant for P_B = 40 Torr, and the core region flow is dominant for P_B = 200 Torr. The flow in the vapor phase shows a three-dimensional convective flow structure with a single cell (unicellular) for P_B = 40 Torr and 200 Torr, exhibits an asymmetrical flow with respect to the x, y central axis under the horizontally oriented configuration with an aspect ratio (length-to-width) of 3 and linear conducting walls. The critical temperature difference between the source and crystal region is about 30 K. The total molar flux of Hg₂Br₂ increases with the temperature difference until the total molar flux reaches the critical value. At the critical total molar flux, the total molar flux abruptly decreases.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.12
• /
• pp.889-898
• /
• 2022

In this paper, parallel computing method on the three-dimensional electromagnetic field is proposed. The present electromagnetic scattering analysis is conducted based on the time-harmonic vector wave equation and the finite element method. The edge-based element and 2^nd -order absorbing boundary condition are used. Parallelization of the elemental numerical integration and the matrix assemblage is accomplished by allocating the partitioned finite element subdomain for each processor. The graph partitioning library, METIS, is employed for the subdomain generation. The large sparse matrix computation is conducted by MUMPS, which is the parallel computing library based on the multi-frontal method. The accuracy of the present program is validated by the comparison against the Mie-series analytical solution and the results by ANSYS HFSS. In addition, the scalability is verified by measuring the speed-up in terms of the number of processors used. The present electromagnetic scattering analysis is performed for a perfect electric conductor sphere, isotropic/anisotropic dielectric sphere, and the missile configuration. The algorithm of the present program will be applied to the finite element and tearing method, aiming for the further extended parallel computing performance.