Tissue Optical Clearing Devices(TOCDs)have been shown to increase light transmission through mechanically compressed regions of naturally turbid biological tissues.We hypothesize that zones of high compressive strain ...Tissue Optical Clearing Devices(TOCDs)have been shown to increase light transmission through mechanically compressed regions of naturally turbid biological tissues.We hypothesize that zones of high compressive strain induced by TOCD pins produce localized water displacement and reversible changes in tissue optical properties.In this paper,we demonstrate a novel combined mechanical finite element model and optical Monte Carlo model which simulates TOCD pin compression of an ex vivo porcine skin sample and modified spatial photon fluence distributions within the tissue.Results of this simulation qualitatively suggest that light transmission through the skin can be significantly affected by changes in compressed tissue geometry as well as concurrent changes in tissue optical properties.The development of a comprehensive multi-domain model of TOCD application to tissues such as skin could ultimately be used as a framework for optimizing future design of TOCDs.展开更多
A long slope consisting of spatially random soils is a common geographical feature. This paper examined the necessity of three-dimensional(3 D) analysis when dealing with slope with full randomness in soil properties....A long slope consisting of spatially random soils is a common geographical feature. This paper examined the necessity of three-dimensional(3 D) analysis when dealing with slope with full randomness in soil properties. Although 3 D random finite element analysis can well reflect the spatial variability of soil properties, it is often time-consuming for probabilistic stability analysis. For this reason, we also examined the least advantageous(or most pessimistic) cross-section of the studied slope. The concept of"most pessimistic" refers to the minimal cross-sectional average of undrained shear strength. The selection of the most pessimistic section is achievable by simulating the undrained shear strength as a 3 D random field. Random finite element analysis results suggest that two-dimensional(2 D) plane strain analysis based the most pessimistic cross-section generally provides a more conservative result than the corresponding full 3 D analysis. The level of conservativeness is around 15% on average. This result may have engineering implications for slope design where computationally tractable 2 D analyses based on the procedure proposed in this study could ensure conservative results.展开更多
This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced ...This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced by combination of Proper Orthogonal Decomposition(POD)and the Radial Basis Function(RBF).In order to avoid re-meshing and retain the geometric exactness,isogeometric boundary element method(IGABEM)is employed for simulation,in which the Non-Uniform Rational B-splines(NURBS)are employed for representing the crack surfaces and discretizing dual boundary integral equations.The stress intensity factors(SIFs)are extracted by M integral method.The numerical examples simulate several cracked structures with various uncertain parameters such as load effects,materials,geometric dimensions,and the results are verified by comparison with the analytical solutions.展开更多
This paper deals with the procedure and methodology which can be used to select the optimal treatment and disposal technology of municipal solid waste (MSW), and to provide practical and effective technical support ...This paper deals with the procedure and methodology which can be used to select the optimal treatment and disposal technology of municipal solid waste (MSW), and to provide practical and effective technical support to policy-making, on the basis of study on solid waste management status and development trend in China and abroad. Focusing on various treatment and disposal technologies and processes of MSW, this study established a Monte-Carlo mathematical model of cost minimization for MSW handling subjected to environmental constraints. A new method of element stream (such as C, H, O, N, S) analysis in combination with economic stream analysis of MSW was developed. By following the streams of different treatment processes consisting of various techniques from generation, separation, transfer, transport, treatment, recycling and disposal of the wastes, the element constitution as well as its economic distribution in terms of possibility functions was identified. Every technique step was evaluated economically. The Mont-Carlo method was then conducted for model calibration. Sensitivity analysis was also carried out to identify the most sensitive factors. Model calibration indicated that landfill with power generation of landfill gas was economically the optimal technology at the present stage under the condition of more than 58% of C, H, O, N, S going to landfill. Whether or not to generate electricity was the most sensitive factor. If landfilling cost increases, MSW separation treatment was recommended by screening first followed with incinerating partially and composting partially with residue landfilling. The possibility of incineration model selection as the optimal technology was affected by the city scale. For big cities and metropolitans with large MSW generation, possibility for constructing large-scale incineration facilities increases, whereas, for middle and small cities, the effectiveness of incinerating waste decreases.展开更多
To investigate the strong random nature of the geometric interfaces between soil and rock, a rock-soil slope is considered as a two-phase random medium. A nonlinear translation of a Gaussian field is utilized to simul...To investigate the strong random nature of the geometric interfaces between soil and rock, a rock-soil slope is considered as a two-phase random medium. A nonlinear translation of a Gaussian field is utilized to simulate the two-phase random media, such that the soil(or rock) volume fraction and the inclination of the soil layer can be examined. The finite element method with random media incorporated as the material properties is used to determine the factor of safety of the rock-soil slope. Monte-Carlo simulations are used to estimate the statistical characteristics of the factor of safety. The failure mode of the rock-soil slope is examined by observing the maximum principal plastic strain at incipient slope failure. It is found that the critical surface of a rock-soil slope is fairly irregular, and it significantly differs from that of a pure soil slope. The factor of safety is sensitive to the soil volume faction, but it is predictable. The average factor of safety could be well predicted by the weighted harmonic average between the strength of soil and rock; the prediction model is practical and simple. Parametric studies on the inclination of the soil layer demonstrate that the most instable scenario occurs when the slope angle is consistent with the inclination of the soil layer.展开更多
Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafros...Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafrost embankment design in frozen ground regions. The prediction may be further complicated due to the inherent uncertainties of material properties. Hence, stochastic analyses should be conducted. Firstly, Karhunen-Loeve expansion is applied to attain the random fields for hydraulic and thermal conductions. Next, the mixed-form modified Richards equation for mass transfer (i.e., mass equation) and the heat transport equation for heat transient flow in a variably saturated frozen soil are combined into one equation with temperature unknown. Furthermore, the finite element formulation for the coupled thermal-hydraulic fields is derived. Based on the random fields, the stochastic finite element analyses on stability of embankment are carried out. Numerical results show that stochastic analyses of embankment stability may provide a more rational picture for the distribution of factors of safety (FOS), which is definitely useful for embankment design in frozen ground regions.展开更多
The paper presents a method of numerical solution of the Schrodinger equation, which combines the finite-difference and Monte-Carlo approaches. The resulting method was effective and economical and, to a certain exten...The paper presents a method of numerical solution of the Schrodinger equation, which combines the finite-difference and Monte-Carlo approaches. The resulting method was effective and economical and, to a certain extent, not improved, <em>i</em>.<em>e</em>. optimal. The method itself is formalized as an algorithm for the numerical solution of the Schrodinger equation for a molecule with an arbitrary number of quantum particles. The method is presented and simultaneously illustrated by examples of solving the one-dimensional and multidimensional Schrodinger equation in such problems: linear one-dimensional oscillator, hydrogen atom, ion and hydrogen molecule, water, benzene and metallic hydrogen.展开更多
Identification results of water quality model parameter directly affect the accuracy of water quality numerical simulation. To overcome the difficulty of parameter identification caused by the measurement’s uncertain...Identification results of water quality model parameter directly affect the accuracy of water quality numerical simulation. To overcome the difficulty of parameter identification caused by the measurement’s uncertainty, a new method which is the coupling of Finite Difference Method and Markov Chain Monte Carlo is developed to identify the parameters of water quality model in this paper. Taking a certain long distance open channel as an example, the effects to the results of parameters identification with different noise are discussed under steady and un-steady non-uniform flow scenarios. And also this proposed method is compared with finite difference method and Nelder Mead Simplex. The results show that it can give better results by the new method. It has good noise resistance and provides a new way to identify water quality model parameters.展开更多
We propose a multiscale multilevel Monte Carlo(MsMLMC)method to solve multiscale elliptic PDEs with random coefficients in the multi-query setting.Our method consists of offline and online stages.In the offline stage,...We propose a multiscale multilevel Monte Carlo(MsMLMC)method to solve multiscale elliptic PDEs with random coefficients in the multi-query setting.Our method consists of offline and online stages.In the offline stage,we construct a small number of reduced basis functions within each coarse grid block,which can then be used to approximate the multiscale finite element basis functions.In the online stage,we can obtain the multiscale finite element basis very efficiently on a coarse grid by using the pre-computed multiscale basis.The MsMLMC method can be applied to multiscale RPDE starting with a relatively coarse grid,without requiring the coarsest grid to resolve the smallestscale of the solution.We have performed complexity analysis and shown that the MsMLMC offers considerable savings in solving multiscale elliptic PDEs with random coefficients.Moreover,we provide convergence analysis of the proposed method.Numerical results are presented to demonstrate the accuracy and efficiency of the proposed method for several multiscale stochastic problems without scale separation.展开更多
基金funded in part by NIH BRP grant#1R01HL098912-01by a fellowship for William CVogt provided by the Institute for Critical Technology and Applied Science(ICTAS)Doctoral Scholars Program at Virginia Tech.
文摘Tissue Optical Clearing Devices(TOCDs)have been shown to increase light transmission through mechanically compressed regions of naturally turbid biological tissues.We hypothesize that zones of high compressive strain induced by TOCD pins produce localized water displacement and reversible changes in tissue optical properties.In this paper,we demonstrate a novel combined mechanical finite element model and optical Monte Carlo model which simulates TOCD pin compression of an ex vivo porcine skin sample and modified spatial photon fluence distributions within the tissue.Results of this simulation qualitatively suggest that light transmission through the skin can be significantly affected by changes in compressed tissue geometry as well as concurrent changes in tissue optical properties.The development of a comprehensive multi-domain model of TOCD application to tissues such as skin could ultimately be used as a framework for optimizing future design of TOCDs.
基金supported by the Key Research&Development Plan Science and Technology Cooperation Programme of Hainan Province,China(Grant No.ZDYF2016226)the National Natural Science Foundation of China(Grant Nos.51879203,51808421)
文摘A long slope consisting of spatially random soils is a common geographical feature. This paper examined the necessity of three-dimensional(3 D) analysis when dealing with slope with full randomness in soil properties. Although 3 D random finite element analysis can well reflect the spatial variability of soil properties, it is often time-consuming for probabilistic stability analysis. For this reason, we also examined the least advantageous(or most pessimistic) cross-section of the studied slope. The concept of"most pessimistic" refers to the minimal cross-sectional average of undrained shear strength. The selection of the most pessimistic section is achievable by simulating the undrained shear strength as a 3 D random field. Random finite element analysis results suggest that two-dimensional(2 D) plane strain analysis based the most pessimistic cross-section generally provides a more conservative result than the corresponding full 3 D analysis. The level of conservativeness is around 15% on average. This result may have engineering implications for slope design where computationally tractable 2 D analyses based on the procedure proposed in this study could ensure conservative results.
基金The authors thank the financial support of National Natural Science Foundation of China(NSFC)under Grant(Nos.51904202,11902212,11901578).
文摘This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced by combination of Proper Orthogonal Decomposition(POD)and the Radial Basis Function(RBF).In order to avoid re-meshing and retain the geometric exactness,isogeometric boundary element method(IGABEM)is employed for simulation,in which the Non-Uniform Rational B-splines(NURBS)are employed for representing the crack surfaces and discretizing dual boundary integral equations.The stress intensity factors(SIFs)are extracted by M integral method.The numerical examples simulate several cracked structures with various uncertain parameters such as load effects,materials,geometric dimensions,and the results are verified by comparison with the analytical solutions.
基金Project Supported by Tsinghua Research Foundation (No. Jc2003010).
文摘This paper deals with the procedure and methodology which can be used to select the optimal treatment and disposal technology of municipal solid waste (MSW), and to provide practical and effective technical support to policy-making, on the basis of study on solid waste management status and development trend in China and abroad. Focusing on various treatment and disposal technologies and processes of MSW, this study established a Monte-Carlo mathematical model of cost minimization for MSW handling subjected to environmental constraints. A new method of element stream (such as C, H, O, N, S) analysis in combination with economic stream analysis of MSW was developed. By following the streams of different treatment processes consisting of various techniques from generation, separation, transfer, transport, treatment, recycling and disposal of the wastes, the element constitution as well as its economic distribution in terms of possibility functions was identified. Every technique step was evaluated economically. The Mont-Carlo method was then conducted for model calibration. Sensitivity analysis was also carried out to identify the most sensitive factors. Model calibration indicated that landfill with power generation of landfill gas was economically the optimal technology at the present stage under the condition of more than 58% of C, H, O, N, S going to landfill. Whether or not to generate electricity was the most sensitive factor. If landfilling cost increases, MSW separation treatment was recommended by screening first followed with incinerating partially and composting partially with residue landfilling. The possibility of incineration model selection as the optimal technology was affected by the city scale. For big cities and metropolitans with large MSW generation, possibility for constructing large-scale incineration facilities increases, whereas, for middle and small cities, the effectiveness of incinerating waste decreases.
基金supported by the International Science and Technology Cooperation Programme of Hainan Province,China (Grant No.ZDYF2016226)the National Natural Science Foundation of China(Grant No.51879203)
文摘To investigate the strong random nature of the geometric interfaces between soil and rock, a rock-soil slope is considered as a two-phase random medium. A nonlinear translation of a Gaussian field is utilized to simulate the two-phase random media, such that the soil(or rock) volume fraction and the inclination of the soil layer can be examined. The finite element method with random media incorporated as the material properties is used to determine the factor of safety of the rock-soil slope. Monte-Carlo simulations are used to estimate the statistical characteristics of the factor of safety. The failure mode of the rock-soil slope is examined by observing the maximum principal plastic strain at incipient slope failure. It is found that the critical surface of a rock-soil slope is fairly irregular, and it significantly differs from that of a pure soil slope. The factor of safety is sensitive to the soil volume faction, but it is predictable. The average factor of safety could be well predicted by the weighted harmonic average between the strength of soil and rock; the prediction model is practical and simple. Parametric studies on the inclination of the soil layer demonstrate that the most instable scenario occurs when the slope angle is consistent with the inclination of the soil layer.
基金supported by the National 973 Project of China (No. 2012CB026104)the National Natural Science Foundation of China (No. 51378057)
文摘Prediction on the coupled thermal-hydraulic fields of embankment and cutting slopes is essential to the assessment on evolution of melting zone and natural permafrost table, which is usually a key factor for permafrost embankment design in frozen ground regions. The prediction may be further complicated due to the inherent uncertainties of material properties. Hence, stochastic analyses should be conducted. Firstly, Karhunen-Loeve expansion is applied to attain the random fields for hydraulic and thermal conductions. Next, the mixed-form modified Richards equation for mass transfer (i.e., mass equation) and the heat transport equation for heat transient flow in a variably saturated frozen soil are combined into one equation with temperature unknown. Furthermore, the finite element formulation for the coupled thermal-hydraulic fields is derived. Based on the random fields, the stochastic finite element analyses on stability of embankment are carried out. Numerical results show that stochastic analyses of embankment stability may provide a more rational picture for the distribution of factors of safety (FOS), which is definitely useful for embankment design in frozen ground regions.
文摘The paper presents a method of numerical solution of the Schrodinger equation, which combines the finite-difference and Monte-Carlo approaches. The resulting method was effective and economical and, to a certain extent, not improved, <em>i</em>.<em>e</em>. optimal. The method itself is formalized as an algorithm for the numerical solution of the Schrodinger equation for a molecule with an arbitrary number of quantum particles. The method is presented and simultaneously illustrated by examples of solving the one-dimensional and multidimensional Schrodinger equation in such problems: linear one-dimensional oscillator, hydrogen atom, ion and hydrogen molecule, water, benzene and metallic hydrogen.
文摘Identification results of water quality model parameter directly affect the accuracy of water quality numerical simulation. To overcome the difficulty of parameter identification caused by the measurement’s uncertainty, a new method which is the coupling of Finite Difference Method and Markov Chain Monte Carlo is developed to identify the parameters of water quality model in this paper. Taking a certain long distance open channel as an example, the effects to the results of parameters identification with different noise are discussed under steady and un-steady non-uniform flow scenarios. And also this proposed method is compared with finite difference method and Nelder Mead Simplex. The results show that it can give better results by the new method. It has good noise resistance and provides a new way to identify water quality model parameters.
基金partially supported by the Hong Kong Ph D Fellowship Schemesupported by the Hong Kong RGC General Research Funds(Projects 27300616,17300817,and 17300318)+2 种基金National Natural Science Foundation of China(Project 11601457)Seed Funding Programme for Basic Research(HKU)Basic Research Programme(JCYJ20180307151603959)of the Science,Technology and Innovation Commission of Shenzhen Municipality。
文摘We propose a multiscale multilevel Monte Carlo(MsMLMC)method to solve multiscale elliptic PDEs with random coefficients in the multi-query setting.Our method consists of offline and online stages.In the offline stage,we construct a small number of reduced basis functions within each coarse grid block,which can then be used to approximate the multiscale finite element basis functions.In the online stage,we can obtain the multiscale finite element basis very efficiently on a coarse grid by using the pre-computed multiscale basis.The MsMLMC method can be applied to multiscale RPDE starting with a relatively coarse grid,without requiring the coarsest grid to resolve the smallestscale of the solution.We have performed complexity analysis and shown that the MsMLMC offers considerable savings in solving multiscale elliptic PDEs with random coefficients.Moreover,we provide convergence analysis of the proposed method.Numerical results are presented to demonstrate the accuracy and efficiency of the proposed method for several multiscale stochastic problems without scale separation.