GLOBAL SOLUTIONS TO 1D COMPRESSIBLE NAVIER-STOKES/ALLEN-CAHN SYSTEM WITH DENSITY-DEPENDENT VISCOSITY AND FREE-BOUNDARY

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.

Comparison of perfectly matched layer and multi-transmitting formula artificial boundary condition based on hybrid finite element formulation

The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.

Dirac method for nonlinear and non-homogenous boundary value problems of plates

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.

Grain boundary engineering for enhancing intergranular damage resistance of ferritic/martensitic steel P92

Ferritic/martensitic(F/M)steel is widely used as a structural material in thermal and nuclear power plants.However,it is susceptible to intergranular damage,which is a critical issue,under service conditions.In this study,to improve the resistance to intergranular damage of F/M steel,a thermomechanical process(TMP)was employed to achieve a grain boundary engineering(GBE)microstructure in F/M steel P92.The TMP,including cold-rolling thickness reduction of 6%,9%,and 12%,followed by austenitization at 1323 K for 40 min and tempering at 1053 K for 45 min,was applied to the as-received(AR)P92 steel.The prior austenite grain(PAG)size,prior austenite grain boundary character distribution(GBCD),and connectivity of prior austenite grain boundaries(PAGBs)were investigated.Compared to the AR specimen,the PAG size did not change significantly.The fraction of coincident site lattice boundaries(CSLBs,3≤Σ≤29)and Σ3^(n) boundaries along PAGBs decreased with increasing reduction ratio because the recrystallization fraction increased with increasing reduction ratio.The PAGB connectivity of the 6%deformed specimen slightly deteriorated compared with that of the AR specimen.Moreover,potentiodynamic polarization studies revealed that the intergranular damage resistance of the studied steel could be improved by increasing the fraction of CSLBs along the PAGBs,indicating that the TMP,which involves low deformation,could enhance the intergranular damage resistance.

A High-Accuracy Curve Boundary Recognition Method Based on the Lattice Boltzmann Method and Immersed Moving Boundary Method

Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

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.

Experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering boundary effect

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.

Development of D_α band symmetrical visible optical diagnostic for boundary reconstruction on EAST tokamak

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.

CONVEXITY OF THE FREE BOUNDARY FOR AN AXISYMMETRIC INCOMPRESSIBLE IMPINGING JET

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.

Two Stages Segmentation Algorithm of Breast Tumor in DCE-MRI Based on Multi-Scale Feature and Boundary Attention Mechanism

Nuclearmagnetic resonance imaging of breasts often presents complex backgrounds.Breast tumors exhibit varying sizes,uneven intensity,and indistinct boundaries.These characteristics can lead to challenges such as low accuracy and incorrect segmentation during tumor segmentation.Thus,we propose a two-stage breast tumor segmentation method leveraging multi-scale features and boundary attention mechanisms.Initially,the breast region of interest is extracted to isolate the breast area from surrounding tissues and organs.Subsequently,we devise a fusion network incorporatingmulti-scale features and boundary attentionmechanisms for breast tumor segmentation.We incorporate multi-scale parallel dilated convolution modules into the network,enhancing its capability to segment tumors of various sizes through multi-scale convolution and novel fusion techniques.Additionally,attention and boundary detection modules are included to augment the network’s capacity to locate tumors by capturing nonlocal dependencies in both spatial and channel domains.Furthermore,a hybrid loss function with boundary weight is employed to address sample class imbalance issues and enhance the network’s boundary maintenance capability through additional loss.Themethod was evaluated using breast data from 207 patients at RuijinHospital,resulting in a 6.64%increase in Dice similarity coefficient compared to the benchmarkU-Net.Experimental results demonstrate the superiority of the method over other segmentation techniques,with fewer model parameters.

Effects of Fe solid solute on grain boundaries of bi-crystal Cu: A molecular dynamics simulation

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.

Modeling the Interaction between Vacancies and Grain Boundaries during Ductile Fracture

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.

Delineation of urban growth boundary based on FLUS model under the perspective of land use evaluation in hilly mountainous areas

With rapid economic development,the size of urban land in China is expanding dramatically.The Urban Growth Boundary(UGB)is an expandable spatial boundary for urban construction in a certain period in order to control the urban sprawl.Reasonable delineation of UGB can inhibit the disorderly spread of urban space and guide the normal development of the city.It is of practical significance for the construction of green urban space.The study utilizes GIS technology to establish a land construction suitability evaluation system for Nankang city,which is experiencing rapid urban expansion,and outlines the preliminary UGB under the future land use simulation(FLUS)model.At the same time,considering the coupled coordination of"Production-Living-Ecological Space",and based on the suitability evaluation,we revised the preliminary UGB by combining the advantages of the patch-generating land use simulation(PLUS)model and the convex hull model to delineate the final UGB.The results show that:1)the comprehensive score of the evaluation of the suitability of the construction of land from high to low shows the distribution of the center of the city to the surrounding circle type spread,the center of the city has the highest suitability score.The results of convex hull model show that the urban expansion type of Nankang is epitaxial.In the future,the urban expansion will mainly occur in the northern part of the city.The PLUS model predicts an increase of 3359.97 hm^(2)of construction land in Nankang by 2035,of which 2022.97 hm^(2)is urban construction land.2)The FLUS model has a prediction accuracy of 86.3%and delineates a preliminary UGB area of 9215.07 hm^(2).3)We used the results of the construction suitability evaluation,PLUS model simulation results,and convex hull model predictions to revise the originally delineated UGB.The final delineated UGB area is 8895.67 hm^(2)and it is capable of meeting the future development of the study area.The results of the delineation can promote sustainable urban development,and the delineation methodology can provide a reference basis for the preparation of territorial spatial planning.

Analytical solutions of turbulent boundary layer beneath forward-leaning waves

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 Boundary Element Method for Ordinary State-Based Peridynamics

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.

Multi-scale physics-informed neural networks for solving high Reynolds number boundary layer flows based on matched asymptotic expansions

Multi-scale system remains a classical scientific problem in fluid dynamics,biology,etc.In the present study,a scheme of multi-scale Physics-informed neural networks is proposed to solve the boundary layer flow at high Reynolds numbers without any data.The flow is divided into several regions with different scales based on Prandtl's boundary theory.Different regions are solved with governing equations in different scales.The method of matched asymptotic expansions is used to make the flow field continuously.A flow on a semi infinite flat plate at a high Reynolds number is considered a multi-scale problem because the boundary layer scale is much smaller than the outer flow scale.The results are compared with the reference numerical solutions,which show that the msPINNs can solve the multi-scale problem of the boundary layer in high Reynolds number flows.This scheme can be developed for more multi-scale problems in the future.

Enriched Constant Elements in the Boundary Element Method for Solving 2D Acoustic Problems at Higher Frequencies

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.

Generalized nth-Order Perturbation Method Based on Loop Subdivision Surface Boundary Element Method for Three-Dimensional Broadband Structural Acoustic Uncertainty Analysis

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.

CAUCHY TYPE INTEGRALS AND A BOUNDARY VALUE PROBLEM IN A COMPLEX CLIFFORD ANALYSIS

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.

Effect of boundary conditions on shakedown analysis of heterogeneous materials

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.