• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.8-17
• /
• 2004

A Semi-Lagrangian method based on CIP(Cubic Interpolated Pseudoparticle)method is proposed and it is applied to solve the two dimensional advection equation. Especially the attentions are given to settle the pole problem and to enhance the accuracy in solving the advection equation on the spherical coordinate system. Tn this algorithm, the CU method is employed as the Semi-Lagrangian method and extended to the spherical coordinate system. To enhance the accuracy of the solution, the spatial discretization is made by CIP method. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional rotating cosine bell problem and the frontogenesis problem are simulated by the present scheme. As results, it is confirmed that the present scheme gives an accurate solution and settles the pole problem in the advection equation on the sphere.

• Journal of Environmental Science International
• /
• v.11 no.9
• /
• pp.907-914
• /
• 2002

Various Eulerian-Lagrangian models for the one-dimensional longitudinal dispersion equation in nonuniform flow were studied comparatively. In the models studied, the transport equation was decoupled into two component parts by the operator-splitting approach; one part is governing advection and the other is governing dispersion. The advection equation has been solved by using the method of characteristics following fluid particles along the characteristic line and the results were interpolated onto an Eulerian grid on which the dispersion equation was solved by Crank-Nicholson type finite difference method. In the solution of the advection equation, Lagrange fifth, cubic spline, Hermite third and fifth interpolating polynomials were tested by numerical experiment and theoretical error analysis. Among these, Hermite interpolating polynomials are generally superior to Lagrange and cubic spline interpolating polynomials in reducing both dissipation and dispersion errors.

• Water for future
• /
• v.27 no.2
• /
• pp.155-166
• /
• 1994

Various Eulerian-Lagrangian numerical models for the one-dimensional longitudinal dispersion equation are studied comparatively. In the model studied, the transport equation is decoupled into two component parts by the operator-splitting approach ; one part governing adveciton and the other dispersion. The advection equation has been solved using the method of characteristics following fluid particles along the characteristic line and the results are interpolated onto an Eulerian grid on which the dispersion equation is solved by Crank-Nicholson type finite difference method. In solving the advection equation, various interpolation schemes are tested. Among those, Hermite interpolation polynomials are superior to Lagrange interpolation polynomials in reducing dissipation and dispersion errors in the simulation.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1993.07a
• /
• pp.151-157
• /
• 1993

• Journal of applied mathematics & informatics
• /
• v.18 no.1_2
• /
• pp.329-337
• /
• 2005

This paper is to establish a functional equation satisfied by the generating function for counting rooted cubic c-nets and then to determine the parametric expressions of the equation directly. Meanwhile, the explicit formulae for counting rooted cubic c-nets are derived immediately by employing Lagrangian inversion with one or two parameters. Both of them are summation-free and in which one is just an answer to the open problem (8.6.5) in [1].

• Journal of Mechanical Science and Technology
• /
• v.18 no.8
• /
• pp.1470-1478
• /
• 2004

In this paper, a new 'presumed' Probability Density Function (PDF) approach coupled with a Lagrangian tracking method is proposed for turbulent combustion modeling. The test and the investigation of the model are conducted by comparing the model results with DNS data for a premixed flame subjected in a decaying turbulent field. The newly constructed PDF, which incorporates the instantaneous chemical reaction term, demonstrates consistent improvement over conventional assumed PDF models. It has been found that the time evolution of the mean scalar, the variance and the mean reaction rate are strongly influenced by a parameter deduced by a Lagrangian equation which takes into account explicitly the local reaction rate. Tests have been performed for a moderate Damkohler number, and it is expected the model may cover a broader range of Damkohler number. The comparison with the DNS data demonstrates that the proposed model may be promising and affordable for implementation in a moment-equation solver.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.5
• /
• pp.41-48
• /
• 2002

Numerical experiments are conducted using a three-dimensional baroclinic equation model and a Lagrangian method for clarifying the effect of th settling velocity on the suspended solids diffusion caused by the dredging and the reclamination works. Diffusion characteristics of the neutral particles and the weighting particles is experimented by the Lagrangian particles trajectory model, The results show that the diffusion characteristics of the suspended solids is effected by the settling velocity classified by the particles size in the density layered semi-closed bay. To estimate exactly the diffusion characteristics of the suspended solids and the contaminant with weight the three-dimensional baroclinic equation model and the three-dimensional Lagrangian particles trajectory model considering the settling velocity of the particle in the density layered semi-closed bay must be used.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.43-48
• /
• 2003

When the Eulerian-Lagrangian method is used to analyze the particle laden two-phase flow, a large number of particles should be used to obtain statistically meaningful solutions. Then it takes too much time to track the particles and to average the particle properties in the numerical analysis of two-phase flow. The purpose of this paper is to reduce the computation time by means of a set of particle gird separate to the flow grid. Particle motion equation here is the simplified B-B-O equation, which is integrated to get the particle trajectories. Particle turbulent dispersion, wall collision, and wall roughness effects are considered but the two-way coupling effects between gas and particles are neglected. Particle laden 2-D channel flow is solved and it is shown that the computational efficiency is indeed improved by using the current method

• Water for future
• /
• v.28 no.4
• /
• pp.125-136
• /
• 1995

ELM, a characteristic line based method, was applied to advection-dispersion equation, and the results obtained were compared with those of Eulerian schemes(Stone-Brian and QUICKEST). The calculation methods consisted of Lagrangian interpolation scheme and cubic spline interpolation scheme for the advection calculation, and the Crank-Nicholson scheme for the dispersion calculation. The results of numerical methods were as follows: (1) for Gaussian hill: ELM, using Lagrangian interpolation scheme, gave the most accurate computational result, ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical damping for Peclet number 50. Stone-Brian scheme gave phase shift introduced in the numerical solution for Peclet number 10 and 50. (2) for advanced front: All schemes gave accurate computational results for Peclet number 1 and 4. ELM, Lagrangian interpolation scheme, and Stone,Brian scheme gave dissipation error and ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical oscillation for Peclet number 50.

• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4395-4407
• /
• 2023

Bubble columns are widely encountered in several industries, especially in the field of nuclear safety. The Eulerian-Eulerian and the Eulerian-Lagrangian methods are commonly used to investigate bubble columns. Eulerian approaches require additional tasks such as strict volume conservation at the interface and a predefined well-structured grid. In contrast, the Lagrangian approach can be easily implemented. Hence, we introduce a fully Lagrangian approach for the simulation of bubble columns using the discrete bubble model (DBM) and moving particle semi-implicit (MPS) methods. Additionally, we propose a rigorous method to estimate the volume fraction accurately, and verified it through experimental data and analytical results. The MPS method was compared with the experimental data of Dambreak. The DBM was verified by analyzing the terminal velocity of a single bubble for each bubble size. It agreed with the analytical results for each of the four drag correlations. Additionally, the improved method for calculating the volume fraction showed agreement with the Ergun equation for the pressure drop in a packed bed. The implemented MPS-DBM was used to simulate the bubble column, and the results were compared with the experimental results. We demonstrated that the MPS-DBM was in quantitative agreement with the experimental data.