A level set method, the TVD scheme of second-order upwind procedure coupledwith flux-limiter, and SIMPLE algorithm were incorporated to simulate the flow and interfacialmotion of immiscible two-fluids with large densi...A level set method, the TVD scheme of second-order upwind procedure coupledwith flux-limiter, and SIMPLE algorithm were incorporated to simulate the flow and interfacialmotion of immiscible two-fluids with large density ratio and viscosity ratios, large topologydistortion and surface tension. As a numerical example axi-symmetric rising bubbles wereinvestigated. It is found that the method is numerically stable and has good convergence propertyand the results are in good agreement with other works.展开更多
Presents a method of proof which improves the estimates of entropy production for general total variation diminishing (TVD) schemes. Elements of the general theory of TVD schemes; Basis for obtaining the entropy inequ...Presents a method of proof which improves the estimates of entropy production for general total variation diminishing (TVD) schemes. Elements of the general theory of TVD schemes; Basis for obtaining the entropy inequality of a class of second order resolution-TVD schemes for strict convex conservation laws; Definition of discrete entropy inequality.展开更多
We present a computational gas dynamics method based on the Spectral Deferred Corrections (SDC) time integration technique and the Piecewise Parabolic Method (PPM) finite volume method. The PPM framework is used to de...We present a computational gas dynamics method based on the Spectral Deferred Corrections (SDC) time integration technique and the Piecewise Parabolic Method (PPM) finite volume method. The PPM framework is used to define edge-averaged quantities, which are then used to evaluate numerical flux functions. The SDC technique is used to integrate solution in time. This kind of approach was first taken by Anita et al in [1]. However, [1] is problematic when it is implemented to certain shock problems. Here we propose significant improvements to [1]. The method is fourth order (both in space and time) for smooth flows, and provides highly resolved discontinuous solutions. We tested the method by solving variety of problems. Results indicate that the fourth order of accuracy in both space and time has been achieved when the flow is smooth. Results also demonstrate the shock capturing ability of the method.展开更多
In recent ten years high resolution difference schenies for the computation of thefull unsteady Eulerian system of equations for invisid compressible gas finds celebratedprogress. This paper tests furtherly, by a comp...In recent ten years high resolution difference schenies for the computation of thefull unsteady Eulerian system of equations for invisid compressible gas finds celebratedprogress. This paper tests furtherly, by a complex two-dimensional unsteady problem,four recent schemes. to them attentions are paid. The test problem is the initial stageof a two-dimensional diffraction and reflection of a plane shock wave, impinging on arectangular obstacle. At whose top side there are two sharp corners, near which flow.parameters finds severe variation. There is occurrence of expansion fan with a centerand also concentrated vortices. To simulate them well, the schemes should have goodadaptivity. The special shock Mach number M,=2.068 is so chosen, that at this M,the partical velocity behind impinging shock in fixed coordinate system is just equal tothe speed of sound there, this condition also occurs along a curve in the region ofexpansion fan with a center at the corner. This can clarify the computational featureof different schemes in case,when one of the eigenvalues is just zero. Zero eigenvaluemay spoil some schemes locally. Graphical visualization of the computational resultsmay, show features of the tested schemes about the shock wave resolution, schemeviscosity, expansion wave and the ability. to simulate the process of the generation ofunsteadv concentrated vortex.展开更多
Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though loo...Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though look simple are numerically challenging to study, due to the loss of numerical stability. Over the years, various methodologies have been developed to overcome this numerical limitation. In spite of this, numerical solutions are considered distant from accuracy, as first-order upwind-differencing scheme （UDS） is often employed for improving the stability of algorithm. To elude this effect, some works been reported in the past, where high-resolution-schemes （HRS） were employed and Deborah number was varied. However, these works are limited to creeping flows and do not detail any information on the numerical stability of HRS. Hence, this article presents the numerical study of high shearing contraction flows, where stability of HRS are addressed in reference to fluid elasticity. Results suggest that all I-IRS show some order of undue oscillations in flow variable profiles, measured along vertical lines placed near contraction region in the upstream section of domain, at varied elasticity number E ~ 5. Furthermore, by E, a clear relationship between numerical stability of HRS and E was obtained, which states that the order of undue oscillations in flow variable profiles is directly proportional to E.展开更多
文摘A level set method, the TVD scheme of second-order upwind procedure coupledwith flux-limiter, and SIMPLE algorithm were incorporated to simulate the flow and interfacialmotion of immiscible two-fluids with large density ratio and viscosity ratios, large topologydistortion and surface tension. As a numerical example axi-symmetric rising bubbles wereinvestigated. It is found that the method is numerically stable and has good convergence propertyand the results are in good agreement with other works.
基金The project supported partly by National Natural Science Foundation No.19901031, State Major KeyProject for Basic Research.
文摘Presents a method of proof which improves the estimates of entropy production for general total variation diminishing (TVD) schemes. Elements of the general theory of TVD schemes; Basis for obtaining the entropy inequality of a class of second order resolution-TVD schemes for strict convex conservation laws; Definition of discrete entropy inequality.
文摘We present a computational gas dynamics method based on the Spectral Deferred Corrections (SDC) time integration technique and the Piecewise Parabolic Method (PPM) finite volume method. The PPM framework is used to define edge-averaged quantities, which are then used to evaluate numerical flux functions. The SDC technique is used to integrate solution in time. This kind of approach was first taken by Anita et al in [1]. However, [1] is problematic when it is implemented to certain shock problems. Here we propose significant improvements to [1]. The method is fourth order (both in space and time) for smooth flows, and provides highly resolved discontinuous solutions. We tested the method by solving variety of problems. Results indicate that the fourth order of accuracy in both space and time has been achieved when the flow is smooth. Results also demonstrate the shock capturing ability of the method.
文摘In recent ten years high resolution difference schenies for the computation of thefull unsteady Eulerian system of equations for invisid compressible gas finds celebratedprogress. This paper tests furtherly, by a complex two-dimensional unsteady problem,four recent schemes. to them attentions are paid. The test problem is the initial stageof a two-dimensional diffraction and reflection of a plane shock wave, impinging on arectangular obstacle. At whose top side there are two sharp corners, near which flow.parameters finds severe variation. There is occurrence of expansion fan with a centerand also concentrated vortices. To simulate them well, the schemes should have goodadaptivity. The special shock Mach number M,=2.068 is so chosen, that at this M,the partical velocity behind impinging shock in fixed coordinate system is just equal tothe speed of sound there, this condition also occurs along a curve in the region ofexpansion fan with a center at the corner. This can clarify the computational featureof different schemes in case,when one of the eigenvalues is just zero. Zero eigenvaluemay spoil some schemes locally. Graphical visualization of the computational resultsmay, show features of the tested schemes about the shock wave resolution, schemeviscosity, expansion wave and the ability. to simulate the process of the generation ofunsteadv concentrated vortex.
文摘Viscoelastic fluids due to their non-linear nature play an important role in process and polymer industries. These non-linear characteristics of fluid, influence final outcome of the product. Such processes though look simple are numerically challenging to study, due to the loss of numerical stability. Over the years, various methodologies have been developed to overcome this numerical limitation. In spite of this, numerical solutions are considered distant from accuracy, as first-order upwind-differencing scheme （UDS） is often employed for improving the stability of algorithm. To elude this effect, some works been reported in the past, where high-resolution-schemes （HRS） were employed and Deborah number was varied. However, these works are limited to creeping flows and do not detail any information on the numerical stability of HRS. Hence, this article presents the numerical study of high shearing contraction flows, where stability of HRS are addressed in reference to fluid elasticity. Results suggest that all I-IRS show some order of undue oscillations in flow variable profiles, measured along vertical lines placed near contraction region in the upstream section of domain, at varied elasticity number E ~ 5. Furthermore, by E, a clear relationship between numerical stability of HRS and E was obtained, which states that the order of undue oscillations in flow variable profiles is directly proportional to E.
