In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interfac...In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.展开更多
A generalized upwind scheme with fractional steps for 3-D mathematical models of convection dominating groundwater quality is suggested. The mass transport equation is split into a convection equation and a dispersive...A generalized upwind scheme with fractional steps for 3-D mathematical models of convection dominating groundwater quality is suggested. The mass transport equation is split into a convection equation and a dispersive equation. The generalized upwind scheme is used to solve the convection equation and the finite element method is used to compute the dispersive equation.These procedures which not only overcome the phenomenon of the negative concentration and numerical dispersion appear frequently with normal FEM or FDM to solve models of convection dominating groundwater transport but also avoid the step for computing each node velocity give a more suitable method to calculate the concentrations of the well points.展开更多
Microbial cultures are comprised of heterogeneous cells that differ according to their size and intracellular concentrations of DNA, proteins and other constituents. Because of the included level of details, multi-var...Microbial cultures are comprised of heterogeneous cells that differ according to their size and intracellular concentrations of DNA, proteins and other constituents. Because of the included level of details, multi-variable cell population balance models (PBMs) offer the most general way to describe the complicated phenomena associated with cell growth, substrate consumption and product formation. For that reason, solving and understanding of such models are essential to predict and control cell growth in the processes of biotechnological interest. Such models typically consist of a partial integro-differential equation for describing cell growth and an ordinary integro-differential equation for representing substrate consumption. However, the involved mathematical complexities make their numerical solutions challenging for the given numerical scheme. In this article, the central upwind scheme is applied to solve the single-variate and bivariate cell population balance models considering equal and unequal partitioning of cellular materials. The validity of the developed algorithms is verified through several case studies. It was found that the suggested scheme is more reliable and effective.展开更多
In this paper, we discuss the uniform convergence of the simple upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh for solving a singularly perturbed Robin boundary value problem, and investigate the m...In this paper, we discuss the uniform convergence of the simple upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh for solving a singularly perturbed Robin boundary value problem, and investigate the midpoint upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh to achieve better uniform convergence. The elaborate ε-uniform pointwise estimates are proved by using the comparison principle and barrier functions. The numerical experiments support the theoretical results for the schemes on the meshes.展开更多
The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase...The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.展开更多
A parallelized upwind flux splitting scheme for supersonic reacting flows on hybrid meshes is presented. The complexity of super/hyper-sonic combustion flows makes it necessary to establish solvers with higher resolut...A parallelized upwind flux splitting scheme for supersonic reacting flows on hybrid meshes is presented. The complexity of super/hyper-sonic combustion flows makes it necessary to establish solvers with higher resolution and efficiency for multi-component Euler/N-S equations. Hence, a spatial second-order van Leer type flux vector splitting scheme is established by introducing auxiliary points in interpolation, and a domain decomposition method used on unstructured hybrid meshes for obtaining high calculating efficiency. The numerical scheme with five-stage Runge-Kutta time step method is implemented to the simulation of combustion flows, including the supersonic hydrogen/air combustion and the normal injection of hydrogen into reacting flows. Satisfying results are obtained compared with limited references.展开更多
We study the L^l-error estimates for the upwind scheme to the linear advection equations with a piecewise constant coefficients modeling linear waves crossing interfaces. Here the interface condition is immersed into ...We study the L^l-error estimates for the upwind scheme to the linear advection equations with a piecewise constant coefficients modeling linear waves crossing interfaces. Here the interface condition is immersed into the upwind scheme. We prove that, for initial data with a bounded variation, the numerical solution of the immersed interface upwind scheme converges in L^l-norm to the differential equation with the corresponding interface condition. We derive the one-halfth order L^l-error bounds with explicit coefficients following a technique used in [25]. We also use some inequalities on binomial coefficients proved in a consecutive paper [32].展开更多
Nowadays,the upwind schemes are in a rapid development to capture shock accurately.However,these upwind schemes’properties at low speeds,such as their reconstruction scheme dependencies,grid dependencies,and Mach num...Nowadays,the upwind schemes are in a rapid development to capture shock accurately.However,these upwind schemes’properties at low speeds,such as their reconstruction scheme dependencies,grid dependencies,and Mach number dependencies,are concerned by few people.In this paper,a systematic study on their low speeds’issues is conducted.Through a series of tests,we can find that most parameter-free upwind schemes,widely used in practice today,are not applicable to low speeds’simulations.In contrast,SLAU and SLAU2 can give reliable results.Also,the upwind scheme’s influence on the accuracy is stronger than the reconstruction scheme’s influence at low speeds.展开更多
Aim To construct a third order upwind scheme for convection equation. Methods Upwind Lagrange interpolation was used. Results and Conclusion The schemes L p stability for p∈ is proved. Numerical exam...Aim To construct a third order upwind scheme for convection equation. Methods Upwind Lagrange interpolation was used. Results and Conclusion The schemes L p stability for p∈ is proved. Numerical examples show that performance of the third order upwind scheme is better than that of most second order schemes.展开更多
The importance of eliminating errors in grid-metric evaluation for highorder difference schemes has been widely recognized in recent years,and it is known from the proof by Vinokur and Yee(NASA TM 209598,2000)that whe...The importance of eliminating errors in grid-metric evaluation for highorder difference schemes has been widely recognized in recent years,and it is known from the proof by Vinokur and Yee(NASA TM 209598,2000)that when conservative derivations of grid metric are used by Thomas,Lombard and Neier(AIAA J.,1978,17(10)and J.Spacecraft and rocket,1990,27(2)),errors caused by metric evaluation could be eliminated by linear schemes when flux splitting is not considered.According to the above achievement,central schemes without the use of flux splitting could fulfill the requirement of error elimination.Difficulties will arise for upwind schemes to attain the objective when the splitting is considered.In this study,further investigations are made on three aspects:Firstly,an idea of central scheme decomposition is introduced,and the procedure to derive the central scheme is proposed to evaluate grid metrics only.Secondly,the analysis has been made on the requirement of flux splitting to acquire free-stream preservation,and a Lax-Friedrichs-type splitting scheme is proposed as an example.Discussions about current study with that by Nonomura et al.(Computers and Fluids,2015,107)have been made.Thirdly,for halfnode-or mixed-type schemes,interpolations should be used to derive variables at half nodes.The requirement to achieve metric identity on this situation is analyzed and an idea of directionally consistent interpolation is proposed,which is manifested to be indispensable to avoid violations of metric identity and to eliminate metric-caused errors thereafter.Two numerical problems are tested,i.e.,the free-stream and vortex preservation onwavy,largely randomized and triangular-like grids.Numerical results validate aforementioned theoretical outcomes.展开更多
A bounded high order upwind scheme is presented for the modified Burgers’equation by using the normalized-variable formulation in the finite volume framework.The characteristic line of the present scheme in the norma...A bounded high order upwind scheme is presented for the modified Burgers’equation by using the normalized-variable formulation in the finite volume framework.The characteristic line of the present scheme in the normalizedvariable diagram is designed on the Hermite polynomial interpolation.In order to suppress unphysical oscillations,the present scheme respects both the TVD(total variational diminishing)constraint and CBC(convection boundedness criterion)condition.Numerical results demonstrate the present scheme possesses good robustness and high resolution for the modified Burgers’equation.展开更多
In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial deriv...In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.展开更多
WT5,5”BX] A new class of numerical schemes is proposed to solve convection diffusion equations by combining the upwind technique and the method of operator splitting. For every time step, the multi dimensional approx...WT5,5”BX] A new class of numerical schemes is proposed to solve convection diffusion equations by combining the upwind technique and the method of operator splitting. For every time step, the multi dimensional approximation is performed in several independent directions alternatively, while the upwind technique is applied to treat the convection term in every individual direction. This scheme possesses maximum principle. Stability and convergence are analysed by energy method.[WT5,5”HZ]展开更多
An upwind weak Galerkin finite element scheme was devised and analyzed in this article for convection-dominated Oseen equations. The numerical algorithm was based on the weak Galerkin method enhanced by upwind stabili...An upwind weak Galerkin finite element scheme was devised and analyzed in this article for convection-dominated Oseen equations. The numerical algorithm was based on the weak Galerkin method enhanced by upwind stabilization. The resulting finite element scheme uses equal-order, say k, polynomial spaces on each element for the velocity and the pressure unknowns. With finite elements of order k≥1, the numerical solutions are proved to converge at the rate of O(h^(k+1/2)) in an energy-like norm for convection-dominated Oseen equations. Numerical results are presented to demonstrate the accuracy and effectiveness of the upwind weak Galerkin scheme.展开更多
The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the ...The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the bow shock.To improve the accuracy of aero-heating predictions,three multi-dimensional approaches on unstructured grids are coupled in our Reynolds-averaged Navier-Stokes(RANS)solver,including multi-dimensional upwind flux reconstruction(MUP),multi-dimensional limiter(MLP-u2)and multi-dimensional gradient reconstruction(MLR).The coupled multi-dimensional RANS solver is validated by several typical verification and validation(V&V)cases,including hypersonic flows over a cylinder,a blunt biconic,and a double-ellipsoid,with commonly used prism/tetrahedral hybrid grids.Finally,the coupled multi-dimensional solver is applied to simulating the heat flux distribution over a 3D engineering configuration,i.e.a Hermes-like space shuttle model.The obtained numerical results are compared with experimental data.The predicted results demonstrate that the coupled multi-dimensional approach has a good prediction capability on aerodynamic heating over a wide range of complex engineering configurations.展开更多
Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind sch...Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.展开更多
A fully discrete energy stability analysis is carried out for linear advection-diffusion problems discretized by generalized upwind summation-by-parts(upwind gSBP)schemes in space and implicit-explicit Runge-Kutta(IME...A fully discrete energy stability analysis is carried out for linear advection-diffusion problems discretized by generalized upwind summation-by-parts(upwind gSBP)schemes in space and implicit-explicit Runge-Kutta(IMEX-RK)schemes in time.Hereby,advection terms are discretized explicitly,while diffusion terms are solved implicitly.In this context,specific combinations of space and time discretizations enjoy enhanced stability properties.In fact,if the first-and second-derivative upwind gSBP operators fulfill a compatibility condition,the allowable time step size is independent of grid refinement,although the advective terms are discretized explicitly.In one space dimension it is shown that upwind gSBP schemes represent a general framework including standard discontinuous Galerkin(DG)schemes on a global level.While previous work for DG schemes has demonstrated that the combination of upwind advection fluxes and the central-type first Bassi-Rebay(BR1)scheme for diffusion does not allow for grid-independent stable time steps,the current work shows that central advection fluxes are compatible with BR1 regarding enhanced stability of IMEX time stepping.Furthermore,unlike previous discrete energy stability investigations for DG schemes,the present analysis is based on the discrete energy provided by the corresponding SBP norm matrix and yields time step restrictions independent of the discretization order in space,since no finite-element-type inverse constants are involved.Numerical experiments are provided confirming these theoretical findings.展开更多
In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the...In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the non physical oscillation by means of the group velocity control. The scheme is used to simulate the interactions of shock density stratified interface and the disturbed interface developing to vortex rollers. Numerical results are satisfactory.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.11925204).
文摘In this paper,we develop a fourth-order conservative wavelet-based shock-capturing scheme.The scheme is constructed by combining a wavelet collocation upwind method with the monotonic tangent of hyperbola for interface capturing(THINC)technique.We employ boundary variation diminishing(BVD)reconstruction to enhance the scheme’s effectiveness in handling shocks.First,we prove that wavelet collocation upwind schemes based on interpolating wavelets can be reformulated into a conservative form within the framework of wavelet theory,forming the foundation of the proposed scheme.The new fourthorder accurate scheme possesses significantly better spectral resolution than the fifth-and even seventh-order WENO-Z(weighted essentially non-oscillatory)schemes over the entire wave-number range.Moreover,the inherent low-pass filtering property of the wavelet bases allows them to filter high-frequency numerical oscillations,endowing the wavelet upwind scheme with robustness and accuracy in solving problems under extreme conditions.Notably,due to the wavelet multiresolution approximation,the proposed scheme possesses a distinctive shape-preserving property absent in the WENO-Z schemes and the fifth-order schemes with BVD reconstruction based on polynomials.Furthermore,compared to the fifth-order scheme with BVD reconstruction based on polynomials—which is significantly superior to the WENO schemes—the proposed scheme further enhances the ability to capture discontinuities.
文摘A generalized upwind scheme with fractional steps for 3-D mathematical models of convection dominating groundwater quality is suggested. The mass transport equation is split into a convection equation and a dispersive equation. The generalized upwind scheme is used to solve the convection equation and the finite element method is used to compute the dispersive equation.These procedures which not only overcome the phenomenon of the negative concentration and numerical dispersion appear frequently with normal FEM or FDM to solve models of convection dominating groundwater transport but also avoid the step for computing each node velocity give a more suitable method to calculate the concentrations of the well points.
文摘Microbial cultures are comprised of heterogeneous cells that differ according to their size and intracellular concentrations of DNA, proteins and other constituents. Because of the included level of details, multi-variable cell population balance models (PBMs) offer the most general way to describe the complicated phenomena associated with cell growth, substrate consumption and product formation. For that reason, solving and understanding of such models are essential to predict and control cell growth in the processes of biotechnological interest. Such models typically consist of a partial integro-differential equation for describing cell growth and an ordinary integro-differential equation for representing substrate consumption. However, the involved mathematical complexities make their numerical solutions challenging for the given numerical scheme. In this article, the central upwind scheme is applied to solve the single-variate and bivariate cell population balance models considering equal and unequal partitioning of cellular materials. The validity of the developed algorithms is verified through several case studies. It was found that the suggested scheme is more reliable and effective.
文摘In this paper, we discuss the uniform convergence of the simple upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh for solving a singularly perturbed Robin boundary value problem, and investigate the midpoint upwind scheme on the Shishkin mesh and the Bakhvalov-Shishkin mesh to achieve better uniform convergence. The elaborate ε-uniform pointwise estimates are proved by using the comparison principle and barrier functions. The numerical experiments support the theoretical results for the schemes on the meshes.
文摘The upwind scheme is very important in the numerical approximation of some problems such as the convection dominated problem, the two-phase flow problem, and so on. For the fractional flow formulation of the two-phase flow problem, the Penalty Discontinuous Galerkin (PDG) methods combined with the upwind scheme are usually used to solve the phase pressure equation. In this case, unless the upwind scheme is taken into consideration in the velocity reconstruction, the local mass balance cannot hold exactly. In this paper, we present a scheme of velocity reconstruction in some H(div) spaces with considering the upwind scheme totally. Furthermore, the different ways to calculate the nonlinear coefficients may have distinct and significant effects, which have been investigated by some authors. We propose a new algorithm to obtain a more effective and stable approximation of the coefficients under the consideration of the upwind scheme.
文摘A parallelized upwind flux splitting scheme for supersonic reacting flows on hybrid meshes is presented. The complexity of super/hyper-sonic combustion flows makes it necessary to establish solvers with higher resolution and efficiency for multi-component Euler/N-S equations. Hence, a spatial second-order van Leer type flux vector splitting scheme is established by introducing auxiliary points in interpolation, and a domain decomposition method used on unstructured hybrid meshes for obtaining high calculating efficiency. The numerical scheme with five-stage Runge-Kutta time step method is implemented to the simulation of combustion flows, including the supersonic hydrogen/air combustion and the normal injection of hydrogen into reacting flows. Satisfying results are obtained compared with limited references.
基金supported in part by the Knowledge Innovation Project of the Chinese Academy of Sciences Nos. K5501312S1 and K5502212F1, and NSFC grant No. 10601062supported in part by NSF grant Nos. DMS-0305081 and DMS-0608720, NSFC grant No. 10228101 and NSAF grant No. 10676017
文摘We study the L^l-error estimates for the upwind scheme to the linear advection equations with a piecewise constant coefficients modeling linear waves crossing interfaces. Here the interface condition is immersed into the upwind scheme. We prove that, for initial data with a bounded variation, the numerical solution of the immersed interface upwind scheme converges in L^l-norm to the differential equation with the corresponding interface condition. We derive the one-halfth order L^l-error bounds with explicit coefficients following a technique used in [25]. We also use some inequalities on binomial coefficients proved in a consecutive paper [32].
基金supported by the National Basic Research Program of China("973" Project)(Grant No.2009CB724104)
文摘Nowadays,the upwind schemes are in a rapid development to capture shock accurately.However,these upwind schemes’properties at low speeds,such as their reconstruction scheme dependencies,grid dependencies,and Mach number dependencies,are concerned by few people.In this paper,a systematic study on their low speeds’issues is conducted.Through a series of tests,we can find that most parameter-free upwind schemes,widely used in practice today,are not applicable to low speeds’simulations.In contrast,SLAU and SLAU2 can give reliable results.Also,the upwind scheme’s influence on the accuracy is stronger than the reconstruction scheme’s influence at low speeds.
文摘Aim To construct a third order upwind scheme for convection equation. Methods Upwind Lagrange interpolation was used. Results and Conclusion The schemes L p stability for p∈ is proved. Numerical examples show that performance of the third order upwind scheme is better than that of most second order schemes.
基金This work is sponsored by the National Science Foundation of China under the Grant Number 11272037 and 91541105is also partially supported by National Key Basic Research and Development 973 Program of China under Grant Number 2014CB744100.
文摘The importance of eliminating errors in grid-metric evaluation for highorder difference schemes has been widely recognized in recent years,and it is known from the proof by Vinokur and Yee(NASA TM 209598,2000)that when conservative derivations of grid metric are used by Thomas,Lombard and Neier(AIAA J.,1978,17(10)and J.Spacecraft and rocket,1990,27(2)),errors caused by metric evaluation could be eliminated by linear schemes when flux splitting is not considered.According to the above achievement,central schemes without the use of flux splitting could fulfill the requirement of error elimination.Difficulties will arise for upwind schemes to attain the objective when the splitting is considered.In this study,further investigations are made on three aspects:Firstly,an idea of central scheme decomposition is introduced,and the procedure to derive the central scheme is proposed to evaluate grid metrics only.Secondly,the analysis has been made on the requirement of flux splitting to acquire free-stream preservation,and a Lax-Friedrichs-type splitting scheme is proposed as an example.Discussions about current study with that by Nonomura et al.(Computers and Fluids,2015,107)have been made.Thirdly,for halfnode-or mixed-type schemes,interpolations should be used to derive variables at half nodes.The requirement to achieve metric identity on this situation is analyzed and an idea of directionally consistent interpolation is proposed,which is manifested to be indispensable to avoid violations of metric identity and to eliminate metric-caused errors thereafter.Two numerical problems are tested,i.e.,the free-stream and vortex preservation onwavy,largely randomized and triangular-like grids.Numerical results validate aforementioned theoretical outcomes.
基金This work is supported by the National Natural Science Foundation of China(Grant No.11061021)Key Project of Chinese Ministry of Education(12024)+1 种基金the Scientific Research Projection of Higher Schools of Inner Mongolia(NJ10016,NJ10006,NJZZ12011)the National Natural Science Foundation of Inner Mongolia Province(2011BS0102).
文摘A bounded high order upwind scheme is presented for the modified Burgers’equation by using the normalized-variable formulation in the finite volume framework.The characteristic line of the present scheme in the normalizedvariable diagram is designed on the Hermite polynomial interpolation.In order to suppress unphysical oscillations,the present scheme respects both the TVD(total variational diminishing)constraint and CBC(convection boundedness criterion)condition.Numerical results demonstrate the present scheme possesses good robustness and high resolution for the modified Burgers’equation.
文摘In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.
文摘WT5,5”BX] A new class of numerical schemes is proposed to solve convection diffusion equations by combining the upwind technique and the method of operator splitting. For every time step, the multi dimensional approximation is performed in several independent directions alternatively, while the upwind technique is applied to treat the convection term in every individual direction. This scheme possesses maximum principle. Stability and convergence are analysed by energy method.[WT5,5”HZ]
基金supported by Postdoctoral Research Foundation of China(Grant No.2022M721054)the Natural Science Foundation of Henan Province of China(Grant No.222300420213).
文摘An upwind weak Galerkin finite element scheme was devised and analyzed in this article for convection-dominated Oseen equations. The numerical algorithm was based on the weak Galerkin method enhanced by upwind stabilization. The resulting finite element scheme uses equal-order, say k, polynomial spaces on each element for the velocity and the pressure unknowns. With finite elements of order k≥1, the numerical solutions are proved to converge at the rate of O(h^(k+1/2)) in an energy-like norm for convection-dominated Oseen equations. Numerical results are presented to demonstrate the accuracy and effectiveness of the upwind weak Galerkin scheme.
基金the National Key Research&Development Program of China(2016YFB020071)the National Natural Science Foundation of China(Grants 11532016 and 11702315).
文摘The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the bow shock.To improve the accuracy of aero-heating predictions,three multi-dimensional approaches on unstructured grids are coupled in our Reynolds-averaged Navier-Stokes(RANS)solver,including multi-dimensional upwind flux reconstruction(MUP),multi-dimensional limiter(MLP-u2)and multi-dimensional gradient reconstruction(MLR).The coupled multi-dimensional RANS solver is validated by several typical verification and validation(V&V)cases,including hypersonic flows over a cylinder,a blunt biconic,and a double-ellipsoid,with commonly used prism/tetrahedral hybrid grids.Finally,the coupled multi-dimensional solver is applied to simulating the heat flux distribution over a 3D engineering configuration,i.e.a Hermes-like space shuttle model.The obtained numerical results are compared with experimental data.The predicted results demonstrate that the coupled multi-dimensional approach has a good prediction capability on aerodynamic heating over a wide range of complex engineering configurations.
基金supported by the Department of Science & Technology, Government of India under research grant SR/S4/MS:318/06.
文摘Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.
文摘A fully discrete energy stability analysis is carried out for linear advection-diffusion problems discretized by generalized upwind summation-by-parts(upwind gSBP)schemes in space and implicit-explicit Runge-Kutta(IMEX-RK)schemes in time.Hereby,advection terms are discretized explicitly,while diffusion terms are solved implicitly.In this context,specific combinations of space and time discretizations enjoy enhanced stability properties.In fact,if the first-and second-derivative upwind gSBP operators fulfill a compatibility condition,the allowable time step size is independent of grid refinement,although the advective terms are discretized explicitly.In one space dimension it is shown that upwind gSBP schemes represent a general framework including standard discontinuous Galerkin(DG)schemes on a global level.While previous work for DG schemes has demonstrated that the combination of upwind advection fluxes and the central-type first Bassi-Rebay(BR1)scheme for diffusion does not allow for grid-independent stable time steps,the current work shows that central advection fluxes are compatible with BR1 regarding enhanced stability of IMEX time stepping.Furthermore,unlike previous discrete energy stability investigations for DG schemes,the present analysis is based on the discrete energy provided by the corresponding SBP norm matrix and yields time step restrictions independent of the discretization order in space,since no finite-element-type inverse constants are involved.Numerical experiments are provided confirming these theoretical findings.
基金NKBRSF CG 1990 3 2 80 5 National Natural Science F oundation of China !( No.5 98760 0 2 )
文摘In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the non physical oscillation by means of the group velocity control. The scheme is used to simulate the interactions of shock density stratified interface and the disturbed interface developing to vortex rollers. Numerical results are satisfactory.
