• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1248-1253
• /
• 2003

Sensitivity information has been used for linearization of nonlinear functions in optimization. Basically, sensitivity is a derivative of a function with respect to a design variable. Design sensitivity is repeatedly calculated in optimization. Since sensitivity calculation is extremely expensive, there are studies to directly use the sensitivity in the design process. When a small design change is required, an engineer makes design changes by considering the sensitivity information. Generally, the current process is performed one-by-one for design variables. Methods to exploit the sensitivity information are developed. When a designer wants to change multiple variables with some relationship, the directional derivative can be utilized. In this case, the first derivative can be calculated. Only small design changes can be made from the first derivatives. Orthogonal arrays can be used for moderate changes of multiple variables. Analysis of Variance is carried out to find out the regional influence of variables. A flow is developed for efficient use of the methods. The sensitivity information is calculated by finite difference method. Various examples are solved to evaluate the proposed algorithm.

• The Journal of Engineering Geology
• /
• v.18 no.4
• /
• pp.493-499
• /
• 2008

In order to evaluate the lateral behavior characteristics of single drilled shaft embedded in granite gneiss, a lateral load test was conducted in field. Horizontal displacement according to lateral load were measured along with the depth by an inclinometer installed in the shaft. In this study, We have evaluated horizontal displacement characteristics comparing the measurement values with calculating results by theoretical formula. Based on the comparison, the Chang's method was similar with the measurement values even though it was slightly underestimated. However, the finite analysis method and p-y method was overestimated, especially on the upper part of the ground.

• Nuclear Engineering and Technology
• /
• v.18 no.2
• /
• pp.71-77
• /
• 1986

A one dimensional neutron kinetics program, BIK which is applicable to the safety analyses of PWR's is developed to analyze the reactor core in axial dimension. The BIK employs the finite difference technique in space and $\theta$-time integration method in time. Detailed models for the Doppler and moderator feedbacks and control rod motion are included. The benchmark of the nuclear model is carried out through the ANL benchmark problem and the time dependent nuclear power change in the rod ejection accident of KNU1 is calculated by BIK code. The results indicate that the BIK can predict the neutron dynamics with fair accuracy within the limits of one dimensional analysis and it is useful for the safety analyses of PWR's.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.17 no.4
• /
• pp.227-236
• /
• 1997

Quantitative analysis and imaging of elastic wave propagation are very important for the materials evaluation as well as flaw detection. The elastic wave propagation in an anisotropic media is more complex, and analysis and imaging become essential for flaw detection and materials evaluation. In the anisotropic media, the wave velocity is dependent on the propagation direction. In addition, the direction of group velocity is different from that of phase velocity, the direction of energy flow is not same as the propagation direction of wavefront (beam skewing effect). Especially, this effect becomes critical for the large anisotropic media such as fiber composite materials, and the results using elastic waves for those materials have to be analyzed considering the wave propagation mechanism. Since the analytical approach for the wave propagation in the anisotropic materials is limited, the numerical analysis such as finite difference method (FDM) have been used for these case. Therefore, 2-dimensional FDM program for the elastic wave propagation is developed, and wave propagation in anisotropic media are simulated.

• Journal of Korean Society of Steel Construction
• /
• v.14 no.4
• /
• pp.489-497
• /
• 2002

Nonlinear behavior and large deformation cannot be analyzed using techniques based on linear theory. Nonetheless, they are emerging as gradually huge and complex structures. In addition, the optimum design of structure is necessary in the development of high-performance computation and numerical methods. as well as stricter design-criterion. Therefore, the structural problems in engineering that are limited to the linear region must be extended to the nonlinear region. Likewise, structural behavior must be accurately analyzed. In turn, this requires considering the expected problems beforehand. Only then can an efficient, economical, and optimized structure be designed. This paper presents the solution of the geometrical nonlinear problem of anisotropic cylindrical shell. The characteristics of the geometrical nonlinear behavior of anisotropic circular cylindrical shells may vary according to several causes. e.g., change of fibers, curvature in the circumferential direction, subtended angle, aspect, etc. Parametric studies were conducted to determine the effect of factors on the large deflection behavior of laminated shells, with interesting observations.

• Journal of the Korean Association of Geographic Information Studies
• /
• v.2 no.2
• /
• pp.45-58
• /
• 1999

Digital map is utilized for an effective display of the distributions of the hydrogeological variables such as water table height(hydraulic head) and log-transmissivity(lnT) in north Pohang, KyungBug. Specifically the geostatistical method kriging is used to construct the distributions in an unconfined aquifer from a finite set of measured data. The experimental variograms for both the head and lnT suggest spherical models with nugget of 0 and range of 6km. The kriged results by using these variograms show that the head decreases primarily from the west to the east with a large peak in the north-western part and lnT is at the maximum level in the central part with outwardly decreasing trend. The constructed delineation is also used to calculate the flow field in the region. Finite differences with second order consistency are used to calculate the fluxes in the east(x) and north(y) across a vertical cross-section of unit width and height equal to the thickness of the wet zone in the aquifer. It is demonstrated that the flow is dominantly in the east with diverging trend on the eastern hillside of the water table peak. A few convergent spots also appear.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.3
• /
• pp.23-27
• /
• 2008

The new conjunctive surface-subsurface flow model at a large scale was developed by using a 1-D Diffusion Wave (DW) model for surface flow interacting with the 3-D Volume Averaged Soil-moisture Transport (VAST) model for subsurface flow for the comprehensive terrestrial water and energy predictions in Land Surface Models (LSMs). A selection of numerical implementation schemes is employed for each flow component. The 3-D VAST model is implemented using a time splitting scheme applying an explicit method for lateral flow after a fully implicit method for vertical flow. The 1-D DW model is then solved by MacCormack finite difference scheme. This new conjunctive flow model is substituted for the existing 1-D hydrologic scheme in Common Land Model (CLM), one of the state-of-the-art LSMs. The new conjunctive flow model coupled to CLM is tested for a study domain around the Ohio Valley. The simulation results show that the interaction between surface flow and subsurface flow associated with the flow routing scheme matches the runoff prediction with the observations more closely in the new coupled CLM simulations. This improved terrestrial hydrologic module will be coupled to the Climate extension of the next-generation Weather Research and Forecasting (CWRF) model for advanced regional, continental, and global hydroclimatological studies and the prevention of disasters caused by climate changes.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.111-123
• /
• 2007

Unlike the conventional reverse time migration method which is implemented by simply extrapolating wavefield in reverse time, this paper presents a derivation of another reverse time migration operator as the exact adjoint of the presumed forward wavefield extrapolation operator. The adjoint operator is obtained by formulating the forward time extrapolation operator in an explicit matrix equation form and then taking the adjoint to this matrix equation followed by determining the corresponding operator. The reverse time migration operator as the exact adjoint to the implied forward operator can be used not only as a migration algorithm but also as an adjoint operator which is required in the imaging through an inversion such as least-squares migration.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.430-437
• /
• 2010

In this paper, we develop a shape design optimization method for thermo-elastoplasticity problems that is applicable to the welding or thermal deformation problems of ship structures. Shell elements and a programming language APDL in a commercial finite element analysis code, ANSYS, are employed in the shape optimization. The point of developed method is to determine the design parameters such that the deformed shape after welding fits very well to a desired design. The geometric parameters of surfaces are selected as the design parameters. The modified method of feasible direction (MMFD) and finite difference sensitivity are used for the optimization algorithm. Two numerical examples demonstrate that the developed shape design method is applicable to existing hull structures and effective for the structural design of ships.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.