Since the West Pingdingshan Section in Chaohu was proposed as the candidate of the Global Stratotype Section and Point of the Induan-Olenekian boundary in 2003, the Lower Triassic of Chaohu has been extensively studie...Since the West Pingdingshan Section in Chaohu was proposed as the candidate of the Global Stratotype Section and Point of the Induan-Olenekian boundary in 2003, the Lower Triassic of Chaohu has been extensively studied. Based on the studies on the Lower Triassic of Chaohu, (1) a continuous conodont zonation is established, which has become an important reference for Lower Triassic stratigraphic correlation over the world; (2) the First Appearance Datum of conodont Neospathodus waageni was suggested and has been basically accepted as the primary marker to define the InduanOienekian boundary; (3) a characteristic Lower Triassic excursion of carbon isotopes was brought to light and has been proven to be not only an excellent index for the stratigraphic correlation but also a unique indication for the perturbation of ecological environments in the aftermath of the end-Permian mass extinction; (4) a magnetostratigraphic sequence is constituted with a certain biostratigraphic control in the low-latitude region and it presents an important correlation to the Boreal sequence; (5) a cyclostratigraphic study provides an alternative method to constrain the age of the chronostratigraphic units; and (6) a scheme of the Olenekian subdivision is recently suggested to define the boundary between the Smithian and Spathian Substages. In addition, Chaohu is also the type locality of the Chaohuan Stage, the upper stage of the Lower Triassic in the China Chronostratigraphic System.Thus, the Lower Triassic of Chaohu is not only a classic sequence in South China, but also a key reference sequence to the investigation of the corresponding stratigraphy and geological events over the world. The recent achievements are viewed here for an overall understanding of the sequence. Then the current situation of the Induan-Olenekian and Smithian-Spathian boundaries is discussed to provide a reference for later works.展开更多
The Lower Triassic in Chaohu (巢湖) area, Anhui (安徽) Province, China, is well developed and its sequence is typical in South China. After a brief introduction of the Induan-Olenekian boundary of Chaohu, this art...The Lower Triassic in Chaohu (巢湖) area, Anhui (安徽) Province, China, is well developed and its sequence is typical in South China. After a brief introduction of the Induan-Olenekian boundary of Chaohu, this article presents some new data on conodonts. More than ten times of conodont samplings and investigations have recovered thousands of conodont specimens, which are especially rich in the Induan-Olenekian boundary strata at the West Pingdingshan Section in Chaohu City, Anhui Province. The most distinctive forms are the conodonts of the Neospathodus dieneri group and N. waageni group. The first occurrence of N. waageni eowaageni, which is regarded as the indicator of the Induan-Olenekian boundary, is situated at 40.49 m above the base of Yinkeng (殷坑) Formation. Some key conodonts and seven new specimens are introduced.展开更多
This paper proposes a scheme for the definition of the Lower Triassic Induan Olenekian boundary (IOB) based on investigation of sections in Chaohu, Anhui Province, China as well as data accumulated from other studies...This paper proposes a scheme for the definition of the Lower Triassic Induan Olenekian boundary (IOB) based on investigation of sections in Chaohu, Anhui Province, China as well as data accumulated from other studies elsewhere. The conodont Neospathodus waageni is suggested as the index fossil of the boundary. According to the FAD of N. waageni , the IOB is at the base of bed 25 2 of the West Pingdingshan Section in Chaohu, 42.19 m above the Permian Triassic boundary, and it is slightly higher than the base of the Flemingites Euflemingites Ammonoid Zone at the section.展开更多
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.展开更多
Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the...Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the formation energy,excess volume,dislocations and melting behaviors of GBs in CuFe alloys.It is illustrated that Fe solute affects the structural stability of Cu GBs substantially,the formation energy of GBs is reduced,but the thickness and melting point of GBs are increased,that is,the structural stability of Cu GBs is significantly improved owing to the Fe solutes.A strong scaling law exists between the formation energy,excess volume,thickness and melting point of GBs.Therefore,Fe solid solute plays an important role in the characteristics of GBs in bi-crystal Cu.展开更多
The experimental results in previous studies have indicated that during the ductile fracture of pure metals,vacancies aggregate and form voids at grain boundaries.However,the physical mechanism underlying this phenome...The experimental results in previous studies have indicated that during the ductile fracture of pure metals,vacancies aggregate and form voids at grain boundaries.However,the physical mechanism underlying this phenomenon remains not fully understood.This study derives the equilibrium distribution of vacancies analytically by following thermodynamics and the micromechanics of crystal defects.This derivation suggests that vacancies cluster in regions under hydrostatic compression to minimize the elastic strain energy.Subsequently,a finite element model is developed for examining more general scenarios of interaction between vacancies and grain boundaries.This model is first verified and validated through comparison with some available analytical solutions,demonstrating consistency between finite element simulation results and analytical solutions within a specified numerical accuracy.A systematic numerical study is then conducted to investigate the mechanism that might govern the micromechanical interaction between grain boundaries and the profuse vacancies typically generated during plastic deformation.The simulation results indicate that the reduction in total elastic strain energy can indeed drive vacancies toward grain boundaries,potentially facilitating void nucleation in ductile fracture.展开更多
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.展开更多
In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Mill...In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.展开更多
Clifford analysis is an important branch of modern analysis;it has a very important theoretical significance and application value,and its conclusions can be applied to the Maxwell equation,Yang-Mill field theory,quan...Clifford analysis is an important branch of modern analysis;it has a very important theoretical significance and application value,and its conclusions can be applied to the Maxwell equation,Yang-Mill field theory,quantum mechanics and value problems.In this paper,we first give the definition of a quasi-Cauchy type integral in complex Clifford analysis,and get the Plemelj formula for it.Second,we discuss the H?lder continuity for the Cauchy-type integral operators with values in a complex Clifford algebra.Finally,we prove the existence of solutions for a class of linear boundary value problems and give the integral representation for the solution.展开更多
The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is nece...The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.展开更多
For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geomet...For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.展开更多
The helium bubbles induced by 14 MeV neutron irradiation can cause intergranular fractures in reduced activation ferritic martensitic steel,which is a candidate structural material for fusion reactors.In order to eluc...The helium bubbles induced by 14 MeV neutron irradiation can cause intergranular fractures in reduced activation ferritic martensitic steel,which is a candidate structural material for fusion reactors.In order to elucidate the susceptibility of different grain boundaries(GBs)to helium-induced embrittlement,the tensile fracture processes of 10 types of GBs with and without helium bubbles in body-centered cubic(bcc)iron at the relevant service temperature of 600 K were investigated via molecular dynamics methods.The results indicate that in the absence of helium bubbles,the GBs studied here can be classified into two distinct categories:brittle GBs and ductile GBs.The atomic scale analysis shows that the plastic deformation of ductile GB at high temperatures originates from complex plastic deformation mechanisms,including the Bain/Burgers path phase transition and deformation twinning,in which the Bain path phase transition is the most dominant plastic deformation mechanism.However,the presence of helium bubbles severely inhibits the plastic deformation channels of the GBs,resulting in a significant decrease in elongation at fractures.For bubble-decorated GBs,the ultimate tensile strength increases with the increase in the misorientation angle.Interestingly,the coherent twin boundary∑3{112}was found to maintain relatively high fracture strength and maximum failure strain under the influence of helium bubbles.展开更多
We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and...We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.展开更多
基金supported by the Ministry of Science and Technology project (no. 2006FY120300-11)Ministry of Education project (no. 200804910503)National Natural Science Foundation of China (nos. 40830212, 40921062, 40972003)
文摘Since the West Pingdingshan Section in Chaohu was proposed as the candidate of the Global Stratotype Section and Point of the Induan-Olenekian boundary in 2003, the Lower Triassic of Chaohu has been extensively studied. Based on the studies on the Lower Triassic of Chaohu, (1) a continuous conodont zonation is established, which has become an important reference for Lower Triassic stratigraphic correlation over the world; (2) the First Appearance Datum of conodont Neospathodus waageni was suggested and has been basically accepted as the primary marker to define the InduanOienekian boundary; (3) a characteristic Lower Triassic excursion of carbon isotopes was brought to light and has been proven to be not only an excellent index for the stratigraphic correlation but also a unique indication for the perturbation of ecological environments in the aftermath of the end-Permian mass extinction; (4) a magnetostratigraphic sequence is constituted with a certain biostratigraphic control in the low-latitude region and it presents an important correlation to the Boreal sequence; (5) a cyclostratigraphic study provides an alternative method to constrain the age of the chronostratigraphic units; and (6) a scheme of the Olenekian subdivision is recently suggested to define the boundary between the Smithian and Spathian Substages. In addition, Chaohu is also the type locality of the Chaohuan Stage, the upper stage of the Lower Triassic in the China Chronostratigraphic System.Thus, the Lower Triassic of Chaohu is not only a classic sequence in South China, but also a key reference sequence to the investigation of the corresponding stratigraphy and geological events over the world. The recent achievements are viewed here for an overall understanding of the sequence. Then the current situation of the Induan-Olenekian and Smithian-Spathian boundaries is discussed to provide a reference for later works.
基金the National Natural Science Foundation of China (Nos. 40621002, 90714010, 40574028)the State Key Laboratory of Geological Processes and Mineral Resources (GPMR2007).
文摘The Lower Triassic in Chaohu (巢湖) area, Anhui (安徽) Province, China, is well developed and its sequence is typical in South China. After a brief introduction of the Induan-Olenekian boundary of Chaohu, this article presents some new data on conodonts. More than ten times of conodont samplings and investigations have recovered thousands of conodont specimens, which are especially rich in the Induan-Olenekian boundary strata at the West Pingdingshan Section in Chaohu City, Anhui Province. The most distinctive forms are the conodonts of the Neospathodus dieneri group and N. waageni group. The first occurrence of N. waageni eowaageni, which is regarded as the indicator of the Induan-Olenekian boundary, is situated at 40.49 m above the base of Yinkeng (殷坑) Formation. Some key conodonts and seven new specimens are introduced.
基金ThisresearchissupportedbytheNationalNaturalScienceFoundationofChina (No .40 0 72 0 11) theMinistryofScienceandTechnology (No .2 0 0 1DEA2 0 0 2 0 )andtheChinese"973Program" (No .G2 0 0 0 0 7770 5 ) .
文摘This paper proposes a scheme for the definition of the Lower Triassic Induan Olenekian boundary (IOB) based on investigation of sections in Chaohu, Anhui Province, China as well as data accumulated from other studies elsewhere. The conodont Neospathodus waageni is suggested as the index fossil of the boundary. According to the FAD of N. waageni , the IOB is at the base of bed 25 2 of the West Pingdingshan Section in Chaohu, 42.19 m above the Permian Triassic boundary, and it is slightly higher than the base of the Flemingites Euflemingites Ammonoid Zone at the section.
基金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.
基金supported by National Key Research and Development Program of China(No.2021YFB3400800)National Natural Science Foundation of China(Grant No.52271136,51901177)Natural Science Foundation of Shaanxi Province(No.2021JC-06,2019TD-020).
文摘Grain boundaries(GBs)play a crucial role on the structural stability and mechanical properties of Cu and its alloys.In this work,molecular dynamics(MD)simulations are employed to study the effects of Fe solutes on the formation energy,excess volume,dislocations and melting behaviors of GBs in CuFe alloys.It is illustrated that Fe solute affects the structural stability of Cu GBs substantially,the formation energy of GBs is reduced,but the thickness and melting point of GBs are increased,that is,the structural stability of Cu GBs is significantly improved owing to the Fe solutes.A strong scaling law exists between the formation energy,excess volume,thickness and melting point of GBs.Therefore,Fe solid solute plays an important role in the characteristics of GBs in bi-crystal Cu.
基金supported by the National Key Research and Development Program of China under Grant No.2023YFB3712401the National Natural Science Foundation of China under Grant Nos.12102254 and 12327802.
文摘The experimental results in previous studies have indicated that during the ductile fracture of pure metals,vacancies aggregate and form voids at grain boundaries.However,the physical mechanism underlying this phenomenon remains not fully understood.This study derives the equilibrium distribution of vacancies analytically by following thermodynamics and the micromechanics of crystal defects.This derivation suggests that vacancies cluster in regions under hydrostatic compression to minimize the elastic strain energy.Subsequently,a finite element model is developed for examining more general scenarios of interaction between vacancies and grain boundaries.This model is first verified and validated through comparison with some available analytical solutions,demonstrating consistency between finite element simulation results and analytical solutions within a specified numerical accuracy.A systematic numerical study is then conducted to investigate the mechanism that might govern the micromechanical interaction between grain boundaries and the profuse vacancies typically generated during plastic deformation.The simulation results indicate that the reduction in total elastic strain energy can indeed drive vacancies toward grain boundaries,potentially facilitating void nucleation in ductile fracture.
基金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.
基金sponsored by the Graduate Student Research and Innovation Fund of Xinyang Normal University under No.2024KYJJ012.
文摘In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.
基金supported by the NSF of Hebei Province(A2022208007)the NSF of China(11571089,11871191)the NSF of Henan Province(222300420397)。
文摘Clifford analysis is an important branch of modern analysis;it has a very important theoretical significance and application value,and its conclusions can be applied to the Maxwell equation,Yang-Mill field theory,quantum mechanics and value problems.In this paper,we first give the definition of a quasi-Cauchy type integral in complex Clifford analysis,and get the Plemelj formula for it.Second,we discuss the H?lder continuity for the Cauchy-type integral operators with values in a complex Clifford algebra.Finally,we prove the existence of solutions for a class of linear boundary value problems and give the integral representation for the solution.
基金Project supported by the National Natural Science Foundation of China (Nos. 52075070 and12302254)the Dalian City Supports Innovation and Entrepreneurship Projects for High-Level Talents (No. 2021RD16)the Liaoning Revitalization Talents Program (No. XLYC2002108)。
文摘The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.
文摘For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.
基金supported by the National Natural Science Foundation of China(Nos.12175231 and 11805131)Anhui Natural Science Foundation of China(No.2108085J05)+1 种基金the National Key Research and Development Plan of China(No.2018YFE0307101)the Collaborative Innovation Program of the Hefei Science Center,CAS(Nos.2021HSC-CIP020 and 2022HSC-CIP009)。
文摘The helium bubbles induced by 14 MeV neutron irradiation can cause intergranular fractures in reduced activation ferritic martensitic steel,which is a candidate structural material for fusion reactors.In order to elucidate the susceptibility of different grain boundaries(GBs)to helium-induced embrittlement,the tensile fracture processes of 10 types of GBs with and without helium bubbles in body-centered cubic(bcc)iron at the relevant service temperature of 600 K were investigated via molecular dynamics methods.The results indicate that in the absence of helium bubbles,the GBs studied here can be classified into two distinct categories:brittle GBs and ductile GBs.The atomic scale analysis shows that the plastic deformation of ductile GB at high temperatures originates from complex plastic deformation mechanisms,including the Bain/Burgers path phase transition and deformation twinning,in which the Bain path phase transition is the most dominant plastic deformation mechanism.However,the presence of helium bubbles severely inhibits the plastic deformation channels of the GBs,resulting in a significant decrease in elongation at fractures.For bubble-decorated GBs,the ultimate tensile strength increases with the increase in the misorientation angle.Interestingly,the coherent twin boundary∑3{112}was found to maintain relatively high fracture strength and maximum failure strain under the influence of helium bubbles.
基金supported by the National Natural Science Foundation of China(11571132,12301542)the Natural Science Foundation of Hubei(2022CFB725)the Natural Science Foundation of Yichang(A23-2-027)。
文摘We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.