We investigate the following inverse problem:starting from the acoustic wave equation,reconstruct a piecewise constant passive acoustic source from a single boundary temporal measurement without knowing the speed of s...We investigate the following inverse problem:starting from the acoustic wave equation,reconstruct a piecewise constant passive acoustic source from a single boundary temporal measurement without knowing the speed of sound.When the amplitudes of the source are known a priori,we prove a unique determination result of the shape and propose a level set algorithm to reconstruct the singularities.When the singularities of the source are known a priori,we show unique determination of the source amplitudes and propose a least-squares fitting algorithm to recover the source amplitudes.The analysis bridges the low-frequency source inversion problem and the inverse problem of gravimetry.The proposed algorithms are validated and quantitatively evaluated with numerical experiments in 2D and 3D.展开更多
This paper proposes a new level-set-based shape recovery approach that can be applied to a wide range of binary tomography reconstructions.In this technique,we derive generic evolution equations for shape reconstructi...This paper proposes a new level-set-based shape recovery approach that can be applied to a wide range of binary tomography reconstructions.In this technique,we derive generic evolution equations for shape reconstruction in terms of the underlying level-set parameters.We show that using the appropriate basis function to parameterize the level-set function results in an optimization problem with a small number of parameters,which overcomes many of the problems associated with the traditional level-set approach.More concretely,in this paper,we use Gaussian functions as a basis function placed at sparse grid points to represent the parametric level-set function and provide more flexibility in the binary representation of the reconstructed image.In addition,we suggest a convex optimization method that can overcome the problem of the local minimum of the cost function by successfully recovering the coefficients of the basis function.Finally,we illustrate the performance of the proposed method using synthetic images and real X-ray CT projection data.We show that the proposed reconstruction method compares favorably to various state-of-the-art reconstruction techniques for limited-data tomography,and it is also relatively stable in the presence of modest amounts of noise.Furthermore,the shape representation using a compact Gaussian radial basis function works well.展开更多
A novel approach to extract flame fronts, which is called the conditioned level-set method with block division （CLSB）, has been developed. Based on a two-phase level-set formulation, the conditioned initialization a...A novel approach to extract flame fronts, which is called the conditioned level-set method with block division （CLSB）, has been developed. Based on a two-phase level-set formulation, the conditioned initialization and region-lock optimiza-tion appear to be beneficial to improve the efficiency and accuracy of the flame contour identification. The original block- division strategy enables the approach to be unsupervised by calculating local self-adaptive threshold values autonomously before binarization. The CLSB approach has been applied to deal with a large set of experimental data involving swirl- stabilized premixed combustion in diluted regimes operating at atmospheric pressures. The OH-PLIF measurements have been carried out in this framework. The resulting images are, thus, featured by lower signal-to-noise ratios （SNRs） than the ideal image; relatively complex flame structures lead to significant non-uniformity in the OH signal intensity; and, the mag- nitude of the maximum OH gradient observed along the flame front can also vary depending on flow or local stoichiometry. Compared with other conventional edge detection operators, the CLSB method demonstrates a good ability to deal with the OH-PLIF images at low SNR and with the presence of a multiple scales of both OH intensity and OH gradient. The robustness to noise sensitivity and intensity inhomogeneity has been evaluated throughout a range of experimental images of diluted flames, as well as against a circle test as Ground Truth （GT）.展开更多
This paper presents an efficient liver-segmentation system developed by combining three ideas under the operations of a level-set method and consequent processes. First, an effective initial process creates mask and s...This paper presents an efficient liver-segmentation system developed by combining three ideas under the operations of a level-set method and consequent processes. First, an effective initial process creates mask and seed regions. The mask regions assist in prevention of leakage regions due to an overlap of gray-intensities between liver and another soft-tissue around ribs and verte-brae. The seed regions are allocated inside the liver to measure statistical values of its gray-intensities. Second, we introduce liver-corrective images to represent statistical regions of the liver and preserve edge information. These images help a geodesic active contour (GAC) to move without obstruction from high level of image noises. Lastly, the computation time in a level-set based on reaction-diffusion evolution and the GAC method is reduced by using a concept of multi-resolution. We applied the proposed system to 40 sets of 3D CT-liver data, which were acquired from four patients (10 different sets per patient) by a 4D-CT imaging system. The segmentation results showed 86.38% ± 4.26% (DSC: 91.38% ± 2.99%) of similarities to outlines of manual delineation provided by a radiologist. Meanwhile, the results of liver segmentation only using edge images presented 79.17% ± 5.15% or statistical regions showed 74.04% ± 9.77% of similarities.展开更多
The separation-of-variable(SOV)methods,such as the improved SOV method,the variational SOV method,and the extended SOV method,have been proposed by the present authors and coworkers to obtain the closed-form analytica...The separation-of-variable(SOV)methods,such as the improved SOV method,the variational SOV method,and the extended SOV method,have been proposed by the present authors and coworkers to obtain the closed-form analytical solutions for free vibration and eigenbuckling of rectangular plates and circular cylindrical shells.By taking the free vibration of rectangular thin plates as an example,this work presents the theoretical framework of the SOV methods in an instructive way,and the bisection–based solution procedures for a group of nonlinear eigenvalue equations.Besides,the explicit equations of nodal lines of the SOV methods are presented,and the relations of nodal line patterns and frequency orders are investigated.It is concluded that the highly accurate SOV methods have the same accuracy for all frequencies,the mode shapes about repeated frequencies can also be precisely captured,and the SOV methods do not have the problem of missing roots as well.展开更多
This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydra...This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.展开更多
To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
In the independent electro-hydrogen system(IEHS)with hybrid energy storage(HESS),achieving optimal scheduling is crucial.Still,it presents a challenge due to the significant deviations in values ofmultiple optimizatio...In the independent electro-hydrogen system(IEHS)with hybrid energy storage(HESS),achieving optimal scheduling is crucial.Still,it presents a challenge due to the significant deviations in values ofmultiple optimization objective functions caused by their physical dimensions.These deviations seriously affect the scheduling process.A novel standardization fusion method has been established to address this issue by analyzing the variation process of each objective function’s values.The optimal scheduling results of IEHS with HESS indicate that the economy and overall energy loss can be improved 2–3 times under different optimization methods.The proposed method better balances all optimization objective functions and reduces the impact of their dimensionality.When the cost of BESS decreases by approximately 30%,its participation deepens by about 1 time.Moreover,if the price of the electrolyzer is less than 15¥/kWh or if the cost of the fuel cell drops below 4¥/kWh,their participation will increase substantially.This study aims to provide a more reasonable approach to solving multi-objective optimization problems.展开更多
In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least s...In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.展开更多
Understanding the wind power potential of a site is essential for designing an optimal wind power conditioning system. The Weibull distribution and wind speed extrapolation methods are powerful mathematical tools for ...Understanding the wind power potential of a site is essential for designing an optimal wind power conditioning system. The Weibull distribution and wind speed extrapolation methods are powerful mathematical tools for efficiently predicting the frequency distribution of wind speeds at a site. Hourly wind speed and direction data were collected from the National Aeronautics and Space Administration (NASA) website for the period 2013 to 2023. MATLAB software was used to calculate the distribution parameters using the graphical method and to plot the corresponding curves, while WRPLOTView software was used to construct the wind rose. The average wind speed obtained is 3.33 m/s and can reach up to 5.71 m/s at a height of 100 meters. The wind energy is estimated to be 1315.30 kWh/m2 at a height of 100 meters. The wind rose indicates the prevailing winds (ranging from 3.60 m/s to 5.70 m/s) in the northeast-east direction.展开更多
Modal parameters can accurately characterize the structural dynamic properties and assess the physical state of the structure.Therefore,it is particularly significant to identify the structural modal parameters accordi...Modal parameters can accurately characterize the structural dynamic properties and assess the physical state of the structure.Therefore,it is particularly significant to identify the structural modal parameters according to the monitoring data information in the structural health monitoring(SHM)system,so as to provide a scientific basis for structural damage identification and dynamic model modification.In view of this,this paper reviews methods for identifying structural modal parameters under environmental excitation and briefly describes how to identify structural damages based on the derived modal parameters.The paper primarily introduces data-driven modal parameter recognition methods(e.g.,time-domain,frequency-domain,and time-frequency-domain methods,etc.),briefly describes damage identification methods based on the variations of modal parameters(e.g.,natural frequency,modal shapes,and curvature modal shapes,etc.)and modal validation methods(e.g.,Stability Diagram and Modal Assurance Criterion,etc.).The current status of the application of artificial intelligence(AI)methods in the direction of modal parameter recognition and damage identification is further discussed.Based on the pre-vious analysis,the main development trends of structural modal parameter recognition and damage identification methods are given to provide scientific references for the optimized design and functional upgrading of SHM systems.展开更多
Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
The prediction of the fracture plane orientation in fatigue is a scientific topic and remains relevant for every type of material. However, in this work, we compared the orientation of the fracture plane obtained expe...The prediction of the fracture plane orientation in fatigue is a scientific topic and remains relevant for every type of material. However, in this work, we compared the orientation of the fracture plane obtained experimentally through tests on specimens under multiaxial loading with that calculated by the variance method. In the statistical approach criteria, several methods have been developed but we have presented only one method, namely the variance method using the equivalent stress. She assumes that the fracture plane orientation is the one on which the variance of the equivalent stress is maximum. Three types of equivalent stress are defined for this method [1]: normal stress, shear stress and combined normal and shear stress. The results obtained were compared with experimental results for multiaxial cyclic stress states, and it emerges that the variance method for the case of combined loading is conservative as it gives a better prediction of the fracture plane.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using in...The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.展开更多
Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extens...Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extensive attention due to their advantages of small size,high multiplexing efficiency,convenient mass production,and low cost.An RGB beam combiner based on directional couplers is designed,with a core-cladding relative refractive index difference of 0.75%.The RGB beam combiner is optimized from the perspective of parameter optimization.Using the beam propagation method(BPM),the relationship between the performance of the RGB beam combiner and individual parameters is studied,achieving preliminary optimization of the device’s performance.The key parameters of the RGB beam combiner are optimized using the entropy weight-technique for order preference by similarity to an ideal solution TOPSIS method,establishing the optimal parameter scheme and further improving the device’s performance indicators.The results show that after optimization,the multiplexing efficiencies for red,green,and blue lights,as well as the average multiplexing efficiency,reached 99.17%,99.76%,96.63%and 98.52%,respectively.The size of the RGB beam combiner is 4.768 mm×0.062 mm.展开更多
Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump mate...Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.展开更多
Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, et...Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, etc. Therefore, it is of great significance to study fractional-order time-delay differential equations. In this paper, we discuss a finite volume element method for a class of fractional-order neutral time-delay differential equations. By introducing an intermediate variable, the fourth-order problem is transformed into a system of equations consisting of two second-order partial differential equations. The L1 formula is used to approximate the time fractional order derivative terms, and the finite volume element method is used in space. A fully discrete format of the equations is established, and we prove the existence, uniqueness, convergence and stability of the solution. Finally, the validity of the format is verified by numerical examples.展开更多
基金partially supported by the NSF(Grant Nos.2012046,2152011,and 2309534)partially supported by the NSF(Grant Nos.DMS-1715178,DMS-2006881,and DMS-2237534)+1 种基金NIH(Grant No.R03-EB033521)startup fund from Michigan State University.
文摘We investigate the following inverse problem:starting from the acoustic wave equation,reconstruct a piecewise constant passive acoustic source from a single boundary temporal measurement without knowing the speed of sound.When the amplitudes of the source are known a priori,we prove a unique determination result of the shape and propose a level set algorithm to reconstruct the singularities.When the singularities of the source are known a priori,we show unique determination of the source amplitudes and propose a least-squares fitting algorithm to recover the source amplitudes.The analysis bridges the low-frequency source inversion problem and the inverse problem of gravimetry.The proposed algorithms are validated and quantitatively evaluated with numerical experiments in 2D and 3D.
基金This work was supported by JST-CREST Grant Number JPMJCR1765,Japan.
文摘This paper proposes a new level-set-based shape recovery approach that can be applied to a wide range of binary tomography reconstructions.In this technique,we derive generic evolution equations for shape reconstruction in terms of the underlying level-set parameters.We show that using the appropriate basis function to parameterize the level-set function results in an optimization problem with a small number of parameters,which overcomes many of the problems associated with the traditional level-set approach.More concretely,in this paper,we use Gaussian functions as a basis function placed at sparse grid points to represent the parametric level-set function and provide more flexibility in the binary representation of the reconstructed image.In addition,we suggest a convex optimization method that can overcome the problem of the local minimum of the cost function by successfully recovering the coefficients of the basis function.Finally,we illustrate the performance of the proposed method using synthetic images and real X-ray CT projection data.We show that the proposed reconstruction method compares favorably to various state-of-the-art reconstruction techniques for limited-data tomography,and it is also relatively stable in the presence of modest amounts of noise.Furthermore,the shape representation using a compact Gaussian radial basis function works well.
文摘A novel approach to extract flame fronts, which is called the conditioned level-set method with block division （CLSB）, has been developed. Based on a two-phase level-set formulation, the conditioned initialization and region-lock optimiza-tion appear to be beneficial to improve the efficiency and accuracy of the flame contour identification. The original block- division strategy enables the approach to be unsupervised by calculating local self-adaptive threshold values autonomously before binarization. The CLSB approach has been applied to deal with a large set of experimental data involving swirl- stabilized premixed combustion in diluted regimes operating at atmospheric pressures. The OH-PLIF measurements have been carried out in this framework. The resulting images are, thus, featured by lower signal-to-noise ratios （SNRs） than the ideal image; relatively complex flame structures lead to significant non-uniformity in the OH signal intensity; and, the mag- nitude of the maximum OH gradient observed along the flame front can also vary depending on flow or local stoichiometry. Compared with other conventional edge detection operators, the CLSB method demonstrates a good ability to deal with the OH-PLIF images at low SNR and with the presence of a multiple scales of both OH intensity and OH gradient. The robustness to noise sensitivity and intensity inhomogeneity has been evaluated throughout a range of experimental images of diluted flames, as well as against a circle test as Ground Truth （GT）.
文摘This paper presents an efficient liver-segmentation system developed by combining three ideas under the operations of a level-set method and consequent processes. First, an effective initial process creates mask and seed regions. The mask regions assist in prevention of leakage regions due to an overlap of gray-intensities between liver and another soft-tissue around ribs and verte-brae. The seed regions are allocated inside the liver to measure statistical values of its gray-intensities. Second, we introduce liver-corrective images to represent statistical regions of the liver and preserve edge information. These images help a geodesic active contour (GAC) to move without obstruction from high level of image noises. Lastly, the computation time in a level-set based on reaction-diffusion evolution and the GAC method is reduced by using a concept of multi-resolution. We applied the proposed system to 40 sets of 3D CT-liver data, which were acquired from four patients (10 different sets per patient) by a 4D-CT imaging system. The segmentation results showed 86.38% ± 4.26% (DSC: 91.38% ± 2.99%) of similarities to outlines of manual delineation provided by a radiologist. Meanwhile, the results of liver segmentation only using edge images presented 79.17% ± 5.15% or statistical regions showed 74.04% ± 9.77% of similarities.
基金supported by the National Natural Science Foundation of China(12172023).
文摘The separation-of-variable(SOV)methods,such as the improved SOV method,the variational SOV method,and the extended SOV method,have been proposed by the present authors and coworkers to obtain the closed-form analytical solutions for free vibration and eigenbuckling of rectangular plates and circular cylindrical shells.By taking the free vibration of rectangular thin plates as an example,this work presents the theoretical framework of the SOV methods in an instructive way,and the bisection–based solution procedures for a group of nonlinear eigenvalue equations.Besides,the explicit equations of nodal lines of the SOV methods are presented,and the relations of nodal line patterns and frequency orders are investigated.It is concluded that the highly accurate SOV methods have the same accuracy for all frequencies,the mode shapes about repeated frequencies can also be precisely captured,and the SOV methods do not have the problem of missing roots as well.
基金Funded by the Natural Science Foundation of Jiangsu Province(No.BK20241529)China Postdoctoral Science Foundation(No.2024M750736)。
文摘This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金sponsored by R&D Program of Beijing Municipal Education Commission(KM202410009013).
文摘In the independent electro-hydrogen system(IEHS)with hybrid energy storage(HESS),achieving optimal scheduling is crucial.Still,it presents a challenge due to the significant deviations in values ofmultiple optimization objective functions caused by their physical dimensions.These deviations seriously affect the scheduling process.A novel standardization fusion method has been established to address this issue by analyzing the variation process of each objective function’s values.The optimal scheduling results of IEHS with HESS indicate that the economy and overall energy loss can be improved 2–3 times under different optimization methods.The proposed method better balances all optimization objective functions and reduces the impact of their dimensionality.When the cost of BESS decreases by approximately 30%,its participation deepens by about 1 time.Moreover,if the price of the electrolyzer is less than 15¥/kWh or if the cost of the fuel cell drops below 4¥/kWh,their participation will increase substantially.This study aims to provide a more reasonable approach to solving multi-objective optimization problems.
基金supported by the Graduate Student Scientific Research Innovation Project through Research Innovation Fund for Graduate Students in Shanxi Province(Project No.2024KY648).
文摘In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.
文摘Understanding the wind power potential of a site is essential for designing an optimal wind power conditioning system. The Weibull distribution and wind speed extrapolation methods are powerful mathematical tools for efficiently predicting the frequency distribution of wind speeds at a site. Hourly wind speed and direction data were collected from the National Aeronautics and Space Administration (NASA) website for the period 2013 to 2023. MATLAB software was used to calculate the distribution parameters using the graphical method and to plot the corresponding curves, while WRPLOTView software was used to construct the wind rose. The average wind speed obtained is 3.33 m/s and can reach up to 5.71 m/s at a height of 100 meters. The wind energy is estimated to be 1315.30 kWh/m2 at a height of 100 meters. The wind rose indicates the prevailing winds (ranging from 3.60 m/s to 5.70 m/s) in the northeast-east direction.
基金supported by the Innovation Foundation of Provincial Education Department of Gansu(2024B-005)the Gansu Province National Science Foundation(22YF7GA182)the Fundamental Research Funds for the Central Universities(No.lzujbky2022-kb01)。
文摘Modal parameters can accurately characterize the structural dynamic properties and assess the physical state of the structure.Therefore,it is particularly significant to identify the structural modal parameters according to the monitoring data information in the structural health monitoring(SHM)system,so as to provide a scientific basis for structural damage identification and dynamic model modification.In view of this,this paper reviews methods for identifying structural modal parameters under environmental excitation and briefly describes how to identify structural damages based on the derived modal parameters.The paper primarily introduces data-driven modal parameter recognition methods(e.g.,time-domain,frequency-domain,and time-frequency-domain methods,etc.),briefly describes damage identification methods based on the variations of modal parameters(e.g.,natural frequency,modal shapes,and curvature modal shapes,etc.)and modal validation methods(e.g.,Stability Diagram and Modal Assurance Criterion,etc.).The current status of the application of artificial intelligence(AI)methods in the direction of modal parameter recognition and damage identification is further discussed.Based on the pre-vious analysis,the main development trends of structural modal parameter recognition and damage identification methods are given to provide scientific references for the optimized design and functional upgrading of SHM systems.
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
文摘The prediction of the fracture plane orientation in fatigue is a scientific topic and remains relevant for every type of material. However, in this work, we compared the orientation of the fracture plane obtained experimentally through tests on specimens under multiaxial loading with that calculated by the variance method. In the statistical approach criteria, several methods have been developed but we have presented only one method, namely the variance method using the equivalent stress. She assumes that the fracture plane orientation is the one on which the variance of the equivalent stress is maximum. Three types of equivalent stress are defined for this method [1]: normal stress, shear stress and combined normal and shear stress. The results obtained were compared with experimental results for multiaxial cyclic stress states, and it emerges that the variance method for the case of combined loading is conservative as it gives a better prediction of the fracture plane.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金funded by Vice Chancellor of Research at Shiraz University(grant 3GFU2M1820).
文摘The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.
基金Project(52175445)supported by the National Natural Science Foundation of ChinaProject(2022JJ30743)supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(2023GK2024)supported by the Key Research and Development Program of Hunan Province,ChinaProject(2023ZZTS0391)supported by the Fundamental Research Funds for the Central Universities of China。
文摘Red-green-blue(RGB)beam combiners are widely used in scenarios such as augmented reality/virtual reality(AR/VR),laser projection,biochemical detection,and other fields.Optical waveguide combiners have attracted extensive attention due to their advantages of small size,high multiplexing efficiency,convenient mass production,and low cost.An RGB beam combiner based on directional couplers is designed,with a core-cladding relative refractive index difference of 0.75%.The RGB beam combiner is optimized from the perspective of parameter optimization.Using the beam propagation method(BPM),the relationship between the performance of the RGB beam combiner and individual parameters is studied,achieving preliminary optimization of the device’s performance.The key parameters of the RGB beam combiner are optimized using the entropy weight-technique for order preference by similarity to an ideal solution TOPSIS method,establishing the optimal parameter scheme and further improving the device’s performance indicators.The results show that after optimization,the multiplexing efficiencies for red,green,and blue lights,as well as the average multiplexing efficiency,reached 99.17%,99.76%,96.63%and 98.52%,respectively.The size of the RGB beam combiner is 4.768 mm×0.062 mm.
基金the financial support provided by MHRD,Govt.of IndiaCoal India Limited for providing financial assistance for the research(Project No.CIL/R&D/01/73/2021)the partial financial support provided by the Ministry of Education,Government of India,under SPARC project(Project No.P1207)。
文摘Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.
文摘Fractional-order time-delay differential equations can describe many complex physical phenomena with memory or delay effects, which are widely used in the fields of cell biology, control systems, signal processing, etc. Therefore, it is of great significance to study fractional-order time-delay differential equations. In this paper, we discuss a finite volume element method for a class of fractional-order neutral time-delay differential equations. By introducing an intermediate variable, the fourth-order problem is transformed into a system of equations consisting of two second-order partial differential equations. The L1 formula is used to approximate the time fractional order derivative terms, and the finite volume element method is used in space. A fully discrete format of the equations is established, and we prove the existence, uniqueness, convergence and stability of the solution. Finally, the validity of the format is verified by numerical examples.
