• Journal of Korean Society on Water Environment
• /
• v.23 no.5
• /
• pp.683-688
• /
• 2007

Water flow into unsaturated soils is most often modeled by Richards equation consisting of the mass conservation law and Darcy's law. Three standard forms of Richards equation are presented as the head (${\Psi}$)-based form, the moisture content (${\theta}$) based form, and the mixed form. Numerical solutions of these partial differential equations with highly nonlinear terms can cause poor results along with significant mass balance errors. The numerical solution based on the mixed form of Richards equation is known that the mass is perfectly conserved without any additional computational efforts. The aim of this study is to develop fully implicit numerical scheme of Richards equation for one-dimensional vertical unsaturated flow in homogeneous soils using the finite difference approximation, and then to perform sensitivity analysis of infiltration to the variations in the unsaturated soil properties and to different soil types.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.779-787
• /
• 2007

This paper presents a highly efficient moving least squares finite difference method (MLS FDM) for a heat transfer problem of bi-material with interfacial boundary. The MLS FDM directly discretizes governing differential equations based on a node set without a grid structure. In the method, difference equations are constructed by the Taylor polynomial expanded by moving least squares method. The wedge function is designed on the concept of hyperplane function and is embedded in the derivative approximation formula on the moving least squares sense. Thus interfacial singular behavior like normal derivative jump is naturally modeled and the merit of MLS FDM in fast derivative computation is assured. Numerical experiments for heat transfer problem of bi-material with different heat conductivities show that the developed method achieves high efficiency as well as good accuracy in interface problems.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.106-112
• /
• 2004

A kinetic theory analysis is made of low-speed gas flows around a micro-plate. The Boltzmann equation simplified by a collision model is solved by means of a finite difference approximation with the discrete ordinate method. The method does not suffer from statistical noise which is common in particle based methods and requires much less amount of computational effort. Calculations are made for flows around a micro-scale flat plate with a finite length of 20 microns. The method is assessed by comparing the results with those from several different methods and available experimental data.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.1-7
• /
• 2004

A kinetic theory analysis is made of low-speed rarefied gas flows around a flat plate. The Boltzmann equation simplified by a collision model is solved by means of a finite difference approximation with the discrete ordinate method. The method does not suffer from statistical noise which is common in particle based methods and requires much less amount of computational effort. Calculations are made for flows around a micro-scale flat plate with a finite length of 20 microns. The method is assessed by comparing the results with those from several different methods and available experimental data.

• Bulletin of the Korean Mathematical Society
• /
• v.58 no.5
• /
• pp.1109-1127
• /
• 2021

We study the long-term dynamics for a family of incompressible three-dimensional Leray-α-like models that employ the spectral fractional Laplacian operators. This family of equations interpolates between incompressible hyperviscous Navier-Stokes equations and the Leray-α model when varying two nonnegative parameters 𝜃₁ and 𝜃₂. We prove the existence of a finite-dimensional global attractor for the continuous semigroup associated to these models. We also show that an operator which projects the weak solution of Leray-α-like models into a finite-dimensional space is determining if it annihilates the difference of two "nearby" weak solutions asymptotically, and if it satisfies an approximation inequality.

• Korea-Australia Rheology Journal
• /
• v.21 no.1
• /
• pp.47-58
• /
• 2009

We have investigated the transient behavior of 1D fully developed Poiseuille viscoelastic flow under finite pressure gradient described by the Oldroyd-B and Leonov constitutive equations. For analysis we employ a simple $2^{nd}$ order discretization scheme such as central difference for space and the Crank-Nicolson for time approximation. For the analysis of the Oldroyd-B model, we also apply the analytical solution, which is obtained again in this work in terms of elementary solution procedure simpler than the previous one (Waters and King, 1970). Both models demonstrate qualitatively similar solutions, but their eventual steady flowrate exhibits noticeable difference due to the absence or presence of shear thinning behavior. In the inertialess flow, the flowrate instantaneously attains a large value corresponding to the Newtonian creeping flow and then decreases to its steady value when the applied pressure gradient is low. However with finite liquid density the flow field shows severe fluctuation even accompanying reversals of flow directions. As the assigned pressure gradient increases, the flowrate achieves its steady value significantly higher than its value during oscillations after quite long period of time. We have also illustrated comparison between 1D and 2D results and possible mechanism of complex 2D flow rearrangement employing a previous solution of [mite element computation. In addition, we discuss some mathematical points regarding missing boundary conditions in 2D modeling due to the change of the type of differential equations when varying from inertialess to inertial flow.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.491-496
• /
• 2003

In sensitivity analysis, semi-analytical method(SAM) reveals severe inaccuracy problem when relatively large rigid body motions are identified for individual elements. Recently such errors of SAM resulted by the finite difference scheme have been improved by the separation of rigid body mode. But the eigenvalue should be obtained first before the sensitivity analysis is performed and it takes much time in the case that large system is considered. In the present study, by constructing a reduced one from the original system, iterative method combined with mode decomposition technique is proposed to compute reliable semi-analytical design sensitivities. The sensitivity analysis is performed by the eigenvector acquired from the reduced system. The error of SAM caused by difference scheme is alleviated by Von Neumann series approximation.

• Korea-Australia Rheology Journal
• /
• v.16 no.4
• /
• pp.227-233
• /
• 2004

Simplified model of slide-fed curtain coating flow has been developed and tested in this study. It rests on the sheet profile equations for curtain thickness in curtain flow and its trajectory derived by the integral momentum balance approach of Higgins and Scriven (1979) and Kistler (1983). It also draws on the film profile equation of film thickness variation in flow down a slide. The equations have been solved in finite difference approximation by Newton iteration with continuation. The results show that how inertia (Rey­nolds number), surface tension (capillary number), inclination angle of the slide, and air pressure difference across the curtain affect sheet trajectory and thickness profile. It has been revealed that approximate models can be useful to easily analyze coating flow dynamics without complex computations, giving qualitative agreement with full theory and with experiment.

• Proceedings of the KSME Conference
• /
• 2001.06a
• /
• pp.686-691
• /
• 2001

Structural optimization often require the evaluation of design sensitivities. The Semi Analytic method(SAM) is popular for shape optimization because this method has several advantages. But when relatively large rigid body motions are identified for individual elements, the SA method shows severe inaccuracy. In this paper, the improvement of design sensitivities corresponding to the rigid body mode is evaluated by exact differentiation of the rigid body modes. Moreover, the error of the SA method caused by numerical difference scheme is alleviated by using a series approximation for the sensitivity derivatives and considering the higher order terms.