• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.945-956
• /
• 1999

This work reports a set of approximate analytical solutions describing the initial transient process of close-contact melting between a rectangular parallelepiped solid and a flat plate on which either constant temperature or constant heat flux is imposed. Not only relative motion of the solid block tangential to the heating plate, but also the density difference between the solid and liquid phase is incorporated in the model. The thin film approximation reduces the force balance between the solid weight and liquid pressure, and the energy balance at the melting front into a simultaneous ordinary differential equation system. The normalized model equations admit compactly expressed analytical solutions which include the already approved two-dimensional solutions as a subset. In particular, the normalized liquid film thickness is independent of all pertinent parameters, thereby facilitating to define the transition period of close-contact melting. A unique behavior of the solid descending velocity due to the density difference is also resolved by the present solution. A new geometric function which alone represents the three-dimensional effect is introduced, and its properties are clarified. One of the representative results is that heat transfer is at least enhanced at the expense of the increase in friction as the cross-sectional shape deviates from the square under the same contact area.

• Journal of The Korean Astronomical Society
• /
• v.38 no.2
• /
• pp.59-66
• /
• 2005

A new ion transport code for planetary ionospheric studies has been developed with consideration of velocity differences among ion species involving ion-ion collision. Most of previous planetary ionosphere models assumed that ions diffuse through non-moving ion and neutral background in order to consolidate continuity and momentum equations for ions into a simple set of diffusion equations. The simplification may result in unreliable density profiles of ions at high altitudes where ion velocities are fast and their velocity differences are significant enough to cause inaccuracy when computing ion-ion collision. A new code solves explicitly one-dimensional continuity and momentum equations for ion densities and velocities by utilizing divided Jacobian matrices in matrix inversion necessary to the Newton iteration procedure. The code has been applied to Martian nightside ionosphere models, as an example computation. The computed density profiles of $O^+,\;OH^+$, and $HCO^+$ differ by more than a factor of 2 at altitudes higher than 200 km from a simple diffusion model, whereas the density profile of the dominant ion, $O_2^+$, changes little. Especially, the density profile of $HCO^+$ is reduced by a factor of about 10 and its peak altitude is lowered by about 40 km relative to a simple diffusion model in which $HCO^+$ ions are assumed to diffuse through non-moving ion background, $O_2^+$. The computed effects of the new code on the Martian nightside models are explained readily in terms of ion velocities that were solved together with ion densities, which were not available from diffusion models. The new code should thus be expected as a significantly improved tool for planetary ionosphere modelling.

• Geomechanics and Engineering
• /
• v.31 no.1
• /
• pp.53-69
• /
• 2022

In densely built areas, the development of underground transportation systems often involves twin excavations, which are sometimes unavoidably constructed adjacent to existing piled foundations. Because soil stiffness degrades with induced stress release and shear strain during excavation, it is vital to investigate the piled raft responses to subsequent excavation after the first tunnel in a twin-excavation system. The effects of deep excavations on existing piled foundations have been extensively investigated, but the influence of twin excavations on a piled raft is seldom reported in the literature. In this study, three-dimensional numerical analyses were carried out to investigate the influence of sand density on an existing piled raft (with a working load on top of the raft) due to twin excavations. A wide range of relative density (D_r) from loosest (30%), loose to medium (50% and 70%), and densest (90%) were selected to investigate the effects on settlement and load transfer mechanism of the piled raft during twin excavations. An advanced hypoplastic sand model (which can capture small-strain stiffness and stress-state dependent dilatancy of sand) was adopted. The model parameters are calibrated against centrifuge test results in sand reported in the literature. From the computed results, it is found that twin excavations in loose sand (D_r=30%) caused the most significant settlement. This is because of the higher stiffness of denser sand (D_r=90%) than that of loose sand. In contrast, a much larger tilting (maximum magnitude=0.18%) was computed in dense sand than in loose sand after the completion of the first excavation. As far as the load transfer mechanism along the piles is concerned, an upward load transfer to mobilize shaft resistance is observed in loose sand. On the contrary, a downward load transfer is observed in dense sand.

• Communications for Statistical Applications and Methods
• /
• v.25 no.2
• /
• pp.131-157
• /
• 2018

We propose a robust event date model to estimate the date of a target event by a combination of individual dates obtained from archaeological artifacts assumed to be contemporaneous. These dates are affected by errors of different types: laboratory and calibration curve errors, irreducible errors related to contaminations, and taphonomic disturbances, hence the possible presence of outliers. Modeling based on a hierarchical Bayesian statistical approach provides a simple way to automatically penalize outlying data without having to remove them from the dataset. Prior information on individual irreducible errors is introduced using a uniform shrinkage density with minimal assumptions about Bayesian parameters. We show that the event date model is more robust than models implemented in BCal or OxCal, although it generally yields less precise credibility intervals. The model is extended in the case of stratigraphic sequences that involve several events with temporal order constraints (relative dating), or with duration, hiatus constraints. Calculations are based on Markov chain Monte Carlo (MCMC) numerical techniques and can be performed using ChronoModel software which is freeware, open source and cross-platform. Features of the software are presented in Vibet et al. (ChronoModel v1.5 user's manual, 2016). We finally compare our prior on event dates implemented in the ChronoModel with the prior in BCal and OxCal which involves supplementary parameters defined as boundaries to phases or sequences.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1165-1173
• /
• 1997

The steady state close-contact melting phenomenon occurring between a phase change material and an isothermally heated flat plate with relative motion is investigated analytically, in which the effects of vertical convection in the liquid film and solid-liquid density difference are incorporated simultaneously. Not only the scale analysis is conducted to estimate a priori qualitative dependence of system variables on characteristic parameters, but also an analytical solution to a set of simplified model equations is obtained to specify the effects under consideration. These two results are consistent with each other, in that the vertical convection affects both the solid descending velocity and the film thickness, and that the density difference alters only the solid descending velocity. While the effect of vertical convection can be characterized conveniently by a newly introduced temperature gradient factor which asymptotically approaches the unity/zero with decreasing/increasing the Stefan number, that of density difference is represented by the liquid-to-solid density ratio. It is shown that the solid descending velocity depends linearly on the density ratio, and that the ratios of solid descending velocity, film thickness and friction coefficient to the conduction solution are proportional to 3/4, 1/4 and -1/4 powers of the temperature gradient factor, respectively. Also, established is the fact that the effect of convection can be legitimately neglected in the analysis for the range of the Stefan number less than 0.1.

• Journal of Industrial Technology
• /
• v.24 no.B
• /
• pp.107-116
• /
• 2004

This thesis is results of experimental works on the behavior of the cut-and-cover tunnel. Centrifuge model tests were performed to simulate the behavior of the cut-and-cover tunnels having cross sections of national road and subway tunnels. Model experiments were carried out with changing the cut slope and the slope of filling ground surface. Displacements of tunnel lining resulted from artificially accelerated gravitational force up to 40g of covered material used in model tests, were measured during centrifuge model tests. In model tests, Jumunjin Standard Sand with the relative density of 80 % and the zinc plates were used for the covered material and the flexible tunnel lining, respectively. Basic soil property tests were performed to obtain it's the property of Jumumjin Standard Sand. Shear strength parameters of Jumunjin Standard Sand were obtained by performing the triaxial compression tests. Direct shear tests were also carried out to find the mechanical properties of the interface between the lining and the covered material. Compared results model tests estimation with respect to displacements of the lining.

• Journal of Industrial Technology
• /
• v.23 no.A
• /
• pp.119-130
• /
• 2003

This thesis is results of experimental works on the behavior of the cut-and-cover tunnel. Centrifuge model tests were performed to simulate the behavior of the cut-and-cover tunnels having cross sections of national road and subway tunnels. Model experiments were carried out with changing the cut slope and the slope of filling ground surface. Displacements of tunnel lining resulted from artificially accelerated gravitational force up to 40g of covered material used in model tests, were measured during centrifuge model tests. In model tests, Jumunjin Standard Sand with the relative density of 80 % and the zinc plates were used for the covered material and the flexible tunnel lining, respectively. Basic soil property tests were performed to obtain it's the property of Jumumjin Standard Sand Shear strength parameters of Jumunjin Standard Sand were obtained by performing the triaxial compression tests. Direct shear tests were also carried out to find the mechanical properties of the interface between the lining and the covered material. Compared results model tests estimation with respect to displacements of the lining.

• Journal of Industrial Technology
• /
• v.18
• /
• pp.7-16
• /
• 1998

This paper is an experimental study of investigating the reinforcing effect and the behavior of cohesionless slope installed with reticulated root pils. Reduced scale model tests with plane strain conditions were performed to study the behavior of the strip footing located on the surface of cohesionless slopes reinforced with root piles. Model tests were carried out with Jumunjin Standard Sand of 45% relative density prepared by raining method to have an uniform slope foundation during tests. Slope of model foundation was 1 : 1.5 and a rigid model slop. Parametric model tests were performed with changing location of model footing, arrangements of root piles and angles of pile installation. On the other hands, the technique with camera shooting was used to monitor sliding surface formed with discontinuty of dyed sand prepared during formation o foudation. From test results, parameters affecting the behavior of model footing were analyzed qualitatively to evaluate their effects on the characteristic of load - settlement, ultimate bearing capacity of model footing and failure mechanism based on the formation of failure surface.

• Geomechanics and Engineering
• /
• v.10 no.1
• /
• pp.109-124
• /
• 2016

In this paper, process of dynamic powder compaction is investigated experimentally using impact of drop hammer and die tube. A series of test is performed using aluminum powder with different grain size. The energy of compaction of powder is determined by measuring height of hammer and the results presented in term of compact density and rupture stress. This paper also presents a mathematical modeling using experimental data and neural network. The purpose of this modeling is to display how the variations of the significant parameters changes with the compact density and rupture stress. The closed-form obtained model shows very good agreement with experimental results and it provides a way of studying and understanding the mechanics of dynamic powder compaction process. In the considered energy level (from 733 to 3580 J), the relative density is varied from 63.89% to 87.41%, 63.93% to 91.52%, 64.15% to 95.11% for powder A, B and C respectively. Also, the maximum rupture stress are obtained for different types of powder and the results shown that the rupture stress increases with increasing energy level and grain size.

• Wind and Structures
• /
• v.26 no.3
• /
• pp.147-161
• /
• 2018

For wind engineering applications downbursts are, presently, almost exclusively modeled, both experimentally and numerically, as transient impinging momentum jets (IJ), even though that model contains none of the physics of real events. As a result, there is no connection between the IJ-simulated downburst wind fields and the conditions of formation of the event. The cooling source (CS) model offers a significant improvement since it incorporates the negative buoyancy forcing and baroclinic vorticity generation that occurs in nature. The present work aims at using large-scale numerical simulation of downburst-producing thunderstorms to develop a simpler model that replicates some of the key physics whilst maintaining the relative simplicity of the IJ model. Using an example of such a simulated event it is found that the non-linear scaling of the velocity field, based on the peak potential temperature (and, hence, density) perturbation forcing immediately beneath the storm cloud, produces results for the radial location of the peak radial outflow wind speeds near the ground, the magnitude of that peak and the time at which the peak occurs that match well (typically within 5%) of those produced from a simple axi-symmetric constant-density dense source simulation. The evolution of the downdraft column within the simulated thunderstorm is significantly more complex than in any axi-symmetric model, with a sequence of downdraft winds that strengthen then weaken within a much longer period (>17 minutes) of consistently downwards winds over almost all heights up to at least 2,500 m.