• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Anatomy & Biological Anthropology
• /
• v.31 no.4
• /
• pp.121-131
• /
• 2018

Stature is one of the unique biological properties of a person, which can be used for identification of the individual. In this regard, statures are estimated for the unknown victims from crimes and disasters. However, the accuracy of estimates may be compromised by inappropriate methodologies and/or practices of stature estimation. Discussed in this study are the methodological issues related to the current practices of forensic anthropological stature estimation in Korea, followed by suggestions to enhance the accuracy of the stature estimates. Summaries of forensic anthropological examinations for 560 skeletal remains, which were conducted at the National Forensic Service (NFS), were reviewed. Mr. Yoo Byung-eun's case is utilized as an example of the NFS's practices. To estimate Mr. Yoo's stature, Trotter's (1970) femur equation was applied even though the fibula equation of a lower standard error was available. In his case report, the standard error associated with the equation (${\pm}3.8cm$) was interpreted as an 'error range', which gave a hasty impression that the prediction interval is that narrow. Also, stature shrinkage by aging was not considered, so the estimated stature in Mr. Yoo's case report should be regarded as his maximum living stature, rather than his stature-at-death. Lastly, applying Trotter's (1970) White female equations to Korean female remains is likely to underestimate their statures. The anatomical method will enhance the accuracy of stature estimates. However, in cases that the anatomical method is not feasible, the mathematical method based on Korean samples should be considered. Since 1980's, effort has been made to generate stature estimation equations using Korean samples. Applying the equations based on Korean samples to Korean skeletal remains will enhance the accuracy of the stature estimates, which will eventually increase the likelihood of successful identification of the unknowns.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.3
• /
• pp.187-197
• /
• 2010

Carbon dioxide Capture and Storage(CCS) is regarded as one of the most promising options to response climate change. CCS is a three-stage process consisting of the capture of carbon dioxide($CO_2$), the transport of $CO_2$ to a storage location, and the long term isolation of $CO_2$ from the atmosphere for the purpose of carbon emission mitigation. Up to now, process design for this $CO_2$ marine geological storage has been carried out mainly on pure $CO_2$. Unfortunately the $CO_2$ mixture captured from the power plants and steel making plants contains many impurities such as $N_2$, $O_2$, Ar, $H_2O$, $SO_2$, $H_2S$. A small amount of impurities can change the thermodynamic properties and then significantly affect the compression, purification, transport and injection processes. In order to design a reliable $CO_2$ marine geological storage system, it is necessary to analyze the impact of these impurities on the whole CCS process at initial design stage. The purpose of the present paper is to compare and analyse the relevant physical property models including BWRS, PR, PRBM, RKS and SRK equations of state, and NRTL-RK model which are crucial numerical process simulation tools. To evaluate the predictive accuracy of the equation of the state for $CO_2-SO_2$ mixture, we compared numerical calculation results with reference experimental data. In addition, optimum binary parameter to consider the interaction of $CO_2$ and $SO_2$ molecules was suggested based on the mean absolute percent error. In conclusion, we suggest the most reliable physical property model with optimized binary parameter in designing the $CO_2-SO_2$ mixture marine geological storage process.

• Journal of the Korean Society of Food Science and Nutrition
• /
• v.26 no.4
• /
• pp.704-713
• /
• 1997

This study was conducted to describe qualitatively the protein metabolism of pigs during growth depending on the feed protein quality and to describe quantitatively amino acids requirements, using a simulation model. The used model has a non-linear structure. In the used model, the protein utilization system of a pig, which is in the non-steady-state, is described with 15 flux equations and 11 differential equations and is composed with two compartments. Protein deposition(g/day) of pigs on the 30th, 60th, 90th, and 120th day of feeding duration with three-quality protein, beginning with body weight 20kg, were calculated according to the empirical model, PAF(the product of amino acid functions) of Menke, and was used as object function for the simulation. The mean of relative difference between the simulated protein deposition and PAF calculated values, lied in a range of 8.8%. The simulated protein deposition showed different behavior according to feed protein quality. In the high-quality protein, it showed paraboloidal form with extending growth simulation up to 150eh day. So the maximum of protein deposition was acquired on the 105th day of simulate growth time and then it decreased fast. In the low-quality protein, this form of protein deposition in the course of simulated growth did not appear until 150th day. The simulated protein mass also showed a difference in accordance with feed protein quality. The difference was small on the 30th day of simulated growth, but with duration of the simulated growth it was larger. On the 150th day the simulated protein deposition of high quality protein was 1.5 times higher as compared to the low-quality protein. The simulated protein synthesis and break-down rates(g/day) in the whole body showed a parallel behavior in the course of growth, according to feed protein quality. It was found that the improvement of feed protein quality increased protein deposition in the whole body through a increase of both protein synthesis and breakdown during growth. Also protein deposition efficiency, which was calculated from simulated protein deposition and protein synthesis, showed a difference in dependence on the protein qualify of feed protein. The protein deposition efficiency was higher in pigs fed with high quality protein, especially at the simulation time 30th day. But this phenomena disappeared with growth, so on the 150th day of growth, the protein deposition of the high feed protein quality was lowest among the three different quality of feed protein. The simulated total requirement of the 10 essential amino acids for the growth of pigs was 28.1(g/100g protein), similar to NRC. The requirement of lysine was 4.2(g/100g protein).

• Elastomers and Composites
• /
• v.55 no.4
• /
• pp.306-313
• /
• 2020

To investigate the reinforcing effects of functional fillers in nitrile rubber (NBR) materials, high-structure carbon black (HS45), coated calcium carbonate (C-CaCO3), silica (200MP), and multi-walled carbon nanotubes (MWCNTs) were used as functional filler, and carbon black (SRF) as a common filler were used for oil-resistant rubber. The curing and mechanical properties of HS45-, 200MP-, and MWCNT-filled NBR compounds were improved compared to those of the SRF-filled NBR compound. The reinforcing effect also increased with a decrease in the particle size of the fillers. The C-CaCO3-filled NBR compound exhibited no reinforcing effect with increasing filler concentration because of their large primary particle size (2 ㎛). The reinforcing behavior based on 100% modulus of the functional filler based NBR compounds was compared by using several predictive equation models. The reinforcing behavior of the C-CaCO3-filled NBR compound was in accordance with the Smallwood-Einstein equation whereas the 200MP- and MWCNT-filled NBR compounds fitted well with the modified Guth-Gold (m-Guth-Gold) equation. The SRF- and HS45-filled NBR compounds exhibited reinforcing behavior in accordance with the Guth-Gold and m-Guth-Gold equations, respectively, at a low filler content. However, the values of reinforcement parameter (100Mf/100Mu) of the SRF- and HS45-filled NBR compounds were higher than those determined by the predictive equation model at a high filler content. Because the chains of SRF composed of spherical filler particles are similarly changed to rod-like filler particles embedded in a rubber matrix and the reinforcement parameter rapidly increased with a high content of HS45, the higher-structured filler. The reinforcing effectiveness of the functional fillers was numerically evaluated on the basis of the effectiveness index (��SRF/��f) determined by the ratio of the volume fraction of the functional filler (��f) to that of the SRF filler (��SRF) at three unit of reinforcing parameter (100Mf/100Mu). On the basis of their effectiveness index, MWCNT-, 200MP-, and HS45-filled compounds showed higher reinforcing effectiveness of 420%, 70%, and 20% than that of SRF-filled compound, respectively whereas C-CaCO3-filled compound exhibited lower reinforcing effectiveness of -50% than that of SRF-filled compound.

• Geomechanics and Engineering
• /
• v.23 no.2
• /
• pp.151-163
• /
• 2020

Grouting method is an effective way of reinforcing cracked rock masses and plugging water gushing. Current grouting diffusion models are generally developed for horizontal cracks, which is contradictory to the fact that the crack generally occurs in rock masses with irregular spatial distribution characteristics in real underground environments. To solve this problem, this study selected a cement-sodium silicate slurry (C-S slurry) generally used in engineering as a fast-curing grouting material and regarded the C-S slurry as a Bingham fluid with time-varying viscosity for analysis. Based on the theory of fluid mechanics, and by simultaneously considering the deadweight of slurry and characteristics of non-uniform spatial distribution of viscosity of fast-curing grouts, a theoretical model of slurry diffusion in an oblique crack in rock masses at constant grouting rate was established. Moreover, the viscosity and pressure distribution equations in the slurry diffusion zone were deduced, thus quantifying the relationship between grouting pressure, grouting time, and slurry diffusion distance. On this basis, by using a 3-d finite element program in multi-field coupled software Comsol, the numerical simulation results were compared with theoretical calculation values, further verifying the effectiveness of the theoretical model. In addition, through the analysis of two engineering case studies, the theoretical calculations and measured slurry diffusion radius were compared, to evaluate the application effects of the model in engineering practice. Finally, by using the established theoretical model, the influence of cracking in rock masses on the diffusion characteristics of slurry was analysed. The results demonstrate that the inclination angle of the crack in rock masses and azimuth angle of slurry diffusion affect slurry diffusion characteristics. More attention should be paid to the actual grouting process. The results can provide references for determining grouting parameters of fast-curing grouts in engineering practice.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.6
• /
• pp.436-449
• /
• 2013

In this paper, a numerical analysis is carried out to study the characteristics of supercavitating flows and the drag of relatively simple two-dimensional and axisymmetric bodies which can be used for supercavity generation device, cavitator, of a high-speed underwater vehicle. In order to investigate the suitability of numerical models, cavity flows around the hemispherical head form and two-dimensional wedge are calculated with combinations of three turbulence models(standard $k-{\epsilon}$, realizable $k-{\epsilon}$, Reynolds stress) and two cavitation models(Schnerr-Sauer, Zwart-Gerber-Belamri). From the results, it is confirmed that the calculated cavity flow is more affected by the turbulence model than the cavitation model. For the calculation of steady state cavity flows, the convergence in case of the realizable $k-{\epsilon}$ model is better than the other turbulence models. The numerical result of the Schnerr-Sauer cavitation model is changed less by turbulence model and more robust than the Zwart-Gerber-Belamri model. Thus the realizable $k-{\epsilon}$ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate supercavitating flows around disks, two dimensional $10^{\circ}$ and $30^{\circ}$ wedges. In case of the disk, the cavitation number dependences of the cavity size and the drag coefficient predicted are similar to either experimental data or Reichardt's semi-empirical equations, but the drag coefficient is overestimated about 3% higher than the Reichardt's equation. In case of the wedges, the cavitation number dependences of the cavity size are similar to experimental data and Newman's linear theory, and the agreement of the cavity length predicted and Newman's linear theory becomes better as decreasing cavitation number. However, the drag coefficients of wedges agree more with experimental data than those of Newman's analytic solution. The cavitation number dependences of the drag coefficients of both the disk and the wedge appear linear and simple formula for estimating the drag of supercavitating disks and wedges are suggested. Consequently, the CFD scheme of this study can be applied for numerical analysis of supercavitating flows of the cavitator and the cavitator design.

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

In the development of engineering designs, decisions are required irrespective of the state of completeness and quality of information, and are formulated under conditions of uncertainty. Furthermore, under conditions of uncertainty the design invokes risks. Thus, in the design of the structures, the currently used deterministic design method does not provide a realistic assessment of the actual safety or the reliability of the structures. It is desirable that decisions required in The process of the design invariably must be made based on the reliability analysis. Properties of soil material are subject to more uncertainty than those of other structural material. In the field of soil mechanics and foundation engineering, it needed to develop reliability-based design methods. In order to simplify the reliability analysis or the reliability-based design process of the structures associated with the active earth pressure, it is necessary to find the variation and the distribution type of the active earth pressure calculated from the basic properties of soils. Monte Carlo simulation is performed to obtain the relationship between the variation of the active earth pressure for cohessionless soils calculated by using Rankine formula and the basic soil properties and the distribution type of the earth pressure. A series of regression equations obtained by utilizing the multi-linear regression analysis is suggested in this paper and the sensitivity of the basic soil properties to the variation of The earth pressure is investigated. The type of distribution of the active earth pressure was found to be the beta distribution in most cases or to be very similar to the beta distribution, if the basic soil variables are normally distributed.

• Steel and Composite Structures
• /
• v.36 no.5
• /
• pp.553-568
• /
• 2020

Steel-concrete composite structures have been extensively used in building, bridges, and other civil engineering infrastructure. Shear stud connectors between steel and concrete are essential in composite members to guarantee the effectiveness of their behavior in terms of strength and deformability. This study focuses on investigating the shear stiffness of headed studs embedded in several types of concrete with wide range of compressive strength, and their effects on the elastic behavior of steel-concrete composite girders were evaluated. Firstly, totally 206 monotonic push-out tests from the literature were reviewed to investigate the shear stiffness of headed studs embedded in various types of concrete (NC, HPC, UHPC etc.). Shear stiffness of studs is defined as the secant stiffness of the load-slip curve at 0.5V_u, and a formulation for predicting defined shear stiffness in elastic state was proposed, indicating that the stud diameter and the elastic modulus of steel and concrete are the main factors. And the shear stiffness predicted by the new formula agree well with test results for studs with a diameter ranging from 10 to 30 mm in the concrete with compressive strength ranging from 22.0 to 200.0MPa. Then, the effects of shear stiffness on the elastic behaviors of composite girders with different sizes and under different loading conditions were analyzed, the equations for calculating the stress and deformation of simply supported composite girders considering the influence of connection's shear stiffness were derived under different loading conditions using classical linear partial-interaction theory. As the increasing of shear stiffness, the stress and deflection at the most unfavorable section under partial connected condition tend to be those under full connected condition, but the approaching speed decreases gradually. Finally, the connector's shear stiffness was recommended for fully connection in composite girders with different dimensions under different loading conditions. The findings from present study may provide a reference for the prediction of shear stiffness for headed studs and the elastic design of steel-concrete composite girder.

• 제어로봇시스템학회:학술대회논문집
• /
• 1991.10b
• /
• pp.1639-1645
• /
• 1991

Because of the important role LD converters play in the production of high quality steel, various dynamic models have been attempted in the past by many researchers not only to understand the complex chemical reactions that take place in the converter process but also to assist the converter operation itself using computers. And yet no single dynamic model was found to be completely satisfactory because of the complexity involved with the process. The process indeed involves dynamic energy and mass balances at high temperatures accompanied by complex chemical reactions and transport phenomena in the molten state. In the present study, a mathematical model describing the dynamic behavior of LD converter process has been developed. The dynamic model describes the time behavior of the temperature and the concentrations of chemical species in the hot metal bath and slag. The analysis was greatly facilitated by dividing the entire process into three zones according to the physical boundaries and reaction mechanisms. These three zones were hot metal (zone 1), slag (zone 2) and emulsion (zone 3) zones. The removal rate of Si, C, Mn and P and the rate of Fe oxidation in the hot metal bath, and the change of composition in the slag were obtained as functions of time, operating conditions and kinetic parameters. The temperature behavior in the metal bath and the slag was also obtained by considering the heat transfer between the mixing and the slag zones and the heat generated from chemical reactions involving oxygen blowing. To identify the unknown parameters in the equations and simulate the dynamic model, Hooke and Jeeves parttern search and Runge-Kutta integration algorithm were used. By testing and fitting the model with the data obtained from the operation of POSCO #2 steelmaking plant, the dynamic model was able to predict the characteristics of the main components in the LD converter. It was possible to predict the optimum CO gas recovery by computer simulation