Method“Monte Carlo”in healthcare

In public health,simulation modeling stands as an invaluable asset,enabling the evaluation of new systems without their physical implementation,experimentation with existing systems without operational adjustments,and testing system limits without real-world repercussions.In simulation modeling,the Monte Carlo method emerges as a powerful yet underutilized tool.Although the Monte Carlo method has not yet gained widespread prominence in healthcare,its technological capabilities hold promise for substantial cost reduction and risk mitigation.In this review article,we aimed to explore the transformative potential of the Monte Carlo method in healthcare contexts.We underscore the significance of experiential insights derived from simulated experimentation,especially in resource-constrained scenarios where time,financial constraints,and limited resources necessitate innovative and efficient approaches.As public health faces increasing challenges,incorporating the Monte Carlo method presents an opportunity for enhanced system construction,analysis,and evaluation.

Neumann stochastic finite element method for calculating temperature field of frozen soil based on random field theory

To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory （MATLAB） sottware, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by compar- ing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatial- ly random fields.

Monte Carlo Simulation of Fractures Using Isogeometric Boundary Element Methods Based on POD-RBF

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.

PARAMETRIC STUDY OF TISSUE OPTICAL CLEARING BY LOCALIZED MECHANICAL COMPRESSION USING COMBINED FINITE ELEMENT AND MONTE CARLO SIMULATION

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.

Interval Analysis of the Finite Element Method for Stochastic Structures

A random parameter can be transformed into an interval number in the structural analysis with the concept of the confidence interval. Hence, analyses of uncertain structural systems can be used in the traditional FEM software. In some cases, the amount of solutions in stochastic structures is nearly as many as that in the traditional structural problems. In addition, a new method to evaluate the failure probability of structures is presented for the needs of the modern engineering design.

Stochastic Finite Element Method for Mechanical Vibration Based on Conjugate Gradient(CG)

When material properties, geometry parameters and applied loads are assumed to be stochastic, the vibration equation of a system is transformed to static problem by using Newmark method. In order to improve the computational efficiency and to save storage, the Conjugate Gradient （CG） method is presented. The CG is an effective method for solving a large system of linear equations and belongs to the method of iteration with rapid convergence and high precision. An example is given and calculated results are compared to validate the proposed methods.

A METHOD OF SOLVING THE FUZZY FINITE ELEMENT EQUATIONS IN MONOSOURCE FUZZY NUMBERS

For same cases the rules of monosource fuzzy numbers con be used into the solution of fuzzy stochastic finite element equations in engineering. This method can reduce the computing quantity of the solution. It can be proved that the amount of the solution is nearly as much as that with the general stochastic finite element method (SFEM). In addition, a new method to appreciate the structural fuzzy failure probability is presented for the needs of the modem engineering design.

Probabilistic stability analyses of undrained slopes by 3D random fields and finite element methods

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.

Approach of technical decision-making by element flow analysis and Monte-Carlo simulation of municipal solid waste stream

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.

Probabilistic analysis of embankment slope stability in frozen ground regions based on random finite element method

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.

Solving the Schrodinger Equation on the Basis of Finite-Difference and Monte-Carlo Approaches

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, i.e. 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.

Sequential Monte Carlo Method Toward Online RUL Assessment with Applications

Online assessment of remaining useful life(RUL) of a system or device has been widely studied for performance reliability, production safety, system conditional maintenance, and decision in remanufacturing engineering. However,there is no consistency framework to solve the RUL recursive estimation for the complex degenerate systems/device.In this paper, state space model(SSM) with Bayesian online estimation expounded from Markov chain Monte Carlo(MCMC) to Sequential Monte Carlo(SMC) algorithm is presented in order to derive the optimal Bayesian estimation.In the context of nonlinear & non-Gaussian dynamic systems, SMC(also named particle filter, PF) is quite capable of performing filtering and RUL assessment recursively. The underlying deterioration of a system/device is seen as a stochastic process with continuous, nonreversible degrading. The state of the deterioration tendency is filtered and predicted with updating observations through the SMC procedure. The corresponding remaining useful life of the system/device is estimated based on the state degradation and a predefined threshold of the failure with two-sided criterion. The paper presents an application on a milling machine for cutter tool RUL assessment by applying the above proposed methodology. The example shows the promising results and the effectiveness of SSM and SMC online assessment of RUL.

RELIABILITY OF ELASTO-PLASTIC STRUCTURE USING FINITE ELEMENT METHOD

A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases.

Identification of Water Quality Model Parameter Based on Finite Difference and Monte Carlo

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.

Analysis of the elemental spectral characteristics of single elemental capture spectrum log using Monte-Carlo simulation

In this study,the gamma-ray spectrum of single elemental capture spectrum log was simulated.By numerical simulation we obtain a single-element neutron capture gamma spectrum.The neutron and photon transportable processes were simulated using the Monte Carlo N-Particle Transport Code System(MCNP),where an Am–Be neutron source generated the neutrons and thermal neutron capture reactions with the stratigraphic elements.The characteristic gamma rays and the standard gamma spectra were recorded,from analyzing of the characteristic spectra analysis we obtain the ten elements in the stratum,such as Si,Ca,Fe,S,Ti,Al,K,Na,Cl,and Ba.Comparing with single elemental capture gamma spectrum of Schlumberger,the simulated characteristic peak and the spectral change results are in good agreement with Schlumberger.The characteristic peak positions observed also consistent with the data obtained from the National Nuclear Data Center of the International Atomic Energy Agency.The neutron gamma spectrum results calculated using this simple method have practical applications.They also serve as an reference for data processing using other types of element logging tools.

High-Order Spectral Stochastic Finite Element Analysis of Stochastic Elliptical Partial Differential Equations

This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finite element formulation of an elliptic partial differential equation having stochastic coefficients. Deriving this spectral stochastic finite element formulation couples a two-dimensional deterministic finite element formulation of an elliptic partial differential equation with generalized polynomial chaos expansions of stochastic coefficients. Further inspection of the performance of resulting spectral stochastic finite element formulation with adopting linear and quadratic (9-node or 8-node) quadrilateral elements finds that more accurate standard deviations of unknowns are surprisingly predicted using quadratic quadrilateral elements, especially under high autocorrelation function values of stochastic coefficients. In addition, creating spectral stochastic finite element results using quadratic quadrilateral elements is not unacceptably time-consuming. Therefore, this study concludes that adopting high-order elements can be a lower-cost method to improve the performance of spectral stochastic finite element method.

A Multiscale Multilevel Monte Carlo Method for Multiscale Elliptic PDEs with Random Coefficients

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.

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of this paper.Firstly,we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space,which is competitive for high-dimensional random inputs.Secondly,the a priori error estimates are derived for the state,the co-state and the control variables.Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.

Evaluation the Price of Multi-Asset Rainbow Options Using Monte Carlo Method

Solution of the system stochastic differential equations in multi dimensional case using Monte Carlo method had many useful features in compare with the other computational methods. One of them is the solution of boundary value problems to be found at just one point, if required (with associated saving in computation), whereas deterministic methods necessarily find the solution at large number of points simultaneously. This property can be particularly useful in problems such option pricing, where the value of an option is required only at the time of striking, and for the state of the market at that time. In this work we consider a European multi-asset options which mathematically described by the system of stochastic differential equations. We will apply Monte Carlo method for the solution of that system which is the price of Multi-asset rainbow options.

Stochastic Finite Element Analysis of Plate and Shell Construction

The response of random plate and shell construction is analyzed with the stochastic finite element method (SFEM). Random material properties and geometric dimensions of construction are involved in this paper. A simplified isoparametric local average model is used to describe the random field. Numerical results of the examples indicate that the approach presented herein is an economical and efficient solution for such an analysis compared with Monte Carlo simulation (MCS).