Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Latti...Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.展开更多
Recently,a Schwarz crystal structure with curved grain boundaries(GBs)constrained by twin-boundary(TB)networks was discovered in nanocrystalline Cu through experiments and atomistic simulations.Nanocrystalline Cu with...Recently,a Schwarz crystal structure with curved grain boundaries(GBs)constrained by twin-boundary(TB)networks was discovered in nanocrystalline Cu through experiments and atomistic simulations.Nanocrystalline Cu with nanosized Schwarz crystals exhibited high strength and excellent thermal stability.However,the grainsize effect and associated deformation mechanisms of Schwarz nanocrystals remain unknown.Here,we performed large-scale atomistic simulations to investigate the deformation behaviors and grain-size effect of nanocrystalline Cu with Schwarz crystals.Our simulations showed that similar to regular nanocrystals,Schwarz nanocrystals exhibit a strengthening-softening transition with decreasing grain size.The critical grain size in Schwarz nanocrystals is smaller than that in regular nanocrystals,leading to a maximum strength higher than that of regular nanocrystals.Our simulations revealed that the softening in Schwarz nanocrystals mainly originates from TB migration(or detwinning)and annihilation of GBs,rather than GB-mediated processes(including GB migration,sliding and diffusion)dominating the softening in regular nanocrystals.Quantitative analyses of simulation data further showed that compared with those in regular nanocrystals,the GB-mediated processes in Schwarz nanocrystals are suppressed,which is related to the low volume fraction of amorphous-like GBs and constraints of TB networks.The smaller critical grain size arises from the suppression of GB-mediated processes.展开更多
To deal with the staircase approximation problem in the standard finite-difference time-domain(FDTD) simulation,the two-dimensional boundary condition equations(BCE) method is proposed in this paper.In the BCE met...To deal with the staircase approximation problem in the standard finite-difference time-domain(FDTD) simulation,the two-dimensional boundary condition equations(BCE) method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.展开更多
In Li and Ren(Int.J.Numer.Methods Fluids 70:742–763,2012),a high-order k-exact WENO finite volume scheme based on secondary reconstructions was proposed to solve the two-dimensional time-dependent Euler equations in ...In Li and Ren(Int.J.Numer.Methods Fluids 70:742–763,2012),a high-order k-exact WENO finite volume scheme based on secondary reconstructions was proposed to solve the two-dimensional time-dependent Euler equations in a polygonal domain,in which the high-order numerical accuracy and the oscillations-free property can be achieved.In this paper,the method is extended to solve steady state problems imposed in a curved physical domain.The numerical framework consists of a Newton type finite volume method to linearize the nonlinear governing equations,and a geometrical multigrid method to solve the derived linear system.To achieve high-order non-oscillatory numerical solutions,the classical k-exact reconstruction with k=3 and the efficient secondary reconstructions are used to perform the WENO reconstruction for the conservative variables.The non-uniform rational B-splines(NURBS)curve is used to provide an exact or a high-order representation of the curved wall boundary.Furthermore,an enlarged reconstruction patch is constructed for every element of mesh to significantly improve the convergence to steady state.A variety of numerical examples are presented to show the effectiveness and robustness of the proposed method.展开更多
A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the...A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.展开更多
This paper proposes and applies a method to sort two-dimensional control points of triangular Bezier surfaces in a row vector. Using the property of bivariate Jacobi basis functions, it further presents two algorithms...This paper proposes and applies a method to sort two-dimensional control points of triangular Bezier surfaces in a row vector. Using the property of bivariate Jacobi basis functions, it further presents two algorithms for multi-degree reduction of triangular Bezier surfaces with constraints, providing explicit degree-reduced surfaces. The first algorithm can obtain the explicit representation of the optimal degree-reduced surfaces and the approximating error in both boundary curve constraints and corner constraints. But it has to solve the inversion of a matrix whose degree is related with the original surface. The second algorithm entails no matrix inversion to bring about computational instability, gives stable degree-reduced surfaces quickly, and presents the error bound. In the end, the paper proves the efficiency of the two algorithms through examples and error analysis.展开更多
For steady Euler equations in complex boundary domains,high-order shockcapturing schemes usually suffer not only from the difficulty of steady-state convergence but also from the problem of dealing with physical bound...For steady Euler equations in complex boundary domains,high-order shockcapturing schemes usually suffer not only from the difficulty of steady-state convergence but also from the problem of dealing with physical boundaries on Cartesian grids to achieve uniform high-order accuracy.In this paper,we utilize a fifth-order finite difference hybrid WENO scheme to simulate steady Euler equations,and the same fifth-order WENO extrapolation methods are developed to handle the curved boundary.The values of the ghost points outside the physical boundary can be obtained by applying WENO extrapolation near the boundary,involving normal derivatives acquired by the simplified inverse Lax-Wendroff procedure.Both equivalent expressions involving curvature and numerical differentiation are utilized to transform the tangential derivatives along the curved solid wall boundary.This hybrid WENO scheme is robust for steady-state convergence and maintains high-order accuracy in the smooth region even with the solid wall boundary condition.Besides,the essentially non-oscillation property is achieved.The numerical spectral analysis also shows that this hybrid WENO scheme has low dispersion and dissipation errors.Numerical examples are presented to validate the high-order accuracy and robust performance of the hybrid scheme for steady Euler equations in curved domains with Cartesian grids.展开更多
High-temperature mechanical properties of medium-entropy carbide ceramics have attracted significant attention.Tailoring the microstructure is an effective way to improve these high-temperature mechanical properties,w...High-temperature mechanical properties of medium-entropy carbide ceramics have attracted significant attention.Tailoring the microstructure is an effective way to improve these high-temperature mechanical properties,which can be affected by the evolution of the enthalpy and entropy,as well as by lattice distortion and sluggish diffusion.In this study,the effects of equiatomic Zr/(Ti,Nb)substitution(Zr content of 10-40 at%)on the microstructure and high-temperature strength of(Ti,Zr,Nb)C medium-entropy ceramics were investigated.The grain size of the(Ti,Zr,Nb)C medium-entropy ceramics was refined from 9.4±3.7 to 1.1±0.4μm with an increase in the Zr content from 10.0 to 33.3 at%.A further increase in the Zr content to 40 at%resulted in a slight increase in the grain size.At 1900℃,the(Ti,Zr,Nb)C medium-entropy ceramics with the Zr contents of 33.3 and 40 at%exhibited ultra-high flexural strengths of 875±43 and 843±71 MPa,respectively,which were higher than those of the transition metal carbides previously reported under similar conditions.Furthermore,relatively smooth grain boundaries,which were detected at a test temperature of 1000℃,transformed into curved and serrated boundaries as the temperature increased to 1900℃,which may be considered the primary reason for the improved high-temperature flexural strength.The associated mechanism was analyzed and discussed in detail.展开更多
基金WJD,JYZ,CLC,ZX,and ZGY were supported by the National Natural Science Foundation of China(Grant Number 51705143)the Education Department of Hunan Province(Grant Number 22B0464)the Postgraduate Scientific Research Innovation Project of Hunan Province(Grant Number QL20230249).
文摘Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.
基金the financial support from National Natural Science Foundation of China (Grants Nos.12325203,91963117,and 11921002)。
文摘Recently,a Schwarz crystal structure with curved grain boundaries(GBs)constrained by twin-boundary(TB)networks was discovered in nanocrystalline Cu through experiments and atomistic simulations.Nanocrystalline Cu with nanosized Schwarz crystals exhibited high strength and excellent thermal stability.However,the grainsize effect and associated deformation mechanisms of Schwarz nanocrystals remain unknown.Here,we performed large-scale atomistic simulations to investigate the deformation behaviors and grain-size effect of nanocrystalline Cu with Schwarz crystals.Our simulations showed that similar to regular nanocrystals,Schwarz nanocrystals exhibit a strengthening-softening transition with decreasing grain size.The critical grain size in Schwarz nanocrystals is smaller than that in regular nanocrystals,leading to a maximum strength higher than that of regular nanocrystals.Our simulations revealed that the softening in Schwarz nanocrystals mainly originates from TB migration(or detwinning)and annihilation of GBs,rather than GB-mediated processes(including GB migration,sliding and diffusion)dominating the softening in regular nanocrystals.Quantitative analyses of simulation data further showed that compared with those in regular nanocrystals,the GB-mediated processes in Schwarz nanocrystals are suppressed,which is related to the low volume fraction of amorphous-like GBs and constraints of TB networks.The smaller critical grain size arises from the suppression of GB-mediated processes.
基金Project supported by the National Natural Science Foundation of China(Grant No.51025622)
文摘To deal with the staircase approximation problem in the standard finite-difference time-domain(FDTD) simulation,the two-dimensional boundary condition equations(BCE) method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.
基金the Scientific Research Fund of Beijing Normal University(Grant No.28704-111032105)the Start-up Research Fund from BNU-HKBU United International College(Grant No.R72021112)+2 种基金The research of Guanghui Hu was partially supported by the FDCT of the Macao S.A.R.(0082/2020/A2)the National Natural Science Foundation of China(Grant Nos.11922120,11871489)the Multi-Year Research Grant(2019-00154-FST)of University of Macao,and a Grant from Department of Science and Technology of Guangdong Province(2020B1212030001).
文摘In Li and Ren(Int.J.Numer.Methods Fluids 70:742–763,2012),a high-order k-exact WENO finite volume scheme based on secondary reconstructions was proposed to solve the two-dimensional time-dependent Euler equations in a polygonal domain,in which the high-order numerical accuracy and the oscillations-free property can be achieved.In this paper,the method is extended to solve steady state problems imposed in a curved physical domain.The numerical framework consists of a Newton type finite volume method to linearize the nonlinear governing equations,and a geometrical multigrid method to solve the derived linear system.To achieve high-order non-oscillatory numerical solutions,the classical k-exact reconstruction with k=3 and the efficient secondary reconstructions are used to perform the WENO reconstruction for the conservative variables.The non-uniform rational B-splines(NURBS)curve is used to provide an exact or a high-order representation of the curved wall boundary.Furthermore,an enlarged reconstruction patch is constructed for every element of mesh to significantly improve the convergence to steady state.A variety of numerical examples are presented to show the effectiveness and robustness of the proposed method.
基金Supported by the National Natural Science Foundation of China (11071216 and 11101361)
文摘A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.
基金Supported by the National Natural Science Foundation of China (6087311160933007)
文摘This paper proposes and applies a method to sort two-dimensional control points of triangular Bezier surfaces in a row vector. Using the property of bivariate Jacobi basis functions, it further presents two algorithms for multi-degree reduction of triangular Bezier surfaces with constraints, providing explicit degree-reduced surfaces. The first algorithm can obtain the explicit representation of the optimal degree-reduced surfaces and the approximating error in both boundary curve constraints and corner constraints. But it has to solve the inversion of a matrix whose degree is related with the original surface. The second algorithm entails no matrix inversion to bring about computational instability, gives stable degree-reduced surfaces quickly, and presents the error bound. In the end, the paper proves the efficiency of the two algorithms through examples and error analysis.
基金Research supported by NSFC grant No.12271498National Key R&D Program of China No.2022YFA1005202/2022YFA1005200.
文摘For steady Euler equations in complex boundary domains,high-order shockcapturing schemes usually suffer not only from the difficulty of steady-state convergence but also from the problem of dealing with physical boundaries on Cartesian grids to achieve uniform high-order accuracy.In this paper,we utilize a fifth-order finite difference hybrid WENO scheme to simulate steady Euler equations,and the same fifth-order WENO extrapolation methods are developed to handle the curved boundary.The values of the ghost points outside the physical boundary can be obtained by applying WENO extrapolation near the boundary,involving normal derivatives acquired by the simplified inverse Lax-Wendroff procedure.Both equivalent expressions involving curvature and numerical differentiation are utilized to transform the tangential derivatives along the curved solid wall boundary.This hybrid WENO scheme is robust for steady-state convergence and maintains high-order accuracy in the smooth region even with the solid wall boundary condition.Besides,the essentially non-oscillation property is achieved.The numerical spectral analysis also shows that this hybrid WENO scheme has low dispersion and dissipation errors.Numerical examples are presented to validate the high-order accuracy and robust performance of the hybrid scheme for steady Euler equations in curved domains with Cartesian grids.
基金Financial supports from the National Natural Science Foundation of China(Nos.52172076,52032001,11575275,and 52102081)the State Key Laboratory of High Performance Ceramics and Superfine Microstructure are greatly appreciated.
文摘High-temperature mechanical properties of medium-entropy carbide ceramics have attracted significant attention.Tailoring the microstructure is an effective way to improve these high-temperature mechanical properties,which can be affected by the evolution of the enthalpy and entropy,as well as by lattice distortion and sluggish diffusion.In this study,the effects of equiatomic Zr/(Ti,Nb)substitution(Zr content of 10-40 at%)on the microstructure and high-temperature strength of(Ti,Zr,Nb)C medium-entropy ceramics were investigated.The grain size of the(Ti,Zr,Nb)C medium-entropy ceramics was refined from 9.4±3.7 to 1.1±0.4μm with an increase in the Zr content from 10.0 to 33.3 at%.A further increase in the Zr content to 40 at%resulted in a slight increase in the grain size.At 1900℃,the(Ti,Zr,Nb)C medium-entropy ceramics with the Zr contents of 33.3 and 40 at%exhibited ultra-high flexural strengths of 875±43 and 843±71 MPa,respectively,which were higher than those of the transition metal carbides previously reported under similar conditions.Furthermore,relatively smooth grain boundaries,which were detected at a test temperature of 1000℃,transformed into curved and serrated boundaries as the temperature increased to 1900℃,which may be considered the primary reason for the improved high-temperature flexural strength.The associated mechanism was analyzed and discussed in detail.