• Journal of Mechanical Science and Technology
• /
• v.19 no.4
• /
• pp.1027-1035
• /
• 2005

Thermophoresis in dense gases is studied by using a multi-scale approach and Born- Yvon­Green (BYG) equation. The problem of a particle movement in an ambient dense gas under temperature gradient is divided into inter and outer ones. The pressure gradient in the inner region is obtained from the solutions of BYG equation. The velocity profile is derived from the conservation equations and calculated using the pressure gradient, which provides the particle velocity in the outer problem. It is shown that the temperature gradient applied to the quiescent ambient gas induces some pressure gradient and thus flow tangential to the particle surface in the interfacial region. The mechanism that induces the flow may be the dominant source of the thermophretic particle movement in dense gases. It is also shown that the particle velocity has a nonlinear relationship with the applied temperature gradient and decreases with increasing temperature.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.130-133
• /
• 2007

A two-phase method in CFD has been developed and is applied to model the cavitation flow. The governing equation system is two-phase Navier-Stokes equation, comprised of the mixture mass, momentum and liquid-phase mass equation. It employs an implicite, dual time, preconditioned algorithm using finite difference scheme in curvilineal coordinates and Chien ${\kappa}-{\varepsilon}$ turbulence equation. The experimental cavitating flows around ogive-cylinder and venturi type objects are employed to test the solver. To prove the capabilities of the solver, several three-dimentional examples are presented.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.11 no.2
• /
• pp.157-163
• /
• 2003

This work was undertaken for the numerical analysis of defrosting phenomena of automobile windshield. To analyze the defrost, the flow and temperature field of cabin interior, heat transfer through the windshield glass, and phase change of the frost should be analyzed simultaneously. The flow field was obtained by solving the 3-D unsteady Navier-Stokes equation and the temperature field was computed by energy equation. The phase-change process of Stefan problem was solved by enthalpy method. For code validation, the temperature field of the driven cavity was calculated. The result of calculation shows a good agreement with the other numerical results. Then, the present code was applied to the defrosting analysis of a real automobile and, also, a good agreement with experiment was obtained.

• Journal of Korean Society of Water and Wastewater
• /
• v.21 no.4
• /
• pp.421-428
• /
• 2007

Many engineering problems on the pipeline flow use continuity, energy, friction loss head equation. To calculate friction loss head in a pipeline, Darcy-Weisbach and many average velocity equations can be used and Hazen-Williams equation is used frequently in the pipe network for the water supply systems. Darcy-Weisbach equation is a general one acquired from applying Bernoulli's equation in the pipeline flow and Hazen-Williams equation is a experimental one in case that pipe velocity is below 3m/sec and pipe diameter is over 50mm. In this study, comparing Darcy-Weisbach with Hazen-Williams equation, relation f and C that are expressed as roughness coefficients of those equations is explained. Next, head losses calculated from using those equations are compared and those are applied in realistic pipelines. Comparing f with C, the f is decreasing linearly according to increase of the Reynolds number Re and increasing in case the C is decreasing. additionally, the C is increasing up to a point and then is decreasing according to increase of the Re. Next, the C is increasing and Re's range for increase of the C lengthens in case of decreasing of the relative roughness ${\varepsilon}/d$. Comparing head losses acquired from the two equations, head loss appears large in case that the C is decreasing and the ${\varepsilon}/d$ is increasing. additionally, Head loss calculated by the Darcy-Weisbach equation varies larger than one by Hazen-Williams equation in regard of the Re. Next, change aspect of head loss acquired by the C is distinguished more clearly than the one by the ${\varepsilon}/d$.

• Journal of Energy Engineering
• /
• v.7 no.1
• /
• pp.65-72
• /
• 1998

This paper presents a fast calculation method for evaluating of voltage stability margin (MW) using the line flow equation in polar form. Here, Line flow equations $(P_{ij},\;Q_{ij}$ are comprised of state variable, $V_i,\;{\Delta}_i,\;V_j$ and ${Delta}_j$, and line parameter, r and x. using the feature of polar coordinate, these becomes one equation with two variables, $V_j,;V_j$. Moreover, if bus j is slack or generator bus, which is specified voltage magnitude in load flow calculation, it becomes one equation with one variable $V_ i $ that is, may be formulated with the second-order equation for $V^2_i$. Therefore, multiple load flow solutions may be obtained with simple computation. The obtained load flow multiple solutions are used for evaluating of voltage stability through sensitivity analysis or its closeness. Also, the method is proposed to calculate for voltage stability margin considering shunt capacitor, which is important element for evaluating of voltage stability. The proposed method was validated to sample systems.

• 한국연소학회:학술대회논문집
• /
• 1998.10a
• /
• pp.139-154
• /
• 1998

Under certain circumstance, premixed turbulent flame can be treated as wrinkled thin laminar flame and its motion in a hydrodynamic flow field has been investigated by employing G-equation. Past studies on G-equation successfully described certain aspects of laminar flame propagation such as effects of stretch on flame speed. In those studies, flames were regarded as a passive interface that does not influence the flow field. The experimental evidences, however, indicate that flow field can be significantly modified by the propagation of flames through the volume expansion of burned gas. In the present study, a new method to be used with G -equation is described to include the effect of volume expansion in the flame dynamics. The effect of volume expansion on the flow field is approximated by Biot-Savart law. The newly developed model is validated by comparison with existing analytical solutions of G -equation to predict flames propagating in hydrodynamic flow field without volume expansion. To further investigate the influence of volume expansion, present method was applied to initially wrinkled or planar flame propagating in an imposed velocity field and the average flame speed was evaluated from the ratio of flame surface area and projected area of unburned stream channel. It was observed that the initial wrinkling of flame cannot sustain itself without velocity disturbance and wrinkled structure decays into planar flame as the flame propagates. The rate of decay of the structure increased with volume expansion. The asymptotic change in the average burning speed occurs only with disturbed velocity field. Because volume expansion acts directly on the velocity field, the average burning speed is affected at all time when its effect is included. With relatively small temperature ratio of 3, the average flame speed increased 10%. The combined effect of volume expansion and flame stretch is also considered and the result implied that the effect of stretch is independent of volume release.

• Geomechanics and Engineering
• /
• v.24 no.4
• /
• pp.307-321
• /
• 2021

In an earlier publication (Serrano et al. 2014), the theoretical basis for evaluating the shear strength in rock joints was presented and used to derive an equation that governs the relationship between tangential and normal stresses on the joint during slippage between the joint faces. In this paper, the theoretical equation is applied to two non-linear failure criteria by using non-associated flow laws, including the modified Hoek and Brown and modified Mohr-Coulomb equations. The theoretical model considers the geometric dilatancy, the instantaneous friction angle, and a parameter that considers joint surface roughness as dependent variables. This model uses a similar equation structure to the empirical law that was proposed by Barton in 1973. However, a good correlation with the empirical values and, therefore, Barton's equation is necessary to incorporate a non-associated flow law that governs breakage processes in rock masses and becomes more significant in highly fractured media, which can be induced in a rock joint. A linear law of dilatancy is used to assess the importance of the non-associated flow to obtain very close values for different roughness states, so the best results are obtained for null material dilatancy, which considers significant changes that correspond to soft rock masses or altered zones of weakness.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.59 no.4
• /
• pp.97-107
• /
• 2017

This study is to determine the coefficients of regression equations and to select the optimal regression equation in the LOADEST model after classifying the whole study period into 5 flow conditions for 16 watersheds located in the Nakdonggang waterbody. The optimized coefficients of regression equations were derived using the gradient descent method as a learning method in Tensorflow which is the engine of machine-learning method. In South Korea, the variability of streamflow is relatively high, and rainfall is concentrated in summer that can significantly affect the characteristic analysis of pollutant loads. Thus, unlike the previous application of the LOADEST model (adjusting whole study period), the study period was classified into 5 flow conditions to estimate the optimized coefficients and regression equations in the LOADEST model. As shown in the results, the equation #9 which has 7 coefficients related to flow and seasonal characteristics was selected for each flow condition in the study watersheds. When compared the simulated load (SS) to observed load, the simulation showed a similar pattern to the observation for the high flow condition due to the flow parameters related to precipitation directly. On the other hand, although the simulated load showed a similar pattern to observation in several watersheds, most of study watersheds showed large differences for the low flow conditions. This is because the pollutant load during low flow conditions might be significantly affected by baseflow or point-source pollutant load. Thus, based on the results of this study, it can be found that to estimate the continuous pollutant load properly the regression equations need to be determined with proper coefficients based on various flow conditions in watersheds. Furthermore, the machine-learning method can be useful to estimate the coefficients of regression equations in the LOADEST model.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.1
• /
• pp.26-32
• /
• 2012

The 3-D unsteady viscous flow around an impulsively started rotating marine propeller is simulated using VIC(Vortex-In-Cell) method which is adequate to analyze the strong vortical flow around complicatedly-shaped body. The computational procedure is governed by the vorticity transport equation in Lagrangian form. In order to solve the equation, a regular grid which is independent to the shape of a body is introduced and each term of the equation is evaluated numerically on the grid by applying immersed boundary concept. In this paper, the overall algorithm including the formulation of governing equations and boundary conditions is described and some computational results are presented with discussing their physical validity.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.293-298
• /
• 2005

We calculate the coordinates of an axisymmetric nozzle with a central body. This nozzle ensures a transonic flow with a plane sound surface, which is orthogonal to the symmetry axis and has a wall kink at the sonic point, The Chaplygin transformation in the subsonic part of the flow leads the Dirichlet problem for a system of nonlinear equations. The definition domain of the solution in the velocity-hodograph plane is taken as a rectangle. This enables one to obtain the nozzle with a monotonic distribution of velocity along its subsonic part. In the nonlinear differential equation, the linear Chaplygin operator for plane flows is separated, which allows the iterative calculation of the solution. The supersonic part of the nozzle is calculated under the assumption that the flow at the nozzle exit is uniform and parallel to the symmetry axis; i.e., the supersonic jet outflows to the submerged space with the same pressure. The calculation is performed by the characteristic method. The exact solution of Tricomi equation for near-sonic flows with the straight sonic line is used to 'move away' the sound plane. The velocity distribution alone the supersonic part of the nozzle is also monotonic, which ensures the absence of the boundary-layer separation and, therefore, the adequacy of the ideal-gas model. calculations show that the flow in the supersonic part of the nozzle is continuous (compression shocks are absent)