• Economic and Environmental Geology
• /
• v.49 no.1
• /
• pp.31-41
• /
• 2016

A program code was developed to calculate block hydraulic conductivity of the 2-D DFN(discrete fracture network) system based on equivalent pipe network, and implemented to examine the effect of joint orientation distribution on the hydraulic characteristics of fractured rock masses through numerical experiments. A rock block of size $32m{\times}32m$ was used to generate the DFN systems using two joint sets with fixed input parameters of joint frequency and gamma distributed joint size, and various normal distributed joint trend. DFN blocks of size $20m{\times}20m$ were selected from center of the $32m{\times}32m$ blocks to avoid boundary effect. Twelve fluid flow directions were chosen every $30^{\circ}$ starting at $0^{\circ}$. The directional block conductivity including the theoretical block conductivity, principal conductivity tensor and average block conductivity were estimated for generated 180 2-D DFN blocks. The effect of joint orientation distribution on block hydraulic conductivity and chance for the equivalent continuum behavior of the 2-D DFN system were found to increase with the decrease of mean intersection angle of the two joint sets. The effect of variability of joint orientation on block hydraulic conductivity could not be ignored for the DFN having low intersection angle between two joint sets.

• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.373-384
• /
• 2000

The goal of this paper is to develop a computational method to predict cancellous bone density distributions based upon continuum levels of volumetric strain. Volumetric strain is defined as the summation of normal strains, excluding shear strains, within an elastic range of loadings. Volumetric strain at a particular location in a cancellous structure changes with changes of the boundary conditions (prescribed displacements, tractions, and pressure). This change in the volumetric strain is postulated to predict the adaptive change in the bone apparent density. This bone remodeling theory based on volumetric strain is then used with the finite element method to compute the apparent density distribution for cancellous bone in both lumbar spine and proximal femur using an iterative algorithm, considering the dead zone of strain stimuli. The apparent density distribution of cancellous bone predicted by this method has the same pattern as experimental data reported in the literature (Wolff 1892, Keller et al. 1989, Cody et al. 1992). The resulting bone apparent density distributions predict Young's modulus and strength distributions throughout cancellous bone in agreement with the literature (Keller et al. 1989, Carter and Hayes 1977). The method was convergent and sensitive to changes in boundary conditions. Therefore, the computational algorithm of the present study appears to be a useful approach to predict the apparent density distribution of cancellous bone (i.e. a numerical approximation for Wolff's Law)

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.4
• /
• pp.425-431
• /
• 2015

Based on a bond-based peridynamics theory for dynamic crack propagation problems, this paper presents a design sensitivity analysis and optimization method. Peridynamics has a peculiar advantage over the existing continuum theory in the mathematical modelling of problems where discontinuities arise. For the design optimization of the crack propagation problems, a non-shape design sensitivity is derived using the adjoint variable method. The obtained adjoint sensitivity of displacement and strain energy turns out to be very accurate and efficient compared to the finite different sensitivity. The obtained design sensitivities are futher utilized to optimally control the position of bifurcation point in the design optimization of crack propagation in a plate under tension. A numerical experiment demonstrates that the optimal distribution of material density could delay the position of bifurcation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.1
• /
• pp.82-87
• /
• 1987

In this study a method of shape optimization is applied to two dimensional heat transfer system. For this the optimization problem is defined in a functional form including cost, constraints and the system governing equation. Then the material derivative concept in continuum mechanics and the adjoint variable method are employed for the shape design sensitivity analysis. With the sensitivity analysis results, an optimum is sought with the gradient projection optimization algorithm. The two dimensional isoparametric finite elements are used for accurate analysis and sensitivity calculations. The above method is employed to find the boundary shape to achieve a desired temperature distribution along a segment of the boundary subject to the maximum area constraint.

• Tunnel and Underground Space
• /
• v.26 no.4
• /
• pp.283-292
• /
• 2016

A computer program code was developed to estimate hydraulic behavior of the 2-D connected pipe network system, and implemented to evaluate the effect of joint aperture on hydraulic parameters of fractured rock masses through numerical experiments. A total of 216 stochastic 2-D DFN(discrete fracture network) blocks of $20m{\times}20m$ were prepared using two joint sets with fixed input parameters of joint orientation, frequency and size distribution. Two different cases of joint aperture variation are considered in this study. The hydraulic parameters were estimated for generated 2-D DFN blocks. The hydraulic anisotropy and the chance for equivalent continuum behavior of the DFN system were found to depend on the variability of joint aperture.

• Steel and Composite Structures
• /
• v.28 no.5
• /
• pp.589-606
• /
• 2018

The previous studies reflected the significant effect of neutral-axis position and coupling of in-plane and out-of-plane displacements on behavior of functionally graded (FG) nanobeams. In thin FG beam, this coupling can be eliminated by a proper choice of the reference axis. In shear deformable FG nanobeam, not only this coupling can't be eliminated but also the position of neutral-axis is dependent on through-thickness distribution of shear strain. For the first time, in this paper it is avoided to guess a shear strain shape function and the exact shape function and consequently the exact position of neutral axis for arbitrary gradation of higher order nanobeam are obtained. This paper presents new methodology based on differential transform and collocation methods to solve coupled partial differential equations of motion without any simplifications. Using exact position of neutral axis and higher order beam kinematics as well as satisfying equilibrium equations and traction-free conditions without shear correction factor requirement yields to better results in comparison to the previously published results in literature. The classical rule of mixture and Mori-Tanaka homogenization scheme are considered. The Eringen's nonlocal continuum theory is applied to capture the small scale effects. For the first time, the dependency of exact position of neutral axis on length to thickness ratio is investigated. The effects of small scale, length to thickness ratio, Poisson's ratio, inhomogeneity of materials and various end conditions on vibration and buckling of local and nonlocal FG beams are investigated. Moreover, the effect of axial load on natural frequencies of the first modes is examined. After degeneration of the governing equations, the exact new formulas for homogeneous nanobeams are computed.

• Journal of The Korean Astronomical Society
• /
• v.45 no.4
• /
• pp.93-99
• /
• 2012

We present a 1.1mm emission map of the OMC1 region observed with AzTEC, a new large-format array composed of 144 silicon-nitride micromesh bolometers, that was in use at the James Clerk Maxwell Telescope (JCMT). These AzTEC observations reveal dozens of cloud cores and a tail of filaments in a manner that is almost identical to the submillimeter continuum emission of the entire OMC1 region at 450 and $850{\mu}m$. We perform Fourier analysis of the image with a modified periodogram and the density power spectrum, which provides the distribution of the length scale of the structures, is determined. The expected value of the periodogram converges to the resulting power spectrum in the mean squared sense. The present analysis reveals that the power spectrum steepens at relatively smaller scales. At larger scales, the spectrum flattens and the power law becomes shallower. The power spectra of the 1.1mm emission show clear deviations from a single power law. We find that at least three components of power law might be fitted to the calculated power spectrum of the 1.1mm emission. The slope of the best fit power law, ${\gamma}{\approx}-2.7$ is similar to those values found in numerical simulations. The effect of beam size and the noise spectrum on the shape and slope of the power spectrum are also included in the present analysis. The slope of the power law changes significantly at higher spatial frequency as the beam size increases.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.2
• /
• pp.75.2-75.2
• /
• 2014

CS molecule as an important tracer for studying inward motions in dense cores is known to be adsorbed onto dusts in cold (T~10K) dense cores, resulting in its significant depletion in the central region of the cores which may hamper a proper study of kinematics stage of star formation. In this study we choose five 'evolved' dense starless cores, L1544, L1552, L1689B, L694-2 and L1197, to investigate how depletion of CS molecule is significant and how the molecule differentiates depending on the evolutional status of the dense cores, by using a rare isotopomer C34S. We performed mapping observations in C34S (J=2-1) and N2H+ (J=1-0) with Nobeyama 45 m telescope, and compared $850{\mu}m$ continuum data as a reference of the density distribution of the dense cores. Our data confirm the claim that CS molecule generally depletes out in the central region in dense starless cores, while N2H+ keeps abundant as they get evolved. All of integrated intensity maps show 'semi-ring-like' depletion holes in CS, and all of abundance radial profiles show decrease toward center. The CS depletion and molecular chemical differentiation seems to depend on the evolutional status in dense cores. The evolved cores shows low abundance at both central and outer regions, implying that in the case of highly evolved cores CS freeze-out occurs over the most area of the cores.

• The Bulletin of The Korean Astronomical Society
• /
• v.42 no.1
• /
• pp.54.1-54.1
• /
• 2017

The UKIRT Wide-field Infrared survey for Fe+(UWIFE) is an unbiased survey of the first Galactic quadrant, with narrow-band filter centered on $1.644{\mu}m$. This survey covers $7^{\circ}$ < l < $62^{\circ}$ and |b| < $1.5^{\circ}$, where active interaction of stars and interstellar medium is expected. With median seeing of 0.8 arcsec, 5 - sigma detection limit of 18.7 mag and surface brightness limit of $8.1{\times}10^{-20}W\;m^{-2}arcsec^{-2}$, this survey gives an opportunity to statistically study Galactic [Fe II] - emitting sources for the first time. In order to identify Ionized Fe Objects (IFOs) in survey area systematically, we conducted visual inspection and automatic detection simultaneously. Total of ~300 extended IFOs are identified, most of them are found out to be part of supernova remnants (SNRs), young stellar objects, HII regions and planetary nebulae. The majority of IFOs are new discoveries which reveal shocked structures in high-extinction region. Spatial distribution of IFOs suggest that they trace Galactic structure. As a part of spectroscopic follow-up, we observed SNR candidate IFO J183740.829-061452.41 with IGRINS (Immersion Grating Infrared Spectrograph, Yuk+2010), mounted on 2.7m Harlan Smith telescope. This unknown arc-like, 6'-long IFO is coincident with inner part of radio continuum loop G25.8+0.2, which has been known as HII region. However, interior of this radio shell is filled with diffuse soft X-ray emission, and possible association of hard X-ray pulsar / pulsar wind nebula makes the nature of the IFO unclear. The H and K-band 2D spectrum shows shock-ionized [Fe II] filaments, which is apart from photoionized HII filaments. In this presentation we present basic statistics of newly identified IFOs, as well as the follow-up study of IFO J183740.829-061452.41.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.7 no.4
• /
• pp.335-341
• /
• 2005

In construction of a large cross section NATM tunnel, to keep the tunnel face stability by the ground itself bench cut method is commonly used. In order to necessity of partial face excavation method, we have to look for more enhanced method that can maintain better stress intensity. This paper presents a stress distribution of the Center Diaphragm Method from the partial face excavation methods, with the numerical analysis, and induced the optimal face distance, which is minimizing stress concentration and the optimal excavation step. Commerical 3 dimensional continuum analyzing FLAC-3D Ver. 2.1 program is used for the analysis. Analyses were performed to investigate ground behavior for tunnels with variable bench-length varying from 2m to 40m.