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.

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.

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.

Semi-analytical solution for drained expansion analysis of a hollow cylinder of critical state soils

The expansion of a thick-walled hollow cylinder in soil is of non-self-similar nature that the stress/deformation paths are not the same for different soil material points.As a result,this problem cannot be solved by the common self-similar-based similarity techniques.This paper proposes a novel,exact solution for rigorous drained expansion analysis of a hollow cylinder of critical state soils.Considering stress-dependent elastic moduli of soils,new analytical stress and displacement solutions for the nonself-similar problem are developed taking the small strain assumption in the elastic zone.In the plastic zone,the cavity expansion response is formulated into a set of first-order partial differential equations(PDEs)with the combination use of Eulerian and Lagrangian descriptions,and a novel solution algorithm is developed to efficiently solve this complex boundary value problem.The solution is presented in a general form and thus can be useful for a wide range of soils.With the new solution,the non-self-similar nature induced by the finite outer boundary is clearly demonstrated and highlighted,which is found to be greatly different to the behaviour of cavity expansion in infinite soil mass.The present solution may serve as a benchmark for verifying the performance of advanced numerical techniques with critical state soil models and be used to capture the finite boundary effect for pressuremeter tests in small-sized calibration chambers.

Solution Domain Boundary Analysis Method and Its Application in Parameter Spaces of Nonlinear Gear System

Mastering the influence laws of parameters on the solution structure of nonlinear systems is the basis of carrying out vibration isolation and control.Many researches on solution structure and bifurcation phenomenon in parameter spaces are carried out broadly in many fields,and the research on nonlinear gear systems has attracted the attention of many scholars.But there is little study on the solution domain boundary of nonlinear gear systems.For a periodic non-autonomous nonlinear dynamic system with several control parameters,a solution domain boundary analysis method of nonlinear systems in parameter spaces is proposed,which combines the cell mapping method based on Poincaré point mapping in phase spaces with the domain decomposition technique of parameter spaces.The cell mapping is known as a global analysis method to analyze the global behavior of a nonlinear dynamic system with finite dimensions,and the basic idea of domain decomposition techniques is to divide and rule.The method is applied to analyze the solution domain boundaries in parameter spaces of a nonlinear gear system.The distribution of different period domains,chaos domain and the domain boundaries between different period domains and chaotic domain are obtained in control parameter spaces constituted by meshing damping ratio with excitation frequency,fluctuation coefficient of meshing stiffness and average exciting force respectively by calculation.The calculation results show that as the meshing damping increases,the responses of the system change towards a single motion,while the variations of the excitation frequency,meshing stiffness and exciting force make the solution domain presenting diversity.The proposed research contribution provides evidence for vibration control and parameter design of the gear system,and confirms the validity of the solution domain boundary analysis method.

Time-Domain Analysis of Underground Station-Layered Soil Interaction Based on High-Order Doubly Asymptotic Transmitting Boundary

Based on the modified scale boundary finite element method and continued fraction solution,a high-order doubly asymptotic transmitting boundary(DATB)is derived and extended to the simulation of vector wave propagation in complex layered soils.The high-order DATB converges rapidly to the exact solution throughout the entire frequency range and its formulation is local in the time domain,possessing high accuracy and good efficiency.Combining with finite element method,a coupled model is constructed for time-domain analysis of underground station-layered soil interaction.The coupled model is divided into the near and far field by the truncated boundary,of which the near field is modelled by FEM while the far field is modelled by the high-order DATB.The coupled model is implemented in an open source finite element software,OpenSees,in which the DATB is employed as a super element.Numerical examples demonstrate that results of the coupled model are stable,accurate and efficient compared with those of the extended mesh model and the viscous-spring boundary model.Besides,it has also shown the fitness for long-time seismic response analysis of underground station-layered soil interaction.Therefore,it is believed that the coupled model could provide a new approach for seismic analysis of underground station-layered soil interaction and could be further developed for engineering.

An equivalent-boundary method for the shell analysis of buried pipelines under fault movement

A new shell finite element method (FEM) model with an equivalent boundary is presented for estimating the re- sponse of a buried pipeline under large fault movement. The length of affected pipeline under fault movement is usually too long for a shell-mode calculation because of the limitation of memory and time of computers. In this study, only the pipeline segment near fault is modeled with plastic shell elements to study the local buckling and the large section deformation in pipe. The material property of pipe segment far away from the fault is considered as elastic, and nonlinear spring elements at equivalent boundaries are obtained and applied to two ends of shell model. Compared with the fixed-boundary shell model, the shell model with an equivalent boundary proposed by the study can remarkably reduce the needed memory and calculating time.

Isogeometric Boundary Element Analysis for 2D Transient Heat Conduction Problem with Radial Integration Method

This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.

Ship Hull Slamming Analysis with Nonlinear Boundary Element Method

A boundary element method is developed for calculating the flare ship hull slamming problem. The nonlinear free surface elevation and the linear element assumption are employed. The method has been verified by comparisons with results for the water entry of wedges with various deadrise angles. Numerical results show that the pressure distribution varies greatly with the ship hull with different curvilinear equations, and the slamming features are also different; From the numerical simulation, the authors found that the structural damage of the flare hull might be caused by the increasing hydrodynamic pressure over an extensive area on the flare when the upper part of the flare comes into contact with water.

Hydraulic Fracture Parameter Inversion Method for Shale GasWells Based on Transient Pressure-Drop Analysis during Hydraulic Fracturing Shut-in Period

Horizontal well drilling and multi-stage hydraulic fracturing are key technologies for the development of shale gas reservoirs.Instantaneous acquisition of hydraulic fracture parameters is crucial for evaluating fracturing effectiveness,optimizing processes,and predicting gas productivity.This paper establishes a transient flow model for shale gas wells based on the boundary element method,achieving the characterization of stimulated reservoir volume for a single stage.By integrating pressure monitoring data following the pumping shut-in period of hydraulic fracturing for well testing interpretation,a workflow for inverting fracture parameters of shale gas wells is established.This new method eliminates the need for prolonged production testing and can interpret parameters of individual hydraulic fracture segments,offering significant advantages over the conventional pressure transient analysismethod.The practical application of thismethodology was conducted on 10 shale gaswellswithin the Changning shale gas block of Sichuan,China.The results show a high correlation between the interpreted single-stage total length and surface area of hydraulic fractures and the outcomes of gas production profile tests.Additionally,significant correlations are observed between these parameters and cluster number,horizontal stress difference,and natural fracture density.This demonstrates the effectiveness of the proposed fracture parameter inversion method and the feasibility of field application.The findings of this study aim to provide solutions and references for the inversion of fracture parameters in shale gas wells.

Axisymmetric Powell-Eyring fluid flow with convective boundary condition: optimal analysis

The effects of axisymmetric flow of a Powell-Eyring fluid over an impermeable radially stretching surface are presented. Characteristics of the heat transfer process are analyzed with a more realistic condition named the convective boundary condition. Governing equations for the flow problem are derived by the boundary layer approximations. The modeled highly coupled partial differential system is converted into a system of ordinary differential equations with acceptable similarity transformations. The convergent series solutions for the resulting system are constructed and analyzed. Optimal values are obtained and presented in a numerical form using an optimal homotopy analysis method （OHAM）. The rheological characteristics of different parameters of the velocity and temperature profiles are presented graphically. Tabular variations of the skin friction coefficient and the Nusselt number are also calculated. It is observed that the temperature distribution shows opposite behavior for Prandtl and Biot numbers. Furthermore, the rate of heating/cooling is higher for both the Prandtl and Biot numbers.

Optimal Homotopy Analysis of MHD Natural Convection Flow of Thixotropic Fluid under Subjection of Thermal Stratification: Boundary Layer Analysis

In this article, the model of a non-Newtonian fluid (Thixotropic) flow past a vertical surface in the presence of exponential space and temperature dependent heat source in a thermally stratified medium is studied. It is assumed that free convection is induced by buoyancy and exponentially decaying internal heat source across the space. The dynamic viscosity is taken to be constant and thermal conductivity of this particular fluid model is assumed to vary linearly with temperature. Thermal stratification has been properly incorporated into the governing equation so that its effect can be revealed and properly reported. The governing partial differential equations describing the model are transformed and parameterized to a system of non-linear ordinary differential equation using similarity transformations. Approximate analytic solutions were obtained by adopting Optimal Homotopy Analysis Method (OHAM). The results show that for both cases of non-Newtonian parameters (Thixotropic) (K1=K2=0?& K1=K2=1.0), increasing stratification parameters, relate to decreasing in the heat energy entering into the fluid region and thus reducing the temperature of the Thixotropic fluid as it flows.

Resolving Domain Integral Issues in Isogeometric Boundary Element Methods via Radial Integration:A Study of Thermoelastic Analysis

The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.

Boundary Element Analysis (Laplace Transform Solution) of Groundwater Unsteady Flow to a Multiple Well System in a Confined Aquifer

The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.

Second-Order Adjoint Sensitivity Analysis Methodology for Computing Exactly Response Sensitivities to Uncertain Parameters and Boundaries of Linear Systems: Mathematical Framework

This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2nd-CASAM)” for the efficient and exact computation of 1st- and 2nd-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (i.e., model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2nd-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1st-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2nd-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (N2/2 + 3 N/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2nd-LASS requires very little additional effort beyond the construction of the 1st-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1st-LASS and 2nd-LASS requires the same computational solvers as needed for solving (i.e., “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1st-LASS and the 2nd-LASS. Since neither the 1st-LASS nor the 2nd-LASS involves any differentials of the operators underlying the original system, the 1st-LASS is designated as a “first-level” (as opposed to a “first-order”) adjoint sensitivity system, while the 2nd-LASS is designated as a “second-level” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2nd-LASS that involve the imprecisely known boundary parameters. Notably, the 2nd-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.

AES Analysis of Ce and P Cosegregation at Grain Boundaries

The chemical state of grain boundary cosegregation of Ce and P in Fe-P-Ce alloy system was studied by means of Auger Electron Spectroscopy(AES).The Auger peaks of Ce segregated at grain boundaries are found within the range of 60～180 eV.By comparing with the Auger spectra of the Fe-Ce-P intermetallic compound, it is supposed that there is a two-dimensional interfacial phase at grain boundaries with Ce and P cosegregation which is similar to the structure of the Fe-Ce-P compound.

STOCHASTIC BOUNDARY ELEMENT METHODS FOR 3－D PROBLEMS WITH CENTRIFUGAL FORCES AND RELIABILITY ANALYSIS

A stochastic boundary element method (SBEM) is developed for 3-Dproblems with body forces. The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables, considering the strengthlimit, rotation speeds and material density to be the fundamental stochastic variables.The method developed is applied to analyzing the strength reliability of the turbo diskof an aero-engine.

LAPLACE TRANSFORM-BOUNDARY ELEMENT COUPLING METHOD FOR VISCOUS FLUIDSTRUCTURE IMPACT ANALYSIS

Based on linearized 2-D Navier-Stokes equation, a Laplacetransform-boundary element coupling method for viscousfluid-structure impact analysis is proposed. Under assumption ofincompressibility for the fluid, the corresponding equivalentboundary integral equation in terms of the potential function andstream function is first established by Lamb's transform in theLaplace transform domain.

THE STOCHASTIC BOUNDARY ELEMENT METHOD IN STATISTICAL ANALYSIS OF MODERATELY THICK PLATES

In this paper a stochastic boundary element method (SEEM) is developed to analyze moderately thick plates with random material parameters and random thickness. Based on the Taylor series expansion, the boundary integration equations concerning the mean and deviation of the generalized displacements are derived, respectively. It is found that the randomness of material parameters is equivalent to a random load, so the mean and covariance matrices of unknown generalized boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of generalized displacements and forces at internal points can also be obtained. A numerical example has been worked out with the method proposed and necessary analysis is made for the results.

Bifurcation Analysis for a Free Boundary Problem Modeling Growth of Solid Tumor with Inhibitors

This paper is concerned with the bifurcation analysis for a free boundary problem modeling the growth of solid tumor with inhibitors.In this problem,surface tension coefficient plays the role of bifurcation parameter,it is proved that there exists a sequence of the nonradially stationary solutions bifurcate from the radially symmetric stationary solutions.Our results indicate that the tumor grown in vivo may have various shapes.In particular,a tumor with an inhibitor is associated with the growth of protrusions.