A discontinuity-capturing scheme of finite element method(FEM)is proposed.The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unstea...A discontinuity-capturing scheme of finite element method(FEM)is proposed.The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unsteady flows,which exhibits the capability of capturing the shock waves and/or thin shear layers accurately in an unsteady viscous flow at high Reynolds number. In particular,a new testing variable,i.e.,the disturbed kinetic energy E,is suggested and used in the adaptive mesh computation,which is universally applicable to the capturing of both shock waves and shear layers in the inviscid flow and viscous flow at high Reynolds number.Based on several calculated examples,this approach has been proved to be effective and efficient for the calculations of compressible and incompressible flows.展开更多
In this paper,we present a discontinuity and cusp capturing physicsinformed neural network(PINN)to solve Stokes equations with a piecewiseconstant viscosity and singular force along an interface.We first reformulate t...In this paper,we present a discontinuity and cusp capturing physicsinformed neural network(PINN)to solve Stokes equations with a piecewiseconstant viscosity and singular force along an interface.We first reformulate the governing equations in each fluid domain separately and replace the singular force effect with the traction balance equation between solutions in two sides along the interface.Since the pressure is discontinuous and the velocity has discontinuous derivatives across the interface,we hereby use a network consisting of two fully-connected sub-networks that approximate the pressure and velocity,respectively.The two sub-networks share the same primary coordinate input arguments but with different augmented feature inputs.These two augmented inputs provide the interface information,so we assume that a level set function is given and its zero level set indicates the position of the interface.The pressure sub-network uses an indicator function as an augmented input to capture the function discontinuity,while the velocity sub-network uses a cusp-enforced level set function to capture the derivative discontinuities via the traction balance equation.We perform a series of numerical experiments to solve two-and three-dimensional Stokes interface problems and perform an accuracy comparison with the augmented immersed interface methods in literature.Our results indicate that even a shallow network with a moderate number of neurons and sufficient training data points can achieve prediction accuracy comparable to that of immersed interface methods.展开更多
Based on newly developed weight-based smoothness detectors and non-linear interpolations designed to capture discontinuities for the multiderivative com-bined dissipative compact scheme(MDCS),hybrid linear and nonline...Based on newly developed weight-based smoothness detectors and non-linear interpolations designed to capture discontinuities for the multiderivative com-bined dissipative compact scheme(MDCS),hybrid linear and nonlinear interpolations are proposed to form hybrid MDCS.These detectors are derived from the weights used for the nonlinear interpolations and can provide suitable switches between the linear and the nonlinear schemes to realize the characteristics for the hybrid MDCS of capturing discontinuities and maintaining high resolution in the region without large discontinuities.To save computational cost,the nonlinear scheme with characteris-tic decomposition is only applied in the detected discontinuities region by specially designed hybrid strategy.Typical tests show that the hybrid MDCS is capable of cap-turing discontinuities and maintaining high resolution power for the smooth region at the same time.With the satisfaction of the geometric conservative law(GCL),the MDCS is further applied on curvilinear mesh to present its promising capability of handling pragmatic simulations.展开更多
基金The project supported by the National Natural Science Foundation of China (10125210),the Hundred-Talent Programme of the Chinese Academy of Sciences and the Innovation Project of the Chinese Academy of Sciences (KJCX-SW-L04,KJCX2-SW-L2)
文摘A discontinuity-capturing scheme of finite element method(FEM)is proposed.The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unsteady flows,which exhibits the capability of capturing the shock waves and/or thin shear layers accurately in an unsteady viscous flow at high Reynolds number. In particular,a new testing variable,i.e.,the disturbed kinetic energy E,is suggested and used in the adaptive mesh computation,which is universally applicable to the capturing of both shock waves and shear layers in the inviscid flow and viscous flow at high Reynolds number.Based on several calculated examples,this approach has been proved to be effective and efficient for the calculations of compressible and incompressible flows.
基金supports by National Science and Technology Council,Taiwan,under research grants 111-2115-M-390-002 and 110-2115-M-A49-011-MY3,respectively.
文摘In this paper,we present a discontinuity and cusp capturing physicsinformed neural network(PINN)to solve Stokes equations with a piecewiseconstant viscosity and singular force along an interface.We first reformulate the governing equations in each fluid domain separately and replace the singular force effect with the traction balance equation between solutions in two sides along the interface.Since the pressure is discontinuous and the velocity has discontinuous derivatives across the interface,we hereby use a network consisting of two fully-connected sub-networks that approximate the pressure and velocity,respectively.The two sub-networks share the same primary coordinate input arguments but with different augmented feature inputs.These two augmented inputs provide the interface information,so we assume that a level set function is given and its zero level set indicates the position of the interface.The pressure sub-network uses an indicator function as an augmented input to capture the function discontinuity,while the velocity sub-network uses a cusp-enforced level set function to capture the derivative discontinuities via the traction balance equation.We perform a series of numerical experiments to solve two-and three-dimensional Stokes interface problems and perform an accuracy comparison with the augmented immersed interface methods in literature.Our results indicate that even a shallow network with a moderate number of neurons and sufficient training data points can achieve prediction accuracy comparable to that of immersed interface methods.
基金supported by the National Key Research and Development Plan(grant No.2016YFB0200700)the National Natural Science Foundation of China(grant Nos.11372342,11572342,and 11672321)the National Key Project GJXM92579.
文摘Based on newly developed weight-based smoothness detectors and non-linear interpolations designed to capture discontinuities for the multiderivative com-bined dissipative compact scheme(MDCS),hybrid linear and nonlinear interpolations are proposed to form hybrid MDCS.These detectors are derived from the weights used for the nonlinear interpolations and can provide suitable switches between the linear and the nonlinear schemes to realize the characteristics for the hybrid MDCS of capturing discontinuities and maintaining high resolution in the region without large discontinuities.To save computational cost,the nonlinear scheme with characteris-tic decomposition is only applied in the detected discontinuities region by specially designed hybrid strategy.Typical tests show that the hybrid MDCS is capable of cap-turing discontinuities and maintaining high resolution power for the smooth region at the same time.With the satisfaction of the geometric conservative law(GCL),the MDCS is further applied on curvilinear mesh to present its promising capability of handling pragmatic simulations.