Monte Carlo simulation of network formation based on structural fragments in epoxy-anhydride systems

A method combining the Monte Carlo technique and the simple fragment approach has been developed for simulating network formation in amine-catalysed epoxy-anhydride systems. The method affords a detailed insight into the nature and composition of the network, showing the distribution of various fragments. It has been used to characterize the network formation in the reaction of the diglycidyl ester of isophthalic acid with hexahydrophthalic anhydride, catalysed by benzyldimethylamine. Pre-gel properties like number and weight distributions and average molecular weights have been calculated as a function of epoxy conversion, leading to a prediction of the gel-point conversion. Analysis of the simulated network further yields other characteristic properties such as concentration of crosslink points, distribution and concentration of elastically active chains, average molecular weight between crosslinks, sol content and mass fraction of pendent chains. A comparison has been made of the properties obtained through simulation with those predicted by the fragment approach alone, which, however, gives only average properties. The Monte Carlo simulation results clearly show that loops and other cyclic structures occur in the gel. This may account for the differences observed between the results of the simulation and the fragment model in the post-gel phase. Copyright (C) 1996 Elsevier Science Ltd.

An improved Monte Carlo scheme for simulation of synthesis of nanoparticles in reverse micelles

An improved Monte Carlo technique is presented in this work to simulate nanoparticle formation through a micellar route. The technique builds on the simulation technique proposed by Bandyopadhyaya et al. (Langmuir 2000, 16, 7139) which is general and rigorous but at the same time very computation intensive, so much so that nanoparticle formation in low occupancy systems cannot be simulated in reasonable time. In view of this, several strategies, rationalized by simple mathematical analyses, are proposed to accelerate Monte Carlo simulations. These are elimination of infructuous events, removal of excess reactant postreaction, and use of smaller micelle population a large number of times. Infructuous events include collision of an empty micelle with another empty one or with another one containing only one molecule or only a solid particle. These strategies are incorporated in a new simulation technique which divides the entire micelle population in four classes and shifts micelles from one class to other as the simulation proceeds. The simulation results, throughly tested using chi-square and other tests, show that the predictions of the improved technique remain unchanged, but with more than an order of magnitude decrease in computational effort for some of the simulations reported in the literature. A post priori validation scheme for the correctness of the simulation results has been utilized to propose a new simulation strategy to arrive at converged simulation results with near minimum computational effort.

Simulation-driven design of helmets for head impact protection using an injury-based criterion

The current paper suggests a new procedure for designing helmets for head impact protection for users such as motorcycle riders. According to the approach followed here, a helmet is mounted on a featureless Hybrid 3 headform that is used in assessing vehicles for compliance to the FMVSS 201 regulation in the USA for upper interior head impact safety. The requirement adopted in the latter standard, i.e. not exceeding a threshold HIC(d) limit of 1000, is applied in the present study as a likely criterion for adjudging the efficacy of helmets. An impact velocity of 6 m/s (13.5 mph) for the helmet-headform system striking a rigid target can probably be acceptable for ascertaining a helmet's effectiveness as a countermeasure for minimizing the risk of severe head injury. The proposed procedure is demonstrated with the help of a validated LS-DYNA model of a featureless Hybrid 3 headform in conjunction with a helmet model comprising an outer polypropylene shell to the inner surface of which is bonded a protective polyurethane foam padding of a given thickness. Based on simulation results of impact on a rigid surface, it appears that a minimum foam padding thickness of 40 mm is necessary for obtaining an acceptable value of HIC(d).

Simulation of length-preserving motions of flexible one dimensional oObjects using optimization

During the motion of one dimensional flexible objects such as ropes, chains, etc., the assumption of constant length is realistic. Moreover,their motion appears to be naturally minimizing some abstract distance measure, wherein the disturbance at one end gradually dies down along the curve defining the object. This paper presents purely kinematic strategies for deriving length-preserving transformations of flexible objects that minimize appropriate ‘motion’. The strategies involve sequential and overall optimization of the motion derived using variational calculus. Numerical simulations are performed for the motion of a planar curve and results show stable converging behavior for single-step infinitesimal and finite perturbations 1 as well as multi-step perturbations. Additionally, our generalized approach provides different intuitive motions for various problem-specific measures of motion, one of which is shown to converge to the conventional tractrix-based solution. Simulation results for arbitrary shapes and excitations are also included.

Some new results on signal reconstruction from Fourier transform magnitude spectrum

We address the problem of signal reconstruction from Fourier transform magnitude spectrum. The problem arises in many real-world scenarios where magnitude-only measurements are possible, but it is required to construct a complex-valued signal starting from those measurements. We present some new general results in this context and show that the previously known results on minimum-phase rational transfer functions, and recoverability of minimum-phase functions from magnitude spectrum, form special cases of the results reported in this paper. Some simulation results are also provided to demonstrate the practical feasibility of the reconstruction methodology.

Computational fluid dynamics simulation of open ended thermoacoustic prime mover with straight and tapered resonator

A thermoacoustic refrigerator driven by a thermoacoustic primemover is an effective way to produce durable and long lasting refrigeration due to high reliability, no exotic materials, and no moving parts. Resonator geometry is also one of the important factors that influence the performance of a thermoacoustic prime mover, namely, frequency. Computational fluid dynamics simulation of performance comparison of thermoacoustic prime mover with a straight and tapered resonator is chosen for the present study under an identical stack condition with the air as a working fluid. The frequency and pressure amplitude of oscillations obtained from simulation results were found to be more in the tapered resonator than the straight resonator. Apart from computational fluid dynamics simulation, the simulation studies have also been conducted using design environment for low-amplitude thermoacoustic energy conversion, which predicts the performance of thermoacoustic primemover comparatively well. Simulation results from computational fluid dynamics and design environment for low-amplitude thermoacoustic energy conversion were compared and found to be matching well, representing the good validity of computational fluid dynamics modeling.

Experimental and simulation studies on the performance of standing wave thermoacoustic prime mover for pulse tube refrigerator

The thermoacoustic prime mover (TAPM) has gained considerable attention as a pressure wave generator to drive pulse tube refrigerator (PTR) due to no moving parts, reasonable efficiency, use of environmental friendly working fluids etc. To drive PTCs, lower frequencies (f) with larger pressure amplitudes (Delta P) are essential, which are affected by geometric and operating parameters of TAPM as well as working fluids. For driving PTRs, a twin standing wave TAPM is built and studied by using different working fluids such as helium, argon, nitrogen and their binary mixtures. Simulation results of DeltaEc are compared with experimental data wherever possible. DeltaEc predicts slightly increased resonance frequencies, but gives larger Delta P and lower temperature difference Delta T across stack. High mass number working fluid leads to lower frequency with larger Delta P, but higher Delta T. Studies indicate that the binary gas mixture of right composition with lower Delta T can be arrived at to drive TAPM of given geometry. (C) 2013 Elsevier Ltd and IIR. All rights reserved.

A Gibbs-ensemble based technique for Monte Carlo simulation of electric double layer capacitors (EDLC) at constant voltage

Current methods for molecular simulations of Electric Double Layer Capacitors (EDLC) have both the electrodes and the electrolyte region in a single simulation box. This necessitates simulation of the electrode-electrolyte region interface. Typical capacitors have macroscopic dimensions where the fraction of the molecules at the electrode-electrolyte region interface is very low. Hence, large systems sizes are needed to minimize the electrode-electrolyte region interfacial effects. To overcome these problems, a new technique based on the Gibbs Ensemble is proposed for simulation of an EDLC. In the proposed technique, each electrode is simulated in a separate simulation box. Application of periodic boundary conditions eliminates the interfacial effects. This in addition to the use of constant voltage ensemble allows for a more convenient comparison of simulation results with experimental measurements on typical EDLCs. (C) 2014 AIP Publishing LLC.

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.

Computationally efficient models for simulation of non-ideal DC-DC converters operating in continuous and discontinuous conduction modes

This paper discusses dynamic modeling of non-isolated DC-DC converters (buck, boost and buck-boost) under continuous and discontinuous modes of operation. Three types of models are presented for each converter, namely, switching model, average model and harmonic model. These models include significant non-idealities of the converters. The switching model gives the instantaneous currents and voltages of the converter. The average model provides the ripple-free currents and voltages, averaged over a switching cycle. The harmonic model gives the peak to peak values of ripple in currents and voltages. The validity of all these models is established by comparing the simulation results with the experimental results from laboratory prototypes, at different steady state and transient conditions. Simulation based on a combination of average and harmonic models is shown to provide all relevant information as obtained from the switching model, while consuming less computation time than the latter.

Atomistic Simulation of Stacked Nucleosome Core Particles: Tail Bridging, the H4 Tail, and Effect of Hydrophobic Forces

We report the first atomistic simulation of two stacked nucleosome core particles (NCPs), with an aim to understand, in molecular detail, how they interact, the effect of salt concentration, and how different histone tails contribute to their interaction, with a special emphasis on the H4 tail, known to have the largest stabilizing effect on the NCP-NCP interaction. We do not observe specific K16-mediated interaction between the H4 tail and the H2A-H2B acidic patch, in contrast with the findings from crystallographic studies, but find that the stacking was stable even in the absence of this interaction. We perform simulations with the H4 tail (partially/completely) removed and find that the region between LYS-16 and LYS-20 of the H4 tail holds special importance in mediating the inter-NCP interaction. Performing similar tail-clipped simulations with the H3 tail removed, we compare the roles of the H3 and H4 tails in maintaining the stacking. We discuss the relevance of our simulation results to the bilayer and other liquid-crystalline phases exhibited by NCPs in vitro and, through an analysis of the histone-histone interface, identify the interactions that could possibly stabilize the inter-NCP interaction in these columnar mesophases. Through the mechanical disruption of the stacked nucleosome system using steered molecular dynamics, we quantify the strength of inter-NCP stacking in the presence and absence of salt. We disrupt the stacking at some specific sites of internucleosomal tail-DNA contact and perform a comparative quantification of the binding strengths of various tails in stabilizing the stacking. We also examine how hydrophobic interactions may contribute to the overall stability of the stacking and find a marked difference in the role of hydrophobic forces as compared with electrostatic forces in determining the stability of the stacked nucleosome system.

Efficient simulation and rendering of realistic motion of one-dimensional flexible objects

In gross motion of flexible one-dimensional (1D) objects such as cables, ropes, chains, ribbons and hair, the assumption of constant length is realistic and reasonable. The motion of the object also appears more natural if the motion or disturbance given at one end attenuates along the length of the object. In an earlier work, variational calculus was used to derive natural and length-preserving transformation of planar and spatial curves and implemented for flexible 1D objects discretized with a large number of straight segments. This paper proposes a novel idea to reduce computational effort and enable real-time and realistic simulation of the motion of flexible 1D objects. The key idea is to represent the flexible 1D object as a spline and move the underlying control polygon with much smaller number of segments. To preserve the length of the curve to within a prescribed tolerance as the control polygon is moved, the control polygon is adaptively modified by subdivision and merging. New theoretical results relating the length of the curve and the angle between the adjacent segments of the control polygon are derived for quadratic and cubic splines. Depending on the prescribed tolerance on length error, the theoretical results are used to obtain threshold angles for subdivision and merging. Simulation results for arbitrarily chosen planar and spatial curves whose one end is subjected to generic input motions are provided to illustrate the approach. (C) 2016 Elsevier Ltd. All rights reserved.

Energy conversion in shape memory alloy heat engine - Part II: Simulation

In this paper, a theoretical model proposed in Part I (Zhu et al., 2001a) is used to simulate the behavior of a twin crank NiTi SMA spring based heat engine, which has been experimentally studied by Iwanaga et al. (1988). The simulation results are compared favorably with the measurements. It is found that (1) output torque and heat efficiency decrease as rotation speed increase; (2) both output torque and output power increase with the increase of hot water temperature; (3) at high rotation speed, higher water temperature improves the heat efficiency. On the contrary, at low rotation speed, lower water temperature is more efficient; (4) the effects of initial spring length may not be monotonic as reported. According to the simulation, output torque, output power and heat efficiency increase with the decrease of spring length only in the low rotation speed case. At high rotation speed, the result might be on the contrary.

Molecular Dynamics Simulation of Peeling a DNA Molecule on Substrate

Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.

Numerical Simulation of Nox Formation in Coal Combustion with Inlet Natural Gas Burning

A full two-fluid model of reacting gas-particle flows and coal combustion is used to simulate coal combustion with and without inlet natural gas added in the inlet. The simulation results for the case without natural gas burning is in fair agreement with the experimental results reported in references. The simulation results of different natural gas adding positions indicate that the natural gas burning can form lean oxygen combustion enviroment at the combustor inlet region and the NOz concentration is reduced. The same result can be obtained from chemical equilibrium analysis.