• Geomechanics and Engineering
• /
• v.36 no.2
• /
• pp.131-143
• /
• 2024

Shield tunneling construction commonly crosses underground pipelines in urban areas, resulting in soil loss and followed deformation of grounds and pipelines nearby, which may threaten the safe operation of shield tunneling. This paper investigated the pipeline deformation caused by double curved shield tunnels in soil-rock composite stratum in Nanjing, China. The stratum settlement equation was modified to consider the double shield tunneling. Moreover, a three dimensional finite element model was established to explore the effects of hard-layer ratio, tunnel curvature radius, pipeline buried depth and other influencing factors. The results indicate the subsequent shield tunnel would cause secondary disturbance to the soil around the preceding tunnel, resulting in increased pipeline and ground surface settlement above the preceding tunnel. The settlement and stress of the pipeline increased gradually as buried depth of the pipeline increased or the hard-layer ratio (the ratio of hard-rock layer thickness to shield tunnel diameter within the range of the tunnel face) decreased. The modified settlement calculation equation was consistent with the measured data, which can be applied to the settlement calculation of ground surface and pipeline settlement. The modified coefficients a and b ranged from 0.45 to 0.95 and 0.90 to 1.25, respectively. Moreover, the hard-layer ratio had the most significant influence on the pipeline settlement, but the tunnel curvature radius and the included angle between pipeline and tunnel axis played a dominant role in the scope of the pipeline settlement deformation.

• The Journal of the Institute of Internet, Broadcasting and Communication
• /
• v.9 no.5
• /
• pp.139-146
• /
• 2009

Parametric array for audio applications is analyzed by numerical modeling and analytical approximation. The nonlinear wave equations are used to provide design guidelines for the audio parametric array. A time domain finite difference code that accurately solves the KZK (Khokhlov-Zabolotskaya-Kuznetsov) nonlinear parabolic wave equation is used to predict the response of the parametric array. The time domain code relates the source size and the carrier frequency to the audible signal response including the output level and beamwidth to considering the implementation issues for audio applications of the parametric array, the emphasis is given to the frequency response and distortion. We use the time domain code to find out the optimal parameters that will help produce the parametric array with highest achievable output in terms of the average power within the demodulated signal. Parameters such as primary input frequency, audio source radius and the modulation method are given utmost importance. The output effect of those parameters are demonstrated through the numerical simulation.

• Journal of The Korean Astronomical Society
• /
• v.35 no.4
• /
• pp.159-174
• /
• 2002

Cosmological hydrodynamic simulations of large scale structure in the universe have shown that accretion shocks and merger shocks form due to flow motions associated with the gravitational collapse of nonlinear structures. Estimated speed and curvature radius of these shocks could be as large as a few 1000 km/s and several Mpc, respectively. According to the diffusive shock acceleration theory, populations of cosmic-ray particles can be injected and accelerated to very high energy by astrophysical shocks in tenuous plasmas. In order to explore the cosmic ray acceleration at the cosmic shocks, we have performed nonlinear numerical simulations of cosmic ray (CR) modified shocks with the newly developed CRASH (Cosmic Ray Amr SHock) numerical code. We adopted the Bohm diffusion model for CRs, based on the hypothesis that strong Alfven waves are self-generated by streaming CRs. The shock formation simulation includes a plasma-physics-based 'injection' model that transfers a small proportion of the thermal proton flux through the shock into low energy CRs for acceleration there. We found that, for strong accretion shocks, CRs can absorb most of shock kinetic energy and the accretion shock speed is reduced up to $20\%$, compared to pure gas dynamic shocks. For merger shocks with small Mach numbers, however, the energy transfer to CRs is only about $10-20\%$ with an associated CR particle fraction of $10^{-3}$. Nonlinear feedback due to the CR pressure is insignificant in the latter shocks. Although detailed results depend on models for the particle diffusion and injection, these calculations show that cosmic shocks in large scale structure could provide acceleration sites of extragalactic cosmic rays of the highest energy.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.6
• /
• pp.111-120
• /
• 1999

Leachate flow behavior due to intermediate cover soil of low hydraulic conductivity and the applicability of pumping method for reducing the leachate level in the landfill are analyzed with the numerical flow model, MODFLOW. Using the hydraulic conductivity and storativity data obtained from the field pumping and slug tests(Jang and Cho, 1999), the hydraulic condition within the landfill is validated. The optimum rate of pumping, the radius of influence, and the efficiency of horizontal drain are analyzed for reducing the leachate level in the landfill. From the results of the analyses, the barrier effect that the buried cover soil of low hydraulic conductivity prevents the vertical movement of leachate flow through the cover soil, which is found from the in-situ geotechnical studies(Jang and Cho, 1999), is identified again. Also, the installation of horizontal drains to the pumping well can increase the pumping rate from 120 ton/day per a well to 300 ton/day. The length of horizontal drain did not influence significantly on the drawdown-time curve of leachate in the landfill.

• Atmosphere
• /
• v.23 no.1
• /
• pp.1-11
• /
• 2013

The angular momentum transport of an idealized axisymmetric vortex in the developing stage was investigated using the Weather Research and Forecast (WRF) model. The balanced axisymmetric vortex was constructed based on an empirical function for tangential wind, and the temperature, geopotential, and surface pressure were obtained from the balanced equation. The numerical simulation was carried out for 6 days on the f-plane with the Sea Surface Temperature (SST) set as constant. The weak vortex at initial time was intensified with time, and reached the strength of tropical cyclone in a couple of days. The Absolute Angular Momentum (AAM) was transported along with the secondary circulation of the vortex. Total AAM integrated over a cylinder of radius of 2000 km decreased with simulation time, but total kinetic energy increased rapidly. From the budget analysis, it was found that the surface friction is mainly responsible for the decrease of total AAM. Also, contribution of the surface friction to the AAM loss was about 90% while that of horizontal advection was as small as 8%. The trajectory of neutral numerical tracers following the secondary circulation was presented for the Lagrangian viewpoint of the transports of absolute angular momentum. From the analysis using the trajectory of tracers it was found that the air parcel was under the influence of the surface friction continuously until it leaves the boundary layer near the core. Then the air parcel with reduced amount of angular momentum compared to its original amount was transported from boundary layer to upper level of the vortex and contributed to form the anti-cyclone. These results suggest that the tropical cyclone loses angular momentum as it develops, which is due to the dissipation of angular momentum by the surface friction.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.5
• /
• pp.471-482
• /
• 2016

In this study, we focused on a speed-up of KIESSI-3D program, which is based on FE-IE techniques, by introducing a p-version dynamic infinite element method. In order to evaluate performance of the KIESSI-3D, numerical analyses for eight real-scale SSI problems are carried out. We considered three types of KIESSI-3D numerical models whose radii of near-field soil region($r_0$)are 1.2, 1.5, and 3.0 times of basemat radius of structure(R). In addition, SSI analyses using the SASSI2010 program are carried out used for comparison of accuracy and runtime against those of the KIESSI-3D. Numerical results show that the KIESSI-3D model of $r_0=1.2R$ is enough to give accurate solution. In view of the computing speed, the new KIESSI-3D was up to 25 times faster than the old KIESSI-3D.

• The Journal of the Petrological Society of Korea
• /
• v.26 no.1
• /
• pp.83-91
• /
• 2017

In order to determine the runout range of pyroclastic density currents on Jeju island, lava dome collapse on 8 locations of outer rim of Baekrokdam crater were simulated by TITAN2D numerical simulation program. We set parameters as internal friction angle as $30^{\circ}$ and bed friction angle as $20^{\circ}$ to control velocity of currents occurred by lava dome collapse. Then we set the height and radius of lava dome, initial speed of collapse and simulation times. And we carried out numerical simulations for a total of 96 scenarios. The result shows that the maximum runout distance was 13.4 km in case of lava dome collapse. This study can be used database for manufacturing of hazard map to minimize damages caused by pyroclastic density currents occurred on Jeju island.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.1 s.24
• /
• pp.29-38
• /
• 2007

The research on the track performance and stability of the tilting-train was performed and the settlement of the roadbed was estimated as the tilting train was being operated on the rail joint under the allowable velocity subjected to the track performance and the stability of the tilting-train. Since the impact on the continuous welded rail (CWR) induced by the tilting-train loading is different from the impact on the rail joint, it needs to investigate the settlement of the roadbed beneath the CWR. In this study, when the tilting-train is being operated on the CWR under the allowable velocity subjected to the track performance and the stability of the tilting-train, the settlement and bearing capacity of the roadbed beneath the CWR have been evaluated using numerical analysis and compared with those beneath the rail joint. The numerical results show that the settlements of the roadbed beneath CWR and rail joint are amount to 71.2% and 88.8% of the allowable settlement, respectively. And the stresses are amount to 10.4% and 12.1% of the allowable bearing capacity, respectively.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.43 no.4
• /
• pp.469-476
• /
• 2023

Soil plasticity models for cyclic and dynamic loads are essential in non-linear numerical analysis of geotechnical structures. While a single yield surface model shows a linear behavior for cyclic loads, J₂-bounding surface plasticity model with zero elastic region can effectively simulate a nonlinearity of the ground response with the same material properties. The radius of the yield surface inside the boundary surface converged to 0 to make the elastic region disappear, and plastic hardening modulus and dilatancy define plastic strain increment. This paper presents the stress-strain incremental equation of the developed model, and derives plastic hardening modulus for the hyperbolic model. The comparative analyses of the triaxial compression test and the shallow foundation under the cyclic load can show stable numerical convergence, consistency with the theoretical solution, and hysteresis behavior. In addition, plastic hardening modulus for the modified hyperbolic function is presented, and a methodology to estimate model variables conforming 1D equivalent linear model is proposed for numerical modeling of the multi-dimensional behavior of the ground.

• Geomechanics and Engineering
• /
• v.34 no.4
• /
• pp.397-408
• /
• 2023

Pipe pile walls are commonly used as retaining structures for excavation projects, particularly in densely populated coastal cities such as Hong Kong. Pipe pile walls are preferred in reclaimed land due to their cost-effectiveness and convenience for installation. However, the pre-bored piling techniques used to install pipe piles can cause significant ground disturbance, posing risks to nearby sensitive structures. This study reports a well-documented case history in a reclamation site, and it was found that pipe piling could induce ground settlement of up to 100 mm. Statutory design submissions in Hong Kong typically specify a ground settlement alarm level of 10 mm, which is significantly lower than the actual settlement observed in this study. In addition, lateral soil movement of approximately 70 mm was detected in the marine deposit. The lateral soil displacement in the marine deposit was found to be up to 3.4 and 3.1 times that of sand fill and CDG, respectively, mainly due to the relatively low stiffness of the marine deposit. Based on the monitoring data and site-investigation data, a 3D numerical analysis was established to back-analyze soil movements due to the installation of the pipe pile wall. The comparison between measured and computed results indicates that the equivalent ground loss ratio is 20%, 40%, and 20% for the fill, marine deposit and CDG, respectively. The maximum ground settlement increases with an increase in the ground loss ratio of the marine deposit, whereas the associated influence radius remains stationary at 1.2 times the pipe pile wall depth (H). The maximum ground settlement increases rapidly when the thickness of marine deposit is less than 0.32H, particularly for the ground loss ratio of larger than 40%. This study provides new insights into the pipe piling construction in reclamation sites.