• Journal of Korean Society of Forest Science
• /
• v.88 no.3
• /
• pp.320-326
• /
• 1999

In this study, 15 survival functions in integral and difference forms for forest plantation were derived based on assumptions for the number of surviving trees and the differential forms of the mortality rate model. Then, performance of the models was evaluated by fitting to remeasurement data of unthinned white pine(Pinus strobes) forest plantation. As a result, three equations associated with a power function of age, $t^{\beta}$, are somewhat more suitable for describing the effect of self-thinning over time. On the other hand, a general survival function for Japanese larch(Larix leptolepts) forest plantation was derived in order to exam the effect of site quality on self-thinning procedures. The results indicate that the $N_{min}$ is negatively correlated with site index and, even though the same initial stand density was assumed, the survival function curves differ in shapes associated with site index values.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.1
• /
• pp.53-60
• /
• 2005

The use of pipe-bends brings about non-uniform flows at the exit of them due to the velocity difference between inner and outer flows inside the bend. These phenomena may cause turbopump of satellite launch vehicle to run off-design and reduce its efficiency, and also introduce unstable influx of propellants to engine manifold after passing through a turbopump. In order to improve the uniformity of flow at the bend exit, certain turning vanes are set up in the bend pipe normally. Correspondingly the design is an $90^{\circ}\;and\;45^{\circ}$ bend pipes that incorporate with the maximum three turning vanes. All designs were analyzed with numerical analysis by solving the Navier-Stokes equations in three dimensions in case of each respective fuel and oxidizer. Evaluations of the vaned pipe bends designs were accomplished by the velocity magnitude distributions and the predicted pressure drops. We could find that the more vaned bend pipe and larger angle pipe under consideration effectively, the more uniform velocity magnitude of the bend and pressure losses.

• Korean Journal of Community Nutrition
• /
• v.8 no.6
• /
• pp.993-1000
• /
• 2003

The purpose of this study was to provide baseline data for revising the recommended energy intake for Korean adults. We recruited 290 adults so as to determine their resting energy expenditure (REE) and energy intake. The REE was measured by indirect calorimetry. We also calculated the REE from prediction equations formulated by World Health Organization (WHO), The energy intake for two consecutive days was assessed using the 24 hour recall method. The body weight, lean body mass (LBM) and percentage body fat were measured using the INBODY 3.0 system. We compared the results of three age groups ; 20 to 29 years,30 to 49 years and 50 years or more. The average energy intake of each age group was below the 7th Korean Recommended Dietary Allowances (RDA). The average energy intake was lower in the older age groups. However, no difference was observed among the age groups when the energy intake as a percentage of the Korean RDA was compared. Our measurement of the REE was higher than the REE calculated by the WHO's method. Correlation coefficients between the measured and the calculated values of REE for all age groups showed significant correlations (r=0.475-0.672). As the ages of all the subjects increased, the REE/kg of body weight decreased. There were no significant differences in the REE / kg of the LBM between the different age groups; however, the REE/kg of the LBM was higher in the female group than in the male group. Negative correlations of the REE with the age (r=-0.242) and body fat ratio (r=-0.313) were observed; positive correlations of the REE with the BMI (r=0.265), height (r=0.570), weight (r=0.562) and LBM (r=0.586) were also found (p<0.01).

• The Korean Journal of Ecology
• /
• v.18 no.1
• /
• pp.189-201
• /
• 1995

In order to clarify regeneration processes and mechanisms of the disturbed Pinus densiflora forest, responses of Pinus densiflora to gap formed by disturbance were analysed by growth of saplings and mature and growth equations were obtained from branch growth of mature trees and height growth of saplings, and age distribution of saplings and young trees recruited within gap was analysed in relation to gap age. Height growth of saplings within gaps was accelerated after gap formation. Such abrupt increases of growth of saplings after the gap formation might be resulted in the difference of growth of saplings between gap and non-gap areas. In fact, height and diameter of saplings in the central part of gap were larger than those of saplings in marginal parts of gap and non-gap area. However, density of saplings was not different in both parts. In addition, growth of annual rings of mature trees bordering on gap also increased after gap formation. Branch growth of mature trees bodering on gap was 6.3 - 6.5 cm /year and the mean radius of gaps created by death of only one canopy tree was about 3 m. Therefore, for those gaps to be closed by branch growth it will take 46 years. Growth of saplings within gap showed exponential equation. Fifty years will be required for the saplings to enter the forest canopy by the exponential growth equation. Therefore, gap created by only one tree might be closed by branch growth of surrounding canopy trees in advance of being done by height growth of saplings. But gaps created by death of trees more than 2 will be closed by the growth of saplings. Among the regenerating saplings and young trees within gaps, individuals established in advance of gap formation were more than those established after the gap formation. From these results, it was assumed that the disturbed Pinus densiflora forests in these sites were regenerated by height growth of saplings recruited in advance of gap formation.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.60 no.1
• /
• pp.142-149
• /
• 2011

This paper presents a new state of charge(SOC) estimation of large capacity of Li-ion battery (LIB) based on the multiple model adaptive estimation(MMAE) method. We first introduce an equivalent circuit model of LIB. The relationship between the terminal voltage and the open circuit voltage(OCV) is nonlinear and may vary depending on the changes of temperature and C-rate. In this paper, such behaviors are described as a set of multiple linear time invariant impedance models. Each model is identified at a temperature and a C-rate. These model set must be obtained a priori for a given LIB. It is shown that most of impedances can be modeled by first-order and second-order transfer functions. For the real time estimation, we transform the continuous time models into difference equations. Subsequently, we construct the model banks in the manner that each bank consists of four adjacent models. When an operating point of cell temperature and current is given, the corresponding model bank is directly determined so that it is included in the interval generated by four operating points of the model bank. The MMAE of SOC at an arbitrary operating point (T $^{\circ}C$, $I_{bat}$[A]) is performed by calculating a linear combination of voltage drops, which are obtained by four models of the selected model bank. The demonstration of the proposed method is shown through simulations using DUALFOIL.

• Journal of applied mathematics & informatics
• /
• v.11 no.1_2
• /
• pp.305-316
• /
• 2003

Recently, an iterative algorithm for finding the interior eigenvalues of a definite matrix by CG-type method has been proposed. This method compares to the inverse power method. The given matrices A, and B are assumed to be large and sparse, and SPD( Symmetric Positive Definite) The CG scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for large sparse eigenproblems for smallest eigenvalue. Also, it is very amenable to parallel computations, like the CG method for the linear systems. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. But for parallel computations we need to find an efficient parallel preconditioner. Our candidates we ILU(0) in the wave-front order, ILU(0) in the multi-coloring order, Point-SSOR(Symmetric Successive Overrelaxation), and Multi-Color Block SSOR preconditioner. Wavefront order is a simple way to increase parallelism in the natural order, and Multi-coloring realizes a parallelism of order(N), where N is the order of the matrix. Another choice is the Multi-Color Block SSOR(Symmetric Successive OverRelaxation) preconditioning. Block SSOR is a symmetric preconditioner which is expected to minimize the interprocessor communication due to the blocking. We implemented the results on the CRAY-T3E with 128 nodes. The MPI (Message Passing Interface) library was adopted for the interprocessor communications. The test problem was drawn from the discretizations of partial differential equations by finite difference methods. The results show that for small number of processors Multi-Color ILU(0) has the best performance, while for large number of processors Multi-Color Block SSOR performs the best.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.1
• /
• pp.1-7
• /
• 2014

The present study numerically investigated the effect of the advance ratio on the wake characteristics of the marine propeller in the propeller open water test. Therefore, a wide range of the advance ratio(0.2${\kappa}-{\omega}$SST Model are considered. The three-dimensional vortical structures of tip vortices are visualized by the swirl strength, resulting in fast decay of the tip vortices with increasing the advance ratio. Furthermore, to better understanding of the wake evolution, the contraction ratio of the slip stream for different advance ratios is extracted from the velocity fields. Consequently, the slip stream contraction ratio decreases with increasing the advance ratio and successively the difference of the slip stream contraction ratio between J=0.2 and J=0.8 is about 0.1R.

• ETRI Journal
• /
• v.21 no.1
• /
• pp.1-27
• /
• 1999

We present a new description of envelope-function equation of the superlattice (SL). The SL wave function and corresponding effective-mass equation are formulated in terms of a linear combination of Bloch states of the constituent material with smaller band gap. In this envelope-function formalism, we review the fundamental concept on the motion of a wave packet in the SL structure subjected to steady and uniform electric fields F. The review confirms that the average of SL crystal momentums K = ($k_x,k_y,q$), where ($K_x,k_y$) are bulk inplane wave vectors and q SL wave vector, included in a wave packet satisfies the equation of motion = $_0+Ft/h$; and that the velocity and acceleration theorems provide the same type of group velocity and definition of the effective mass tensor, respectively, as in the Bulk. Finally, Schlosser and Marcus's method for the band theory of metals has been by Altarelli to include the interface-matching condition in the variational calculation for the SL structure in the multi-band envelope-function approximation. We re-examine this procedure more thoroughly and present variational equations in both general and reduced forms for SLs, which agrees in form with the proposed envelope-function formalism. As an illustration of the application of the present work and also for a brief investigation of effects of band-parameter difference on the subband energy structure, we calculate by the proposed variational method energies of non-strained $GaAs/Al_{0.32}Ga_{0.68}As$ and strained $In_{0.63}Ga_{0.37}As/In_{0.73}Ga_{0.27}As_{0.58}P_{0.42}SLs$ with well/barrier widths of $60{\AA}/500{\AA}$ and 30${\AA}/30{\AA}$, respectively.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.

• Structural Engineering and Mechanics
• /
• v.12 no.6
• /
• pp.657-668
• /
• 2001

Fiber reinforced cementitious composites are nowadays widely applied in civil engineering. The postcracking performance of this material depends on the interaction between a steel fiber, which is obliquely across a crack, and its surrounding matrix. While the partly debonded steel fiber is subjected to pulling out from the matrix and simultaneously subjected to transverse force, it may be modelled as a Bernoulli-Euler beam partly supported on an elastic foundation with non-linearly varying modulus. The fiber bridging the crack may be cut into two parts to simplify the problem (Leung and Li 1992). To obtain the transverse displacement at the cut end of the fiber (Fig. 1), it is convenient to directly solve the corresponding differential equation. At the first glance, it is a classical beam on foundation problem. However, the differential equation is not analytically solvable due to the non-linear distribution of the foundation stiffness. Moreover, since the second order deformation effect is included, the boundary conditions become complex and hence conventional numerical tools such as the spline or difference methods may not be sufficient. In this study, moment equilibrium is the basis for formulation of the fundamental differential equation for the beam (Timoshenko 1956). For the cantilever part of the beam, direct integration is performed. For the non-linearly supported part, a transformation is carried out to reduce the higher order differential equation into one order simultaneous equations. The Runge-Kutta technique is employed for the solution within the boundary domain. Finally, multi-dimensional optimization approaches are carefully tested and applied to find the boundary values that are of interest. The numerical solution procedure is demonstrated to be stable and convergent.