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.展开更多
The aim of this study is to numerically investigate the impact of boundary slip on electroosmotic flow(EOF) in curved rectangular microchannels. Navier slip boundary conditions were employed at the curved microchannel...The aim of this study is to numerically investigate the impact of boundary slip on electroosmotic flow(EOF) in curved rectangular microchannels. Navier slip boundary conditions were employed at the curved microchannel walls. The electric potential distribution was governed by the Poisson–Boltzmann equation, whereas the velocity distribution was determined by the Navier–Stokes equation. The finite-difference method was employed to solve these two equations. The detailed discussion focuses on the impact of the curvature ratio, electrokinetic width, aspect ratio and slip length on the velocity. The results indicate that the present problem is strongly dependent on these parameters. The results demonstrate that by varying the dimensionless slip length from 0.001 to 0.01 while maintaining a curvature ratio of 0.5 there is a twofold increase in the maximum velocity. Moreover, this increase becomes more pronounced at higher curvature ratios. In addition, the velocity difference between the inner and outer radial regions increases with increasing slip length. Therefore, the incorporation of the slip boundary condition results in an augmented velocity and a more non-uniform velocity distribution. The findings presented here offer valuable insights into the design and optimization of EOF performance in curved hydrophobic microchannels featuring rectangular cross-sections.展开更多
Cross-flows around two,three and four circular cylinders in tandem,side-by-side,isosceles triangle and square arrangements are simulated using the incompressible lattice Boltzmann method with a second-order accurate c...Cross-flows around two,three and four circular cylinders in tandem,side-by-side,isosceles triangle and square arrangements are simulated using the incompressible lattice Boltzmann method with a second-order accurate curved boundary condition at Reynolds number 200 and the cylinder center-to-center transverse or/and longitudinal spacing 1.5D,where D is the identical circular cylinder diameter.The wake patterns,pressure and force distributions on the cylinders and mechanism of flow dynamics are investigated and compared among the four cases.The results also show that flows around the three or four cylinders significantly differ from those of the two cylinders in the tandem and side-by-side arrangements although there are some common features among the four cases due to their similarity of structures,which are interesting,complex and useful for practical applications.This study provides a useful database to validate the simplicity,accuracy and robustness of the Lattice Boltzmann method.展开更多
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.展开更多
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.展开更多
The impact of boundary layer suction on the aerodynamic performance of a high-turning compressor cascade was numerically simulated and discussed.The aerodynamic performance of a curved and a straight cascade with and ...The impact of boundary layer suction on the aerodynamic performance of a high-turning compressor cascade was numerically simulated and discussed.The aerodynamic performance of a curved and a straight cascade with and without boundary layer suction were comparatively studied at several suction flow rates.The results showed that boundary layer suction dramatically improved the flow behavior within the flow passage.Moreover,higher loading over the whole blade height,lower total pressure loss,and higher passage throughflow were achieved with a relatively small amount of boundary layer removal.The integration of curved blade and boundary layer suction contributed to better aerodynamic performance than the cascades with only curved blade or boundary layer suction used,and the more favorable effect resulted from the weakening of the three dimensional effects of the boundary layer close to the endwalls.展开更多
The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding freque...The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...展开更多
基金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 Natural Science Foundation of Inner Mongolia of China(Grant No.2021BS01008)the Program for Innovative Research Team in Universities of Inner Mongolia Autonomous Region(Grant No.NMGIRT2323)the Scientific Research Funding Project for introduced high level talents of IMNU(Grant No.2020YJRC014)。
文摘The aim of this study is to numerically investigate the impact of boundary slip on electroosmotic flow(EOF) in curved rectangular microchannels. Navier slip boundary conditions were employed at the curved microchannel walls. The electric potential distribution was governed by the Poisson–Boltzmann equation, whereas the velocity distribution was determined by the Navier–Stokes equation. The finite-difference method was employed to solve these two equations. The detailed discussion focuses on the impact of the curvature ratio, electrokinetic width, aspect ratio and slip length on the velocity. The results indicate that the present problem is strongly dependent on these parameters. The results demonstrate that by varying the dimensionless slip length from 0.001 to 0.01 while maintaining a curvature ratio of 0.5 there is a twofold increase in the maximum velocity. Moreover, this increase becomes more pronounced at higher curvature ratios. In addition, the velocity difference between the inner and outer radial regions increases with increasing slip length. Therefore, the incorporation of the slip boundary condition results in an augmented velocity and a more non-uniform velocity distribution. The findings presented here offer valuable insights into the design and optimization of EOF performance in curved hydrophobic microchannels featuring rectangular cross-sections.
基金supported by the Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education
文摘Cross-flows around two,three and four circular cylinders in tandem,side-by-side,isosceles triangle and square arrangements are simulated using the incompressible lattice Boltzmann method with a second-order accurate curved boundary condition at Reynolds number 200 and the cylinder center-to-center transverse or/and longitudinal spacing 1.5D,where D is the identical circular cylinder diameter.The wake patterns,pressure and force distributions on the cylinders and mechanism of flow dynamics are investigated and compared among the four cases.The results also show that flows around the three or four cylinders significantly differ from those of the two cylinders in the tandem and side-by-side arrangements although there are some common features among the four cases due to their similarity of structures,which are interesting,complex and useful for practical applications.This study provides a useful database to validate the simplicity,accuracy and robustness of the Lattice Boltzmann method.
基金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.
基金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.
基金Sponsored by the National Natural Science Foundation of China(Grant No.50646021).
文摘The impact of boundary layer suction on the aerodynamic performance of a high-turning compressor cascade was numerically simulated and discussed.The aerodynamic performance of a curved and a straight cascade with and without boundary layer suction were comparatively studied at several suction flow rates.The results showed that boundary layer suction dramatically improved the flow behavior within the flow passage.Moreover,higher loading over the whole blade height,lower total pressure loss,and higher passage throughflow were achieved with a relatively small amount of boundary layer removal.The integration of curved blade and boundary layer suction contributed to better aerodynamic performance than the cascades with only curved blade or boundary layer suction used,and the more favorable effect resulted from the weakening of the three dimensional effects of the boundary layer close to the endwalls.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20060056036)
文摘The lattice Boltzmann method (LBM) is employed to simulate the uniform flow past a circular cylinder. The performance of the two-dimensional LBM model on the prediction of force coefficients and vortex shedding frequency is investigated. The local grid refinement technique and second-order boundary condition for curved walls are applied in the calculations. It is found that the calculated vortex shedding frequency, drag coefficient and lift coefficient are consistent with experimental results at Reynolds nu...
