Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation o...Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.展开更多
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.展开更多
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.展开更多
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.展开更多
In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery...Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
This paper presents a comparative study of a meshless level-set method in the simulation of sloshing flows. The numerical moving particle semi-implicit （MPS） method and a grid based schemes of the MPS and level-set ...This paper presents a comparative study of a meshless level-set method in the simulation of sloshing flows. The numerical moving particle semi-implicit （MPS） method and a grid based schemes of the MPS and level-set methods are outlined and two violent sloshing cases are considered. The computed results are compared with the corresponding experimental data for validation. The impact pressure and the deformations of free surface induced by sloshing are comparatively analyzed, and are in good agreement with experimental ones. Results show that both the MPS and level-set methods are good tools for simulation of violent sloshing flows. However, the second pressure peaks as well as breaking and splashing of free surface by the MPS method are captured better than by the level-set method.展开更多
We present and compare twoflexible and effective methodologies to predict disturbance zones ahead of underground tunnels by using elastic full-waveform inversion.One methodology uses a linearized,iterative approach bas...We present and compare twoflexible and effective methodologies to predict disturbance zones ahead of underground tunnels by using elastic full-waveform inversion.One methodology uses a linearized,iterative approach based on misfit gradients computed with the adjoint method while the other uses iterative,gradient-free unscented Kalmanfiltering in conjunction with a level-set representation.Whereas the former does not involve a priori assumptions on the distribution of elastic properties ahead of the tunnel,the latter intro-duces a massive reduction in the number of explicit model parameters to be inverted for by focusing on the geometric form of potential disturbances and their average elastic properties.Both imaging methodologies are validated through successful reconstructions of simple disturbances.As an application,we consider an elastic multiple disturbance scenario.By using identical synthetic time-domain seismo-grams as test data,we obtain satisfactory,albeit different,reconstruction results from the two inversion methodologies.The computa-tional costs of both approaches are of the same order of magnitude,with the gradient-based approach showing a slight advantage.The model parameter space reduction approach compensates for this by additionally providing a posteriori estimates of model parameter uncertainty.展开更多
A new Semi-Lagrangian scheme is proposed to discretize the surface convection-diffusion equation.The other involved equations including the the levelset convection equation,the re-initialization equation and the exten...A new Semi-Lagrangian scheme is proposed to discretize the surface convection-diffusion equation.The other involved equations including the the levelset convection equation,the re-initialization equation and the extension equation are also solved by S-L schemes.The S-L method removes both the CFL condition and the stiffness caused by the surface Laplacian,allowing larger time step than the Eulerian method.The method is extended to the block-structured adaptive mesh.Numerical examples are given to demonstrate the efficiency of the S-L method.展开更多
Semi-Lagrangian(S-L)methods have no CFL stability constraint,and are more stable than the Eulerian methods.In the literature,the S-L method for the levelset re-initialization equation was complicated,which may be unne...Semi-Lagrangian(S-L)methods have no CFL stability constraint,and are more stable than the Eulerian methods.In the literature,the S-L method for the levelset re-initialization equation was complicated,which may be unnecessary.Since the re-initialization procedure is auxiliary,we propose to use the first-order S-L scheme coupled with a projection technique to improve the accuracy at the grid points just adjacent to the interface.Standard second-order S-L method is used for evolving the level-set convection equation.The implementation is simple,including on the block-structured adaptive mesh.The efficiency of the S-L method is demonstrated by extensive numerical examples including passive convection of interfaces with corners/kinks/large deformation under given velocity fields,a geometrical flow with topological changes,simulations of bubble/droplet dynamics in incompressible twophase flows.In terms of accuracy it is comparable to the other existing methods.展开更多
This paper presents an approach to model the solvent-excluded surface(SES)of 3D protein molecular structures using the geometric PDE-based level-set method.The level-set method embeds the shape of 3D molecular objects...This paper presents an approach to model the solvent-excluded surface(SES)of 3D protein molecular structures using the geometric PDE-based level-set method.The level-set method embeds the shape of 3D molecular objects as an isosurface or level set corresponding to some isovalue of a scattered dense scalar field,which is saved as a discretely-sampled,rectilinear grid,i.e.,a volumetric grid.Our level-set model is described as a class of tri-cubic tensor product B-spline implicit surface with control point values that are the signed distance function.The geometric PDE is evolved in the discrete volume.The geometric PDE we use is the mean curvature specified flow,which coincides with the definition of the SES and is geometrically intrinsic.The technique of speeding up is achieved by use of the narrow band strategy incorporated with a good initial approximate construction for the SES.We get a very desirable approximate surface for the SES.展开更多
We present a general framework for a higher-order spline level-set （HLS） method and apply this to biomolecule surfaces construction. Starting from a first order energy functional, we obtain a general level set formu...We present a general framework for a higher-order spline level-set （HLS） method and apply this to biomolecule surfaces construction. Starting from a first order energy functional, we obtain a general level set formulation of geometric partial differential equation, and provide an efficient approach to solving this partial differential equation using a C2 spline basis. We also present a fast cubic spline interpolation algorithm based on convolution and the Z-transform, which exploits the local relationship of interpolatory cubic spline coefficients with respect to given function data values. One example of our HLS method is demonstrated their individual atomic coordinates which is the construction of biomolecule and solvated radii as prerequisites. surfaces （an implicit solvation interface） with展开更多
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.展开更多
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.展开更多
This study investigated the physicochemical properties,enzyme activities,volatile flavor components,microbial communities,and sensory evaluation of high-temperature Daqu(HTD)during the maturation process,and a standar...This study investigated the physicochemical properties,enzyme activities,volatile flavor components,microbial communities,and sensory evaluation of high-temperature Daqu(HTD)during the maturation process,and a standard system was established for comprehensive quality evaluation of HTD.There were obvious changes in the physicochemical properties,enzyme activities,and volatile flavor components at different storage periods,which affected the sensory evaluation of HTD to a certain extent.The results of high-throughput sequencing revealed significant microbial diversity,and showed that the bacterial community changed significantly more than did the fungal community.During the storage process,the dominant bacterial genera were Kroppenstedtia and Thermoascus.The correlation between dominant microorganisms and quality indicators highlighted their role in HTD quality.Lactococcus,Candida,Pichia,Paecilomyces,and protease activity played a crucial role in the formation of isovaleraldehyde.Acidic protease activity had the greatest impact on the microbial community.Moisture promoted isobutyric acid generation.Furthermore,the comprehensive quality evaluation standard system was established by the entropy weight method combined with multi-factor fuzzy mathematics.Consequently,this study provides innovative insights for comprehensive quality evaluation of HTD during storage and establishes a groundwork for scientific and rational storage of HTD and quality control of sauce-flavor Baijiu.展开更多
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.展开更多
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.展开更多
Soil improvement is one of the most important issues in geotechnical engineering practice.The wide application of traditional improvement techniques(cement/chemical materials)are limited due to damage ecological en-vi...Soil improvement is one of the most important issues in geotechnical engineering practice.The wide application of traditional improvement techniques(cement/chemical materials)are limited due to damage ecological en-vironment and intensify carbon emissions.However,the use of microbially induced calcium carbonate pre-cipitation(MICP)to obtain bio-cement is a novel technique with the potential to induce soil stability,providing a low-carbon,environment-friendly,and sustainable integrated solution for some geotechnical engineering pro-blems in the environment.This paper presents a comprehensive review of the latest progress in soil improvement based on the MICP strategy.It systematically summarizes and overviews the mineralization mechanism,influ-encing factors,improved methods,engineering characteristics,and current field application status of the MICP.Additionally,it also explores the limitations and correspondingly proposes prospective applications via the MICP approach for soil improvement.This review indicates that the utilization of different environmental calcium-based wastes in MICP and combination of materials and MICP are conducive to meeting engineering and market demand.Furthermore,we recommend and encourage global collaborative study and practice with a view to commercializing MICP technique in the future.The current review purports to provide insights for engineers and interdisciplinary researchers,and guidance for future engineering applications.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.12372200 and 12072242).
文摘Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.
基金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.
基金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.
文摘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 in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
基金Supported by the Natural Science Foundation of Guangxi Province(Grant Nos.2023GXNSFAA026067,2024GXN SFAA010521)the National Natural Science Foundation of China(Nos.12361079,12201149,12261026).
文摘Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金supported by the National Natural Science Foundation of China(Grant Nos.51379125,51411130131 and 11272120)the National Key Basic Research Development of China(973 Program,Grant No.2013CB036103)+1 种基金the High Te-chnology of Marine Research Project of the Ministry of Indu-stry and the Information Technology of Chinathe Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(Grant No.2013022)
文摘This paper presents a comparative study of a meshless level-set method in the simulation of sloshing flows. The numerical moving particle semi-implicit （MPS） method and a grid based schemes of the MPS and level-set methods are outlined and two violent sloshing cases are considered. The computed results are compared with the corresponding experimental data for validation. The impact pressure and the deformations of free surface induced by sloshing are comparatively analyzed, and are in good agreement with experimental ones. Results show that both the MPS and level-set methods are good tools for simulation of violent sloshing flows. However, the second pressure peaks as well as breaking and splashing of free surface by the MPS method are captured better than by the level-set method.
文摘We present and compare twoflexible and effective methodologies to predict disturbance zones ahead of underground tunnels by using elastic full-waveform inversion.One methodology uses a linearized,iterative approach based on misfit gradients computed with the adjoint method while the other uses iterative,gradient-free unscented Kalmanfiltering in conjunction with a level-set representation.Whereas the former does not involve a priori assumptions on the distribution of elastic properties ahead of the tunnel,the latter intro-duces a massive reduction in the number of explicit model parameters to be inverted for by focusing on the geometric form of potential disturbances and their average elastic properties.Both imaging methodologies are validated through successful reconstructions of simple disturbances.As an application,we consider an elastic multiple disturbance scenario.By using identical synthetic time-domain seismo-grams as test data,we obtain satisfactory,albeit different,reconstruction results from the two inversion methodologies.The computa-tional costs of both approaches are of the same order of magnitude,with the gradient-based approach showing a slight advantage.The model parameter space reduction approach compensates for this by additionally providing a posteriori estimates of model parameter uncertainty.
基金This work is partially supported by National natural science fund of China(No.91430213 and No.11571293)Hunan Provincial Innovation Foundation for Postgraduate(No.CX2015B208).
文摘A new Semi-Lagrangian scheme is proposed to discretize the surface convection-diffusion equation.The other involved equations including the the levelset convection equation,the re-initialization equation and the extension equation are also solved by S-L schemes.The S-L method removes both the CFL condition and the stiffness caused by the surface Laplacian,allowing larger time step than the Eulerian method.The method is extended to the block-structured adaptive mesh.Numerical examples are given to demonstrate the efficiency of the S-L method.
基金This work is partially supported by National natural science fund of China(No.91430213 and No.11571293)Hunan Provincial Innovation Foundation for Postgraduate(No.CX2015B208)。
文摘Semi-Lagrangian(S-L)methods have no CFL stability constraint,and are more stable than the Eulerian methods.In the literature,the S-L method for the levelset re-initialization equation was complicated,which may be unnecessary.Since the re-initialization procedure is auxiliary,we propose to use the first-order S-L scheme coupled with a projection technique to improve the accuracy at the grid points just adjacent to the interface.Standard second-order S-L method is used for evolving the level-set convection equation.The implementation is simple,including on the block-structured adaptive mesh.The efficiency of the S-L method is demonstrated by extensive numerical examples including passive convection of interfaces with corners/kinks/large deformation under given velocity fields,a geometrical flow with topological changes,simulations of bubble/droplet dynamics in incompressible twophase flows.In terms of accuracy it is comparable to the other existing methods.
基金We acknowledge the support from:NSFC grant 10701071Key Laboratory of Computational&Stochastic Mathematics and Their Applications at Universities of Hunan Province,MOE(Ministry of Education)Tier II project T207N2202IDM project NRF2007IDM-IDM002-010 and NTU project SUG 20/07.
文摘This paper presents an approach to model the solvent-excluded surface(SES)of 3D protein molecular structures using the geometric PDE-based level-set method.The level-set method embeds the shape of 3D molecular objects as an isosurface or level set corresponding to some isovalue of a scattered dense scalar field,which is saved as a discretely-sampled,rectilinear grid,i.e.,a volumetric grid.Our level-set model is described as a class of tri-cubic tensor product B-spline implicit surface with control point values that are the signed distance function.The geometric PDE is evolved in the discrete volume.The geometric PDE we use is the mean curvature specified flow,which coincides with the definition of the SES and is geometrically intrinsic.The technique of speeding up is achieved by use of the narrow band strategy incorporated with a good initial approximate construction for the SES.We get a very desirable approximate surface for the SES.
基金Bajaj is supported in part by NSF of USA under Grant No. CNS-0540033NIH under Grant Nos. P20-RR020647, R01- EB00487, R01-GM074258, R01-GM07308.+2 种基金Xu and Zhang are supported by the National Natural Science Foundation of China under Grant No. 60773165the National Basic Research 973 Program of China under Grant No. 2004CB318000. Zhang is also supported by Beijing Educational Committee Foundation under Grant No. KM200811232009.
文摘We present a general framework for a higher-order spline level-set （HLS） method and apply this to biomolecule surfaces construction. Starting from a first order energy functional, we obtain a general level set formulation of geometric partial differential equation, and provide an efficient approach to solving this partial differential equation using a C2 spline basis. We also present a fast cubic spline interpolation algorithm based on convolution and the Z-transform, which exploits the local relationship of interpolatory cubic spline coefficients with respect to given function data values. One example of our HLS method is demonstrated their individual atomic coordinates which is the construction of biomolecule and solvated radii as prerequisites. surfaces （an implicit solvation interface） with
基金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 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.
文摘This study investigated the physicochemical properties,enzyme activities,volatile flavor components,microbial communities,and sensory evaluation of high-temperature Daqu(HTD)during the maturation process,and a standard system was established for comprehensive quality evaluation of HTD.There were obvious changes in the physicochemical properties,enzyme activities,and volatile flavor components at different storage periods,which affected the sensory evaluation of HTD to a certain extent.The results of high-throughput sequencing revealed significant microbial diversity,and showed that the bacterial community changed significantly more than did the fungal community.During the storage process,the dominant bacterial genera were Kroppenstedtia and Thermoascus.The correlation between dominant microorganisms and quality indicators highlighted their role in HTD quality.Lactococcus,Candida,Pichia,Paecilomyces,and protease activity played a crucial role in the formation of isovaleraldehyde.Acidic protease activity had the greatest impact on the microbial community.Moisture promoted isobutyric acid generation.Furthermore,the comprehensive quality evaluation standard system was established by the entropy weight method combined with multi-factor fuzzy mathematics.Consequently,this study provides innovative insights for comprehensive quality evaluation of HTD during storage and establishes a groundwork for scientific and rational storage of HTD and quality control of sauce-flavor Baijiu.
基金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(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 National Natural Science Foundation of China(No.41962016)the Natural Science Foundation of NingXia(Nos.2023AAC02023,2023A1218,and 2021AAC02006).
文摘Soil improvement is one of the most important issues in geotechnical engineering practice.The wide application of traditional improvement techniques(cement/chemical materials)are limited due to damage ecological en-vironment and intensify carbon emissions.However,the use of microbially induced calcium carbonate pre-cipitation(MICP)to obtain bio-cement is a novel technique with the potential to induce soil stability,providing a low-carbon,environment-friendly,and sustainable integrated solution for some geotechnical engineering pro-blems in the environment.This paper presents a comprehensive review of the latest progress in soil improvement based on the MICP strategy.It systematically summarizes and overviews the mineralization mechanism,influ-encing factors,improved methods,engineering characteristics,and current field application status of the MICP.Additionally,it also explores the limitations and correspondingly proposes prospective applications via the MICP approach for soil improvement.This review indicates that the utilization of different environmental calcium-based wastes in MICP and combination of materials and MICP are conducive to meeting engineering and market demand.Furthermore,we recommend and encourage global collaborative study and practice with a view to commercializing MICP technique in the future.The current review purports to provide insights for engineers and interdisciplinary researchers,and guidance for future engineering applications.
