Experimentally validated Raman Monte Carlo simulation for a cuboid object to obtain Raman spectroscopic signatures for hidden material

In conventional Raman spectroscopic measurements of liquids or surfaces the preferred geometry for detection of the Raman signal is the backscattering (or reflection) mode. For non-transparent layered materials, sub-surface Raman signals have been retrieved using spatially offset Raman spectroscopy (SORS), usually with light collection in the same plane as the point of excitation. However, as a result of multiple scattering in a turbid medium, Raman photons will be emitted in all directions. In this study, Monte Carlo simulations for a three-dimensional layered sample with finite geometry have been performed to confirm the detectability of Raman signals at all angles and at all sides of the object. We considered a non-transparent cuboid container (high density polyethylene) with explosive material (ammonium nitrate) inside. The simulation results were validated with experimental Raman intensities. Monte Carlo simulation results reveal that the ratio of sub-surface to surface signals improves at geometries other than backscattering. In addition, we demonstrate through simulations the effects of the absorption and scattering coefficients of the layers, and that of the diameter of the excitation beam. The advantage of collecting light from all possible 4 angles, over other collection modes, is that this technique is not geometry specific and molecular identification of layers underneath non-transparent surfaces can be obtained with minimal interference from the surface layer. To what extent all sides of the object will contribute to the total signal will depend on the absorption and scattering coefficients and the physical dimensions. Copyright (c) 2015 John Wiley & Sons, Ltd.

Role of molecular reorientation in the vitrification of molecular liquids: A Monte Carlo study of the liquid and glassy states of 2,2-dimethylbutane

The liquid and the glassy phases of 2,2-dimethylbutane have been investigated by isothermal isobaric ensemble Monte Carlo simulation. Thermodynamic Properties and radial distribution functions for both the liquid and the glass have been obtained. The radial distribution functions have been classified into three types based on the accessibility of the group. It has been shown that the structure of the Iiquid and the glass can be understood in terms of the above classification of the radial distribution functions. Molecular reorientation plays an important role in the structural rearrangement accompanying glass formation. As much as 35% of the contribution to the increase in the intermolecular interaction energy on vitrification is due to the reorientation of the neighbouring pairs of molecules. The observed changes in the dimerisation energy and the bonding energy distribution function are consistent with the observed structural changes.

Nonequilibrium-phase transition and 'specific-heat' singularity in the kinetic Ising model: A Monte Carlo study

The nonequilibrium-phase transition has been studied by Monte Carlo simulation in a ferromagnetically interacting (nearest-neighbour) kinetic Ising model in presence of a sinusoidally oscillating magnetic field. The ('specific-heat') temperature derivative of energies (averaged over a full cycle of the oscillating field) diverge near the dynamic transition point.

Prediction of the stiffness of nanoclay-polypropylene composites using a monte carlo finite element analysis approach

The present work deals with the prediction of stiffness of an Indian nanoclay-reinforced polypropylene composite (that can be termed as a nanocomposite) using a Monte Carlo finite element analysis (FEA) technique. Nanocomposite samples are at first prepared in the laboratory using a torque rheometer for achieving desirable dispersion of nanoclay during master batch preparation followed up with extrusion for the fabrication of tensile test dog-bone specimens. It has been observed through SEM (scanning electron microscopy) images of the prepared nanocomposite containing a given percentage (3–9% by weight) of the considered nanoclay that nanoclay platelets tend to remain in clusters. By ascertaining the average size of these nanoclay clusters from the images mentioned, a planar finite element model is created in which nanoclay groups and polymer matrix are modeled as separate entities assuming a given homogeneous distribution of the nanoclay clusters. Using a Monte Carlo simulation procedure, the distribution of nanoclay is varied randomly in an automated manner in a commercial FEA code, and virtual tensile tests are performed for computing the linear stiffness for each case. Values of computed stiffness modulus of highest frequency for nanocomposites with different nanoclay contents correspond well with the experimentally obtained measures of stiffness establishing the effectiveness of the present approach for further applications.

Calculation of three-phase methane-ethane binary clathrate hydrate phase equilibrium from Monte Carlo molecular simulations

Methane and ethane are the simplest hydrocarbon molecules that can form clathrate hydrates. Previous studies have reported methods for calculating the three-phase equilibrium using Monte Carlo simulation methods in systems with a single component in the gas phase. Here we extend those methods to a binary gas mixture of methane and ethane. Methane-ethane system is an interesting one in that the pure components form sII clathrate hydrate whereas a binary mixture of the two can form the sII clathrate. The phase equilibria computed from Monte Carlo simulations show a good agreement with experimental data and are also able to predict the sI-sII structural transition in the clathrate hydrate. This is attributed to the quality of the TIP4P/Ice and TRaPPE models used in the simulations. (C) 2014 Elsevier B.V. All rights reserved.

Probabilistic Seismic Design of Pile Foundations in Non Liquefiable Soil by Response Spectrum Approach

The behavior of pile foundations in non liquefiable soil under seismic loading is considerably influenced by the variability in the soil and seismic design parameters. Hence, probabilistic models for the assessment of seismic pile design are necessary. Deformation of pile foundation in non liquefiable soil is dominated by inertial force from superstructure. The present study considers a pseudo-static approach based on code specified design response spectra. The response of the pile is determined by equivalent cantilever approach. The soil medium is modeled as a one-dimensional random field along the depth. The variability associated with undrained shear strength, design response spectrum ordinate, and superstructure mass is taken into consideration. Monte Carlo simulation technique is adopted to determine the probability of failure and reliability indices based on pile failure modes, namely exceedance of lateral displacement limit and moment capacity. A reliability-based design approach for the free head pile under seismic force is suggested that enables a rational choice of pile design parameters.

Continuum percolation with discs having a distribution of radii

It is shown that for continuum percolation with overlapping discs having a distribution of radii, the net areal density of discs at percolation threshold depends non-trivially on the distribution, and is not bounded by any finite constant. Results of a Monte Carlo simulation supporting the argument are presented.

Damage assessment of structures with uncertainty by using mode-shape curvatures and fuzzy logic

Uncertainties associated with the structural model and measured vibration data may lead to unreliable damage detection. In this paper, we show that geometric and measurement uncertainty cause considerable problem in damage assessment which can be alleviated by using a fuzzy logic-based approach for damage detection. Curvature damage factor (CDF) of a tapered cantilever beam are used as damage indicators. Monte Carlo simulation (MCS) is used to study the changes in the damage indicator due to uncertainty in the geometric properties of the beam. Variation in these CDF measures due to randomness in structural parameter, further contaminated with measurement noise, are used for developing and testing a fuzzy logic system (FLS). Results show that the method correctly identifies both single and multiple damages in the structure. For example, the FLS detects damage with an average accuracy of about 95 percent in a beam having geometric uncertainty of 1 percent COV and measurement noise of 10 percent in single damage scenario. For multiple damage case, the FLS identifies damages in the beam with an average accuracy of about 94 percent in the presence of above mentioned uncertainties. The paper brings together the disparate areas of probabilistic analysis and fuzzy logic to address uncertainty in structural damage detection.

An imprecise fuzzy risk approach for water quality management of a river system

Uncertainty plays an important role in water quality management problems. The major sources of uncertainty in a water quality management problem are the random nature of hydrologic variables and imprecision (fuzziness) associated with goals of the dischargers and pollution control agencies (PCA). Many Waste Load Allocation (WLA)problems are solved by considering these two sources of uncertainty. Apart from randomness and fuzziness, missing data in the time series of a hydrologic variable may result in additional uncertainty due to partial ignorance. These uncertainties render the input parameters as imprecise parameters in water quality decision making. In this paper an Imprecise Fuzzy Waste Load Allocation Model (IFWLAM) is developed for water quality management of a river system subject to uncertainty arising from partial ignorance. In a WLA problem, both randomness and imprecision can be addressed simultaneously by fuzzy risk of low water quality. A methodology is developed for the computation of imprecise fuzzy risk of low water quality, when the parameters are characterized by uncertainty due to partial ignorance. A Monte-Carlo simulation is performed to evaluate the imprecise fuzzy risk of low water quality by considering the input variables as imprecise. Fuzzy multiobjective optimization is used to formulate the multiobjective model. The model developed is based on a fuzzy multiobjective optimization problem with max-min as the operator. This usually does not result in a unique solution but gives multiple solutions. Two optimization models are developed to capture all the decision alternatives or multiple solutions. The objective of the two optimization models is to obtain a range of fractional removal levels for the dischargers, such that the resultant fuzzy risk will be within acceptable limits. Specification of a range for fractional removal levels enhances flexibility in decision making. The methodology is demonstrated with a case study of the Tunga-Bhadra river system in India.

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation-based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1-D as well as the 3-D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1-D problems.

Optimum Design for External Seismic Stability of Geosynthetic Reinforced Soil Walls: Reliability Based Approach

In this paper, an analytical study considering the effect of uncertainties in the seismic analysis of geosynthetic-reinforced soil (GRS) walls is presented. Using limit equilibrium method and assuming sliding wedge failure mechanism, analysis is conducted to evaluate the external stability of GRS walls when subjected to earthquake loads. Target reliability based approach is used to estimate the probability of failure in three modes of failure, viz., sliding, bearing, and eccentricity failure. The properties of reinforced backfill, retained backfill, foundation soil, and geosynthetic reinforcement are treated as random variables. In addition, the uncertainties associated with horizontal seismic acceleration and surcharge load acting on the wall are considered. The optimum length of reinforcement needed to maintain the stability against three modes of failure by targeting various component and system reliability indices is obtained. Studies have also been made to study the influence of various parameters on the seismic stability in three failure modes. The results are compared with those given by first-order second moment method and Monte Carlo simulation methods. In the illustrative example, external stability of the two walls, Gould and Valencia walls, subjected to Northridge earthquake is reexamined.

Nonlinear vibration analysis of composite laminated and sandwich plates with random material properties

Nonlinear vibration analysis is performed using a C-0 assumed strain interpolated finite element plate model based on Reddy's third order theory. An earlier model is modified to include the effect of transverse shear variation along the plate thickness and Von-Karman nonlinear strain terms. Monte Carlo Simulation with Latin Hypercube Sampling technique is used to obtain the variance of linear and nonlinear natural frequencies of the plate due to randomness in its material properties. Numerical results are obtained for composite plates with different aspect ratio, stacking sequence and oscillation amplitude ratio. The numerical results are validated with the available literature. It is found that the nonlinear frequencies show increasing non-Gaussian probability density function with increasing amplitude of vibration and show dual peaks at high amplitude ratios. This chaotic nature of the dispersion of nonlinear eigenvalues is also r

Optimum Design for External Seismic Stability of Geosynthetic Reinforced Soil Walls: Reliability Based Approach

In this paper, an analytical study considering the effect of uncertainties in the seismic analysis of geosynthetic-reinforced soil (GRS) walls is presented. Using limit equilibrium method and assuming sliding wedge failure mechanism, analysis is conducted to evaluate the external stability of GRS walls when subjected to earthquake loads. Target reliability based approach is used to estimate the probability of failure in three modes of failure, viz., sliding, bearing, and eccentricity failure. The properties of reinforced backfill, retained backfill, foundation soil, and geosynthetic reinforcement are treated as random variables. In addition, the uncertainties associated with horizontal seismic acceleration and surcharge load acting on the wall are considered. The optimum length of reinforcement needed to maintain the stability against three modes of failure by targeting various component and system reliability indices is obtained. Studies have also been made to study the influence of various parameters on the seismic stability in three failure modes. The results are compared with those given by first-order second moment method and Monte Carlo simulation methods. In the illustrative example, external stability of the two walls, Gould and Valencia walls, subjected to Northridge earthquake is reexamined.

Spectral element approach to wave propagation in uncertain beam structures

This paper presents a study on the uncertainty in material parameters of wave propagation responses in metallic beam structures. Special effort is made to quantify the effect of uncertainty in the wave propagation responses at high frequencies. Both the modulus of elasticity and the density are considered uncertain. The analysis is performed using a Monte Carlo simulation (MCS) under the spectral finite element method (SEM). The randomness in the material properties is characterized by three different distributions, the normal, Weibull and extreme value distributions. Their effect on wave propagation in beams is investigated. The numerical study shows that the CPU time taken for MCS under SEM is about 48 times less than for MCS under a conventional one-dimensional finite element environment for 50 kHz loading. The numerical results presented investigate effects of material uncertainties on high frequency modes. A study is performed on the usage of different beam theories and their uncertain responses due to dynamic impulse load. These studies show that even for a small coefficient of variation, significant changes in the above parameters are noticed. A number of interesting results are presented, showing the true effects of uncertainty response due to dynamic impulse load.

Hysteresis in model spin system

The aim of this paper is to construct a nonequilibrium statistical‐mechanics theory to study hysteresis in ferromagnetic systems. We study the hysteretic response of model spin systems to periodic magnetic fields H(t) as a function of the amplitude H0 and frequency Ω. At fixed H0, we find conventional, squarelike hysteresis loops at low Ω, and rounded, roughly elliptical loops at high Ω, in agreement with experiments. For the O(N→∞), d=3, (Φ2)2 model with Langevin dynamics, we find a novel scaling behavior for the area A of the hysteresis loop, of the form A∝H0.660Ω0.33. We carry out a Monte Carlo simulation of the hysteretic response of the two‐dimensional, nearest‐neighbor, ferromagnetic Ising model. These results agree qualitatively with the results obtained for the O(N) model.