The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive comp...The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.展开更多
This paper is concerned with the Navier-Stokes/Allen-Cahn system,which is used to model the dynamics of immiscible two-phase flows.We consider a 1D free boundary problem and assume that the viscosity coefficient depen...This paper is concerned with the Navier-Stokes/Allen-Cahn system,which is used to model the dynamics of immiscible two-phase flows.We consider a 1D free boundary problem and assume that the viscosity coefficient depends on the density in the form ofη(ρ)=ρ^(α).The existence of unique global H^(2m)-solutions(m∈N)to the free boundary problem is proven for when 0<α<1/4.Furthermore,we obtain the global C^(∞)-solutions if the initial data is smooth.展开更多
The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundar...The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundary value problem of rectangular plates is proposed. The key concept behind this method is to transform the nonlinear or non-homogeneous part on the boundary into a lateral force within the governing function by the Dirac operator, which linearizes and homogenizes the original boundary, allowing one to employ the modal superposition method for obtaining solutions to reconstructive governing equations. Once projected into the modal space, the harmonic balance method(HBM) is utilized to solve coupled ordinary differential equations(ODEs)of truncated systems with nonlinearity. To validate the convergence and accuracy of the proposed Dirac method, the results of typical examples, involving nonlinearly restricted boundaries, moment excitation, and displacement excitation, are compared with those of the differential quadrature element method(DQEM). The results demonstrate that when dealing with nonlinear boundaries, the Dirac method exhibits more excellent accuracy and convergence compared with the DQEM. However, when facing displacement excitation, there exist some discrepancies between the proposed approach and simulations;nevertheless, the proposed method still accurately predicts resonant frequencies while being uniquely capable of handling nonuniform displacement excitations. Overall, this methodology offers a convenient way for addressing nonlinear and non-homogenous plate boundaries.展开更多
In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be r...In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.展开更多
The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundar...The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.展开更多
This paper is devoted to the study of the shape of the free boundary for a threedimensional axisymmetric incompressible impinging jet.To be more precise,we will show that the free boundary is convex to the fluid,provi...This paper is devoted to the study of the shape of the free boundary for a threedimensional axisymmetric incompressible impinging jet.To be more precise,we will show that the free boundary is convex to the fluid,provided that the uneven ground is concave to the fluid.展开更多
To investigate the potential of utilizing visible spectral imaging for controlling the plasma boundary shape during stable operation of plasma in future tokamak, a D_α band symmetric visible light diagnostic system w...To investigate the potential of utilizing visible spectral imaging for controlling the plasma boundary shape during stable operation of plasma in future tokamak, a D_α band symmetric visible light diagnostic system was designed and implemented on the Experimental Advanced Superconducting Tokamak(EAST). This system leverages two symmetric optics for joint plasma imaging. The optical system exhibits a spatial resolution less than 2 mm at the poloidal cross-section, distortion within the field of view below 10%, and relative illumination of 91%.The high-quality images obtained enable clear observation of both the plasma boundary position and the characteristics of components within the vacuum vessel. Following system calibration and coordinate transformation, the image coordinate boundary features are mapped to the tokamak coordinate system. Utilizing this system, the plasma boundary was reconstructed, and the resulting representation showed alignment with the EFIT(Equilibrium Fitting) results. This underscores the system's superior performance in boundary reconstruction applications and provides a diagnostic foundation for boundary shape control based on visible spectral imaging.展开更多
As a typical nonlinear wave,forward-leaning waves can be frequently encountered in the near-shore areas,which can impact coastal sediment transport significantly.Hence,it is of significance to describe the characteris...As a typical nonlinear wave,forward-leaning waves can be frequently encountered in the near-shore areas,which can impact coastal sediment transport significantly.Hence,it is of significance to describe the characteristics of the boundary layer beneath forward-leaning waves accurately,especially for the turbulent boundary layer.In this work,the linearized turbulent boundary layer model with a linear turbulent viscosity coefficient is applied,and the novel expression of the near-bed orbital velocity that has been worked out by the authors for forward-leaning waves of arbitrary forward-leaning degrees is further used to specify the free stream boundary condition of the bottom boundary layer.Then,a variable transformation is found so as to make the equation of the turbulent boundary layer model be solved analytically through a modified Bessel function.Consequently,an explicit analytical solution of the turbulent boundary layer beneath forward-leaning waves is derived by means of variable separation and variable transformation.The analytical solutions of the velocity profile and bottom shear stress of the turbulent boundary layer beneath forward-leaning waves are verified by comparing the present analytical results with typical experimental data available in the previous literature.展开更多
The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic...The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.展开更多
The boundary element method(BEM)is a popular method for solving acoustic wave propagation problems,especially those in exterior domains,owing to its ease in handling radiation conditions at infinity.However,BEM models...The boundary element method(BEM)is a popular method for solving acoustic wave propagation problems,especially those in exterior domains,owing to its ease in handling radiation conditions at infinity.However,BEM models must meet the requirement of 6–10 elements per wavelength,using the conventional constant,linear,or quadratic elements.Therefore,a large storage size of memory and long solution time are often needed in solving higher-frequency problems.In this work,we propose two new types of enriched elements based on conventional constant boundary elements to improve the computational efficiency of the 2D acoustic BEM.The first one uses a plane wave expansion,which can be used to model scattering problems.The second one uses a special plane wave expansion,which can be used tomodel radiation problems.Five examples are investigated to showthe advantages of the enriched elements.Compared with the conventional constant elements,the new enriched elements can deliver results with the same accuracy and in less computational time.This improvement in the computational efficiency is more evident at higher frequencies(with the nondimensional wave numbers exceeding 100).The paper concludes with the potential of our proposed enriched elements and plans for their further improvement.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.12072217 and 42077254)the Natural Science Foundation of Hunan Province,China(Grant No.2022JJ30567).
文摘The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.
基金supported by the Key Project of the NSFC(12131010)the NSFC(11771155,12271032)+1 种基金the NSF of Guangdong Province(2021A1515010249,2021A1515010303)supported by the NSFC(11971179,12371205)。
文摘This paper is concerned with the Navier-Stokes/Allen-Cahn system,which is used to model the dynamics of immiscible two-phase flows.We consider a 1D free boundary problem and assume that the viscosity coefficient depends on the density in the form ofη(ρ)=ρ^(α).The existence of unique global H^(2m)-solutions(m∈N)to the free boundary problem is proven for when 0<α<1/4.Furthermore,we obtain the global C^(∞)-solutions if the initial data is smooth.
基金Project supported by the National Natural Science Foundation of China (No. 12002195)the National Science Fund for Distinguished Young Scholars (No. 12025204)the Program of Shanghai Municipal Education Commission (No. 2019-01-07-00-09-E00018)。
文摘The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundary value problem of rectangular plates is proposed. The key concept behind this method is to transform the nonlinear or non-homogeneous part on the boundary into a lateral force within the governing function by the Dirac operator, which linearizes and homogenizes the original boundary, allowing one to employ the modal superposition method for obtaining solutions to reconstructive governing equations. Once projected into the modal space, the harmonic balance method(HBM) is utilized to solve coupled ordinary differential equations(ODEs)of truncated systems with nonlinearity. To validate the convergence and accuracy of the proposed Dirac method, the results of typical examples, involving nonlinearly restricted boundaries, moment excitation, and displacement excitation, are compared with those of the differential quadrature element method(DQEM). The results demonstrate that when dealing with nonlinear boundaries, the Dirac method exhibits more excellent accuracy and convergence compared with the DQEM. However, when facing displacement excitation, there exist some discrepancies between the proposed approach and simulations;nevertheless, the proposed method still accurately predicts resonant frequencies while being uniquely capable of handling nonuniform displacement excitations. Overall, this methodology offers a convenient way for addressing nonlinear and non-homogenous plate boundaries.
基金supported by the National Natural Science Foundation of China (No.12172154)the 111 Project (No.B14044)+1 种基金the Natural Science Foundation of Gansu Province (No.23JRRA1035)the Natural Science Foundation of Anhui University of Finance and Economics (No.ACKYC20043).
文摘In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.
基金The research work described herein was funded by the National Nature Science Foundation of China(Grant No.41877213).This financial support is gratefully acknowledged.
文摘The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.
基金supported in part by the National Natural Science Foundation of China(12101088)the Natural Science Foundation of Sichuan Province(2022NSFSC1858)。
文摘This paper is devoted to the study of the shape of the free boundary for a threedimensional axisymmetric incompressible impinging jet.To be more precise,we will show that the free boundary is convex to the fluid,provided that the uneven ground is concave to the fluid.
基金supported by the National MCF Energy R&D Program of China (Nos. 2018YFE0302103 and 2018YFE 0302100)National Natural Science Foundation of China (Nos. 12205195 and 11975277)。
文摘To investigate the potential of utilizing visible spectral imaging for controlling the plasma boundary shape during stable operation of plasma in future tokamak, a D_α band symmetric visible light diagnostic system was designed and implemented on the Experimental Advanced Superconducting Tokamak(EAST). This system leverages two symmetric optics for joint plasma imaging. The optical system exhibits a spatial resolution less than 2 mm at the poloidal cross-section, distortion within the field of view below 10%, and relative illumination of 91%.The high-quality images obtained enable clear observation of both the plasma boundary position and the characteristics of components within the vacuum vessel. Following system calibration and coordinate transformation, the image coordinate boundary features are mapped to the tokamak coordinate system. Utilizing this system, the plasma boundary was reconstructed, and the resulting representation showed alignment with the EFIT(Equilibrium Fitting) results. This underscores the system's superior performance in boundary reconstruction applications and provides a diagnostic foundation for boundary shape control based on visible spectral imaging.
基金Project supported by the National Key R&D Program of China(No.2022YFC3204303)the National Natural Science Foundation of China(Nos.12202503,12132018,and 52394254)。
文摘As a typical nonlinear wave,forward-leaning waves can be frequently encountered in the near-shore areas,which can impact coastal sediment transport significantly.Hence,it is of significance to describe the characteristics of the boundary layer beneath forward-leaning waves accurately,especially for the turbulent boundary layer.In this work,the linearized turbulent boundary layer model with a linear turbulent viscosity coefficient is applied,and the novel expression of the near-bed orbital velocity that has been worked out by the authors for forward-leaning waves of arbitrary forward-leaning degrees is further used to specify the free stream boundary condition of the bottom boundary layer.Then,a variable transformation is found so as to make the equation of the turbulent boundary layer model be solved analytically through a modified Bessel function.Consequently,an explicit analytical solution of the turbulent boundary layer beneath forward-leaning waves is derived by means of variable separation and variable transformation.The analytical solutions of the velocity profile and bottom shear stress of the turbulent boundary layer beneath forward-leaning waves are verified by comparing the present analytical results with typical experimental data available in the previous literature.
基金supported by the National Key R&D Program of China(2020YFA0710500).
文摘The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.
基金the National Natural Science Foundation of China(https://www.nsfc.gov.cn/,Project No.11972179)the Natural Science Foundation of Guangdong Province(http://gdstc.gd.gov.cn/,No.2020A1515010685)the Department of Education of Guangdong Province(http://edu.gd.gov.cn/,No.2020ZDZX2008).
文摘The boundary element method(BEM)is a popular method for solving acoustic wave propagation problems,especially those in exterior domains,owing to its ease in handling radiation conditions at infinity.However,BEM models must meet the requirement of 6–10 elements per wavelength,using the conventional constant,linear,or quadratic elements.Therefore,a large storage size of memory and long solution time are often needed in solving higher-frequency problems.In this work,we propose two new types of enriched elements based on conventional constant boundary elements to improve the computational efficiency of the 2D acoustic BEM.The first one uses a plane wave expansion,which can be used to model scattering problems.The second one uses a special plane wave expansion,which can be used tomodel radiation problems.Five examples are investigated to showthe advantages of the enriched elements.Compared with the conventional constant elements,the new enriched elements can deliver results with the same accuracy and in less computational time.This improvement in the computational efficiency is more evident at higher frequencies(with the nondimensional wave numbers exceeding 100).The paper concludes with the potential of our proposed enriched elements and plans for their further improvement.