Flux Limiter Lattice Boltzmann for Compressible Flows
|
|
- Leonard Charles
- 5 years ago
- Views:
Transcription
1 Commun. Theor. Phys. 56 ( Vol. 56, No. 2, August 15, 2011 Flux Limiter Lattice Boltzmann for Compressible Flows CHEN Feng (íô, 1 XU Ai-Guo (ÆÇÁ, 2, ZHANG Guang-Cai (¾, 2 and LI Ying-Jun (Ó 1, 1 State Key Laboratory for GeoMechanics and Deep Underground Engineering, China University of Mining and Technology (Beijing, Beijing , China 2 National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, P.O. Box , Beijing , China (Received December 27, 2010; revised manuscript received May 11, 2011 Abstract In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann (LB model for shacked compressible flows. The proposed flux limiter scheme is efficient in decreasing the artificial oscillations and numerical diffusion around the interface. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. Numerical simulations for the Richtmyer Meshkov instability show that with the new model the computed interfaces are smoother and more consistent with physical analysis. The growth rates of bubble and spike present a satisfying agreement with the theoretical predictions and other numerical simulations. PACS numbers: j, y, Dd Key words: lattice Boltzmann method, flux limiter technique, compressible flows, multiple-relaxation-time, Richtmyer Meshkov instability 1 Introduction In the past few decades, the lattice Boltzmann (LB method [1 2] has become a powerful tool for computational fluids. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. The starting point of LB method is the evolution equation for the distribution function f i, f i t + v f i iα = Ω i, (1 x α where v iα is the discrete particle velocity, Ω i is the collision operator. The equation (1 may be solved with an appropriate finite difference scheme defined on the lattice. Generally, the forward Euler method is used to compute the time derivative f i / t. Appropriate spatial difference scheme is essential to reduce spurious behaviour around the interface, and the choice requires more consideration. The upwind difference scheme leads to stable simulation, but it gives very smeared solution with excessive numerical diffusion; [3 4] the higher order schemes, such as Lax Wendroff or Beam Warming scheme, are more accurate, but the wiggle phenomenon introduces unphysical oscillations in the density profile. [5 6] One strategy to overcome these problems is the hybrid method, flux limiter scheme. [7 11] At discontinuities the flux limiter is close to 0, so the first order upwind scheme is applied; while for the smooth areas the flux limiter approaches 1, a second order or second order approximation term is used. The flux limiter technique has been introduced into the Single-Relaxation-Time (SRT LB simulations by Sofonea, et al. [7 8] They showed that such a scheme can successfully improve the numerical accuracy. However, to date, the flux limiter scheme has not been incorporated into the Multiple-Relaxation-Time (MRT LB method for compressible flows. In order to simulate the multiphase flows under shock wave action, in this work, we propose a flux limiter scheme with the splitting technique [12] for an MRT LB model. The rest of the paper is organized as follows: Details of our finite difference scheme based on the flux limiter and splitting techniques are given in Sec. 2. In Sec. 3 we present some simulation results with the new model, and compare the simulation results with the theoretical predictions. A brief conclusion is given in Sec Numerical Schemes Recently, we proposed a new MRT Finite Difference (FD lattice Boltzmann model for compressible flows. [13 14] This MRT model has better numerical stability and accuracy than its single-relaxation-time counterpart. In this MRT model the evolution of distribution function f i is governed by the following equation: where f i t + v f i iα = M 1 il Ŝ lk ( x ˆf k α eq ˆf k, (2 Supported by the Science Foundation of Laboratory of Computational Physics, Science Foundation of China Academy of Engineering Physics under Grant Nos. 2009A , 2009B , National Basic Research Program of China under Grant No. 2007CB815105, National Natural Science Foundation of China under Grant Nos , , and , and the Fundamental Research Funds for the Central Universities under Grant No. 2010YS03 Corresponding author, Xu Aiguo@iapcm.ac.cn Corresponding author, lyj@aphy.iphy.ac.cn c 2011 Chinese Physical Society and IOP Publishing Ltd
2 334 Communications in Theoretical Physics Vol. 56 cyc : (±1, 0, for 1 i 4, cyc : (±6, 0, for 5 i 8, (v i1, v i2 = 2(±1, ±1, for 9 i 12, 3 2 (±1, ±1, for 13 i 16. (3 Details of transformation matrix M and the equilibrium distribution function to the original publication. [13] ˆf eq i in the moment space are referred Fig. 1 Characteristic lines and corresponding projections in the x and y directions for the following distribution functions: (a is for f 1(x, t; (b is for f 9(x, t. In this paper, we propose a flux limiter scheme corresponding to the MRT model. The splitting technique is used in composing the flux limiter. Here, the splitting technique means that the one-dimensional flux limiter scheme is used in each direction (x or y for a multidimensional flux limiter scheme. Figure 1 shows the characteristic lines in the flux limiter scheme and corresponding projections in x and y directions for two distribution functions f 1 (x, t and f 9 (x, t, respectively. (J 1 x and (J 1 y are corresponding projections of node J 1 in the x and y directions. Let f n i,j be the value of the distribution function at time t at the node J along the direction i, we rewrite the evolution of f i at node J at time step t + dt as follows, where f n+1 i,j = fi,j n dt A i dx [F i,j+1/2 n x Fi,J 1/2 n x] dt A i dy [F i,j+1/2 n y Fi,J 1/2 n y] dtm 1 il Ŝ lk ( ˆf k,j n n,eq ˆf k,j, (4 1, for 1 i 4, 1/6, for 5 i 8, A i = 1/ (5 2, for 9 i 12, 2/3, for 13 i 16. Fi,J+1/2 n x (Fi,J 1/2 n x and Fi,J+1/2 n y (Fi,J 1/2 n y are x and y components of the outgoing (incoming flux at node J along the direction i, F n i,j+1/2 x = f n i (ix, iy ( A i dx F n i,j 1/2 x = f n i (ix A i v ix, iy Fi,J+1/2 n y = fi n (ix, iy + 1 ( 2 A i dy F n i,j 1/2 y = f n i (ix, iy A i v iy [f n i (ix + A iv ix, iy f n i (ix, iy]ψ x(ix, iy, (6a ( A i dx [f n i (ix, iy f n i (ix A i v ix, iy]ψ x (ix A i v ix, iy, (6b [f n i (ix, iy + A iv iy f n i (ix, iy]ψ y(ix, iy, (6c ( A i dy [f n i (ix, iy f n i (ix, iy A i v iy ]ψ y (ix, iy A i v iy. (6d The sufficient condition for the discrete space with second order accuracy and TVD property is that the limiter lies in the limited area, where the upper contour is the superbee limiter, and the lower contour is the minmod limiter. The
3 No. 2 Communications in Theoretical Physics 335 closer the chosen limiter to the upper contour, the smaller the dissipation caused by the limiter is, and the higher the resolution is, but the lower the stability and the convergence are. On the contrary, the closer the chosen limiter to the lower contour, the bigger the dissipation, and the lower the resolution, but the higher the stability and convergence. Based on an overall consideration of the resolution and stability, we adopt the Monitorized Central Difference (MCD flux limiter [15] in this paper. The flux limiters ψ x and ψ y are expressed as follows 0, θi n(ix, iy x 0, 2θi n ψ x (ix, iy = (ix, iy x, 0 θi n(ix, iy x 1 3, (1 + θi n (ix, iy 1 x/2, 3 θn i (ix, iy (7a x 3, 2, 3 θi n(ix, iy x, 0, θi n(ix, iy y 0, 2θi n ψ y (ix, iy = (ix, iy y, 0 θi n(ix, iy y 1 3, (1 + θi n (ix, iy 1 y/2, 3 θn i (ix, iy (7b y 3, 2, 3 θi n(ix, iy y, where the smoothness functions are θi n (ix, iy x = fn i (ix, iy fn i (ix A iv ix, iy fi n(ix + A iv ix, iy fi n, (ix, iy (8a θi n (ix, iy y = fn i (ix, iy fn i (ix, iy A iv iy fi n(ix, iy + A iv iy fi n. (ix, iy (8b The Lax Wendroff scheme is recovered for the flux limiter ψ x = ψ y = 1. The first order upwind scheme is recovered when ψ x = ψ y = 0. 3 Simulation Results The Richtmyer Meshkov (RM instability has been studied extensively in the linear and nonlinear regimes. Before going to LB simulations, we first summarize briefly some relevant results. For the model summary given below, it is useful to define some key quantities appearing in the models: a is the mixing layer amplitude, a b and a s are the bubble and spike amplitudes, k = 2π/λ is the wave number, λ is the initial perturbation wavelength, u is the velocity change across the interface, A 1 is the post-shock Atwood number, a 1 represents the post-shock initial perturbation amplitude. Zhang and Sohn [16] developed a model for the growth of RM unstable interface in the case of light-heavy transition. They proposed da dt = v k 2 v 0 a 1 t + max[0, (ka 1 2 (A ](kv 0 t 2 (9a da b dt = da dt A 1 kv0t k 2 v 0 a 1 t + 4[(ka /3(1 A 2 1 ](kv 0t 2, (9b da s dt = da dt + A 1 kv0 2t 1 + 2k 2 v 0 a 1 t + 4[(ka /3(1 A 2 1 ](kv 0t 2, (9c for the overall growth rate and the growth rates of bubble and spike, respectively. Here v 0 = k ua 1 a 1. Goncharov [17] extended the two-dimensional Layzer model to the A 1 1 case and obtained Sadot et al. [18] proposed the growth rate of the bubble and spike, da b dt = 3 + A 1 3(1 + A 1 kt. (10 da b/s dt = (1 + kv 0 tv (1 ± A 1 kv 0 t + [(1 ± A 1 /(1 + A 1 ][(kv 0 t 2 /2πC], (11 where C = 1/(3π for A 1 0.5, C = 1/(2π for A 1 0, and the upper (+ and lower ( signs in ± denote the bubble and spike, respectively. Sohn [9] extended the Zufiria model to arbitrary Atwood numbers, and the bubble velocity is da b dt = [ 3 + A1 3(1 + A 1 1 q + 2A 1 3(1 + A 1 q 2 ] 1 kt, (12 where q = q(a 1 is the root of the cubic polynomial (3 A 1 q 3 (21 + 9A 1 q 2 + (3 + 15A 1 q 4A 1 = 0.
4 336 Communications in Theoretical Physics Vol. 56 Here, the simulation results with the new model are compared to the results from Lax Wendroff scheme simulations, and to the predictions of models mentioned above. The following is the initial condition of our simulations: An incident shock wave, with the Mach number 1.2, traveling from the left side, hits an interface with sinusoidal perturbation. The initial macroscopic quantities are as follows: (ρ, u 1, u 2, p l = ( , , 0, , (ρ, u 1, u 2, p m = (1, 0, 0, 1, (ρ, u 1, u 2, p r = (5.04, 0, 0, 1, (13 where ρ is the density, u 1 and u 2 are the x and y components of fluid velocity, p is pressure, and the subscripts l, m, r indicate the left, middle, right regions of the whole domain. The initial sinusoidal perturbation at the interface is: x = 0.25 Nx dx cos(20πy, where the cycle of initial perturbation is 0.1, the amplitude is 0.008, Nx is grid number in the x direction, and dx is grid size. The following boundary conditions are imposed: (i inflow boundary is applied at the left side, (ii periodic boundary conditions are applied at the top and bottom boundaries, and (iii extrapolation technique is adopted at the right boundary. γ = 1.4 in the whole domain. Fig. 2 Comparison of the simulation results with flux limiter scheme and Lax Wendroff scheme at t = 0.1 and t = 0.3, respectively: dotted lines are for flux limiter scheme, solid lines are for Lax Wendroff scheme. Fig. 3 Comparison of the simulation results with flux limiter scheme and Lax Wendroff scheme: amplitude and growth rate changes with time. Figure 2 shows a comparison of the simulation results with flux limiter scheme and Lax Wendroff scheme at t = 0.1 and t = 0.3, respectively. The dotted lines correspond to simulation results with the flux limiter scheme, and solid lines represent the simulation results with Lax Wendroff scheme. The abscissa is for ix, and the vertical axis is for iy, where ix and iy are the indexes of lattice node in the x- and y-directions. The parameters are dx = dy = 0.001, dt = 10 5, s 5 = 10 3, and 10 5 for the other relaxation parameters. We adopt the same parameters in the following simulations. Figure 3 shows the changes of perturbation amplitudes and growth rates with time simulated with the two different schemes. From these figures, we can find that the flux limiter scheme is efficient in decreasing the oscillations, and the computed interfaces are more smooth and consistent with physical analysis. The smooth interface is beneficial to accurately
5 No. 2 Communications in Theoretical Physics 337 estimate the disturbance amplitude. Fig. 4 Snapshots of RM instability: density contours at t = 0, t = 0.05, t = 0.3, t = 0.5, t = 0.7, t = 1.1, respectively. From black to white the gray level corresponds to the increase of values. In Fig. 4, we show the density contours simulated with the flux limiter scheme at t = 0, t = 0.05, t = 0.3, t = 0.5, t = 0.7, t = 1.1, respectively. In the first snapshot corresponding to the initial condition (13, the gray levels on the left side of the interface disturbance are not obvious, because the density difference of these regions is too small compared with the value of right region. From the density distribution, we can clearly find the locations of the interface and transmission wave. At t = 0.7, the transmission wave has been out of the computational domain. In the last snapshot (t = 1.1, we also find a discontinuity near ix = 470, which is a reflected shock wave from the right boundary. It should be noted that the extrapolation technique is very difficult to ensure complete nonreflective outflow. But, since the reflection wave has not reach the location of interface, so the calculation results are still valid and impervious. How to construct an appropriate and complete non-reflective outflow boundary condition for this model is still an open problem for us. In the impact process, the heavy and light fluids gradually penetrate into each other, the light fluid rises to form a bubble and the heavy fluid falls to generate a spike. At time t = 1.1, the interface disturbance has obtained full development. At the top of the spike, due to the Kelvin Helmholtz instability, a mushroom shape is formed eventually. Fig. 5 Simulation results of flux limiter scheme: bubble and spike. (a, (b, (c, (d are the bubble amplitude, spike amplitude, bubble and spike growth rates, respectively. In Fig. 5, we present the amplitudes and growth rates of bubble and spike, where Figs. 5(a and 5(b describe the changes of bubble and spike amplitudes with time, Figs. 5(c and 5(d show the comparisons of bubble and spike growth rates with theoretical predictions, respectively. The location of the shocked, unperturbed interface is used as a reference for the determination of the bubble and spike amplitudes. The distance from the unperturbed interface to the bubble tip is the bubble amplitude a b, while the distance from the unperturbed interface to the spike tip is the spike amplitude a s. When the shock wave passes through the initial interface disturbance area, due
6 338 Communications in Theoretical Physics Vol. 56 to the compression, the interface produces a small deformation. The bubble and spike amplitudes reduce slightly, as shown in Figs. 5(a and 5(b. In Fig. 5(c, we compare our simulated bubble growth rate with the predictions of Zhang Sohn model (Eq. (9b, Goncharov model (Eq. (10, model by Sadot et al. (Eq. (11 and Sohn Zufiria model (Eq. (12. The bubble growth rates predicted by the model of Sadot et al. and the Sohn-Zufiria model are in very close agreement in the late time and are closest to the simulation bubble growth rate. The Goncharov model overpredicts the bubble growth rate, and the Zhang Sohn model underpredicts the bubble growth rate. In Figs. 5(d, we compare the simulated spike growth rate with the predictions of Zhang Sohn model (Eq. (11 and the model by Sadot et al. (Eq. (11. The model by Sadot et al. overpredicts the spike growth rate, and the Zhang Sohn model underpredicts the spike growth rate. The agreement between the theories and the simulations for the bubble is better than that for the spike. And it is consistent with a comment proposed by Zhang Sohn, the prediction for the spike growth rate becomes less accurate than that for the bubble at later times. [20] The simulation results are also accordant with those by other numerical methods. [21] 4 Conclusions We propose a flux limiter scheme with splitting technique for a multiple-relaxation-time lattice Boltzmann method. The new model works for both high speed and low speed thermal compressible flows under shack wave reaction. It is shown that, the new scheme can efficiently decrease the unphysical oscillations and numerical diffusion around the interfaces. Some interface problems which are difficult to handle at the macroscopic level can be modeled more naturally in the LB method. The new model was applied to simulate a traditional difficult problem, the Richtmyer Meshkov instability. The obtained growth rates of bubble and spike are compared with predictions of the Zhang-Sohn model, Goncharov model, Sohn Zufiria model and the model by Sadot et al. The numerical results are shown to be in good agreement with those from analytic calculations and other macroscopic approaches. Since the modeling of LB method includes intrinsically non-equilibrium effects, it has more potential advantages than the macroscopic methods in studying the non-equilibrium behaviors of fluid instability problems, such as the Richtmyer Meshkov mixing. References [1] S. Succi, The Lattice Boltzmann Equation for Fluid Dynamics and Beyond, Oxford University Press, New York (2001. [2] Ai-Guo Xu, G. Gonnella, and A. Lamura, Phys. Rev. E 67 ( ; Phys. Rev. E 74 ( ; Physica A 331 ( ; Physica A 344 ( [3] V. Sofonea and R.F. Sekerka, J. Comp. Phys. 184 ( ; Int. J. Mod. Phys. C 16 ( [4] A. Cristea and V. Sofonea, Int. J. Mod. Phys. C 14 ( [5] E.F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics, Springer Verlag, Berlin (1999. [6] S. Teng, Y. Chen, and H. Ohashi, Int. J. Heat and Fluid Flow 21 ( [7] A. Cristea and V. Sofonea, Proceedings of the Romanian Academy, Series A 4 ( ; Central European J. Phys. 2 ( [8] V. Sofonea, A. Lamura, G. Gonnella, and A. Cristea, Phys. Rev. E 70 ( [9] V. Sofonea, J. Comput. Phys. 228 ( [10] G. Gonnella, A. Lamura, and V. Sofonea, Eur. Phys. J. Special Topics 171 ( [11] A. Cristea, Int. J. Mod. Phys. C 17 ( [12] S.P. Spekreijse, Math. Comput. 49 ( [13] F. Chen, A.G. Xu, G.C. Zhang, Y.J. Li, and S. Succi, Europhys. Lett. 90 ( [14] F. Chen, A.G. Xu, G.C. Zhang, and Y.J. Li, Commun. Theor. Phys. 55 ( [15] E. F. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics, Second Edition, Springer Verlag, Berlin (1999. [16] Q. Zhang and S. Sohn, Phys. Fluids 9 ( [17] V.N. Goncharov, Phys. Rev. Lett. 88 ( [18] O. Sadot, L. Erez, U. Alon, D. Oron, and L.A. Levin, Phys. Rev. Lett. 80 ( [19] S. Sohn, Phys. Rev. E 67 ( [20] Q. Zhang and S. Sohn, Appl. Math. Lett. 10 ( [21] M. Latini, O. Schilling, and W.S. Don, Phys. Fluids 19 (
A Vortex Model for Studying the Effect of Shock Proximity on Richtmyer-Meshkov Instability at High Mach Number
A Vortex Model for Studying the Effect of Shock Proximity on Richtmyer-Meshkov Instability at High Mach Number H. F. Robey, S. G. Glendinning, J. A. Greenough, & S. V. Weber Lawrence Livermore National
More information470 SADOT ET AL. Vol. 127
THE ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES, 127:469È473, 2000 April ( 2000. The American Astronomical Society. All rights reserved. Printed in U.S.A. STUDIES ON THE NONLINEAR EVOLUTION OF THE RICHTMYER-MESHKOV
More informationOn The Collapse of a Gas Cavity by an Imploding Molten Lead Shell and Richtmyer-Meshkov Instability
On The Collapse of a Gas Cavity by an Imploding Molten Lead Shell and Richtmyer-Meshkov Instability Victoria Suponitsky, Sandra Barsky, and Aaron Froese General Fusion Inc., 8-368 Bonneville Place, Burnaby,
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More informationPHYSICS OF FLUIDS 19,
PHYSICS OF FLUIDS 19, 04104 007 High-resolution simulations and modeling of reshocked single-mode Richtmyer-Meshkov instability: Comparison to experimental data and to amplitude growth model predictions
More informationLos Alamos National Laboratory Hydrodynamic Methods Applications and Research 1 LA-UR
Rayleigh-Taylor instability is generated when a heavy fluid sits above a lighter fluid in a gravitational field. The flow behavior is described in terms of bubbles of light fluid rising into the heavier
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationLattice Boltzmann approach to liquid - vapour separation
Lattice Boltzmann approach to liquid - vapour separation T.Biciușcă 1,2, A.Cristea 1, A.Lamura 3, G.Gonnella 4, V.Sofonea 1 1 Centre for Fundamental and Advanced Technical Research, Romanian Academy Bd.
More informationLattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel
Commun. Theor. Phys. 56 (2011) 756 760 Vol. 56, No. 4, October 15, 2011 Lattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel QIU Bing ( ), 1, TAN Hui-Li ( Û), 2 and
More informationNumerical investigation of shock wave dense particles cloud interaction
25 th ICDERS August 2 7, 215 Leeds, UK Numerical investigation of shock wave dense particles cloud interaction Pavel S. Utkin Moscow Institute of Physics and Technology Dolgoprudny, Moscow region, Russia
More informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More information68 Guo Wei-Bin et al Vol. 12 presented, and are thoroughly compared with other numerical data with respect to the Strouhal number, lift and drag coeff
Vol 12 No 1, January 2003 cfl 2003 Chin. Phys. Soc. 1009-1963/2003/12(01)/0067-08 Chinese Physics and IOP Publishing Ltd Lattice-BGK simulation of a two-dimensional channel flow around a square cylinder
More informationThermal lattice Bhatnagar-Gross-Krook model for flows with viscous heat dissipation in the incompressible limit
PHYSICAL REVIEW E 70, 066310 (2004) Thermal lattice Bhatnagar-Gross-Krook model for flows with viscous heat dissipation in the incompressible limit Yong Shi, T. S. Zhao,* Z. L. Guo Department of Mechanical
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
More informationNoise Shielding Using Acoustic Metamaterials
Commun. Theor. Phys. (Beijing, China) 53 (2010) pp. 560 564 c Chinese Physical Society and IOP Publishing Ltd Vol. 53, No. 3, March 15, 2010 Noise Shielding Using Acoustic Metamaterials LIU Bin ( Ê) and
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationSimulation of Kelvin Helmholtz Instability with Flux-Corrected Transport Method
Commun. Theor. Phys. (Beijing, China) 51 (2009) pp. 909 913 c Chinese Physical Society and IOP Publishing Ltd Vol. 51, No. 5, May 15, 2009 Simulation of Kelvin Helmholtz Instability with Flux-Corrected
More informationApplication of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationViscous non-linear theory of Richtmyer-Meshkov Instability. Abstract
Viscous non-linear theory of Richtmyer-Meshkov Instability Pierre Carles and Stéphane Popinet Laboratoire de Modélisation en Mécanique, Université Pierre et Marie Curie, Case 162, 4 place Jussieu, 75252
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationStructure Formation and Particle Mixing in a Shear Flow Boundary Layer
Structure Formation and Particle Mixing in a Shear Flow Boundary Layer Matthew Palotti palotti@astro.wisc.edu University of Wisconsin Center for Magnetic Self Organization Ellen Zweibel University of Wisconsin
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More information550 XU Hai-Bo, WANG Guang-Rui, and CHEN Shi-Gang Vol. 37 the denition of the domain. The map is a generalization of the standard map for which (J) = J
Commun. Theor. Phys. (Beijing, China) 37 (2002) pp 549{556 c International Academic Publishers Vol. 37, No. 5, May 15, 2002 Controlling Strong Chaos by an Aperiodic Perturbation in Area Preserving Maps
More informationNew Feedback Control Model in the Lattice Hydrodynamic Model Considering the Historic Optimal Velocity Difference Effect
Commun. Theor. Phys. 70 (2018) 803 807 Vol. 70, No. 6, December 1, 2018 New Feedback Control Model in the Lattice Hydrodynamic Model Considering the Historic Optimal Velocity Difference Effect Guang-Han
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More information2 The incompressible Kelvin-Helmholtz instability
Hydrodynamic Instabilities References Chandrasekhar: Hydrodynamic and Hydromagnetic Instabilities Landau & Lifshitz: Fluid Mechanics Shu: Gas Dynamics 1 Introduction Instabilities are an important aspect
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationPredictor-Corrector Finite-Difference Lattice Boltzmann Schemes
Applied Mathematical Sciences, Vol. 9, 2015, no. 84, 4191-4199 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.5138 Predictor-Corrector Finite-Difference Lattice Boltzmann Schemes G. V.
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationLattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocks
Lattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocs Maysam Saidi, Hassan Basirat Tabrizi, Member, IAENG, and Reza Sepahi Samian Abstract Lattice Boltzmann method
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationThe Discontinuous Galerkin (dg) method is a popular numerical method that has received much attention
Fluid Dynamics and Co-located Conferences June 24-27, 2013, San Diego, CA 21st AIAA Computational Fluid Dynamics Conference AIAA 2013-2595 Discontinuous Galerkin method for multifluid Euler equations Marc
More informationSolution of the Hirota Equation Using Lattice-Boltzmann and the Exponential Function Methods
Advanced Studies in Theoretical Physics Vol. 11, 2017, no. 7, 307-315 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2017.7418 Solution of the Hirota Equation Using Lattice-Boltzmann and the
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationSecond Order Positive Schemes by means of Flux Limiters for the Advection Equation
Second Order Positive Schemes by means of Flux Limiters for the Advection Equation Riccardo Fazio and Alessandra Jannelli Abstract In this paper, we study first and second order positive numerical methods
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationA gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions
Acta Mech. Sin. 2017) 334):733 741 DOI 10.1007/s10409-017-0695-2 RESEARCH PAPER A gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions Xiaodong Ren 1
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationA General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations
A General Technique for Eliminating Spurious Oscillations in Conservative Schemes for Multiphase and Multispecies Euler Equations Ronald P. Fedkiw Xu-Dong Liu Stanley Osher September, 2 Abstract Standard
More informationCascade of Random Rotation and Scaling in a Shell Model Intermittent Turbulence
Commun. Theor. Phys. (Beijing, China) 46 (2006) pp. 1074 1080 c International Academic Publishers Vol. 46, No. 6, December 15, 2006 Cascade of Random Rotation and Scaling in a Shell Model Intermittent
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationPropagation of Lorentz Gaussian Beams in Strongly Nonlocal Nonlinear Media
Commun. Theor. Phys. 6 04 4 45 Vol. 6, No., February, 04 Propagation of Lorentz Gaussian Beams in Strongly Nonlocal Nonlinear Media A. Keshavarz and G. Honarasa Department of Physics, Faculty of Science,
More informationDouble-distance propagation of Gaussian beams passing through a tilted cat-eye optical lens in a turbulent atmosphere
Double-distance propagation of Gaussian beams passing through a tilted cat-eye optical lens in a turbulent atmosphere Zhao Yan-Zhong( ), Sun Hua-Yan( ), and Song Feng-Hua( ) Department of Photoelectric
More informationOn the Early Development of Dispersion in Flow through a Tube with Wall Reactions
World Academy of Science, Engineering and Technology 9 7 On the Early Development of Dispersion in Flow through a Tube with Wall Reactions M. W. Lau, and C. O. Ng Abstract This is a study on numerical
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationThe Equivalent Differential Equation of the Gas-Kinetic BGK Scheme
The Equivalent Differential Equation of the Gas-Kinetic BGK Scheme Wei Gao Quanhua Sun Kun Xu Abstract The equivalent partial differential equation is derived for the gas-kinetic BGK scheme (GKS) [Xu J.
More informationSchemes for Mixture Modeling
DEPARTMENT OF ENERGETICS Consistency of Multiplerelaxation-time Lattice Boltzmann Schemes for Mixture Modeling Invited Talk (Prof. Taku Ohwada) Department of Aeronautics and Astronautics Graduate School
More informationEffects of Atomic Mixing in Inertial Confinement Fusion by Multifluid Interpenetration Mix Model
Commun. Theor. Phys. (Beijing, China) 52 (2009) pp. 1102 1106 c Chinese Physical Society and IOP Publishing Ltd Vol. 52, No. 6, December 15, 2009 Effects of Atomic Mixing in Inertial Confinement Fusion
More informationPalindromic Discontinuous Galerkin Method
Palindromic Discontinuous Galerkin Method David Coulette, Emmanuel Franck, Philippe Helluy, Michel Mehrenberger, Laurent Navoret To cite this version: David Coulette, Emmanuel Franck, Philippe Helluy,
More informationDetonation Simulations with a Fifth-Order TENO Scheme
Commun. Comput. Phys. doi: 10.4208/cicp.OA-2018-0008 Vol. x, No. x, pp. 1-37 xxx 201x Detonation Simulations with a Fifth-Order TENO Scheme Haibo Dong 1, Lin Fu 2, Fan Zhang 3,, Yu Liu 4 and Jun Liu 1
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
More information3.4. Monotonicity of Advection Schemes
3.4. Monotonicity of Advection Schemes 3.4.1. Concept of Monotonicity When numerical schemes are used to advect a monotonic function, e.g., a monotonically decreasing function of x, the numerical solutions
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationAdiabatic piston - Draft for MEMPHYS
diabatic piston - Draft for MEMPHYS 4th ovember 2006 Christian GRUBER (a) and nnick LESE (b) (a) Institut de héorie des Phénomènes Physiques École Polytechnique Fédérale de Lausanne, CH-05 Lausanne, Switzerland
More informationOn the lattice Boltzmann method for multiphase flows
Center for Turbulence Research Annual Research Briefs 2009 377 On the lattice Boltzmann method for multiphase flows By S. H. Kim AND H. Pitsch 1. Motivation and objectives The lattice Boltzmann (LB) method
More informationA New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent Sources
Commun. Theor. Phys. Beijing, China 54 21 pp. 1 6 c Chinese Physical Society and IOP Publishing Ltd Vol. 54, No. 1, July 15, 21 A New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent
More informationComparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models
Comparison of Numerical Solutions for the Boltzmann Equation and Different Moment Models Julian Koellermeier, Manuel Torrilhon October 12th, 2015 Chinese Academy of Sciences, Beijing Julian Koellermeier,
More information2 In this paper, we give a brief summary of the front tracking algorithm for axisymmetric ow. We report simulations of spherical shock refraction by a
A Fast Algorithm for Moving Interface Problems S. Dutta 1,J.Glimm 1 2,J.W.Grove 3,D.H.Sharp 4, and Y. Zhang 1 1 Department of Applied Mathematics and Statistics, University at Stony Brook, Stony Brook,
More information2 Three Dimensional Finite Volume Gas-Kinetic Schemes
Gas-Kinetic Finite Volume Methods K. Xu, L. Martinelli, A. Jameson Department of M.A.E, Princeton University, Prificeton N.J. 08544, USA 1 Introduction Gas-kinetic schemes developed from the BGK model
More informationNew Homoclinic and Heteroclinic Solutions for Zakharov System
Commun. Theor. Phys. 58 (2012) 749 753 Vol. 58, No. 5, November 15, 2012 New Homoclinic and Heteroclinic Solutions for Zakharov System WANG Chuan-Jian ( ), 1 DAI Zheng-De (à ), 2, and MU Gui (½ ) 3 1 Department
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationGeneralized Local Equilibrium in the Cascaded Lattice Boltzmann Method. Abstract
Accepted for publication on Physical Review E (R), code ERR1034 Generalized Local Equilibrium in the Cascaded Lattice Boltzmann Method Pietro Asinari Department of Energetics, Politecnico di Torino, Corso
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationGas-Kinetic Relaxation (BGK-Type) Schemes for the Compressible Euler Equations
Gas-Kinetic Relaxation (BGK-Type) Schemes for the Compressible Euler Equations Iizln Xu* and Antony ~arneson~ Department of Mechanical and Aerospace Engineering Princeton University Princeton, New Jersey
More informationThe First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
Commun. Theor. Phys. (Beijing, China 35 (2 pp. 42 46 c International Academic Publishers Vol. 35, No. 4, April 5, 2 The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
More informationComputational issues and algorithm assessment for shock/turbulence interaction problems
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln NASA Publications National Aeronautics and Space Administration 2007 Computational issues and algorithm assessment for shock/turbulence
More informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationAdvection algorithms II. Flux conservation, subgrid models and flux limiters
Chapter 4 Advection algorithms II. Flu conservation, subgrid models and flu limiters In this chapter we will focus on the flu-conserving formalism of advection algorithms, and we shall discuss techniues
More informationEffects of ion temperature on electrostatic solitary structures in nonthermal plasmas
PHYSICAL REVIEW E VOLUME 55, NUMBER FEBRUARY 1997 Effects of ion temperature on electrostatic solitary structures in nonthermal plasmas A. A. Mamun Department of Physics, Jahangirnagar University, Savar,
More informationDynamical behaviour of a controlled vibro-impact system
Vol 17 No 7, July 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(07)/2446-05 Chinese Physics B and IOP Publishing Ltd Dynamical behaviour of a controlled vibro-impact system Wang Liang( ), Xu Wei( ), and
More informationCoupled double-distribution-function lattice Boltzmann method for the compressible Navier-Stokes equations
PHYSICAL REVIEW E 76, 5675 7 Coupled double-distribution-function lattice Boltzmann method for the compressible Navier-Stokes equations Q. Li, Y. L. He,* Y. Wang, and W. Q. Tao State Key Laboratory of
More informationNon-linear Methods for Scalar Equations
Non-linear Methods for Scalar Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 3, 04 / 56 Abstract
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationStable Propagating Waves and Wake Fields in Relativistic Electromagnetic Plasma
Commun. Theor. Phys. (Beijing, China) 49 (2008) pp. 753 758 c Chinese Physical Society Vol. 49, No. 3, March 15, 2008 Stable Propagating Waves and Wake Fields in Relativistic Electromagnetic Plasma XIE
More informationProlongation structure for nonlinear integrable couplings of a KdV soliton hierarchy
Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy Yu Fa-Jun School of Mathematics and Systematic Sciences, Shenyang Normal University, Shenyang 110034, China Received
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationTime-delay feedback control in a delayed dynamical chaos system and its applications
Time-delay feedback control in a delayed dynamical chaos system and its applications Ye Zhi-Yong( ), Yang Guang( ), and Deng Cun-Bing( ) School of Mathematics and Physics, Chongqing University of Technology,
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationFFT-LB Modeling of Thermal Liquid-Vapor System
Commun. Theor. Phys. 57 (2012) 681 694 Vol. 57, No. 4, April 15, 2012 FFT-LB Modeling of Thermal Liquid-Vapor System GAN Yan-Biao ( Á), 1,2,3 XU Ai-Guo (ÆÇÁ), 2, ZHANG Guang-Cai (¾), 2 and LI Ying-Jun
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationGas Dynamics Equations: Computation
Title: Name: Affil./Addr.: Gas Dynamics Equations: Computation Gui-Qiang G. Chen Mathematical Institute, University of Oxford 24 29 St Giles, Oxford, OX1 3LB, United Kingdom Homepage: http://people.maths.ox.ac.uk/chengq/
More informationARTICLE IN PRESS Mathematical and Computer Modelling ( )
Mathematical and Computer Modelling Contents lists available at ScienceDirect Mathematical and Computer Modelling ournal homepage: wwwelseviercom/locate/mcm Total variation diminishing nonstandard finite
More information