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.展开更多
Accurate prediction of coal and gas outburst(CGO)hazards is paramount in gas disaster prevention and control.This paper endeavors to overcome the constraints posed by traditional prediction indexes when dealing with C...Accurate prediction of coal and gas outburst(CGO)hazards is paramount in gas disaster prevention and control.This paper endeavors to overcome the constraints posed by traditional prediction indexes when dealing with CGO incidents under low gas pressure conditions.In pursuit of this objective,we have studied and established a mechanical model of the working face under abnormal stress and the excitation energy conditions of CGO,and proposed a method for predicting the risk of CGO under abnormal stress.On site application verification shows that when a strong outburst hazard level prediction is issued,there is a high possibility of outburst disasters occurring.In one of the three locations where we predicted strong outburst hazards,a small outburst occurred,and the accuracy of the prediction was higher than the traditional drilling cuttings index S and drilling cuttings gas desorption index q.Finally,we discuss the mechanism of CGO under the action of stress anomalies.Based on the analysis of stress distribution changes and energy accumulation characteristics of coal under abnormal stress,this article believes that the increase in outburst risk caused by high stress abnormal gradient is mainly due to two reasons:(1)The high stress abnormal gradient leads to an increase in the plastic zone of the coal seam.After the working face advances,it indirectly leads to an increase in the gas expansion energy that can be released from the coal seam before reaching a new stress equilibrium.(2)Abnormal stress leads to increased peak stress of coal body in front of working face.When coal body in elastic area transforms to plastic area,its failure speed is accelerated,which induces accelerated gas desorption and aggravates the risk of outburst.展开更多
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.展开更多
Legendre polynomial method is well-known in modeling acoustic wave characteristics.This method uses for the mechanical displacements a single polynomial expansion over the entire sandwich layers.This results in a limi...Legendre polynomial method is well-known in modeling acoustic wave characteristics.This method uses for the mechanical displacements a single polynomial expansion over the entire sandwich layers.This results in a limitation in the accuracy of the field profile restitution.Thus,it can deal with the guided waves in layered sandwich only when the material properties of adjacent layers do not change significantly.Despite the great efforts regarding this issue in the literature,there remain open questions.One of them is:“what is the exact threshold of contrasting material properties of adjacent layers for which this polynomial method cannot correctly restitute the roots of guided waves?”We investigated this numerical issue using the calculated guided phase velocities in 0°/φ/0°-carbon fibre reinforced plastics(CFRP)sandwich plates with gradually increasing angleφ.Then,we approached this numerical problem by varying the middle layer thickness h90°for the 0°/90°/0°-CFRP sandwich structure,and we proposed an exact thickness threshold of the middle layer for the Legendre polynomial method limitations.We showed that the polynomial method fails to calculate the quasi-symmetric Lamb mode in 0°/φ/0°-CFRP whenφ>25°.Moreover,we introduced a new Lamb mode so-called minimum-group-velocity that has never been addressed in literature.展开更多
This study proposes an alternative calculation mode for stresses on the slip surface(SS).The calculation of the normal stress(NS)on the SS involves examining its composition and expanding its unknown using the Taylor ...This study proposes an alternative calculation mode for stresses on the slip surface(SS).The calculation of the normal stress(NS)on the SS involves examining its composition and expanding its unknown using the Taylor series.This expansion enables the reasonable construction of a function describing the NS on the SS.Additionally,by directly incorporating the nonlinear Generalized Hoke-Brown(GHB)strength criterion and utilizing the slope factor of safety(FOS)definition,a function of the shear stress on the SS is derived.This function considers the mutual feedback mechanism between the NS and strength parameters of the SS.The stress constraints conditions are then introduced at both ends of the SS based on the spatial stress relation of one point.Determining the slope FOS and stress solution for the SS involves considering the mechanical equilibrium conditions and the stress constraint conditions satisfied by the sliding body.The proposed approach successfully simulates the tension-shear stress zone near the slope top and provides an intuitive description of the concentration effect of compression-shear stress of the SS near the slope toe.Furthermore,compared to other methods,the present method demonstrates superior processing capabilities for the embedded nonlinear GHB strength criterion.展开更多
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.展开更多
The LiMnPO4/C composite material was synthesized via a sol-gel method based on the citric acid. The X-ray diffraction （XRD）, scanning electron microscopy （SEM） and electrochemical performance tests were adopted to...The LiMnPO4/C composite material was synthesized via a sol-gel method based on the citric acid. The X-ray diffraction （XRD）, scanning electron microscopy （SEM） and electrochemical performance tests were adopted to characterize the properties of LiMnPO4/C. The XRD studies show that the pure olivine phase LiMnPO4 can be obtained at a low temperature of 500 °C. The SEM analyses illustrate that the citric acid used as the chelating reagent and carbon source can restrain the particle size of LiMnPO4/C well. The LiMnPO4/C sample synthesized at 500 °C for 10 h performs the highest initial discharge capacity of 122.6 mA-h/g, retaining 112.4 mA-h/g over 30 cycles at 0.05C rate. The citric acid based sol-gel method is favor to obtain the high electrochemical performance of LiMnPO4/C.展开更多
[Objective] This study was conducted to establish a microbial limit test method for Bupi Qiangli Ointment. [Method] The conventional method and medium dilution method were used for bacterial, mold and yeast counting i...[Objective] This study was conducted to establish a microbial limit test method for Bupi Qiangli Ointment. [Method] The conventional method and medium dilution method were used for bacterial, mold and yeast counting in sample recovery test. [Result] The medium dilution method （1：10 test solution, 0.5 ml/plate） could effectively eliminate the inhibition effect of the Bupi Qiangli Ointment, and the recovery of Staphylococcus aureus was greater than 70% in the 3 batches of samples; and the conventional method exhibited the recoveries of E. coil, Bacillus subtilis, Candida albicans and Aspergillus greater than 70% in the 3 batches of samples. [Conclusion] Due to Bupi Qiangli Ointment has strongly antibacterial effect on Staphylococcus au- reus, the medium dilution method was used for bacterial counting, and the conventional method was used for mold and yeast counting; and the conventional method was used for controlled bacterium examination of E. coll.展开更多
The degree of numerical linear independence is proposed and discussed. Based on this linear independence theory, a modified limited memory BFGS method is deve loped. Similar to the standard limited memory method, thi...The degree of numerical linear independence is proposed and discussed. Based on this linear independence theory, a modified limited memory BFGS method is deve loped. Similar to the standard limited memory method, this new method determines the new update by applying the updating formula m times to an initial positive diagonal matrix using the m previous pairs of the change in iteration and gradient. Besides the most recent pair of the change, which guarantees the quadratic termination, the choice of the other ( m -1) pairs of the change in the new method is dependent on the degree of numerical linear independence of previous search directions. In addition, the numerical linear independence theory is further discussed and the computation of the degree of linear independence is simplified. Theoretical and numerical results show that this new modified method improves efficiently the standard limited memory method.展开更多
Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured...Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.展开更多
A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio （EBR）,and the elastic modulus reduction method （EMRM） is proposed for limit load evaluation of frame structures. The EBR...A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio （EBR）,and the elastic modulus reduction method （EMRM） is proposed for limit load evaluation of frame structures. The EBR is defined employing the generalized yield criterion,and the reference EBR is determined by introducing the extrema and the degree of uniformity of EBR in the structure. The elastic modulus in the element with an EBR greater than the reference one is reduced based on the linear elastic finite element analysis and the equilibrium of strain energy. The lower-bound of limit-loads of frame structures are analyzed and the numerical example demonstrates the flexibility,accuracy and effciency of the proposed method.展开更多
Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,...Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,the total solution is the sum of the internally and externally induced parts.For the internally induced part,three numerical schemes（grid-staggering,local-nesting and piecewise continuous integration） are designed to deal with the singularity of the Green＇s function encountered in numerical calculations.For the externally induced part,by setting the velocity potential（or streamfunction） component to zero,the other component of the solution can be computed in two ways：（1） Solve for the density function from its boundary integral equation and then construct the solution from the boundary integral of the density function.（2） Use the Cauchy integral to construct the solution directly.The boundary integral can be discretized on a uniform grid along the boundary.By using local-nesting（or piecewise continuous integration）,the scheme is refined to enhance the discretization accuracy of the boundary integral around each corner point（or along the entire boundary）.When the domain is not free of data holes,the total solution contains a data-hole-induced part,and the Cauchy integral method is extended to construct the externally induced solution with irregular external and internal boundaries.An automated algorithm is designed to facilitate the integrations along the irregular external and internal boundaries.Numerical experiments are performed to evaluate the accuracy and efficiency of each scheme relative to others.展开更多
An analysis of the radar backscattering from the ocean surface covered by oil spill is presented using a mi- crowave scattering model and Monte-Carlo simulation. In the analysis, a one-dimensional rough sea sur- face ...An analysis of the radar backscattering from the ocean surface covered by oil spill is presented using a mi- crowave scattering model and Monte-Carlo simulation. In the analysis, a one-dimensional rough sea sur- face is numerically generated with an ocean waveheight spectrum for a given wind velocity. A two-layered medium is then generated by adding a thin oil layer on the simulated rough sea surface. The electric fields backscattered from the sea surface with two-layered medium are computed with the method of moments （MoM）, and the backscattering coefficients are statistically obtained with N independent samples for each oil-spilled surface using the Monte-Carlo technique for various conditions of surface roughness, oil-layer thickness, frequency, polarization and incidence angle. The numerical simulation results are compared with theoretical models for clean sea surfaces and SAR images of an off-spilled sea surface caused by the Hebei （Hebei province, China） Spirit oil tanker in 2007. Further, conditions for better oil spill extraction are sought by the numerical simulation on the effects of wind speed and oil-layer thickness at different inci- dence angles on the backscattering coefficients.展开更多
We study how to use the SR1 update to realize minimization methods for problems where the storage is critical. We give an update formula which generates matrices using information from the last m iterations. The numer...We study how to use the SR1 update to realize minimization methods for problems where the storage is critical. We give an update formula which generates matrices using information from the last m iterations. The numerical tests show that the method is efficent.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(52174162)the Fundamental Research Funds for the Central Universities(FRF-TP-20-002A3).
文摘Accurate prediction of coal and gas outburst(CGO)hazards is paramount in gas disaster prevention and control.This paper endeavors to overcome the constraints posed by traditional prediction indexes when dealing with CGO incidents under low gas pressure conditions.In pursuit of this objective,we have studied and established a mechanical model of the working face under abnormal stress and the excitation energy conditions of CGO,and proposed a method for predicting the risk of CGO under abnormal stress.On site application verification shows that when a strong outburst hazard level prediction is issued,there is a high possibility of outburst disasters occurring.In one of the three locations where we predicted strong outburst hazards,a small outburst occurred,and the accuracy of the prediction was higher than the traditional drilling cuttings index S and drilling cuttings gas desorption index q.Finally,we discuss the mechanism of CGO under the action of stress anomalies.Based on the analysis of stress distribution changes and energy accumulation characteristics of coal under abnormal stress,this article believes that the increase in outburst risk caused by high stress abnormal gradient is mainly due to two reasons:(1)The high stress abnormal gradient leads to an increase in the plastic zone of the coal seam.After the working face advances,it indirectly leads to an increase in the gas expansion energy that can be released from the coal seam before reaching a new stress equilibrium.(2)Abnormal stress leads to increased peak stress of coal body in front of working face.When coal body in elastic area transforms to plastic area,its failure speed is accelerated,which induces accelerated gas desorption and aggravates the risk of outburst.
基金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 National Natural Science Foundation of China(Grant No.12102131).
文摘Legendre polynomial method is well-known in modeling acoustic wave characteristics.This method uses for the mechanical displacements a single polynomial expansion over the entire sandwich layers.This results in a limitation in the accuracy of the field profile restitution.Thus,it can deal with the guided waves in layered sandwich only when the material properties of adjacent layers do not change significantly.Despite the great efforts regarding this issue in the literature,there remain open questions.One of them is:“what is the exact threshold of contrasting material properties of adjacent layers for which this polynomial method cannot correctly restitute the roots of guided waves?”We investigated this numerical issue using the calculated guided phase velocities in 0°/φ/0°-carbon fibre reinforced plastics(CFRP)sandwich plates with gradually increasing angleφ.Then,we approached this numerical problem by varying the middle layer thickness h90°for the 0°/90°/0°-CFRP sandwich structure,and we proposed an exact thickness threshold of the middle layer for the Legendre polynomial method limitations.We showed that the polynomial method fails to calculate the quasi-symmetric Lamb mode in 0°/φ/0°-CFRP whenφ>25°.Moreover,we introduced a new Lamb mode so-called minimum-group-velocity that has never been addressed in literature.
基金Project(52278380)supported by the National Natural Science Foundation of ChinaProject(2023JJ30670)supported by the National Science Foundation of and Technology Major Project of Hunan Province,China。
文摘This study proposes an alternative calculation mode for stresses on the slip surface(SS).The calculation of the normal stress(NS)on the SS involves examining its composition and expanding its unknown using the Taylor series.This expansion enables the reasonable construction of a function describing the NS on the SS.Additionally,by directly incorporating the nonlinear Generalized Hoke-Brown(GHB)strength criterion and utilizing the slope factor of safety(FOS)definition,a function of the shear stress on the SS is derived.This function considers the mutual feedback mechanism between the NS and strength parameters of the SS.The stress constraints conditions are then introduced at both ends of the SS based on the spatial stress relation of one point.Determining the slope FOS and stress solution for the SS involves considering the mechanical equilibrium conditions and the stress constraint conditions satisfied by the sliding body.The proposed approach successfully simulates the tension-shear stress zone near the slope top and provides an intuitive description of the concentration effect of compression-shear stress of the SS near the slope toe.Furthermore,compared to other methods,the present method demonstrates superior processing capabilities for the embedded nonlinear GHB strength criterion.
基金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.
基金Project (0991025) supported by Natural Science Foundation of Guangxi, ChinaProject (51164007) supported by the National Natural Science Foundation of ChinaProject (201101ZD008) supported by Educational Commission of Guangxi, China
文摘The LiMnPO4/C composite material was synthesized via a sol-gel method based on the citric acid. The X-ray diffraction （XRD）, scanning electron microscopy （SEM） and electrochemical performance tests were adopted to characterize the properties of LiMnPO4/C. The XRD studies show that the pure olivine phase LiMnPO4 can be obtained at a low temperature of 500 °C. The SEM analyses illustrate that the citric acid used as the chelating reagent and carbon source can restrain the particle size of LiMnPO4/C well. The LiMnPO4/C sample synthesized at 500 °C for 10 h performs the highest initial discharge capacity of 122.6 mA-h/g, retaining 112.4 mA-h/g over 30 cycles at 0.05C rate. The citric acid based sol-gel method is favor to obtain the high electrochemical performance of LiMnPO4/C.
基金Supported by Traditional Chinese Medicine and Ethnic Minority Medicine Technology Research Subject of Guizhou Administrative Bureau(QZYY2013-70,QZYY-2014-026)Joint Application of Basic Research Fund of Science and Technology Department o Yunnan Province-Guiyang University of Chinese Medicine(QKHLHZ[2014]7323)~~
文摘[Objective] This study was conducted to establish a microbial limit test method for Bupi Qiangli Ointment. [Method] The conventional method and medium dilution method were used for bacterial, mold and yeast counting in sample recovery test. [Result] The medium dilution method （1：10 test solution, 0.5 ml/plate） could effectively eliminate the inhibition effect of the Bupi Qiangli Ointment, and the recovery of Staphylococcus aureus was greater than 70% in the 3 batches of samples; and the conventional method exhibited the recoveries of E. coil, Bacillus subtilis, Candida albicans and Aspergillus greater than 70% in the 3 batches of samples. [Conclusion] Due to Bupi Qiangli Ointment has strongly antibacterial effect on Staphylococcus au- reus, the medium dilution method was used for bacterial counting, and the conventional method was used for mold and yeast counting; and the conventional method was used for controlled bacterium examination of E. coll.
文摘The degree of numerical linear independence is proposed and discussed. Based on this linear independence theory, a modified limited memory BFGS method is deve loped. Similar to the standard limited memory method, this new method determines the new update by applying the updating formula m times to an initial positive diagonal matrix using the m previous pairs of the change in iteration and gradient. Besides the most recent pair of the change, which guarantees the quadratic termination, the choice of the other ( m -1) pairs of the change in the new method is dependent on the degree of numerical linear independence of previous search directions. In addition, the numerical linear independence theory is further discussed and the computation of the degree of linear independence is simplified. Theoretical and numerical results show that this new modified method improves efficiently the standard limited memory method.
基金financially supported by the National Natural Science Foundation of China(No.41774125)Key Program of National Natural Science Foundation of China(No.41530320)+1 种基金the Key National Research Project of China(Nos.2016YFC0303100 and 2017YFC0601900)the Strategic Priority Research Program of Chinese Academy of Sciences Pilot Special(No.XDA 14020102)
文摘Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.
基金supported by the National Natural Science Foundation of China (No. 50768001)the Foundation of New Century Excellent Talents in University (No. NCET-04-0834)the Guangxi Natural Science Foundation (No. 0728026)
文摘A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio （EBR）,and the elastic modulus reduction method （EMRM） is proposed for limit load evaluation of frame structures. The EBR is defined employing the generalized yield criterion,and the reference EBR is determined by introducing the extrema and the degree of uniformity of EBR in the structure. The elastic modulus in the element with an EBR greater than the reference one is reduced based on the linear elastic finite element analysis and the equilibrium of strain energy. The lower-bound of limit-loads of frame structures are analyzed and the numerical example demonstrates the flexibility,accuracy and effciency of the proposed method.
基金supported by the Office of Naval Research (Grant No.N000141010778) to the University of Oklahomathe National Natural Sciences Foundation of China (Grant Nos. 40930950,41075043,and 4092116037) to the Institute of Atmospheric Physicsprovided by NOAA/Office of Oceanic and Atmospheric Research under NOAA-University of Oklahoma Cooperative Agreement No. (NA17RJ1227),U.S. Department of Commerce
文摘Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,the total solution is the sum of the internally and externally induced parts.For the internally induced part,three numerical schemes（grid-staggering,local-nesting and piecewise continuous integration） are designed to deal with the singularity of the Green＇s function encountered in numerical calculations.For the externally induced part,by setting the velocity potential（or streamfunction） component to zero,the other component of the solution can be computed in two ways：（1） Solve for the density function from its boundary integral equation and then construct the solution from the boundary integral of the density function.（2） Use the Cauchy integral to construct the solution directly.The boundary integral can be discretized on a uniform grid along the boundary.By using local-nesting（or piecewise continuous integration）,the scheme is refined to enhance the discretization accuracy of the boundary integral around each corner point（or along the entire boundary）.When the domain is not free of data holes,the total solution contains a data-hole-induced part,and the Cauchy integral method is extended to construct the externally induced solution with irregular external and internal boundaries.An automated algorithm is designed to facilitate the integrations along the irregular external and internal boundaries.Numerical experiments are performed to evaluate the accuracy and efficiency of each scheme relative to others.
基金The Project "Development of Korea Operational Oceanographic System (PM57041)" funded by the Ministry of Land, Transport and Maritime Affairs of Korean Governmentthe Project "Cooperation on the Development of Basic Technologies for the Yellow Sea and East China Sea Operational Oceanographic System (YOOS)" funded by CKJORC and the Basic Research Projects (PE98731, PG47770 and PE98732) of the Korea Institute Ocean Science and Technologysupport by the PASCO Corporation,Japan is also apreciated
文摘An analysis of the radar backscattering from the ocean surface covered by oil spill is presented using a mi- crowave scattering model and Monte-Carlo simulation. In the analysis, a one-dimensional rough sea sur- face is numerically generated with an ocean waveheight spectrum for a given wind velocity. A two-layered medium is then generated by adding a thin oil layer on the simulated rough sea surface. The electric fields backscattered from the sea surface with two-layered medium are computed with the method of moments （MoM）, and the backscattering coefficients are statistically obtained with N independent samples for each oil-spilled surface using the Monte-Carlo technique for various conditions of surface roughness, oil-layer thickness, frequency, polarization and incidence angle. The numerical simulation results are compared with theoretical models for clean sea surfaces and SAR images of an off-spilled sea surface caused by the Hebei （Hebei province, China） Spirit oil tanker in 2007. Further, conditions for better oil spill extraction are sought by the numerical simulation on the effects of wind speed and oil-layer thickness at different inci- dence angles on the backscattering coefficients.
文摘We study how to use the SR1 update to realize minimization methods for problems where the storage is critical. We give an update formula which generates matrices using information from the last m iterations. The numerical tests show that the method is efficent.
