This paper presents a gas-kinetic theory based multidimensional high-order method for the compressible Naiver–Stokes solutions. In our previous study, a spatially and temporally dependent third-order flux scheme with...This paper presents a gas-kinetic theory based multidimensional high-order method for the compressible Naiver–Stokes solutions. In our previous study, a spatially and temporally dependent third-order flux scheme with the use of a third-order gas distribution function is employed.However, the third-order flux scheme is quite complicated and less robust than the second-order scheme. In order to reduce its complexity and improve its robustness, the secondorder flux scheme is adopted instead in this paper, while the temporal order of method is maintained by using a two stage temporal discretization. In addition, its CPU cost is relatively lower than the previous scheme. Several test cases in two and three dimensions, containing high Mach number compressible flows and low speed high Reynolds number laminar flows, are presented to demonstrate the method capacity.展开更多
A subspace expanding technique(SET) is proposed to efficiently discover and find all zeros of nonlinear functions in multi-degree-of-freedom(MDOF) engineering systems by discretizing the space into smaller subdomains,...A subspace expanding technique(SET) is proposed to efficiently discover and find all zeros of nonlinear functions in multi-degree-of-freedom(MDOF) engineering systems by discretizing the space into smaller subdomains, which are called cells. The covering set of the cells is identified by parallel calculations with the root bracketing method. The covering set can be found first in a low-dimensional subspace, and then gradually extended to higher dimensional spaces with the introduction of more equations and variables into the calculations. The results show that the proposed SET is highlyefficient for finding zeros in high-dimensional spaces. The subdivision technique of the cell mapping method is further used to refine the covering set, and the obtained numerical results of zeros are accurate. Three examples are further carried out to verify the applicability of the proposed method, and very good results are achieved. It is believed that the proposed method will significantly enhance the ability to study the stability, bifurcation,and optimization problems in complex MDOF nonlinear dynamic systems.展开更多
We report the numerical observation of discrete spatial solitons in a periodically poled lithium niobate waveguide array by applying an electrical field through electro-optical effect. We show that discrete spatial so...We report the numerical observation of discrete spatial solitons in a periodically poled lithium niobate waveguide array by applying an electrical field through electro-optical effect. We show that discrete spatial soliton can be controlled by applied voltage in the periodically poled lithium niobate.展开更多
We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusi...We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.展开更多
In this article, the fmite element solution of quasi-three-dimensional (quasi-3-D) groundwater flow was mathematically analyzed. The research shows that the spurious oscillation solution to the Finite Element Model ...In this article, the fmite element solution of quasi-three-dimensional (quasi-3-D) groundwater flow was mathematically analyzed. The research shows that the spurious oscillation solution to the Finite Element Model (FEM) is the results choosing the small time step △t or the large element size L and using the non-diagonal storage matrix. The mechanism for this phenomenon is explained by the negative weighting factor of implicit part in the discretized equations. To avoid spurious oscillation solution, the criteria on the selection of △t and L for quasi-3-D groundwater flow simulations were identified. An application example of quasi-3-D groundwater flow simulation was presented to verify the criteria. The results indicate that temporal discretization scale has significant impact on the spurious oscillations in the finite-element solutions, and the spurious oscillations can be avoided in solving practical quasi-3-D groundwater flow problems if the criteria are satisfied.展开更多
Focuses on a study which discreted Ginzburg-Landau-BBM equations with periodic initial boundary value conditions by the finite difference method in spatial direction. Background on the discretization of the equations ...Focuses on a study which discreted Ginzburg-Landau-BBM equations with periodic initial boundary value conditions by the finite difference method in spatial direction. Background on the discretization of the equations and the priori estimates; Existence of the attractors for the discrete system; Estimates of the upper bounds of Hausdorff and fractal dimensions for the attractors.展开更多
Global linear instability analysis is a powerful tool for the complex flow diagnosis.However,the methods used in the past would generally suffer from some dis-advantages,either the excessive computational resources fo...Global linear instability analysis is a powerful tool for the complex flow diagnosis.However,the methods used in the past would generally suffer from some dis-advantages,either the excessive computational resources for the low-order methods or the tedious mathematical derivations for the high-order methods.The present work proposed a CFD-aided Galerkin methodology which combines the merits from both the low-order and high-order methods,where the expansion on proper basis func-tions is preserved to ensure a small matrix size,while the differentials,incompressibility constraints and boundary conditions are realized by applying the low-order linearized Navier-Stokes equation solvers on the basis functions on a fine grid.Several test cases have shown that the new method can get satisfactory results for one-dimensional,two-dimensional and three-dimensional flow problems and also for the problems with complex geometries and boundary conditions.展开更多
文摘This paper presents a gas-kinetic theory based multidimensional high-order method for the compressible Naiver–Stokes solutions. In our previous study, a spatially and temporally dependent third-order flux scheme with the use of a third-order gas distribution function is employed.However, the third-order flux scheme is quite complicated and less robust than the second-order scheme. In order to reduce its complexity and improve its robustness, the secondorder flux scheme is adopted instead in this paper, while the temporal order of method is maintained by using a two stage temporal discretization. In addition, its CPU cost is relatively lower than the previous scheme. Several test cases in two and three dimensions, containing high Mach number compressible flows and low speed high Reynolds number laminar flows, are presented to demonstrate the method capacity.
基金the National Natural Science Foundation of China (Nos. 11702213,11772243,11572215,and 11332008)the Natural Science Foundation of Shaanxi Province of China(No. 2018JQ1061)。
文摘A subspace expanding technique(SET) is proposed to efficiently discover and find all zeros of nonlinear functions in multi-degree-of-freedom(MDOF) engineering systems by discretizing the space into smaller subdomains, which are called cells. The covering set of the cells is identified by parallel calculations with the root bracketing method. The covering set can be found first in a low-dimensional subspace, and then gradually extended to higher dimensional spaces with the introduction of more equations and variables into the calculations. The results show that the proposed SET is highlyefficient for finding zeros in high-dimensional spaces. The subdivision technique of the cell mapping method is further used to refine the covering set, and the obtained numerical results of zeros are accurate. Three examples are further carried out to verify the applicability of the proposed method, and very good results are achieved. It is believed that the proposed method will significantly enhance the ability to study the stability, bifurcation,and optimization problems in complex MDOF nonlinear dynamic systems.
基金This work was supported by the National Natural Science Foundation of China (No. 60007001)the Foundation for Development of Science and Technology of Shanghai (No. OOJC14027)
文摘We report the numerical observation of discrete spatial solitons in a periodically poled lithium niobate waveguide array by applying an electrical field through electro-optical effect. We show that discrete spatial soliton can be controlled by applied voltage in the periodically poled lithium niobate.
基金Open Access funding enabled and organized by Projekt DEAL.
文摘We develop error-control based time integration algorithms for compressible fluid dynam-ics(CFD)applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime.Focusing on discontinuous spectral element semidis-cretizations,we design new controllers for existing methods and for some new embedded Runge-Kutta pairs.We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice.We compare a wide range of error-control-based methods,along with the common approach in which step size con-trol is based on the Courant-Friedrichs-Lewy(CFL)number.The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances,while additionally providing control of the temporal error at tighter tolerances.The numerical examples include challenging industrial CFD applications.
文摘In this article, the fmite element solution of quasi-three-dimensional (quasi-3-D) groundwater flow was mathematically analyzed. The research shows that the spurious oscillation solution to the Finite Element Model (FEM) is the results choosing the small time step △t or the large element size L and using the non-diagonal storage matrix. The mechanism for this phenomenon is explained by the negative weighting factor of implicit part in the discretized equations. To avoid spurious oscillation solution, the criteria on the selection of △t and L for quasi-3-D groundwater flow simulations were identified. An application example of quasi-3-D groundwater flow simulation was presented to verify the criteria. The results indicate that temporal discretization scale has significant impact on the spurious oscillations in the finite-element solutions, and the spurious oscillations can be avoided in solving practical quasi-3-D groundwater flow problems if the criteria are satisfied.
基金Project supported by the National Natural Science Poundation of China (No. 19861004).
文摘Focuses on a study which discreted Ginzburg-Landau-BBM equations with periodic initial boundary value conditions by the finite difference method in spatial direction. Background on the discretization of the equations and the priori estimates; Existence of the attractors for the discrete system; Estimates of the upper bounds of Hausdorff and fractal dimensions for the attractors.
基金supported by the National Science Foundation of China(NSFC Grants No.11822208,No.11772297,No.91752201,No.91752202)。
文摘Global linear instability analysis is a powerful tool for the complex flow diagnosis.However,the methods used in the past would generally suffer from some dis-advantages,either the excessive computational resources for the low-order methods or the tedious mathematical derivations for the high-order methods.The present work proposed a CFD-aided Galerkin methodology which combines the merits from both the low-order and high-order methods,where the expansion on proper basis func-tions is preserved to ensure a small matrix size,while the differentials,incompressibility constraints and boundary conditions are realized by applying the low-order linearized Navier-Stokes equation solvers on the basis functions on a fine grid.Several test cases have shown that the new method can get satisfactory results for one-dimensional,two-dimensional and three-dimensional flow problems and also for the problems with complex geometries and boundary conditions.