• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.11 no.1
• /
• pp.95-105
• /
• 2007

A mathematical modeling for the drying process of hygroscopic porous media, such as wood, has been developed in the past decades. The governing equations for wood drying consist of three conservation equations with respect to the three state variables, moisture content, temperature and air density. They are involving simultaneous, highly coupled heat and mass transfer phenomena. In recent, the equations were extended to account for material heterogeneity through the density of the wood and via the density variation of the material process, capillary pressure, absolute permeability, bound water diffusivity and effective thermal conductivity. In this paper, we investigate the drying behavior for the three primary variables of the drying process in terms of control volume finite element method to the heterogeneous transport model on one-dimensional grid.

• Journal of the korean Society of Automotive Engineers
• /
• v.4 no.2
• /
• pp.58-68
• /
• 1982

It is well known to diesel engineers that the heat release pattern is one of the most important factors affecting engine performance. Thorough research in heat release pattern has materially helped the progress in high-speed diesel engine development . This paper is based on the research conducted at KAIST and Daewoo Heavy Industry last year. The purpose of this paper is to determine the heat release pattern in combustion chamber of MAN M type, the famous low-noise engine. Thermodynamic cycle simulation was performed using Whitehous-Way's heat release pattern with modified coefficients and Annand's heat transfer model. Instantaneous temperature and pressure of gas in cylinder could be determined by the numerical solution of simultaneous equation of mass conservation, equation of energy conservation, and state equation of ideal gas. Calculated results were compared with measured values in some details emphasizing upon the factors affecting rate of heat release. The agreement was fairly good and revealed why M type should have lower burning velocity at the early part of combustion in spite of high injection rate. Additional results by parametric studies were given in relation to fuel injection conditions for further application to engine development.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2003.09a
• /
• pp.497-504
• /
• 2003

An optimal design method for hybrid structural control system of building structures subject to earthquake excitation is presented in this paper. Designing a hybrid structural control system nay be defined as a process that optimizes the capacities and configuration of passive and active control systems as well as structural members. The optimal design proceeds by formulating the optimization problem via a multi-stage goal programming technique and, then, by finding reasonable solution to the optimization problem by means of a goal-updating genetic algorithm. The process of the integrated optimization design is illustrated by a numerical simulation of a nine-story building structure subject to earthquake excitation. The effectiveness of the proposed method is demonstrated by comparing the optimally designed results with those of a hybrid structural control system where structural members, passive and active control systems are uniformly distributed.

• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.829-862
• /
• 2004

This paper presents a method for the nonlinear analysis of beam elements subjected to the cyclical combined actions of torsion, biaxial flexure and axial forces based on an extension of the disturbed compression field (DSFM). The theoretical model is based on a hybrid formulation between the full rotation of the cracks model and the fixed direction of the cracking model. The described formulation, which treats cracked concrete as an orthotropic material, includes a new approach for the evaluation of the re-orientation of both the compression field and the deformation field by removing the restriction of their coincidence. A new equation of congruence permits evaluating the deformation of the middle line. The problem consists in the solution of coupled nonlinear simultaneous equations expressing equilibrium, congruence and the constitutive laws. The proposed method makes it possible to determine the deformations of the beam element according to the external stresses applied.

• 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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.1
• /
• pp.13-22
• /
• 2013

This paper deals with a novel method for numerical analyses of the tapered geometrical non-linear beam with three unknown parameters, subjected a floating point load. The beams with hinged-movable end constraint are chosen as the objective beam. Cross sections of the beam whose flexural rigidities are functionally varied with the axial coordinate. The first order simultaneous differential equations governing the elastica of such beam are derived on the basis of the Bernoulli-Euler beam theory. A novel numerical method for solving these equations is developed by using the iteration technique. The processes of the solution method are extensively discussed through a typical numerical example. For validating theories developed herein, laboratory scaled experiments are conducted.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.1
• /
• pp.34-44
• /
• 1992

A numerical model using Hybrid Element Method(HEM) is presented for the prediction of wave agitations in a harbor which are induced by the intrusion and transformation of incident short-period waves. A linear mild-slope equation including bottom friction is used as the governing equation and a partial absorbing boundary condition is used on solid boundaries. Functional derived in the present paper is based on the Chen and Mei(1974)'s concept which uses finite element net in the inner region and analytical solution of Helmholtz equation in the outer region. Final simultaneous equations are solved using the Gaussian Elimination Method. The model appears to be reasonably good from the comparison of numerical calculation with hydraulic experimental results of short-wave diffraction through a breakwater gap(Pos and Kilner, 1987). The problem of requring large computational memory could be overcome using 8-noded isoparametric elements.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.185-192
• /
• 1993

A numerical model has been used for the prediction of wave agitations in a harbor which are induced by the intrusion and transformation of incident waves. Based on linear wave theory a mild-slope equation has been used. A partial absorbing boundary condition has been used on solid boundary. Functional has been derived following Chen and Mei(l974)'s technique based on Hybrid Element Method which uses finite discretisation in the inner region and analytical solution of Helmholtz equation in the outer region. Final simultaneous equation has been solved using the Gaussian Elimination Method. Helmholtz natural period and second peak period of seiche in Donghae Harbor coincide very well with the results from numerical calculation. Computed amplification factors show good agreement, especially when the reflection coefficient on solid boundary is 0.99, with those of measurements.

• Journal of Korean Society of Industrial and Systems Engineering
• /
• v.38 no.4
• /
• pp.80-87
• /
• 2015

The proliferation of the Internet and communication technologies and applications, besides the conventional retailers, has led to a new form of distribution channel, namely home sopping through the telephone, TV, catalog or the Internet. The conventional and new distribution channels have different transaction costs perceived by the consumers in the following perspectives: the accessibility to the product information, the traffic cost and the opportunity cost for the time to visit the store, the possibility of 'touch and feel' to test the quality of the product, the delivery time and the concern for the security for the personal information. Difference in the transaction costs between the distribution channels results in the different selling prices even for the same product. Moreover, distribution channels with different selling prices necessarily result in different business surpluses. In this paper, we study the multilateral bargaining strategy of a manufacturer who sells a product through multiple distribution channels with different transaction costs. We first derive the Nash equilibrium solutions for both simultaneous and sequential bargaining games. The numerical analyses for the Nash equilibrium solutions show that the optimal bargaining strategy of the manufacturer heavily depends not only on the degree of competition between the distribution channels but on the difference of the business surpluses of the distribution channels. First, it is shown that there can be four types of locally optimal bargaining strategies if we assume the market powers of the manufacturer over the distribution channels can be different. It is also shown that, among the four local optimal bargaining strategies, simultaneous bargaining with the distribution channels is the most preferred bargaining strategy for the manufacturer.

• Clean Technology
• /
• v.23 no.1
• /
• pp.42-53
• /
• 2017

Numerical analysis was done to evaluate the fluid distribution inside of the mixing quencher to increase the reaction efficiency of the aqueous hydrogen peroxide solution in the scrubbing column which is used for simultaneous desulfurization and denitrification. Effective injection of the aqueous hydrogen peroxide ($H_2O_2$) solution in the mixing quencher has major effects for improving the reaction efficiency in the scrubbing column by enhancing the mixing of the aqueous $H_2O_2$ solution with the exhaust gas. The current study is to optimize the array of nozzles and the spray angles of the aqueous $H_2O_2$ solution in the mixing quencher by using the computational method. Main concerns of the analysis are how to enhance the uniformity of the $H_2O_2$ concentration distribution in the internal flow. Numerical analysis was done to check the distribution of the internal flow in the mixing quencher in terms of RMS values of the $H_2O_2$ concentration at the end of quencher. The concentration distribution of $H_2O_2$ at the end of is evaluated with respect to the different array of the nozzle pipes and the nozzle tip angles, and we also analyzed the turbulence formation and fluid mixing in the zone. The effect of the spray angle was evaluated with respect to the mixing efficiency in different flow directions. The optimized mixing quencher had the nozzle array at location of 0.3 m from the inlet duct surface and the spray angle is $15^{\circ}$ with the co-current flow. The RMS value of the $H_2O_2$ concentration at the end of the mixing quencher was 12.4%.