• International Journal of Aeronautical and Space Sciences
• /
• v.15 no.3
• /
• pp.232-240
• /
• 2014

Transonic flow past a NASA SC(2)-0710 airfoil with deployments of a spoiler up to $6^{\circ}$ was studied numerically. We consider angles of attack from $-0.6^{\circ}$ to $0.6^{\circ}$ and free-stream Mach numbers from 0.81 to 0.86. Solutions of the unsteady Reynolds-averaged Navier-Stokes equations were obtained with a finite-volume solver using several turbulence models. Both stationary and time-dependent deployments of the spoiler were examined. The study revealed the existence of narrow bands of the Mach number, angle of attack, and spoiler deflection angle, in which the flow was extremely sensitive to small perturbations. Simulations of 3D flow past a swept wing confirmed the flow sensitivity to small perturbations of boundary conditions.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.5
• /
• pp.434-441
• /
• 2016

Numerical simulations of 1D Burgers equation and 2D sloshing problem were carried out to study numerical discrepancy between serial and parallel computations. The numerical domain was decomposed into 2 and 4 subdomains for parallel computations with message passing interface. The numerical solution of Burgers equation disclosed that fully explicit boundary conditions used on subdomains of parallel computation was responsible for the numerical discrepancy of transient solution between serial and parallel computations. Two dimensional sloshing problems in a rectangular domain were solved using OpenFOAM. After a lapse of initial transient time sloshing patterns of water were significantly different in serial and parallel computations although the same numerical conditions were given. Based on the histograms of pressure measured at two points near the wall the statistical characteristics of numerical solution was not affected by the number of subdomains as much as the transient solution was dependent on the number of subdomains.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.1
• /
• pp.49-62
• /
• 1994

An analysis method for the time-dependent nonlinear analysis of segmentally erected planar prestressed concrete cable-stayed bridges was described. To account for the time-dependent effects, load history, creep, shrinkage. aging of concrete and relaxation of prestress were considered. Changes in boundary conditions and loads, installing and removing frame elements, stressing, restressing and removing cables and prestressing tendons were incorporated for modeling segmental erection operations. One typical example on segmentally erected prestressed concrete cable-stayed bridge was presented to illustrate the analysis method. Results of this example show that it is important to follow the development of stresses and deformations at all stages of construction to predict the true response of the bridge through its various load history.

• Steel and Composite Structures
• /
• v.32 no.2
• /
• pp.239-252
• /
• 2019

Since conical sandwich shells are important structures in the modern industries, in this paper, for the first time, vibration behavior of the truncated conical sandwich shells which include temperature dependent porous FG face sheets and temperature dependent homogeneous core in various thermal conditions are investigated. A high order theory of sandwich shells which modified by considering the flexibility of the core and nonlinear von Karman strains are utilized. Power law rule which modified by considering the two types of porosity volume fractions are applied to model the functionally graded materials. By utilizing the Hamilton's energy principle, and considering the in-plane and thermal stresses in the face-sheets and the core, the governing equations are obtained. A Galerkin procedure is used to solve the equations in a simply supported boundary condition. Uniform, linear and nonlinear temperature distributions are used to model the effect of the temperature changing in the sandwich shell. To verify the results of this study, they are compared with FEM results obtained by Abaqus software and for special cases with the results in literatures. Eigen frequencies variations are surveyed versus the temperature changing, geometrical effects, porosity, and some others in the numerical examples.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2003.05a
• /
• pp.427-432
• /
• 2003

An analytical method to predict the flexural behavior of composite girder is presented in which the early-age properties of concrete are specified including maturing of elastic modulus, creep and shrinkage. The time dependent constitutive relation accounting for the early-age concrete properties is derived in an incremental format by expanding the total form of stress-strain relation by the first order Taylor series with respect to the reference time. The sectional analysis calculates the axial and curvature strains based on the force and moment equilibriums. The deflection curve of the box 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 girders which is a composite bridge girder filled with concrete at the bottom of the steel box in the negative moment region. The one dimensional finite element analysis results are compared with those of the three dimensional finite element analysis and the analytical method based on the sectional analysis. Close agreement is observed among the three methods.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.328-337
• /
• 2005

This work reports results of our study on the dynamic responses of the buried pipelines both along the axial and the transverse directions under various boundary end conditions. We have considered three cases, i.e., the free ends, the fixed ends, and the fixed-free ends for the axial direction, and three more cases including the guided ends, the simply supported ends, and the supported-guided ends for the transverse direction. In order to investigate the effect of the boundary end conditions for the dynamic responses of the buried pipeline, we have devised a computer program to find the solutions of the formulae on the dynamic responses (displacements, axial strains, and bending strains) under the various boundary end conditions considered in this study. The dynamic behavior of the buried pipelines for the forced vibration is found to exhibit two different forms, a transient response and a steady state response, depending on the time before and after the transfer of a seismic wave on the end of the buried pipeline. The former is identified by a slight change in its behavior before the sinusoidal-shaped seismic wave travels along the whole length of the pipeline whereas the latter by the complete form of a sinusoidal wave when the wave travels throughout the pipeline. The transient response becomes insignificant as the wave speed increases. We have observed a resonance when the mode wavelength matches the wavelength of the seismic wave, where the mode number(k) of resonance for the axial direction is found to be $\overline{\omega}/{\pi}V+1/2$ for the fixed-free ends, $\overline{\omega}/{\pi}V+1$ for the free ends, and $\overline{\omega}/{\pi}V$ for the fixed ends, respectively. By adding 10 more modes to the mode number(k) of resonance, we were able to study all the dynamic responses of the buried pipeline for the axial direction. On the other hand, we have not been able to observe a resonance in the analysis for the transverse direction, because the dynamic responses are found to vanish after the seventh mode. From the results of the dynamic responses at the many points of the pipeline, we have found that the responses appeared to be dependent critically on the boundary end conditions. Such effects are found to be most prominent especially for the maximum values of the displacement and the strain and its position.

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

This work reports results of our study on the dynamic responses of the buried pipelines both along the axial and the transverse directions under various boundary end conditions. In order to investigate the effect of the boundary end conditions for the dynamic responses of the buried pipeline, we have devised a computer program to find the solutions of the formulae on the dynamic responses (displacements, axial strains, and bending strains) under the various boundary end conditions considered in this study, The dynamic behavior of the buried pipelines for the forced vibration is found to exhibit two different forms, a transient response and a steady state response, depending on the time before and after the transfer of a seismic wave on the end of the buried pipeline. We have observed a resonance when the mode wavelength matches the wavelength of the seismic wave, where the mode number(k) of resonance f3r the axial direction. On the other hand, we have not been able to observe a resonance in the analysis of the transverse direction, because the dynamic responses are found to vanish after the seventh mode. From the results of the dynamic responses at many points of the pipeline, we have found that the responses appeared to be dependent critically on the boundary end conditions. Such effects are found to be most prominent especially for the maximum values of the displacement, the strain and its position.

• Journal of the Korean GEO-environmental Society
• /
• v.13 no.3
• /
• pp.39-51
• /
• 2012

The time-dependent characteristics of mudstone in this study were analyzed by performing creep tests for the mudstone in Pohang. The creep tests were conducted on various boundary conditions to observe the time-dependent behavior of mudstone and the results were compared with values predicted by established visco-elastic model and empirical equations. As a result of analysis for empirical equations proposed by Griggs(1936), Cottrell(1952) and Singh(1975), the creep constants generally tend to increase as increasing the stress level and water content. Moreover, the values predicted by Singh's equation were well fitted for the test results. Therefore, it is expected to be reasonable that creep behavior of mudstone is analyzed using Singh's equation. As a result of analysis for a Burger's model, the predicted values were well fitted for the test results. Therefore, it is expected be reasonable that the creep behavior of mudstone by a rheological model is analyzed using Burger's model.

• KIEAE Journal
• /
• v.16 no.1
• /
• pp.81-87
• /
• 2016

Purpose: Heat transfers phenomena are described by the second order partial differential equation and its boundary conditions. In a three-dimensional structure of a building, the heat transfer phenomena generally include more than one material, and thus, become complicate. The analytic solutions are useful to understand heat transfer phenomena, but they can hardly be applied in engineering or design problems. Engineers and designers have generally been forced to use numerical methods providing reliable results. Finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains. Method: To obtain an numerical solution, a discretization method, which approximates the differential equations, and the interpolation methods for temperature and heat flux between two or more materials are required. The discretization methods are applied to small domains in space and time, and these numerical solutions form the descretized equations provide approximated solutions in both space and time. The accuracy of numerical solutions is dependent on the quality of discretizations and size of cells used. The higher accuracy, the higher numerical resources are required. The balance between the accuracy and difficulty of the numerical methods is critical for the success of the numerical analysis. A simple and easy interpolation methods among multiple materials are developed. The linear equations are solved with the BiCGSTAB being a effective matrix solver. Result: This study provides an overview of discretization methods, boundary interface, and matrix solver for the 3-dimensional numerical heat transfer including two materials.

• The Bulletin of The Korean Astronomical Society
• /
• v.40 no.1
• /
• pp.65.2-65.2
• /
• 2015

We developed a three-dimensional magnetohydrodynamic (MHD) code to reproduce the structure of a solar wind, the properties of a coronal mass ejection (CME) and the interaction between them. This MHD code is based on the finite volume method incorporating total variation diminishing (TVD) scheme with an unstructured grid system. In particular, this grid system can avoid the singularity at the north and south poles and relax tight CFL conditions around the poles, both of which would arise in a spherical coordinate system (Tanaka 1994). In this model, we first apply an MHD tomographic method (Hayashi et al. 2003) to interplanetary scintillation (IPS) observational data and derive a solar wind from the physical values obtained at 50 solar radii away from the Sun. By comparing the properties of this solar wind to observational data obtained near the Earth orbit, we confirmed that our model captures the velocity, temperature and density profiles of a solar wind near the Earth orbit. We then insert a spheromak-type CME (Kataoka et al. 2009) into the solar wind to reproduce an actual CME event. This has been done by introducing a time-dependent boundary condition to the inner boundary of our simulation domain. On the basis of a comparison between a simulated CME and observations near the Earth, we discuss the physics involved in an ICME interacting with a solar wind.