We theoretically investigate the high-order harmonic generation(HHG)of defect-free solids by solving the timedependent Schrodinger equation(TDSE).The results show that the harmonic intensity can be enhanced,harmonic o...We theoretically investigate the high-order harmonic generation(HHG)of defect-free solids by solving the timedependent Schrodinger equation(TDSE).The results show that the harmonic intensity can be enhanced,harmonic order can be extended,and modulation near the cutoff order becomes smaller for the second plateau by increasing the time delay.These effects are due to an increase of the electron population in higher energy bands,where the larger band gap allows electrons to release more energy,and the long electronic paths are suppressed.Additionally,we also investigate the HHG of defective solids by Bohmian trajectories(BT).It is found that the harmonic intensity of the second plateau can be further enhanced.Simultaneously,cutoff order is also extended due to Bohmian particles moving farther away from the defective zone.展开更多
A new oxidative N-heterocyclic carbene(NHC)-catalyzed high-order[7+3]annulation reaction ofγ-indolyl phenols as 1,7-dinucleophiles andα,β-alkynals with the aid of Sc(OTf)_(3)is reported,enabling the highly regiosel...A new oxidative N-heterocyclic carbene(NHC)-catalyzed high-order[7+3]annulation reaction ofγ-indolyl phenols as 1,7-dinucleophiles andα,β-alkynals with the aid of Sc(OTf)_(3)is reported,enabling the highly regioselective access to unprecedented polyarene-fused ten-membered lactams bearing a bridged aryl-aryl-indole scaffold in moderate to good yields.This protocol demonstrates a broad substrate scope,good compatibility with substituents and complete regioselectivity,providing an organocatalytic modular synthetic strategy for creating medium-sized lactams.展开更多
We investigate theoretically the effects of chirped laser pulses on high-order harmonic generation(HHG)from solids.We find that the harmonic spectra display redshifts for the driving laser pulses with negative chirp a...We investigate theoretically the effects of chirped laser pulses on high-order harmonic generation(HHG)from solids.We find that the harmonic spectra display redshifts for the driving laser pulses with negative chirp and blueshifts for those with positive chirp,which is due to the change in the instantaneous frequency of the driving laser for different chirped pulses.The analysis of crystal-momentum-resolved(k-resolved)HHG reveals that the frequency shifts are equal for the harmonics generated by different crystal momentum channels.The frequency shifts in the cutoff region are larger than those in the plateau region.With the increase of the absolute value of the chirp parameters,the frequency shifts of HHG become more significant,leading to the shifts from odd-to even-order harmonics.We also demonstrate that the frequency shifts of harmonic spectra are related to the duration of the chirped laser field,but are insensitive to the laser intensity and dephasing time.展开更多
Thermal expansion is crucial for various industrial processes and is increasingly the focus of research endeavors aimed at improving material performance.However,it is the continuous advancements in first-principles c...Thermal expansion is crucial for various industrial processes and is increasingly the focus of research endeavors aimed at improving material performance.However,it is the continuous advancements in first-principles calculations that have enabled researchers to understand the microscopic origins of thermal expansion.In this study,we propose a coefficient of thermal expansion(CTE)calculation scheme based on self-consistent phonon theory,incorporating the fourth-order anharmonicity.We selected four structures(Si,CaZrF_(6),SrTiO_(3),NaBr)to investigate high-order anharmonicity’s impact on their CTEs,based on bonding types.The results indicate that our method goes beyond the second-order quasi-harmonic approximation and the third-order perturbation theory,aligning closely with experimental data.Furthermore,we observed that an increase in the ionicity of the structures leads to a more pronounced influence of high-order anharmonicity on CTE,with this effect primarily manifesting in variations of the Grüneisen parameter.Our research provides a theoretical foundation for accurately predicting and regulating the thermal expansion behavior of materials.展开更多
An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into ...An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into an equivalent system,and the k-order backward differentiation formula(BDF k)and central difference formula are used to discretize the temporal and spatial derivatives,respectively.Different from the traditional discrete method that adopts full implicit or full explicit for the nonlinear integral terms,the proposed scheme is based on the SAV idea and can be treated semi-implicitly,taking into account both accuracy and effectiveness.Numerical results are presented to demonstrate the high-order convergence(up to fourth-order)of the developed schemes and it is computationally efficient in long-time computations.展开更多
We present a comprehensive study on the role of various excited states in high-order harmonic generation of hydrogen atoms driven by a long-wavelength(1500 nm)laser field.By numerically solving the time-dependent Schr...We present a comprehensive study on the role of various excited states in high-order harmonic generation of hydrogen atoms driven by a long-wavelength(1500 nm)laser field.By numerically solving the time-dependent Schrodinger equation(TDSE)and performing a time-frequency analysis,we investigate the influence of individual excited states on the harmonic spectrum.Our results reveal that the 2s excited state primarily contributes to the enhancement of high-energy harmonic yields by facilitating long electron trajectories,while the 2p excited state predominantly suppresses harmonic yields in the lower-energy region(20th-50th orders)by altering the contributions of electron trajectories.Our results highlight the critical role of the excited states in the HHG process,even at longer laser wavelengths.展开更多
We performed real-time and real-space numerical simulations of high-order harmonic generation in the threedimensional structured molecule methane(CH_(4)) using time-dependent density functional theory. By irradiating ...We performed real-time and real-space numerical simulations of high-order harmonic generation in the threedimensional structured molecule methane(CH_(4)) using time-dependent density functional theory. By irradiating the methane molecule with an elliptically polarized laser pulse polarized in the x–y plane, we observed significant even-order harmonic emission in the z-direction. By analyzing the electron dynamics in the electric field and the multi-orbital effects of the molecule, we revealed that electron recombination near specific atoms in methane is the primary source of highorder harmonic generation in the z-direction. Furthermore, we identified the dominant molecular orbitals responsible for the enhancement of harmonics in this direction and demonstrated the critical role played by multi-orbital effects in this process.展开更多
The high-speed development of space defense technology demands a high state estimation capacity for spacecraft tracking methods.However,reentry flight is accompanied by complex flight environments,which brings to the ...The high-speed development of space defense technology demands a high state estimation capacity for spacecraft tracking methods.However,reentry flight is accompanied by complex flight environments,which brings to the uncertain,complex,and strongly coupled non-Gaussian detection noise.As a result,there are several intractable considerations on the problem of state estimation tasks corrupted by complex non-Gaussian outliers for non-linear dynamics systems in practical application.To address these issues,a new iterated rational quadratic(RQ)kernel high-order unscented Kalman filtering(IRQHUKF)algorithm via capturing the statistics to break through the limitations of the Gaussian assumption is proposed.Firstly,the characteristic analysis of the RQ kernel is investigated in detail,which is the first attempt to carry out an exploration of the heavy-tailed characteristic and the ability on capturing highorder moments of the RQ kernel.Subsequently,the RQ kernel method is first introduced into the UKF algorithm as an error optimization criterion,termed the iterated RQ kernel-UKF(RQ-UKF)algorithm by derived analytically,which not only retains the high-order moments propagation process but also enhances the approximation capacity in the non-Gaussian noise problem for its ability in capturing highorder moments and heavy-tailed characteristics.Meanwhile,to tackle the limitations of the Gaussian distribution assumption in the linearization process of the non-linear systems,the high-order Sigma Points(SP)as a subsidiary role in propagating the state high-order statistics is devised by the moments matching method to improve the RQ-UKF.Finally,to further improve the flexibility of the IRQ-HUKF algorithm in practical application,an adaptive kernel parameter is derived analytically grounded in the Kullback-Leibler divergence(KLD)method and parametric sensitivity analysis of the RQ kernel.The simulation results demonstrate that the novel IRQ-HUKF algorithm is more robust and outperforms the existing advanced UKF with respect to the kernel method in reentry vehicle tracking scenarios under various noise environments.展开更多
With the availability of high-performance computing technology and the development of advanced numerical simulation methods, Computational Fluid Dynamics (CFD) is becoming more and more practical and efficient in engi...With the availability of high-performance computing technology and the development of advanced numerical simulation methods, Computational Fluid Dynamics (CFD) is becoming more and more practical and efficient in engineering. As one of the high-precision representative algorithms, the high-order Discontinuous Galerkin Method (DGM) has not only attracted widespread attention from scholars in the CFD research community, but also received strong development. However, when DGM is extended to high-speed aerodynamic flow field calculations, non-physical numerical Gibbs oscillations near shock waves often significantly affect the numerical accuracy and even cause calculation failure. Data driven approaches based on machine learning techniques can be used to learn the characteristics of Gibbs noise, which motivates us to use it in high-speed DG applications. To achieve this goal, labeled data need to be generated in order to train the machine learning models. This paper proposes a new method for denoising modeling of Gibbs phenomenon using a machine learning technique, the zero-shot learning strategy, to eliminate acquiring large amounts of CFD data. The model adopts a graph convolutional network combined with graph attention mechanism to learn the denoising paradigm from synthetic Gibbs noise data and generalize to DGM numerical simulation data. Numerical simulation results show that the Gibbs denoising model proposed in this paper can suppress the numerical oscillation near shock waves in the high-order DGM. Our work automates the extension of DGM to high-speed aerodynamic flow field calculations with higher generalization and lower cost.展开更多
The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass cons...The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.展开更多
In this paper,a high-order scheme based on the lattice Boltzmann flux solver(LBFS)is proposed to simulate viscous compressible flows.The flux reconstruction(FR)approach is adopted to implement the spatial discretizati...In this paper,a high-order scheme based on the lattice Boltzmann flux solver(LBFS)is proposed to simulate viscous compressible flows.The flux reconstruction(FR)approach is adopted to implement the spatial discretization.The LBFS is employed to compute the inviscid flux by using the local reconstruction of the lattice Boltzmann equation solutions from macroscopic flow variables.Meanwhile,a switch function is used in LBFS to adjust the magnitude of the numerical viscosity.Thus,it is more beneficial to capture both strong shock waves and thin boundary layers.Moreover,the viscous flux is computed according to the local discontinuous Galerkin method.Some typical compressible viscous problems,including manufactured solution case,lid-driven cavity flow,supersonic flow around a cylinder and subsonic flow over a NACA0012 airfoil,are simulated to demonstrate the accuracy and robustness of the proposed FR-LBFS.展开更多
The lattice Boltzmann method （LBM） is coupled with the multiple-relaxation- time （MRT） collision model and the three-dimensional 19-discrete-velocity （D3Q19） model to resolve intermittent behaviors on small scal...The lattice Boltzmann method （LBM） is coupled with the multiple-relaxation- time （MRT） collision model and the three-dimensional 19-discrete-velocity （D3Q19） model to resolve intermittent behaviors on small scales in isotropic turbulent flows. The high- order scaling exponents of the velocity structure functions, the probability distribution functions of Lagrangian accelerations, and the local energy dissipation rates are investi- gated. The self-similarity of the space-time velocity structure functions is explored using the extended self-similarity （ESS） method, which was originally developed for velocity spatial structure functions. The scaling exponents of spatial structure functions at up to ten orders are consistent with the experimental measurements and theoretical results, implying that the LBM can accurately resolve the intermittent behaviors. This valida~ tion provides a solid basis for using the LBM to study more complex processes that are sensitive to small scales in turbulent flows, such as the relative dispersion of pollutants and mesoscale structures of preferential concentration of heavy particles suspended in turbulent flows.展开更多
In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solut...In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.展开更多
基金supported by the Natural Science Foundation of Jilin Province of China(Grant No.20230101014JC)the Fundamental Research Funds for the Central Universities(Grant No.2572021BC05)the National Natural Science Foundation of China(Grant No.12374265)。
文摘We theoretically investigate the high-order harmonic generation(HHG)of defect-free solids by solving the timedependent Schrodinger equation(TDSE).The results show that the harmonic intensity can be enhanced,harmonic order can be extended,and modulation near the cutoff order becomes smaller for the second plateau by increasing the time delay.These effects are due to an increase of the electron population in higher energy bands,where the larger band gap allows electrons to release more energy,and the long electronic paths are suppressed.Additionally,we also investigate the HHG of defective solids by Bohmian trajectories(BT).It is found that the harmonic intensity of the second plateau can be further enhanced.Simultaneously,cutoff order is also extended due to Bohmian particles moving farther away from the defective zone.
基金National Natural Science Foundation of China(Nos.21971090 and 22271123)the NSF of Jiangsu Province(No.BK20230201)+1 种基金the Natural Science Foundation of Jiangsu Education Committee(No.22KJB150024)the Natural Science Foundation of Jiangsu Normal University(No.21XSRX010)。
文摘A new oxidative N-heterocyclic carbene(NHC)-catalyzed high-order[7+3]annulation reaction ofγ-indolyl phenols as 1,7-dinucleophiles andα,β-alkynals with the aid of Sc(OTf)_(3)is reported,enabling the highly regioselective access to unprecedented polyarene-fused ten-membered lactams bearing a bridged aryl-aryl-indole scaffold in moderate to good yields.This protocol demonstrates a broad substrate scope,good compatibility with substituents and complete regioselectivity,providing an organocatalytic modular synthetic strategy for creating medium-sized lactams.
基金Project supported by the Natural Science Foundation of Jilin Province of China(Grant No.20230101014JC)the National Natural Science Foundation of China(Grant No.12374265)。
文摘We investigate theoretically the effects of chirped laser pulses on high-order harmonic generation(HHG)from solids.We find that the harmonic spectra display redshifts for the driving laser pulses with negative chirp and blueshifts for those with positive chirp,which is due to the change in the instantaneous frequency of the driving laser for different chirped pulses.The analysis of crystal-momentum-resolved(k-resolved)HHG reveals that the frequency shifts are equal for the harmonics generated by different crystal momentum channels.The frequency shifts in the cutoff region are larger than those in the plateau region.With the increase of the absolute value of the chirp parameters,the frequency shifts of HHG become more significant,leading to the shifts from odd-to even-order harmonics.We also demonstrate that the frequency shifts of harmonic spectra are related to the duration of the chirped laser field,but are insensitive to the laser intensity and dephasing time.
基金Project supported by the National Natural Science Foundation of China(Grant No.62125402).
文摘Thermal expansion is crucial for various industrial processes and is increasingly the focus of research endeavors aimed at improving material performance.However,it is the continuous advancements in first-principles calculations that have enabled researchers to understand the microscopic origins of thermal expansion.In this study,we propose a coefficient of thermal expansion(CTE)calculation scheme based on self-consistent phonon theory,incorporating the fourth-order anharmonicity.We selected four structures(Si,CaZrF_(6),SrTiO_(3),NaBr)to investigate high-order anharmonicity’s impact on their CTEs,based on bonding types.The results indicate that our method goes beyond the second-order quasi-harmonic approximation and the third-order perturbation theory,aligning closely with experimental data.Furthermore,we observed that an increase in the ionicity of the structures leads to a more pronounced influence of high-order anharmonicity on CTE,with this effect primarily manifesting in variations of the Grüneisen parameter.Our research provides a theoretical foundation for accurately predicting and regulating the thermal expansion behavior of materials.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12001210 and 12261103)the Natural Science Foundation of Henan(Grant No.252300420308)the Yunnan Fundamental Research Projects(Grant No.202301AT070117).
文摘An efficient and accurate scalar auxiliary variable(SAV)scheme for numerically solving nonlinear parabolic integro-differential equation(PIDE)is developed in this paper.The original equation is first transformed into an equivalent system,and the k-order backward differentiation formula(BDF k)and central difference formula are used to discretize the temporal and spatial derivatives,respectively.Different from the traditional discrete method that adopts full implicit or full explicit for the nonlinear integral terms,the proposed scheme is based on the SAV idea and can be treated semi-implicitly,taking into account both accuracy and effectiveness.Numerical results are presented to demonstrate the high-order convergence(up to fourth-order)of the developed schemes and it is computationally efficient in long-time computations.
基金supported by the Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi。
文摘We present a comprehensive study on the role of various excited states in high-order harmonic generation of hydrogen atoms driven by a long-wavelength(1500 nm)laser field.By numerically solving the time-dependent Schrodinger equation(TDSE)and performing a time-frequency analysis,we investigate the influence of individual excited states on the harmonic spectrum.Our results reveal that the 2s excited state primarily contributes to the enhancement of high-energy harmonic yields by facilitating long electron trajectories,while the 2p excited state predominantly suppresses harmonic yields in the lower-energy region(20th-50th orders)by altering the contributions of electron trajectories.Our results highlight the critical role of the excited states in the HHG process,even at longer laser wavelengths.
基金Project supported by the National Natural Science Foundation of China (Grant No. 12204214)the National Key Research and Development Program of China (Grant No. 2022YFE0134200)+1 种基金the Fundamental Research Funds for the Central Universities (Grant No. GK202207012), QCYRCXM-2022-241the Guangdong Basic and Applied Basic Research Foundation (Grant No. 2025A1515011117)。
文摘We performed real-time and real-space numerical simulations of high-order harmonic generation in the threedimensional structured molecule methane(CH_(4)) using time-dependent density functional theory. By irradiating the methane molecule with an elliptically polarized laser pulse polarized in the x–y plane, we observed significant even-order harmonic emission in the z-direction. By analyzing the electron dynamics in the electric field and the multi-orbital effects of the molecule, we revealed that electron recombination near specific atoms in methane is the primary source of highorder harmonic generation in the z-direction. Furthermore, we identified the dominant molecular orbitals responsible for the enhancement of harmonics in this direction and demonstrated the critical role played by multi-orbital effects in this process.
基金supported by the National Natural Science Foundation of China under Grant No.12072090.
文摘The high-speed development of space defense technology demands a high state estimation capacity for spacecraft tracking methods.However,reentry flight is accompanied by complex flight environments,which brings to the uncertain,complex,and strongly coupled non-Gaussian detection noise.As a result,there are several intractable considerations on the problem of state estimation tasks corrupted by complex non-Gaussian outliers for non-linear dynamics systems in practical application.To address these issues,a new iterated rational quadratic(RQ)kernel high-order unscented Kalman filtering(IRQHUKF)algorithm via capturing the statistics to break through the limitations of the Gaussian assumption is proposed.Firstly,the characteristic analysis of the RQ kernel is investigated in detail,which is the first attempt to carry out an exploration of the heavy-tailed characteristic and the ability on capturing highorder moments of the RQ kernel.Subsequently,the RQ kernel method is first introduced into the UKF algorithm as an error optimization criterion,termed the iterated RQ kernel-UKF(RQ-UKF)algorithm by derived analytically,which not only retains the high-order moments propagation process but also enhances the approximation capacity in the non-Gaussian noise problem for its ability in capturing highorder moments and heavy-tailed characteristics.Meanwhile,to tackle the limitations of the Gaussian distribution assumption in the linearization process of the non-linear systems,the high-order Sigma Points(SP)as a subsidiary role in propagating the state high-order statistics is devised by the moments matching method to improve the RQ-UKF.Finally,to further improve the flexibility of the IRQ-HUKF algorithm in practical application,an adaptive kernel parameter is derived analytically grounded in the Kullback-Leibler divergence(KLD)method and parametric sensitivity analysis of the RQ kernel.The simulation results demonstrate that the novel IRQ-HUKF algorithm is more robust and outperforms the existing advanced UKF with respect to the kernel method in reentry vehicle tracking scenarios under various noise environments.
基金co-supported by the Aeronautical Science Foundation of China(Nos.2018ZA52002,2019ZA052011).
文摘With the availability of high-performance computing technology and the development of advanced numerical simulation methods, Computational Fluid Dynamics (CFD) is becoming more and more practical and efficient in engineering. As one of the high-precision representative algorithms, the high-order Discontinuous Galerkin Method (DGM) has not only attracted widespread attention from scholars in the CFD research community, but also received strong development. However, when DGM is extended to high-speed aerodynamic flow field calculations, non-physical numerical Gibbs oscillations near shock waves often significantly affect the numerical accuracy and even cause calculation failure. Data driven approaches based on machine learning techniques can be used to learn the characteristics of Gibbs noise, which motivates us to use it in high-speed DG applications. To achieve this goal, labeled data need to be generated in order to train the machine learning models. This paper proposes a new method for denoising modeling of Gibbs phenomenon using a machine learning technique, the zero-shot learning strategy, to eliminate acquiring large amounts of CFD data. The model adopts a graph convolutional network combined with graph attention mechanism to learn the denoising paradigm from synthetic Gibbs noise data and generalize to DGM numerical simulation data. Numerical simulation results show that the Gibbs denoising model proposed in this paper can suppress the numerical oscillation near shock waves in the high-order DGM. Our work automates the extension of DGM to high-speed aerodynamic flow field calculations with higher generalization and lower cost.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11862003 and 81860635)the Key Project of the Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2017GXNSFDA198038)+3 种基金the Project of Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2018GXNSFAA281302)the Project for Promotion of Young and Middle-aged Teachers’Basic Scientific Research Ability in Guangxi Universities,China(Grant No.2019KY0084)the“Bagui Scholar”Teams for Innovation and Research Project of Guangxi Zhuang Autonomous Region,Chinathe Graduate Innovation Program of Guangxi Normal University,China(Grant No.JXYJSKT-2019-007)。
文摘The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.
基金supported by the National Natural Science Foundation of China(No.12072158)the Natural Science Foundation of Jiangsu Province(No.BK20191271)+1 种基金the Research Fund of Key Laboratory of Computational AerodynamicsAVIC Aerodynamics Research Institute(No.YL2022XFX0402)。
文摘In this paper,a high-order scheme based on the lattice Boltzmann flux solver(LBFS)is proposed to simulate viscous compressible flows.The flux reconstruction(FR)approach is adopted to implement the spatial discretization.The LBFS is employed to compute the inviscid flux by using the local reconstruction of the lattice Boltzmann equation solutions from macroscopic flow variables.Meanwhile,a switch function is used in LBFS to adjust the magnitude of the numerical viscosity.Thus,it is more beneficial to capture both strong shock waves and thin boundary layers.Moreover,the viscous flux is computed according to the local discontinuous Galerkin method.Some typical compressible viscous problems,including manufactured solution case,lid-driven cavity flow,supersonic flow around a cylinder and subsonic flow over a NACA0012 airfoil,are simulated to demonstrate the accuracy and robustness of the proposed FR-LBFS.
基金Project supported by the Science Challenge Program(No.TZ2016001)the National Natural Science Foundation of China(Nos.11472277,11572331,11232011,and 11772337)+2 种基金the Strategic Priority Research Program,Chinese Academy of Sciences(CAS)(No.XDB22040104)the Key Research Program of Frontier Sciences,CAS(No.QYZDJ-SSW-SYS002)the National Basic Research Program of China(973 Program)(No.2013CB834100)
文摘The lattice Boltzmann method （LBM） is coupled with the multiple-relaxation- time （MRT） collision model and the three-dimensional 19-discrete-velocity （D3Q19） model to resolve intermittent behaviors on small scales in isotropic turbulent flows. The high- order scaling exponents of the velocity structure functions, the probability distribution functions of Lagrangian accelerations, and the local energy dissipation rates are investi- gated. The self-similarity of the space-time velocity structure functions is explored using the extended self-similarity （ESS） method, which was originally developed for velocity spatial structure functions. The scaling exponents of spatial structure functions at up to ten orders are consistent with the experimental measurements and theoretical results, implying that the LBM can accurately resolve the intermittent behaviors. This valida~ tion provides a solid basis for using the LBM to study more complex processes that are sensitive to small scales in turbulent flows, such as the relative dispersion of pollutants and mesoscale structures of preferential concentration of heavy particles suspended in turbulent flows.
基金supported by the National Natural Science Foundation of China(No.12002290)。
文摘In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.
