This paper develops the high-order accurate entropy stable finite difference schemes for one-and two-dimensional special relativistic hydrodynamic equations.The schemes are built on the entropy conservative flux and t...This paper develops the high-order accurate entropy stable finite difference schemes for one-and two-dimensional special relativistic hydrodynamic equations.The schemes are built on the entropy conservative flux and the weighted essentially non-oscillatory(WENO)technique as well as explicit Runge-Kutta time discretization.The key is to technically construct the affordable entropy conservative flux of the semi-discrete second-order accurate entropy conservative schemes satisfying the semi-discrete entropy equality for the found convex entropy pair.As soon as the entropy conservative flux is derived,the dissipation term can be added to give the semidiscrete entropy stable schemes satisfying the semi-discrete entropy inequality with the given convex entropy function.The WENO reconstruction for the scaled entropy variables and the high-order explicit Runge-Kutta time discretization are implemented to obtain the fully-discrete high-order entropy stable schemes.Several numerical tests are conducted to validate the accuracy and the ability to capture discontinuities of our entropy stable schemes.展开更多
Accurate and efficient analysis of the coupled electroelastic behavior of piezoelectric structures is a challenging task in the community of computational mechanics.During the past few decades,the method of fundamenta...Accurate and efficient analysis of the coupled electroelastic behavior of piezoelectric structures is a challenging task in the community of computational mechanics.During the past few decades,the method of fundamental solutions(MFS)has emerged as a popular and well-established meshless boundary collocation method for the numerical solution of many engineering applications.The classical MFS formulation,however,leads to dense and non-symmetric coefficient matrices which will be computationally expensive for large-scale engineering simulations.In this paper,a localized version of the MFS(LMFS)is devised for electroelastic analysis of twodimensional(2D)piezoelectric structures.In the LMFS,the entire computational domain is divided into a set of overlapping small sub-domains where the MFS-based approximation and the moving least square(MLS)technique are employed.Different to the classical MFS,the LMFS will produce banded and sparse coefficient matrices which makes the method very attractive for large-scale simulations.Preliminary numerical experiments illustrate that the present LMFM is very promising for coupled electroelastic analysis of piezoelectric materials.展开更多
In this paper,based on the imaginary time gradient flow model in the density functional theory,a scalar auxiliary variable(SAV)method is developed for the ground state calculation of a given electronic structure syste...In this paper,based on the imaginary time gradient flow model in the density functional theory,a scalar auxiliary variable(SAV)method is developed for the ground state calculation of a given electronic structure system.To handle the orthonormality constraint on those wave functions,two kinds of penalty terms are introduced in designing the modified energy functional in SAV,i.e.,one for the norm preserving of each wave function,another for the orthogonality between each pair of different wave functions.A numerical method consisting of a designed scheme and a linear finite element method is used for the discretization.Theoretically,the desired unconditional decay of modified energy can be obtained from our method,while computationally,both the original energy and modified energy decay behaviors can be observed successfully from a number of numerical experiments.More importantly,numerical results show that the orthonormality among those wave functions can be automatically preserved,without explicitly preserving orthogonalization operations.This implies the potential of our method in large-scale simulations in density functional theory.展开更多
A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navie...A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navier-Stokes(N-S)equations and lattice Boltzmann equation(LBE).The macroscopic differential equations are discretized by the finite volume method,where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers.The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.LBFS is validated by its application to simulate the viscous decaying vortex flow,the driven cavity flow,the viscous flow past a circular cylinder,and the inviscid flow past a circular cylinder.The obtained numerical results compare very well with available data in the literature,which show that LBFS has the second order of accuracy in space,and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.展开更多
In Zhu,Wang and Gao(SIAM J.Sci.Comput.,43(2021),pp.A3009–A3031),we proposed a new framework of troubled-cell indicator(TCI)using K-means clustering and the numerical results demonstrate that it can detect the trouble...In Zhu,Wang and Gao(SIAM J.Sci.Comput.,43(2021),pp.A3009–A3031),we proposed a new framework of troubled-cell indicator(TCI)using K-means clustering and the numerical results demonstrate that it can detect the troubled cells accurately using the KXRCF indication variable.The main advantage of this TCI framework is its great potential of extensibility.In this follow-up work,we introduce three more indication variables,i.e.,the TVB,Fu-Shu and cell-boundary jump indication variables,and show their good performance by numerical tests to demonstrate that the TCI framework offers great flexibility in the choice of indication variables.We also compare the three indication variables with the KXRCF one,and the numerical results favor the KXRCF and the cell-boundary jump indication variables.展开更多
The energy gradient method has been proposed with the aim of better understanding the mechanism of flow transition from laminar flow to turbulent flow.In this method,it is demonstrated that the transition to turbulenc...The energy gradient method has been proposed with the aim of better understanding the mechanism of flow transition from laminar flow to turbulent flow.In this method,it is demonstrated that the transition to turbulence depends on the relative magnitudes of the transverse gradient of the total mechanical energy which amplifies the disturbance and the energy loss from viscous friction which damps the disturbance,for given imposed disturbance.For a given flow geometry and fluid properties,when the maximum of the function K(a function standing for the ratio of the gradient of total mechanical energy in the transverse direction to the rate of energy loss due to viscous friction in the streamwise direction)in the flow field is larger than a certain critical value,it is expected that instability would occur for some initial disturbances.In this paper,using the energy gradient analysis,the equation for calculating the energy gradient function K for plane Couette flow is derived.The result indicates that K reaches the maximum at the moving walls.Thus,the fluid layer near the moving wall is the most dangerous position to generate initial oscillation at sufficient high Re for given same level of normalized perturbation in the domain.The critical value of K at turbulent transition,which is observed from experiments,is about 370 for plane Couette flow when two walls move in opposite directions(anti-symmetry).This value is about the same as that for plane Poiseuille flow and pipe Poiseuille flow(385-389).Therefore,it is concluded that the critical value of K at turbulent transition is about 370-389 for wall-bounded parallel shear flows which include both pressure(symmetrical case)and shear driven flows(anti-symmetrical case).展开更多
The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying...The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying solutions of weakly singular Volterra integro-differential equations are sufficiently smooth.We provide a rigorous error analysis for the spectral methods,which shows that both the errors of approximate solutions and the errors of approximate derivatives of the solutions decay exponentially in L^(∞)-norm and weighted L^(2)-norm.The numerical examples are given to illustrate the theoretical results.展开更多
This paper proposes a high order deep domain decomposition method(HOrderDeepDDM)for solving high-frequency interface problems,which combines high order deep neural network(HOrderDNN)with domain decomposition method(DD...This paper proposes a high order deep domain decomposition method(HOrderDeepDDM)for solving high-frequency interface problems,which combines high order deep neural network(HOrderDNN)with domain decomposition method(DDM).The main idea of HOrderDeepDDM is to divide the computational domain into some sub-domains by DDM,and apply HOrderDNNs to solve the high-frequency problem on each sub-domain.Besides,we consider an adaptive learning rate annealing method to balance the errors inside the sub-domains,on the interface and the boundary during the optimization process.The performance of HOrderDeepDDM is evaluated on high-frequency elliptic and Helmholtz interface problems.The results indicate that:HOrderDeepDDM inherits the ability of DeepDDM to handle discontinuous interface problems and the power of HOrderDNN to approximate high-frequency problems.In detail,HOrderDeepDDMs(p>1)could capture the high-frequency information very well.When compared to the deep domain decomposition method(DeepDDM),HOrderDeepDDMs(p>1)converge faster and achieve much smaller relative errors with the same number of trainable parameters.For example,when solving the high-frequency interface elliptic problems in Section 3.3.1,the minimum relative errors obtained by HOrderDeepDDMs(p=9)are one order of magnitude smaller than that obtained by DeepDDMs when the number of the parameters keeps the same,as shown in Fig.4.展开更多
In this paper spectral Galerkin approximation of optimal control problem governed by fractional elliptic equation is investigated.To deal with the nonlocality of fractional Laplacian operator the Caffarelli-Silvestre ...In this paper spectral Galerkin approximation of optimal control problem governed by fractional elliptic equation is investigated.To deal with the nonlocality of fractional Laplacian operator the Caffarelli-Silvestre extension is utilized.The first order optimality condition of the extended optimal control problem is derived.A spectral Galerkin discrete scheme for the extended problem based on weighted Laguerre polynomials is developed.A priori error estimates for the spectral Galerkin discrete scheme is proved.Numerical experiments are presented to show the effectiveness of our methods and to verify the theoretical findings.展开更多
Recently,the use of porous materials has grown widely in many structures,such as beams,plates,and shells.The characteristics of porous materials change in the thickness direction by different functions.This study has ...Recently,the use of porous materials has grown widely in many structures,such as beams,plates,and shells.The characteristics of porous materials change in the thickness direction by different functions.This study has investigated the free vibration analysis of a sandwich porous elliptical micro-shell with a magneto-rheological fluid(MRF)core for the first time.Initially,we examined the displacement of the middle layer’s macro-and micro-components,using Love’s shell theory.Next,we used the modified couple stress theory(MCST)to obtain the strain and symmetrical curvature tensors for the three layers.The Hamilton’s principle was implemented to derive the equations of motion.We also used the Galerkin’s method to solve the equations of motion,resulting in a system of equations in the form of a linear eigenvalue problem.By solving the governing equations,we obtained the various natural frequencies and loss factors of the elliptical micro-shell,and compared them with the results in earlier studies.Lastly,we investigated the effects of thickness,porosity distribution pattern,aspect ratio,length scale parameter,and magnetic field intensity on the natural frequency and loss factor of the micro-shell.The data accuracy was validated by comparing them with those of reputable previous articles.展开更多
We propose a numerical method to solve the Monge-Ampère equation which admits a classical convex solution.The Monge-Ampère equation is reformulated into an equivalent first-order system.We adopt a novel reco...We propose a numerical method to solve the Monge-Ampère equation which admits a classical convex solution.The Monge-Ampère equation is reformulated into an equivalent first-order system.We adopt a novel reconstructed discontinuous approximation space which consists of piecewise irrotational polynomials.This space allows us to solve the first-order system in two sequential steps.In the first step,we solve a nonlinear system to obtain the approximation to the gradient.A Newton iteration is adopted to handle the nonlinearity of the system.The approximation to the primitive variable is obtained from the approximate gradient by a trivial least squares finite element method in the second step.Numerical examples in both two and three dimensions are presented to show an optimal convergence rate in accuracy.It is interesting to observe that the approximation solution is piecewise convex.Particularly,with the reconstructed approximation space,the proposed method numerically demonstrates a remarkable robustness.The convergence of the Newton iteration does not rely on the initial values.The dependence of the convergence on the penalty parameter in the discretization is also negligible,in comparison to the classical discontinuous approximation space.展开更多
In this paper,we propose a variational multiscale method(VMM)for the stationary incompressible magnetohydrodynamics equations.This method is defined by large-scale spaces for the velocity field and the magnetic field,...In this paper,we propose a variational multiscale method(VMM)for the stationary incompressible magnetohydrodynamics equations.This method is defined by large-scale spaces for the velocity field and the magnetic field,which aims to solve flows at high Reynolds numbers.We provide a new VMM formulation and prove its stability and convergence.Finally,some numerical experiments are presented to indicate the optimal convergence of our method.展开更多
The Richtmyer-Meshkov instability of interfaces separating elastic-plastic materials from vacuum is investigated by numerical simulation using a multi-material solid mechanics algorithm based on an Eulerian framework....The Richtmyer-Meshkov instability of interfaces separating elastic-plastic materials from vacuum is investigated by numerical simulation using a multi-material solid mechanics algorithm based on an Eulerian framework.The research efforts are directed to reveal the influence of the initial perturbation and material strength on the deformation of the perturbed interface impacted by an initial shock.By varying the initial amplitude(kx0)of the perturbed interface and the yield stress(sY),three typical modes of interface deformation have been identified as the broken mode,the stable mode and the oscillating mode.For the broken mode,the interface width(i.e.,the bubble position with respect to that of the spike)increases continuously resulting in a final separation of the spike from the perturbed interface.For the stable mode,the interface width grows to saturation and then maintains a nearly constant value in the long term.For the oscillating mode,the wavy-like interface moving forward obtains an aperiodic oscillation of small amplitude,namely,the interface width varies in time slightly around zero.The intriguing difference of the typical modes is interpreted qualitatively by comparing the early-stage wave motion and the commensurate pressure and effective stress.Further,the subsequent interface deformation is illustrated quantitatively via the time series of the interface positions and velocities of these three typical modes.展开更多
The background numerical noise#0 is determined by the maximum of truncation error and round-off error.For a chaotic system,the numerical error#(t)grows exponentially,say,#(t)=#0exp(kt),where k>0 is the so-called no...The background numerical noise#0 is determined by the maximum of truncation error and round-off error.For a chaotic system,the numerical error#(t)grows exponentially,say,#(t)=#0exp(kt),where k>0 is the so-called noise-growing exponent.This is the reason why one can not gain a convergent simulation of chaotic systems in a long enough interval of time by means of traditional algorithms in double precision,since the background numerical noise#0 might stop decreasing because of the use of double precision.This restriction can be overcome by means of the clean numerical simulation(CNS),which can decrease the background numerical noise#0 to any required tiny level.A lot of successful applications show the novelty and validity of the CNS.In this paper,we further propose some strategies to greatly increase the computational efficiency of the CNS algorithms for chaotic dynamical systems.It is highly suggested to keep a balance between truncation error and round-off error and besides to progressively enlarge the background numerical noise#0,since the exponentially increasing numerical noise#(t)is much larger than it.Some examples are given to illustrate the validity of our strategies for the CNS.展开更多
This investigation focused on the influence of the radial inertia effect on the propagation behavior of stress waves in thin-walled tubes subjected to combined longitudinal and torsional impact loads.Generalized chara...This investigation focused on the influence of the radial inertia effect on the propagation behavior of stress waves in thin-walled tubes subjected to combined longitudinal and torsional impact loads.Generalized characteristics theory was used to analyze the main features of the characteristic wave speeds and simple wave solutions in thin-walled tubes.The incremental elastic-plastic constitutive relations described by the rate-independent plasticity were adopted,and the finite difference method was used to investigate the evolution and propagation behaviors of combined elasticplastic stress waves in thin-walled tubes when the radial inertial effect was considered.The numerical results were compared with those obtained when the radial inertia effect was not considered.The results showed that the speed of the coupled stress wave increased when the radial inertia effect was considered.The hardening modulus of the material in the plastic stage had a greater impact on the coupled slow waves than on the coupled fast waves.展开更多
In this paper,we develop a direct-forcing immersed boundary projection method for simulating the dynamics in thermal fluid-solid interaction problems.The underlying idea of the method is that we treat the solid as mad...In this paper,we develop a direct-forcing immersed boundary projection method for simulating the dynamics in thermal fluid-solid interaction problems.The underlying idea of the method is that we treat the solid as made of fluid and introduce two virtual forcing terms.First,a virtual fluid force distributed only on the solid region is appended to the momentum equation to make the region behave like a real solid body and satisfy the prescribed velocity.Second,a virtual heat source located inside the solid region near the boundary is added to the energy transport equation to impose the thermal boundary condition on the solid boundary.We take the implicit second-order backward differentiation to discretize the time variable and employ the Choi-Moin projection scheme to split the coupled system.As for spatial discretization,second-order centered differences over a staggered Cartesian grid are used on the entire domain.The advantages of this method are its conceptual simplicity and ease of implementation.Numerical experiments are performed to demonstrate the high performance of the proposed method.Convergence tests show that the spatial convergence rates of all unknowns seem to be super-linear in the 1-norm and 2-norm while less than linear in the maximum norm.展开更多
This paper develops and analyzes interior penalty discontinuous Galerkin(IPDG)method by patch reconstruction technique for Helmholtz problems.The technique achieves high order approximation by locally solving a discre...This paper develops and analyzes interior penalty discontinuous Galerkin(IPDG)method by patch reconstruction technique for Helmholtz problems.The technique achieves high order approximation by locally solving a discrete least-squares over a neighboring element patch.We prove a prior error estimates in the L 2 norm and energy norm.For each fixed wave number k,the accuracy and efficiency of the method up to order five with high-order polynomials.Numerical examples are carried out to validate the theoretical results.展开更多
In this paper,the finite-time stability and instability are studied for nonlinear impulsive systems.There are mainly four concerns.1)For the system with stabilizing impulses,a Lyapunov theorem on global finite-time st...In this paper,the finite-time stability and instability are studied for nonlinear impulsive systems.There are mainly four concerns.1)For the system with stabilizing impulses,a Lyapunov theorem on global finite-time stability is presented.2)When the system without impulsive effects is globally finite-time stable(GFTS)and the settling time is continuous at the origin,it is proved that it is still GFTS over any class of impulse sequences,if the mixed impulsive jumps satisfy some mild conditions.3)For systems with destabilizing impulses,it is shown that to be finite-time stable,the destabilizing impulses should not occur too frequently,otherwise,the origin of the impulsive system is finite-time instable,which are formulated by average dwell time(ADT)conditions respectively.4)A theorem on finite-time instability is provided for system with stabilizing impulses.For each GFTS theorem of impulsive systems considered in this paper,the upper boundedness of settling time is given,which depends on the initial value and impulsive effects.Some numerical examples are given to illustrate the theoretical analysis.展开更多
In this paper,we propose an approach for constructing conservative and maximum-principle-preserving finite volume schemes by using the method of undetermined coefficients,which depend nonlinearly on the linear non-con...In this paper,we propose an approach for constructing conservative and maximum-principle-preserving finite volume schemes by using the method of undetermined coefficients,which depend nonlinearly on the linear non-conservative onesided fluxes.In order to facilitate the derivation of expressions of these undetermined coefficients,we explicitly provide a simple constriction condition with a scaling parameter.Such constriction conditions can ensure the final schemes are exact for linear solution problems and may induce various schemes by choosing different values for the parameter.In particular,when this parameter is taken to be 0,the nonlinear terms in our scheme degenerate to a harmonic average combination of the discrete linear fluxes,which has often been used in a variety of maximum-principle-preserving finite volume schemes.Thus our method of determining the coefficients of the nonlinear terms is more general.In addition,we prove the convergence of the proposed schemes by using a compactness technique.Numerical results demonstrate that our schemes can preserve the conservation property,satisfy the discrete maximum principle,possess a second-order accuracy,be exact for linear solution problems,and be available for anisotropic problems on distorted meshes.展开更多
Stochastic temperature distribution should be carefully inspected in the thermal-failure design of heterogeneous solids with unexpected random energy excitations.Stochastic multiscale modeling for these problems invol...Stochastic temperature distribution should be carefully inspected in the thermal-failure design of heterogeneous solids with unexpected random energy excitations.Stochastic multiscale modeling for these problems involve multiscale and highdimensional uncertain thermal parameters,which remains limitation of prohibitive computation.In this paper,we propose a multi-modes based constrained energy minimization generalized multiscale finite element method(MCEM-GMsFEM),which can transform the original stochastic multiscale model into a series of recursive multiscale models sharing the same deterministic material parameters by multiscale analysis.Thus,MCEM-GMsFEM reveals an inherent low-dimensional representation in random space,and is designed to effectively reduce the complexity of repeated computation of discretized multiscale systems.In addition,the convergence analysis is established,and the optimal error estimates are derived.Finally,several typical random fluctuations on multiscale thermal conductivity are considered to validate the theoretical results in the numerical examples.The numerical results indicate that the multi-modes multiscale approach is a robust integrated method with the excellent performance.展开更多
基金supported by the Special Project on High-performance Computing under the National Key R&D Program(No.2016YFB0200603)Science Challenge Project(No.TZ2016002)the National Natural Science Foundation of China(Nos.91630310 and 11421101),and High-Performance Computing Platform of Peking University.
文摘This paper develops the high-order accurate entropy stable finite difference schemes for one-and two-dimensional special relativistic hydrodynamic equations.The schemes are built on the entropy conservative flux and the weighted essentially non-oscillatory(WENO)technique as well as explicit Runge-Kutta time discretization.The key is to technically construct the affordable entropy conservative flux of the semi-discrete second-order accurate entropy conservative schemes satisfying the semi-discrete entropy equality for the found convex entropy pair.As soon as the entropy conservative flux is derived,the dissipation term can be added to give the semidiscrete entropy stable schemes satisfying the semi-discrete entropy inequality with the given convex entropy function.The WENO reconstruction for the scaled entropy variables and the high-order explicit Runge-Kutta time discretization are implemented to obtain the fully-discrete high-order entropy stable schemes.Several numerical tests are conducted to validate the accuracy and the ability to capture discontinuities of our entropy stable schemes.
基金supported by the National Natural Science Foundation of China(Nos.11872220,12111530006)the Natural Science Foundation of Shandong Province of China(Nos.ZR2021JQ02,2019KJI009)the Key Laboratory of Road Construction Technology and Equipment(Chang’an University,No.300102251505).
文摘Accurate and efficient analysis of the coupled electroelastic behavior of piezoelectric structures is a challenging task in the community of computational mechanics.During the past few decades,the method of fundamental solutions(MFS)has emerged as a popular and well-established meshless boundary collocation method for the numerical solution of many engineering applications.The classical MFS formulation,however,leads to dense and non-symmetric coefficient matrices which will be computationally expensive for large-scale engineering simulations.In this paper,a localized version of the MFS(LMFS)is devised for electroelastic analysis of twodimensional(2D)piezoelectric structures.In the LMFS,the entire computational domain is divided into a set of overlapping small sub-domains where the MFS-based approximation and the moving least square(MLS)technique are employed.Different to the classical MFS,the LMFS will produce banded and sparse coefficient matrices which makes the method very attractive for large-scale simulations.Preliminary numerical experiments illustrate that the present LMFM is very promising for coupled electroelastic analysis of piezoelectric materials.
基金The first author would like to thank the support from the UM-Funded PhD Assistantship from University of MacaoThe second author was partially supported by Macao Young Scholar Program(AM201919)+5 种基金excellent youth project of Hunan Education Department(19B543)Hunan National Applied Mathematics Center of Hunan Provincial Science and Technology Department(2020ZYT003)The third author would like to thank financial support from National Natural Science Foundation of China(Grant Nos.11922120,11871489)FDCT of Macao SAR(Grant No.0082/2020/A2)University of Macao(Grant No.MYRG2020-00265-FST)Guangdong-Hong Kong-Macao Joint Laboratory for Data-Driven Fluid Mechanics and Engineering Applications(Grant No.2020B1212030001).
文摘In this paper,based on the imaginary time gradient flow model in the density functional theory,a scalar auxiliary variable(SAV)method is developed for the ground state calculation of a given electronic structure system.To handle the orthonormality constraint on those wave functions,two kinds of penalty terms are introduced in designing the modified energy functional in SAV,i.e.,one for the norm preserving of each wave function,another for the orthogonality between each pair of different wave functions.A numerical method consisting of a designed scheme and a linear finite element method is used for the discretization.Theoretically,the desired unconditional decay of modified energy can be obtained from our method,while computationally,both the original energy and modified energy decay behaviors can be observed successfully from a number of numerical experiments.More importantly,numerical results show that the orthonormality among those wave functions can be automatically preserved,without explicitly preserving orthogonalization operations.This implies the potential of our method in large-scale simulations in density functional theory.
文摘A lattice Boltzmann flux solver(LBFS)is presented in this work for simulation of incompressible viscous and inviscid flows.The new solver is based on Chapman-Enskog expansion analysis,which is the bridge to link Navier-Stokes(N-S)equations and lattice Boltzmann equation(LBE).The macroscopic differential equations are discretized by the finite volume method,where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers.The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh,tie-up of mesh spacing and time interval,limitation to viscous flows.LBFS is validated by its application to simulate the viscous decaying vortex flow,the driven cavity flow,the viscous flow past a circular cylinder,and the inviscid flow past a circular cylinder.The obtained numerical results compare very well with available data in the literature,which show that LBFS has the second order of accuracy in space,and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.
基金We thank the anonymous reviewers and the editor for their valuable comments and suggestions.The research of Z.Gao is partially supported by the National Key R&D Program of China(No.2021YFF0704002)The four authors,Z.Wang,Z.Gao,H.Wang and H.Zhu,want to acknowledge the funding support by NSFC grant No.11871443+3 种基金The research of Z.Wang and H.Zhu is also partially sponsored by NUPTSF(Grant No.NY220040)Natural Science Foundation of Jiangsu Province of China(No.BK20191375)Postgraduate Research&Practice Innovation Program of Jiangsu Province under Grant No.KYCX200787The research of Q.Zhang is partially supported by NSFC grant No.12071214.
文摘In Zhu,Wang and Gao(SIAM J.Sci.Comput.,43(2021),pp.A3009–A3031),we proposed a new framework of troubled-cell indicator(TCI)using K-means clustering and the numerical results demonstrate that it can detect the troubled cells accurately using the KXRCF indication variable.The main advantage of this TCI framework is its great potential of extensibility.In this follow-up work,we introduce three more indication variables,i.e.,the TVB,Fu-Shu and cell-boundary jump indication variables,and show their good performance by numerical tests to demonstrate that the TCI framework offers great flexibility in the choice of indication variables.We also compare the three indication variables with the KXRCF one,and the numerical results favor the KXRCF and the cell-boundary jump indication variables.
文摘The energy gradient method has been proposed with the aim of better understanding the mechanism of flow transition from laminar flow to turbulent flow.In this method,it is demonstrated that the transition to turbulence depends on the relative magnitudes of the transverse gradient of the total mechanical energy which amplifies the disturbance and the energy loss from viscous friction which damps the disturbance,for given imposed disturbance.For a given flow geometry and fluid properties,when the maximum of the function K(a function standing for the ratio of the gradient of total mechanical energy in the transverse direction to the rate of energy loss due to viscous friction in the streamwise direction)in the flow field is larger than a certain critical value,it is expected that instability would occur for some initial disturbances.In this paper,using the energy gradient analysis,the equation for calculating the energy gradient function K for plane Couette flow is derived.The result indicates that K reaches the maximum at the moving walls.Thus,the fluid layer near the moving wall is the most dangerous position to generate initial oscillation at sufficient high Re for given same level of normalized perturbation in the domain.The critical value of K at turbulent transition,which is observed from experiments,is about 370 for plane Couette flow when two walls move in opposite directions(anti-symmetry).This value is about the same as that for plane Poiseuille flow and pipe Poiseuille flow(385-389).Therefore,it is concluded that the critical value of K at turbulent transition is about 370-389 for wall-bounded parallel shear flows which include both pressure(symmetrical case)and shear driven flows(anti-symmetrical case).
基金This work is supported by the Foundation for Talent Introduction of Guangdong Provincial University,Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)National Science Foundation of China(10971074).
文摘The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying solutions of weakly singular Volterra integro-differential equations are sufficiently smooth.We provide a rigorous error analysis for the spectral methods,which shows that both the errors of approximate solutions and the errors of approximate derivatives of the solutions decay exponentially in L^(∞)-norm and weighted L^(2)-norm.The numerical examples are given to illustrate the theoretical results.
基金supported partly by National Key R&D Program of China(grants Nos.2019YFA0709600 and 2019YFA0709602)National Natural Science Foundation of China(grants Nos.11831016 and 12101609)the Innovation Foundation of Qian Xuesen Laboratory of Space Technology。
文摘This paper proposes a high order deep domain decomposition method(HOrderDeepDDM)for solving high-frequency interface problems,which combines high order deep neural network(HOrderDNN)with domain decomposition method(DDM).The main idea of HOrderDeepDDM is to divide the computational domain into some sub-domains by DDM,and apply HOrderDNNs to solve the high-frequency problem on each sub-domain.Besides,we consider an adaptive learning rate annealing method to balance the errors inside the sub-domains,on the interface and the boundary during the optimization process.The performance of HOrderDeepDDM is evaluated on high-frequency elliptic and Helmholtz interface problems.The results indicate that:HOrderDeepDDM inherits the ability of DeepDDM to handle discontinuous interface problems and the power of HOrderDNN to approximate high-frequency problems.In detail,HOrderDeepDDMs(p>1)could capture the high-frequency information very well.When compared to the deep domain decomposition method(DeepDDM),HOrderDeepDDMs(p>1)converge faster and achieve much smaller relative errors with the same number of trainable parameters.For example,when solving the high-frequency interface elliptic problems in Section 3.3.1,the minimum relative errors obtained by HOrderDeepDDMs(p=9)are one order of magnitude smaller than that obtained by DeepDDMs when the number of the parameters keeps the same,as shown in Fig.4.
基金supported by the National Natural Science Foundation of China Project(Nos.12071402,11931003,12261131501,and 11971276)the Project of Scientific Research Fund of the Hunan Provincial Science and Technology Department(No.2022RC3022).
文摘In this paper spectral Galerkin approximation of optimal control problem governed by fractional elliptic equation is investigated.To deal with the nonlocality of fractional Laplacian operator the Caffarelli-Silvestre extension is utilized.The first order optimality condition of the extended optimal control problem is derived.A spectral Galerkin discrete scheme for the extended problem based on weighted Laguerre polynomials is developed.A priori error estimates for the spectral Galerkin discrete scheme is proved.Numerical experiments are presented to show the effectiveness of our methods and to verify the theoretical findings.
文摘Recently,the use of porous materials has grown widely in many structures,such as beams,plates,and shells.The characteristics of porous materials change in the thickness direction by different functions.This study has investigated the free vibration analysis of a sandwich porous elliptical micro-shell with a magneto-rheological fluid(MRF)core for the first time.Initially,we examined the displacement of the middle layer’s macro-and micro-components,using Love’s shell theory.Next,we used the modified couple stress theory(MCST)to obtain the strain and symmetrical curvature tensors for the three layers.The Hamilton’s principle was implemented to derive the equations of motion.We also used the Galerkin’s method to solve the equations of motion,resulting in a system of equations in the form of a linear eigenvalue problem.By solving the governing equations,we obtained the various natural frequencies and loss factors of the elliptical micro-shell,and compared them with the results in earlier studies.Lastly,we investigated the effects of thickness,porosity distribution pattern,aspect ratio,length scale parameter,and magnetic field intensity on the natural frequency and loss factor of the micro-shell.The data accuracy was validated by comparing them with those of reputable previous articles.
基金This research was supported by the National Natural Science Foundation in China(Nos.12201442,and 11971041).
文摘We propose a numerical method to solve the Monge-Ampère equation which admits a classical convex solution.The Monge-Ampère equation is reformulated into an equivalent first-order system.We adopt a novel reconstructed discontinuous approximation space which consists of piecewise irrotational polynomials.This space allows us to solve the first-order system in two sequential steps.In the first step,we solve a nonlinear system to obtain the approximation to the gradient.A Newton iteration is adopted to handle the nonlinearity of the system.The approximation to the primitive variable is obtained from the approximate gradient by a trivial least squares finite element method in the second step.Numerical examples in both two and three dimensions are presented to show an optimal convergence rate in accuracy.It is interesting to observe that the approximation solution is piecewise convex.Particularly,with the reconstructed approximation space,the proposed method numerically demonstrates a remarkable robustness.The convergence of the Newton iteration does not rely on the initial values.The dependence of the convergence on the penalty parameter in the discretization is also negligible,in comparison to the classical discontinuous approximation space.
基金supported by the National Natural Science Foundation of China(Nos.12071404,11971410,12261131501 and 12026254)Young Elite Scientist Sponsorship Program by CAST(No.2020QNRC001)+2 种基金Key Project of Scientific Research Project of Hunan Provincial Department of Education(No.22A0136)International Scientific and Technological Innovation Cooperation Base of Hunan Province for Computational Science(No.2018WK4006)Postgraduate Scientific Research Innovation Project of Hunan Province(No.CX20210612).
文摘In this paper,we propose a variational multiscale method(VMM)for the stationary incompressible magnetohydrodynamics equations.This method is defined by large-scale spaces for the velocity field and the magnetic field,which aims to solve flows at high Reynolds numbers.We provide a new VMM formulation and prove its stability and convergence.Finally,some numerical experiments are presented to indicate the optimal convergence of our method.
基金supported by the National Natural Science Foundation of China(Nos.12172353,92052301,11621202 and 12202436)the Science Challenge Project(No.TZ2016001)the National Science Foundation(No.CBET0755269).
文摘The Richtmyer-Meshkov instability of interfaces separating elastic-plastic materials from vacuum is investigated by numerical simulation using a multi-material solid mechanics algorithm based on an Eulerian framework.The research efforts are directed to reveal the influence of the initial perturbation and material strength on the deformation of the perturbed interface impacted by an initial shock.By varying the initial amplitude(kx0)of the perturbed interface and the yield stress(sY),three typical modes of interface deformation have been identified as the broken mode,the stable mode and the oscillating mode.For the broken mode,the interface width(i.e.,the bubble position with respect to that of the spike)increases continuously resulting in a final separation of the spike from the perturbed interface.For the stable mode,the interface width grows to saturation and then maintains a nearly constant value in the long term.For the oscillating mode,the wavy-like interface moving forward obtains an aperiodic oscillation of small amplitude,namely,the interface width varies in time slightly around zero.The intriguing difference of the typical modes is interpreted qualitatively by comparing the early-stage wave motion and the commensurate pressure and effective stress.Further,the subsequent interface deformation is illustrated quantitatively via the time series of the interface positions and velocities of these three typical modes.
基金supported by National Natural Science Foundation of China(No.12272230)Shanghai Pilot Program for Basic Research-Shanghai Jiao Tong University(No.21TQ1400202).
文摘The background numerical noise#0 is determined by the maximum of truncation error and round-off error.For a chaotic system,the numerical error#(t)grows exponentially,say,#(t)=#0exp(kt),where k>0 is the so-called noise-growing exponent.This is the reason why one can not gain a convergent simulation of chaotic systems in a long enough interval of time by means of traditional algorithms in double precision,since the background numerical noise#0 might stop decreasing because of the use of double precision.This restriction can be overcome by means of the clean numerical simulation(CNS),which can decrease the background numerical noise#0 to any required tiny level.A lot of successful applications show the novelty and validity of the CNS.In this paper,we further propose some strategies to greatly increase the computational efficiency of the CNS algorithms for chaotic dynamical systems.It is highly suggested to keep a balance between truncation error and round-off error and besides to progressively enlarge the background numerical noise#0,since the exponentially increasing numerical noise#(t)is much larger than it.Some examples are given to illustrate the validity of our strategies for the CNS.
文摘This investigation focused on the influence of the radial inertia effect on the propagation behavior of stress waves in thin-walled tubes subjected to combined longitudinal and torsional impact loads.Generalized characteristics theory was used to analyze the main features of the characteristic wave speeds and simple wave solutions in thin-walled tubes.The incremental elastic-plastic constitutive relations described by the rate-independent plasticity were adopted,and the finite difference method was used to investigate the evolution and propagation behaviors of combined elasticplastic stress waves in thin-walled tubes when the radial inertial effect was considered.The numerical results were compared with those obtained when the radial inertia effect was not considered.The results showed that the speed of the coupled stress wave increased when the radial inertia effect was considered.The hardening modulus of the material in the plastic stage had a greater impact on the coupled slow waves than on the coupled fast waves.
基金the Ministry of Science and Technology of Taiwan under grants MOST 107-2115-M-035-007-MY2(C.-S.You)MOST 106-2115-M-005-005-MY2(P.-W.Hsieh)MOST 106-2115-M-008-014-MY2(S.-Y.Yang)。
文摘In this paper,we develop a direct-forcing immersed boundary projection method for simulating the dynamics in thermal fluid-solid interaction problems.The underlying idea of the method is that we treat the solid as made of fluid and introduce two virtual forcing terms.First,a virtual fluid force distributed only on the solid region is appended to the momentum equation to make the region behave like a real solid body and satisfy the prescribed velocity.Second,a virtual heat source located inside the solid region near the boundary is added to the energy transport equation to impose the thermal boundary condition on the solid boundary.We take the implicit second-order backward differentiation to discretize the time variable and employ the Choi-Moin projection scheme to split the coupled system.As for spatial discretization,second-order centered differences over a staggered Cartesian grid are used on the entire domain.The advantages of this method are its conceptual simplicity and ease of implementation.Numerical experiments are performed to demonstrate the high performance of the proposed method.Convergence tests show that the spatial convergence rates of all unknowns seem to be super-linear in the 1-norm and 2-norm while less than linear in the maximum norm.
基金the National Key Research and Development Program of China(No.2020YFA0714200)the Science and Technology Major Project of Hubei Province under Grant No.2021AAA010+1 种基金the National Science Foundation of China(Nos.12125103 and 12071362)the Natural Science Foundation of Hubei Province(No.2019CFA007).
文摘This paper develops and analyzes interior penalty discontinuous Galerkin(IPDG)method by patch reconstruction technique for Helmholtz problems.The technique achieves high order approximation by locally solving a discrete least-squares over a neighboring element patch.We prove a prior error estimates in the L 2 norm and energy norm.For each fixed wave number k,the accuracy and efficiency of the method up to order five with high-order polynomials.Numerical examples are carried out to validate the theoretical results.
基金National Natural Science Foundation of China(No.61807017)the National Natural Science Foundation of China(Nos.12171122,11771128)+3 种基金Shenzhen Science and Technology Program(Grant No.RCJC20210609103755110)Fundamental Research Project of Shenzhen(No.JCYJ20190806143201649)Project(HIT.NSRIF.2020056)Supported by Natural Scientific Research Innovation Foundation in Harbin Institute of TechnologyResearch start-up fund Foundation in Harbin Institute of Technology(No.20190019)。
文摘In this paper,the finite-time stability and instability are studied for nonlinear impulsive systems.There are mainly four concerns.1)For the system with stabilizing impulses,a Lyapunov theorem on global finite-time stability is presented.2)When the system without impulsive effects is globally finite-time stable(GFTS)and the settling time is continuous at the origin,it is proved that it is still GFTS over any class of impulse sequences,if the mixed impulsive jumps satisfy some mild conditions.3)For systems with destabilizing impulses,it is shown that to be finite-time stable,the destabilizing impulses should not occur too frequently,otherwise,the origin of the impulsive system is finite-time instable,which are formulated by average dwell time(ADT)conditions respectively.4)A theorem on finite-time instability is provided for system with stabilizing impulses.For each GFTS theorem of impulsive systems considered in this paper,the upper boundedness of settling time is given,which depends on the initial value and impulsive effects.Some numerical examples are given to illustrate the theoretical analysis.
基金This work of J.Lv was partially supported by the Natural Science Foundation of Jilin province grant(No.20200201259JC)the National Natural Science Foundation of China grant(No.12271209)This work of Z.Sheng was partially supported by the National Natural Science Foundation of China grant(No.12071045).
文摘In this paper,we propose an approach for constructing conservative and maximum-principle-preserving finite volume schemes by using the method of undetermined coefficients,which depend nonlinearly on the linear non-conservative onesided fluxes.In order to facilitate the derivation of expressions of these undetermined coefficients,we explicitly provide a simple constriction condition with a scaling parameter.Such constriction conditions can ensure the final schemes are exact for linear solution problems and may induce various schemes by choosing different values for the parameter.In particular,when this parameter is taken to be 0,the nonlinear terms in our scheme degenerate to a harmonic average combination of the discrete linear fluxes,which has often been used in a variety of maximum-principle-preserving finite volume schemes.Thus our method of determining the coefficients of the nonlinear terms is more general.In addition,we prove the convergence of the proposed schemes by using a compactness technique.Numerical results demonstrate that our schemes can preserve the conservation property,satisfy the discrete maximum principle,possess a second-order accuracy,be exact for linear solution problems,and be available for anisotropic problems on distorted meshes.
基金the Natural Science Foundation of Shanghai(No.21ZR1465800)the Science Challenge Project(No.TZ2018001)+2 种基金the Interdisciplinary Project in Ocean Research of Tongji University,the Aeronautical Science Foundation of China(No.2020001053002)the National Key R&D Program of China(No.2020YFA0713603)the Fundamental Research Funds for the Central Universities.
文摘Stochastic temperature distribution should be carefully inspected in the thermal-failure design of heterogeneous solids with unexpected random energy excitations.Stochastic multiscale modeling for these problems involve multiscale and highdimensional uncertain thermal parameters,which remains limitation of prohibitive computation.In this paper,we propose a multi-modes based constrained energy minimization generalized multiscale finite element method(MCEM-GMsFEM),which can transform the original stochastic multiscale model into a series of recursive multiscale models sharing the same deterministic material parameters by multiscale analysis.Thus,MCEM-GMsFEM reveals an inherent low-dimensional representation in random space,and is designed to effectively reduce the complexity of repeated computation of discretized multiscale systems.In addition,the convergence analysis is established,and the optimal error estimates are derived.Finally,several typical random fluctuations on multiscale thermal conductivity are considered to validate the theoretical results in the numerical examples.The numerical results indicate that the multi-modes multiscale approach is a robust integrated method with the excellent performance.
