We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C^(1,α) estimates across th...We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C^(1,α) estimates across the discontinuity surfaces and provide an example to illustrate the issue regarding the regularity at the corners.展开更多
The cohesive zone model(CZM)has been used widely and successfully in fracture propagation,but some basic problems are still to be solved.In this paper,artificial compliance and discontinuous force in CZM are investiga...The cohesive zone model(CZM)has been used widely and successfully in fracture propagation,but some basic problems are still to be solved.In this paper,artificial compliance and discontinuous force in CZM are investigated.First,theories about the cohesive element(local coordinate system,stiffness matrix,and internal nodal force)are presented.The local coordinate system is defined to obtain local separation;the stiffness matrix for an eight-node cohesive element is derived from the calculation of strain energy;internal nodal force between the cohesive element and bulk element is obtained from the principle of virtual work.Second,the reason for artificial compliance is explained by the effective stiffnesses of zero-thickness and finite-thickness cohesive elements.Based on the effective stiffness,artificial compliance can be completely removed by adjusting the stiffness of the finite-thickness cohesive element.This conclusion is verified from 1D and 3D simulations.Third,three damage evolution methods(monotonically increasing effective separation,damage factor,and both effective separation and damage factor)are analyzed.Under constant unloading and reloading conditions,the monotonically increasing damage factor method without discontinuous force and healing effect is a better choice than the other two methods.The proposed improvements are coded in LS-DYNA user-defined material,and a drop weight tear test verifies the improvements.展开更多
Crack propagation in brittle material is not only crucial for structural safety evaluation,but also has a wideranging impact on material design,damage assessment,resource extraction,and scientific research.A thorough ...Crack propagation in brittle material is not only crucial for structural safety evaluation,but also has a wideranging impact on material design,damage assessment,resource extraction,and scientific research.A thorough investigation into the behavior of crack propagation contributes to a better understanding and control of the properties of brittle materials,thereby enhancing the reliability and safety of both materials and structures.As an implicit discrete elementmethod,the Discontinuous Deformation Analysis(DDA)has gained significant attention for its developments and applications in recent years.Among these developments,the particle DDA equipped with the bonded particle model is a powerful tool for predicting the whole process of material from continuity to failure.The primary objective of this research is to develop and utilize the particle DDAtomodel and understand the complex behavior of cracks in brittle materials under both static and dynamic loadings.The particle DDA is applied to several classical crack propagation problems,including the crack branching,compact tensile test,Kalthoff impact experiment,and tensile test of a rectangular plate with a hole.The evolutions of cracks under various stress or geometrical conditions are carefully investigated.The simulated results are compared with the experiments and other numerical results.It is found that the crack propagation patterns,including crack branching and the formation of secondary cracks,can be well reproduced.The results show that the particle DDA is a qualified method for crack propagation problems,providing valuable insights into the fracture mechanism of brittle materials.展开更多
We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separat...We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.展开更多
Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating land...Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating landslides.However,DDA fails to accurately capture the degradation in shear strength of rock joints commonly observed in high-speed landslides.In this study,DDA is modified by incorporating simplified joint shear strength degradation.Based on the modified DDA,the kinematics of the Baige landslide that occurred along the Jinsha River in China on 10 October 2018 are reproduced.The violent starting velocity of the landslide is considered explicitly.Three cases with different violent starting velocities are investigated to show their effect on the landslide movement process.Subsequently,the landslide movement process and the final accumulation characteristics are analyzed from multiple perspectives.The results show that the violent starting velocity affects the landslide motion characteristics,which is found to be about 4 m/s in the Baige landslide.The movement process of the Baige landslide involves four stages:initiation,high-speed sliding,impact-climbing,low-speed motion and accumulation.The accumulation states of sliding masses in different zones are different,which essentially corresponds to reality.The research results suggest that the modified DDA is applicable to similar high-level rock landslides.展开更多
This paper reviews the adaptive sparse grid discontinuous Galerkin(aSG-DG)method for computing high dimensional partial differential equations(PDEs)and its software implementation.The C++software package called AdaM-D...This paper reviews the adaptive sparse grid discontinuous Galerkin(aSG-DG)method for computing high dimensional partial differential equations(PDEs)and its software implementation.The C++software package called AdaM-DG,implementing the aSG-DG method,is available on GitHub at https://github.com/JuntaoHuang/adaptive-multiresolution-DG.The package is capable of treating a large class of high dimensional linear and nonlinear PDEs.We review the essential components of the algorithm and the functionality of the software,including the multiwavelets used,assembling of bilinear operators,fast matrix-vector product for data with hierarchical structures.We further demonstrate the performance of the package by reporting the numerical error and the CPU cost for several benchmark tests,including linear transport equations,wave equations,and Hamilton-Jacobi(HJ)equations.展开更多
In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the diver...In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.展开更多
In this paper,we propose a novel Local Macroscopic Conservative(LoMaC)low rank tensor method with discontinuous Galerkin(DG)discretization for the physical and phase spaces for simulating the Vlasov-Poisson(VP)system....In this paper,we propose a novel Local Macroscopic Conservative(LoMaC)low rank tensor method with discontinuous Galerkin(DG)discretization for the physical and phase spaces for simulating the Vlasov-Poisson(VP)system.The LoMaC property refers to the exact local conservation of macroscopic mass,momentum,and energy at the discrete level.The recently developed LoMaC low rank tensor algorithm(arXiv:2207.00518)simultaneously evolves the macroscopic conservation laws of mass,momentum,and energy using the kinetic flux vector splitting;then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables.This paper is a generalization of our previous work,but with DG discretization to take advantage of its compactness and flexibility in handling boundary conditions and its superior accuracy in the long term.The algorithm is developed in a similar fashion as that for a finite difference scheme,by observing that the DG method can be viewed equivalently in a nodal fashion.With the nodal DG method,assuming a tensorized computational grid,one will be able to(i)derive differentiation matrices for different nodal points based on a DG upwind discretization of transport terms,and(ii)define a weighted inner product space based on the nodal DG grid points.The algorithm can be extended to the high dimensional problems by hierarchical Tucker(HT)decomposition of solution tensors and a corresponding conservative projection algorithm.In a similar spirit,the algorithm can be extended to DG methods on nodal points of an unstructured mesh,or to other types of discretization,e.g.,the spectral method in velocity direction.Extensive numerical results are performed to showcase the efficacy of the method.展开更多
In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al...In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.展开更多
In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal e...In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.展开更多
This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis techniq...This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.展开更多
Due to the limited uplink capability in heterogeneousnetworks (HetNets), the decoupled uplinkand downlink access (DUDA) mode has recently beenproposed to improve the uplink performance. In thispaper, the random discon...Due to the limited uplink capability in heterogeneousnetworks (HetNets), the decoupled uplinkand downlink access (DUDA) mode has recently beenproposed to improve the uplink performance. In thispaper, the random discontinuous transmission (DTX)at user equipment (UE) is adopted to reduce the interferencecorrelation across different time slots. By utilizingstochastic geometry, we analytically derive themean local delay and energy efficiency (EE) of an uplinkHetNet with UE random DTX scheme under theDUDA mode. These expressions are further approximatedas closed forms under reasonable assumptions.Our results reveal that under the DUDA mode, there isan optimal EE with respect to mute probability underthe finite local delay constraint. In addition, with thesame finite mean local delay as under the coupled uplinkand downlink access (CUDA) mode, the HetNetsunder the DUDA mode can achieve a higher EE witha lower mute probability.展开更多
In this paper,we develop an entropy-conservative discontinuous Galerkin(DG)method for the shallow water(SW)equation with random inputs.One of the most popular methods for uncertainty quantifcation is the generalized P...In this paper,we develop an entropy-conservative discontinuous Galerkin(DG)method for the shallow water(SW)equation with random inputs.One of the most popular methods for uncertainty quantifcation is the generalized Polynomial Chaos(gPC)approach which we consider in the following manuscript.We apply the stochastic Galerkin(SG)method to the stochastic SW equations.Using the SG approach in the stochastic hyperbolic SW system yields a purely deterministic system that is not necessarily hyperbolic anymore.The lack of the hyperbolicity leads to ill-posedness and stability issues in numerical simulations.By transforming the system using Roe variables,the hyperbolicity can be ensured and an entropy-entropy fux pair is known from a recent investigation by Gerster and Herty(Commun.Comput.Phys.27(3):639–671,2020).We use this pair and determine a corresponding entropy fux potential.Then,we construct entropy conservative numerical twopoint fuxes for this augmented system.By applying these new numerical fuxes in a nodal DG spectral element method(DGSEM)with fux diferencing ansatz,we obtain a provable entropy conservative(dissipative)scheme.In numerical experiments,we validate our theoretical fndings.展开更多
We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations ...We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.展开更多
In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-depe...In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.展开更多
Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to r...Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.展开更多
Magnesium alloys are lightweight materials with great potential,and plasma electrolytic oxidation(PEO)is effective surface treatment for necessary improvement of corrosion resistance of magnesium alloys.However,the∼1...Magnesium alloys are lightweight materials with great potential,and plasma electrolytic oxidation(PEO)is effective surface treatment for necessary improvement of corrosion resistance of magnesium alloys.However,the∼14µm thick and rough PEO protection layer has inferior wear resistance,which limits magnesium alloys as sliding or reciprocating parts,where magnesium alloys have special advantages by their inherent damping and denoising properties and attractive light-weighting.Here a novel super wear-resistant coating for magnesium alloys was achieved,via the discontinuous sealing(DCS)of a 1.3µm thick polytetrafluoroethylene(PTFE)polymer layer with an initial area fraction(A_(f))of 70%on the necessary PEO protection layer by selective spraying,and the wear resistance was exceptionally enhanced by∼5500 times in comparison with the base PEO coating.The initial surface roughness(Sa)under PEO+DCS(1.54µm)was imperfectly 59%higher than that under PEO and conventional continuous sealing(CS).Interestingly,DCS was surprisingly 20 times superior for enhancing wear resistance in contrast to CS.DCS induced nano-cracks that splitted DCS layer into multilayer nano-blocks,and DCS also provided extra space for the movement of nano-blocks,which resulted in rolling friction and nano lubrication.Further,DCS promoted mixed wear of the PTFE polymer layer and the PEO coating,and the PTFE layer(HV:6 Kg·mm^(−2),A_(f):92.2%)and the PEO coating(HV:310 Kg·mm^(−2),A_(f):7.8%)served as the soft matrix and the hard point,respectively.Moreover,the dynamic decrease of Sa by 29%during wear also contributed to the super wear resistance.The strategy of depositing a low-frictional discontinuous layer on a rough and hard layer or matrix also opens a window for achieving super wear-resistant coatings in other materials.展开更多
In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy o...In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.展开更多
Accurate wave propagation simulation in anisotropic media is important for forward modeling, migration and inversion. In this study, the weighted Runge-Kutta discontinuous Galerkin (RKDG) method is extended to solve t...Accurate wave propagation simulation in anisotropic media is important for forward modeling, migration and inversion. In this study, the weighted Runge-Kutta discontinuous Galerkin (RKDG) method is extended to solve the elastic wave equations in 2D transversely isotropic media. The spatial discretization is based on the numerical flux discontinuous Galerkin scheme. An explicit weighted two-step iterative Runge-Kutta method is used as time-stepping algorithm. The weighted RKDG method has good flexibility and applicability of dealing with undulating geometries and boundary conditions. To verify the correctness and effectiveness of this method, several numerical examples are presented for elastic wave propagations in vertical transversely isotropic and tilted transversely isotropic media. The results show that the weighted RKDG method is promising for solving wave propagation problems in complex anisotropic medium.展开更多
In automotive industries,panel acoustic contribution analysis(PACA)is used to investigate the contributions of the body panels to the acoustic pressure at a certain point of interest.Currently,PACA is implementedmostl...In automotive industries,panel acoustic contribution analysis(PACA)is used to investigate the contributions of the body panels to the acoustic pressure at a certain point of interest.Currently,PACA is implementedmostly by either experiment-based methods or traditional numerical methods.However,these schemes are effort-consuming and inefficient in solving engineering problems,thereby restraining the further development of PACA in automotive acoustics.In this work,we propose a PACA scheme using discontinuous isogeometric boundary element method(IGABEM)to build an easily implementable and efficient method to identify the relative acoustic contributions of each automotive body panel.Discontinuous IGABEMis more accurate and converges faster than continuous BEM and IGABEM in the interior sound pressure evaluation of automotive compartments.In this work,a contribution ratio is defined to estimate the relative acoustic contribution of the structure panels;it can be calculated by reusing the coefficient matrix that has already been generated in the sound pressure evaluation process.The utilization of the parallel technique enables the proposed method to be more efficient than conventional methods;it is validated in two numerical examples,including a car passenger compartment subjected to realistic boundary conditions.A sound pressure response experiment based on a steel box is conducted to verify the accuracy of the interior sound pressure calculation using discontinuous IGABEM.This work is expected to promote the practical process of IGABEM for application in automotive acoustic problems.展开更多
基金supported by National Natural Science Foundation of China(12061080,12161087 and 12261093)the Science and Technology Project of the Education Department of Jiangxi Province(GJJ211601)supported by National Natural Science Foundation of China(11871305).
文摘We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C^(1,α) estimates across the discontinuity surfaces and provide an example to illustrate the issue regarding the regularity at the corners.
文摘The cohesive zone model(CZM)has been used widely and successfully in fracture propagation,but some basic problems are still to be solved.In this paper,artificial compliance and discontinuous force in CZM are investigated.First,theories about the cohesive element(local coordinate system,stiffness matrix,and internal nodal force)are presented.The local coordinate system is defined to obtain local separation;the stiffness matrix for an eight-node cohesive element is derived from the calculation of strain energy;internal nodal force between the cohesive element and bulk element is obtained from the principle of virtual work.Second,the reason for artificial compliance is explained by the effective stiffnesses of zero-thickness and finite-thickness cohesive elements.Based on the effective stiffness,artificial compliance can be completely removed by adjusting the stiffness of the finite-thickness cohesive element.This conclusion is verified from 1D and 3D simulations.Third,three damage evolution methods(monotonically increasing effective separation,damage factor,and both effective separation and damage factor)are analyzed.Under constant unloading and reloading conditions,the monotonically increasing damage factor method without discontinuous force and healing effect is a better choice than the other two methods.The proposed improvements are coded in LS-DYNA user-defined material,and a drop weight tear test verifies the improvements.
基金supported by the National Natural Science Foundation of China(Grant No.42372310).
文摘Crack propagation in brittle material is not only crucial for structural safety evaluation,but also has a wideranging impact on material design,damage assessment,resource extraction,and scientific research.A thorough investigation into the behavior of crack propagation contributes to a better understanding and control of the properties of brittle materials,thereby enhancing the reliability and safety of both materials and structures.As an implicit discrete elementmethod,the Discontinuous Deformation Analysis(DDA)has gained significant attention for its developments and applications in recent years.Among these developments,the particle DDA equipped with the bonded particle model is a powerful tool for predicting the whole process of material from continuity to failure.The primary objective of this research is to develop and utilize the particle DDAtomodel and understand the complex behavior of cracks in brittle materials under both static and dynamic loadings.The particle DDA is applied to several classical crack propagation problems,including the crack branching,compact tensile test,Kalthoff impact experiment,and tensile test of a rectangular plate with a hole.The evolutions of cracks under various stress or geometrical conditions are carefully investigated.The simulated results are compared with the experiments and other numerical results.It is found that the crack propagation patterns,including crack branching and the formation of secondary cracks,can be well reproduced.The results show that the particle DDA is a qualified method for crack propagation problems,providing valuable insights into the fracture mechanism of brittle materials.
基金supported by the National Natural Science Foundation of China(11871218,12071298)in part by the Science and Technology Commission of Shanghai Municipality(21JC1402500,22DZ2229014)。
文摘We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.
基金supported by the National Natural Science Foundations of China(grant numbers U22A20601 and 52209142)the Opening fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(grant number SKLGP2022K018)+1 种基金the Science&Technology Department of Sichuan Province(grant number 2023NSFSC0284)the Science and Technology Major Project of Tibetan Autonomous Region of China(grant number XZ202201ZD0003G)。
文摘Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating landslides.However,DDA fails to accurately capture the degradation in shear strength of rock joints commonly observed in high-speed landslides.In this study,DDA is modified by incorporating simplified joint shear strength degradation.Based on the modified DDA,the kinematics of the Baige landslide that occurred along the Jinsha River in China on 10 October 2018 are reproduced.The violent starting velocity of the landslide is considered explicitly.Three cases with different violent starting velocities are investigated to show their effect on the landslide movement process.Subsequently,the landslide movement process and the final accumulation characteristics are analyzed from multiple perspectives.The results show that the violent starting velocity affects the landslide motion characteristics,which is found to be about 4 m/s in the Baige landslide.The movement process of the Baige landslide involves four stages:initiation,high-speed sliding,impact-climbing,low-speed motion and accumulation.The accumulation states of sliding masses in different zones are different,which essentially corresponds to reality.The research results suggest that the modified DDA is applicable to similar high-level rock landslides.
基金supported by the NSF grant DMS-2111383Air Force Office of Scientific Research FA9550-18-1-0257the NSF grant DMS-2011838.
文摘This paper reviews the adaptive sparse grid discontinuous Galerkin(aSG-DG)method for computing high dimensional partial differential equations(PDEs)and its software implementation.The C++software package called AdaM-DG,implementing the aSG-DG method,is available on GitHub at https://github.com/JuntaoHuang/adaptive-multiresolution-DG.The package is capable of treating a large class of high dimensional linear and nonlinear PDEs.We review the essential components of the algorithm and the functionality of the software,including the multiwavelets used,assembling of bilinear operators,fast matrix-vector product for data with hierarchical structures.We further demonstrate the performance of the package by reporting the numerical error and the CPU cost for several benchmark tests,including linear transport equations,wave equations,and Hamilton-Jacobi(HJ)equations.
基金supported by the NSFC Grant 11901555,12271499the Cyrus Tang Foundationsupported by the NSFC Grant 11871448 and 12126604.
文摘In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.
基金supported by the NSF(Grant Nos.the NSF-DMS-1818924 and 2111253)the Air Force Office of Scientific Research FA9550-22-1-0390 and Department of Energy DE-SC0023164+1 种基金supported by the NSF(Grant Nos.NSF-DMS-1830838 and NSF-DMS-2111383)the Air Force Office of Scientific Research FA9550-22-1-0390.
文摘In this paper,we propose a novel Local Macroscopic Conservative(LoMaC)low rank tensor method with discontinuous Galerkin(DG)discretization for the physical and phase spaces for simulating the Vlasov-Poisson(VP)system.The LoMaC property refers to the exact local conservation of macroscopic mass,momentum,and energy at the discrete level.The recently developed LoMaC low rank tensor algorithm(arXiv:2207.00518)simultaneously evolves the macroscopic conservation laws of mass,momentum,and energy using the kinetic flux vector splitting;then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables.This paper is a generalization of our previous work,but with DG discretization to take advantage of its compactness and flexibility in handling boundary conditions and its superior accuracy in the long term.The algorithm is developed in a similar fashion as that for a finite difference scheme,by observing that the DG method can be viewed equivalently in a nodal fashion.With the nodal DG method,assuming a tensorized computational grid,one will be able to(i)derive differentiation matrices for different nodal points based on a DG upwind discretization of transport terms,and(ii)define a weighted inner product space based on the nodal DG grid points.The algorithm can be extended to the high dimensional problems by hierarchical Tucker(HT)decomposition of solution tensors and a corresponding conservative projection algorithm.In a similar spirit,the algorithm can be extended to DG methods on nodal points of an unstructured mesh,or to other types of discretization,e.g.,the spectral method in velocity direction.Extensive numerical results are performed to showcase the efficacy of the method.
基金supported by the National Science Foundation grant DMS-1818998.
文摘In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.
基金supported by the NSF under Grant DMS-1818467Simons Foundation under Grant 961585.
文摘In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.
基金supported by the National Natural Science Foundation of China(Grant Nos.11871428 and 12071214)the Natural Science Foundation for Colleges and Universities of Jiangsu Province of China(Grant No.20KJB110011)+1 种基金supported by the National Science Foundation(Grant No.DMS-1620335)and the Simons Foundation(Grant No.637716)supported by the National Natural Science Foundation of China(Grant Nos.11871428 and 12272347).
文摘This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.
基金supported in part by the National Key R&D Program of China under Grant 2021YFB 2900304the Shenzhen Science and Technology Program under Grants KQTD20190929172545139 and ZDSYS20210623091808025.
文摘Due to the limited uplink capability in heterogeneousnetworks (HetNets), the decoupled uplinkand downlink access (DUDA) mode has recently beenproposed to improve the uplink performance. In thispaper, the random discontinuous transmission (DTX)at user equipment (UE) is adopted to reduce the interferencecorrelation across different time slots. By utilizingstochastic geometry, we analytically derive themean local delay and energy efficiency (EE) of an uplinkHetNet with UE random DTX scheme under theDUDA mode. These expressions are further approximatedas closed forms under reasonable assumptions.Our results reveal that under the DUDA mode, there isan optimal EE with respect to mute probability underthe finite local delay constraint. In addition, with thesame finite mean local delay as under the coupled uplinkand downlink access (CUDA) mode, the HetNetsunder the DUDA mode can achieve a higher EE witha lower mute probability.
文摘In this paper,we develop an entropy-conservative discontinuous Galerkin(DG)method for the shallow water(SW)equation with random inputs.One of the most popular methods for uncertainty quantifcation is the generalized Polynomial Chaos(gPC)approach which we consider in the following manuscript.We apply the stochastic Galerkin(SG)method to the stochastic SW equations.Using the SG approach in the stochastic hyperbolic SW system yields a purely deterministic system that is not necessarily hyperbolic anymore.The lack of the hyperbolicity leads to ill-posedness and stability issues in numerical simulations.By transforming the system using Roe variables,the hyperbolicity can be ensured and an entropy-entropy fux pair is known from a recent investigation by Gerster and Herty(Commun.Comput.Phys.27(3):639–671,2020).We use this pair and determine a corresponding entropy fux potential.Then,we construct entropy conservative numerical twopoint fuxes for this augmented system.By applying these new numerical fuxes in a nodal DG spectral element method(DGSEM)with fux diferencing ansatz,we obtain a provable entropy conservative(dissipative)scheme.In numerical experiments,we validate our theoretical fndings.
文摘We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.
文摘In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.
基金The authors would like to thank the National Natural Science Foundation of China(Grant Nos.51879184 and 52079091)for funding this work.
文摘Since its introduction,discontinuous deformation analysis(DDA)has been widely used in different areas of rock mechanics.By dividing large blocks into subblocks and introducing artificial joints,DDA can be applied to rock fracture simulation.However,parameter calibration,a fundamental issue in discontinuum methods,has not received enough attention in DDA.In this study,the parameter calibration of DDA for intact rock is carefully studied.To this end,a subblock DDA with Voronoi tessellation is presented first.Then,a modified contact constitutive law is introduced,in which the tensile and shear meso-strengths are modified to be independent of the bond lengths.This improvement can prevent the unjustified preferential failure of short edges.A method for imposing confining pressure is also introduced.Thereafter,sensitivity analysis is performed to investigate the influence of the calculated parameters and meso-parameters on the mechanical properties of modeled rock.Based on the sensitivity analysis,a unified calibration procedure is suggested for both cases with and without confining pressure.Finally,the calibration procedure is applied to two examples,including a biaxial compression test.The results show that the proposed Voronoi-based DDA can simulate rock fracture with and without confining pressure very well after careful parameter calibration.
基金This work was financially supported by the Jiangsu Distinguished Professor Project,the Innovate UK(Project reference:10004694)the National Key R&D Program of China 2021YFB3401200.The Experimental Techniques Centre at Brunel University London and Nanjing University of Aeronautics and Astronautics are acknowledged.The authors also acknowledge the characterization facility at Shanghai Jiao Tong University,Central South University,University of Birmingham and University of Lille.
文摘Magnesium alloys are lightweight materials with great potential,and plasma electrolytic oxidation(PEO)is effective surface treatment for necessary improvement of corrosion resistance of magnesium alloys.However,the∼14µm thick and rough PEO protection layer has inferior wear resistance,which limits magnesium alloys as sliding or reciprocating parts,where magnesium alloys have special advantages by their inherent damping and denoising properties and attractive light-weighting.Here a novel super wear-resistant coating for magnesium alloys was achieved,via the discontinuous sealing(DCS)of a 1.3µm thick polytetrafluoroethylene(PTFE)polymer layer with an initial area fraction(A_(f))of 70%on the necessary PEO protection layer by selective spraying,and the wear resistance was exceptionally enhanced by∼5500 times in comparison with the base PEO coating.The initial surface roughness(Sa)under PEO+DCS(1.54µm)was imperfectly 59%higher than that under PEO and conventional continuous sealing(CS).Interestingly,DCS was surprisingly 20 times superior for enhancing wear resistance in contrast to CS.DCS induced nano-cracks that splitted DCS layer into multilayer nano-blocks,and DCS also provided extra space for the movement of nano-blocks,which resulted in rolling friction and nano lubrication.Further,DCS promoted mixed wear of the PTFE polymer layer and the PEO coating,and the PTFE layer(HV:6 Kg·mm^(−2),A_(f):92.2%)and the PEO coating(HV:310 Kg·mm^(−2),A_(f):7.8%)served as the soft matrix and the hard point,respectively.Moreover,the dynamic decrease of Sa by 29%during wear also contributed to the super wear resistance.The strategy of depositing a low-frictional discontinuous layer on a rough and hard layer or matrix also opens a window for achieving super wear-resistant coatings in other materials.
文摘In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.
基金supported by the National Natural Science Foundation of China(Grant Nos.41974114,41604105)the Fundamental Research Funds for the Central Universities(2020YQLX01)+1 种基金supported in part by the Project of Cultivation for Young Top-notch Talents of Beijing Municipal Institutions under Grant BPHR202203047in part by the Young Elite Scientists Sponsorship Program by BAST.
文摘Accurate wave propagation simulation in anisotropic media is important for forward modeling, migration and inversion. In this study, the weighted Runge-Kutta discontinuous Galerkin (RKDG) method is extended to solve the elastic wave equations in 2D transversely isotropic media. The spatial discretization is based on the numerical flux discontinuous Galerkin scheme. An explicit weighted two-step iterative Runge-Kutta method is used as time-stepping algorithm. The weighted RKDG method has good flexibility and applicability of dealing with undulating geometries and boundary conditions. To verify the correctness and effectiveness of this method, several numerical examples are presented for elastic wave propagations in vertical transversely isotropic and tilted transversely isotropic media. The results show that the weighted RKDG method is promising for solving wave propagation problems in complex anisotropic medium.
基金funded by the National Natural Science Foundation of China (Grant No.52175111)。
文摘In automotive industries,panel acoustic contribution analysis(PACA)is used to investigate the contributions of the body panels to the acoustic pressure at a certain point of interest.Currently,PACA is implementedmostly by either experiment-based methods or traditional numerical methods.However,these schemes are effort-consuming and inefficient in solving engineering problems,thereby restraining the further development of PACA in automotive acoustics.In this work,we propose a PACA scheme using discontinuous isogeometric boundary element method(IGABEM)to build an easily implementable and efficient method to identify the relative acoustic contributions of each automotive body panel.Discontinuous IGABEMis more accurate and converges faster than continuous BEM and IGABEM in the interior sound pressure evaluation of automotive compartments.In this work,a contribution ratio is defined to estimate the relative acoustic contribution of the structure panels;it can be calculated by reusing the coefficient matrix that has already been generated in the sound pressure evaluation process.The utilization of the parallel technique enables the proposed method to be more efficient than conventional methods;it is validated in two numerical examples,including a car passenger compartment subjected to realistic boundary conditions.A sound pressure response experiment based on a steel box is conducted to verify the accuracy of the interior sound pressure calculation using discontinuous IGABEM.This work is expected to promote the practical process of IGABEM for application in automotive acoustic problems.