We show how to obtain a fast component-by-component construction algorithm for higher order polynomial lattice rules. Such rules are useful for multivariate quadrature of high-dimensional smooth functions over the unit cube as they achieve the near optimal order of convergence. The main problem addressed in this paper is to find an efficient way of computing the worst-case error. A general algorithm is presented and explicit expressions for base~2 are given. To obtain an efficient component-by-component construction algorithm we exploit the structure of the underlying cyclic group. We compare our new higher order multivariate quadrature rules to existing quadrature rules based on higher order digital nets by computing their worst-case error. These numerical results show that the higher order polynomial lattice rules improve upon the known constructions of quasi-MonteCarlo rules based on higher order digital nets.
Much of the current technology on 3D-prestack depth migration is based on the Kirchhoff method and the traveltime information needed for this purpose is obtained by the ray-theory approximation. The main motivation driving these options is to reduce computations but it can also have an adverse effect on the quality of results. Experience with real data indicates that the corresponding poststack data migrated by the reverse-time method often yields superior results. Considerable computational efficiency can be obtained by employing nonuniform coarse grid spacing during the reverse-time migration. We followed this strategy for developing a parallel scheme for doing 3D- prestack depth migration. However, we worked only on a small subset of shot gathers whose determined by the Quasi-MonteCarlo method. The results were indistinguishable from those obtained by using a much larger population of uniformly spaced shot gathers. We present examples of ...
Objective: To estimate exposure dose of personnel around patients by Monte Carlo method and integration method. Methods: Exposures were estimated by a Monte Carlo practical program with Visual Basic 6.0 and integration method using 'pen-and-paper'. Results: Exposures for rectangle and ellipse sources were calculated. The difference between different methods for various sources were 0.88% and 0.61%, respectively. Conclusion: The results estimated by Monte Carlo method are close to those of integration method. It is illustrated that doses of other people estimated by Monte Carlo method are significant
The results of this study show that Monte Carlo simulations of diffusion in homogeneous turbulence can be formulated in terms of the Eulerian space-time velocity autocorrelation function. Numerical results obtained using this approach agree with results obtained by Taylor (1921) using the Lagrangian autocorrelation function. Use was made of the equivalence of the Lagrangian and Eulerian Monte Carlo approaches to derive analytic relations between the Lagrangian integral time scale and the Eulerian integral space and time scales. These analytic results were verified by comparison to Monte Carlo simulations and to other theoretical results. They are in general agreement with many existing theories and semi-empirical relations.
This paper reports on Monte Carlo path integral simulations of a Bianchi Type IX cosmology which yield a wave function which reflects the known chaotic dynamics of the classical system and predicts a significant probability for a correlation between large universe volume and large anisotropy.
SummaryThis paper examines the sensitivity and uncertainty of evaporation estimates from Nasser Lake, south of Egypt using data collected at the floating station Raft, 2km upstream of the Aswan High Dam, and 12 evaporation models. Based on their approach the evaporation models were ranked in five classes, respectively the Bowen ratio energy budget (BREB) approach, three energy-aerodynamic based methods, a mass transfer method, two temperature and five solar radiation-temperature approaches. State dependent parameter models were used to estimate the sensitivity indices of the different evaporation estimation methods at daily and monthly frequencies. A drastic reduction in the cost of the sensitivity analysis was obtained when coupled with efficient sampling methods, producing quasi-Monte Ca...
An integrated beam optics-nuclear processes framework is essential for accurate simulation of fragment separator beam dynamics. The code COSY INFINITY provides powerful differential algebraic methods for modeling and beam dynamics simulations in absence of beam-material interactions. However, these interactions are key for accurately simulating the dynamics of heavy ion fragmentation and fission. We have developed an extended version of the code that includes these interactions, and a set of new tools that allow efficient and accurate particle transport: by transfer map in vacuum and by Monte Carlo methods in materials. The new framework is presented, along with several examples from a preliminary layout of a fragment separator for a facility for rare isotope beams.
Organ and tissue doses and integral doses in patients due to diagnostic x-ray examinations of the chest and of the head will be presented. These data have been obtained by measuring the dose distributions in man phantoms by LiF-dosimetry. The data will be compared with data obtained by Monte-Carlo calculations. To obtain average organ doses for a special type of examination for the population or for a specified subpopulation in the variation of the exposure, data due to the different techniques used by the different radiologists has to be known. Such data will be presented. The results of the investigation will be discussed in connection with the concept of the 'somatic significant.
This lectures aim at giving graduate students an introduction to a working knowledge of path integral methods in a wide variety of fields in physics. Consequently, the the lecture notes are organized in three main parts dealing with non-relativistic quantum mechanics, many-body physics and field theory. In the first part the basic concepts of path integrals are developed in the usual heuristic, non-mathematical way followed by the standard examples of quadratic Lagrangians for which the path integrals can be solved exactly. Applications include semi-classical expansions, scattering problems and the representation of Green functions as path integrals. In the last chapter of this part it is shown how (euclidean) path integrals can be treated numerically by Monte-Carlo methods with a program for the anharmonic oscillator as an explicit example. The second part ...
The integrity and safety of beam-like structures are dependent in part on their boundary conditions which can vary with time due to damage or aging. Structural health monitoring of such structures should therefore include attention to boundary conditions. Where the boundary conditions can be represented by a lumped spring then the identification of associated stiffness parameter values may be a means to quantifying the integrity of the support. This paper investigates such a method for identifying the equivalent translational and rotational stiffness of a constrained tapered beam-like structure. An analytical model of a beam of tapered width and thickness is adopted as a simplified representation of a tower-like structure. The model is used to explore in what scenarios natural frequencies and/or nodal points might be sufficiently sensitive to changes in support conditions to be measurable indicators of damage. The method is evaluated by Monte ...
In this paper we consider the adsorption of argon on the surface of graphitized thermal carbon black and in slit pores at temperatures ranging from subcritical to supercritical conditions by the method of grand canonical Monte Carlo simulation. Attention is paid to the variation of the adsorbed density when the temperature crosses the critical point. The behavior of the adsorbed density versus pressure (bulk density) shows interesting behavior at temperatures in the vicinity of and those above the critical point and also at extremely high pressures. Isotherms at temperatures greater than the critical temperature exhibit a clear maximum, and near the critical temperature this maximum is a very sharp spike. Under the supercritical conditions and very high pressure the excess of adsorbed density decreases towards zero value for a graphite surface, while for slit pores negative excess density is possible at extremely high pressures. For imperfect pores (defined as ...
A new property which relies on the linear programming (LP) decoder, the approximate maximum-likelihood certificate (AMLC), is introduced. When using the belief propagation decoder, this property is a measure of how close the decoded codeword is to the LP solution. Using upper bounding techniques, it is demonstrated that the conditional frame error probability given that the AMLC holds is, with some degree of confidence, below a threshold. In channels with low noise, this threshold is several orders of magnitude lower than the simulated frame error rate, and our bound holds with very high degree of confidence. In contrast, showing this error performance by simulation would require very long Monte Carlo runs. When the AMLC holds, our approach thus provides the decoder with extra error detection capability, which is especially important in applications requiring high data integrity.
The effects of strong Coulomb correlations in dense three-dimensional electron-hole plasmas are studied by means of unbiased direct path integral Monte Carlo simulations. The formation and dissociation of bound states, such as excitons and bi-excitons is analyzed and the density-temperature region of their appearance is identified. At high density, the Mott transition to the fully ionized metallic state (electron-hole liquid) is detected. Particular attention is paid to the influence of the hole to electron mass ratio $M$ on the properties of the plasma. Above a critical value of about M=80 formation of a hole Coulomb crystal was recently verified [Phys. Rev. Lett. {\\bf 95}, 235006 (2005)] which is supported by additional results. Results are related to the excitonic phase diagram of intermediate valent Tm[Se,Te], where large values of $M$ have been observed experimentally.
This paper presents a trial of probabilistic simulation for performance assessment of high-level waste (HLW) disposal using the Monte Carlo method. Based on the Reference Case conceptual model in H12, a new integrated simulation system which allowed us to rapidly evaluate the effect of data uncertainty was developed. The doses to hypothetical exposure group were compared with the results of H12 that were performed by using a point-wise approach, in terms of maximum total dose. This study showed that H12 results were consistent with results of probabilistic simulation and also showed that transmissivity had a strong influence for the uncertainty of the system performance in all simulation time. (author)
Event generators that handle neutrino-nucleon interaction have been developed for the FLUKA code [1]. In earlier FLUKA versions only quasi-elastic (QEL) interactions were included, and the code relied on external event generators for the resonance (RES) and deep inelastic scattering (DIS). The new DIS+RES event generator is fully integrated in FLUKA and uses the same hadronization routines as those used for simulating hadron-nucleon interactions. Nuclear effects in neutrino-nucleus interactions are simulated within the same framework as in the FLUKA hadron-nucleus interaction model (PEANUT), thus profiting from its detailed physics modelling and longstanding benchmarking. The generators are available in the standard FLUKA distribution. They are presently under development and several improvements are planned to be implemented. The physics relevant to the neutrino-nucleon interactions and the results of comparisons with experimental data are discussed.
A clean benchmark experiment on beryllium was performed with D-T neutrons at the FNS facility of the Japan Atomic Energy Agency. The main objective was to verify the integral data related to the tritium production on lithium isotopes. Tritium production rates, as well as activation reaction rates were measured inside the beryllium assembly that was shaped as a pseudo-cylindrical slab with an area-equivalent diameter of 628 mm and a thickness of 355 mm. Experimental results were analyzed with a three-dimensional Monte Carlo transport code MCNP-4C and FENDL/MC-2.0, JENDL-3.2/3.3 neutron transport libraries. Evaluation of reaction rates was based on the cross section data taken from the JENDL Dosimetry File and ENDF B-VI data libraries. Analysis shows that all calculation combinations (transport and activation cross section libraries) used for evaluation of reaction rates give data that is agreeable with measured values within 10%.
A probability-based approach is presented as the integration of probabilistic methods and deterministic modelling based on the finite element method. An existing finite element software package was linked to an existing probabilistic package to analyse the complex mechanics that occur during the transient non-linear analysis of impact problems. This methodology is applied to a pipe whip analysis of a group-distribution-header, which results from a guillotine break, and subsequent impact with the adjacent building wall; this is a postulated accident for the Ignalina Nuclear Power Plant RBMK-1500 reactors. The uncertainties of material properties, component geometry data and loads were taken into consideration. The probabilities of failure of the impacted header and of the header support-wall were estimated given uncertainties in material properties, geometrical parameters and loading. The software ProFES was used for the probabilistic analysis and the finite element ...
The epithermal cross section shielding methods used in the lattice physics code EPRI-CELL (E-C) have been extensively studied to determine its major approximations and to examine the sensitivity of computed results to these approximations. The study has resulted in several improvements in the original methodology. These include: treatment of the external moderator source with intermediate resonance (IR) theory, development of a new Dancoff factor expression to account for clad interactions, development of a new method for treating resonance interference, and application of a generalized least squares method to compute best-estimate values for the Bell factor and group-dependent IR parameters. The modified E-C code with its new ENDF/B-V cross section library is tested for several numerical benchmark problems. Integral parameters computed by EC are compared with those obtained with point-cross section Monte Carlo calculations, and E-C fine group ...
This paper presents a new technique--Integrated Bayesian Uncertainty Estimator (IBUNE) to account for the major uncertainties of hydrologic rainfall-runoff predictions explicitly. The uncertainties from the input (forcing) data--mainly the precipitation observations and from the model parameters are reduced through a Monte Carlo Markov Chain (MCMC) scheme named Shuffled Complex Evolution Metropolis (SCEM) algorithm which has been extended to include a precipitation error model. Afterwards, the Bayesian Model Averaging (BMA) scheme is employed to further improve the prediction skill and uncertainty estimation using multiple model output. A series of case studies using three rainfall-runoff models to predict the streamflow in the Leaf River basin, Mississippi are used to examine the necessity and usefulness of this technique. The results suggests that ignoring either input forcings error or model structural uncertainty will lead to unrealistic ...
Brachytherapy is used for treating certain types of cancer by inserting radioactive sources into tumours. CDTN/CNEN is developing brachytherapy seeds to be used mainly in prostate cancer treatment. Dose calculations play a very significant role in the characterization of the developed seeds. The current state-of-the-art of computation dosimetry relies on Monte Carlo methods using, for instance, MCNP codes. However, deterministic calculations have some advantages, as, for example, short computer time to find solutions. This paper presents a software developed to calculate doses in a two-dimensional space surrounding the seed, using a deterministic algorithm. The analysed seeds consist of capsules similar to IMC6711 (OncoSeed), that are commercially available. The exposure rates and absorbed doses are computed using the Sievert integral and the Meisberger third order polynomial, respectively. The software also allows the isodose visualization at ...
In clinical applications, it is very important to accurately determine the range of a proton beam in a patient. Recently, a prototype prompt-gamma scanning system, called a prompt-gamma scanner (PGS), has been constructed at Hanyang University in Korea with the collaboration of the Korean national cancer center (NCC), based on a series of Monte Carlo simulations. The PGS system is composed of a prompt-gamma collimator, a CsI(Tl) scintillation detector, a multi-channel analyzer (MCA), a precision movement system, and an integrated readout system. The prompt-gamma collimator is designed to effectively shield high-energy neutrons from the proton beam passage and subsequent capture gammas. The PGS system can be used to scan the distribution of the prompt gammas in the patient from the proton beam passage. The distribution of the prompt gammas can then be used to determine the proton beam range or the distal dose edge in the patient. In this study, ...
In the convolution/superposition method of photon beam dose calculations, inhomogeneities are usually handled by using some form of scaling involving the relative electron densities of the inhomogeneities. In this paper the accuracy of density scaling as applied to primary electrons generated in photon interactions is examined. Monte Carlo calculations are compared with density scaling calculations for air and cork slab inhomogeneities. For individual primary photon kernels as well as for photon interactions restricted to a thin layer, the results can differ significantly, by up to 50%, between the two calculations. However, for realistic photon beams where interactions occur throughout the whole irradiated volume, the discrepancies are much less severe. The discrepancies for the kernel calculation are attributed to the scattering characteristics of the electrons and the consequent oversimplified modeling used in the density scaling method. A technique called the ...
Spatially resolved studies of star forming regions show that the assumption of spherical geometry is not realistic in most cases, with a major complication posed by the gas being ionised by multiple non-centrally located stars or star clusters. We try to isolate the effects of multiple non-centrally located stars on the temperature and ionisation structure of HII regions, via the construction of 3D photoionisation models using the 3D Monte Carlo photoionisation code MOCASSIN. We find that the true temperature fluctuations due to the stellar distribution (as opposed to the large-scale temperature gradients due to other gas properties) are small in all cases and not a significant cause of error in metallicity studies. Strong emission lines from HII regions are often used to study the metallicity of star-forming regions. We compare integrated emission line spectra from our models and quantify any systematic errors caused by the simplifying ...
Radiation damage and its recovery in the light emittance efficiency and the attenuation length of barium fluoride (BaF_2) scintillator have been investigated. The light yield and transmittance of small samples of BaF_2 scintillator were measured after #gamma#-irradiation from 0.5x10"4 to 1.1x10"5 Gy for deterioration, and after sunlight exposure for recovery. Suspension of deterioration was observed both in light yield and in attenuation length at an integrated dose of #gamma#-rays of about 10"4 Gy. Fairly good and quick recovery of the deteriorated BaF_2 scintillator was obtained by sunlight exposure. Using a Monte Carlo simulation method, the dependence of the light emittance efficiency on #gamma#-irradiation and sunlight exposure was studied. It has been found that the light emittance efficiency, as well as the attenuation length, is influenced by #gamma#-irradiation and sunlight exposure. (orig.).
The O({alpha}) electroweak radiative corrections to {gamma}{gamma}{yields}WW{yields}4f within the electroweak standard model are calculated in double-pole approximation (DPA). Virtual corrections are treated in DPA, leading to a classification into factorizable and non-factorizable contributions, and real-photonic corrections are based on complete lowest-order matrix elements for {gamma}{gamma}{yields}4f+{gamma}. Soft and collinear singularities appearing in the virtual and real corrections are combined alternatively in two different ways, namely by using the dipole subtraction method or by applying phase-space slicing. The radiative corrections are implemented in a Monte Carlo generator called Coffer {gamma}{gamma} - the computer code can be obtained from the authors upon request - which optionally includes anomalous triple and quartic gauge-boson couplings in addition and performs a convolution over realistic spectra of the photon beams. A detailed survey of ...
The O(#alpha#) electroweak radiative corrections to #gamma##gamma##->#WW#->#4f within the electroweak standard model are calculated in double-pole approximation (DPA). Virtual corrections are treated in DPA, leading to a classification into factorizable and non-factorizable contributions, and real-photonic corrections are based on complete lowest-order matrix elements for #gamma##gamma##->#4f+#gamma#. Soft and collinear singularities appearing in the virtual and real corrections are combined alternatively in two different ways, namely by using the dipole subtraction method or by applying phase-space slicing. The radiative corrections are implemented in a Monte Carlo generator called Coffer #gamma##gamma# - the computer code can be obtained from the authors upon request - which optionally includes anomalous triple and quartic gauge-boson couplings in addition and performs a convolution over realistic spectra of the photon beams. A detailed survey of ...
In most conceptual fusion power reactor designs, it is proposed to use beryllium as a neutron multiplier in the blanket. Detailed chemical composition of beryllium is necessary for evaluation of the tritium breeding ratio, and estimating the activation and transmutation of beryllium in the fusion reactor. In the present report, special attention was paid to a detailed analysis of impurities in beryllium, relevant to the tritium breeding ratio evaluation. Two different methods were used for the study of impurities: an analysis of the local sample by the ICP-MS method, and an integral analysis of the beryllium assembly, using the pulsed neutron method. The latter method was proposed as the most effective way of analyzing the integral effect to impurities in beryllium on production of the tritium on the lithium-6. The evaluation of the integral effect was based on time behaviour observations of the thermal neutron flux, ...
For heavy-ion beam driven inertial fusion ''liquid-protected'' reactor designs such as HYLIFE-II, a mixture of molten salts made of F{sup 10}, Li{sup -6}, Li{sup 7} and Be{sup 9} (called flibe) allows small chambers and final-focus magnets closer to the target with superconducting coils suffering higher radiation damage, though they can stand only a certain amount of energy deposited before quenching. This work has been primarily focusing on verifying that total energy deposited by fusion neutrons and induced gamma rays remain under such limit values and the final purpose is the optimization of the shielding of the magnetic lens system from the points of view of the geometrical configuration and of the physical nature of the materials adopted. The system is analyzed in terms of six geometrical models going from simplified up to much more realistic representations of a system of 192 beam lines, each focused by six magnets. A 3-D transport calculation ...
Volumetric properties of liquid mixtures of neon and hydrogen have been calculated using path integral hybrid Monte Carlo simulations. Realistic potentials have been used for the three interactions involved. Molar volumes and excess volumes of these mixtures have been evaluated for various compositions at 29 and 31.14 K, and 30 atm. Significant quantum effects are observed in molar volumes. Quantum simulations agree well with experimental molar volumes. Calculated excess volumes agree qualitatively with experimental values. However, contrary to the existing understanding that large positive deviations from ideal mixtures are caused due to quantum effects in Ne{endash}H{sub 2} mixtures, both classical as well as quantum simulations predict the large positive deviations from ideal mixtures. Further investigations using two other Ne{endash}H{sub 2} potentials of Lennard{endash}Jones (LJ) type show that excess volumes are very sensitive to the ...
Volumetric properties of liquid mixtures of neon and hydrogen have been calculated using path integral hybrid Monte Carlo simulations. Realistic potentials have been used for the three interactions involved. Molar volumes and excess volumes of these mixtures have been evaluated for various compositions at 29 and 31.14 K, and 30 atm. Significant quantum effects are observed in molar volumes. Quantum simulations agree well with experimental molar volumes. Calculated excess volumes agree qualitatively with experimental values. However, contrary to the existing understanding that large positive deviations from ideal mixtures are caused due to quantum effects in Ne - H_2 mixtures, both classical as well as quantum simulations predict the large positive deviations from ideal mixtures. Further investigations using two other Ne - H_2 potentials of Lennard - Jones (LJ) type show that excess volumes are very sensitive to the cross-interaction potential. ...
Recent efforts to identify methods of modeling pitting corrosion damage of high-level radioactive-waste containers are described. The need to develop models that can provide information useful to higher level system performance assessment models is emphasized, and examples of how this could be accomplished are described. Work to date has focused upon physically-based phenomenological stochastic models of pit initiation and growth. These models may provide a way to distill information from mechanistic theories in a way that provides the necessary information to the less detailed performance assessment models. Monte Carlo implementations of the stochastic theory have resulted in simulations that are, at least qualitatively, consistent with a wide variety of experimental data. The effects of environment on pitting corrosion have been included in the model using a set of simple phenomenological equations relating the parameters of the stochastic model to key ...
Details of the interaction of photons with tissue phantoms are elucidated using Monte Carlo simulations. In particular, photon sampling volumes and photon pathlengths are determined for a variety of scattering and absorption parameters. The Monte Carlo simulations are specifically designed to model light delivery and collection geometries relevant to clinical applications of optical biopsy techniques. The Monte Carlo simulations assume that light is delivered and collected by two, nearly-adjacent optical fibers and take into account the numerical aperture of the fibers as well as reflectance and refraction at interfaces between different media. To determine the validity of the Monte Carlo simulations for modeling the interactions between the photons and the tissue phantom in these geometries, the simulations were compared to measurements of aqueous suspensions of polystyrene microspheres in the ...
Monte Carlo methods have been used to compute k{sub eff} and the fundamental mode eigenfunction of critical systems since the 1950s. While such calculations have become routine using standard codes such as MCNP and SCALE/KENO, there still remain 3 concerns that must be addressed to perform calculations correctly: convergence of k{sub eff} and the fission distribution, bias in k{sub eff} and tally results, and bias in statistics on tally results. This paper provides a review of the fundamental problems inherent in Monte Carlo criticality calculations. To provide guidance to practitioners, suggested best practices for avoiding these problems are discussed and illustrated by examples.
Monte Carlo criticality calculations have been performed for over 50 years for reactor physics and criticality safety applications. With today's faster computers, these calculations are being carried out to greater precision (smaller uncertainties) in keff, and detailed distributions of power and reaction rates are being computed routinely. This paper provides a review of the fundamental theory of Monte Carlo criticality calculations, with guidance on practical methods for: (1) assuring convergence of both keff and the source distribution, (2) minimizing the bias in keff and reaction rate distributions, and (3) dealing with the under-prediction bias in uncertainties for keff and reaction rate distributions. (authors)
The mass generation in the Standard Model of Particles Physics relies on a spontaneous symmetry breaking mechanism. Its implementation is recalled, along with its constraints, both theoretical (Naturalness, Stability, Triviality, Unitarity) and experimental (limits of direct and indirect searches, prospects). Calculation techniques for observables evaluation in Perturbative Field Theory are described, particularly Helicity Amplitude method, which is given in details: fermions and vector bosons, massless and massive. Monte-Carlointegration, and structure functions approximations (which allows non-perturbative calculations) are also detailed. With these tools, a process giving to Physics beyond the Standard Model is studied: it leads to an experimental prediction for the LEP collision ring, taking the classical background into account. Technical aspects of a future photon linear collider are reviewed. The production of heavy vector bosons, ...
Single-crystalline nanostructures often exhibit gradients of surface (and/or interface) curvature that emerge from fabrication and growth processes or from thermal fluctuations. Thus, the system-inherent capillary force can initiate morphological transformations during further processing steps or during operation at elevated temperature. Therefore and because of the ongoing miniaturization of functional structures which causes a general rise in surface-to-volume ratios, solid-state capillary phenomena will become increasingly important: On the one hand diffusion-mediated capillary processes can be of practical use in view of non-conventional nanostructure fabrication methods based on self-organization mechanisms, on the other hand they can destroy the integrity of nanostructures which can go along with the failure of functionality. Additionally, capillarity-induced shape transformations are effected and can thereby be controlled by applied fields and forces (guided ...
An embodiment of a mercury-free fluorescent lamp combines a low pressure rare gas discharges with a phosphor having a quantum efficiency grater than one. The choice of the rare gas depends on a number of factors, one of which is the resonance transition energy. Less demand is placed the quantum efficiency of the phosphor for a lower energy resonance photon. Xenon has the lowest energy resonance transition of the stable rare gases at 8.5 eV (147 nm) and thus is a good candidate to study. The usefulness of a xenon-based discharge depends on the radiant emittance of the discharge at the resonance wavelength of 147 nm. The radiant emittance from a low pressure xenon positive column discharge is measured using two independent techniques. The first relies on the measurement of the resonance level density using absorption techniques. The effective decay rate of the resonance level is calculated using radiation trapping theory. The product of this density and trapped decay rate yields the ...
A high resolution Compton telescope has been identified by the Gamma Ray Astronomy Program Working Group (GRAPWG) as the highest priority major mission in gamma ray astrophysics following GLAST. This mission should provide 25-100 times improved sensitivity, relative to CGRO and INTEGRAL, for MeV gamma ray lines. It must have good performance for narrow and broad lines and for discrete and diffuse emissions. Several instrumental approaches are being pursued to achieve these goals. We discuss issues relating to this mission including alternative detector concepts, instrumental configurations, and background reduction techniques. We have pursued the development of position-sensitive solid-state detectors (Ge, Si) for a high spectral resolution Compton telescope mission. A #approx#1 m2 germanium Compton telescope of position-sensitive germanium detectors was the basis for one of the GRAPWG concepts. Preliminary Monte Carlo estimates for the ...
In this paper we investigate the effects of surface mediation on the adsorption behavior of argon at different temperatures on homogeneous graphitized thermal carbon black and on heterogeneous nongraphitized carbon black surface. The grand canonical Monte Carlo (GCMC) simulation is used to study the adsorption, and its performance is tested against a number of experimental data on graphitized thermal carbon black (which is known to be highly homogeneous) that are available in the literature. The surface-mediation effect is shown to be essential in the correct description of the adsorption isotherm because without accounting for that effect the GCMC simulation results are always greater than the experimental data in the region where the monolayer is being completed. This is due to the overestimation of the fluid-fluid interaction between particles in the first layer close to the solid surface. It is the surface mediation that reduces this fluid-fluid interaction in ...
In the fields of medical imaging, geophysical well logging, and industrial radiography, it is often of interest to characterize the spatially distributed sensitivities of neutron and gamma-ray measurement devices to the physical properties of the materials being examined. For instance, one may wish to know how the count rate in a detector varies in response to small changes in the local density of the irradiated object as a function of position. Experimental determination of such sensitivity functions is often impractical. Consequently, we have developed a general three-dimensional Monte Carlo numerical technique that allows us to directly compute the differential sensitivity of an arbitrary integral response parameter, such as a time- or energy-discriminated count rate, with respect to the spatial distribution of macroscopic cross sections and sources in the irradiated medium. Sensitivities to object density, porosity, etc., can easily be ...
We argue that the tachyons which exist in some string models are a signal of deconfinement. We compare the deconfinement temperature obtained from Monte Carlo calculations with some string models. (orig.).
We argue that the tachyons which exist in some string models are a signal of deconfinement. We compare the deconfinement temperature obtained from Monte Carlo calculations with some string models. (orig.).
Ray tracing Monte Carlo calculations have been carried out using an existing Monte Carlo radiation transport code to obtain estimates of the microsave power exiting the torus coupling links in EPT microwave manifolds. The microwave power loss and polarization at surface reflections were accounted for by treating the microwaves as plane waves reflecting off plane surfaces. Agreement on the order of 10% was obtained between the measured and calculated output power distribution for an existing EBT-S toroidal manifold. A cost effective iterative procedure utilizing the Monte Carlo history data was implemented to predict design changes which could produce increased manifold efficiency and improved output power uniformity.
ambassador in Mexico, Carlos Pascual, does he still have your confidence after the WikiLeaks scandals? (In Spanish.) SECRETARY CLINTON: Well, I deeply regret the WikiLeaks...
The Monte Carlo method is used to analyse the dose fall-off at the exit surface of a megavoltage photon beam. The convolution/superposition method of dose calculation using Monte-Carlo-generated homogeneous photon kernels is shown to be in error for exit dose calculation. Instead, photon kernels that incorporate modelling of the exit surface were generated, also using Monte Carlo, to analyse the problem, and the calculated dose fall-off using these kernels agrees well with measured data. In addition, the physics underlying the characteristics of the dose fall-off is analysed based on complete Monte Carlo modelling. Practical improvements to the convolution/superposition method are suggested. (Author).
The Monte Carlo simulation method is applied to calculate the nucleon transfer part of the imaginary optical-model potential in colliding heavy ions. In the process of transfer of nucleons from one nucleus to the other, forward and backward neutrons and protons are treated separately by taking into consideration the effect of the Coulomb forces and the potential energy surfaces (driving forces). The Pauli exclusion principle is also incorporated properly in the calculation.
A thermal neutron imaging facility for computed tomography and real-time neutron radiography is being developed at the University of Texas at Austin. The TRIGA reactor is a graphite-reflected Mark It pool-type research reactor. The neutron imaging facility will use beam port, which is at one end of a through part. Monte Carlo calculations were used to design the neutron collimator for this facility.
Monte Carlo calculations performed to validate observed capture reaction rates of "1"2"9I and "2"3"7Np were in a reasonable agreement with experimental data obtained by the Joint Institute for Nuclear Research. The most significant discrepancy between the measurements and calculations (within a factor of 3) was observed for the NpO_2 sample and Pb target at 3.7 GeV proton energy
The Fermi National Laboratory (Fermilab) operates the Tevatron proton-antiproton collider at a center-of-mass energy of {radical}s = 1.96 TeV, the is therefore the only collider which is today able to produce the heaviest known particle, the top quark. The top quark was discovered at the Tevatron by the CDF and D0 collaborations in 1995. At the Tevatron, most top quarks are produced via the strong interaction, whereby quark-antiquark annihilation dominates with 85%, and gluon fusion contributes with 15%. Considering next-to-leading order (NLO) contributions in the cross section of top-antitop quark production, leads to a slight positive asymmetry in the differential distribution of the production angle {alpha} of the top quarks. This asymmetry is due to the interference of certain NLO contributions. The charge asymmetry A in the cosine of {alpha} is predicted [14] to amount to 4-6%. Information about the partonic rest frame, necessary for a measurement of A in the observable cos ...
Monte Carlo neutron transport methods can be used to verify the applicability of point kinetics for safety analysis of nuclear reactors. KENO-NR was used to obtain the transfer function of the Advanced Neutron Source reactor and the time delay between the core power production and the external detectors, a parameter of interest to the safety systems design. The good agreement between the Monte Carlo generated transfer function and the point kinetics transfer function validates that the uncommon ANS geometry does not preclude the use of point kinetics in the frequency range that was investigated. Various features of the power spectral densities also demonstrated the applicability of point kinetics. The time delay was obtained from the cross-power spectral density (CPSD) and is {approximately}15 ms. These analyses show that frequency analysis can be used experimentally to investigate the validity of the use of point kinetics models in critical ...
During the last few decades, accuracy in photon and electron radiotherapy has increased substantially. This is partly due to enhanced linear accelerator technology, providing more flexibility in field definition (e.g. the usage of computer-controlled dynamic multileaf collimators), which led to intensity modulated radiotherapy (IMRT). Important improvements have also been made in the treatment planning process, more specifically in the dose calculations. Originally, dose calculations relied heavily on analytic, semi-analytic and empirical algorithms. The more accurate convolution/superposition codes use pre-calculated Monte Carlo dose "kernels" partly accounting for tissue density heterogeneities. It is generally recognized that the Monte Carlo method is able to increase accuracy even furt...
An accurate neutronics calculation of a local thermal environment within a fast reactor presents a major challenge. A method was previously described that used Monte Carlo techniques within a macrocell to make accurate and reasonably efficient design calculations for such an environment. This method is now being further optimized for the calculation of "2"3"8Pu production in the Fast Flux Test Facility (FFTF). Here, it is not only important to determine the "2"3"8Pu production from neutron capture in "2"3"7Np, but also to calculate the production of the contamination isotope "2"3"6Pu from high-energy (n,2n) and (#gamma#,n) reactions. The power of the Monte Carlo method to automatically include geometry and energy self-shielding is retained by optimization using fission neutron source biasing in both space and energy.
An ytterbium-169 high dose rate brachytherapy source, distinguished by an intensity-weighted average photon energy of 92.7 keV and a 32.015#+-#0.009 day half-life, is characterized in terms of the updated AAPM Task Group Report No. 43 specifications using the MCNP5 Monte Carlo computer code. In accordance with these specifications, the investigation included Monte Carlo simulations both in water and air with the in-air photon spectrum filtered to remove low-energy photons below 10 keV. TG-43 dosimetric data including S_K, D(r,#theta#), #LAMBDA#, g_L(r), F(r,#theta#), #phi#_a_n(r), and #phi#_a_n were calculated and statistical uncertainties in these parameters were derived and calculated in the appendix.
Monte Carlo calculations using the codes PENELOPE and GEANT4 have been performed to characterize the dosimetric properties of monoenergetic photon point sources in water. The dose rate in water has been calculated for energies of interest in brachytherapy, ranging between 10 keV and 2 MeV. A comparison of the results obtained using the two codes with the available data calculated with other Monte Carlo codes is carried out. A chi2-like statistical test is proposed for these comparisons. PENELOPE and GEANT4 show a reasonable agreement for all energies analyzed and distances to the source larger than 1 cm. Significant differences are found at distances from the source up to 1 cm. A similar situation occurs between PENELOPE and EGS4.
In areas where radiation dose monitoring has not been performed, it is essential to use material available in the environment be able to rapidly assess doses to individuals for immediate emergency medical care or for general estimation of the radiological consequences. It was shown that certain types of telephone cards containing microchips have the potential to be used as individual radiation dosimeters in emergency situations to detect doses over 250 mGy by luminescence measurements. In order to understand the dosimetric properties of chip cards, the components obtained from INFINIEON Company at various stages of production were used for luminescence measurements. It is found that the protecting layer used above the chips so called 'globe top' is the main source of radiation induced signal in chip cards. The globe top produced by INFINIEON at that stage is found to contain SiO2 and Epoxy. In order to improve the dosimetric properties of the chip cards, the raw material of the globe ...
A constructive method is proposed to solve a class of massive Feynman integrals within the Negative Dimensional Integration Method. The results are compared with the ones given by known positive D models and expressed in terms of dimensionally and analytically regularized integrals 3 refs.
A Monte Carlo simulation of the vacuum Bianchi type-IX (mixmaster) cosmology yields a significant correlation between large universe volume and high anisotropy. An analog of the model's chaotic classical behavior is seen in the break up of the universe wave function at large volume into fingers in the corners of the minisuperspace anisotropy potential.
The reference medium for brachytherapy dose measurements is water. Accuracy of dose measurements of brachytherapy sources is critically dependent on precise measurement of the source–detector...Full Text Available
Physicists grouped around a screen in the LEP control room at the strat-up of LEP on 14 July 1989. The emotion of the moment is clear. Carlo Rubbia, Director-General of CERN at the time, is in the centre and on his left, Herwig Schopper, former Director-General of the Organization.
The Feynmann and other fluctuation techniques have been shown to be useful for determining the multiplication of subcritical systems. The moments of the counting distribution from neutron detectors is analyzed to yield the multiplication value. We present the methodology and some selected applications and results and comparisons with Monte Carlo calculations.
Initial results obtained with a Monte Carlo computer program designed to link initial physical events in irradiated liquid water with subsequent chemical and biological events are presented. 10 refs., 4 figs., 3 tabs.
Accurate dose calculation in radiation treatment planning is most important for successful treatment. Since human body is composed of various materials and not an ideal shape, it is not easy to calculate the accurate effective dose in the patients. Many methods have been proposed to solve inhomogeneity and surface contour problems. Monte Carlo simulations are regarded as the most accurate method, but it is not appropriate for routine planning because it takes so much time. Pencil beam kernel based convolution/superposition methods were also proposed to correct those effects. Nowadays, many commercial treatment planning systems have adopted this algorithm as a dose calculation engine. The purpose of this study is to verify the accuracy of the dose calculated from pencil beam kernel based treatment planning system comparing to Monte Carlo simulations and measurements especially in inhomogeneous region. Home-made inhomogeneous phantom, Helax-TMS ...
BackgroundPlotless density estimators are those that are based on distance measures rather than counts per unit area (quadrats or plots) to estimate the density of some usually stationary...Full Text Available
Integrating phytotechnologies with energy crop production for biofuels, bioenergy, and bioproducts. In: Sixth International Phytotechnologies Conference ; ...
We review the development and application of kinetic Monte Carlo simulations to investigate defect and dopant diffusion in ion implanted silicon. In these type of Monte Carlo models, defects and dopants are treated at the atomic scale, and move according to reaction rates given as input principles. These input parameters can be obtained from first principles calculations and/or empirical molecular dynamics simulations, or can be extracted from fits to experimental data. Time and length scales differing several orders of magnitude can be followed with this method, allowing for direct comparison with experiments. The different approaches are explained and some results presented.
A method, which estimates the lightning performance of high voltage transmission lines based on the Monte-Carlo simulation technique, is described in this paper. The average number of faults which occur in a transmission line, dividing them in single-phase and three-phase faults, as well as the average grounding resistances of the transmission lines are calculated. The method is applied, on several operating Greek transmission lines, showing good correlation between predicted and field observation results. The proposed method can be used as a useful tool in the design of electric power systems, aiding in the right insulation dimensioning of a transmission line.(author)
This work reports on the use of three state-of-the-art Monte Carlo codes (MCNPX, PENELOPE, FLUKA) in the efficiency calibration of a Broad-Energy Germanium (BEGe) detector. Initial discrepancies found between the experimental and computational efficiency values are related to the poor knowledge of some physical parameters of the detector (dead-layers, crystal dimensions, etc.). As a consequence, a sensitivity analysis was carried out. Each parameter was systematically analyzed, and an accurate model of the detector was determined. The obtained results are consistent, allowing this model to be used in computational efficiency calibrations of the equipment at stake.
The Monte Carlo ray-tracing method is applied and coupled with optical properties to predict the radiation performance of solar concentrator/cavity receiver systems. Several different cavity geometries are compared on the radiation performance. A flux density distribution measurement system for dish parabolic concentrators is developed. The contours of the flux distribution for target placements at different distances from the dish vertex of a solar concentrator are taken by using an indirect method with a Lambert and a charge coupled device (CCD) camera. Further, the measured flux distributions are compared with a Monte Carlo-predicted distribution. The results can be a valuable reference for the design and assemblage of the solar collector system.
A Monte Carlo method called transient curve fitting method was developed to calculate a eigenvalues by first simulating the existing neutrons and precursors in the system, then calculating the Eigen-distribution of neutron flux and calculating the a eigenvalues using the transient results based on the Eigen-distribution by the code TMCC. The results of this method are tested by calculating Godiva Benchmark problems and they agree well with the benchmark results. Then the reasonable results of Subcritical Facility in Tsinghua University are given by TMCC, and the results are compared with the experimental results measured by Rossi-a method. Even in the deep subcritical cases, the method can give results consistent with experimental results. (authors)
The global effect of activation cross-section uncertainties on calculated radiological quantities is investigated for the first time using a methodology based on Monte Carlo random sampling. The method is applied to the calculation of the uncertainty in the contact dose rate from the gunite shielding of the National Ignition Facility chamber after 30 yr of pulsed irradiation. Some critical cross section contributing significantly to the overall uncertainty are identified. By a reasonable reduction of the uncertainty in those cross sections, the accuracy in the calculated total contact dose rate is greatly improved.
This paper presents calculations of the {sup 252}Cf-source-driven noise analysis measurements for subcritical highly enriched uranium metal cylinders using the Monte Carlo code MCNP-DSP. This code directly calculates the noise analysis data from the {sup 252}Cf- source-driven noise analysis method for both neutron and gamma ray detectors. Direct calculation of experimental observables by the Monte Carlo method allows for the benchmarking of the calculational model and the cross sections and for determining the bias in the calculation.
During the last few decades, accuracy in photon and electron radiotherapy has increased substantially. This is partly due to enhanced linear accelerator technology, providing more flexibility in field definition (e.g. the usage of computer-controlled dynamic multileaf collimators), which led to intensity modulated radiotherapy (IMRT). Important improvements have also been made in the treatment planning process, more specifically in the dose calculations. Originally, dose calculations relied heavily on analytic, semi-analytic and empirical algorithms. The more accurate convolution/superposition codes use pre-calculated Monte Carlo dose 'kernels' partly accounting for tissue density heterogeneities. It is generally recognized that the Monte Carlo method is able to increase accuracy even further. Since the second half of the 1990s, several Monte Carlo dose engines for radiotherapy treatment planning have been introduced. To ...
We perform a recursive reduction of one-loop n-point rank R tensor Feynman integrals [in short: (n,R)-integrals] for n{<=}6 with R{<=}n by representing (n,R)-integrals in terms of (n,R-1)- and (n-1,R-1)-integrals. We use the known representation of tensor integrals in terms of scalar integrals in higher dimension, which are then reduced by recurrence relations to integrals in generic dimension. With a systematic application of metric tensor representations in terms of chords, and by decomposing and recombining these representations, we find the recursive reduction for the tensors. The procedure represents a compact, sequential algorithm for numerical evaluations of tensor Feynman integrals appearing in next-to-leading order contributions to massless and massive three- and four-particle production at LHC and ...
This document presents the integrated reactors concepts by a presentation of four reactors: PIUS, SIR, IRIS and CAREM. The core conception, the operating, the safety, the economical aspects and the possible users are detailed. From the performance of the classical integrated PWR, the necessity of new innovative fuels utilization, the research of a simplified design to make easier the safety and the KWh cost decrease, a new integrated reactor is presented: SCAR 600. (A.L.B.)
Two integrated hepatitis B virus (HBV) DNA molecules were cloned from two primary hepatocellular carcinomas each containing only a single integration. One integration (C3) contained a single linear...Full Text Available
The advent of sophisticated automation equipment and computer hardware and software is changing the way manufacturing is carried out. To compete in the global marketplace, manufacturing companies must integrate these new technologies into their factories. In addition, they must integrate the planning, control, and data management methodologies needed to make effective use of these technologies. This paper provides an overview of recent approaches to achieving this enterprise integration. It then describes, using simulation as a particular example, a new tool`s perspective of enterprise integration.
A scaling Monte Carlo method has been developed to calculate diffuse reflectance from multilayered media with a wide range of optical properties in the ultraviolet-visible wavelength range. This multilayered scaling method employs the photon trajectory information generated from a single baseline Monte Carlo simulation of a homogeneous medium to scale the exit distance and exit weight of photons for a new set of optical properties in the multilayered medium. The scaling method is particularly suited to simulating diffuse reflectance spectra or creating a Monte Carlo database to extract optical properties of layered media, both of which are demonstrated in this paper. Particularly, it was found that the root-mean-square error (RMSE) between scaled diffuse reflectance, for which the anisotropy factor and refractive index in the baseline simulation were, respectively, 0.9 and 1.338, and independently simulated diffuse ...
Molecular simulation aims at simulating particles in interaction, describing a physico-chemical system. When considering Markov Chain Monte Carlo sampling in this context, we often meet the same problem of statistical efficiency as with Molecular Dynamics for the simulation of complex molecules (polymers for example). The search for a correct sampling of the space of possible configurations with respect to the Boltzmann-Gibbs distribution is directly related to the statistical efficiency of such algorithms (i.e. the ability of rapidly providing uncorrelated states covering all the configuration space). We investigated how to improve this efficiency with the help of Artificial Evolution (AE). AE algorithms form a class of stochastic optimization algorithms inspired by Darwinian evolution. Efficiency measures that can be turned into efficiency criteria have been first searched before identifying parameters that could be optimized. Relative frequencies for each type ...
International standards and guidelines for calibrating high-dose dosimetry systems to be used in industrial radiation processing recommend that dose-rate effects on dosimeters be evaluated under conditions of use. This is important when the irradiation relies on high-current electron accelerators, which usually provide very high dose-rates. However, most dosimeter calibration facilities use low-intensity gamma radiation or low-current electron accelerators, which deliver comparatively low dose-rates. Because of issues of thermal conductivity and response, portable calorimeters cannot be practically used with high-current accelerators, where product conveyor speeds under an electron beam can exceed several meters per second and the calorimeter is not suitable for use with product handling systems. As an alternative, Monte Carlo calculations can give theoretical estimates of the absorbed dose in materials with flat or complex configurations such that the results are ...
BackgroundIntegration of medical care across clinicians and settings could enhance the quality of care for patients. To date, there is limited data on the levels of integration in...Full Text Available
In this paper, the following topics are studied: uranium 238 effective integral; thermal range uranium 238 capture cross section; Americium 242 m capture cross section. The mentioned examples show that differential and integral experiments are both useful...
The discovery that hepatitis B virus (HBV) integrates into host chromosomes raises the question of whether such viral DNA integration correlates directly with the activation of specific oncogenes or...Full Text Available
The process of ensuring the quality of computer codes can be very time consuming and expensive. The Oak Ridge National Laboratory (ORNL) Monte Carlo codes all predate the existence of quality assurance (QA) standards and configuration control. The number of person-years and the amount of money spent on code development make it impossible to adhere strictly to all the current requirements. At ORNL, the Nuclear Engineering Applications Section of the Computing Applications Division is responsible for the development, maintenance, and application of the Monte Carlo codes MORSE and KENO. The KENO code is used for doing criticality analyses; the MORSE code, which has two official versions, CGA and SGC, is used for radiation transport analyses. Because KENO and MORSE were very thoroughly checked out over the many years of extensive use both in the United States and in the international community, the existing codes were open-quotes baselined.close ...
TART97 is a coupled neutron-photon, 3 dimensional, combinatorial geometry, time dependent Monte Carlo transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART97 is also incredibly fast: if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy. TART 97 is distributed on CD. This CD contains on-line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART97 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART97 and ist data files.
This paper reports on three-dimensional continuous-energy coupled neutron-gamma Monte Carlo models of the Advanced Neutron Source (ANS) final preconceptual and conceptual reference core designs developed using the Monte Carlo Neutron and Photon transport code (MCNP) Version 3b. These models contain the reactor core with control rods, the heavy water reflector tank with shutdown rods and some beam tubes, and the outer light water pool. Eighty homogenized fuel zones per fuel element are used to represent the radial and axial {sup 235}U fuel distribution. These models are the most sophisticated, physically accurate reactor physics models of the ANS currently available. The use of MCNP methods and applications to the ANS are demonstrated. Beam tube studies, coolant voiding studies, and many criticality studies have already been performed, as have studies with variance reduction techniques. In comparison with deterministic methods, MCNP proves ...
The use of computer codes based on Monte Carlo methods to perform criticality calculations has become common-place. Although results frequently published in the literature report calculated k_e_f_f values to four decimal places, people who use the codes in their everyday work say that they only believe the first two decimal places of any result. The lack of confidence in the computed k_e_f_f values may be due to the tendency of the reported standard deviation to underestimate errors associated with the Monte Carlo process. The standard deviation as reported by the codes is the standard deviation of the mean of the k_e_f_f values for individual generations in the computer simulation, not the standard deviation of the computed k_e_f_f value compared with the physical system. A more subtle problem with the standard deviation of the mean as reported by the codes is that all the k_e_f_f values from the separate generations are not statistically ...
The technical feasibility of the solar carbothermal reduction of ZnO has been successfully demonstrated in a pilot plant. The economics of this process is addressed by means of a single sensitivity analysis and a Monte-Carlo risk analysis. A medium-term and a long-term scenario have been investigated, each for a 5 and a 30 MW{sub th} plant. For a discount rate of 15% the zinc production costs vary between 482 and 245 $/t for the medium-term scenario and between 312 and 146 $/t for the long-term scenario, respectively. These costs do not account for the zinc oxide input material. In addition, a risk analysis was conducted for the 30 MW{sub th} long-term scenario. For each input parameter, a probability distribution was estimated and the probability distribution of the zinc production cost was calculated by means of a Monte-Carlo method. The expected mean zinc production costs vary from 95 $/t for a discount rate of 0%-286 $/t for a discount rate ...
A novel direct calculation method of response matrices on heterogeneous lattices by using the Monte Carlo method is proposed. These direct response matrices (DRMs) can be used in core calculations in place of the conventional homogenized lattice constants. The DRMs are formalized by four sub response matrices (sub-RMs) in order to respond to a core eigenvalue, k; thus the DRMs can be re-evaluated on each outer iteration in the core calculations. The sub-RMs can be evaluated by analyzing each neutron's trajectory from ordinary lattice calculations with the Monte Carlo code. Since these sub-RMs are calculated directly under an actual complex assembly geometry, i,e., without a homogenization process, intra-assembly heterogeneous effects can be reflected on global partial current balance calculations. With using two of the sub-RMs, which deal with neutron production probabilities for each fuel pin, and the obtained partial current balance, ...
The High Flux Isotope Reactor (HFIR) located at Oak Ridge National Laboratory is one of the world's most powerful research reactors. In 1996, one year after the demise of the Advanced Neutron Source Project, the U.S. Department of Energy embarked on an aggressive program to upgrade the neutron scattering facilities at the HFIR. These upgrades, which are now in progress, include the installation of larger beam tubes, a high-performance hydrogen cold source, and additional neutron guides and neutron scattering instruments. An extensive analysis effort was performed over the past 4 yr to support the design of the modified beamlines and new user facilities and to assess the impact of the upgrades on the integrity of the existing reactor system. The results of three of these analyses are summarized here. Specifically, results are presented for analyses related to the design of the new cold neutron source (CNS), the assessment of beam tube changes on the anticipated ...
Lentiviral vectors enter cells with high efficiency and deliver stable transduction through integration into host chromosomes, but their preference for integration within actively transcribing genes...Full Text Available
Here we describe the fabrication, optimization, and application of a microfluidic device that integrates microdialysis (MD) sampling, microchip electrophoresis (ME), and electrochemical detection...Full Text Available
We write a space-time Feynman Path Integral representation for scattered wave fields from a weakly/compact supported anisotropic non-homogeneity. (author)
PurposeDescription and analysis of the effects and side-effects of integrated mental health care in the Netherlands.Context of caseDue...Full Text Available
We demonstrate the use of a variational method to determine a quantitative lower bound on the rate of convergence of Markov Chain Monte Carlo (MCMC) algorithms as a function of the target density and proposal density. The bound relies on approximating the second largest eigenvalue in the spectrum of the MCMC operator using a variational principle and the approach is applicable to problems with continuous state spaces. We apply the method to one dimensional examples with Gaussian and quartic target densities, and we contrast the performance of the basic Metropolis-Hastings algorithms with a ``smart'' variant that incorporates gradient information into the trial moves. We find that the variational method agrees quite closely with numerical simulations. We also see that the smart MCMC algorithm often fails to converge geometrically in the tails of the target density except in the simplest case we examine, and even then care must be taken to choose the appropriate ...
This book examines the applications of Monte Carlo (MC) calculations in therapeutic nuclear medicine, from basic principles to computer implementations of software packages and their applications in radiation dosimetry and treatment planning. It is written for nuclear medicine physicists and physicians as well as radiation oncologists, and can serve as a supplementary text for medical imaging, radiation dosimetry and nuclear engineering graduate courses in science, medical and engineering faculties. With chapters is written by recognised authorities in that particular field, the book covers the entire range of MC applications in therapeutic medical and health physics, from its use in imaging prior to therapy to dose distribution modelling targeted radiotherapy. The contributions discuss the fundamental concepts of radiation dosimetry, radiobiological aspects of targeted radionuclide therapy and the various components and steps required for implementing a dose ...
This study presents results from Monte Carlo radiation transport calculations aimed at characterizing a novel methodology being developed to detect partial defects in Pressurized Water Reactor (PWR) spent fuel assemblies (SFAs). The methodology uses a combination of measured neutron and gamma fields inside a spent fuel assembly in an in-situ condition where no movement of the fuel assembly is required. Previous studies performed on single isolated assemblies resulted in a unique base signature that would change when some of the fuel in the assembly is replaced with dummy fuel. These studies indicate that this signature is still valid in the in-situ condition enhancing the prospect of building a practical tool, Partial Defect Detector (PDET), which can be used in the field for partial defect detection.
To represent a heterogeneous unsaturated fractured rock by its homogeneous equivalent, Monte Carlo simulations are used to obtain upscaled (effective) flow properties. In this study, we present a numerical procedure for upscaling the van Genuchten parameters of unsaturated fractured rocks by conducting Monte Carlo simulations of the unsaturated flow in a domain under gravity-dominated regime. The simulation domain can be chosen as the scale of block size in the field-scale modeling. The effective conductivity is computed from the steady-state flux at the lower boundary and plotted as a function of the averaging pressure head or saturation over the domain. The scatter plot is then fitted using van Genuchten model and three parameters, i.e., the saturated conductivity K{sub s}, the air-entry parameter {alpha}, the pore-size distribution parameter n, corresponding to this model are considered as the effective K{sub s}, effective {alpha}, and ...
An algorithm for creating synthetic telescope images of Smoothed Particle Hydrodynamics (SPH) density fields is presented, which utilises the adaptive nature of the SPH formalism in full. The imaging process uses Monte Carlo Radiative Transfer (MCRT) methods to model the scattering and absorption of photon packets in the density field, which then exit the system and are captured on a pixelated image plane, creating a 2D image (or a 3D datacube, if the photons are also binned by their wavelength). The algorithm is implemented on the density field directly: no gridding of the field is required, allowing the density field to be described to an identical level of accuracy as the simulations that generated it. Some applications of the method to star and planet formation simulations are presented to illustrate the advantages of this new technique, and suggestions as to how this framework could support a Radiative Equilibrium algorithm are also given as an indication for ...
A three-dimensional, continuous-energy, coupled neutron-gamma Monte Carlo model of the advanced neutron source (ANS) final preconceptual reference core design has been developed using MCNP Version 3b. This model contains the reactor core with control rods, the heavy-water reflector tank with shutdown rods and representative beam tubes, and the outer light-water poor. Eighty homogenized fuel zones per fuel element are used to represent the radical and axial {sup 235}U fuel grading. This model is the most sophisticated, physically accurate reactor physics model of the ANS currently available. The purpose of this summary is to demonstrate the MCNP methods and applications for the ANS.
The MCNPE-BO and MCNP4 Monte Carlo electron-photon codes were used to calculate the dose equivalent per unit fluence at various depths in tissue-equivalent slab phantoms for broad parallel beams of monoenergetic electrons with energies from 50 keV to 10 MeV. The study was carried out in the framework of the activities of a ICRP/ICRU Joint Task Group with the support of EURADOS WG4 (Numerical Dosimetry). Some preliminary results and comparisons as well as a general discussion on the performances of the codes are presented, demonstrating quite a satisfactory agreement among the results obtained using the two codes and those of other authors. (author).
The liquid/vapor coexistence density, the partial vapor pressure, and the heat of vaporization were calculated using Gibbs ensemble Monte Carlo simulation techniques. Long-range interactions such as charge-charge, charge-dipole, and dipole-dipole were evaluated using Ewald summation techniques. A polarizable potential model was used to describe the water-water interactions (Dang and Chang, J. Chem. Phys. 106, 8149, 1997). The model yields good agreement with the corresponding experimental data in the lower temperature region and moderate agreement in the higher temperature region. The critical temperature and density were estimated to be 565 K and 0.28 g/cm3.
The aim of this work was to compare dose calculation algorithm results at orthovoltage energies for a phantom composed of a bone slab in water. The calculation methods investigated were: no correction, ETAR, Batho, convolution/superposition and Monte Carlo. All algorithms calculated depth dose curves in a water phantom within 4% of experiment. However in the bone phantom, differences of over 40% between the No Correction / ETAR / Batho / Convolution and Monte Carlo results in the 1 cm thick bone slab were observed. These differences are predominantly because the algorithms do not account for the differing atomic number of the bone compared to water The increased dose to bone and the tissue adjacent to the bone interface should be considered when treating with orthovoltage photons. Copyright (1998) Australasian Physical and Engineering Sciences in Medicine
The Monte-Carlo method and experimental methods were used to determine the neutron fluxes in the irradiation channels of the Ghana Research Reactor -1. The MCNP5 code was used for this purpose to simulate the radial and axial distribution of the neutron fluxes within all the ten irradiation channels. The results obtained were compared with the experimental results. After the MCNP simulation and experimental procedure, it was observed that axially, the fluxes rise to a peak before falling and then finally leveling out. Axially and radially, it was also observed that the fluxes in the centre of the channels were lower than on the sides. Radially, the fluxes dip in the centre while it increases steadily towards the sides of the channels. The results have shown that there are flux variations within the irradiation channels both axially and radially. (au)
HANARO fuel test loop will be equipped in HANARO to obtain the development betterment of advanced fuel and materials through the irradiation test. The object of this study is to evaluate the structural integrity of FTL in-pool piping by investigating a dynamic analysis of the loop containing a postulated rupture section. The method to perform the dynamic analysis and structural integrity evaluation caused by the pipe whip in water environment can be a reference for a similar structural integrity evaluation. (author). 7 refs., 39 tabs., 34 figs.
A method for ranking interval objects is proposed and analyzed; the characteristics of those objects are represented by pessimistic, optimistic, and most probable estimates. The method is based on the approximation of the binary probability preference relation by the binary median preference relation. The method is verified using statistical modeling (the Monte Carlo method). The proposed approach can be used for ranking nonreusable and reusable objects.
The electrogeometric model cannot consider the influence of some parameters on the screening of the wire conductor of the H.V. lines. For analysing the influence of these parameters (the conductor potential, the influence of the distance from phase conductor to earth wire, the influence of the reducing factor of the striking earth distance, the earth slope, etc.), the Monte Carlo method is preferable. The influence of the earth slope on the screening effect of the wire conductor for different conditions is presented in this paper. (author)
Results of the evaluation of the scattered radiation effect on the image of defects in examined products are given. The formula is suggested which characterizes the image of defects on a scintillation screen as well as accumulation factors are calculated by the Monte-Carlo method. Values of fluctuations of absorbed energy in the screen used for the evaluation of sensitivity are experimentally obtained. The effect of high-energy bremsstrahlung on the defect detection is analyzed.
Previous models and evaluation techniques based on the Monte Carlo simulation for the reliability assessment of mixed hydrothermal generation systems are extended. Models for reservoirs on multiple rivers and pumped storage plants are developed. Operating and water management policies are extended to include the relative economical effects of operating alternative sources of generation. These effects and the benefits are described and discussed using an extended IEEE reliability test system. (Author).
A microscopic calculation of the level density is performed, based on a combinatorial evaluation using a realistic single-particle level scheme. This calculation relies on a fast Monte Carlo algorithm, allowing to consider heavy nuclei (i.e., large shell model spaces) which could not be treated previously in combinatorial approaches. An exhaustive comparison of the predicted neutron s-wave resonance spacings with experimental data for a wide range of nuclei is presented.
A Monte Carlo simulation of the vacuum Bianchi type-IX (mixmaster) cosmology yields a significant correlation between large universe volume and high anisotropy. An analog of the model's chaotic classical behavior is seen in the break up of the universe wave function at large volume into fingers in the corners of the minisuperspace anisotropy potential.
We study by means of Quantum Monte Carlo simulations based on the Worm Algorithm the low temperature (down to T = 0.05 K) properties of parahydrogen clusters comprising up to 40 molecules. Three different intermolecular interactions are employed: the Silvera-Goldman, the Buck and the Lennard-Jones potential. Despite important discrepancies observed in the numerical estimates of energy and superfluid fraction, the mechanism by which clusters melt at low T is independent of the particular choice of the potential, whose only effect is to alter the temperature scale.
The Milstar communications satellite system will provide secure antijam communication capabilities for DOD operations into the next century. In order to accomplish this task, the Milstar system will employ precise timekeeping on its satellites and at its ground control stations. The constellation will consist of four satellites in geosynchronous orbit, each carrying a set of four rubidium (Rb) atomic clocks. Several times a day, during normal operation, the Mission Control Element (MCE) will collect timing information from the constellation, and after several days use this information to update the time and frequency of the satellite clocks. The MCE will maintain precise time with a cesium (Cs) atomic clock, synchronized to UTC(USNO) via a GPS receiver. We have developed a Monte Carlo simulation of Milstar's space segment timekeeping. The simulation includes the effects of: uplink/downlink time transfer noise; satellite crosslink time transfer noise; satellite ...
Monte Carlo simulation of beam modifiers such as physical wedges and compensating filters has been performed with a rectilinear voxel geometry module. A modified version of the EGS4/DOSXYZ code has been developed for this purpose. The new implementations have been validated against the BEAM Monte Carlo code using its standard component modules (CMs) in several geometrical conditions. No significant disagreements were found within the statistical errors of 0.5% for photons and 2% for electrons. The clinical applicability and flexibility of the new version of the code has been assessed through an extensive verification versus dosimetric data. Both Varian multi-leaf collimator (MLC) wedges and standard wedges have been simulated and compared against experiments for 6 MV photon beams and different field sizes. Good agreement was found between calculated and measured depth doses and lateral dose profiles along both wedged and unwedged directions for ...
The status of Monte Carlo system for the simulation of Bremsstrahlung in arbitrary decays and for the decay itself of #tau# -lepton is reviewed. During the last #tau# -lepton conference in 2010 several developments of the last two years have been presented: (i) For the TAUOLA Monte Carlo generator of #tau# -lepton decays, automated and simultaneous use of many versions of form-factors for the calculation of optional weights for fits was developed and checked to work in Belle and BaBar software environment. On-going work on alternative parameterizations of hadronic decays is presented too. (ii) The TAUOLA universal interface based on HepMC (the C++ event record) is now public. A similar interface for PHOTOS is now also public. (iii) Extension of PHOTOS Monte Carlo for QED Bremsstrahlung in decays featuring kernels based on complete first order matrix element are gradually becoming widely available thanks to properties of the ...
The purpose of this study was to investigate if interfraction and intrafraction motion in free-breathing and gated lung IMRT can lead to systematic dose differences between 3DCT and 4DCT. Dosimetric...Full Text Available
The purpose of this work was to study and quantify the differences in dose distributions computed with some of the newest dose calculation algorithms available in commercial planning systems. The study was done for clinical cases originally calculated with pencil beam convolution (PBC) where large density inhomogeneities were present. Three other dose algorithms were used: a pencil beam like algorithm, the anisotropic analytic algorithm (AAA), a convolution superposition algorithm, collapsed cone convolution (CCC), and a Monte Carlo program, voxel Monte Carlo (VMC++). The dose calculation algorithms were compared under static field irradiations at 6 MV and 15 MV using multileaf collimators and hard wedges where necessary. Five clinical cases were studied: three lung and two breast cases. We found that, in terms of accuracy, the CCC algorithm performed better overall than AAA compared to VMC++, but AAA remains an attractive option for routine ...
The purpose of this work was to study and quantify the differences in dose distributions computed with some of the newest dose calculation algorithms available in commercial planning systems. The study was done for clinical cases where large density inhomogeneities were present. Three dose algorithms were used: a pencil beam like algorithm, the anisotropic analytic algorithm (AAA), a convolution superposition algorithm, collapsed cone convolution (CCC) and a Monte Carlo program, voxel Monte Carlo (VMC++). The dose calculation algorithms were compared under static field irradiations at 6 MV and 15 MV using multileaf collimators and hard wedges where necessary. Five clinical cases were studied: three lung and two breast cases. We found that the CCC algorithm performed overall better than AAA compared to VMC++, but AAA remains an attractive option for routine use in the clinic due to its short computation times. Dose differences between the ...
The Monte Carlo transport codes LAHET and MCNP were used to calculate energy fluence spectra at three neutron therapy facilities. The results compare very favourably with measured data. Kerma spectra and the ratio of ICRU muscle tissue kerma to A-150 kerma, along the carbon to oxygen kerma ratio, were determined. Absorbed dose rate calculations are in reasonable agreement with measured values. Use of these codes to study modifications to existing therapy beams is briefly discussed. (author).
We present our plans for a Monte-Carlo code simulating all possible combinations of (electromagnetic) interactions between colliding electron, positron, and both high-energy and laser photon beams, based on the ABEL code for beam-beam interaction. The implementation and first results for the laser-e{sup -} interaction are described. ((orig.)).
We present our plans for a Monte-Carlo code simulating all possible combinations of (electromagnetic) interactions between colliding electron, positron, and both high-energy and laser photon beams, based, on the ABEL code for beam-beam interaction. The implementation and first results for the laser-e"- interaction are described.
A microemulsion comprised of water, Brij 78, pentanol and styrene into which PPO and bis-MSB had been dissolved was prepared. Polymerization of the styrene resulted in a suspension of fluor-containing polystyrene nanoparticles (<100 nm). After a concentration step, the aqueous nanosuspension was able to detect {sup 14}C with counting efficiencies over 50% of those of a commercially available scintillation cocktail. Monte Carlo calculations demonstrated that the size and concentration of the nanoparticles were appropriate for optimum detection efficiency.
The theory of neutron moderation and spectroscopy are briefly reviewed, and moderators that are useful for Bonner sphere spectrometers are discussed. The choice of the neutron detector for a Bonner sphere spectrometer is examined. Spectral deconvolution methods are briefly reviewed, including derivative, parametric, quadrature, and Monte Carlo methods. Calibration is then discussed. (LEW)
Depth dose distributions for protons in the energy range from 2 to 600 MeV have been calculated by Monte Carlo and analytical methods. The dependence on energy of Maximum Dose Equivalent (MDE), quality factor and isotropic factor recommended for practical use are presented. The problems connected with regulation of maximum permissible proton fluxes are discussed. (author).
A model has been developed to predict the chemical composition of lakes that form in pyrite-bearing open-pit mines. Pyrite oxidation in wall rock releases sulfuric acid and metals, which mix with ground waters in the lake. Aqueous metal concentrations are calculated using MINTEQA2.
Integrated femtocell/macrocell networks, comprising a conventional cellular network overlaid with femtocells, offer an economically appealing way to improve coverage, quality of service, and access network capacity. The key element to successful femtocells/macrocell integration lies in its self-organizing capability. Provisioning of quality of service is the main technical challenge of the femtocell/macrocell integrated networks, while the main administrative challenge is the choice of the proper evolutionary path from the existing macrocellular networks to the integrated network. In this article, we introduce three integrated network architectures which, while increasing the access capacity, they also reduce the deployment and operational costs. Then, we discuss a number of technical issues, which are key to making such integration a reality, and we offer ...
This paper explores the feasibility of integrating energyefficiency program evaluation with the emerging need for the evaluationof programs from different "energy cultures" (demand response, renewableenergy, and climate change). The paper reviews key features andinformation needs of the energy cultures and critically reviews theopportunities and challenges associated with integrating these withenergy efficiency program evaluation. There is a need to integrate thedifferent policy arenas where energy efficiency, demand response, andclimate change programs are developed, and there are positive signs thatthis integration is starting to occur.
The radiative transport problem is cast in integral form using a transport kernel. The transport kernel has an explicit representation in terms of a Feynman Path Integral over all paths between selected points in a volume. This representation is setup in detail. Numerical evaluation of this Path Integral is formulated with a Frenet-Serret based procedure for generating valid random paths, and with a numerical evaluation of the weight for each valid path. Very early sanity checks of a numerical implementation are reported. Approaches to optimization are identified. (authors)
High-energy ionizing radiation is a prominent modality for the treatment of many cancers. The approaches to electron dose calculation can be categorized into semi-empirical models (e.g. Fermi-Eyges, convolution-superposition) and probabilistic methods (e.g. Monte Carlo). A third approach to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. We derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on Berthon et al (2007 J. Sci. Comput. 31 347-89) that exactly preserves the key properties of the analytical solution on the discrete level. We discuss several test cases taken from the medical physics literature. A test case with an ...
Background We developed a tissue-engineered biphasic cartilage bone substitute construct which has been shown to integrate with host cartilage and differs from autologous osteochondral transfer in which integration with host cartilage does not occur. Questions/purposes (1) Develop a reproducible in vitro model to study the mechanisms regulating tissue-engineered cartilage integration with host cartilage, (2) compare the integrative properties of tissue-engineered cartilage with autologous cartilage and (3) determine if chondrocytes from the in-vitro formed cartilage migrate across the integration site. Methods A biphasic construct was placed into host bovine osteochondral explant and cultured for up to 8 weeks (n?=?6 at each time point). Autologous osteochondral implants served as controls...
Recent supply chain reengineering efforts have focused on integrating firms? production, inventory and replenishment activities with the help of communication networks. While communication networks and supply chain integration facilitate optimization of traditional supply chain functions, they also exacerbate the information security risk: communication networks propagate security breaches from one firm to another, and supply chain integration causes breach on one firm to affect other firms in the supply chain. We study the impact of network security vulnerability and supply chain integration on firms? incentives to invest in information security. We find that even though an increase in either the degree of network vulnerability or the degree of supply chain integration increases the secur...
A nonviral vector for highly efficient site-specific integration would be desirable for many applications in transgenesis, including gene therapy. In this study we directly compared the genomic integration...Full Text Available
Transposons are promising systems for somatic gene integration because they can not only integrate exogenous genes efficiently, but also be delivered to a variety of organs using a range of transfection...Full Text Available
We use functional integral techniques to calculate the scattering amplitude for four open off-shell tachyons in Witten's string field theory and show that the residues of the first three poles agree with those obtained using oscillator methods. (orig.).
SYNOPSISRhomboid proteases are a fascinating class of enzymes that combine a serine protease active site within the core of an integral membrane protein. Despite having key roles...Full Text Available
Hepatitis B virus (HBV) is clearly a factor in the development of hepatocellular carcinoma, but its mechanism of action remains obscure. One possibility is that the HBV integration event alters the...Full Text Available
A theoretical framework centering on four classes of self-referent constructs is offered as a device for integrating the diverse areas constituting medical sociology. Guidance by this framework...Full Text Available
We use functional integral techniques to calculate the scattering amplitude for four open off-shell tachyons in Witten's string field theory and show that the residues of the first three poles agree with those obtained using oscillator methods.
This report highlights the Pipeline integrity Management methods being implemented by gas companies. These aim at maintaining the current high safety level, prevent major hazards, ensure the integrity of the pipeline and protect people and environment in the vicinity of the pipeline in the most cost effective way. It should be noticed that Pipeline Integrity Management aspects, technical and organisational, are included in the more general framework of the Safety Management System. Currently, more and more gas companies implement such a system on the basis of standards like ISO 9000 and so on. In this way, the report shows how practices of Pipeline Integrity Management are continually developing in order to adapt to their environment, and to improve performance. Past experience and imminent developments show that Pipeline Integrity Management is a flexible and efficient approach to ...
On the basis of the resolvent of a simple differential operator, a method for finding approximations to continuous functions is constructed. In this method, both the approximated function and its approximations satisfy the given integral boundary condition.
The developing Integrated Academic Information System (IAIMS) at Columbia-Presbyterian Medical Center provides data sharing links between two separate corporate entities, namely Columbia University...Full Text Available
This quarterly update contains information on the National Bioenergy Center Biochemical Platform Integration Project, R&D progress and related activities.
The integration sites for viral DNA in cells infected with Moloney murine leukemia virus (M-MuLV) were studied by restriction endonuclease cleavage of cellular DNA followed by electrophoresis in agarose...Full Text Available
We have developed miniature (≈1 μm diameter) microcavity surface-plasmon-resonance sensors (MSPRS), integrated them with microfluidics and tested...Full Text Available
We present a methodology and database mediator tool for integrating modern knowledge-based systems, such as the Stanford EON architecture for automated guideline-based decision-support, with legacy...Full Text Available
Objectives: The authors designed and implemented a clinical data mart composed of an integrated information retrieval (IR) and relational database management system (RDBMS).Design:...Full Text Available
The University of Cincinnati (UC) has been active in the National Library of Medicine's Integrated Advanced Information Management Systems (IAIMS) program since IAIMS' inception in 1984. UC received...Full Text Available
The AKR mouse strain is characterized by a high incidence of spontaneous thymic lymphoma that appears in older animals (greater than 6 months of age) and is associated with novel provirus integrations...Full Text Available
The plasmid pE194 (3.7 kilobases) is capable of integrating into the genome of the bacterial host Bacillus subtilis in the absence of the major homology-dependent RecE recombination system. Multiple...Full Text Available
The Kenya Revenue Authority (KRA) has been engaged in revenue administration reform for several years. Among the KRA's impressive, recent accomplishments is the integration of formerly separate departments for income taxes and value added taxes. (VAT) und...
It is usually very complex to design an Integrated Information System in a hospital that includes both in-patient and out-patient medical care data. The difficulty appears because each department...Full Text Available
An integrated system of micro-and mini-computers is described to acquire, analyze, store and report data on the total activity of the heart of a critically ill patient. Real-time beat and rhythm diagnoses...Full Text Available
ObjectiveTo estimate the joint effect of a multifaceted access intervention on primary care physician (PCP) productivity in a large, integrated prepaid group practice.Data...Full Text Available
A new method of random data analysis has been developed with special implications for membrane noise. The integral spectrometer uses overlapping broad-band filters of simple design, whose bandwidth...Full Text Available
The Accelerator Transmutation Waste System is under development at the Los Alamos National Laboratory. The goal is to perform an independent verification of the feasibility of actinide and long-lived fission product burning in this system. The authors' work is divided into five tasks: (a) production of an actinide and long-lived fission product cross section library from JEF 2.2; (b) simulation using MCNP and KENO IV Monte Carlo codes, of the Accelerator Transmutation Waste System configurations existing in literature; (c) validation of HETC Monte Carlo code (production of spallation source); (d) validation of the cross sections by comparison of Keff and reaction rate results, calculated with MNCP and KENO IV, with experimental benchmarks and intercomparison between the authors' calculations of a PWR unit cell and the computations carried out with various codes and cross section libraries (NEACRP critically group data); and (e) simulation of ...
There are two versions of the MORSE multigroup Monte Carlo radiation transport computer code system at Oak Ridge National Laboratory. MORSE-CGA is the most well-known and has undergone extensive use for many years. MORSE-SGC was originally developed in about 1980 in order to restructure the cross-section handling and thereby save storage. However, with the advent of new computer systems having much larger storage capacity, that aspect of SGC has become unnecessary. Both versions use data from multigroup cross-section libraries, although in somewhat different formats. MORSE-SGC is the version of MORSE that is part of the SCALE system, but it can also be run stand-alone. Both CGA and SGC use the Multiple Array System (MARS) geometry package. In the last six months the main focus of the work on these two versions has been on making them operational on workstations, in particular, the IBM RISC 6000 family. A new version of SCALE for workstations is being released to ...
One of the most active areas of physics in the last decades has been that of critical phenomena, and Monte Carlo simulations have played an important role as a guide for the validation and prediction of system properties close to the critical points. The kind of phase transitions occurring for the Betts lattice (lattice constructed removing 1/7 of the sites from the triangular lattice) have been studied before with the Potts model for the values q=3, ferromagnetic and antiferromagnetic regime. Here, we add up to this research line the ferromagnetic case for q=4 and 5. In the first case, the critical exponents are estimated for the second order transition, whereas for the latter case the histogram method is applied for the occurring first order transition. Additionally, Domany's Monte Carlo based clustering technique mainly used to group genes similar in their expression levels is reviewed. Finally, a control theory tool --an adaptive observer-- ...
We investigate the effects of single-particle structure and pairing on the equilibration of positive and negative-parity level densities for the even-even nuclei "5"8","6"2","6"6Fe and "5"8Ni and the odd-A nuclei "5"9Ni and "6"5Fe. Calculations are performed using the shell model Monte Carlo method in the complete fp-gds shell-model space using a pairing+quadrupole type residual interaction. We find for the even-even nuclei that the positive-parity states dominate at low excitation energies due to strong pairing correlations. At excitation energies at which pairs are broken, single-particle structure of these nuclei is seen to play the decisive role for the energy dependence of the ratio of negative-to-positive parity level densities. We also find that equilibration energies are noticeably lower for the odd-A nuclei "5"9Ni and "6"5Fe than for the neighboring even-even nuclei "5"8Ni and "6"6Fe.
We calculate parity-dependent level densities for the even-even isotopes "5"8","6"2","6"6Fe and "5"8Ni and the odd-A nuclei "5"9Ni and "6"5Fe using the shell model Monte Carlo method. We perform these calculations in the complete fp-gds shell-model space using a pairing+quadrupole residual interaction. We find that, due to pairing of identical nucleons, the low-energy spectrum is dominated by positive parity states. Although these pairs break at around the same excitation energy in all nuclei, the energy dependence of the ratio of negative-to-positive parity level densities depends strongly on the particular nucleus of interest. We find equilibration of both parities at noticeably lower excitation energies for the odd-A nuclei "5"9Ni and "6"5Fe than for the neighboring even-even nuclei "5"8Ni and "6"6Fe.
Let $\\mathscr{P}(E)$ be the space of probability measures on a measurable space $(E,\\mathcal{E})$. In this paper we introduce a class of nonlinear Markov chain Monte Carlo (MCMC) methods for simulating from a probability measure $\\pi\\in\\mathscr{P}(E)$. Nonlinear Markov kernels (see [Feynman--Kac Formulae: Genealogical and Interacting Particle Systems with Applications (2004) Springer]) $K:\\mathscr{P}(E)\\times E\\rightarrow\\mathscr{P}(E)$ can be constructed to, in some sense, improve over MCMC methods. However, such nonlinear kernels cannot be simulated exactly, so approximations of the nonlinear kernels are constructed using auxiliary or potentially self-interacting chains. Several nonlinear kernels are presented and it is demonstrated that, under some conditions, the associated approximations exhibit a strong law of large numbers; our proof technique is via the Poisson equation and Foster--Lyapunov conditions. We investigate the performance of our ...
High resolution Monte-Carlo simulations show that the neutron spectrum, fuel burnup and fuel temperature feedback effect of a PWR core loaded with Thoria-based fuel (Th/Pu-O_2) do not significantly differ from the MOX fuelled one due to the similar neutronic characteristics of both fertile materials (Th-232, U-238). The core physics of this fuel variant is characterized by an enhanced moderator/void temperature coefficient (by factor 2.4) and high incineration rate for Pu (approx. 60 %). A PWR core loaded with the Molybdenum-based inert matrix fuel (IMF) - in contrast to MOX-, shows a harder spectrum, resulting in small temperature coefficients of reactivity and particularly in a higher fuel depletion rate as well as an enhanced TRU reduction performance. The incineration of Pu amounts to 46 % resulting, in turn, in generation of minor actinides of about 10 % of the total Pu consumption. The higher excess reactivity resulting from the initial Pu contents is ...
Monte-Carlo simulations of 2-MeV #alpha#-particle channeling in Si_1_-_xSn_x alloys with 0#<=#x#<=#1 have been performed. The simulations are compared with measured channeling-angular scans for strained Si_0_._9_5Sn_0_._0_5 layers grown by molecular beam epitaxy (MBE). Agreement between simulated and measured angular scans can only be achieved if we assume a deviation of the crystal structure from the ideal one. This deviation can be attributed to a mosaic structure in the films and to an atomic-scale distortion of the crystal lattice due to an expected difference in the bond lengths between the Si-Si, Si-Sn and Sn-Sn atoms (such a difference in bond lengths has been observed in the epitaxial Si_1_-_xGe_x system). The contributions from both of these imperfections are estimated and discussed.
A technique capable of measuring the electron energy fluence spectra in a scattering medium was designed. These measurements were performed by setting a bremsstrahlung conversion target on the surface of a phantom, at an intermediate depth, and at a depth equal to electron mean range. The bremsstrahlung produced by the deceleration of electrons in the target was passed through an air channel in the phantom and passed forward by a pinhole collimator into a Na(Tl) detector. The measured pulse height data were unfolded to correct for the distortion of the spectrometer system by using the FORIST unfolding code. The unfolded bremsstrahlung spectra represent the electron energy fluence spectra convolution with the bremsstrahlung produced in the target. To generate the electron energy fluence spectra, the unfolded bremsstrahlung spectra were deconvoluted by using the MAZE2 unfolding code. CYLTRAN, a coupled electron-photon Monte Carlo transport code, was used to produce ...
In this paper, we report kinetic Monte Carlo study on the diffusion behavior of boron in silicon crystal, more particularly on the transient enhanced diffusion (TED) of boron in silicon during implantation and annealing. Firstly, the accuracy of our KMC code was verified by investigating the time evolutionary behavior of interstitial (I) and vacancy (V) when a silicon substrate is implanted with silicon dose with an energy of 10 keV and with a dose of 1 X 10{sup 14} ions/cm{sup 2}. To investigate the influence of native defects (I, V) on boron diffusion, a single and multi boron markers grown by MBE were employed. The simulation results revealed that the precursor of boron cluster (BI{sub 2}) is dominant at the initial stage of annealing, which explains the boron TED phenomenon in terms of the concentration of boron complexes and I, V clusters, respectively. The formation of {l_brace}311{r_brace} defects and dislocation loop were observed from the time-evolutionary ...
In this paper, we report kinetic Monte Carlo study on the diffusion behavior of boron in silicon crystal, more particularly on the transient enhanced diffusion (TED) of boron in silicon during implantation and annealing. Firstly, the accuracy of our KMC code was verified by investigating the time evolutionary behavior of interstitial (I) and vacancy (V) when a silicon substrate is implanted with silicon dose with an energy of 10 keV and with a dose of 1 X 10"1"4 ions/cm"2. To investigate the influence of native defects (I, V) on boron diffusion, a single and multi boron markers grown by MBE were employed. The simulation results revealed that the precursor of boron cluster (BI_2) is dominant at the initial stage of annealing, which explains the boron TED phenomenon in terms of the concentration of boron complexes and I, V clusters, respectively. The formation of #left brace#311#right brace# defects and dislocation loop were observed from the time-evolutionary study ...
Employing a recently developed Monte-Carlo model, we study the fission of 240Pu induced by neutrons with energies from thermal to just below the threshold for second chance fission. Current measurements of the mean number of prompt neutrons emitted in fission, together with less accurate measurements of the neutron energy spectra, place remarkably fine constraints on predictions of microscopic calculations. In particular, the total excitation energy of the nascent fragments must be specified to within 1MeV to avoid disagreement with measurements of the mean neutron multiplicity. The combination of the Monte-Carlo fission model with a statistical likelihood analysis also presents a powerful tool for the evaluation of fission neutron data. Of particular importance is the fission spectrum, which plays a key role in determining reactor criticality. We show that our approach can be used to develop an estimate of the fission spectrum with ...
Monte Carlo calculations using the codes PENELOPE and GEANT4 have been performed to characterize the dosimetric parameters of the new 20 mm long catheter based $^{32}$P beta source manufactured by Guidant Corporation. The dose distribution along the transverse axis and the two dimensional dose rate table have been calculated. Also, the dose rate at the reference point, the radial dose function and the anisotropy function were evaluated according to the adapted TG-60 formalism for cylindrical sources. PENELOPE and GEANT4 codes were first verified against previous results corresponding to the old 27 mm Guidant $^{32}$P beta source. The dose rate at the reference point for the unsheathed 27 mm source in water was calculated to be $0.215 \\pm 0.001$ cGy s$^{-1}$ mCi$^{-1}$, for PENELOPE, and $0.2312 \\pm 0.0008$ cGy s$^{-1}$ mCi$^{-1}$, for GEANT4. For the unsheathed 20 mm source these values were $0.2908 \\pm 0.0009$ cGy s$^{-1}$ mCi$^{-1}$ and $0.311 \\pm 0.001$ cGy ...
In the case wherein nonlinear seismic response analyses are carried out, the response values vary due to the variations in materials and modeling. In this paper, nonlinear analyses of several random variables are carried out using: i. a conventional method; ii. a two-point estimation method (i. and ii. are simplified methods); and iii. Monte Carlo simulation (detailed method) to examine the variability of the response in the excessive nonlinear range for seismic responses of shear walls. The analyses are performed to a PWR-3 loop type reactor building which is one of the most typical reactor buildings in Japan. The variations are considered in specified compressive strength of concrete, concrete damping factor, shear wave velocity of soil and shapes of shear stress-strain relation curves of shear walls. As the results by the two simplified methods closely matched the Monte Carlo simulation results, the appropriateness for applying the ...
Ion implantation and subsequent annealing are essential stages in today's advanced CMOS processing. Although the dopant implanted profile can be accurately predicted by analytical fits calibrated with SIMS profiles, the damage has to be estimated with a binary collision approximation implant simulator. Some models have been proposed, like the '+n', in an attempt to simplify the anneal simulation. We have used the atomistic kinetic Monte Carlo dados to elucidate which are the implant modeling features most relevant in the simulation of transient enhanced diffusion (TED). For the experimental conditions studied we find that the spatial correlation of the I, V Frenkel pairs is not critical in order to yield the correct I supersaturation, that can be simulated just taking into account the net I-V excess distribution. In contrast to, simulate impurity clustering/deactivation when there is an impurity concentration comparable to the net I-V excess, the full I and V ...
This study examines the degree of price-integration of equity indices between the major markets of Africa, namely Morocco, Tunisia, Egypt, Kenya, Nigeria, Namibia and South Africa with the European markets of London and Paris. Vector Autoregressive and Autoregressive Distributed Lag methods reveal that African markets are largely price-segmented. The only markets that are price-integrated have shared economic and financial institutions, such as Namibia and South Africa, and Egypt, Tunisia and France. The evidence suggests that development policy should be focussed on enhancing existing institutions rather than embarking prematurely on regional integration.
An intemted approach to building a networking infrastructure is an absolute necessiry for meeting the multidisciplinary science networking requirements of ...
The structural integrity of the Fuel Test Loop(FTL) in a Korean experimental reactor is evaluated when the FTL, submerged in a water environment, is subjected to a postulated pipe rupture. The analyses are performed under static and dynamic conditions, imposing the thrust force history at each postulated pipe rupture section. Through analysis the following results are found: 1) A double ended guillotine can not be expected based on the toughness of the material, 2) the structural integrity of the chimney surrounding the FTL would not impede the structural integrity by the pipe whip. All analyses are performed by finite element methods.
The conference paper gives a short description of the strategies within Statoil Gas Transport for pipeline condition monitoring and control. This is done for ensuring that the pipeline integrity can be maintained during all phases from construction to operation. Efficient condition monitoring and control can only be performed if the total pipeline integrity is taken into consideration such as design conditions, experience, environmental information, and general improved knowledge or the like for determining the requirement of inspection. Such philosophy will ensure the pipeline integrity better as corrective actions can be taken at a very early stage before any damage to the pipeline system have developed. 1 fig.
... the application of an acousto- optical tunable filter ... Couplers for Large Switch-Array Applications ... Symmetric Integrated Optic X Junction," Electronics ...
... DARPA UUV and the UUV support simulator. ... Modular connections for robot arm manipulators ... Intelligent mobile robots will be performing many ...
Prior to this at TRW, Mr. Schier managed flight software development on the MILSTAR Communications Payload and led integration and verification work on ...
Knowledge integration is the process of incorporating new information into a body of existing knowledge. It involves determining how new and existing knowledge interact and how existing knowledge should be modified to accommodate the new information. KI is a machine learning program that performs knowledge integration. Through actively investigating the interaction of new information with existing knowledge KI is capable of detecting and exploiting a variety of diverse learning opportunities during a single learning episode. Empirical evaluation suggests that KI provides significant assistance to knowledge engineers while integrating new information into a large knowledge base.
... They comprise of a power model for the representation of the electric system, a fluid model that represents the cooling fluid flow, three layers of ...
... to integrate and synchronize SC efforts and capabilities, as noted in the SC Roadmap, the joint community needs to publish SC policy and expand ...
... order probabilistic modelling, multisensor monitoring of ... object recognition, image registration and point ... of the Berkeley Continuous Media Toolkit. ...
Brachytherapy is a radiotherapeutic modality that makes use of radionuclides to deliver a high radiation dose to a well-defined volume while sparing surrounding healthy structures. At the National Institute for Cancer Research of Genova a High Dose Rate remote afterloading system provides Ir(192) endocavitary brachytherapy treatments. We studied the possibility to use the Geant4 Monte Carlo simulation toolkit in brachytherapy for calculation of complex physical parameters, not directly available by experiment al measurements, used in treatment planning dose deposition models.
A calorimetric-time-of-flight (CTOF) technique was used for real-time, high-precision measurement of neutron spectrum at the angle of 175 degrees from the initial proton beam direction, which hits a face plane of a cylindrical iron target of 20 cm in diameter and 25 cm thick. A comparison was performed between the neutron spectra predicted by the MARS and the MCNPX code systems and measured by experiments for 0.4-, 0.6-, 0.8-, 1.0-, and 1.2-GeV protons.
The present paper considers the evaluation of timber structures with the focus on robustness due to connection ductility. The robustness analysis is based on the structural reliability framework applied to a simplified mechanical system. The structural timber system is depicted as a parallel system. An evaluation method of the ductile behaviour is introduced. For different ductile behaviours, the system reliability is estimated based on Monte Carlo simulation. A correlation between the strength of the structural elements is introduced. The results indicate that the reliability of a structural timber system can be significantly increased due to the ductile behavior.
We investigate the formation of nanostructures in 2D strained alloys on face centered cubic (111) surfaces by means of equilibrium Monte Carlo simulations. In the framework of an off-lattice model, we consider one monolayer of two bulk-immiscible adsorbates A and B with negative and positive misfit relative to the substrate, respectively. Simulations show that the adsorbates partly self-organize into island or stripe-like patterns. We show how these structures depend on the relative misfits, interaction, and concentration of components. The morphology is quite different for phase separation and intermixing regimes.
Tests for the applicability of a CFD (Computational Fluid Dynamics) code for simulating activity transport in PNA (Pulsed Neutron Activation) fluid measurements have been performed. The CFD code was combined with a Monte Carlo code used for the calculation of the initial activity distribution. The results from the calculations show that it is possible to use CFD for calculation of the activity distribution in PNA. The mainly qualitative results in this work are encouraging and suggest further work. In the continuation of this work a response function for the gamma detector will be calculated so that a PNA time spectrum can be simulated. A more accurate comparison with experimental data can then be performed
An algorithm is proposed which treats the ECH generated potential in a self-consistent way, by extending the Monte-Carlo Fokker-Planck method used by Murakami [S. Murakami et al., Proc. 17th IAEA Fusion Energy Conference, Yokohama, 1998 (International Atomic Energy Agency, Vienna, in press), paper CN-69/TH2/1]. The additional physics is expected to influence the transport of both thermal and suprathermal electrons in a helical toroidal system. (author)
An important unresolved issue in outcomes analysis for lung complications is the effect of poor or completely lacking heterogeneity corrections in previously archived treatment plans. To estimate this effect, we developed a novel method based on Monte Carlo (MC) dose calculations which can be applied retrospectively to RTOG/AAPM-style archived treatment plans (ATP). We applied this method to 218 archived nonsmall cell lung cancer lung treatment plans that were originally calculated either without heterogeneity corrections or with primitive corrections. To retrospectively specify beam weights and wedges, beams were broken into Monte Carlo-generated beamlets, simulated using the VMC++ code, and mathematical optimization was used to match the archived water-based dose distributions. The derived beam weights (and any wedge effects) were then applied to Monte Carlo beamlets regenerated based on the patient computed tomography ...
By using a model independent Monte Carlo approach, we study the possible structure of charged and neutral lepton mass matrices, under the assumption of an U(2) horizontal symmetry (additional to the usual Standard Model ones) involving the light fermion generations. We assume the most general Majorana mass matrix for neutrinos. We update the results of our previous similar study, by inserting in the analysis the recent KamLAND data, that contributed to find a final solution to the Solar neutrino problem. The introduction of the new experimental data reduce the allowed regions in the nine dimensional space parameters, and show that our procedure gives stable solutions.
An accelerator experiment was performed using a low-energy antiproton beam to measure antiproton detection efficiency of BESS, a balloon-borne spectrometer with a superconducting solenoid. Measured efficiencies showed good agreement with calculated ones derived from the BESS Monte Carlo simulation based on GEANT/GHEISHA. With detailed verification of the BESS simulation, the relative systematic error of detection efficiency derived from the BESS simulation has been determined to be {+-}5%, compared with the previous estimation of {+-}15% which was the dominant uncertainty for measurements of cosmic-ray antiproton flux.
Recent developments in the physical model of 1 MeV to 100 TeV hadron and lepton interactions with nuclei and atoms are described. These include a new nuclear cross section library, a model for soft pion production, the cascade-exciton model, the dual parton model, deuteron-nucleus and neutrino-nucleus interaction models, detailed description of mu, pi and anti p absorption and a unified treatment of muon and charged hadron electromagnetic interactions with matter. New algorithms are implemented into the MARS13(98) Monte Carlo code and benchmarked against experimental data. The code capabilities to simulate cascades and generate a variety of results in complex media have been also enhanced.
multi-parton interactions and underlying event structure, $\\gamma$-mediated processes and luminosity determination, Monte Carlo tuning and even MSSM Higgs discovery in central exclusive production. In this article, the forward detector instrumentation around the CMS interaction point is described and the prospects for diffractive and forward physics using the CMS forward detectors are summarized. In addition, first observation of forward jets as well as early measurements of the forward energy flow in the pseudorapidity range $3.15
The adsorption and corrosion inhibition behavior of three selected thiosemicarbazone derivatives, namely 3-pyridinecarboxaldehyde thiosemicarbazone (META), isonicotinaldehyde thiosemicarbazone (PARA) and 2-pyridinecarboxaldehyde thiosemicarbazone (ORTHO) at the nickel surface were studied electrochemically by Tafel and impedance methods and computationally by carrying out Monte Carlo searches of configurational space on nickel/thiosemicarbazone derivative system. Electrochemical measurements showed that the inhibition efficiency of these compounds increased with increase in their concentration. The recorded inhibition efficiencies of the three tested thiosemicarbazone increase in the order: META>PARA>ORTHO. Polarization studies showed that these compounds act as mixed-type inhibitors for n...
A neutron collimator system is being designed as part of a neutron imaging facility for computed tomography and real-time neutron radiography research at the through beam port of the University of Texas TRIGA reactor. Lack of sufficient information about collimator systems in a through port from the literature necessitated the use of Monte Carlo calculations using the MCNP code 3 to search for optimal design configuration and materials that maximize the thermal neutron intensity at the image plane while minimizing the fast neutrons and gamma radiation.
The CyberKnife Robotic Radiosurgery System (Accuracy Incorporated, Sunnyvale, CA, USA) is used worldwide to treat tumors and neurological disorders anywhere in the body with sub-millimetre beam delivery accuracy. Accuracy has developed a number of new technologies in recent years to enhance the treatment of cancer patients. Such new technologies include a fast Monte Carlo Dose Calculation algorithm, Sequential Optimization dose planning, the IrisTM Variable Aperture Collimator, an 800 MU/min Linear Accelerator, and Optimized Path Traversal. These technologies enable physicists and physicians to plan treatments quickly and easily and deliver them with unrivalled accuracy and precision
Continuum events represent an eminent source of background in any e+e- experiment. As these have a higher branching ratio than BB-bar events (at BaBar this ratio is estimated to about 3.5) or ?+?- events, efficient continuum background suppression is essential in many analyses. Using Artificial Neural Networks and the Nearest Neighbor Method we developed several selectors which, based only on the global event shape variables, efficiently tag BB-bar events and ?+?- events against the continuum background. These selectors could then be combined with the channel specific information in various types of analyses. The study was done using a parametric Monte Carlo.
The diamond norm measures the distance between two quantum channels. From an operational viewpoint, this norm measures how well we can distinguish between two channels by applying them to the input states of arbitrarily large dimensions. In this paper, we show that the diamond norm can be conveniently, and in a physically transparent way, computed by means of a Monte Carlo algorithm based on the Fano representation of quantum states and quantum operations. The effectiveness of this algorithm is illustrated for several single-qubit quantum channels.
An extensive set of benchmark measurement of PDDs and beam profiles was performed in a heterogeneous layer phantom, including a lung equivalent heterogeneity, by means of several detectors and compared against the predicted dose values by different calculation algorithms in two treatment planning systems. PDDs were measured with TLDs, plane parallel and cylindrical ionization chambers and beam profiles with films. Additionally, Monte Carlo simulations by meansof the PENELOPE code were performed. Four different field sizes (10x10, 5x5, 2x2, and1x1 cm"2) and two lung equivalent materials (CIRS, #rho#_e"w=0.195 and St. Bartholomew Hospital, London, #rho#_e"w=0.244-0.322) were studied. The performance of four correction-based algorithms and one based on convolution-superposition was analyzed. The correction-based algorithms were the Batho, the Modified Batho, and the Equivalent TAR implemented in the Cadplan (Varian) treatment planning system and the TMS Pencil Beam ...
We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent results in terms of shock hydrodynamics, cosmic ray spectrum and also escaping flux spectrum and anisotropy. Strong points and limits of each approach are also discussed, as a function of the problem one wants to study.
This handbook breaks down the complex field of nuclear reactor calculations into major steps. Each step presents a detailed analysis of the problems to be solved, the parameters involved, and the elaborate computer programs developed to perform the calculations. This book bridges the gap between nuclear reactor theory and the implementation of that theory, including the problems to be encountered and the level of confidence that should be given to the methods described. Volume II: Monte Carlo Calculations for Nuclear Reactors. In-Core Management of Four Reactor Types. In-Core Management in CANDU-PHW Reactors. Reactor Dynamics. The Theory of Neutron Leakage in Reactor Lattices. Index.
We made Basic designs of neutron guide, choppers, and detectors in order to optimize the design parameters of DC-TOF to be built in the HANARO Cold Neutron Guide Hall. In addition, we calculated the expected performance of DC-TOF using Monte Carlo simulations and evaluated the properties of neutron beam. Based on the results we obtained, we have compared the expected performance of the DC-TOF with those of existing instruments overseas. In conclusion, we believe that we will be able to construct the DC-TOF at HANARO as one of the best instruments of its kinds and it will become an invaluable instrument to researchers in the related field.
An accelerator beam experiment was performed using a low-energy antiproton beam to measure antiproton detection efficiency of the BESS detector. Measured and calculated efficiencies derived from the BESS Monte Carlo simulation based on GRANT/GHEISHA showed good agreement. With detailed verification of the BESS simulation, the relative systematic error of detection efficiency derived from the BESS simulation has been determined to be {+-}5%, compared with the previous estimation of {+-}15% which was the dominant uncertainty for measurements of cosmic-ray antiproton flux. (author)
We have developed a new type personal dosemeter by using a B-10 doped silicon p-n junction detector with a polyethylene radiator and a polyethylene moderator. The purpose of this study was to develop a real time neutron dosemeter with a nearly flat response in the energy range from thermal to 15 MeV and low angular dependence to the incident neutron direction. The neutron response of the dosemeter was obtained with the Monte Carlo calculation and the monoenergetic neutron experiment in a free air field and also under a condition attached on a phantom.
Symplectic integrators designed for simulating soft matter at constant temperature and constant pressure (or constant surface tension) are presented. In addition to the well-known merits of symplectic integrators, such as long time stability, these methods allow the estimation of thermal heat production/absorption and thus entropy differences through a simple function of the time scaling factor ?. The relative entropy estimation around crystal-liquid phase transition is given for systems with soft core repulsive potential.
A renormalization procedure is proposed which applies to lattice Feynman integrals containing zero-mass propagators and is analogous to the BPHZL renormalization procedure for continuum Feynman integrals. The renormalized diagrams are infrared convergent for non-exceptional external momenta, if the vertices of the theory satisfy a general infrared constraint. Under the same conditions as in the massive case, the continuum limit of the renormalized theory exists and is independent of the details of the lattice action.
The aim of this report is to show how Energetic Planning and Territorial Policy should be working together for a better integration of Renewable Energies into Region. This Integration should to contemplate social, economic and environmental aspects of the territory. The report has been classified into 7 items: planning, energetic scenarios, technology transfer for Renewable Energies dissemination, barriers for this dissemination, environmental aspects, European Union Policy and Decision Support Systems (and specially GIS). (Author) 54 refs.
The computer facilities and the computer expertise at certain departments at Aalborg University made these departments natural starting points for integration of Information and Communication Technology (ICT) in their activities. This utilization and integration of ICT was not organized from the top or intermediate levels. It developed from the grassroots. Today Aalborg University as a whole can be said to be a truly ICT integrated institution. The quick development is to a high extent due to a structure with a finely distributed system of equipment, know-how and autonomy at all levels
Accelerated aging tests such as high temperature burn-in, which are in current use on Bendix Kansas City Division's (BKC) purchased small and medium scale integrated circuits, were evaluated to determine if they are effective and necessary to ensure the required reliability. A theoretical analysis, a literature search, and a study of lot acceptance results were used to assess the value of integrated circuit burn-in.
Using the manifestly causal gluon propagator in the light-cone gauge, we evaluate two one-loop Feynman integrals which appear in the computation of the three-gluon vertex correction in the two-component formalism of the Yang-Mills gauge fields. We conclude that they both are finite for {omega} -> 2, where {omega} is the dimensional regulator. They belong to a class of finite one-loop light-cone integrals. (author). 7 refs.
For the Convolution, Superposition and Monte Carlo methods, energy spectrum data is an essential parameter. Especially, by using Monte Carlo method, if the adequate properties of electrons are set, complex energy spectrum data can be calculated by simple input. However, adjustment of the property of electrons incident on the target of linear accelerator is time-consuming procedure. To cover the difference of many accelerators, it is important how we obtain accurate energy spectrum by simple methods. Especially, the total accuracy of dose calculation depends on the agreement of measured and calculated percentage depth dose (PDD) and off-axis ratio (OAR). Thus, a simple method for the determination of the simulation parameters was considered. To reduce the adjustment procedure, we selected values of important three parameters, which were the mean energy of the input electron beam, the energy distribution of the input electron beam and the radial ...
Accurate modeling of system response and scatter distribution is crucial for image reconstruction in emission tomography. Monte Carlo simulations are very well suited to calculate these quantities. However, Monte Carlo simulations are also slow and many simulated counts are needed to provide a sufficiently exact estimate of the detection probabilities. In order to overcome these problems, we propose to split the simulation into two parts, the detection system and the object to be imaged (the patient). A so-called 'virtual boundary' that separates these two parts is introduced. Within the patient, particles are simulated conventionally. Whenever a photon reaches the virtual boundary, its detection probability is calculated analytically by evaluating a multi-dimensional B-spline that depends on the photon position, direction and energy. The unknown B-spline knot values that define this B-spline are fixed by a prior 'pre-' simulation that needs to ...
Brachytherapy refers to short distance treatment of cancer with radiation from small, encapsulated radionuclide sources. This type of treatment is given by placing sources directly into or near the volume to be treated. The dose is then delivered continuously, either over a short period of time (temporary implants) or over the lifetime of the source to a complete decay (permanent implants). This technique is extensively used in the treatments of gynecological, breast, prostate, head and neck, and other soft tissue cancers. The clinical outcomes of this treatment depend on accurate dosimetry and dose delivery. The patient is generally considered to be water equivalent. So the treatment planning system software generally is based on measured or calculated dose distributions in homogenous water medium. In real situation, the treatment volume is not homogenous water medium, but it is a heterogeneous medium consisting of bone, soft tissue, lung, air etc. For better understanding we present ...
A Siemens Mevatron KV2 accelerator installed at the Royal Adelaide Hospital employs cylindrical solid-walled electron cones for some electron collimation. The cones being used at present result in treatment fields that do not always conform with the International Electrotechnical Commission (IEC) Standards (particularly at high energies). The aim of this project was to simulate the existing cones using Monte Carlo methods in order to evaluate potential cone modifications required to overcome the field irregularities. Simulations were performed using the EGS4 (Electron Gamma Shower version 4, distribution II) Monte Carlo code installed on a DEC Alpha workstation at the University of South Australia. To rigorously simulate the existing electron cones it was necessary to also simulate various components within the treatment head of the linear accelerator. Results of simulations for existing cones were found to be consistent with experimental data. ...
For Monte Carlo simulations to be used as an alternative solution to perform scatter correction, accurate modelling of the scanner as well as speed is paramount. General-purpose Monte Carlo packages (Geant4, EGS, MCNP) allow a detailed description of the scanner but are not efficient at simulating voxel-based geometries (patient images). On the other hand, dedicated codes (SimSET, PETSIM) will perform well for voxel-based objects but will be poor in their capacity of simulating complex geometries such as a PET scanner. The approach adopted in this work was to couple a dedicated code (SimSET) with a general-purpose package (Geant4) to have the efficiency of the former and the capabilities of the latter. The combined SimSET+Geant4 code (SimG4) was assessed on the GE Advance PET scanner and compared to the use of SimSET only. A better description of the resolution and sensitivity of the scanner and of the scatter fraction was obtained with SimG4. ...
As for Wolsung NPP unit 1, radiation shielding analysis was performed by using neutron diffusion codes, one-dimensional discrete ordinates code ANISN, and analytical methods. But for Wolsung NPP unit 2, 3, and 4, two-dimensional discrete ordinates code DOT substituted for neutron diffusion codes. In other words, the method of analysis and computer codes used for radiation shielding of CANDU 6 type reactor have been improved. Recently Monte Carlo MCNP code has been widely utilized in the field of radiation physics and other radiation related areas because it can describe an object sophisticately by use of three-dimensional modelling and can adopt continuous energy cross-section library. Nowadays Monte Carlo method has been reported to be competitive to discrete ordinate method in the field of radiation shielding and the former has been known to be superior to the latter for complex geometry problem. However, Monte Carlo ...
This report is composed of the following two parts and appendix. (I) Improvement of the Method for Evaluating Reactivity Based on Monte Carlo Perturbation Theory: Theoretical formulation in Monte Carlo perturbation method had been checked, and then introduced into a calculation code. There are some cases that the results of the change of eigenvalues becomes positive or negative by changing the estimator, and there is no reasonable difference in the results between the conventional method, which does not consider the change of neutron source distribution caused by a perturbation, and the new method, which consider that change. Thus it is still necessary to check the Monte Carlo perturbation code. (II) Improvement of Nodal Transport Method for 3-D Hexagonal Geometry: We can accurately evaluate hexagonal geometry FBR core by nodal transport calculation code for hexagonal-Z geometry named `NSHEX`. However it is also found that ...
The hard-sphere radial distribution functions, g_H_S(r/d,#eta#), for the face-centered cubic and hexagonal close-packed phases have been computed by the Monte Carlo method at nine values of the packing fraction, #eta#[=(#pi#/6)#rho#d"3], ranging from 4% below the melting density to 99% of the close-packed density. The Monte Carlo data are used to improve available analytic expressions for g_H_S(r/d,#eta#). By utilizing the new g_H_S(r/d,#eta#) in the Henderson and Grundke method [J. Chem. Phys. 63, 601 (1975)], we next derive an expression for y_H_S(r/d,#eta#) [=g_H_S(r/d)exp#left brace##beta#V_H_S(r)#right brace#] inside the hard-sphere diameter, d. These expressions are employed in a solid-state perturbation theory [J. Chem. Phys. 84, 4547 (1986)] to compute solid-state and melting properties of the Lennard-Jones and inverse-power potentials. Results are in close agreement with Monte Carlo and lattice-dynamics ...
In the field of reactor and fuel cycle physics, particle transport plays and important role. Neutronic design, operation and evaluation calculations of nuclear system make use of large and powerful computer codes. However, current limitations in terms of computer resources make it necessary to introduce simplifications and approximations in order to keep calculation time and cost within reasonable limits. Two different types of methods are available in these codes. The first one is the deterministic method, which is applicable in most practical cases but requires approximations. The other method is the Monte Carlo method, which does not make these approximations but which generally requires exceedingly long running times. The main motivation of this work is to investigate the possibility of a combined use of the two methods in such a way as to retain their advantages while avoiding their drawbacks. Our work has mainly focused on the speed-up of 3-D continuous ...
In 2002 and again in 2003, an investigative journalist unit at ABC News transported a 6.8 kilogram metallic slug of depleted uranium (DU) via shipping container from Istanbul, Turkey to Brooklyn, NY and from Jakarta, Indonesia to Long Beach, CA. Targeted inspection of these shipping containers by Department of Homeland Security (DHS) personnel, included the use of gamma-ray imaging, portal monitors and hand-held radiation detectors, did not uncover the hidden DU. Monte Carlo analysis of the gamma-ray intensity and spectrum of a DU slug and one consisting of highly-enriched uranium (HEU) showed that DU was a proper surrogate for testing the ability of DHS to detect the illicit transport of HEU. Our analysis using MCNP-5 illustrated the ease of fully shielding an HEU sample to avoid detection. The assembly of an Improvised Nuclear Device (IND) -- a crude atomic bomb -- from sub-critical pieces of HEU metal was then examined via Monte Carlo ...
In this paper, we studied vapor-liquid equilibria (VLE) and adsorption of ethylene on graphitized thermal carbon black and in slit pores whose walls are composed of graphene layers. Simple models of a one-center Lennard-Jones (LJ) potential and a two-center united atom (UA)-LJ potential are investigated to study the impact of the choice of potential models in the description of VLE and adsorption behavior. Here, we used a Monte Carlo simulation method with grand canonical Monte Carlo (GCMC) and Gibbs ensemble Monte Carlo ensembles. The one-center potential model cannot describe adequately the VLE over the practical range of temperature from the triple point to the critical point. On the other hand, the two-center potential model (Wick et al. J. Phys. Chem. B 2000, 104, 8008-8016) performs well in the description of VLE (saturated vapor and liquid densities and vapor pressure) over the wide range of temperature. This UA-LJ ...
Using Monte Carlo simulation and the convolution/superposition algorithm, this work examines percent depth dose curves of the central axis in an acrylic phantom (20x20x20 cm"3) with variously sized air cavities (20x20x1.0, 20x20x2.0, 20x20x3.0, 20x20x4.0 and 20x20x4.95 cm"3 for study of longitudinal electron disequilibrium (ED) and 3.6x3.6x4.95, 4.5x4.5x4.95, 5.4x5.4x4.95 and 20x20x4.95 cm"3 for study of lateral ED). Radiochromic film samples are also measured to verify the Monte Carlo results. The Monte Carlo simulation is performed using OMEGA/BEAM and DOSXYZ codes, and the convolution/superposition calculation relies on an ADAC commercial treatment planning system. Underestimating the dose kernel expansion leads to overestimating the dose of what was found in the air cavity of ED using the convolution/superposition algorithm. Consequently, the dose in the rebuild-up region is influenced. The influenced region is on the ...
Summary 1.-The modules of plants experience different levels of resources, and clonal plants can integrate resource heterogeneity. However, no studies have tested whether the benefits of clonal integration depend on patterns of heterogeneity in multiple resources, like high levels of above-ground and below-ground resources in the same patches (-coincident patchiness-) or different patches (-reciprocal patchiness-). 2.-We hypothesized that the benefits of clonal integration can vary depending on whether patchiness is reciprocal or coincident, and that clonal species experience greater benefits from integration when qualitative patterns of resource heterogeneity are more like those likely to occur in their habitats. To test these hypotheses, we grew pairs of connected ramets of Cynodon dacty...
In this paper, we show how to improve and extend the integration by fractional expansion technique (IBFE) by applying it to certain families of scalar massive Feynman diagrams. The strategy is based on combining this method together with the integration by parts technique. In particular, we want to calculate certain Feynman diagrams which have a triangle loop as a subgraph. The main idea is to use the integration by parts technique in this subgraph in order to simplify the topology of the original diagram in which it is immersed, using then, in a second step, the IBFE technique. The result we have obtained, after the application of both techniques, represents a simplification in the complexity of the solution, compared with having used only the IBFE technique.
The integration of European energy markets is a key goal of EU energy policy, and has also been the focal point of many scientific studies in recent years. International markets for coal, oil, natural gas and electricity have previously been investigated in order to determine the extent of the respective markets. This study enhances this field of research to bioenergy markets. Price series data and time series econometrics are used to determine whether residential sector wood pellet markets of Austria, Germany and Sweden are integrated. The results of the econometric tests show that the German and Austrian markets can be considered to be integrated, whereas the Swedish market is separate from the other two countries. Although increased internationalization of wood pellet markets is likely to contribute to European price convergence and market integration, this process is far from completed. (author)
...3 MtC, with 78% of these emissions from integrated steel making (26.9 MtC) and the rest (22%)...6 MtC, with 71% of these emissions from integrated steel making (21.9 MtC) and the rest (29%)...1 million ton of carbon (MtC) in 1994 (3.9 MtC from integrated steel making, and 2.2 MtC from secondary ...7 MtC in 2002 (3.7 MtC from integrated steel making, and 2.0 MtC from secondary steel making), corresponding ...
Unitarity cuts are widely used in analytic computation of loop amplitudes in gauge theories such as QCD. We expand upon the technique introduced in hep-ph/0503132 to carry out any finite unitarity cut integral. This technique naturally separates the contributions of bubble, triangle and box integrals in one-loop amplitudes and is not constrained to any particular helicity configurations. Loop momentum integration is reduced to a sequence of algebraic operations. We discuss the extraction of the residues at higher-order poles. Additionally, we offer concise algebraic formulas for expressing coefficients of three-mass triangle integrals. As an application, we compute all remaining coefficients of bubble and triangle integrals for nonsupersymmetric six-gluon amplitudes.
This studys purpose was to explore the relationship between creativity and integrity in the ideation phase of feature filmmaking. Integrity refers to ones self-defined ability to maintain authenticity and moral autonomy while preserving ones sense of membership and loyalty to the team or organization. When team members choose elements for the screenplays story that they feel will attract the ideal audience, the dynamic tension between creativity and integrity is most apparent. The forces at play during this phase of work yielded the research question: In what ways and to what degree do screenwriters andor studio executives feel their personal integrity is in alignment with the creative process of feature film development? Several concepts from the literature formed the design around which ...
The accuracy of dose computation within the lungs depends strongly on the performance of the calculation algorithm in regions of electronic disequilibrium that arise near tissue inhomogeneities with large density variations. There is a lack of data evaluating the performance of highly developed analytical dose calculation algorithms compared to Monte Carlo computations in a clinical setting. We compared full Monte Carlo calculations (performed by our Monte Carlo dose engine MCDE) with two different commercial convolution/superposition (CS) implementations (Pinnacle-CS and Helax-TMS's collapsed cone model Helax-CC) and one pencil beam algorithm (Helax-TMS's pencil beam model Helax-PB) for 10 intensity modulated radiation therapy (IMRT) lung cancer patients. Treatment plans were created for two photon beam qualities (6 and 18 MV). For each dose calculation algorithm, patient, and beam quality, the following set of clinically ...
A path integral evaluation of the Green's function for the hydrogen atom initiated by Duru and Kleinert is studied by recognizing it as a special case of the general treatment of the separable Hamiltonian of Liouville type. The basic dynamical principle involved is identified as Jacobi's principle of least action for given energy which is reparametrization invariant, and thus the appearance of a gauge freedom is naturally understood. The separation of variables in the operator formalism corresponds to a choice of gauge in the path integral, and the Green's function is shown to be gauge independent if the operator ordering is properly taken into account. Unlike the conventional Feynman path integral, which deals with a space-time picture of particle motion, the path integral on the basis of Jacobi's principle sums over orbits in space. We illustrate these properties by evaluating an exact path ...
Throughout the 1990's to year 2000, energy integration, and particularly natural gas integration, was seen as a major goal in the southern cone of South America. The regional perspective for energy was related to an even more challenging objective - Latin American economic integration, which was quickly moving forward by the constitution of Mercosur, the free trade zone built by Argentina, Brazil, Uruguay and Paraguay, and also having Bolivia and Chile as special partners. Mercosur was getting stronger and becoming a successful political project. Trade conflicts existed, but they were considered normal and nobody would disagree about the viability of Mercosur. By the beginning of the 21. century, the economic integration process in South America was already advancing timidly as compared to the impetus from the beginning of 1990's. Mercosur started facing complex and almost insurmountable challenges ...
Daylight responsive dimming systems have been used in few buildings to date because they require improvements to improve reliability. The key underlying factor contributing to poor performance is the variability of the ratio of the photosensor signal to daylight workplane illuminance in accordance with sun position, sky condition, and fenestration condition. Therefore, this paper describes the integrated systems between automated roller shade systems and daylight responsive dimming systems with an improved closed-loop proportional control algorithm, and the relative performance of the integrated systems and single systems. The concept of the improved closed-loop proportional control algorithm for the integrated systems is to predict the varying correlation of photosensor signal to daylight workplane illuminance according to roller shade height and sky conditions for improvement of the system accuracy. In this study, the ...
The Project ``Middle School Integrated Science, Mathematics and Technology Curriculum`` had two goals: (1) to survey the literature of energy education; and (2) to develop a theme for a possible integrated middle school energy based curriculum. We aimed to respond to the challenge of developing thematic integrated curricula as advocated by the NSTA, AAAS and other organizations analyzing the future of American science and mathematics education. The survey of middle school energy curriculum materials has been completed. A list of the resources surveyed are included in this report. Though many energy based curriculum materials have been produced, none of them appears to be broadly disseminated throughout the country. Some energy based curriculum materials are far less well developed than others. We found that an integrated set of modular materials concerning the energy based theme of light and optics does ...
A case study is presented where an integrated, ecologically relevant exposure assessment is presented for a polluted lowland river. Using partial least squares regression of latent structures (PLS), an analysis of the impact of two effluents on physico-chemical water quality measures, macroinvertebrate and diatom communities, and in situ bioassay responses with four different test species are combined into an integrative exposure assessment. Bioassays focussed on growth and condition related endpoints, because they are key functional processes of organisms and populations. Integrating these multiple lines of evidence, we were able to discriminate among the impact of both effluents, link changes in physico-chemical water quality with bioassay endpoints and ecological quality of the ecosystem, and address the importance of integrating all information into one exposure assessment framework. The bioassays ...
Integration of positron emission tomography (PET) and magnetic resonance (MR) has become a topic of increasing interest to the imaging community over the past two years. In this text, the authors attempt to distinguish facts from fiction concerning such integrated systems. Analysis of existing information of combined imaging on existing brain PET/MR systems and imaging experience with PET-computed tomography (CT) is reviewed. Various types of system integration of PET and MR are discussed with completely independent systems on one hand and completely integrated systems with the possibility of simultaneous data acquisition on the other hand. Furthermore, it is discussed, what simultaneous data acquisition with nuclear imaging systems combined with MR or CT really means, as technical simultaneity may not be relevant in light of the pharmacokinetics of the nuclear tracers used. The authors conclude that ...
The effective integration of technology into the teaching and learning of mathematics remains one of the critical challenges facing contemporary tertiary mathematics. This article reports on some significant findings of a wider study investigating the use of technology in undergraduate mathematics. It first discusses a taxonomy developed to describe and compare technology use within individual courses and departments that identifies a complex range of factors, summarized under six defining characteristics of an integrated technology mathematics curriculum (ITMC). An instrument for a simple comparison of technology use employing the elements of this taxonomy is provided. It then presents evidence gathered from an observational study of technology implementation at The University of Auckland...
The two-dimensional WKB model has been the basis for several investigations of the quasi-Landau hydrogenic spectrum. Whereas other authors have used numerical integration, it is shown that the results can be obtained analytically in terms of elliptic integrals. The later are easily generated by even programmable pocket calculators, from which numerical results - which are in good agreement with experiments - are easily obtained. A further advantage of using elliptic integrals is that several strong-field problems can be shown to be special cases of a general potential whose energy and spacing are expressible in terms of them.
To improve operational safety of the fuel channel for Wolsung units the recommended primary heat transport system (PHTS) operating guidelines were developed, which provides a safe operation window in terms of pressure and temperature to maintain fuel channel integrity, which is similar to pressure-temperature curves used for primary pressure boundary integrity in PWRs and deterministic Leak-Before Break (LBB) analyses were conducted using recent understanding of delayed hydride cracking (DHC) mechanism. In addition, a comparative assessment of LBB analysis was made using the procedures of Fitness For Service Guidelines(FFSG) for pressure tubes rev. 0 and rev. 1.
Ascochyta blight causes significant yield loss in pulse crops worldwide. Integrated disease management is essential to take advantage of cultivars with partial resistance to this disease. The most effective practices, established by decades of research, use a combination of disease-free seed, destruction or avoidance of inoculum sources, manipulation of sowing dates, seed and foliar fungicides, and cultivars with improved resistance. An understanding of the pathosystems and the inter-relationship between host, pathogen and the environment is essential to be able to make correct decisions for disease control without compromising the agronomic or economic ideal. For individual pathosystems, some components of the integrated management principles may need to be given greater consideration tha...
The goal of the book is to summarize those methods for evaluating Feynman integrals that have been developed over a span of more than fifty years. The book characterizes the most powerful methods and illustrates them with numerous examples starting from very simple ones and progressing to nontrivial examples. The book demonstrates how to choose adequate methods and combine evaluation methods in a non-trivial way. The most powerful methods are characterized and then illustrated through numerous examples. This is an updated textbook version of the previous book (Evaluating Feynman integrals, STMP 211) of the author. (orig.)
This study does not recommend the use of a pressure vessel made of integral stainless steel, due to the failure mechanisms under stress corrosion assisted by chlorides or polythionic acid. Are presented case studies of literature and analysis of reports of proceedings of RLAM reactors, showing that the materials produced by bimetallic clad overlay or are more appropriate, in terms of integrity, for use in equipment that the internal environment requires austenitic stainless steel specification.
Air conditioning is getting maximum attention these days when building new offices. For the new office building of the assurance company Zwitserleven in Amstelveen, Netherlands, a new air conditioning concept is applied, based on four existing air conditioning methods: cold storage and heat storage in the soil, heat pump, double-wall facades with integrated climate control, and radiant heating combined with a cooling system in the ceiling. The integration of these methods not only provides complete climate comfort for the employees but also substantial energy savings. 3 figs.,1 ills.
Despite the fact that there is still a potential to improve the efficiency of classic brown coal power plants by various measures from 36% to approx. 40-41%, the essential big increase of the efficiency to values of approx. 46% can only be achieved by the transition to combined gas and steam turbine processes. This means for the solid fuel brown coal that a gasification is integrated into the power plant process. The combined cycle power plant with an integrated brown coal gasification represents this new power plant concept. The applied High Temperature Winkler method (HTW-method) is an advancement of the Winkler method that works at environmental pressure and was developed by the Rheinbraun AG. (orig./DG).
System identification is a necessity in control theory. Classical control theory usually considers processes with integer order transfer functions. Real processes are usually of fractional order as opposed to the ideal integral order models. A simple and elegant scheme is presented for approximation of such a real world fractional order process by an ideal integral order model. A population of integral order process models is generated and updated by PSO technique, the fitness function being the sum of squared deviations from the set of observations obtained from the actual fractional order process. Results show that the proposed scheme offers a high degree of accuracy.
In this paper we present a new, accurate form of the heat balance integral method, termed the combined integral method (CIM). The application of this method to Stefan problems is discussed. For simple test cases the results are compared with exact and asymptotic limits. In particular, it is shown that the CIM is more accurate than the second order, large Stefan number, perturbation solution for a wide range of Stefan numbers. In the initial examples it is shown that the CIM reduces the standard problem, consisting of a PDE defined over a domain specified by an ODE, to the solution of one or two algebraic equations. The latter examples, where the boundary temperature varies with time, reduce to a set of three first order ODEs.
Cancer is a severe threat to human health. Early detection is considered the best way to increase the chance for survival. While the traditional cancer detection method, biopsy, is invasive, noninvasive optical diagnostic techniques are revolutionizing the way that cancer is diagnosed. Reflectance spectroscopy is one of these optical spectroscopy techniques showing promise as a diagnostic tool for pre-cancer detection. When a neoplasia occurs in tissue, morphologic and biochemical changes happen in the tissue, which in turn results in the change of optical properties and reflectance spectroscopy. Therefore, a pre-cancer can be detected by extracting optical properties from reflectance spectroscopy. This dissertation described the construction of a fiberoptic based reflectance system and the development of a series of modeling studies. This research is aimed at establishing an improved understanding of the optical properties of mucosal tissues by analyzing reflectance signals at ...
The goal of this work is to present a systematic Monte Carlo validation study on the clinical implementation of the enhanced dynamic wedges (EDWs) into the Pinnacle"3 (Philips Medical Systems, Fitchburg, WI) treatment planning system (TPS) and QA procedures for patient plan verification treated with EDWs. Modeling of EDW beams in the Pinnacle"3 TPS, which employs a collapsed-cone convolution superposition (CCCS) dose model, was based on a combination of measured open-beam data and the 'Golden Segmented Treatment Table' (GSTT) provided by Varian for each photon beam energy. To validate EDW models, dose profiles of 6 and 10 MV photon beams from a Clinac 2100 C/D were measured in virtual water at depths from near-surface to 30 cm for a wide range of field sizes and wedge angles using the Profiler 2 (Sun Nuclear Corporation, Melbourne, FL) diode array system. The EDW output factors (EDWOFs) for square fields from 4 to 20 cm wide were measured in virtual water using a ...
integrates multiple fisheries information systems to obtain a view of U.S. freshwater fish distribution. Information about NFDI National Framework for Early Detection, Rapid...
Advances in molecular technologies challenge the different concepts of causality in biology, epidemiology and multistage mathematical models. The lack of integration of the different aspects of causality...Full Text Available
... The acousto-optic device have a 30 MHz 1 ... coherent systems including compact non-coherent optical ... a relatively simple phase switching approach. ...
... communication with the divisions of the Seventh Army ... the 1st Infantry Division could not ... commanding general of the 82nd Airborne Division, landed ...
the southern Arabian Gulf region left its signature on the heterogeneous aerosol .... Arabian Gulf region, since large differences in ? may be caused by ...
As organizations merge or collaborate closely, an important question is how their existing software assets should be handled. If these previously separate organizations are in the same business domain - they might even have been competitors - it is likely that they have developed similar software systems. To rationalize, these existing software assets should be integrated, in the sense that similar features should be implemented only once. The integration can be achieved in different ways. Success of it involves properly managing challenges such as making as well founded decisions as early as possible, maintaining commitment within the organization, managing the complexities of distributed teams, and synchronizing the integration efforts with concurrent evolution of the existing systems.Th...
Described as the body's largest organ, the skin is strategically located at the interface with the external environment where it has evolved to detect, integrate and respond to a diverse range...Full Text Available
A device made of amorphous silicon which exhibits inductive properties at certain voltage biases and in certain frequency ranges in described. Devices of the type described can be made in integrated circuit form.
The recent developments in the mobile technology (mobile phones, middleware) created a need for new methods of protecting the code transmitted through the network. The proposed mechanisms not only secure the compiled program, but also the data, that can be gathered during its "journey". The oldest and the simplest methods are more concentrated on integrity of the code itself and on the detection of unauthorized manipulation. Other, more advanced proposals protect not only the code but also the execution state and the collected data. The paper is divided into two parts. The first one is mostly devoted to different methods of securing the code and protecting its integrity; starting from watermarking and fingerprinting, up to methods designed specially for mobile agent systems: encrypted function, cryptographic traces, time limited black-box security, chained-MAC protocol, publicly-verifiable chained digital signatures The second part presents new ...
screens of fine wire was placed between the gun barrel and the panel to determine whether the gun performance was repeatable. The transducers were placed ...
BackgroundMost reports of scientific misconduct have been focused on principal investigators and other scientists (e.g., biostatisticians) involved in the research...Full Text Available
A quantum mechanical analysis of the guided light in integrated photonics waveguides is presented. The analysis is made starting from one-dimensional (1D) guided vector modes by taking into account the modal orthonormalization property on a cross section of an optical waveguide, the vector structure of the guided optical modes and the reversal-time symmetry in order to quantize the 1D vector modes and to derive the quantum momentum operator and the Heisenberg equations. The results provide a quantum-consistent formulation of the linear and nonlinear quantum light propagations as a function of forward and backward creation and annihilation operators in integrated photonics. As an illustration, an application to an integrated nonlinear directional coupler is given, that is, both the nonlinear momentum and the Heisenberg equations of the nonlinear coupler are derived.
This article argues that anti-corruption agencies at the local-level have been successful in a way that can be evaluated and emulated. A related contention is that corruption control is most effective when the central public integrity agency is part of both a local anticorruption network and a local public management network. Quite reasonably, the international anti-corruption project has focused most time and energy on advocating and assessing efforts made to ensure public integrity at the national level. Baseline studies by scholars and supra-national integrity nongovernmental organizations (NGOs) identify the form of corruption control (if any) adopted by the central government. Key considerations in assessing the status of national anti-corruption agencies (ACAs) are the ones mentioned...
This paper presents a set of heuristics useful for optimal design of heat-exchange networks and integrated heat and power systems. These heuristics consolidate other approaches reported in the literature in a formalism useful for design purposes.
Jul 7, 2008 ... HOUSTON -- Following a detailed, integrated assessment, NASA selected ... in the center that provides a 360-degree view around the station. ...
Global supply chain practices and their effects have received considerable attention over the last two decades. In the recent past, the need for integration across supply chains has been identified as a key for effective and efficient operations of supply chains. This is observed with the increasing trend of collaborative partnerships among supply chain partners. This paper presents an integrated approach for manufacturing and distribution networks within the supply chain system of a global car company. The paper shows that the integration of manufacturing and distribution networks creates the environment for effective planning of many components and execution/follow-up of those plans. These components include materials, resources, operations/activities, suppliers and customers. The main f...
... for team sizes ranging from 2 to 6 robots. ... intercep- tion mission, a robot must overtake a ... of unmanned underwater vehicles (UUV), unmanned aerial ...
Hexavalent chromium in corrosion preventive coating 64 ...To prevent corrosion of rims, weights can be provided with a coating. ...integrated clips brought up technical problems due to corrosion
Genetics, Genomics, and Molecular Biology USGS scientists develop and integrate new genetic and molecular techniques into systematic analyses to describe individuals and populations of fish .....
We used neuropsychological tasks to investigate integrity of brain circuits linking orbitofrontal cortex and amygdala (orbitofrontal-amygdala), and dorsolateral prefrontal cortex and hippocampus...Full Text Available
This Contractor Performance Plan is a detailed execution plan covering fiscal years 2007 and 2008 that implements the Integration Project Baseline in a more aggressive cost and schedule approach.
BackgroundSenescence is integral to the flowering plant life-cycle. Senescence-like processes occur also in non-angiosperm land plants, algae and photosynthetic prokaryotes. Increasing...Full Text Available
During the last years, large-scale simulations of realistic physical environments which support the interaction of multiple participants over the Internet have become increasingly available and economically significant, most notably in the computer gaming industry. Such systems, commonly called networked virtual environments (NVEs), are usually based on a client-server architecture where for performance reasons and bandwidth restrictions, the simulation is partially deferred to the clients. This inevitable architectural choice renders the simulation vulnerable to attacks against the semantic integrity of the simulation: malicious clients may attempt to compromise the physical and logical laws governing the simulation, or to alter the causality of events a posteriori. In this paper, we initiate the systematic study of semantic integrity in NVEs from a security point of view. We argue that naive policies to enforce semantic ...
Emergency Core Cooling System (ECCS) is a safety critical system provided to mitigate the consequence of Loss of Coolant Accident (LOCA) in PHWR. Unlike 220MWe, all header injection has been introduced in 540MWe to simplify the logic. ECCS Integrated Test is schematic approach to establish that ECC system will behave as per design intent during actual LOCA condition. Objective of ECCS Integrated test is to ascertain that various ECC system components operate as intended in design. Additionally, the various system resistances which form the input to LOCA analysis are validated. This test has been carried out by creating actual LOCA during cold and pressurised condition of PHT system to establish all phases of injection with overlap. This paper discusses the results obtained during the Integrated Test and comparison with the prediction during the commissioning of first unit of 540 MWe. (author)
Advances in modern neuroscience require the identification of principles that connect different levels of experimental analysis, from molecular mechanisms to explanations of cellular functions,...Full Text Available
Integration of Eddy Current probes with MAUS scanning system has been demonstrated to produce scan rates of 36 sqft/hour. System will be designed to ...
The objectives of the Phase-I research address the development of mechanics- based structural integrity analysis methodologies for composite sandwich panels ...
assessment.Both frustumswerein excellentcondition.No missingTPSor debonds/ unbondswere ...... 2) The average time for a 360 degree rotation was 29.9 seconds. ...
The use of repeated expressions to establish coreference allows an investigation of the relationship between basic processes of word recognition and higher-level...Full Text Available
The pleiotropic actions of neuromodulators on pre- and postsynaptic targets present challenges to disentangling the mechanisms underlying regulation of synaptic transmission. Within the striatum,...Full Text Available
Purpose - The article discusses how care programmes and integrated care pathways can be linked, finding ways to improve healthcare process professional and logistical quality from a supply chain and a network point-of-view. Design/methodology/approach - The authors argue that owing to cost containment goals and increasing healthcare demand, healthcare services systems are challenged to improve service quality, whilst at the same time finding ways to improve delivery processes. It explores if the combination of two instruments, care programmes and integrated care pathways, can meet both goals. This combination is illustrated by an example from the Institute of Mental Health Care Eindhoven en de Kempen. Findings - Analysis suggests that care programmes can be combined with integrated care pa...
... being designed for projected Minotaur IV launch load environments, with design objectives of light weight, integrated vibration isolation, low shock ...
Bone integrity and mineral status were studied with a noninvasive method in uremic patients with severe secondary hyperparathyroidism undergoing maintenance hemodialysis.
We identify new organelles associated with the vacuolar system in plant cells. These organelles are defined biochemically by their internal content of three integral membrane proteins: a chimeric reporter...Full Text Available
... Development for Envir Activities Materials Use: Sci, Eng, & Society Synopsis of Program: This ... for Environmental Activities (IDEA); Materials Use: Science, Engineering, & Society (MUSES). In all ...
... reliability of fiber optics in preference over space-based systems. ... STC web site at www.aiaa.org/tc ... and improve lethalities of US weapon systems 3 ...
ObjectiveTo describe the development and assessment of monographs as an assignment to incorporate evidence-based medicine (EBM) and pharmacoeconomic principles into a third-year...Full Text Available
flow separation and vortices shed from the airplane forebody when flying at ...... AGARD Conference Proceedings 548: Technologies for Highly Maneuverable ...
... HAL and VAL limits with the specified PMD of 10-7, the GPS range-rate error must be monitored to assure that the inertial navigation error is not ...
A numerical and experimental investigation is carried out in a solar thermochemical reactor for the thermal dissociation of ZnO at 2000 K using concentrated solar energy. The reactor consists of a cavity-receiver lined with ZnO particles and directly exposed to high-flux irradiation. A transient heat transfer model is formulated to link the rate of radiation, convection, and conduction heat transfer to the reaction kinetics. The radiosity and Monte Carlo methods are applied to obtain the distribution of net radiative fluxes at the internal surfaces of the reactor cavity and at the surface of the ZnO bed. Validation is accomplished in terms of the calculated and measured transient temperature profiles and chemical reaction rates.
Solar cavity receiver plays a dominant role in the light-heat conversion. Its performance can directly affect the efficiency of the whole power generation system. A combined calculation method for evaluating the thermal performance of the solar cavity receiver is raised in this paper. This method couples the Monte-Carlo method, the correlations of the flow boiling heat transfer, and the calculation of air flow field. And this method can ultimately figure out the surface heat flux inside the cavity, the wall temperature of the boiling tubes, and the heat loss of the solar receiver with an iterative solution. With this method, the thermal performance of a solar cavity receiver, a saturated steam receiver, is simulated under different wind environments. The highest wall temperature of the boi...
A large database (507 station-years) of daily suspended particulate matter (SPM) concentration and discharge data from 36 stations on river basins ranging from 600?km2 to 600,000?km2 in size (USA and Europe) was collected to assess the effects of SPM transport regime on bias and imprecision of flux estimates when using infrequent surveys and the discharge-weighted mean concentration method. By extracting individual SPM concentrations and corresponding discharge values from the database, sampling frequencies from 12 to 200 per year were simulated using Monte Carlo techniques. The resulting estimates of yearly SPM fluxes were compared to reference fluxes derived from the complete database. For each station and given frequency, bias was measured by the median of relative errors between estima...
We present results from a kinetic Monte Carlo simulation of boron transient enhanced diffusion (TED) in silicon. Our approach avoids the use of phenomenological fits to experimental data by using a complete and self-consistent set of values for defect and dopant energetics derived mostly from {ital ab initio} calculations. The results predict that, during annealing of 40 keV B-implanted Si at 800{degree}C, there exists a time window during which all the implanted boron atoms are substitutional. At earlier or later times, the interactions between free silicon self-interstitials and boron atoms drive the growth of boron clusters and result in an inactive boron fraction. The results show that the majority of boron TED takes place during the growth period of interstitial clusters and not during their dissolution. {copyright} {ital 1998 American Institute of Physics.}
We present results from a kinetic Monte Carlo simulation of boron transient enhanced diffusion (TED) in silicon. Our approach avoids the use of phenomenological fits to experimental data by using a complete and self-consistent set of values for defect and dopant energetics derived mostly from ab initio calculations. The results predict that, during annealing of 40 keV B-implanted Si at 800 degree C, there exists a time window during which all the implanted boron atoms are substitutional. At earlier or later times, the interactions between free silicon self-interstitials and boron atoms drive the growth of boron clusters and result in an inactive boron fraction. The results show that the majority of boron TED takes place during the growth period of interstitial clusters and not during their dissolution. copyright 1998 American Institute of Physics.
We study the finite temperature electroweak transition with non-perturbative lattice Monte Carlo simulations. We find that it is of first order, at least for Higgs masses up to 80 GeV. The critical temperature of the phase transition is found to be smaller than that determined by a 2-loop renormalization group improved effective potential. The jump of the order parameter at the critical temperature is considerably larger than the perturbative value. By comparing lattice data and perturbation theory, we demonstrate that the latter, for the computation of the vacuum expectation value of the Higgs field v(T) in the broken phase at given temperature, converges quite well, provided v(T)/T>1. An upper bound on the Higgs mass necessary for electroweak baryogenesis in the light of the lattice data is briefly discussed.
The effects of pore curvature and surface heterogeneity on the adsorption of water on a graphitic surface at 298 K were investigated using a Grand Canonical Monte Carlo (GCMC) simulation. Slit and cylindrical pores are used to study the curvature effects. To investigate the surface heterogeneity the functional group and the structural defect on the surface were specifically considered. The hydroxyl group (OH) is used as a model for the functional group and the water potential model proposed by Muller et al. is used to calculate the water interaction. For the homogeneous cylinder, the pore filling occurs at a pressure lower than the saturation pressure of the water model, while it is greater in the case of homogeneous slit pore. The size of hysteresis loop is more sensitive to the length of...
We point out that the occurrence of a large-distance ''Coulomb'' term in the static quark-antiquark potential is related to tachyons in the underlying string model. Thus, the tachyon contents determine the coefficient of the ''Coulomb'' term. For a scalar string, one has Luescher's result -(d-2) /24R, for a Neveu-Schwarz string one obtains -(d-2) /16R, whereas for the Ramond string (and strings without tachyons) the ''Coulomb'' term vanishes. Monte Carlo data for QCD prefers a ''Coulomb'' term. This is not inconsistent with the tachyonic origin of this term, since the large-distance behaviour of string models is consistent even with tachyons. Also, critical dimensions are not necessary at large distances.
We study spectroscopically determined iron abundances of 642 solar-type stars to search for the signature of accreted iron-rich material. We find that the metallicity [Fe/H] of a subset of 466 main sequence stars, when plotted as a function of stellar mass, mimics the pattern seen in lithium abundances in open clusters. Using Monte Carlo models we find that, on average, these stars have accreted about 0.4 Earth masses of iron while on the main sequence. A much smaller sample of 19 stars in the Hertzsprung gap, which are slightly evolved and whose convection zones are significantly more massive, have lower average [Fe/H], and their metallicity shows no clear variation with stellar mass. These findings suggest that terrestrial-type material is common around solar type stars.
We have studied, by a Monte Carlo simulation code ACAT-DIFFUSE, the fluence-dependence of the amount of retained helium atoms in tungsten carbide at room temperature under helium ion bombardment. The retention behavior may be understood qualitatively in terms of irradiation-dependent diffusion coefficient assumed and range. The emission processes from tungsten carbide under helium ion irradiation derived were compared with each other. We have discussed the retention curves for incident energy of 5 keV at incident angles of 0deg and 80deg and of 500 eV at 0deg. The energy spectra of helium atoms reflected from tungsten carbide for incident energy of 500 eV at 0deg and 80deg were compared with those from graphite and tungsten. (author)
The difficulty and unstability of the parameter estimation of mass action model for radioimmunoassay are the major obstacles to its application in routine work. A unique robust curve fitting method for the mass action model was introduced. large numbers of simulated data were generated with Monte Carlo method to compare the accuracy and robustness of the robust fitting and the least squre fitting (Marquardt method) of the mass action model. The results demonstrated that the new robusst method significantly increased the accuracy and reliability of data processing with the mass action model. It also markedly reduced the influence of outliers in the fitting of RIA data, and it was always 'convergent'. This robust method could also be applied to other mathematical models which can be transformed into high degree equations and no initial estimates are needed for curve fitting.
In this report we present an update on the results from the Amanda-B10 detector which operated in 1997 at depths of 1500 to 2000 meters in the deep Antarctic ice. The goal of Amanda project is to search for extra-terrestrial neutrinos. As a precursor to such a search we have studied atmospheric neutrinos which act as a calibration source for the detector. The observation of atmospheric neutrinos at a rate consistent with Monte-Carlo predictions establishes Amanda-B10 as a neutrino telescope. The Amanda-B10 data has been searched for evidence of several classes of neutrinos and for magnetic monopoles. Searches for a diffuse high energy neutrino flux and for neutrinos in coincidence with gamma-ray bursts have been conducted. Preliminary data analyses show no excess of neutrinos has been found. (A.C.)
New values of the factors required to convert the reading of a radiation dosemeter calibrated in terms of air kerma (or exposure) into absorbed dose to water for medium-energy x-radiation were given in a code of practice published by the International Atomic Energy Agency in 1987. These are not considered to possess sufficient support from other sources. It is therefore recommended that the F-factors given or codes of practice should continue to be used. Values of backscatter factors for low-energy x-radiation (below 140 kV) in Supplement 17 of the British Journal of Radiology appear to be inaccurate. New values based on Monte Carlo calculations, and supported by new experimental data, are given for use in radiotherapy. (author).
New values of the factors required to convert the reading of a radiation dosemeter calibrated in terms of air kerma (or exposure) into absorbed dose to water for medium-energy x-radiation were given in a code of practice published by the International Atomic Energy Agency in 1987. These are not considered to possess sufficient support from other sources. It is therefore recommended that the F-factors given or codes of practice should continue to be used. Values of backscatter factors for low-energy x-radiation (below 140 kV) in Supplement 17 of the British Journal of Radiology appear to be inaccurate. New values based on Monte Carlo calculations, and supported by new experimental data, are given for use in radiotherapy. (author).
The present authors have previously developed a new method for Probabilistic Fracture Mechanics (PFM), which they call Recursive Distribution (RD) method. The method is based on the construction of the Lebesgue-Stieltjes measure through a deterministic mapping defining a crack growth process. In the present paper, its theoretical background is first discussed, and the Lebesgue decomposition of the measure is given. Then a numerical example of a Light Water Reactor (LWR)`s piping problem is solved by the present method, and the results are compared with those of the Monte Carlo (MC) method. In addition to leakage probability, a variation in stress cycles of the marginal distribution of an aspect ratio of a semi-elliptical surface crack is calculated, which will be used in a study on LBB evaluation.
An energy transfer and conversion model for high-temperature solar cavity receivers has been developed using the transport behaviour of solar radiation as described by the spectral radiative exchange factors. A Monte-Carlo ray-tracing method coupled with optical properties was adopted, to predict radiation characteristics of the solar collector system by calculating radiative exchange factors. A cavity receiver with a plano-convexo quartz window was proposed, based upon the directional characteristics of the focal flux and the redistribution effect of the quartz window. Parametric studies on the windowed receiver provided a more uniform flux distribution, higher efficiency and lower loss than the windowless receivers. The predicted results serve as a design reference for the solar receiver...
The Front End Enclosure (FEE) of the Linac Coherent Light Source (LCLS) is a shielding housing located between the electron dump area and the first experimental hutch. The upstream part of the FEE hosts the commissioning diagnostics for the FEL beam. In the downstream part of the FEE, two sets of grazing incidence mirror and several collimators are used to direct the beam to one of the experimental stations and reduce the bremsstrahlung background and the hard component of the spontaneous radiation spectrum. This paper addresses the beam loss assumptions and radiation sources entering the FEE used for the design of the FEE shielding using the Monte-Carlo code FLUKA. The beam containment system prevents abnormal levels of radiations inside the FEE and ensures that the beam remains in its intended path is also described.
The general objective of this workshop is to investigate and discuss methods by which uncertainties in mass balance models for toxics in the Great Lakes may be reduced. As described by the workshop prospectus, this paper is focused on problems of reducing (and quantifying) uncertainty as they relate to in situ field observations/system response measurements for the establishment of initial conditions, boundary conditions, calibration/confirmation data sets, and model post-audit data sets.'' I have taken this description to refer not only to the evaluation of uncertainty in the field observations themselves, but also to the uncertainty associated the analyses of in situ observations as they interact in the overall modeling process. Thus, I will be concerned here with quantification and reduction of uncertainty both (1) as they may be applied to descriptions of the system that is being modeled and (2) as they may be associated with model simulations.
We report on the implementation of an interface between the SANC generator framework for Drell-Yan hard processes, which includes next-to-leading order electroweak (NLO EW) corrections, and the Herwig++ and Pythia8 QCD parton shower Monte Carlos. A special aspect of this implementation is that the initial-state shower evolution in both shower generators has been augmented to handle the case of an incoming photon-in-a-proton, diagrams for which appear at the NLO EW level. The difference between shower algorithms leads to residual differences in the relative corrections of 2-3% in the p_T(mu) distributions at p_T(mu)>~50 GeV (where the NLO EW correction itself is of order 10%).
Monte Carlo calculations of proton irradiation of permanent magnets for the TRADE experiment have been performed. An irradiation dose of about 4´106 Gy/yr/mA has been estimated due to beam losses in normal operating conditions. Existing experimental results indicate that this irradiation level may induce a considerable demagnetization: in fact, a dose of 6´107 Gy induces a remanence loss of 0.3 % on samples of Sm2Co17 magnets, which are the most resistant type. More detailed calculations with the final design of the magnets and of the beam line are suggested, to determine if the irradiation levels allowed a reliable operation of the permanent magnets for the entire duration of the TRADE experiment. Damage and gas production rates have also been calculated; the values obtained are very low, thus confirming that the demagnetization process is in great part reversible.
Drug particle size distribution has a profound impact to the content uniformity in low-dose solid drug products. We derived theoretically the skewness of potency distribution as a function of particle size distribution and target dose. It was demonstrated that both skewness and coefficient of variation diverge simultaneously with inverse square root of the target dose. This scaling relation was observed in recent experiment and was verified by Monte Carlo (MC) simulation, which was employed for the first time to solve for the full potency distribution from a random retrieving model. When tested against the criteria from USP uniformity of dosage units, MC simulation showed a striking anisotropic distribution of the data. This suggests a full-scale consideration of the potency distribu...
Analyses of pressure vessel damage resulting from neutron irradiation have primarily relied on two-dimensional transport calculations and a spatial-synthesis methodology to accommodate three-dimensional effects in the results of two two-dimensional calculations. In this paper, the authors report on calculations made on the Pool Critical Assembly (PCA) Benchmark, Configuration 12/13, using the three-dimensional, continuous energy Monte Carlo transport code, MCNP, and the three-dimensional, multigroup, diffusion accelerated discrete ordinates transport code THREEDANT. Neutron fluxes and activation rates as determined from these two calculations are compared to each other and to experimental results in the literature. The authors also draw some conclusions on the value of 3D calculations on the interpretation of experimental results.
To establish the phylogeographic relationships in rabies viruses in Brazil, we studied a dataset retrieved from GenBank consisting of 71 genetic sequences from the coding region of the N gene of rabies viruses isolated in dogs over a period of 22?years. The Bayesian Markov chain Monte Carlo method available in the BEAST package was used with the GTR+G+?4 evolutionary model in conjunction with the relaxed uncorrelated lognormal molecular clock model and an exponential growth tree prior. A discrete phylogeographic diffusion model was also analyzed using a standard continuous-time Markov chain viewed with Google Earth to provide a spatial projection of the diffusion of genetic lineages based on their phylogeographic relationships. The topology of the time and substitution phylogenetic trees a...
The basic relationships of the convolution/superposition dose calculation technique are reviewed, and a modelling technique that can be used for obtaining a satisfactory beam model for a commercially available convolution/superposition-based treatment planning system is described. A fluence energy spectrum for a Co-60 treatment machine obtained from a Monte Carlo simulation was used for modelling the fluence spectrum for a Rokus M machine. In order to achieve this model we measured the depth dose distribution and the dose profiles with a Wellhofer dosimetry system. The primary fluence was iteratively modelled by comparing the computed depth dose curves and beam profiles with the depth dose curves and crossbeam profiles measured in a water phantom. The objective of beam modelling is to build a model of the primary fluence that the patient is exposed to, which can then be used for the calculation of the dose deposited in the patient. (authors)
Using a Markov Chain Monte Carlo optimization algorithm and a computer simulation, I find the passenger ordering which minimizes the time required to board the passengers onto an airplane. The model that I employ assumes that the time that a passenger requires to load his or her luggage is the dominant contribution to the time needed to completely fill the aircraft. The optimal boarding strategy may reduce the time required to board and airplane by over a factor of four and possibly more depending upon the dimensions of the aircraft. I explore some features of the optimal boarding method and discuss practical modifications to the optimal. Finally, I mention some of the benefits that could come from implementing an improved passenger boarding scheme.
The possibility of radio galaxies being random sample of otherwise normal elliptical galaxies is tested. Starting with the observed optical luminosity functions for elliptical galaxies, it is shown that the probability of an elliptical forming a radio source is a continuous, increasing function of optical luminosity, precisely proportional to square of the optical luminosity of the galaxy. Once the probability function is fixed, the luminosity function of normal elliptical galaxies is used as input for Monte Carlo simulations that reproduce the distribution of radio galaxies in the radio-optical luminosity plane. Our results show that radio galaxies are luminosity biased, but otherwise random sample of elliptical galaxies. This unified view of radio and non-radio ellipticals also explains the well known difference of 0.5 mag in average optical luminosity between FRI and FRII radio galaxies as a simple selection effect. Specifically, FRII appear preferentially in ...
This document reproduces the final project of David San Fabian Ayuso, presented on May 26, 2005, for the obtention of the engineer degree of the Carlos III University of Madrid. A single-phase, isothermal model, including both electron and proton transport, is introduced for the simulation of polymer electrolyte fuel cells (PEM). The model is implemented in the commercial code Fluent 6.0, through the use of UDFs (User Defined Functions). In order to validate the model, a single canal of a PEM monocell is simulated in three dimensions. The obtained result are qualitatively satisfactory. It is observed that it is not essential to solve the current collectors when a monocell is considered (and not a stack). in the present study, the number of nodes is the computational grid appears to be too low in the membrane zone in order to make a complete validation of the model. (Author) 20 refs.
A thermal neutron imaging facility for real-time neutron radiography and computed tomography has recently been developed at the University of Texas TRIGA reactor. Extensive Monte Carlo design calculations were used to determine optimal design parameters of the neutron collimator system to avoid costly trial and error. Thermal neutron flux determined by gold foil activation is 5 {times} 10{sup 6} n/cm{sup 2}{center_dot}s at the primary imaging location with beam size of 22.5 cm in diameter. The collimation ratio can be varied from 125 to 235. The neutron-to-gamma ratio is 7.8 {times} 10{sup 6} n/cm{sup 2}{center_dot}mR. The facility has been tested for radiography and tomography applications and is now fully operational.
The present work summarizes efforts on the simulation of natural convection cooling within the IFMIF target and test cell. The simulations have been performed with the STAR-CD code using the k-#omega# high-Reynolds number turbulence model. A dedicated thermohydraulic model has been devised including Lithium loop components. Nuclear heat production has been calculated by the Monte-Carlo code McDeLicious for different parts of the target and test cell walls and was used as input for the STAR-CD simulations. Helium atmospheres at several pressures from 0.1 to 10"-"5 MPa have been investigated. In order to limit the maximum temperature of the concrete walls to 80 deg. C it was necessary to add thermal insulation layers to the hot Lithium loop surfaces and a conceptual system of two cooling layers in different depths of the concrete walls.
This paper aims to provide a comprehensive review of the various models or simulations for predicting drug release from bulk-degrading systems. A brief description of bulk degradation processes and factors affecting the degradation rate, and consequently the release kinetics, is presented first. Next, several important classical models, often used as the basis for subsequent model development, are discussed. Both mathematical models and Monte-Carlo based simulations have been developed for controlled release from bulk-degrading systems. The mathematical models can be further subdivided into two categories. First, the diffusion-based models whose transport mechanism is mainly governed by diffusion, but with degradation-dependent diffusion coefficients. These are generally simpler and easier...
A model for risk and reliability analysis of complex multifunctional production process systems is presented. The model employs Monte-Carlo and Markov Chain algorithms that uses a weighted index to train and simulate the fuzzy hazard data sets which represents failure outcomes of risk component transient and non-transient systems. Early simulation results shows that hazard rates and the risk of containment loss from typical floating production and storage offloading (FPSO)-Riser system for the risk components in parallel or series increases exponentially with time and decreases as safety ratings fraction increases. The reliability value decreases with time and safety fraction (SFAC) for all fuzzy hazard classifications. The results of the computed mean time before repair (MTBR) show that t...
By using a variational Monte Carlo method, we examine an effective low-energy model for LaFeAsO derived from an ab initio downfolding scheme. We show that quantum and many-body fluctuations near the antiferromagnetic (AF) quantum critical point largely reduce the antiferromagnetic ordered moment. Our derived model not only quantitatively reproduces the small ordered moment in LaFeAsO, but also accounts for the diversity from LaFePO, BaFe_2As_2 to FeTe. Electron correlation is found to determine the observed material dependence. We also find that LaFeAsO is subject to large orbital fluctuations, sandwiched by the AF Mott insulator and weakly correlated metals. The orbital fluctuations and Dirac-cone dispersion hold keys for the diverse magnetic properties. (author)
A successful boron neutron capture treatment (BNCT) of a patient with multiple liver metastases has been first given in Italy, by placing the removed organ into the thermal neutron column of the Triga research reactor of the University of Pavia. In Finland, FiR 1 Triga reactor with an epithermal neutron beam well suited for BNCT has been extensively used to irradiate patients with brain tumors such as glioblastoma and recently also head and neck tumors. In this work we have studied by MCNP Monte Carlo simulations, whether it would be beneficial to treat an isolated liver with epithermal neutrons instead of thermal ones. The results show, that the epithermal field penetrates deeper into the liver and creates a build-up distribution of the boron dose. Our results strongly encourage further studying of irradiation arrangement of an isolated liver with epithermal neutron fields.
The confining effect of a spherical substrate inducing anchoring (normal to the surface) of rod-like liquid crystal molecules contained in a thin film spread over it has been investigated with regard to possible changes in the nature of the isotropic-to-nematic phase transition as the sample is cooled. The focus of these Monte Carlo simulations is to study the competing effects of the homeotropic anchoring due to the surface inducing orientational ordering in the radial direction and the inherent uniaxial order promoted by the intermolecular interactions. By adopting entropic sampling procedure, we could investigate this transition with a high temperature precision, and we studied the effect of the surface anchoring strength on the phase diagram for a specifically chosen geometry. We find that there is a threshold anchoring strength of the surface below which uniaxial nematic phase results, and above which the isotropic fluid cools to a radially ordered nematic ...
We consider low-energy nucleons at next-to-next-to-leading order in lattice chiral effective field theory. Three-body interactions first appear at this order, and we discuss several methods for determining three-body interaction coefficients on the lattice. We compute the energy of the triton and low-energy neutron-deuteron scattering phase shifts in the spin-doublet and spin-quartet channels using Luescher's finite-volume method. In the four-nucleon system we calculate the energy of the {alpha}-particle using auxiliary fields and projection Monte Carlo. (orig.)
Data from Monte Carlo transport codes are used to model radiobiological effects. We previously reported the Fourier analysis of ionization data generated by simulating a 500-keV proton traversing water. Here, we extend Fourier analysis to energy transfer data of another radiation type, a 14-MeV #alpha#-particle. A radiobiological model based on this frequency-domain analysis views cell as an information processing system . It lends itself naturally to traditional engineering analyses. One engineering principle-the output response of a linear system to random signal-is applied here to explain the fact that there is measurable difference in the magnitude of the biological effectiveness when a given biological system is irradiated with two different radiation types of the same Linear Energy Transfer (LET).
The Monte-Carlo transport code MCNP was used to evaluate possible arrangements of cores for the MIT Reactor using monolithic LEU fuel. Plate and moderator thicknesses were varied, and fixed absorbers and inner reflectors added in an effort to maximize available neutron fluxes at in-core and ex-core locations of experimental facilities. Addition of D_2O in the H_2O moderator was also evaluated. Comparisons of the fast, epithermal, and thermal fluxes were made at selected locations. Keff was also evaluated and critical blade heights compared with the existing HEU core. Results indicate that the LEU fluxes could approach HEU values with the use of a fueled in-core experimental facility, a fixed boron absorber spider and an inner beryllium reflector. (author)
A major deficiency of current photon calculation methods that are based on the concept of primary and scatter separation is their inability to handle the condition of electronic disequilibrium. This deficiency is examined and it is shown that the limitation is not inherent in the algorithms themselves but is, at least in part, in the data which the algorithms use. A new concept of primary and scatter separation is developed to cover the condition of electronic disequilibrium. This new concept requires little change to the existing algorithms and only additional data are required, which are generated using Monte Carlo calculation methods. The new concept is tested using programs in the Theratronics Theraplan treatment-planning system, and two calculation examples illustrate the ability to model electron transport and also the improvement over the existing algorithms. Close analogy of the extended concept with the convolution/superposition method of dose calculation ...
Purpose - Burr distribution has been proved to be a useful failure model. It can assume different shapes which allow it to be a good fit for various lifetimes data. Hybrid censoring is an important way of generating lifetimes data. The purpose of this paper is to estimate an unknown parameter of the Burr type XII distribution when data are hybrid censored. Design/methodology/approach - The problem is dealt with through both the classical and Bayesian point of view. Specifically, the methods of estimation used to tackle the problem are maximum likelihood estimation method and Bayesian method. Empirical Bayesian approach is also considered. The performance of all estimates is compared through their mean square error values. The paper employs Monte Carlo simulation to evaluate the mean square...
Incomplete convergence in numerical simulation such as computational physics simulations and/or Monte Carlo simulations can enter into the calculation of the objective function in an optimization problem, producing noise, bias, and topo- graphical inaccuracy in the objective function. These affect accuracy and convergence rate in the optimization problem. This paper is concerned with global searching of a diverse parameter space, graduating to accelerated local convergence to a (hopefully) global optimum, in a framework that acknowledges convergence uncertainty and manages model resolu- tion to efficiently reduce uncertainty in the final optimum. In its own right, the global-to-local optimization engine employed here (devised for noise tolerance) performs better than other classical and contemporary optimization approaches tried individually and in combination on the "industrial" test problem to be presented.
In this paper, the effects of surface mediation on the adsorption isotherm and isosteric heat of adsorption on a graphite surface were investigated, as the surface mediation is known to affect the intermolecular interaction of adsorbed molecules close to the surface. Kim and Steele (Phys. Rev. B 45 (11) (1992) 6226-6233) and others have assumed that the surface mediation is confined only to the first layer. This will be tested in this paper with a combined experimental and Grand Canonical Monte Carlo (GCMC) simulation of adsorption of argon on graphitized thermal carbon black (GTCB) over a range of temperatures (77-95.25K). By matching the simulation results against the experimental data, we have found that the surface mediation is extended up to the fourth layer, rather than only the firs...
The Advanced Photon Source injector synchrotron is a 7-GeV positron machine with a standard alternating gradient lattice. The calculated effect of dipole magnet strength errors on the orbit distortion, simulated by Monte Carlo, was reduced by sorting pairs of magnets having the closest simulated measured strengths to reduce the driving the term of the integer resonance nearest the operating point. This method resulted in a factor of four average reduction in the rms orbit distortion when all 68 magnets were sorted at once. The simulated effect of magnet measurement experimental resolution was found to limit the actual improvement. The {Beta}-beat factors were similarly reduced by sorting the quadrupole magnets according to their gradients.
Dissociative electron capture to CCl_4 is studied by measuring the angular and velocity distribution of Cl"- ions produced in collisions with velocity selected K(np) Rydberg atoms. Analysis of the data using a Monte Carlo collision code that models the detailed kinematics of the reaction indicates that the lifetime of the CCl_4"- intermediate formed by Rydberg electron capture is 7.5 #+-# 2.5 ps and that, upon dissociation, only a small fraction of the excess energy of reaction appears as translational energy of the Cl"- and CCl_3 fragments. The present approach is one of the few experimental techniques yet devised that can probe the lifetimes of collisionally-produced excited states on a picosecond timescale, and demonstrates that Rydberg atoms provide a unique tool with which to investigate the dynamics of dissociative electron attachment. Measurements are being extended to additional species, including CFCl_3.
We investigate the behavior of a circuit QED device when the resonator is initially populated with a mesoscopic coherent field. The strong coupling between the cavity and the qubit produces an entangled state involving mesoscopic quasi-pointer states with respect to cavity dissipation. The overlap of the associated field components results in collapse and revivals for the Rabi oscillation. Although qubit relaxation and dephasing do not preserve these states, a simple analytical description of the dissipative dynamics of the circuit QED device including cavity relaxation as well as qubit dissipation is obtained from the Monte-Carlo approach. Explicit predictions for the spontaneous and induced Rabi oscillation signals are derived and sucessfully compared with exact calculations. We show that these interesting effects could be observed with a 10 photon field in forthcoming circuit QED experiments.
An essential part of fire risk assessment is the analysis of fire hazards and fire propagation. In this work, models and tools for two different aspects of numerical fire simulation have been developed. The primary objectives have been firstly to investigate the possibility of exploiting state-of-the-art fire models within probabilistic fire risk assessments and secondly to develop a computationally efficient solver of thermal radiation for the Fire Dynamics Simulator (FDS) code. In the first part of the work, an engineering tool for probabilistic fire risk assessment has been developed. The tool can be used to perform Monte Carlo simulations of fires and is called the Probabilistic Fire Simulator (PFS). In Monte Carlo simulation, the simulations are repeated multiple times, covering the whole range of variability of the input parameters and thus resulting in a distribution of results covering what can be expected in reality. In practical ...
Obtaining the set of cosmological parameters consistent with observational data is an important exercise in current cosmological research. It involves finding the global maximum of the likelihood function in the multi-dimensional parameter space. Currently sampling based methods, which are in general stochastic in nature, like Markov-Chain Monte Carlo(MCMC), are being commonly used for parameter estimation. The beauty of stochastic methods is that the computational cost grows, at the most, linearly in place of exponentially (as in grid based approaches) with the dimensionality of the search space. MCMC methods sample the full joint probability distribution (posterior) from which one and two dimensional probability distributions, best fit (average) values of parameters and then error bars can be computed. In the present work we demonstrate the application of another stochastic method, named Particle Swarm Optimization (PSO), that is widely used in the field of ...
The spectrum of antiprotons from dark matter annihilation are calculated using the Lund Monte Carlo program, and simple analytic expressions for the spectrum and low-energy antiproton/proton ratio are derived. Comparing the results with recent upper limits on low energy antiprotons, it is concluded that the reported 4-13 GeV antiproton flux cannot be accounted for by dark matter annihilation. The new upper limits do not provide useful constraints on dark matter particles. They restrict the annihilation rate and imply that annihilation gamma ray and e(+) fluxes would be far below the fluxes produced by cosmic-ray collisions. It may be possible to look for a dark matter halo annihilation signal at antiprotons energies below 0.5 GeV, where the flux from cosmic-ray collisions is expected to be negligible. 32 references.
Air toxics from the industrial wastewater treatment plants (IWTPs) impose serious health concerns on its surrounding residential neighborhoods. To address such health concerns, one of the key challenges is to quantify the air emissions from the IWTPs. The objective here is to characterize the air emissions from the IWTPs and quantify its associated uncertainty. An IWTP receiving the wastewaters from an airplane maintenance facility is used for illustration with focus on the quantification of air emissions for benzyl alcohol, phenol, methylene chloride, 2-butanone, and acetone. Two general fate models, i.e., WATER9 and TOXCHEM+V3.0 were used to model the IWTP and quantify the air emissions. Monte Carlo and Bootstrap simulation were used for uncertainty analysis. On average, air emissions fr...
We have studied actinide burn-up in ICF reactor pellets; i.e., 14 MeV neutron fission of the very long-lived actinides that pose storage problems. A major advantage of pellet fuel region burn-up is safety: only milligrams of highly toxic and active material need to be present in the fusion chamber, whereas blanket burn-up requires the continued presence of tons of actinides in a small volume. The actinide data tables required for Monte Carlo calculations of the burn-up of /sup 241/Am and /sup 243/Am are discussed in connection with a study of the sensitivity to cross section uncertainties. More accurate and complete cross sections are required for realistic quantitative calculations.
We present results of multiple-time-scale simulations of 5, 10 and 15 keV low temperature ion implantation of arsenic on silicon (100), followed by high temperature anneals. The simulations start with a molecular dynamics (MD) calculation of the primary state of damage after 10ps. The results are then coupled to a kinetic Monte Carlo (MC) simulation of bulk defect diffusion and clustering. Dose accumulation is achieved considering that at low temperatures the damage produced in the lattice is stable. After the desired dose is accumulated, the system is annealed at 800 degrees C for several seconds. The results provide information on the evolution for the damage microstructure over macroscopic length and time scales and affords direct comparison to experimental results. We discuss the database of inputs to the MC model and how it affects the diffusion process.
We review our recent work on an atomistic approach to the development of predictive process simulation tools. First principles methods, molecular dynamics simulations, and experimental results are used to construct a database of defect and dopant energetics in Si. This is used as input for kinetic Monte Carlo simulations. C and B trapping of the Si self- interstitial is shown to help explain the enormous disparity in its measured diffusivity. Excellent agreement is found between experiments and simulations of transient enhanced diffusion following 20-80 keV B implants into Si, and with those of 50 keV Si implants into complex B-doped structures. Our simulations predict novel behavior of the time evolution of the electrically active B fraction during annealing.
Most active pharmaceutical ingredients (API) exhibit particle size distributions with some degrees of asymmetry deviating from log-normality. A new log-skew-normal (L-S-N) distribution model is proposed for a systematic comparison of the asymmetry effect on content uniformity. The new model originated from the S-N model used by Azzalini gives a close approximation to real API particle size distribution. Monte-Carlo method was employed to simulate the dosage potency distribution. A high risk of over potency is uncovered when either the dose is low or API particle size distribution is positively skewed. This is due to the formation of pseudo heavy tail in potency distribution that decays slower than exponentially. Nomographs of API particle size versus dosage strength were constructed with a...
In this paper, the performance of customer-owned distributed generation (DG) units is quantified from different perspectives through an uncertainty study. A Monte Carlo-based method is applied to assess the stochastic operation of the customer-owned DG units in the power distribution system. Several cases are studied to analyze the impact on system performance of using such generators, with the emphasis on benefits. The results of the studied cases show that proper operation of customer-owned DG units placed close to significant consumption centers offers several benefits which lead to significant energy savings and improvement in the performance indices while maintaining the cost-effectiveness. Furthermore, based on the energy demand, different electricity price scenarios considering a co...
We compare the metallicities of stars with radial velocity planets to the metallicity of a sample of field dwarfs. We confirm recent work indicating that the stars-with-planet sample as a whole is iron rich. However, the lowest mass stars tend to be iron poor, with several having [Fe/H]0.48) that contributes to but does not explain the mass-metallicity trend in the stars-with-planets sample. We use Monte Carlo models to show that adding an average of 6.5 Earth masses of iron to each star can explain both the mass-metallicity and the age-metallicity relations of the stars-with-planets sample. However, for at least one star, HD 38529, there is good evidence that the bulk metallicity is high. We conclude that the observed metallicities and metallicity trends are the result of the interaction of three effects; accretion of about 6 Earth masses of iron rich material, selection effects, and in some cases, high intrinsic metallicity.
We calculate the energy spectra of cosmic rays (CR) and their secondaries produced in a supernova remnant (SNR), taking into account the time-dependence of the SNR shock. We model the trajectories of charged particles as a random walk with a prescribed diffusion coefficient, accelerating the particles at each shock crossing. Secondary production by CRs colliding with gas is included as a Monte Carlo process. We find that SNRs produce less antimatter than suggested previously: The positron/electron ratio and the antiproton/proton ratio are a few percent and few $\\times 10^{-5}$, respectively. Moreover, the obtained positron/electron ratio decreases with energy, while the antiproton/proton ratio rises at most by a factor of two above 10 GeV.
We carried out theoretical analyses based on ab initio calculations that incorporate free energy of the vapor phase in order to determine the initial growth process of cubic GaN on GaN (001)-(4x1). The results suggest that a N-adsorbed structure appears at the initial growth stage and then Ga adsorbs on the N-adsorbed GaN (001)-(4x1) surface. Considering this process, we performed Monte Carlo simulations. The results suggest that the maximum point of Ga coverage after supplying a Formula Not Shown monolayer of atoms shifted toward a Ga-rich condition from V/III=1.0.
It has been shown by recent experiments that bulk immiscible metals (e.g. Ag/Cu, Ag/Co and Au/Ni) can form binary alloys on certain surfaces where the substrate mediates the elastic misfits between the two components, thus relieving the elastic strain in the overlayer. These novel surface alloys exhibit a rich phase structure. We formulate a two-component Frenkel-Kontorova model in one dimension to study surface alloy formation. This model can naturally incorporate dislocation formation that plays a crucial role in determining the actual structure of the system. Using energy minimization calculations we provide a phase diagram in terms of average alloy composition and the energy of mixing. Monte Carlo simulations were also performed to study the structure and interaction of the emerging dislocations.
A complete code based on the finite element method (FEM) is used to study the attenuation of gamma rays, of source energy 9 MeV, in simple homogeneous lead and water systems. The results are compared, in tables and graphs, with those produced by a Monte Carlo method of solution. In the comparison, particular attention is paid to the annihilation gamma rays generated by the pair production process, to energy deposition, and to the leakage spectra. A model, multilayer, slab shield, having a fission neutron source impinging on the left-hand face, is also solved by the FEM code, using data from a coupled neutron/gamma ray cross section library. The resulting neutron and gamma ray flux distributions, and the leakage spectra, are shown graphically. (author).
A classical model is presented for magnetic field-induced Wigner crystallization in electron systems confined within two-dimensional quantum dots. In contrast to other classical models, this one does not treat an electron as a point charge; the electron density is assumed to take a Gaussian form corresponding to the lowest Landau level. Using a Monte Carlo method we have determined the equilibrium configurations as functions of the magnetic field. We have found a classical counterpart of the quantum maximum density droplet (MDD) and studied the breakdown of the MDD into a Wigner molecule as well as the transformations of the Wigner molecule shape induced by the external magnetic field. The phase diagram for the classical Wigner molecules has been presented and its qualitative agreement with previous quantum mechanical calculations has been shown.
The second-order spatial subwavelength interference pattern is observed in a modified Michelson interferometer with single-mode continuous-wave laser beams. By analyzing our subwavelength interference experiment based on Feynman's path integral theory, a unified interpretation for all the second-order subwavelength interference is suggested.
The results for development of methods and computer programs for integrated optimization of parameters of perspective fast reactors are given. The possibilities of the program for the reactor campaign calculation are analysed. This program is based on utilisation of the Bubnov-Galerkin method and Wigner disturbance theory. The possibility of application of approximation methods for the optimization researches is discussed. The results of development of the programs for complex reactor computations with account of control rods system and change of physical parameters in the reactor campaign are discussed. (author).
The thermal decrease of the order parameter can empirically be described by a single T{sup {epsilon}} power law with an exponent {epsilon} which depends on the dimensionality of the magnetic interactions and on whether the spin quantum number is integral or half-integral. We present experimental examples in which the order parameter shows a crossover between different T{sup {epsilon}} power laws as a function of temperature. This indicates that the magnetic interactions can change their dimensionality as a function of temperature. (orig.)
This paper dealt with the identification of possible damaging mechanism of the collector of the WWER 440 steam generator, cracking of primary collectors, corrosion damage of the protective coat of the primary collector circumferential weld, cracking of breathing space in the region of blinding effect by corrosion and strain, leaking of disassembling joint of the primary collector lid and with the integrity of heat exchanging tubes.
Operational Final Test Report for Integrated Systems, Project W-030 (Phase 2 test, RECIRC and HIGH-HEAT Modes). Project W-030 provides a ventilation upgrade for the four Aging Waste Facility tanks, including upgraded vapor space cooling and filtered venting of tanks AY101, Ay102, AZ101, AZ102.
Integrating technology readiness levels (TRL) into the management of engineering projects is critical to the mitigation of risk and improved customer/supplier communications. TRLs provide a common framework and language with which consistent comparisons of different technologies and approaches can be made. At Sandia National Laboratories, where technologies are developed, integrated and deployed into high consequence systems, the use of TRLs may be transformational. They are technology independent and span the full range of technology development including scientific and applied research, identification of customer requirements, modeling and simulation, identification of environments, testing and integration. With this report, we provide a reference set of definitions for TRLs and a brief history of TRLs at Sandia National Laboratories. We then propose and describe two approaches that may be used to ...
Numerical errors of PALLAS calculation due to spatial mesh sizes are examined for a typical deep penetration shielding problem of isotropic incident fission neutrons penetrating a 200-cm-thick water slab. The exponential approximation for the source spatial distribution to solve the transport equation based on the direct integration method is verified to be effective for radiation transport in attenuating medium, while the linear approximation employed in the old PALLAS code is found to produce large errors for calculations with large mesh sizes.
To optimize the emergency care of patients with multiple injuries, a special emergency room with integrated radiological equipment has been developed. It consists of a linear array of a computed tomography system, a mobile C-arm with an integrated scattered radiation grid for conventional films and a C-arm fluorographic system. The intent is to avoid delay caused by transport and repositioning of the patient. It enables all necessary diagnostic procedures to be performed within a minimal amount of time. (orig./MG).
The authors propose an integrated indicator of pump-rod couplings that allows both the couplings and the pump operations to be appraised according to the given formula. (Formula provided). The dynamic relationships of rod operations were determined with nomographs. These relationships involve such factors as: the type and size of the sucker-rod string; the pressure load at the equalizer head and its correlation to threshold pressure at pump discharge; pump diameter; the rod weight and construction; and integral reliability indicators.
This document provides structural load requirements and their basis for maintaining the structural integrity of the Hanford Single-Shell Tanks during waste feed delivery and retrieval activities. The requirements are based on a review of previous requirements and their basis documents as well as load histories with particular emphasis on the proposed lead transfer feed tanks for the privatized vitrification plant.
An integrated mirrorless bistable optical device based on the Mach-Zehnder interferometric optical switch has been proposed and demonstrated experimentally using a Ti:LiNbO3 waveguide. The resulting device is capable of combining more than two of them to realize multifunctional optical devices such as optical multivibrators.
U.S. Department of Energy (DOE) Policy 450.4, Safety Management System Policy commits to institutionalizing an Integrated Safety Management System (ISMS) throughout the DOE complex as a means of accomplishing its missions safely. DOE Acquisition Regulation 970.5204-2 requires that contractors manage and perform work in accordance with a documented safety management system.
The reports in this series are designed primarily to familiarize professional land use and resource planners with the range of possibilities and effective procedures for achieving integrated mining, reclamation, and land use planning. These reports are ba...
The evaluated data testing for Si, Zr, Nb and Fe materials has been performed through the analyses of the neutron leakage spectrum from spherical assemblies measured in integral benchmark experiments with (d,t) and Cf fission neutron sources. Intercomparisons of the calculated from BROND-2, ENDF/B-6 and JENDL-3 data files with experimental results are presented.
We construct the Generalized Monodromy matrix $\\mathcal{\\hat{M}}(\\omega)$ of two dimensional string effective action by introducing the T-duality group properties.The integrability conditions with general solutions depending on spectral parameter are given. This construction is investigated for the exactly solvable Wess, Zumino, Novikov and Witten (WZNW) model in pp-wave Limit when B=0.
The paper describes the Integrated Air Pollution Control System (IAPCS), a computerized program that can be used to estimate the cost and performance of pre-combustion, in situ, and post-combustion air pollution control configurations in pulverized-coal-fired utility boilers of 1...
We report on the calculation of multi-loop Feynman integrals for single-scale problems by means of difference equations in Mellin space. The solution to these difference equations in terms of harmonic sums can be constructed algorithmically over difference fields, the so-called {pi}{sigma}{sup *}-fields. We test the implementation of the Mathematica package Sigma on examples from recent higher order perturbative calculations in Quantum Chromodynamics. (orig.)
Certain problems of gamma-absorption method of measuring the thickness of layers in multi-layer items are analyzed. Two examples of solving the above-mentioned problems have been considered. One of them deals with a two-layer item, its integral thickness being known or measured in advance, the second example is referred to a two-layer item with unknown integral value
In this paper we first establish global pointwise time-space estimates for a class of oscillatory integrals. Then, we use them to establish $L^p-L^q$ estimates for a class of higher order wave-type equations of the form $\\partial_{tt}u+P(D_{x})u=0$, where the symbol $P(\\xi)$ is a real non-degenerate elliptic polynomial of ${\\bf R}^n$.
Comprehensive Structural Integrity is a reference work which covers all activities involved in the assurance of structural integrity. It provides engineers and scientists with an unparalleled depth of knowledge in the disciplines involved. The new online Volume 11 is dedicated to the mechanical characteristics of materials. This paper contains the chapter 11.02 of this volume and is structured as follows: Test techniques; Analysis; Fracture behavior; Fracture toughness tests for nonmetals.
An integrated voltage regulator and power system stabilizer based on one step ahead prediction is presented. The fixed parameter predictive control scheme is shown capable of providing consistently good voltage response and contribution to system damping over a wide range of operating conditions and system configuration. This is in contrast to the case with conventional automatic voltage regulator and power system stabilizer where performance varies widely with the operating situation and can even cause system instability. (author)
We discuss the determination of the infrared singularities of massive one-loop 5-point functions with Mellin-Barnes (MB) representations. Massless internal lines may lead to poles in the {epsilon} expansion of the Feynman diagram, while unresolved massless final state particles give endpoint singularities of the phase space integrals. MB integrals are an elegant tool for their common treatment. An evaluation by taking residues leads to inverse binomial sums. (orig.)
The {sup 252}Cf-source-driven noise analysis method has been used in measurements for subcritical configurations of fissile systems for a variety of applications. Measurements of 25 fissile systems have been performed with a wide variety of materials and configurations. This method has been applied to measurements for (1) initial fuel loading of reactors, (2) quality assurance of reactor fuel elements, (3) fuel preparation facilities, (4) fuel processing facilities, (5) fuel storage facilities, (6) zero-power testing of reactors, and (7) verification of calculational methods for assemblies with the neutron k < l. These previous measurements, performed with a wide variety of multiplying systems, demonstrated the usefulness of the method. The high sensitivity of noise-measured parameters to small changes in fissile systems has been observed in several measurements. This high sensitivity has been evaluated using Monte Carlo neutron transport methods for an array of ...
Within the framework of the EFTTRA cooperation between CEA, ECN, EDF, FZK, IAM and ITU, six metallic {sup 99}Tc rods have been irradiated in the Petten HFR for 193 effective full power days. During this irradiation, more than 6% of the {sup 99}Tc has been transmuted to the stable {sup 100}Ru. At ECN, one of the six rods has been examined in the hot cell laboratory. The ruthenium concentration in the rod measured by Isotope Dilution Mass Spectrometry reaches 6.4% at 5 mm from the bottom of the rod and 6.0% at 5 mm from the top. Also the axial and radial distributions of the ruthenium have been measured by Electron Probe Micro Analysis. The ruthenium concentrations calculated by the three-dimensional Monte Carlo code KENO reach 6.1% at 5 mm from the bottom of the rod and 5.7% at 5 mm from the top. These values are in reasonable agreement with the measured ones. However, the calculated radial distribution of the ruthenium concentration is not in agreement with the ...
Monte Carlo (MC) dose calculations can be accurate but are also computationally intensive. In contrast, convolution superposition (CS) offers faster and smoother results but by making approximations. We investigated MC denoising techniques, which use available convolution superposition results and new noise filtering methods to guide and accelerate MC calculations. Two main approaches were developed to combine CS information with MC denoising. In the first approach, the denoising result is iteratively updated by adding the denoised residual difference between the result and the MC image. Multi-scale methods were used (wavelets or contourlets) for denoising the residual. The iterations are initialized by the CS data. In the second approach, we used a frequency splitting technique by quadrature filtering to combine low frequency components derived from MC simulations with high frequency components derived from CS components. The rationale is to take the scattering ...
The verification of the subcriticality is of utmost importance for the safe transportation and storage of nuclear reactor fuels. Transportation containers and storage facilities are designed such that nuclear fuels remain in a subcritical state. Such designs often involve excess conservatism because of the lack of relevant experimental data to verify the accuracy of Monte Carlo codes used in nuclear criticality safety analyses. A joint experimental research program between Oak Ridge National Laboratory, Westinghouse Safety Management Solutions, Inc., and the University of Missouri was initiated to obtain measured quantities that could be directly related to the subcriticality of simple arrays of Missouri University Research Reactor (MURR) fuel elements. A series of measurement were performed to assess the reactivity of materials such as BORAL, stainless steel, aluminum, and lead that are typically used in the construction of shipping casks. These materials were ...
In this thesis, the potential to discover supersymmetric signatures in both e"+e"- and #gamma##gamma# colliders evaluated with a Monte Carlo analysis, is discussed. The analysis was focused on the detection of muons, essentially. First, we study the detection of #mu#"+_R#mu#"-_R in e"-e"+ collisions, whose purpose is the measurement of the #mu#_R mass. It was found that an uncertainty of 0.11 GeV (stat) can be achieved. Under the assumption of the real conditions of the ILC photon collider, a study covering the detection of #mu#"+_R#mu#"-_R and #mu#"+_L#mu#"-_L pairs for #sq root#(s_e_"-_e_"+)=0.5 and 0.6 TeV was performed. According to the simulation, a statistical error for the branching ratios of #mu#"#+-#_L#->##chi#"0_1#mu#"#+-# of 0.98% and #mu#"#+-#_L#->##chi#"0_2#mu#"#+-# of 3.97% can be reached. In order to judge the Monte Carlo results, we have used a technique based on a multidimensional fit to evaluate the impact of the ...
For the study of radiation biology and its application to radiotherapy, the double differential cross section of electron emission from water vapor induced by 6.0 MeV alpha particle beam is measured. The energy spectra of electrons ranging 7- 10000 eV are detected by the electrostatic analyzer and micro channel plate. The measurements are made at angles between 20 and 160 degrees. With use of this data set, the radial dose distribution in water is calculated by using KURBUC code. It is the Monte Carlo type code of the electron transport process, where the track of the electron is simulated through each individual interactions including elastic scattering, ionization cross section and total excitation cross section in case that electrons with certain energy are put in the liquid-density water. In order to understand the effect of radiation when the particle flux is injected in the human body like radiotherapy using accelerator beam, the dose distribution in the ...
The Anisotropic Analytical Algorithm (AAA) is a new pencil beam convolution/superposition algorithm proposed by Varian for photon dose calculations. The configuration of AAA depends on linear accelerator design and specifications. The purpose of this study was to investigate the accuracy of AAA for an Elekta SL25 linear accelerator for small fields and intensity modulated radiation therapy (IMRT) treatments in inhomogeneous media. The accuracy of AAA was evaluated in two studies. First, AAA was compared both with Monte Carlo (MC) and the measurements in an inhomogeneous phantom simulating lung equivalent tissues and bone ribs. The algorithm was tested under lateral electronic disequilibrium conditions, using small fields (2x2 cm"2). Good agreement was generally achieved for depth dose and profiles, with deviations generally below 3% in lung inhomogeneities and below 5% at interfaces. However, the effects of attenuation and scattering close to the bone ribs were not ...
We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolution codes SSE and BSE from Hurley et al. We describe the modifications that we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared with simulations without any stellar evolution. In particular, we find that the mass loss from the stellar evolution acts as a significant energy production channel simply by reducing the total gravitational binding energy and can significantly ...
A series of first principles molecular dynamics and Monte Carlo simulations were carried out for liquid water to assess the validity and reproducibility of different sampling approaches. These simulations include Car-Parrinello molecular dynamics simulations using the program CPMD with different values of the fictitious electron mass in the microcanonical and canonical ensembles, Born-Oppenheimer molecular dynamics using the programs CPMD and CP2K in the microcanonical ensemble, and Metropolis Monte Carlo using CP2K in the canonical ensemble. With the exception of one simulation for 128 water molecules, all other simulations were carried out for systems consisting of 64 molecules. It is found that the structural and thermodynamic properties of these simulations are in excellent agreement with each other as long as adiabatic sampling is maintained in the Car-Parrinello molecular dynamics simulations either by choosing a sufficiently small ...
We have developed a kinetic Monte Carlo (kMC) simulator that links atomic migration and binding energies determined primarily from first principles calculations to macroscopic phenomena and laboratory time scales. Input for the kMC simulation is obtained from a combination of ab initio planewave pseudopotential calculations, molecular dynamics simulations, and experimental data. The simulator is validated against an extensive series of experimental studies of the diffusion of B spikes in self-implanted Si. The implant energy, dose, and dose rate, as well as the detailed thermal history of the sample, are included. Good agreement is obtained with the experimental data for temperatures between 750 and 950 C and times from 15 to 255 s. At 1050o C we predict too little diffusion after 105 s compared to experiment: apparently, some mechanism which is not adequately represented by our model becomes important at this temperature. Below 1050o C, the kMC simulation produces ...
Currently, x-ray mammography is the method of choice in breast cancer screening programmes. As the mammography technology moves from 2D imaging modalities to 3D, conventional computational phantoms do not have sufficient detail to support the studies of these advanced imaging systems. Studies of these 3D imaging systems call for a realistic and sophisticated computational model of the breast. DeBRa (Detailed Breast model for Radiological studies) is the most advanced, detailed, 3D computational model of the breast developed recently for breast imaging studies. A DeBRa phantom can be constructed to model a compressed breast, as in film/screen, digital mammography and digital breast tomosynthesis studies, or a non-compressed breast as in positron emission mammography and breast CT studies. Both the cranial-caudal and mediolateral oblique views can be modelled. The anatomical details inside the phantom include the lactiferous duct system, the Cooper ligaments and the pectoral muscle. The ...
Adsorption of ethylene and ethane on graphitized thermal carbon black and in slit pores whose walls are composed of graphene layers is studied in detail to investigate the packing efficiency, the two-dimensional critical temperature, and the variation of the isosteric heat of adsorption with loading and temperature. Here we used a Monte Carlo simulation method with a grand canonical Monte Carlo ensemble. A number of two-center Lennard-Jones (LJ) potential models are investigated to study the impact of the choice of potential models in the description of adsorption behavior. We chose two 2C-LJ potential models in our investigation of the (i) UA-TraPPE-LJ model of Martin and Siepmann for ethane and Wick et al. for ethylene and (ii) AUA4-LJ model of Ungerer et al. for ethane and Bourasseau et al. for ethylene. These models are used to study the adsorption of ethane and ethylene on graphitized thermal carbon black. It is found that the solid-fluid ...
The convolution/superposition calculations for radiotherapy dose distributions are traditionally performed by convolving polyenergetic energy deposition kernels with TERMA (total energy released per unit mass) precomputed in each voxel of the irradiated phantom. We propose an alternative method in which the TERMA calculation is replaced by random sampling of photon energy, direction and interaction point. Then, a direction is randomly sampled from the angular distribution of the monoenergetic kernel corresponding to the photon energy. The kernel ray is propagated across the phantom, and energy is deposited in each voxel traversed. An important advantage of the explicit sampling of energy is that spectral changes with depth are automatically accounted for. No spectral or kernel hardening corrections are needed. Furthermore, the continuous sampling of photon direction allows us to model sharp changes in fluence, such as those due to collimator tongue-and-groove. The use of explicit ...
In large samples, the {gamma}-ray count rate of a prompt gamma neutron activation analysis system is a multi-variable function of the elemental dry composition, density, water content and thickness of the material. The experimental calibration curves require tremendous laboratory work, using a great number of standards with well-known compositions. Although a Monte Carlo simulation study does not avoid the experimental calibration work, it reduces the number of experimental calibration standards. This paper is part of a feasibility study for a PGNAA system for on-line continuous characterisation of cement raw material conveyed on a belt (Oliveira, C., Salgado, J. and Carvalho, F. G. (1997) Optimisation of PGNAA instrument design for cement raw materials using the MCNP code. J. Radioanal. Nucl. Chem. 216(2), 191-198; Oliveira, C., Salgado, J., Goncalves, I. F., Carvalho, F. G. and Leitao, F. (1997a) A Monte Carlo study of the influence of ...
Under support from the French Space Agency (CNES), a 3D+t dusty-gas model of Comet 67P/Churyumov-Gerasimenko is being developed, to compute, from the first 2014 Rosetta orbital data, the aerodynamic forces exerted on the Rosetta orbiter and on the descent lander. We report the recently developed dust dynamics part of the code. The multi-species (presently H2O and CO) gas code is optimized in terms of computational speed owing to the use of two complementary methods: (a) 3D+t Direct Simulation Monte Carlo (DSMC) runs in the non-equilibrium regions adjacent to the surface and very distant from it, and (b) solutions of the Navier-Stokes equations in-between. The model is used presently using Lamy et al. (Space Sci. Rev., 2007, 128, 23) coarse information on 67P nucleus shape and rotation, and a range of possible gas production rates Q for the early Rosetta observations at rh 3 AU (Q 1026 - 1027 s-1). In the interim version, simplifying assumptions are made with ...
Particle swarm optimization (PSO) is extensively used for real parameter optimization in diverse fields of study. This paper describes an application of PSO to the problem of designing a fractional-order proportional-integral-derivative (FOPID) controller whose parameters comprise proportionality constant, integral constant, derivative constant, integral order (lambda) and derivative order (delta). The presence of five optimizable parameters makes the task of designing a FOPID controller more challenging than conventional PID controller design. Our design method focuses on minimizing the Integral Time Absolute Error (ITAE) criterion. The digital realization of the deigned system utilizes the Tustin operator-based continued fraction expansion scheme. We carry out a simulation that illustrates the effectiveness of the proposed approach especially for realizing fractional-order plants. This paper also ...
College student persistence and dropout have been studied for decades, but little inquiry has focused on community college or private two-year college students. Although about half of first-time postsecondary students enroll in a two-year college, researchers understand little about why only approximately a quarter of these degree-seekers complete any degree five years after entrance (Rosenbaum, Deil-Amen, & Person, 2006). Models that exist have been generated primarily from research on residential, four-year college students. The present analysis uses rich qualitative data to excavate the potential of Tinto's Theory of Student Departure for understanding the dynamics of persistence for a more nontraditional group--two-year college commuting students. To better reflect how these students successfully cultivate feelings of belonging and competence, reinforce goal commitment, and access valuable social capital, a new conceptual fusion of the formerly distinct processes of academic ...
A review of the status of integral measurements is presented for "2"4"0Pu, "2"4"1Pu, "2"4"2Pu, "2"4"1Am and "2"4"3Am. This review includes integral measurements pertinent to thermal reactor systems, i.e., thermal cross sections and resonance integrals, as well as measurements for fast reactor systems. It appears that for these nuclides the data for thermal reactors are in good shape; however, more work is recommended in defining the branching ratio of the capture cross section of "2"4"1Am to the isomeric and ground states of "2"4"2Am. Also, benchmark irradiation data are needed for cross section data testing using depletion/production codes. For fast reactors, experiments are in progress, in the UK, in France, and also in the US, with partial results available at this time. Fast integral data obtained from these measurements will be very beneficial. The recommendation pertaining to "2"4"1Am and proper ...
Interest in renewable energy systems that employ regenerative components to enable intermittent energy sources to service time varying loads is growing. The overall performance of these systems is intimately linked to the choice of the energy storage media and the efficiency of the energy conversion devices. Recent advances in electrolyser and fuel cell technologies may make small-scale hydrogen buffered renewable energy systems a practical reality in the near term. A research team at the University of Victoria - Institute for Integrated Energy Systems (IESVic) is developing a flexible distributed, laboratory-scale energy system called IRENE (Integrated Renewable Energy Experiment). The initial goal of this research is to explore the potential for constructing a regenerative energy system using early commercial and precommercial hydrogen storage and fuel cell energy conversion devices. Once implemented, research efforts will focus on the ...
Formulation is introduced for discretizing a boundary integral equation into an indirect boundary element method for the solution of 3-dimensional topographic problems. Yokoi and Takenaka propose an analytical solution-capable reference solution (solution for the half space elastic body with flat free surface) to problems of topographic response to seismic motion in a 2-dimensional in-plane field. That is to say, they propose a boundary integral equation capable of effectively suppressing the non-physical waves that emerge in the result of computation in the wake of the truncation of the discretized ground surface making use of the wave field in a semi-infinite elastic body with flat free surface. They apply the proposed boundary integral equation discretized into the indirect boundary element method to solve some examples, and succeed in proving its validity. In this report, the equation is expanded to deal with ...
We derive a full Bern-Kosower-type rule for scalar QED starting from quantum field theory: we derive a set of rules for calculating S-matrix elements for any processes at any order of the coupling constant. A gauge-invariant set of diagrams in general is first written in the world line path-integral expression. Then we integrate over x(#tau#), and the resulting expression is given in terms of a correlation function on the world line left-angle x(#tau#)x(#tau#"')right-angle. Simple rules to decompose the correlation function into basic elements are obtained. A gauge transformation known as the integration by parts technique can be used to reduce the number of independent terms before integration over proper-time variables. The surface terms can be omitted provided the external scalars are on shell. Also, we clarify correspondence to the conventional Feynman rule, which enabled us to avoid any ambiguity ...
The integrity of the cables located in the target room is very important in the Fire PSA, because the CDF and CCDP are changed according to the results of a cable integrity that depends on the surrounding gas temperature. The conservative assumptions used in the Fire PSA typically specify that all of the equipment and cables of a room would fail when a fire happens in the room. But the realistic assessment of a fire risk by using a fire simulation tool has become necessary in the Fire PSA as described in the ANS Fire PRA Standard. This paper evaluates the cable integrity of eight pump rooms in the nuclear power plant by using the CFAST zone fire model. The upper layer gas temperature of each room is estimated, and an analysis based on the results of model simulations is used to judge the cable integrity. According to the analysis results, the integrity of the cable located in the ...
The analysis of oceanographic lidar systems measurements is often carried out with semi-empirical methods, since there is only a rough understanding of the effects of many environmental variables. The development of techniques for interpreting the accuracy of lidar measurements is needed to evaluate the effects of various environmental situations, as well as of different experimental geometric configurations and boundary conditions. A Monte Carlo simulation model represents a tool that is particularly well suited for answering these important questions. The PREMAR-2F Monte Carlo code has been developed taking into account the main molecular and non-molecular components of the marine environment. The laser radiation interaction processes of diffusion, re-emission, refraction and absorption are treated. In particular are considered: the Rayleigh elastic scattering, produced by atoms and molecules with small dimensions with respect to the laser ...
A fully implicit compositional simulator was written to predict the drawdown and buildup responses for a single well in a gas condensate reservoir. The model uses the Redlich-Kwong equation of state to compute fluid properties and to determine phase equilibrium. The accuracy of the model predictions in the two-phase flow regime is tested using the results of a steady-state theory developed here. Two new correlating functions for transient flow - the reservoir integral and the sandface integral - are presented. The reservoir integral is shown to be a good liquid analogue during transient flow (drawdown). The sandface integral correlates with the liquid solution only in the slope. These two theoretical functions are shown to reduce to a computable pseudopressure function by using the results of the steady-state theory. The usefulness of this pseudopressure function for computing the permeability-thickness ...
The infection of cultured human cells with baboon endogenous virus (BEV) frequently leads to an association of viral DNA with a specific genetic locus (termed BEVI, for baboon endogenous virus infection) on chromosome 6. Restriction endonuclease digestion of DNA from BEV-infected human cells and their derived somatic cell clones frequently revealed a common cellular DNA sequence in the proximity of one of the junctions between cellular DNA and the integrated virus. We propose that a short cellular DNA sequence, repeated on chromosome 6 and separated by unique DNA sequences, presents a high-affinity target for the integration of BEV in human cells. PMID:6401843
The Systems Analysis Programs for Hands-on Integrated Reliability Evaluations (SAPHIRE) refers to a set of several microcomputer programs that were developed to create and analyze probabilistic risk assessments (PRAs), primarily for nuclear power plants. The Fault Tree, Event Tree, and Piping and Instrumentation Diagram (FEP) editors allow the user to graphically build and edit fault trees, and event trees, and piping and instrumentation diagrams (P and IDs). The software is designed to enable the independent use of the graphical-based editors found in the Integrated Reliability and Risk Assessment System (IRRAS). FEP is comprised of three separate editors (Fault Tree, Event Tree, and Piping and Instrumentation Diagram) and a utility module. This reference manual provides a screen-by-screen guide of the entire FEP System.
The sequence of events from the assembly of a space nuclear power system to its integration in the Space Shuttle Transportation System (STS) is considered. First, the sequence followed for SNAP-10A, the only free world space reactor electric power system ever launched and operated in space, is reviewed. Before shipment, the SNAP-10A reactor was raised to operating temperature using electrically supplied heat and operated at low power for control calibration. Next we discuss shipment to the launch site, a phase that is critical because of the potential for various accidents. Once the power system arrives at the launch site, the processing sequence is performed. This sequence includes checkout, mating with the payload or upper stage launch vehicle, and integration into the STS.
Results that contradict known physics, data manipulated, lack of vigilance by co-authors, failures in the system for scientific publication... Last September a US Committee of Enquiry unveiled one of the most serious frauds in the history of physics. Over a two year period, a young researcher at Bell Laboratories had published a large number of articles with exciting results for solid state physics, but which, alas, were fraudulent! Obviously a fraud of this magnitude is exceptional. However, it did serve to focus attention on the problem of integrity and honesty in research practices. This subject, crucial to the well-being and credibility of scientific research, will be the central theme of the lecture given by Nicholas Steneck, Professor of History at the University of Michigan. A leading expert on this issue, on which he has published extensively, he is a consultant to the Office of Research Integrity in the US, and has been closely ...
Speech perception integrates auditory and visual information. This is evidenced by the McGurk illusion where seeing the talking face influences the auditory phonetic percept and by the audiovisual detection advantage where seeing the talking face influences the detectability of the acoustic speech signal. Here, we show that identification of phonetic content and detection can be dissociated as speech-specific and non-specific audiovisual integration effects. To this end, we employed synthetically modified stimuli, sine wave speech (SWS), which is an impoverished speech signal that only observers informed of its speech-like nature recognize as speech. While the McGurk illusion only occurred for informed observers, the audiovisual detection advantage occurred for na?ve observers as well. Thi...
We examine moving average (MA) filters for estimating the integrated variance (IV) of a financial asset price in a framework where high-frequency price data are contaminated with market microstructure noise. We show that the sum of squared MA residuals must be scaled to enable a suitable estimator of IV. The scaled estimator is shown to be consistent, first-order efficient, and asymptotically Gaussian distributed about the integrated variance under restrictive assumptions. Under more plausible assumptions, such as time-varying volatility, the MA model is misspecified. This motivates an extensive simulation study of the merits of the MA-based estimator under misspecification. Specifically, we consider nonconstant volatility combined with rounding errors and various forms of dependence between the noise and efficient returns. We benchmark the scaled MA-based estimator to subsample and realized kernel estimators and find that the MA-based ...
Most gas dynamic computations in industrial ducts are done in one dimension with cross-section-averaged Euler equations. This poses a fundamental difficulty as soon as geometrical discontinuities are present. The momentum equation contains a non-conservative term involving a surface pressure integral, responsible for momentum loss. Definition of this integral is very difficult from a mathematical standpoint as the flow may contain other discontinuities (shocks, contact discontinuities). From a physical standpoint, geometrical discontinuities induce multidimensional vortices that modify the surface pressure integral. In the present paper, an improved 1D flow model is proposed. An extra energy (or entropy) equation is added to the Euler equations expressing the energy and turbulent pressure ...
Conventional interconnect and switching technology is rapidly becoming a critical issue in the realization of systems using high speed silicon and GaAs based technologies. In recent years clock speeds and on-chip density for VLSI/VHSIC technology has made packaging these high speed chips extremely difficult. A strong case can be made for using optical interconnects for on-chip/on-wafer, chip-to-chip and board-to-board high speed communications. GaAs integrated optoelectronic circuits (IOC's) are being developed in a number of laboratories for performing Input/Output functions at all levels. In this paper integrated optoelectronic materials, electronics and optoelectronic devices are presented. IOC's are examined from the standpoint of what it takes to fabricate the devices and what performance can be expected.
Knowledge is a success factor for globally acting company networks as a. o. the aim-oriented knowledge transfer between partners is an essential condition for a successful cooperation. Several network specific problems impede however an efficient and effective knowledge management; e.g. the transfer of competition relevant data is a high sensitive theme. The authors describe how a relevant method can be selected, adapted and implemented to the application-specific boundary conditions. They integrate existing attempts processes and methods to establish an effective knowledge development and an efficient knowledge transfer in the network and take into account ''hard'' (IT-technical-oriented solution approaches) as also ''soft'' factors (e.g. cultural and personal aspects). So the authors present unified instruments for an integrated knowledge management in networks. (GL)
This self-evaluation report offers a summary of results from FY2000 actions to achieve Pacific Northwest National Laboratory's strategy and provides an analysis of the state of their self-assessment process. The result of their integrated planning and assessment process identifies Laboratory strengths and opportunities for improvement. Critical elements of that process are included in this report; namely, a high-level summary of external oversight activities, progress against Operations Improvement Initiatives, and a summary of Laboratory strengths and areas for improvement developed by management from across the laboratory. Some key areas targeted for improvement in FY2001 are: systems approach to resource management; information protection; integrated safety management flow-down to the benchtop; cost management; integrated assessment; Price Anderson Amendments Act (PAAA) Program; and travel risk mitigation.
The Fast Flux Test Facility (FFTF) is a liquid-metal-cooled fast reactor, which is designed to test a variety of different structural and fuel materials. A safety analysis is performed for each experiment that is irradiated in FFTF. The FFTF final safety analysis report (FSAR) assumed that all driver fuel assemblies would maintain cladding integrity during normal operations and all design transients. Maintenance of cladding integrity retains three barriers to any fission gas release to the public and also prevents any potential contact between the fuel and coolant. Experiments are, in general, expected to meet the same criterion. Selected experiments can, however, be classified as run-to-cladding-breach experiments (RTCB). The purpose of this paper is to describe alternative acceptance criteria for RTCB experiments that they feel provide protection equivalent to the maintenance of cladding integrity.
Web conferencing tools have entered the mainstream of business applications. Using web conferencing for IEEE conferences has a good potential of adding value to both organizers and participants. Authors propose a concept of Truly Integrated Conference (TIC) according to which a multi-point worldwide-distributed network of conference online authors/participants will enhance the standard (centralized) IEEE conference model, which requires attendance of the participants in person at the main conference location. The concept entails seamless integration of the onsite and online conference systems, including data/presentation, video, audio channels. Benefits and challenges of the TIC concept are analyzed. Requirements to the web conferencing system capable of supporting the TIC conference are presented and reviewed against commercial web conferencing tools. Case study of the IEEE Toronto International Conference ? Science and Technology for ...
Integrated Leakage Rate Test (ILRT) of containment system plays a very important role in safety of a Nuclear Power Plant. Containment system constitutes the last physical barrier to release of radioactivity from the core and is called upon to mitigate the consequences of not only accidents within the design basis, but also some of the highly unlikely severe accidents. Hence, leak tightness of containment becomes uttermost priority for the safety of plant personnel and public. The containment and associated ESFs are tested before the first criticality and there after periodically during service. The pre-operational integrated leakage rate is carried out at LOCA based design pressure, at periodic test pressure and at some intermediate pressure points to assess the leakage characteristics. This paper summarizes the various requirements and activities relevant to the ILRT of the Indian Pressurized Heavy Water Reactor (PHWR) containment system. ...
The conceptual design for MCP to be installed in the integral reactor SMART was carried out. Canned motor pump was adopted in the conceptual design of MCP. Three-dimensional modeling was performed to visualize the conceptual design of the MCP and to check interferences between the parts. The theoretical design procedure for the impeller was developed. The procedures for the flow field and structural analysis of impeller was also developed to assess the design validity and to verify its structural integrity. A computer program to analyze the dynamic characteristics of the rotor shaft of MCP was developed. The rotational speed sensor was designed and its performance test was conducted to verify the possibility of operation. A prototypes of the canned motor was manufactured and tested to confirm the validity of the design concept. The MCP design concept was also investigated for fabricability by establishing the manufacturing procedures. 41 refs., ...
We propose a combined migration velocity analysis and imaging method based on Kirchhoff integral migration and reverse time migration, using the residual curvature analysis and layer stripping strategy to build the velocity model. This method improves the image resolution of Kirchhoff integral migration and reduces the computations of the reverse time migration. It combines the advantages of efficiency and accuracy of the two migration methods. Its application in tunnel seismic prediction shows good results. Numerical experiments show that the imaging results of reverse time migration are better than the imaging results of Kirchhoff integral migration in many aspects of tunnel prediction. Field data show that this method has efficient computations and can establish a reasonable velocity model and a high quality imaging section. Combination with geological information can make an accurate prediction of the front of the ...
The overall aim is to develope programme modules for a Danish integrated computer system for projecting buildings and to ensure that Danish programmes are compatible with international systems. The programmes to be integrated here were AutoCAD and SCRIBE-Models and the simulation programme tsbi3 for the calculation of indoor climate conditions, power demand and energy consumption in buildings. For the calculation of humidity and temperature conditions in constructions a certain amount of integration is carried out with the MATCH programme (Moisture And Temperature calculations for Constructions of Hygroscopic materials). The activities embodied in the pilot projects are described in this working report in addition to individual computer programmes and prototypes, and examples and diagrams are presented. The programmes described deal with handling design elements in relation to AutoCAD designs, and the formatting of a ...
The integrity of the cables located in the target room is very important in the Fire PSA, because the CDF and CCDP are changed according to the results of a cable integrity that depends on the surrounding gas temperature. The conservative assumptions used in the Fire PSA typically specify that all of the equipment and cables of a room would fail when a fire happens in the room. But the realistic assessment of a fire risk by using a fire simulation tool has become necessary in the Fire PSA as described in the ANS Fire PRA Standard. This paper evaluates the cable integrity of eight pump rooms in the nuclear power plant by using the CFAST zone fire model. The upper layer gas temperature of each room is estimated, and an analysis based on the results of model simulations is used to judge the c...
Concrete is the primary material for building envelopes in some parts of the world, and its ability to store heat as well as its dynamic temperature changes will not only affect the deterioration rate of the exterior wall but will also greatly influence the energy efficiency of interior air conditioning. There are many methods for measuring the inner temperature of concrete, but they often have limitations, such as indirect estimation, cable installation requirements, high cost, or heterogeneity of the sample structure. In order to measure the internal temperature of concrete, this study integrated a Radio Frequency Integrated Circuit (RFIC) with a temperature sensor chip and embedded the device in concrete structures. A Smart Temperature Information Material (STIM) was thus developed. Thi...
In this work, we propose a conceptual distributed control framework for electrical grid integrated with distributed renewable energy generation systems in order to enable the development of the so-called ''smart electrical grid''. First, we introduce the key elements and their interactions in the proposed control architecture and discuss the design of the distributed control systems which are able to coordinate their actions to account for optimization considerations on the system operation. Subsequently, we focus on a specific wind/solar energy generation system connected to a reverse osmosis water desalination system and the electrical grid and design two supervisory predictive controllers via model predictive control to operate the integrated system taking into account short-term and lo...
The European Union (EU) market integration is leading to increasingly monopolistic electricity market infrastructures, which has opened a debate on the regulation of these so-called power exchanges. In this paper, we start by stating that there are two types of power exchanges in Europe, i.e. 'merchant' and 'cost-of-service regulated' power exchanges. We then discuss how regulation can be used to better align their incentives with the main power exchange tasks. We conclude that adopting the cost-of-service regulated model for all power exchanges in Europe could be counterproductive in the current context, but that regulation can help ensure that the benefits of the EU market integration materialize. Promising regulatory actions include tempering the reinforced market power of power exchanges, and quality-of-service regulation for the ongoing cooperation among power exchanges to organize trade across borders. - Research highlights: #-># ...
AbstractThe optimal definition of the size, shape and location of gross tumour volume is one of the most important steps in the planning of radiation therapy, and necessitates a proper...Full Text Available
Gas turbine component protection remains essential but environmental regulations emphasize the gas-turbine efficiency. Heating and cooling are features often integrated into the GI air intake packages to provide the optimum operating conditions. This paper describes modern inlet air treatment systems and discusses their effectiveness. (author)
The Brain Database is an information tool designed to aid in the integration of clinical and research results in neuroanatomy and regional biochemistry. It can handle a wide range of data types including...Full Text Available
The commercial potential of the 17 m, 100 kW Darrieus Vertical Axis Wind Turbine (VAWT) is discussed. Long term testing is in progress. As the first commercially adapted Darrieus turbines built by DOE, the superior cost, structural integrity, and output characteristics demonstrated by the prototypes appear particularly promising.
Neuroblasts born in the adult subventricular zone (SVZ) migrate long distances in the rostral migratory stream (RMS) to the olfactory bulbs where they integrate into circuitry as functional interneurons....Full Text Available
The accepted model of retroviral reverse transcription includes a circular DNA intermediate which requires strand displacement synthesis for linearization and creation of an integration-competent, long...Full Text Available
Perturbative approach to two-dimensional gravity and supergravity is considered. One-loop renormalization of the central charge of SL(2,R) Kac-Moody algebra is calculated perturbatively by functional integration and by explicit calculations of the Feynman diagrams. Also the wavefunction renormalization and the anomalous dimensions in the presence of gravity are calculated.
Digital information flow within the intensive care unit (ICU) continues to grow, with advances in technology and computational biology. Recent developments in the integration and archiving of these...Full Text Available
Quality reporting tools, integrated with ambulatory electronic health records, can help clinicians and administrators understand performance, manage populations, and improve quality. Report Central...Full Text Available
Using plasma glutathione S-transferase measurements hepatocellular integrity was assessed in groups of hyperthyroid and hypothyroid patients before and after treatment. Ten of 14 hyperthyroid patients...Full Text Available
Currently several groups are actively researching the integration of a megavoltage teletherapy unit with magnetic resonance (MR) imaging for real-time image-guided radiotherapy. The use of a...Full Text Available
The European demonstration project Hip Hip (house integrated photovoltaic high-tech in public) led the French photovoltaic market to reach a level comparable to those of the other European countries in terms of relevance of implemented solutions and costs. This document presents the best realizations and the experience gained through the Hip Hip project. Its aim is to convince the designers and managers of building projects of the advantages of photovoltaic installations integrated to the structure and connected to the power distribution grid: 1 - presentation of the technology; 2 - status of the Hip Hip demonstration project (goals, results: main innovations and impact on the French market); 3 - different possibilities of integration to the building structure; 4 - examples of projects realized in France in the framework of the Hip Hip project: integration in glass roof, frontage, added elements, ...
We examine how the structure and function of indirect flight muscle (IFM) and the entire flight system of Drosophila melanogaster are affected by phosphorylation of the myosin regulatory light chain...Full Text Available
The intercalated disk (ID) is a complex structure that electromechanically couples adjoining cardiac myocytes into a functional syncitium. The integrity of the disk is essential for normal cardiac function,...Full Text Available
When multiple integrals are approximately evaluated using Korobov cubature formulas, it is necessary to introduce a parameter characterizing the uniform distribution of the grid nodes. A new parameter for Korobov parallelepipedal grids is proposed, and an algorithm for its computation is described.
The long-term goal of nuclear transfer or alternative reprogramming approaches is to create patient-specific donor cells for transplantation therapy, avoiding immunorejection, a major complication in...Full Text Available
Integration of solid-state biosensors and lipid bilayer membranes is important for membrane protein research and drug discovery. In these sensors, it is critical that the solid-state sensing...Full Text Available
BackgroundMacrophages (MΦs) utilize macropinocytosis to integrate immune and metabolic signals in order to initiate an effective immune response. Diabetes is characterized...Full Text Available
We consider the integrable structure of the quantum lattice W_N algebras. We introduce the ultralocal Lax matrix, and show that the Yang-Baxter relation is satisfied with a Z_N invariant R-matrix. (orig.).
Four points are made about globalization and health. First, economic integration is a powerful force for raising the incomes of poor countries. In the past 20 years several large developing countries...Full Text Available
Experiments were conducted with a scale model of an amphibious vehicle towed in calm water. Pressure distribution along the submerged bow and hull bottom was determined for a range of vehicle speed. A waterjet system with a bow intake and an axial flow im...
Although integration generally is considered a critical step in the retrovirus life cycle, it has been reported that visna virus, which causes degenerative neurologic disease in sheep, can productively...Full Text Available
Life history theory proposes that early-life cues induce highly integrated responses in traits associated with energy partitioning, maturation, reproduction, and aging such that the individual phenotype...Full Text Available
This talk sketches out the main parameters (technology, circuit design) that have an influence on the hardening of digital CMOS integrated circuits. For each technology the more common defects are mentioned. General design rules are proposed to prevent or limit those defects. (D.L.). 2 refs., 2 figs.
We consider a symmetric boundary integral formulation associated with a mixed boundary value problem defined on a domain Omega is an element of the set of real numbers(2) with piecewise smooth boundary Gamma. We assume that Omega is mapped onto itself by ...
Few diseases exemplify the integration of research from bench to bedside as well as neonatal lupus, often referred to as a model of passively acquired autoimmunity. In essence, this disease encompasses...Full Text Available
We consider the spin-k/2 XXZ model in the antiferromagnetic regime using the free-field realization of the quantum affine algebra U_q(sl_2) of level k. We give a free-field realization of the type-II q-vertex operator, which describes creation and annihilation of physical particles in the model. By taking a trace of the type-I and type-II q-vertex operators over the irreducible highest-weight representation of U_q(sl_2), we also derive an integral formula for form factors in this model. Investigating the structure of poles, we obtain a residue formula for form factors, which is a lattice analog of the higher-spin extension of Smirnov's formula in the massive integrable quantum field theory. This result as well as the quantum deformation of the Knizhnik-Zamolodchikov equation for form factors shows a deep connection in the mathematical structure of the integrable lattice models and the massive integrable ...
Microelectronic failure analysis (FA) has been an integral part of the development of state-of-the-art integrated circuits. FA of MicroElectroMechanical Systems (MEMS) is moving from its infancy stage to assume an important role in the successful design, fabrication, performance and reliability analysis for this new technology. In previous work, the authors focused on the application of several techniques developed for integrated circuit analysis to an earlier version of a surface micromachined microengine fabricated at Sandia. Recently, they have identified important new failure modes in binary counters that incorporate a newer design of the microengine, using a subset of integrated circuit failure analysis techniques including optical microscopy, focused ion beam (FIB) techniques, atomic force microscopy (AFM), and scanning electron microscopy (SEM). The primary failure mode they have identified is ...
Explicit soliton solutions describing the inelastic interaction between sub- and superluminal particles are found within the framework of a new integrable model of relativistic classical field theory. The corresponding energies are nonnegative irrespective of the choice of reference frame.
BackgroundEscherichia coli strain EL350 contains chromosomally integrated phage lambda Red recombinase genes enabling this strain to be used for modifying the sequence...Full Text Available
An electroencephalogram (EEG) signal is extremely nonstationary, highly composite and very complex, all of which reflects the underlying integral neurodynamics. Understanding the EEG “grammar”,...Full Text Available
A numerical technique is developed for estimating water quality violation frequencies due to pollutant discharges from urban areas during combined sewer overflow events. The first four moments of in-stream pollutant concentration are found by integrating a pollutant loading - wat...
Alveolar enlargement, which is characteristic of bronchopulmonary dysplasia, congenital matrix disorders, and cigarette smoke-induced emphysema, is thought to result from enhanced inflammation and ensuing...Full Text Available
A balanced application of chemiosmosis to plant mitochondria is presented. Electron transport carriers are sequenced to achieve proton translocation across the inner membrane. The processes of proton translocation, use, and release are integrated to yield a description of the effective proton gradient available for metabolite transport and ATP synthesis. 27 references, 6 figures.
We study the integration of multisensory and central input at the level of an identified fly motoneuron, the ventral cervical nerve motoneuron (VCNM) cell, which controls head movements of the animal....Full Text Available
Average human behavior in cue combination tasks is well predicted by Bayesian inference models. As this capability is acquired over developmental timescales, the question arises, how it is learned....Full Text Available
BackgroundFormation of site specific genomic double strand breaks (DSBs), induced by the expression of a pair of engineered zinc-finger nucleases (ZFNs), dramatically increases the...Full Text Available
SummaryFor mammalian somatic cells the importance of microtubule cytoskeleton integrity in interphase cell cycle progression is uncertain. The loss, diminishment, or stabilization...Full Text Available
Rhodococcus fascians is a Gram-positive phytopathogen that induces shooty hyperplasia on its hosts through the secretion of cytokinins. Global transcriptomics using microarrays combined...Full Text Available
The national standards and state derivatives must be reinforced by models of curricular reform. In this paper, ARISE presents one model based on a set of principles--coherence, integration of the sciences, movement from concrete ideas to abstract ones, inquiry, connection and application, sequencing that is responsive to how people learn.
... relative suitability of areas for production of commercial timber was determined based on considerations of operability, proximity ... habitats, and an ecologically-based estimate of sustainable timber pr...
BackgroundThe biopharmaceutical industry requires cell lines to have an optimal proliferation rate and a high integral viable cell number resulting in a maximum volumetric recombinant...Full Text Available
The feasibility of producing a pocket-sized air contaminant monitor based upon a miniature gas chromatograph and an integrated microcomputer has been demonstrated and a practical design for the instrument is presented. This pocket-sized instrument is intended to be used as a pers...