• Wind and Structures
• /
• v.13 no.6
• /
• pp.557-580
• /
• 2010

A new inflow turbulence generation method and a combined dynamic SGS model recently developed by the authors were applied to evaluate the wind effects on 508 m high Taipei 101 Tower. Unlike the majority of the past studies on large eddy simulation (LES) of wind effects on tall buildings, the present numerical simulations were conducted for the full-scale tall building with Reynolds number greater than $10^8$. The inflow turbulent flow field was generated based on the new method called discretizing and synthesizing of random flow generation technique (DSRFG) with a prominent feature that the generated wind velocity fluctuations satisfy any target spectrum and target profiles of turbulence intensity and turbulence integral length scale. The new dynamic SGS model takes both advantages of one-equation SGS model and a dynamic production term without test-filtering operation, which is particular suitable to relative coarse grid situations and high Reynolds number flows. The results of comparative investigations with and without generation of inflow turbulence show that: (1) proper simulation of an inflow turbulent field is essential in accurate evaluation of dynamic wind loads on a tall building and the prescribed inflow turbulence characteristics can be adequately imposed on the inflow boundary by the DSRFG method; (2) the DSRFG can generate a large number of random vortex-like patterns in oncoming flow, leading to good agreements of both mean and dynamic forces with wind tunnel test results; (3) The dynamic mechanism of the adopted SGS model behaves adequately in the present LES and its integration with the DSRFG technique can provide satisfactory predictions of the wind effects on the super-tall building.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.154-160
• /
• 2007

This paper compares the characteristics of three different algorithms for three-dimensional (3D) magnetotelluric (MT) modeling. These methods are developed by Mackie et al. (1994), Sasaki (1999) and Nam et al. (2007). The first and second methods are based on the finite difference method (FDM), while the last one the finite-element method (FEM). MT responses, apparent resistivities and phases, for a COMMEMI 3D-2 model show a good agreement with integral equation solutions and only minor discrepancies are found over the anomalous bodies in the 3D model. The computation time of the two methods based on FDM is short and the static divergence correction contributes to speed up. The FEM modeling scheme is accurate but slow.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.4
• /
• pp.25-34
• /
• 1986

A two-dimensional linear method has been developed for the motion and the second-order steady force arising from the hydrodynamic coupling between waves and currents in the presence of a body of arbitrary shape. Interaction between the incident wave and current in the absence of the body lies in the realm beyond our interest. A Fredholm integral equation of the second kind is employed in association with the Haskind's potential for a steadily moving source of pulsating strength located in or below the free surface. The numerical calculations at the preliminary stage showed a significant fluctuation of the hydrodynamic forces on the surface-piercing body. The problem is approximately solved by using the asymptotic Green function for $U^2{\rightarrow}0$. The original Green function, however, is applied for the fully submerged body. Numerical calculations are made for a submerged and for a half-immersed circular cylinder and extensively for the mid-ship section of a Lewis-form. Some of the results are compared with other analytical results without any available experimental data. The current has strong influence on roll motion near resonance. When the current opposes the waves, the roll response are generally negligible in the low frequency region. The current has strong influence on roll motion near resonance. When the current opposes the wave, the roll response decreases. When the current and wave come from the same direction, the roll response increases significantly, as the current speed increases. The mean drift forces and moment on the submerged body are more affected by current than those on the semi-immersed circular cylinder or on the ship-like section in the encounter frequency domain.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• Polymer(Korea)
• /
• v.27 no.1
• /
• pp.33-39
• /
• 2003

In this study, the bisphenol-based DGEBA/GEBS blend systems were studied in cure kinetics, thermal stabilities, and fracture toughness of the casting specimen. The content of DGEBA/DCEBS was varied in 100 : 0, 90 : 10, 80 : 20, 70 : 30, and 60 : 40 wt%. The cure activation energies ($E_a$) of the blend systems were determined by Ozawa's equation. The thermal stabilities, including initial decomposed temperature (IDT), temperatures of maximum rate of degradation ($T_{max}$), and integral procedural decomposition temperature (IPDT) of the cured specimen were investigated by thermogravimetric analysis (TGA). For the mechanical interfacial properties of the specimens, the critical stress intensity factor ($K_{IC}$) test was performed and their fractured surfaces were examined by using a scanning electron microscope (SEM). As a result, $E_a$, IPDT, and $K_{IC}$ show maximum values in the 20 wt% DGEBS content compared with the neat DGEBA resins. This was probably due to the fact that the elevated networks were farmed by the introduction of sulfonyl groups of the DCEBS resin.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.4
• /
• pp.653-659
• /
• 1989

Elastic-plastic fracture toughness J$_{IC}$ can be used as an effective design criterion in elastic plastic fracture mechanics. In the J$_{IC}$ test methods approved by ASTM and JSME, there are discrepancies such as the definition of J$_{IC}$, the slope of the blunting line, curve fitting method and the measurement of crack extension etc. The objective of this paper is to evaluate the effect of these discrepancies on the determination of J$_{IC}$ values. Fracture toughness tests were performed on A516, SA508 and SCM415 steels, and test results were analyzed according to ASTM E 813-81, ASTM E 813-87 and JSME S 001-1981. Results showed significant differences depending on the analysis methods. Therefore, a conversion equation between two ASTM methods was proposed, and the conversion error was within acceptable range(less then 8.5%)en 8.5%)

• The Proceeding of the Korean Institute of Electromagnetic Engineering and Science
• /
• v.6 no.3
• /
• pp.3-14
• /
• 1995

The broadband input impedance, the input power and the radiation pattern of the monopole antenna attached to the handy phone operated at 800MHz are calculated by using the Finite Difference Time Domain(FDTD) Method. For the FDTD analysis of frequency characteristics of monopole antenna, the handy phone is modeled with the geometry that the monopole antenna is connected to a conducting box, and the modified FDTD algorithm[11] used the thin wire appproximation method and the Maxwell's integral equation from the original Yee algorithm is applied for the analysis of the wire structure. Also, by means of finding the current distribution directly from circumferencial magnetic filelds around the monopole antenna and the conducting box, the radiation pattern is calculated to observe the influence of the conducting box, and is compared with the results of the known mothod for the FDTD calculation of radiation pattern, For the experiments, the handy phone of which full length including antenna is .lambda. $\lambda$/2 is manufactured and we confirm that all computation results are agree well with the mea- sured values.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.18 no.1 s.116
• /
• pp.7-14
• /
• 2007

In this paper, the scan characteristics of phased array antenna consisted of rectangular open-ended waveguide with a triangular grid are investigated. An infinite array structure is analyzed by numerically solving the integral equation for the electric field over the waveguide aperture using waveguide mode function and Floquet mode function. Next, SEP(Scan Element Pattern) and SI(Scan Impedance) characteristics are simulated by CST's MWS(Microwave Studio) and Ansoft's HFSS(High Frequency Structure Simulator) for the finite and infinite array structures. Also, validity of these approaches is verified by comparing the calculated and simulated results with the measured ones for an $8{\times}8$ subarray. Within 10.5 % fractional bandwidth in the X-band, the fabricated subarray showed the flat gain characteristic in the scan range of ${\pm}45^{\circ}C$ in the E-plane(azimuth) and ${\pm}20^{\circ}C$ in the H-plane(elevation), and also showed the return loss characteristic of less than -10 dB.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.1
• /
• pp.44-51
• /
• 1996

The accurate calculation of run-up heights of long waves along the coastline is important in the view of engineering. In this paper the run-up heights of long waves are estimated by using the cnoidal wave theory which also covers both sinusoidal and solitary waves. However, the generation and the calculation of run-up heights of cnoidal waves are difficult both in laboratory and numerical experiments. In this study, the maximum run-up heights of cnoidal waves on steep slopes are computed by using the boundary integral equation model. It has been shown that the run-up heights of cnoidal waves are less than those of solitary waves, while they are larger than those of sinusoidal waves having the same wavelengths and heights. The variation of run-up heights of cnoidal waves is not a monotonic function of the wavelength. However, the run-up heights of cnoidal waves asymptotically approach that of a solitary wave as the wavelength approaches infinity. The calculated run-up heights agreed reasonably with experimental data.

• Journal of the Korean Society of Clothing and Textiles
• /
• v.26 no.5
• /
• pp.606-617
• /
• 2002

The purpose of this study is to develop a decision model that helps manufacturers and retailers determine the optimal timing of markdown in order to maximize their profit. An optimal timing decision model was developed based on three steps; conjoint measurement, scenario analysis and simulation. Data were collected from the sample of 149 out of 170 undergraduate and graduate students in Seoul in 1997. From the Jeans market, 8 brands; Levi's, lee, Guess, Calvin Klein, Pintos, Get used, MFG, and Basic, were selected as competitors for this study. In the conjoint measurement, respondents estimated the level of preference, from 1 to 100, for each item in which brand, price, style, and colors were used to explain product characteristics. Then, in order to reflect competitive situation in Jeans market, four types of scenarios were developed. In each scenario, simulations were applied to decide optimal timing of markdowns that leads to maximal profitability and sales volume. The profit was calculated based on the equation; Profit = Jean's market volume x market share of each brand - cost, where market volume was obtained by integral calculus for market utility function, and market share by logit value of part-worth from the conjoint analysis. For the purpose of the parsimony of the research, costs and the level of markdown were fixed to 30% of the regular price. In results, the optimal timing decision model identified 3 different types of brands. The brands that do not need to take markdown were Ievi's, MFG, and Basic Jeans characterized by the highest brand power and the highest price zone. The brands that needed to take early markdowns were Guess, Lee, Calvin Klein, and Get Used with the intermediate level of brand power and price. The brand that need late markdown was Pintos with the weakest brand power among the competitors and the lowest price. The optimal range of markdown remains for further research.