Pile Load Capacity – Calculation Methods

Directory of Open Access Journals (Sweden)

Wrana Bogumił

2015-12-01

Full Text Available The article is a review of the current problems of the foundation pile capacity calculations. The article considers the main principles of pile capacity calculations presented in Eurocode 7 and other methods with adequate explanations. Two main methods are presented: α – method used to calculate the short-term load capacity of piles in cohesive soils and β – method used to calculate the long-term load capacity of piles in both cohesive and cohesionless soils. Moreover, methods based on cone CPTu result are presented as well as the pile capacity problem based on static tests.

An improved method for calculation of interface pressure force in PLIC-VOF methods

International Nuclear Information System (INIS)

Sefollahi, M.; Shirani, E.

2004-08-01

Conventional methods for the modeling of surface tension force in Piecewise Linear Interface Calculation-Volume of Fluid (PLIC-VOF) methods, such as Continuum Surface Force (CSF), Continuum Surface Stress (CSS) and also Meier's method, convert the surface tension force into a body force. Not only do they include the force in the interfacial cells but also in the neighboring cells. Thus they produce spurious currents. Also the pressure jump, due to the surface tension, is not calculated accurately in these methods. In this paper a more accurate method for the application of interface force in the computational modeling of free surfaces and interfaces which use PLIC-VOF methods is developed. This method is based on the evaluation of the surface tension force only in the interfacial cells and not the neighboring cells. Also the normal and the interface surface area needed for the calculation of the surface tension force is calculated more accurately. The present method is applied to a two-dimensional motionless drop of liquid and a bubble of gas as well as a non-circular two-dimensional drop, which oscillates due to the surface tension force, in an initially stagnant fluid with no gravity force. The results are compared with the results of the cases when CSF, CSS and Meier's methods are used. It is shown that the present method calculates pressure jump at the interface more accurately and produces less spurious currents comparing to CSS an CSF models. (author)

Acceleration methods and models in Sn calculations

International Nuclear Information System (INIS)

Sbaffoni, M.M.; Abbate, M.J.

1984-01-01

In some neutron transport problems solved by the discrete ordinate method, it is relatively common to observe some particularities as, for example, negative fluxes generation, slow and insecure convergences and solution instabilities. The commonly used models for neutron flux calculation and acceleration methods included in the most used codes were analyzed, in face of their use in problems characterized by a strong upscattering effect. Some special conclusions derived from this analysis are presented as well as a new method to perform the upscattering scaling for solving the before mentioned problems in this kind of cases. This method has been included in the DOT3.5 code (two dimensional discrete ordinates radiation transport code) generating a new version of wider application. (Author) [es

New method for calculation of integral characteristics of thermal plumes

DEFF Research Database (Denmark)

Zukowska, Daria; Popiolek, Zbigniew; Melikov, Arsen Krikor

2008-01-01

A method for calculation of integral characteristics of thermal plumes is proposed. The method allows for determination of the integral parameters of plumes based on speed measurements performed with omnidirectional low velocity thermoanemometers. The method includes a procedure for calculation...... of the directional velocity (upward component of the mean velocity). The method is applied for determination of the characteristics of an asymmetric thermal plume generated by a sitting person. The method was validated in full-scale experiments in a climatic chamber with a thermal manikin as a simulator of a sitting...

Pressure algorithm for elliptic flow calculations with the PDF method

Science.gov (United States)

Anand, M. S.; Pope, S. B.; Mongia, H. C.

1991-01-01

An algorithm to determine the mean pressure field for elliptic flow calculations with the probability density function (PDF) method is developed and applied. The PDF method is a most promising approach for the computation of turbulent reacting flows. Previous computations of elliptic flows with the method were in conjunction with conventional finite volume based calculations that provided the mean pressure field. The algorithm developed and described here permits the mean pressure field to be determined within the PDF calculations. The PDF method incorporating the pressure algorithm is applied to the flow past a backward-facing step. The results are in good agreement with data for the reattachment length, mean velocities, and turbulence quantities including triple correlations.

Assessment of seismic margin calculation methods

International Nuclear Information System (INIS)

Kennedy, R.P.; Murray, R.C.; Ravindra, M.K.; Reed, J.W.; Stevenson, J.D.

1989-03-01

Seismic margin review of nuclear power plants requires that the High Confidence of Low Probability of Failure (HCLPF) capacity be calculated for certain components. The candidate methods for calculating the HCLPF capacity as recommended by the Expert Panel on Quantification of Seismic Margins are the Conservative Deterministic Failure Margin (CDFM) method and the Fragility Analysis (FA) method. The present study evaluated these two methods using some representative components in order to provide further guidance in conducting seismic margin reviews. It is concluded that either of the two methods could be used for calculating HCLPF capacities. 21 refs., 9 figs., 6 tabs

Statistical methods for including two-body forces in large system calculations

International Nuclear Information System (INIS)

Grimes, S.M.

1980-07-01

Large systems of interacting particles are often treated by assuming that the effect on any one particle of the remaining N-1 may be approximated by an average potential. This approach reduces the problem to that of finding the bound-state solutions for a particle in a potential; statistical mechanics is then used to obtain the properties of the many-body system. In some physical systems this approach may not be acceptable, because the two-body force component cannot be treated in this one-body limit. A technique for incorporating two-body forces in such calculations in a more realistic fashion is described. 1 figure

ICECON: a computer program used to calculate containment back pressure for LOCA analysis (including ice condenser plants)

International Nuclear Information System (INIS)

1976-07-01

The ICECON computer code provides a method for conservatively calculating the long term back pressure transient in the containment resulting from a hypothetical Loss-of-Coolant Accident (LOCA) for PWR plants including ice condenser containment systems. The ICECON computer code was developed from the CONTEMPT/LT-022 code. A brief discussion of the salient features of a typical ice condenser containment is presented. Details of the ice condenser models are explained. The corrections and improvements made to CONTEMPT/LT-022 are included. The organization of the code, including the calculational procedure, is outlined. The user's manual, to be used in conjunction with the CONTEMPT/LT-022 user's manual, a sample problem, a time-step study (solution convergence) and a comparison of ICECON results with the results of the NSSS vendor are presented. In general, containment pressure calculated with the ICECON code agree with those calculated by the NSSS vendor using the same mass and energy release rates to the containment

Structural system reliability calculation using a probabilistic fault tree analysis method

Science.gov (United States)

Torng, T. Y.; Wu, Y.-T.; Millwater, H. R.

1992-01-01

The development of a new probabilistic fault tree analysis (PFTA) method for calculating structural system reliability is summarized. The proposed PFTA procedure includes: developing a fault tree to represent the complex structural system, constructing an approximation function for each bottom event, determining a dominant sampling sequence for all bottom events, and calculating the system reliability using an adaptive importance sampling method. PFTA is suitable for complicated structural problems that require computer-intensive computer calculations. A computer program has been developed to implement the PFTA.

Sputtering calculations with the discrete ordinated method

International Nuclear Information System (INIS)

Hoffman, T.J.; Dodds, H.L. Jr.; Robinson, M.T.; Holmes, D.K.

1977-01-01

The purpose of this work is to investigate the applicability of the discrete ordinates (S/sub N/) method to light ion sputtering problems. In particular, the neutral particle discrete ordinates computer code, ANISN, was used to calculate sputtering yields. No modifications to this code were necessary to treat charged particle transport. However, a cross section processing code was written for the generation of multigroup cross sections; these cross sections include a modification to the total macroscopic cross section to account for electronic interactions and small-scattering-angle elastic interactions. The discrete ordinates approach enables calculation of the sputtering yield as functions of incident energy and angle and of many related quantities such as ion reflection coefficients, angular and energy distributions of sputtering particles, the behavior of beams penetrating thin foils, etc. The results of several sputtering problems as calculated with ANISN are presented

Comparative Study of the Volumetric Methods Calculation Using GNSS Measurements

Science.gov (United States)

Şmuleac, Adrian; Nemeş, Iacob; Alina Creţan, Ioana; Sorina Nemeş, Nicoleta; Şmuleac, Laura

2017-10-01

This paper aims to achieve volumetric calculations for different mineral aggregates using different methods of analysis and also comparison of results. To achieve these comparative studies and presentation were chosen two software licensed, namely TopoLT 11.2 and Surfer 13. TopoLT program is a program dedicated to the development of topographic and cadastral plans. 3D terrain model, level courves and calculation of cut and fill volumes, including georeferencing of images. The program Surfer 13 is produced by Golden Software, in 1983 and is active mainly used in various fields such as agriculture, construction, geophysical, geotechnical engineering, GIS, water resources and others. It is also able to achieve GRID terrain model, to achieve the density maps using the method of isolines, volumetric calculations, 3D maps. Also, it can read different file types, including SHP, DXF and XLSX. In these paper it is presented a comparison in terms of achieving volumetric calculations using TopoLT program by two methods: a method where we choose a 3D model both for surface as well as below the top surface and a 3D model in which we choose a 3D terrain model for the bottom surface and another 3D model for the top surface. The comparison of the two variants will be made with data obtained from the realization of volumetric calculations with the program Surfer 13 generating GRID terrain model. The topographical measurements were performed with equipment from Leica GPS 1200 Series. Measurements were made using Romanian position determination system - ROMPOS which ensures accurate positioning of reference and coordinates ETRS through the National Network of GNSS Permanent Stations. GPS data processing was performed with the program Leica Geo Combined Office. For the volumetric calculating the GPS used point are in 1970 stereographic projection system and for the altitude the reference is 1975 the Black Sea projection system.

Methods for magnetostatic field calculation

International Nuclear Information System (INIS)

Vorozhtsov, S.B.

1984-01-01

Two methods for magnetostatic field calculation: differential and integrat are considered. Both approaches are shown to have certain merits and drawbacks, choice of the method depend on the type of the solved problem. An opportunity of combination of these tWo methods in one algorithm (hybrid method) is considered

Soil structure interaction calculations: a comparison of methods

Energy Technology Data Exchange (ETDEWEB)

Wight, L.; Zaslawsky, M.

1976-07-22

Two approaches for calculating soil structure interaction (SSI) are compared: finite element and lumped mass. Results indicate that the calculations with the lumped mass method are generally conservative compared to those obtained by the finite element method. They also suggest that a closer agreement between the two sets of calculations is possible, depending on the use of frequency-dependent soil springs and dashpots in the lumped mass calculations. There is a total lack of suitable guidelines for implementing the lumped mass method of calculating SSI, which leads to the conclusion that the finite element method is generally superior for calculative purposes.

Soil structure interaction calculations: a comparison of methods

International Nuclear Information System (INIS)

Wight, L.; Zaslawsky, M.

1976-01-01

Two approaches for calculating soil structure interaction (SSI) are compared: finite element and lumped mass. Results indicate that the calculations with the lumped mass method are generally conservative compared to those obtained by the finite element method. They also suggest that a closer agreement between the two sets of calculations is possible, depending on the use of frequency-dependent soil springs and dashpots in the lumped mass calculations. There is a total lack of suitable guidelines for implementing the lumped mass method of calculating SSI, which leads to the conclusion that the finite element method is generally superior for calculative purposes

Reactor perturbation calculations by Monte Carlo methods

International Nuclear Information System (INIS)

Gubbins, M.E.

1965-09-01

Whilst Monte Carlo methods are useful for reactor calculations involving complicated geometry, it is difficult to apply them to the calculation of perturbation worths because of the large amount of computing time needed to obtain good accuracy. Various ways of overcoming these difficulties are investigated in this report, with the problem of estimating absorbing control rod worths particularly in mind. As a basis for discussion a method of carrying out multigroup reactor calculations by Monte Carlo methods is described. Two methods of estimating a perturbation worth directly, without differencing two quantities of like magnitude, are examined closely but are passed over in favour of a third method based on a correlation technique. This correlation method is described, and demonstrated by a limited range of calculations for absorbing control rods in a fast reactor. In these calculations control rod worths of between 1% and 7% in reactivity are estimated to an accuracy better than 10% (3 standard errors) in about one hour's computing time on the English Electric KDF.9 digital computer. (author)

Development and application of advanced methods for electronic structure calculations

DEFF Research Database (Denmark)

Schmidt, Per Simmendefeldt

. For this reason, part of this thesis relates to developing and applying a new method for constructing so-called norm-conserving PAW setups, that are applicable to GW calculations by using a genetic algorithm. The effect of applying the new setups significantly affects the absolute band positions, both for bulk......This thesis relates to improvements and applications of beyond-DFT methods for electronic structure calculations that are applied in computational material science. The improvements are of both technical and principal character. The well-known GW approximation is optimized for accurate calculations...... of electronic excitations in two-dimensional materials by exploiting exact limits of the screened Coulomb potential. This approach reduces the computational time by an order of magnitude, enabling large scale applications. The GW method is further improved by including so-called vertex corrections. This turns...

Calculation methods for determining dose equivalent

International Nuclear Information System (INIS)

Endres, G.W.R.; Tanner, J.E.; Scherpelz, R.I.; Hadlock, D.E.

1987-11-01

A series of calculations of neutron fluence as a function of energy in an anthropomorphic phantom was performed to develop a system for determining effective dose equivalent for external radiation sources. Critical organ dose equivalents are calculated and effective dose equivalents are determined using ICRP-26 [1] methods. Quality factors based on both present definitions and ICRP-40 definitions are used in the analysis. The results of these calculations are presented and discussed. The effective dose equivalent determined using ICRP-26 methods is significantly smaller than the dose equivalent determined by traditional methods. No existing personnel dosimeter or health physics instrument can determine effective dose equivalent. At the present time, the conversion of dosimeter response to dose equivalent is based on calculations for maximal or ''cap'' values using homogeneous spherical or cylindrical phantoms. The evaluated dose equivalent is, therefore, a poor approximation of the effective dose equivalent as defined by ICRP Publication 26. 3 refs., 2 figs., 1 tab

Development of 3-D FBR heterogeneous core calculation method based on characteristics method

International Nuclear Information System (INIS)

Takeda, Toshikazu; Maruyama, Manabu; Hamada, Yuzuru; Nishi, Hiroshi; Ishibashi, Junichi; Kitano, Akihiro

2002-01-01

A new 3-D transport calculation method taking into account the heterogeneity of fuel assemblies has been developed by combining the characteristics method and the nodal transport method. In the axial direction the nodal transport method is applied, and the characteristics method is applied to take into account the radial heterogeneity of fuel assemblies. The numerical calculations have been performed to verify 2-D radial calculations of FBR assemblies and partial core calculations. Results are compared with the reference Monte-Carlo calculations. A good agreement has been achieved. It is shown that the present method has an advantage in calculating reaction rates in a small region

Statistic method of research reactors maximum permissible power calculation

International Nuclear Information System (INIS)

Grosheva, N.A.; Kirsanov, G.A.; Konoplev, K.A.; Chmshkyan, D.V.

1998-01-01

The technique for calculating maximum permissible power of a research reactor at which the probability of the thermal-process accident does not exceed the specified value, is presented. The statistical method is used for the calculations. It is regarded that the determining function related to the reactor safety is the known function of the reactor power and many statistically independent values which list includes the reactor process parameters, geometrical characteristics of the reactor core and fuel elements, as well as random factors connected with the reactor specific features. Heat flux density or temperature is taken as a limiting factor. The program realization of the method discussed is briefly described. The results of calculating the PIK reactor margin coefficients for different probabilities of the thermal-process accident are considered as an example. It is shown that the probability of an accident with fuel element melting in hot zone is lower than 10 -8 1 per year for the reactor rated power [ru

Radiation transport calculation methods in BNCT

International Nuclear Information System (INIS)

Koivunoro, H.; Seppaelae, T.; Savolainen, S.

2000-01-01

Boron neutron capture therapy (BNCT) is used as a radiotherapy for malignant brain tumours. Radiation dose distribution is necessary to determine individually for each patient. Radiation transport and dose distribution calculations in BNCT are more complicated than in conventional radiotherapy. Total dose in BNCT consists of several different dose components. The most important dose component for tumour control is therapeutic boron dose D B . The other dose components are gamma dose D g , incident fast neutron dose D f ast n and nitrogen dose D N . Total dose is a weighted sum of the dose components. Calculation of neutron and photon flux is a complex problem and requires numerical methods, i.e. deterministic or stochastic simulation methods. Deterministic methods are based on the numerical solution of Boltzmann transport equation. Such are discrete ordinates (SN) and spherical harmonics (PN) methods. The stochastic simulation method for calculation of radiation transport is known as Monte Carlo method. In the deterministic methods the spatial geometry is partitioned into mesh elements. In SN method angular integrals of the transport equation are replaced with weighted sums over a set of discrete angular directions. Flux is calculated iteratively for all these mesh elements and for each discrete direction. Discrete ordinates transport codes used in the dosimetric calculations are ANISN, DORT and TORT. In PN method a Legendre expansion for angular flux is used instead of discrete direction fluxes, land the angular dependency comes a property of vector function space itself. Thus, only spatial iterations are required for resulting equations. A novel radiation transport code based on PN method and tree-multigrid technique (TMG) has been developed at VTT (Technical Research Centre of Finland). Monte Carlo method solves the radiation transport by randomly selecting neutrons and photons from a prespecified boundary source and following the histories of selected particles

Quantum mechanical methods for calculation of force constants

International Nuclear Information System (INIS)

Mullally, D.J.

1985-01-01

The focus of this thesis is upon the calculation of force constants; i.e., the second derivatives of the potential energy with respect to nuclear displacements. This information is useful for the calculation of molecular vibrational modes and frequencies. In addition, it may be used for the location and characterization of equilibrium and transition state geometries. The methods presented may also be applied to the calculation of electric polarizabilities and infrared and Raman vibrational intensities. Two approaches to this problem are studied and evaluated: finite difference methods and analytical techniques. The most suitable method depends on the type and level of theory used to calculate the electronic wave function. Double point displacement finite differencing is often required for accurate calculation of the force constant matrix. These calculations require energy and gradient calculations on both sides of the geometry of interest. In order to speed up these calculations, a novel method is presented that uses geometry dependent information about the wavefunction. A detailed derivation for the analytical evaluation of force constants with a complete active space multiconfiguration self consistent field wave function is presented

Willow growing - Methods of calculation and profitability

International Nuclear Information System (INIS)

Rosenqvist, H.

1997-01-01

The calculation method presented here makes it possible to conduct profitability comparisons between annual and perennial crops and in addition take the planning situation into account. The method applied is a modified total step calculation. The difference between a traditional total step calculation and the modified version is the way in which payments and disbursements are taken into account over a period of several years. This is achieved by combining the present value method and the annuity method. The choice of interest rate has great bearing on the result in perennial calculations. The various components influencing the interest rate are analysed and factors relating to the establishment of the interest rate in different situations are described. The risk factor can be an important variable component of the interest rate calculation. Risk is also addressed from an approach in accordance with portfolio theory. The application of the methods sheds light on the profitability of Salix cultivation from the viewpoint of business economics, and also how different factors influence the profitability of Salix cultivation. Aspects studied are harvesting intervals, the importance of yield level, the competitiveness of Salix versus grain cultivation, the influence of income taxes on profitability etc. Methods for evaluation of activities concerning cultivation of a perennial crop are described and also involve the application of nitrogen fertilization to Salix cultivation. Studies have been performed using these methods to look into nitrogen fertilizer profitability in Salix cultivation during the first rotation period. Nitrogen fertilizer profitability has been investigated involving both production functions and cost calculations, taking the year fertilization into consideration. 72 refs., 2 figs., 52 tabs

Lambda-guided calculation method (LGC method) for xenon/CT CBF

Energy Technology Data Exchange (ETDEWEB)

Sase, Shigeru [Anzai Medical Co., Ltd., Tokyo (Japan); Honda, Mitsuru; Kushida, Tsuyoshi; Seiki, Yoshikatsu; Machida, Keiichi; Shibata, Iekado [Toho Univ., Tokyo (Japan). School of Medicine

2001-12-01

A quantitative CBF calculation method for xenon/CT was developed by logically estimating time-course change rate (rate constant) of arterial xenon concentration from that of end-tidal xenon concentration. A single factor ({gamma}) was introduced to correlate the end-tidal rate constant (Ke) with the arterial rate constant (Ka) in a simplified equation. This factor ({gamma}) is thought to reflect the diffusing capacity of the lung for xenon. When an appropriate value is given to {gamma}, it is possible to calculate the arterial rate constant (Calculated Ka) from Ke. To determine {gamma} for each xenon/CT CBF examination, a procedure was established which utilizes the characteristics of white matter lambda; lambda refers to xenon brain-blood partition coefficient. Xenon/CT studies were performed on four healthy volunteers. Hemispheric CBF values (47.0{+-}9.0 ml/100 g/min) with use of Calculated Ka were close to the reported normative values. For a 27-year-old healthy man, the rate constant for the common carotid artery was successfully measured and nearly equal to Calculated Ka. The authors conclude the method proposed in this work, lambda-guided calculation method, could make xenon/CT CBF substantially reliable and quantitative by effective use of end-tidal xenon. (author)

Lambda-guided calculation method (LGC method) for xenon/CT CBF

International Nuclear Information System (INIS)

Sase, Shigeru; Honda, Mitsuru; Kushida, Tsuyoshi; Seiki, Yoshikatsu; Machida, Keiichi; Shibata, Iekado

2001-01-01

A quantitative CBF calculation method for xenon/CT was developed by logically estimating time-course change rate (rate constant) of arterial xenon concentration from that of end-tidal xenon concentration. A single factor (γ) was introduced to correlate the end-tidal rate constant (Ke) with the arterial rate constant (Ka) in a simplified equation. This factor (γ) is thought to reflect the diffusing capacity of the lung for xenon. When an appropriate value is given to γ, it is possible to calculate the arterial rate constant (Calculated Ka) from Ke. To determine γ for each xenon/CT CBF examination, a procedure was established which utilizes the characteristics of white matter lambda; lambda refers to xenon brain-blood partition coefficient. Xenon/CT studies were performed on four healthy volunteers. Hemispheric CBF values (47.0±9.0 ml/100 g/min) with use of Calculated Ka were close to the reported normative values. For a 27-year-old healthy man, the rate constant for the common carotid artery was successfully measured and nearly equal to Calculated Ka. The authors conclude the method proposed in this work, lambda-guided calculation method, could make xenon/CT CBF substantially reliable and quantitative by effective use of end-tidal xenon. (author)

An assessment of methods of calculating sodium-voiding reactivity in plutonium-fuelled fast reactors

International Nuclear Information System (INIS)

Butland, A.T.D.; Simmons, W.N.; Stevenson, J.M.

1980-01-01

After a survey of the requirements an assessment of the accuracy of calculations of the sodium-void effect using UK methods and data is made on the basis of the following work: (a) The analysis of small and large sodium voids in the MOZART and Zebra 13 small (300 MW(e)) fast reactor mock-ups and the BIZET large fast reactor mock-ups, all of conventional design. The analysis was carried out using the UK FGL5 fine group nuclear data library, the MURAL cell code, whole reactor diffusion theory calculations of the neutron flux and perturbation theory methods. Exact perturbation theory was used in many cases, otherwise first-order perturbation theory calculations were adjusted to give results equivalent to exact perturbation theory. (b) Theoretical studies of some effects, including the following: (i) The effects of extrapolating to fuel operating temperature; (ii) Fuel-cycle and burnup effects, including the gradual replacement through a fuel cycle of control-rod absorption by fission product absorption, the loss of fissile material and the change in fuel nuclide relative composition; (iii) The heterogeneity effects of large fuelled subassemblies in pin geometry. (c) Theoretical studies of approximations in the calculational methods, including the following: (i) The importance in the whole reactor calculation of the energy group structure and the spatial mesh, including comparisons of calculations in two (RZ) and three-dimensional geometry; (ii) The importance of reactor material boundaries in the calculation of resonance shielding effects; (iii) The use of neutron fluxes calculated using neutron diffusion theory rather than transport theory. (author)

Manual method for dose calculation in gynecologic brachytherapy

International Nuclear Information System (INIS)

Vianello, Elizabeth A.; Almeida, Carlos E. de; Biaggio, Maria F. de

1998-01-01

This paper describes a manual method for dose calculation in brachytherapy of gynecological tumors, which allows the calculation of the doses at any plane or point of clinical interest. This method uses basic principles of vectorial algebra and the simulating orthogonal films taken from the patient with the applicators and dummy sources in place. The results obtained with method were compared with the values calculated with the values calculated with the treatment planning system model Theraplan and the agreement was better than 5% in most cases. The critical points associated with the final accuracy of the proposed method is related to the quality of the image and the appropriate selection of the magnification factors. This method is strongly recommended to the radiation oncology centers where are no treatment planning systems available and the dose calculations are manually done. (author)

Advanced numerical methods for three dimensional two-phase flow calculations

Energy Technology Data Exchange (ETDEWEB)

Toumi, I. [Laboratoire d`Etudes Thermiques des Reacteurs, Gif sur Yvette (France); Caruge, D. [Institut de Protection et de Surete Nucleaire, Fontenay aux Roses (France)

1997-07-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses an extension of Roe`s method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations.

Advanced numerical methods for three dimensional two-phase flow calculations

International Nuclear Information System (INIS)

Toumi, I.; Caruge, D.

1997-01-01

This paper is devoted to new numerical methods developed for both one and three dimensional two-phase flow calculations. These methods are finite volume numerical methods and are based on the use of Approximate Riemann Solvers concepts to define convective fluxes versus mean cell quantities. The first part of the paper presents the numerical method for a one dimensional hyperbolic two-fluid model including differential terms as added mass and interface pressure. This numerical solution scheme makes use of the Riemann problem solution to define backward and forward differencing to approximate spatial derivatives. The construction of this approximate Riemann solver uses an extension of Roe's method that has been successfully used to solve gas dynamic equations. As far as the two-fluid model is hyperbolic, this numerical method seems very efficient for the numerical solution of two-phase flow problems. The scheme was applied both to shock tube problems and to standard tests for two-fluid computer codes. The second part describes the numerical method in the three dimensional case. The authors discuss also some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. Such a scheme is not implemented in a thermal-hydraulic computer code devoted to 3-D steady-state and transient computations. Some results obtained for Pressurised Water Reactors concerning upper plenum calculations and a steady state flow in the core with rod bow effect evaluation are presented. In practice these new numerical methods have proved to be stable on non staggered grids and capable of generating accurate non oscillating solutions for two-phase flow calculations

A comparison study of size-specific dose estimate calculation methods

Energy Technology Data Exchange (ETDEWEB)

Parikh, Roshni A. [Rainbow Babies and Children' s Hospital, University Hospitals Cleveland Medical Center, Case Western Reserve University School of Medicine, Department of Radiology, Cleveland, OH (United States); University of Michigan Health System, Department of Radiology, Ann Arbor, MI (United States); Wien, Michael A.; Jordan, David W.; Ciancibello, Leslie; Berlin, Sheila C. [Rainbow Babies and Children' s Hospital, University Hospitals Cleveland Medical Center, Case Western Reserve University School of Medicine, Department of Radiology, Cleveland, OH (United States); Novak, Ronald D. [Rainbow Babies and Children' s Hospital, University Hospitals Cleveland Medical Center, Case Western Reserve University School of Medicine, Department of Radiology, Cleveland, OH (United States); Rebecca D. Considine Research Institute, Children' s Hospital Medical Center of Akron, Center for Mitochondrial Medicine Research, Akron, OH (United States); Klahr, Paul [CT Clinical Science, Philips Healthcare, Highland Heights, OH (United States); Soriano, Stephanie [Rainbow Babies and Children' s Hospital, University Hospitals Cleveland Medical Center, Case Western Reserve University School of Medicine, Department of Radiology, Cleveland, OH (United States); University of Washington, Department of Radiology, Seattle, WA (United States)

2018-01-15

The size-specific dose estimate (SSDE) has emerged as an improved metric for use by medical physicists and radiologists for estimating individual patient dose. Several methods of calculating SSDE have been described, ranging from patient thickness or attenuation-based (automated and manual) measurements to weight-based techniques. To compare the accuracy of thickness vs. weight measurement of body size to allow for the calculation of the size-specific dose estimate (SSDE) in pediatric body CT. We retrospectively identified 109 pediatric body CT examinations for SSDE calculation. We examined two automated methods measuring a series of level-specific diameters of the patient's body: method A used the effective diameter and method B used the water-equivalent diameter. Two manual methods measured patient diameter at two predetermined levels: the superior endplate of L2, where body width is typically most thin, and the superior femoral head or iliac crest (for scans that did not include the pelvis), where body width is typically most thick; method C averaged lateral measurements at these two levels from the CT projection scan, and method D averaged lateral and anteroposterior measurements at the same two levels from the axial CT images. Finally, we used body weight to characterize patient size, method E, and compared this with the various other measurement methods. Methods were compared across the entire population as well as by subgroup based on body width. Concordance correlation (ρ{sub c}) between each of the SSDE calculation methods (methods A-E) was greater than 0.92 across the entire population, although the range was wider when analyzed by subgroup (0.42-0.99). When we compared each SSDE measurement method with CTDI{sub vol,} there was poor correlation, ρ{sub c}<0.77, with percentage differences between 20.8% and 51.0%. Automated computer algorithms are accurate and efficient in the calculation of SSDE. Manual methods based on patient thickness provide

Calculation Method for the Projection of Future Spent Nuclear Fuel Discharges

International Nuclear Information System (INIS)

B. McLeod

2002-01-01

This report describes the calculation method developed for the projection of future utility spent nuclear fuel (SNF) discharges in regard to their timing, quantity, burnup, and initial enrichment. This projection method complements the utility-supplied RW-859 data on historic discharges and short-term projections of SNF discharges by providing long-term projections that complete the total life cycle of discharges for each of the current U.S. nuclear power reactors. The method was initially developed in mid-1999 to update the SNF discharge projection associated with the 1995 RW-859 utility survey (CRWMS M and O 1996). and was further developed as described in Rev. 00 of this report (CRWMS M and O 2001a). Primary input to the projection of SNF discharges is the utility projection of the next five discharges from each nuclear unit, which is provided via the revised final version of the Energy Information Administration (EIA) 1998 RW-859 utility survey (EIA 2000a). The projection calculation method is implemented via a set of Excel 97 spreadsheets. These calculations provide the interface between receipt of the utility five-discharge projections that are provided in the RW-859 survey, and the delivery of projected life-cycle SNF discharge quantities and characteristics in the format requisite for performing logistics analysis to support design of the Civilian Radioactive Waste Management System (CRWMS). Calculation method improvements described in this report include the addition of a reactor-specific maximum enrichment-based discharge burnup limit. This limit is the consequence of the enrichment limit, currently 5 percent. which is imposed as a Nuclear Regulatory Commission (NRC) license condition on nuclear fuel fabrication plants. In addition, the calculation method now includes the capability for projecting future nuclear plant power upratings, consistent with many such recent plant uprates and the prospect of additional future uprates. Finally. this report

Calculation methods for determining dose equivalent

International Nuclear Information System (INIS)

Endres, G.W.R.; Tanner, J.E.; Scherpelz, R.I.; Hadlock, D.E.

1988-01-01

A series of calculations of neutron fluence as a function of energy in an anthropomorphic phantom was performed to develop a system for determining effective dose equivalent for external radiation sources. critical organ dose equivalents are calculated and effective dose equivalents are determined using ICRP-26 methods. Quality factors based on both present definitions and ICRP-40 definitions are used in the analysis. The results of these calculations are presented and discussed

Range calculations using multigroup transport methods

International Nuclear Information System (INIS)

Hoffman, T.J.; Robinson, M.T.; Dodds, H.L. Jr.

1979-01-01

Several aspects of radiation damage effects in fusion reactor neutron and ion irradiation environments are amenable to treatment by transport theory methods. In this paper, multigroup transport techniques are developed for the calculation of particle range distributions. These techniques are illustrated by analysis of Au-196 atoms recoiling from (n,2n) reactions with gold. The results of these calculations agree very well with range calculations performed with the atomistic code MARLOWE. Although some detail of the atomistic model is lost in the multigroup transport calculations, the improved computational speed should prove useful in the solution of fusion material design problems

Calculation of NMR chemical shifts. 7. Gauge-invariant INDO method

Science.gov (United States)

Fukui, H.; Miura, K.; Hirai, A.

A gauge-invariant INDO method based on the coupled Hartree-Fuck perturbation theory is presented and applied to the calculation of 1H and 13C chemical shifts of hydrocarbons including ring compounds. Invariance of the diamagnetic and paramagnetic shieldings with respect to displacement of the coordinate origin is discussed. Comparison between calculated and experimental results exhibits fairly good agreement, provided that the INDO parameters of Ellis et al. (J. Am. Chem. Soc.94, 4069 (1972)) are used with the inclusion of all multicenter one-electron integrals.

An assessment of methods of calculating sodium voiding reactivity in plutonium fuelled fast reactors

International Nuclear Information System (INIS)

Butland, A.T.D.; Simmons, W.N.; Stevenson, J.M.

1979-01-01

After a survey of the requirements an assessment of the accuracy of calculations of the sodium void effect using UK methods and data is made on the basis of the following work. First, the analysis of small and large sodium voids in the MOZART and Zebra 13 small (300 MW(E)) fast reactor mock-ups and the BIZET large fast reactor mock-ups, all of conventional design. The analysis was carried out using the UK FGL5 fine group nuclear data library, the MURAL cell code, whole reactor diffusion theory calculations of the neutron flux and perturbation theory methods. Exact perturbation theory was used in many cases, otherwise first order perturbation theory calculations were adjusted to give results equivalent to exact perturbation theory. Second, theoretical studies of some effects, including, the effects of extrapolating to fuel operating temperatures, fuel cycle and burn-up effects, and the heterogeneity effects of large fuelled subassemblies in pin geometry. Third, theoretical studies of approximations in the calculational methods including, the importance in the whole reactor calculation of the energy group structure and the spatial mesh, the importance of reactor material boundaries in the calculation of resonance shielding effects, and the use of neutron fluxes calculated using neutron diffusion theory rather than transport theory. (U.K.)

Non-iterative method to calculate the periodical distribution of temperature in reactors with thermal regeneration

International Nuclear Information System (INIS)

Sanchez de Alsina, O.L.; Scaricabarozzi, R.A.

1982-01-01

A matrix non-iterative method to calculate the periodical distribution in reactors with thermal regeneration is presented. In case of exothermic reaction, a source term will be included. A computer code was developed to calculate the final temperature distribution in solids and in the outlet temperatures of the gases. The results obtained from ethane oxidation calculation in air, using the Dietrich kinetic data are presented. This method is more advantageous than iterative methods. (E.G.) [pt

Acceleration methods for assembly-level transport calculations

International Nuclear Information System (INIS)

Adams, Marvin L.; Ramone, Gilles

1995-01-01

A family acceleration methods for the iterations that arise in assembly-level transport calculations is presented. A single iteration in these schemes consists of a transport sweep followed by a low-order calculation which is itself a simplified transport problem. It is shown that a previously-proposed method fitting this description is unstable in two and three dimensions. It is presented a family of methods and shown that some members are unconditionally stable. (author). 8 refs, 4 figs, 4 tabs

The New Performance Calculation Method of Fouled Axial Flow Compressor

Directory of Open Access Journals (Sweden)

Huadong Yang

2014-01-01

Full Text Available Fouling is the most important performance degradation factor, so it is necessary to accurately predict the effect of fouling on engine performance. In the previous research, it is very difficult to accurately model the fouled axial flow compressor. This paper develops a new performance calculation method of fouled multistage axial flow compressor based on experiment result and operating data. For multistage compressor, the whole compressor is decomposed into two sections. The first section includes the first 50% stages which reflect the fouling level, and the second section includes the last 50% stages which are viewed as the clean stage because of less deposits. In this model, the performance of the first section is obtained by combining scaling law method and linear progression model with traditional stage stacking method; simultaneously ambient conditions and engine configurations are considered. On the other hand, the performance of the second section is calculated by averaged infinitesimal stage method which is based on Reynolds’ law of similarity. Finally, the model is successfully applied to predict the 8-stage axial flow compressor and 16-stage LM2500-30 compressor. The change of thermodynamic parameters such as pressure ratio, efficiency with the operating time, and stage number is analyzed in detail.

Problems in radiation shielding calculations with Monte Carlo methods

International Nuclear Information System (INIS)

Ueki, Kohtaro

1985-01-01

The Monte Carlo method is a very useful tool for solving a large class of radiation transport problem. In contrast with deterministic method, geometric complexity is a much less significant problem for Monte Carlo calculations. However, the accuracy of Monte Carlo calculations is of course, limited by statistical error of the quantities to be estimated. In this report, we point out some typical problems to solve a large shielding system including radiation streaming. The Monte Carlo coupling technique was developed to settle such a shielding problem accurately. However, the variance of the Monte Carlo results using the coupling technique of which detectors were located outside the radiation streaming, was still not enough. So as to bring on more accurate results for the detectors located outside the streaming and also for a multi-legged-duct streaming problem, a practicable way of ''Prism Scattering technique'' is proposed in the study. (author)

Friction and wear calculation methods

CERN Document Server

Kragelsky, I V; Kombalov, V S

1981-01-01

Friction and Wear: Calculation Methods provides an introduction to the main theories of a new branch of mechanics known as """"contact interaction of solids in relative motion."""" This branch is closely bound up with other sciences, especially physics and chemistry. The book analyzes the nature of friction and wear, and some theoretical relationships that link the characteristics of the processes and the properties of the contacting bodies essential for practical application of the theories in calculating friction forces and wear values. The effect of the environment on friction and wear is a

Eigenvalue translation method for mode calculations

International Nuclear Information System (INIS)

Gerck, E.; Cruz, C.H.B.

1978-11-01

A new method is described for the first few modes calculations in a interferometer that has several advantages over the ALLMAT subroutine, the Prony Method and the Fox and Li Method. In the illustrative results shown for the same cases it can be seen that the eigenvalue translation method is typically 100 fold times faster than the usual Fox and Li Method and 10 times faster than ALLMAT [pt

Multi-reference approach to the calculation of photoelectron spectra including spin-orbit coupling

Energy Technology Data Exchange (ETDEWEB)

Grell, Gilbert; Bokarev, Sergey I., E-mail: sergey.bokarev@uni-rostock.de; Kühn, Oliver [Institut für Physik, Universität Rostock, D-18051 Rostock (Germany); Winter, Bernd; Seidel, Robert [Helmholtz-Zentrum Berlin für Materialien und Energie, Methods for Material Development, Albert-Einstein-Strasse 15, D-12489 Berlin (Germany); Aziz, Emad F. [Helmholtz-Zentrum Berlin für Materialien und Energie, Methods for Material Development, Albert-Einstein-Strasse 15, D-12489 Berlin (Germany); Department of Physics, Freie Universität Berlin, Arnimalle 14, D-14159 Berlin (Germany); Aziz, Saadullah G. [Chemistry Department, Faculty of Science, King Abdulaziz University, 21589 Jeddah (Saudi Arabia)

2015-08-21

X-ray photoelectron spectra provide a wealth of information on the electronic structure. The extraction of molecular details requires adequate theoretical methods, which in case of transition metal complexes has to account for effects due to the multi-configurational and spin-mixed nature of the many-electron wave function. Here, the restricted active space self-consistent field method including spin-orbit coupling is used to cope with this challenge and to calculate valence- and core-level photoelectron spectra. The intensities are estimated within the frameworks of the Dyson orbital formalism and the sudden approximation. Thereby, we utilize an efficient computational algorithm that is based on a biorthonormal basis transformation. The approach is applied to the valence photoionization of the gas phase water molecule and to the core ionization spectrum of the [Fe(H{sub 2}O){sub 6}]{sup 2+} complex. The results show good agreement with the experimental data obtained in this work, whereas the sudden approximation demonstrates distinct deviations from experiments.

Burnup calculations using Monte Carlo method

International Nuclear Information System (INIS)

Ghosh, Biplab; Degweker, S.B.

2009-01-01

In the recent years, interest in burnup calculations using Monte Carlo methods has gained momentum. Previous burn up codes have used multigroup transport theory based calculations followed by diffusion theory based core calculations for the neutronic portion of codes. The transport theory methods invariably make approximations with regard to treatment of the energy and angle variables involved in scattering, besides approximations related to geometry simplification. Cell homogenisation to produce diffusion, theory parameters adds to these approximations. Moreover, while diffusion theory works for most reactors, it does not produce accurate results in systems that have strong gradients, strong absorbers or large voids. Also, diffusion theory codes are geometry limited (rectangular, hexagonal, cylindrical, and spherical coordinates). Monte Carlo methods are ideal to solve very heterogeneous reactors and/or lattices/assemblies in which considerable burnable poisons are used. The key feature of this approach is that Monte Carlo methods permit essentially 'exact' modeling of all geometrical detail, without resort to ene and spatial homogenization of neutron cross sections. Monte Carlo method would also be better for in Accelerator Driven Systems (ADS) which could have strong gradients due to the external source and a sub-critical assembly. To meet the demand for an accurate burnup code, we have developed a Monte Carlo burnup calculation code system in which Monte Carlo neutron transport code is coupled with a versatile code (McBurn) for calculating the buildup and decay of nuclides in nuclear materials. McBurn is developed from scratch by the authors. In this article we will discuss our effort in developing the continuous energy Monte Carlo burn-up code, McBurn. McBurn is intended for entire reactor core as well as for unit cells and assemblies. Generally, McBurn can do burnup of any geometrical system which can be handled by the underlying Monte Carlo transport code

Monte Carlo method for array criticality calculations

International Nuclear Information System (INIS)

Dickinson, D.; Whitesides, G.E.

1976-01-01

The Monte Carlo method for solving neutron transport problems consists of mathematically tracing paths of individual neutrons collision by collision until they are lost by absorption or leakage. The fate of the neutron after each collision is determined by the probability distribution functions that are formed from the neutron cross-section data. These distributions are sampled statistically to establish the successive steps in the neutron's path. The resulting data, accumulated from following a large number of batches, are analyzed to give estimates of k/sub eff/ and other collision-related quantities. The use of electronic computers to produce the simulated neutron histories, initiated at Los Alamos Scientific Laboratory, made the use of the Monte Carlo method practical for many applications. In analog Monte Carlo simulation, the calculation follows the physical events of neutron scattering, absorption, and leakage. To increase calculational efficiency, modifications such as the use of statistical weights are introduced. The Monte Carlo method permits the use of a three-dimensional geometry description and a detailed cross-section representation. Some of the problems in using the method are the selection of the spatial distribution for the initial batch, the preparation of the geometry description for complex units, and the calculation of error estimates for region-dependent quantities such as fluxes. The Monte Carlo method is especially appropriate for criticality safety calculations since it permits an accurate representation of interacting units of fissile material. Dissimilar units, units of complex shape, moderators between units, and reflected arrays may be calculated. Monte Carlo results must be correlated with relevant experimental data, and caution must be used to ensure that a representative set of neutron histories is produced

Methods for calculating nonconcave entropies

International Nuclear Information System (INIS)

Touchette, Hugo

2010-01-01

Five different methods which can be used to analytically calculate entropies that are nonconcave as functions of the energy in the thermodynamic limit are discussed and compared. The five methods are based on the following ideas and techniques: (i) microcanonical contraction, (ii) metastable branches of the free energy, (iii) generalized canonical ensembles with specific illustrations involving the so-called Gaussian and Betrag ensembles, (iv) the restricted canonical ensemble, and (v) the inverse Laplace transform. A simple long-range spin model having a nonconcave entropy is used to illustrate each method

Calculation of an axisymmetric current coil field with the bounding contour integration method

Energy Technology Data Exchange (ETDEWEB)

Telegin, Alexander P.; Klevets, Nickolay I. E-mail: pmsolution@mail.ru

2004-06-01

Method for the economic and stable (in the sense of calculation errors) analysis of an induction of a magnetic field created with axisymmetric coils in arbitrary points of space, including points located inside a coil or on its border, is obtained. The basic idea of the method is to replace a current coil with continuous distribution of current density by magnetization distributed in the volume of the coil and creating the equivalent magnetic field. This allows to use field surface sources at calculation of the fields. Consequently, the range of integration is reduced resulting in reduction of calculation volume by an order in most cases. Besides, the calculation of improper integrals in internal points and on the border is completely excluded.

Calculation of an axisymmetric current coil field with the bounding contour integration method

International Nuclear Information System (INIS)

Telegin, Alexander P.; Klevets, Nickolay I.

2004-01-01

Method for the economic and stable (in the sense of calculation errors) analysis of an induction of a magnetic field created with axisymmetric coils in arbitrary points of space, including points located inside a coil or on its border, is obtained. The basic idea of the method is to replace a current coil with continuous distribution of current density by magnetization distributed in the volume of the coil and creating the equivalent magnetic field. This allows to use field surface sources at calculation of the fields. Consequently, the range of integration is reduced resulting in reduction of calculation volume by an order in most cases. Besides, the calculation of improper integrals in internal points and on the border is completely excluded

Classification of methods for annual energy harvesting calculations of photovoltaic generators

International Nuclear Information System (INIS)

Rus-Casas, C.; Aguilar, J.D.; Rodrigo, P.; Almonacid, F.; Pérez-Higueras, P.J.

2014-01-01

Highlights: • The paper presents a novel classification of methods for annual energy harvesting calculation of grid-connected PV systems. • The methods are classified in direct and indirect methods. • Direct methods directly calculate the energy. Indirect methods calculate the energy from the power. • The classification can help the PV professionals in order to choose the most suitable method for each application. - Abstract: Estimating the energy provided by the generators of grid-connected photovoltaic systems is important in order to analyze their economic viability and supervise their operation. The energy harvesting calculation of a photovoltaic generator is not trivial; there are a lot of methods for this calculation. The aim of this paper is to develop a novel classification of methods for annual energy harvesting calculation of a generator of a grid-connected photovoltaic system. The methods are classified in two groups: (1) those that indirectly calculate the energy, i.e. they first calculate the power and from this, they calculate the energy, and (2) those that directly calculate the energy. Furthermore, the indirect methods are grouped in two categories: those that first calculate the I–V curve of the generator and from this, they calculate the power, and those that directly calculate the power. The study has shown that the existing methods differ in simplicity and accuracy, so that the proposed classification is useful in order to choose the most suitable method for each specific application

SOLAR OPACITY CALCULATIONS USING THE SUPER-TRANSITION-ARRAY METHOD

International Nuclear Information System (INIS)

Krief, M.; Feigel, A.; Gazit, D.

2016-01-01

A new opacity model has been developed based on the Super-Transition-Array (STA) method for the calculation of monochromatic opacities of plasmas in local thermodynamic equilibrium. The atomic code, named STAR (STA-Revised), is described and used to calculate spectral opacities for a solar model implementing the recent AGSS09 composition. Calculations are carried out throughout the solar radiative zone. The relative contributions of different chemical elements and atomic processes to the total Rosseland mean opacity are analyzed in detail. Monochromatic opacities and charge-state distributions are compared with the widely used Opacity Project (OP) code, for several elements near the radiation–convection interface. STAR Rosseland opacities for the solar mixture show a very good agreement with OP and the OPAL opacity code throughout the radiation zone. Finally, an explicit STA calculation was performed of the full AGSS09 photospheric mixture, including all heavy metals. It was shown that, due to their extremely low abundance, and despite being very good photon absorbers, the heavy elements do not affect the Rosseland opacity

SOLAR OPACITY CALCULATIONS USING THE SUPER-TRANSITION-ARRAY METHOD

Energy Technology Data Exchange (ETDEWEB)

Krief, M.; Feigel, A.; Gazit, D., E-mail: menahem.krief@mail.huji.ac.il [The Racah Institute of Physics, The Hebrew University, 91904 Jerusalem (Israel)

2016-04-10

A new opacity model has been developed based on the Super-Transition-Array (STA) method for the calculation of monochromatic opacities of plasmas in local thermodynamic equilibrium. The atomic code, named STAR (STA-Revised), is described and used to calculate spectral opacities for a solar model implementing the recent AGSS09 composition. Calculations are carried out throughout the solar radiative zone. The relative contributions of different chemical elements and atomic processes to the total Rosseland mean opacity are analyzed in detail. Monochromatic opacities and charge-state distributions are compared with the widely used Opacity Project (OP) code, for several elements near the radiation–convection interface. STAR Rosseland opacities for the solar mixture show a very good agreement with OP and the OPAL opacity code throughout the radiation zone. Finally, an explicit STA calculation was performed of the full AGSS09 photospheric mixture, including all heavy metals. It was shown that, due to their extremely low abundance, and despite being very good photon absorbers, the heavy elements do not affect the Rosseland opacity.

Calculation of one-loop anomalous dimensions by means of the background field method

International Nuclear Information System (INIS)

Morozov, A.Yu.

1983-01-01

The knowledge of propagators in background fields makes calculation of anomalous dimensions (AD) straightforward and brief. The paper illustrates this statement by calculation of AD of many spin-zero and one QCD operators up to the eighth dimension included. The method presented does not simplify calculations in case of four-quark operators, therefore these are not discussed. Together with calculational difficulties arising for operators with derivatives this limits capacities of the whole approach and leads to incompleteness of some mixing matrices found in the article

Calculation of radon concentration in water by toluene extraction method

Energy Technology Data Exchange (ETDEWEB)

Saito, Masaaki [Tokyo Metropolitan Isotope Research Center (Japan)

1997-02-01

Noguchi method and Horiuchi method have been used as the calculation method of radon concentration in water. Both methods have two problems in the original, that is, the concentration calculated is changed by the extraction temperature depend on the incorrect solubility data and the concentration calculated are smaller than the correct values, because the radon calculation equation does not true to the gas-liquid equilibrium theory. However, the two problems are solved by improving the radon equation. I presented the Noguchi-Saito equation and the constant B of Horiuchi-Saito equation. The calculating results by the improved method showed about 10% of error. (S.Y.)

Calculation method of passive tracer using regional ocean modeling system (ROMS)

International Nuclear Information System (INIS)

Tubono, Takaki; Tsumune, Daisuke; Misumi, Kazuhiro; Yoshida, Yoshikatsu

2010-01-01

A dose assessments method for packages shipping radioactive materials was developed to calculate concentrations of radioactive materials discharged from the hypothetically sinking point using Regional Ocean Modeling System (ROMS). The programs of ROMS were improved for simultaneously calculating the flow fields around Japan and the concentrations discharged from several bottom points in the continental shelf. The model domain covered the northeast pacific region around Japan, 20degN-65degN, 120degE-165degE, with 1/10 degree horizontal resolution and 30 sigma-levels. The ROMS model captured main features such as Kuroshio currents, including the separation location and the Kuroshio extension, and Oyashio currents. The calculated volume transports through the main straits such as Tokara strait, Tsushima strait and Tsugaru strait were compared well to the measured or estimated that. The calculated concentrations discharged from 12 bottom points showed mixing associated with the eddies of ocean currents. This method is available to evaluate effectively the concentrations of radioactive materials only by setting the subsidence location and the concentrations and volume discharged. (author)

Efficient methods for time-absorption (α) eigenvalue calculations

International Nuclear Information System (INIS)

Hill, T.R.

1983-01-01

The time-absorption eigenvalue (α) calculation is one of the options found in most discrete-ordinates transport codes. Several methods have been developed at Los Alamos to improve the efficiency of this calculation. Two procedures, based on coarse-mesh rebalance, to accelerate the α eigenvalue search are derived. A hybrid scheme to automatically choose the more-effective rebalance method is described. The α rebalance scheme permits some simple modifications to the iteration strategy that eliminates many unnecessary calculations required in the standard search procedure. For several fast supercritical test problems, these methods resulted in convergence with one-fifth the number of iterations required for the conventional eigenvalue search procedure

Simple method for calculating island widths

International Nuclear Information System (INIS)

Cary, J.R.; Hanson, J.D.; Carreras, B.A.; Lynch, V.E.

1989-01-01

A simple method for calculating magnetic island widths has been developed. This method uses only information obtained from integrating along the closed field line at the island center. Thus, this method is computationally less intensive than the usual method of producing surfaces of section of sufficient detail to locate and resolve the island separatrix. This method has been implemented numerically and used to analyze the buss work islands of ATF. In this case the method proves to be accurate to at least within 30%. 7 refs

A method for bubble volume calculating in vertical two-phase flow

International Nuclear Information System (INIS)

Wang, H Y; Dong, F

2009-01-01

The movement of bubble is a basic subject in gas-liquid two-phase flow research. A method for calculating bubble volume which is one of the most important characters in bubble motion research was proposed. A suit of visualized experimental device was designed and set up. Single bubble rising in stagnant liquid in a rectangular tank was studied using the high-speed video system. Bubbles generated by four orifice with different diameter (1mm, 2mm, 3mm, 4mm) were recorded respectively. Sequences of recorded high-speed images were processed by digital image processing method, such as image noise remove, binary image transform, bubble filling, and so on. then, Several parameters could be obtained from the processed image. Bubble area, equivalent diameter, bubble velocity, bubble acceleration are all indispensable in bubble volume calculating. In order to get the force balance equation, forces that work on bubble along vertical direction, including drag force, virtual mass force, buoyancy, gravity and liquid thrust, were analyzed. Finally, the bubble volume formula could be derived from the force balance equation and bubble parameters. Examples were given to shown the computing process and results. Comparison of the bubble volume calculated by geomettic method and the present method have shown the superiority of the proposed method in this paper.

Iterative resonance self-shielding methods using resonance integral table in heterogeneous transport lattice calculations

International Nuclear Information System (INIS)

Hong, Ser Gi; Kim, Kang-Seog

2011-01-01

This paper describes the iteration methods using resonance integral tables to estimate the effective resonance cross sections in heterogeneous transport lattice calculations. Basically, these methods have been devised to reduce an effort to convert resonance integral table into subgroup data to be used in the physical subgroup method. Since these methods do not use subgroup data but only use resonance integral tables directly, these methods do not include an error in converting resonance integral into subgroup data. The effective resonance cross sections are estimated iteratively for each resonance nuclide through the heterogeneous fixed source calculations for the whole problem domain to obtain the background cross sections. These methods have been implemented in the transport lattice code KARMA which uses the method of characteristics (MOC) to solve the transport equation. The computational results show that these iteration methods are quite promising in the practical transport lattice calculations.

Methods for calculation of compensation for expropriation of a foreign investment

Directory of Open Access Journals (Sweden)

Đundić Petar M.

2015-01-01

Full Text Available A legal expropriation of a foreign investment without exception includes an obligation of the host state to pay the investor an appropriate compensation for the taking. Bilateral treaties for protection and encouragement of investments which usually serve as a normative basis for arbitration of expropriation disputes contain provisions instructing contracting states to provide 'adequate', 'just' compensation or to pay the compensation in 'full' or 'genuine' value of the expropriated investment etc. However, the sole existence of the standard for payment of compensation is not by itself enough for establishing the precise sum of compensation in particular cases. This purpose is served by mathematical methods of calculation employed used in judicial and arbitral practice. The paper contains an overview of the most important methods for calculation in practice: the market value method, the discounted cash flow method and the book value method. It also identifies the differences between situations in which different methods are used and explains how the nature of an investment and circumstances of particular case affect the choice of a particular method by the arbitral tribunal.

Comparison study on cell calculation method of fast reactor

International Nuclear Information System (INIS)

Chiba, Gou

2002-10-01

Effective cross sections obtained by cell calculations are used in core calculations in current deterministic methods. Therefore, it is important to calculate the effective cross sections accurately and several methods have been proposed. In this study, some of the methods are compared to each other using a continuous energy Monte Carlo method as a reference. The result shows that the table look-up method used in Japan Nuclear Cycle Development Institute (JNC) sometimes has a difference over 10% in effective microscopic cross sections and be inferior to the sub-group method. The problem was overcome by introducing a new nuclear constant system developed in JNC, in which the ultra free energy group library is used. The system can also deal with resonance interaction effects between nuclides which are not able to be considered by other methods. In addition, a new method was proposed to calculate effective cross section accurately for power reactor fuel subassembly where the new nuclear constant system cannot be applied. This method uses the sub-group method and the ultra fine energy group collision probability method. The microscopic effective cross sections obtained by this method agree with the reference values within 5% difference. (author)

Classical Methods and Calculation Algorithms for Determining Lime Requirements

Directory of Open Access Journals (Sweden)

André Guarçoni

Full Text Available ABSTRACT The methods developed for determination of lime requirements (LR are based on widely accepted principles. However, the formulas used for calculation have evolved little over recent decades, and in some cases there are indications of their inadequacy. The aim of this study was to compare the lime requirements calculated by three classic formulas and three algorithms, defining those most appropriate for supplying Ca and Mg to coffee plants and the smaller possibility of causing overliming. The database used contained 600 soil samples, which were collected in coffee plantings. The LR was estimated by the methods of base saturation, neutralization of Al3+, and elevation of Ca2+ and Mg2+ contents (two formulas and by the three calculation algorithms. Averages of the lime requirements were compared, determining the frequency distribution of the 600 lime requirements (LR estimated through each calculation method. In soils with low cation exchange capacity at pH 7, the base saturation method may fail to adequately supply the plants with Ca and Mg in many situations, while the method of Al3+ neutralization and elevation of Ca2+ and Mg2+ contents can result in the calculation of application rates that will increase the pH above the suitable range. Among the methods studied for calculating lime requirements, the algorithm that predicts reaching a defined base saturation, with adequate Ca and Mg supply and the maximum application rate limited to the H+Al value, proved to be the most efficient calculation method, and it can be recommended for use in numerous crops conditions.

MERSENNE AND HADAMARD MATRICES CALCULATION BY SCARPIS METHOD

Directory of Open Access Journals (Sweden)

N. A. Balonin

2014-05-01

Full Text Available Purpose. The paper deals with the problem of basic generalizations of Hadamard matrices associated with maximum determinant matrices or not optimal by determinant matrices with orthogonal columns (weighing matrices, Mersenne and Euler matrices, ets.; calculation methods for the quasi-orthogonal local maximum determinant Mersenne matrices are not studied enough sufficiently. The goal of this paper is to develop the theory of Mersenne and Hadamard matrices on the base of generalized Scarpis method research. Methods. Extreme solutions are found in general by minimization of maximum for absolute values of the elements of studied matrices followed by their subsequent classification according to the quantity of levels and their values depending on orders. Less universal but more effective methods are based on structural invariants of quasi-orthogonal matrices (Silvester, Paley, Scarpis methods, ets.. Results. Generalizations of Hadamard and Belevitch matrices as a family of quasi-orthogonal matrices of odd orders are observed; they include, in particular, two-level Mersenne matrices. Definitions of section and layer on the set of generalized matrices are proposed. Calculation algorithms for matrices of adjacent layers and sections by matrices of lower orders are described. Approximation examples of the Belevitch matrix structures up to 22-nd critical order by Mersenne matrix of the third order are given. New formulation of the modified Scarpis method to approximate Hadamard matrices of high orders by lower order Mersenne matrices is proposed. Williamson method is described by example of one modular level matrices approximation by matrices with a small number of levels. Practical relevance. The efficiency of developing direction for the band-pass filters creation is justified. Algorithms for Mersenne matrices design by Scarpis method are used in developing software of the research program complex. Mersenne filters are based on the suboptimal by

Calculation Methods for Wallenius’ Noncentral Hypergeometric Distribution

DEFF Research Database (Denmark)

Fog, Agner

2008-01-01

Two different probability distributions are both known in the literature as "the" noncentral hypergeometric distribution. Wallenius' noncentral hypergeometric distribution can be described by an urn model without replacement with bias. Fisher's noncentral hypergeometric distribution...... is the conditional distribution of independent binomial variates given their sum. No reliable calculation method for Wallenius' noncentral hypergeometric distribution has hitherto been described in the literature. Several new methods for calculating probabilities from Wallenius' noncentral hypergeometric...... distribution are derived. Range of applicability, numerical problems, and efficiency are discussed for each method. Approximations to the mean and variance are also discussed. This distribution has important applications in models of biased sampling and in models of evolutionary systems....

Failure Probability Calculation Method Using Kriging Metamodel-based Importance Sampling Method

Energy Technology Data Exchange (ETDEWEB)

Lee, Seunggyu [Korea Aerospace Research Institue, Daejeon (Korea, Republic of); Kim, Jae Hoon [Chungnam Nat’l Univ., Daejeon (Korea, Republic of)

2017-05-15

The kernel density was determined based on sampling points obtained in a Markov chain simulation and was assumed to be an important sampling function. A Kriging metamodel was constructed in more detail in the vicinity of a limit state. The failure probability was calculated based on importance sampling, which was performed for the Kriging metamodel. A pre-existing method was modified to obtain more sampling points for a kernel density in the vicinity of a limit state. A stable numerical method was proposed to find a parameter of the kernel density. To assess the completeness of the Kriging metamodel, the possibility of changes in the calculated failure probability due to the uncertainty of the Kriging metamodel was calculated.

Analyzed method for calculating the distribution of electrostatic field

International Nuclear Information System (INIS)

Lai, W.

1981-01-01

An analyzed method for calculating the distribution of electrostatic field under any given axial gradient in tandem accelerators is described. This method possesses satisfactory accuracy compared with the results of numerical calculation

The pseudo-harmonics method applied to depletion calculation

International Nuclear Information System (INIS)

Silva, F.C. da; Amaral, J.A.C.; Thome, Z.D.

1989-01-01

In this paper, a new method for performing depletion calculations, based on the use of the Pseudo-Harmonics perturbation method, was developed. The fuel burnup was considered as a global perturbation and the multigroup difusion equations were rewriten in such a way as to treat the soluble boron concentration as the eigenvalue. By doing this, the critical boron concentration can be obtained by a perturbation method. A test of the new method was performed for a H 2 O-colled, D 2 O-moderated reactor. Comparison with direct calculation showed that this method is very accurate and efficient. (author) [pt

Electromagnetically induced nuclear beta decay calculated by a Green's function method

International Nuclear Information System (INIS)

Reiss, H.R.

1984-01-01

The transition probability for enhancement of forbidden nuclear beta decay by an applied plane-wave electromagnetic field is calculated in a nonrelativistic spinless approximation by a Green's function method. The calculation involves a stationary-phase approximation. The stationary phase points in the presence of an intense field are located in very different positions than they are in the field-free case. In order-of-magnitude terms, the results are completely consistent with an earlier, much more complete wave-function calculation which includes spin and relativistic effects. Both the present Green's function calculation and the earlier wave function calculation give electromagnetic contributions in first-forbidden nuclear beta decay matrix elements which are of order (R 0 /lambda-dash-bar/sub C/) 2 with respect to allowed decays, where R 0 is the nuclear radius and lambda-dash-bar/sub C/ is the electron Compton wavelength

Simple Calculation Programs for Biology Immunological Methods

Indian Academy of Sciences (India)

First page Back Continue Last page Overview Graphics. Simple Calculation Programs for Biology Immunological Methods. Computation of Ab/Ag Concentration from EISA data. Graphical Method; Raghava et al., 1992, J. Immuno. Methods 153: 263. Determination of affinity of Monoclonal Antibody. Using non-competitive ...

Comparison of different dose calculation methods for irregular photon fields

International Nuclear Information System (INIS)

Zakaria, G.A.; Schuette, W.

2000-01-01

In this work, 4 calculation methods (Wrede method, Clarskon method of sector integration, beam-zone method of Quast and pencil-beam method of Ahnesjoe) are introduced to calculate point doses in different irregular photon fields. The calculations cover a typical mantle field, an inverted Y-field and different blocked fields for 4 and 10 MV photon energies. The results are compared to those of measurements in a water phantom. The Clarkson and the pencil-beam method have been proved to be the methods of equal standard in relation to accuracy. Both of these methods are being distinguished by minimum deviations and applied in our clinical routine work. The Wrede and beam-zone methods deliver useful results to central beam and yet provide larger deviations in calculating points beyond the central axis. (orig.) [de

Consideration of the environmental effects on fatigue behavior of austenitic components. Calculation methods and practical application

International Nuclear Information System (INIS)

Seichter, Johannes; Reese, Sven H.; Klucke, Dietmar

2012-01-01

During the last years environmental effects on the fatigue behavior of nuclear power plant components has worldwide been discussed controversial with respect to the transferability of laboratory data on real components. A publication from Argonne National Laboratory on experimental results concerning environmental effects (air and LWR coolant) on fatigue of austenitic steels included a proposal on calculation methods concerning the lifetime reduction due to environmental effects. This calculation method, i.e. multiplication of the usage factor by a F(en), has been included into the ASME Code, Section III, Division I, as Code Case N-792 (fatigue evaluations including environmental effects). The presented contribution evaluates the practical application of this calculation procedure and demonstrates the determination of the usage factor of an austenitic component under environmental exposure.

A unique manual method for emergency offsite dose calculations

International Nuclear Information System (INIS)

Wildner, T.E.; Carson, B.H.; Shank, K.E.

1987-01-01

This paper describes a manual method developed for performance of emergency offsite dose calculations for PP and L's Susquehanna Steam Electric Station. The method is based on a three-part carbonless form. The front page guides the user through selection of the appropriate accident case and inclusion of meteorological and effluent data data. By circling the applicable accident descriptors, the user circles the dose factors on pages 2 and 3 which are then simply multiplied to yield the whole body and thyroid dose rates at the plant boundary, two, five, and ten miles. The process used to generate the worksheet is discussed, including the method used to incorporate the observed terrain effects on airflow patterns caused by the Susquehanna River Valley topography

Three-dimensional space-charge calculation method

International Nuclear Information System (INIS)

Lysenko, W.P.; Wadlinger, E.A.

1980-09-01

A method is presented for calculating space-charge forces on individual particles in a particle tracing simulation code. Poisson's equation is solved in three dimensions with boundary conditions specified on an arbitrary surface. When the boundary condition is defined by an impressed radio-frequency field, the external electric fields as well as the space-charge fields are determined. A least squares fitting procedure is used to calculate the coefficients of expansion functions, which need not be orthogonal nor individually satisfy the boundary condition

Calculation method for gamma dose rates from Gaussian puffs

Energy Technology Data Exchange (ETDEWEB)

Thykier-Nielsen, S; Deme, S; Lang, E

1995-06-01

The Lagrangian puff models are widely used for calculation of the dispersion of releases to the atmosphere. Basic output from such models is concentration of material in the air and on the ground. The most simple method for calculation of the gamma dose from the concentration of airborne activity is based on the semi-infinite cloud model. This method is however only applicable for puffs with large dispersion parameters, i.e. for receptors far away from the release point. The exact calculation of the cloud dose using volume integral requires large computer time usually exceeding what is available for real time calculations. The volume integral for gamma doses could be approximated by using the semi-infinite cloud model combined with correction factors. This type of calculation procedure is very fast, but usually the accuracy is poor because only a few of the relevant parameters are considered. A multi-parameter method for calculation of gamma doses is described here. This method uses precalculated values of the gamma dose rates as a function of E{sub {gamma}}, {sigma}{sub y}, the asymmetry factor - {sigma}{sub y}/{sigma}{sub z}, the height of puff center - H and the distance from puff center R{sub xy}. To accelerate the calculations the release energy, for each significant radionuclide in each energy group, has been calculated and tabulated. Based on the precalculated values and suitable interpolation procedure the calculation of gamma doses needs only short computing time and it is almost independent of the number of radionuclides considered. (au) 2 tabs., 15 ills., 12 refs.

Calculation method for gamma dose rates from Gaussian puffs

International Nuclear Information System (INIS)

Thykier-Nielsen, S.; Deme, S.; Lang, E.

1995-06-01

The Lagrangian puff models are widely used for calculation of the dispersion of releases to the atmosphere. Basic output from such models is concentration of material in the air and on the ground. The most simple method for calculation of the gamma dose from the concentration of airborne activity is based on the semi-infinite cloud model. This method is however only applicable for puffs with large dispersion parameters, i.e. for receptors far away from the release point. The exact calculation of the cloud dose using volume integral requires large computer time usually exceeding what is available for real time calculations. The volume integral for gamma doses could be approximated by using the semi-infinite cloud model combined with correction factors. This type of calculation procedure is very fast, but usually the accuracy is poor because only a few of the relevant parameters are considered. A multi-parameter method for calculation of gamma doses is described here. This method uses precalculated values of the gamma dose rates as a function of E γ , σ y , the asymmetry factor - σ y /σ z , the height of puff center - H and the distance from puff center R xy . To accelerate the calculations the release energy, for each significant radionuclide in each energy group, has been calculated and tabulated. Based on the precalculated values and suitable interpolation procedure the calculation of gamma doses needs only short computing time and it is almost independent of the number of radionuclides considered. (au) 2 tabs., 15 ills., 12 refs

Quantum Monte Carlo diagonalization method as a variational calculation

International Nuclear Information System (INIS)

Mizusaki, Takahiro; Otsuka, Takaharu; Honma, Michio.

1997-01-01

A stochastic method for performing large-scale shell model calculations is presented, which utilizes the auxiliary field Monte Carlo technique and diagonalization method. This method overcomes the limitation of the conventional shell model diagonalization and can extremely widen the feasibility of shell model calculations with realistic interactions for spectroscopic study of nuclear structure. (author)

Iterative acceleration methods for Monte Carlo and deterministic criticality calculations

Energy Technology Data Exchange (ETDEWEB)

Urbatsch, T.J.

1995-11-01

If you have ever given up on a nuclear criticality calculation and terminated it because it took so long to converge, you might find this thesis of interest. The author develops three methods for improving the fission source convergence in nuclear criticality calculations for physical systems with high dominance ratios for which convergence is slow. The Fission Matrix Acceleration Method and the Fission Diffusion Synthetic Acceleration (FDSA) Method are acceleration methods that speed fission source convergence for both Monte Carlo and deterministic methods. The third method is a hybrid Monte Carlo method that also converges for difficult problems where the unaccelerated Monte Carlo method fails. The author tested the feasibility of all three methods in a test bed consisting of idealized problems. He has successfully accelerated fission source convergence in both deterministic and Monte Carlo criticality calculations. By filtering statistical noise, he has incorporated deterministic attributes into the Monte Carlo calculations in order to speed their source convergence. He has used both the fission matrix and a diffusion approximation to perform unbiased accelerations. The Fission Matrix Acceleration method has been implemented in the production code MCNP and successfully applied to a real problem. When the unaccelerated calculations are unable to converge to the correct solution, they cannot be accelerated in an unbiased fashion. A Hybrid Monte Carlo method weds Monte Carlo and a modified diffusion calculation to overcome these deficiencies. The Hybrid method additionally possesses reduced statistical errors.

Iterative acceleration methods for Monte Carlo and deterministic criticality calculations

International Nuclear Information System (INIS)

Urbatsch, T.J.

1995-11-01

If you have ever given up on a nuclear criticality calculation and terminated it because it took so long to converge, you might find this thesis of interest. The author develops three methods for improving the fission source convergence in nuclear criticality calculations for physical systems with high dominance ratios for which convergence is slow. The Fission Matrix Acceleration Method and the Fission Diffusion Synthetic Acceleration (FDSA) Method are acceleration methods that speed fission source convergence for both Monte Carlo and deterministic methods. The third method is a hybrid Monte Carlo method that also converges for difficult problems where the unaccelerated Monte Carlo method fails. The author tested the feasibility of all three methods in a test bed consisting of idealized problems. He has successfully accelerated fission source convergence in both deterministic and Monte Carlo criticality calculations. By filtering statistical noise, he has incorporated deterministic attributes into the Monte Carlo calculations in order to speed their source convergence. He has used both the fission matrix and a diffusion approximation to perform unbiased accelerations. The Fission Matrix Acceleration method has been implemented in the production code MCNP and successfully applied to a real problem. When the unaccelerated calculations are unable to converge to the correct solution, they cannot be accelerated in an unbiased fashion. A Hybrid Monte Carlo method weds Monte Carlo and a modified diffusion calculation to overcome these deficiencies. The Hybrid method additionally possesses reduced statistical errors

Calculating the albedo characteristics by the method of transmission probabilities

International Nuclear Information System (INIS)

Lukhvich, A.A.; Rakhno, I.L.; Rubin, I.E.

1983-01-01

The possibility to use the method of transmission probabilities for calculating the albedo characteristics of homogeneous and heterogeneous zones is studied. The transmission probabilities method is a numerical method for solving the transport equation in the integrated form. All calculations have been conducted as a one-group approximation for the planes and rods with different optical thicknesses and capture-to-scattering ratios. Above calculations for plane and cylindrical geometries have shown the possibility to use the numerical method of transmission probabilities for calculating the albedo characteristics of homogeneous and heterogeneous zones with high accuracy. In this case the computer time consumptions are minimum even with the cylindrical geometry, if the interpolation calculation of characteristics is used for the neutrons of the first path

Evolution of calculation methods taking into account severe accidents

International Nuclear Information System (INIS)

L'Homme, A.; Courtaud, J.M.

1990-12-01

During the first decade of PWRs operation in France the calculation methods used for design and operation have improved very much. This paper gives a general analysis of the calculation methods evolution in parallel with the evolution of safety approach concerning PWRs. Then a comprehensive presentation of principal calculation tools is presented as applied during the past decade. An effort is done to predict the improvements in near future

Comparison of calculation methods. PME and LTA; Vergelijking methoden energiebesparing. PME en MJA

Energy Technology Data Exchange (ETDEWEB)

Verdonk, M. [Planbureau voor de Leefomgeving PBL, Bilthoven (Netherlands); Boonekamp, P. [ECN Beleidsstudies, Petten (Netherlands)

2009-11-15

Various calculation methods are used to calculate the effects of energy conservation policy. This memo provides insight in the main differences between two frequently used methods, i.e. the second Long Term Agreements Energy Efficiency (LTA2) and the Protocol Monitoring Energy Saving (PME). These insights enable better fine-tuning between PME and the successor to LTA2 (LTA3). The memo starts with a brief characterization of the LTAs, the corresponding calculation methods and PME approach. Next the most striking differences are discussed, including the quantitative effect of the two calculation methods. [Dutch] Er zijn verschillende rekenmethoden in gebruik, waarmee de effecten van het energiebesparingsbeleid kunnen worden berekend. Deze notitie verschaft inzicht in de voornaamste verschillen tussen twee veel gebruikte methoden, namelijk de tweede Meerjarenafspraken energie-efficiency (MJA2) en het Protocol Monitoring Energiebesparing (PME). De inzichten maken een betere afstemming mogelijk tussen PME en de opvolger van MJA2 (MJA3). De notitie start met een korte typering van de MJA's, de bijbehorende rekenmethode en de PME aanpak. Vervolgens worden de meest opvallende verschillen belicht, inclusief het kwantitatieve effect van de twee rekenmethoden.

Comments on Simplified Calculation Method for Fire Exposed Concrete Columns

DEFF Research Database (Denmark)

Hertz, Kristian Dahl

1998-01-01

The author has developed new simplified calculation methods for fire exposed columns. Methods, which are found In ENV 1992-1-2 chapter 4.3 and in proposal for Danish code of Practise DS411 chapter 9. In the present supporting document the methods are derived and 50 eccentrically loaded fire expos...... columns are calculated and compared to results of full-scale tests. Furthermore 500 columns are calculated in order to present each test result related to a variation of the calculation in time of fire resistance....

Calculational methods for lattice cells

International Nuclear Information System (INIS)

Askew, J.R.

1980-01-01

At the current stage of development, direct simulation of all the processes involved in the reactor to the degree of accuracy required is not an economic proposition, and this is achieved by progressive synthesis of models for parts of the full space/angle/energy neutron behaviour. The split between reactor and lattice calculations is one such simplification. Most reactors are constructed of repetitions of similar geometric units, the fuel elements, having broadly similar properties. Thus the provision of detailed predictions of their behaviour is an important step towards overall modelling. We shall be dealing with these lattice methods in this series of lectures, but will refer back from time to time to their relationship with overall reactor calculation The lattice cell is itself composed of somewhat similar sub-units, the fuel pins, and will itself often rely upon a further break down of modelling. Construction of a good model depends upon the identification, on physical and mathematical grounds, of the most helpful division of the calculation at this level

Conventional method for the calculation of the global energy cost of buildings; Methode conventionnelle de calcul du cout global energetique des batiments

Energy Technology Data Exchange (ETDEWEB)

NONE

2002-05-01

A working group driven by Electricite de France (EdF), Chauffage Fioul and Gaz de France (GdF) companies has been built with the sustain of several building engineering companies in order to clarify the use of the method of calculation of the global energy cost of buildings. This global cost is an economical decision help criterion among others. This press kit presents, first, the content of the method (input data, calculation of annual expenses, calculation of the global energy cost, display of results and limitations of the method). Then it fully describes the method and its appendixes necessary for its implementation: economical and financial context, general data of the project in progress, environmental data, occupation and comfort level, variants, investment cost of energy systems, investment cost for the structure linked with the energy system, investment cost for other invariant elements of the structure, calculation of consumptions (space heating, hot water, ventilation), maintenance costs (energy systems, structure), operation and exploitation costs, tariffs and consumption costs and taxes, actualized global cost, annualized global cost, comparison between variants. The method is applied to a council building of 23 flats taken as an example. (J.S.)

Method for calculating anisotropic neutron transport using scattering kernel without polynomial expansion

International Nuclear Information System (INIS)

Takahashi, Akito; Yamamoto, Junji; Ebisuya, Mituo; Sumita, Kenji

1979-01-01

A new method for calculating the anisotropic neutron transport is proposed for the angular spectral analysis of D-T fusion reactor neutronics. The method is based on the transport equation with new type of anisotropic scattering kernels formulated by a single function I sub(i) (μ', μ) instead of polynomial expansion, for instance, Legendre polynomials. In the calculation of angular flux spectra by using scattering kernels with the Legendre polynomial expansion, we often observe the oscillation with negative flux. But in principle this oscillation disappears by this new method. In this work, we discussed anisotropic scattering kernels of the elastic scattering and the inelastic scatterings which excite discrete energy levels. The other scatterings were included in isotropic scattering kernels. An approximation method, with use of the first collision source written by the I sub(i) (μ', μ) function, was introduced to attenuate the ''oscillations'' when we are obliged to use the scattering kernels with the Legendre polynomial expansion. Calculated results with this approximation showed remarkable improvement for the analysis of the angular flux spectra in a slab system of lithium metal with the D-T neutron source. (author)

Factors affecting volume calculation with single photon emission tomography (SPECT) method

International Nuclear Information System (INIS)

Liu, T.H.; Lee, K.H.; Chen, D.C.P.; Ballard, S.; Siegel, M.E.

1985-01-01

Several factors may influence the calculation of absolute volumes (VL) from SPECT images. The effect of these factors must be established to optimize the technique. The authors investigated the following on the VL calculations: % of background (BG) subtraction, reconstruction filters, sample activity, angular sampling and edge detection methods. Transaxial images of a liver-trunk phantom filled with Tc-99m from 1 to 3 μCi/cc were obtained in 64x64 matrix with a Siemens Rota Camera and MDS computer. Different reconstruction filters including Hanning 20,32, 64 and Butterworth 20, 32 were used. Angular samplings were performed in 3 and 6 degree increments. ROI's were drawn manually and with an automatic edge detection program around the image after BG subtraction. VL's were calculated by multiplying the number of pixels within the ROI by the slice thickness and the x- and y- calibrations of each pixel. One or 2 pixel per slice thickness was applied in the calculation. An inverse correlation was found between the calculated VL and the % of BG subtraction (r=0.99 for 1,2,3 μCi/cc activity). Based on the authors' linear regression analysis, the correct liver VL was measured with about 53% BG subtraction. The reconstruction filters, slice thickness and angular sampling had only minor effects on the calculated phantom volumes. Detection of the ROI automatically by the computer was not as accurate as the manual method. The authors conclude that the % of BG subtraction appears to be the most important factor affecting the VL calculation. With good quality control and appropriate reconstruction factors, correct VL calculations can be achieved with SPECT

Method to Calculate Accurate Top Event Probability in a Seismic PSA

Energy Technology Data Exchange (ETDEWEB)

Jung, Woo Sik [Sejong Univ., Seoul (Korea, Republic of)

2014-05-15

ACUBE(Advanced Cutset Upper Bound Estimator) calculates the top event probability and importance measures from cutsets by dividing cutsets into major and minor groups depending on the cutset probability, where the cutsets that have higher cutset probability are included in the major group and the others in minor cutsets, converting major cutsets into a Binary Decision Diagram (BDD). By applying the ACUBE algorithm to the seismic PSA cutsets, the accuracy of a top event probability and importance measures can be significantly improved. ACUBE works by dividing the cutsets into two groups (higher and lower cutset probability groups), calculating the top event probability and importance measures in each group, and combining the two results from the two groups. Here, ACUBE calculates the top event probability and importance measures of the higher cutset probability group exactly. On the other hand, ACUBE calculates these measures of the lower cutset probability group with an approximation such as MCUB. The ACUBE algorithm is useful for decreasing the conservatism that is caused by approximating the top event probability and importance measure calculations with given cutsets. By applying the ACUBE algorithm to the seismic PSA cutsets, the accuracy of a top event probability and importance measures can be significantly improved. This study shows that careful attention should be paid and an appropriate method be provided in order to avoid the significant overestimation of the top event probability calculation. Due to the strength of ACUBE that is explained in this study, the ACUBE became a vital tool for calculating more accurate CDF of the seismic PSA cutsets than the conventional probability calculation method.

Study on calculation methods for the effective delayed neutron fraction

International Nuclear Information System (INIS)

Irwanto, Dwi; Obara, Toru; Chiba, Go; Nagaya, Yasunobu

2011-03-01

The effective delayed neutron fraction β eff is one of the important neutronic parameters from a view point of a reactor kinetics. Several Monte-Carlo-based methods to estimate β eff have been proposed to date. In order to quantify the accuracy of these methods, we study calculation methods for β eff by analyzing various fast neutron systems including the bare spherical systems (Godiva, Jezebel, Skidoo, Jezebel-240), the reflective spherical systems (Popsy, Topsy, Flattop-23), MASURCA-R2 and MASURCA-ZONA2, and FCA XIX-1, XIX-2 and XIX-3. These analyses are performed by using SLAROM-UF and CBG for the deterministic method and MVP-II for the Monte Carlo method. We calculate β eff with various definitions such as the fundamental value β 0 , the standard definition, Nauchi's definition and Meulekamp's definition, and compare these results with each other. Through the present study, we find the following: The largest difference among the standard definition of β eff , Nauchi's β eff and Meulekamp's β eff is approximately 10%. The fundamental value β 0 is quite larger than the others in several cases. For all the cases, Meulekamp's β eff is always higher than Nauchi's β eff . This is because Nauchi's β eff considers the average neutron multiplicity value per fission which is large in the high energy range (1MeV-10MeV), while the definition of Meulekamp's β eff does not include this parameter. Furthermore, we evaluate the multi-generation effect on β eff values and demonstrate that this effect should be considered to obtain the standard definition values of β eff . (author)

Engineering calculations for the Delta S method of solving the orbital allotment problem

Science.gov (United States)

Kohnhorst, P. A.; Levis, C. A.; Walton, E. K.

1987-01-01

The method of calculating single-entry separation requirements for pairs of satellites is extended to include the interference on the top link as well as on the down link. Several heuristic models for analyzing the effects of shaped-beam antenna designs on required satellite separations are introduced and demonstrated with gain contour plots. The calculation of aggregate interference is extended to include the effects of up-link interference. The relationship between the single-entry C/I requirements, used in determining satellite separation constraints for various optimization procedures, and the aggregate C/I values of the resulting solutions is discussed.

Calculation methods for SPF for heat pump systems for comparison, system choice and dimensioning

Energy Technology Data Exchange (ETDEWEB)

Nordman, Roger; Andersson, Kajsa; Axell, Monica; Lindahl, Markus

2010-09-15

In this project, results from field measurements of heat pumps have been collected and summarised. Also existing calculation methods have been compared and summarised. Analyses have been made on how the field measurements compare to existing calculation models for heat pumps Seasonal Performance Factor (SPF), and what deviations may depend on. Recommendations for new calculation models are proposed, which include combined systems (e.g. solar - HP), capacity controlled heat pumps and combined DHW and heating operation

Three-dimensional space charge calculation method

International Nuclear Information System (INIS)

Lysenko, W.P.; Wadlinger, E.A.

1981-01-01

A method is presented for calculating space-charge forces suitable for use in a particle tracing code. Poisson's equation is solved in three dimensions with boundary conditions specified on an arbitrary surface by using a weighted residual method. Using a discrete particle distribution as our source input, examples are shown of off-axis, bunched beams of noncircular crosssection in radio-frequency quadrupole (RFQ) and drift-tube linac geometries

The development and validation of control rod calculation methods

International Nuclear Information System (INIS)

Rowlands, J.L.; Sweet, D.W.; Franklin, B.M.

1979-01-01

Fission rate distributions have been measured in the zero power critical facility, ZEBRA, for a series of eight different arrays of boron carbide control rods. Diffusion theory calculations have been compared with these measurements. The normalised fission rates differ by up to about 30% in some regions, between the different arrays, and these differences are well predicted by the calculations. A development has been made to a method used to produce homogenised cross sections for lattice regions containing control rods. Calculations show that the method also reproduces the reaction rate within the rod and the fission rate dip at the surface of the rod in satisfactory agreement with the more accurate calculations which represent the fine structure of the rod. A comparison between diffusion theory and transport theory calculations of control rod reactivity worths in the CDFR shows that for the standard design method the finite mesh approximation and the difference between diffusion theory and transport theory (the transport correction) tend to cancel and result in corrections to be applied to the standard mesh diffusion theory calculations of about +- 2% or less. This result applies for mesh centred finite difference diffusion theory codes and for the arrays of natural boron carbide control rods for which the calculations were made. Improvements have also been made to the effective diffusion coefficients used in diffusion theory calculations for control rod followers and these give satisfactory agreement with transport theory calculations. (U.K.)

Linear augmented plane wave method for self-consistent calculations

International Nuclear Information System (INIS)

Takeda, T.; Kuebler, J.

1979-01-01

O.K. Andersen has recently introduced a linear augmented plane wave method (LAPW) for the calculation of electronic structure that was shown to be computationally fast. A more general formulation of an LAPW method is presented here. It makes use of a freely disposable number of eigenfunctions of the radial Schroedinger equation. These eigenfunctions can be selected in a self-consistent way. The present formulation also results in a computationally fast method. It is shown that Andersen's LAPW is obtained in a special limit from the present formulation. Self-consistent test calculations for copper show the present method to be remarkably accurate. As an application, scalar-relativistic self-consistent calculations are presented for the band structure of FCC lanthanum. (author)

Introduction to quantum calculation methods in high resolution NMR

International Nuclear Information System (INIS)

Goldman, M.

1996-01-01

New techniques as for instance the polarization transfer, the coherence with several quanta and the double Fourier transformation have appeared fifteen years ago. These techniques constitute a considerable advance in NMR. Indeed, they allow to study more complex molecules than it was before possible. But with these advances, the classical description of the NMR is not enough to understand precisely the physical phenomena induced by these methods. It is then necessary to resort to quantum calculation methods. The aim of this work is to present these calculation methods. After some recalls of quantum mechanics, the author describes the NMR with the density matrix, reviews the main methods of double Fourier transformation and then gives the principle of the relaxation times calculation. (O.M.)

Comparison of calculational methods for EBT reactor nucleonics

International Nuclear Information System (INIS)

Henninger, R.J.; Seed, T.J.; Soran, P.D.; Dudziak, D.J.

1980-01-01

Nucleonic calculations for a preliminary conceptual design of the first wall/blanket/shield/coil assembly for an EBT reactor are described. Two-dimensional Monte Carlo, and one- and two-dimensional discrete-ordinates calculations are compared. Good agreement for the calculated values of tritium breeding and nuclear heating is seen. We find that the three methods are all useful and complementary as a design of this type evolves

Calculation method for gamma-dose rates from spherical puffs

International Nuclear Information System (INIS)

Thykier-Nielsen, S.; Deme, S.; Lang, E.

1993-05-01

The Lagrangian puff-models are widely used for calculation of the dispersion of atmospheric releases. Basic output from such models are concentrations of material in the air and on the ground. The most simple method for calculation of the gamma dose from the concentration of airborne activity is based on semi-infinite cloud model. This method is however only applicable for points far away from the release point. The exact calculation of the cloud dose using the volume integral requires significant computer time. The volume integral for the gamma dose could be approximated by using the semi-infinite cloud model combined with correction factors. This type of calculation procedure is very fast, but usually the accuracy is poor due to the fact that the same correction factors are used for all isotopes. The authors describe a more elaborate correction method. This method uses precalculated values of the gamma-dose rate as a function of the puff dispersion parameter (δ p ) and the distance from the puff centre for four energy groups. The release of energy for each radionuclide in each energy group has been calculated and tabulated. Based on these tables and a suitable interpolation procedure the calculation of gamma doses takes very short time and is almost independent of the number of radionuclides. (au) (7 tabs., 7 ills., 12 refs.)

Simple Calculation Programs for Biology Other Methods

Indian Academy of Sciences (India)

First page Back Continue Last page Overview Graphics. Simple Calculation Programs for Biology Other Methods. Hemolytic potency of drugs. Raghava et al., (1994) Biotechniques 17: 1148. FPMAP: methods for classification and identification of microorganisms 16SrRNA. graphical display of restriction and fragment map of ...

Estimation of subcriticality of TCA using 'indirect estimation method for calculation error'

International Nuclear Information System (INIS)

Naito, Yoshitaka; Yamamoto, Toshihiro; Arakawa, Takuya; Sakurai, Kiyoshi

1996-01-01

To estimate the subcriticality of neutron multiplication factor in a fissile system, 'Indirect Estimation Method for Calculation Error' is proposed. This method obtains the calculational error of neutron multiplication factor by correlating measured values with the corresponding calculated ones. This method was applied to the source multiplication and to the pulse neutron experiments conducted at TCA, and the calculation error of MCNP 4A was estimated. In the source multiplication method, the deviation of measured neutron count rate distributions from the calculated ones estimates the accuracy of calculated k eff . In the pulse neutron method, the calculation errors of prompt neutron decay constants give the accuracy of the calculated k eff . (author)

A finite element method for SSI time history calculation

International Nuclear Information System (INIS)

Ni, X.; Gantenbein, F.; Petit, M.

1989-01-01

The method which is proposed is based on a finite element modelization for the soil and the structure and a time history calculation. It has been developed for plane and axisymmetric geometries. The principle of this method is presented, then applications are given, first to a linear calculation for which results will be compared to those obtained by standard methods. Then results for a non linear behavior are described

7 CFR 51.308 - Methods of sampling and calculation of percentages.

Science.gov (United States)

2010-01-01

..., CERTIFICATION, AND STANDARDS) United States Standards for Grades of Apples Methods of Sampling and Calculation of Percentages § 51.308 Methods of sampling and calculation of percentages. (a) When the numerical... 7 Agriculture 2 2010-01-01 2010-01-01 false Methods of sampling and calculation of percentages. 51...

A tool for standardized collector performance calculations including PVT

DEFF Research Database (Denmark)

Perers, Bengt; Kovacs, Peter; Olsson, Marcus

2012-01-01

A tool for standardized calculation of solar collector performance has been developed in cooperation between SP Technical Research Institute of Sweden, DTU Denmark and SERC Dalarna University. The tool is designed to calculate the annual performance of solar collectors at representative locations...... can be tested and modeled as a thermal collector, when the PV electric part is active with an MPP tracker in operation. The thermal collector parameters from this operation mode are used for the PVT calculations....

A finite element method for SSI time history calculations

International Nuclear Information System (INIS)

Ni, X.M.; Gantenbein, F.; Petit, M.

1989-01-01

The method which is proposed is based on a finite element modelisation for the soil and the structure and a time history calculation. It has been developed for plane and axisymmetric geometries. The principle of this method will be presented, then applications will be given, first to a linear calculation for which results will be compared to those obtained by standard methods. Then results for a non linear behavior will be described

An automated Monte-Carlo based method for the calculation of cascade summing factors

Energy Technology Data Exchange (ETDEWEB)

Jackson, M.J., E-mail: mark.j.jackson@awe.co.uk; Britton, R.; Davies, A.V.; McLarty, J.L.; Goodwin, M.

2016-10-21

A versatile method has been developed to calculate cascade summing factors for use in quantitative gamma-spectrometry analysis procedures. The proposed method is based solely on Evaluated Nuclear Structure Data File (ENSDF) nuclear data, an X-ray energy library, and accurate efficiency characterisations for single detector counting geometries. The algorithm, which accounts for γ–γ, γ–X, γ–511 and γ–e{sup −} coincidences, can be applied to any design of gamma spectrometer and can be expanded to incorporate any number of nuclides. Efficiency characterisations can be derived from measured or mathematically modelled functions, and can accommodate both point and volumetric source types. The calculated results are shown to be consistent with an industry standard gamma-spectrometry software package. Additional benefits including calculation of cascade summing factors for all gamma and X-ray emissions, not just the major emission lines, are also highlighted. - Highlights: • Versatile method to calculate coincidence summing factors for gamma-spectrometry analysis. • Based solely on ENSDF format nuclear data and detector efficiency characterisations. • Enables generation of a CSF library for any detector, geometry and radionuclide. • Improves measurement accuracy and reduces acquisition times required to meet MDA.

Environment-based pin-power reconstruction method for homogeneous core calculations

International Nuclear Information System (INIS)

Leroyer, H.; Brosselard, C.; Girardi, E.

2012-01-01

Core calculation schemes are usually based on a classical two-step approach associated with assembly and core calculations. During the first step, infinite lattice assemblies calculations relying on a fundamental mode approach are used to generate cross-sections libraries for PWRs core calculations. This fundamental mode hypothesis may be questioned when dealing with loading patterns involving several types of assemblies (UOX, MOX), burnable poisons, control rods and burn-up gradients. This paper proposes a calculation method able to take into account the heterogeneous environment of the assemblies when using homogeneous core calculations and an appropriate pin-power reconstruction. This methodology is applied to MOX assemblies, computed within an environment of UOX assemblies. The new environment-based pin-power reconstruction is then used on various clusters of 3x3 assemblies showing burn-up gradients and UOX/MOX interfaces, and compared to reference calculations performed with APOLLO-2. The results show that UOX/MOX interfaces are much better calculated with the environment-based calculation scheme when compared to the usual pin-power reconstruction method. The power peak is always better located and calculated with the environment-based pin-power reconstruction method on every cluster configuration studied. This study shows that taking into account the environment in transport calculations can significantly improve the pin-power reconstruction so far as it is consistent with the core loading pattern. (authors)

An automated Monte-Carlo based method for the calculation of cascade summing factors

Science.gov (United States)

Jackson, M. J.; Britton, R.; Davies, A. V.; McLarty, J. L.; Goodwin, M.

2016-10-01

A versatile method has been developed to calculate cascade summing factors for use in quantitative gamma-spectrometry analysis procedures. The proposed method is based solely on Evaluated Nuclear Structure Data File (ENSDF) nuclear data, an X-ray energy library, and accurate efficiency characterisations for single detector counting geometries. The algorithm, which accounts for γ-γ, γ-X, γ-511 and γ-e- coincidences, can be applied to any design of gamma spectrometer and can be expanded to incorporate any number of nuclides. Efficiency characterisations can be derived from measured or mathematically modelled functions, and can accommodate both point and volumetric source types. The calculated results are shown to be consistent with an industry standard gamma-spectrometry software package. Additional benefits including calculation of cascade summing factors for all gamma and X-ray emissions, not just the major emission lines, are also highlighted.

Statistics of Monte Carlo methods used in radiation transport calculation

International Nuclear Information System (INIS)

Datta, D.

2009-01-01

Radiation transport calculation can be carried out by using either deterministic or statistical methods. Radiation transport calculation based on statistical methods is basic theme of the Monte Carlo methods. The aim of this lecture is to describe the fundamental statistics required to build the foundations of Monte Carlo technique for radiation transport calculation. Lecture note is organized in the following way. Section (1) will describe the introduction of Basic Monte Carlo and its classification towards the respective field. Section (2) will describe the random sampling methods, a key component of Monte Carlo radiation transport calculation, Section (3) will provide the statistical uncertainty of Monte Carlo estimates, Section (4) will describe in brief the importance of variance reduction techniques while sampling particles such as photon, or neutron in the process of radiation transport

Finite difference method calculations of X-ray absorption fine structure for copper

Energy Technology Data Exchange (ETDEWEB)

Bourke, J.D. [School of Physics, University of Melbourne, Parkville, Vic 3010 (Australia); Chantler, C.T. [School of Physics, University of Melbourne, Parkville, Vic 3010 (Australia)]. E-mail: chantler@physics.unimelb.edu.au; Witte, C. [School of Physics, University of Melbourne, Parkville, Vic 3010 (Australia)

2007-01-15

The finite difference method is extended to calculate X-ray absorption fine structure (XAFS) for solid state copper. These extensions include the incorporation of a Monte Carlo frozen phonon technique to simulate the effect of thermal vibrations under a correlated Debye-Waller model, and the inclusion of broadening effects from inelastic processes. Spectra are obtained over an energy range in excess of 300 eV above the K absorption edge-more than twice the greatest energy range previously reported for a solid state calculation using this method. We find this method is highly sensitive to values of the photoelectron inelastic mean free path, allowing us to probe the accuracy of current models of this parameter, particularly at low energies. We therefore find that experimental data for the photoelectron inelastic mean free path can be obtained by this method. Our results compare favourably with high precision measurements of the X-ray mass attenuation coefficient for copper, reaching agreement to within 3%, and improving previous results using the finite difference method by an order of magnitude.

Comparison between ASHRAE and ISO thermal transmittance calculation methods

DEFF Research Database (Denmark)

Blanusa, Petar; Goss, William P.; Roth, Hartwig

2007-01-01

is proportional to the glazing/frame sightline distance that is also proportional to the total glazing spacer length. An example calculation of the overall heat transfer and thermal transmittance (U-value or U-factor) using the two methods for a thermally broken, aluminum framed slider window is presented....... The fenestration thermal transmittance calculations analyses presented in this paper show that small differences exist between the calculated thermal transmittance values produced by the ISO and ASHRAE methods. The results also show that the overall thermal transmittance difference between the two methodologies...... decreases as the total window area (glazing plus frame) increases. Thus, the resulting difference in thermal transmittance values for the two methods is negligible for larger windows. This paper also shows algebraically that the differences between the ISO and ASHRAE methods turn out to be due to the way...

Hybrid SN Laplace Transform Method For Slab Lattice Calculations

International Nuclear Information System (INIS)

Segatto, Cynthia F.; Vilhena, Marco T.; Zani, Jose H.; Barros, Ricardo C.

2008-01-01

In typical lattice cells where a highly absorbing, small fuel element is embedded in the moderator, a large weakly absorbing medium, high-order transport methods become unnecessary. In this paper we describe a hybrid discrete ordinates (S N ) method for slab lattice calculations. This hybrid S N method combines the convenience of a low-order S N method in the moderator with a high-order S N method in the fuel. We use special fuel-moderator interface conditions based on an approximate angular flux interpolation analytical method and the Laplace transform (LTS N ) numerical method to calculate the neutron flux distribution and the thermal disadvantage factor. We present numerical results for a range of typical model problems. (authors)

Direct Discrete Method for Neutronic Calculations

International Nuclear Information System (INIS)

Vosoughi, Naser; Akbar Salehi, Ali; Shahriari, Majid

2002-01-01

The objective of this paper is to introduce a new direct method for neutronic calculations. This method which is named Direct Discrete Method, is simpler than the neutron Transport equation and also more compatible with physical meaning of problems. This method is based on physic of problem and with meshing of the desired geometry, writing the balance equation for each mesh intervals and with notice to the conjunction between these mesh intervals, produce the final discrete equations series without production of neutron transport differential equation and mandatory passing from differential equation bridge. We have produced neutron discrete equations for a cylindrical shape with two boundary conditions in one group energy. The correction of the results from this method are tested with MCNP-4B code execution. (authors)

Full-wave calculation of fast-wave current drive in tokamaks including kparallel upshifts

International Nuclear Information System (INIS)

Jaeger, E.F.; Batchelor, D.B.

1991-01-01

Numerical calculations of fast-wave current drive (FWCD) efficiency have generally been of two types: ray tracing or global wave calculations. Ray tracing shows that the projection of the wave number (k parallel) along the magnetic field can vary greatly over a ray trajectory, particularly when the launch point is above or below the equatorial plane. As the wave penetrates toward the center of the plasma, k parallel increases, causing a decrease in the parallel phase speed and a corresponding decrease in the current drive efficiency, γ. But the assumptions of geometrical optics, namely short wavelength and strong single-pass absorption, are not greatly applicable in FWCD scenarios. Eigenmode structure, which is ignored in ray tracing, can play an important role in determining electric field strength and Landau damping rates. In such cases, a full-wave or global solution for the wave fields is desirable. In full-wave calculations such as ORION k parallel appear as a differential operator (rvec B·∇) in the argument of the plasma dispersion function. Since this leads to a differential system of infinite order, such codes of necessity assume k parallel ∼ k var-phi = const, where k var-phi is the toroidal wave number. Thus, it is not possible to correctly include effects of the poloidal magnetic field on k parallel. The problem can be alleviated by expressing the electric field as a superposition of poloidal modes, in which case k parallel is purely algebraic. This paper describes a new full-wave calculation, Poloidal Ion Cyclotron Expansion Solution, which uses poloidal and toroidal mode expansions to solve the wave equation in general flux coordinates. The calculation includes a full solution for E parallel and uses a reduced-order form of the plasma conductivity tensor to eliminate numerical problems associated with resolution of the very short wavelength ion Bernstein wave

Calculation method of the desiccant amount for the protection of packed electronic equipments and its comparison with other methods currently used

International Nuclear Information System (INIS)

Del Valle, J.C.; Guido, O.O.

1982-01-01

The silica gel characteristics and properties are reviewed from the viewpoint of its use as a desiccant under static conditions. There were developed two calculation methods of the desiccant amount required for the protection from moisture of electronic equipments packed into rigid synthetic foam containers, for a given period and with an adequate safety margin. Besides, MIL specifications methods applicable to this type of packages have been also included. The proposed method considers the depression due to the contraction of the inside air volume during the temperature decrease periods as the only cause of the water vapor ingress into the container, the moisture permeation by diffusion through the closed cell synthetic foams being negligible. Owing to the marked thermal inertia of the package material, the thermal amplitude inside the container is lower than the outside thermal amplitude; therefore the method includes the calculation of the former. The other method is based on a procedure applicable to large storage rooms, briefly described in the literature, and where the moisture entrance is due to the above mentioned cause. A test to verify the reliability of the estimation is proposed. Finally, the amount of silica gel required for the protection of an electronic equipment packed in a given container has been calculated by the aforementioned methods, and the results compared. The value found by the proposed method was over 30% lower than that obtained by the most adequate procedure described in the MIL specification. (author) [es

A New Thermodynamic Calculation Method for Binary Alloys: Part I: Statistical Calculation of Excess Functions

Institute of Scientific and Technical Information of China (English)

无

2002-01-01

The improved form of calculation formula for the activities of the components in binary liquids and solid alloys has been derived based on the free volume theory considering excess entropy and Miedema's model for calculating the formation heat of binary alloys. A calculation method of excess thermodynamic functions for binary alloys, the formulas of integral molar excess properties and partial molar excess properties for solid ordered or disordered binary alloys have been developed. The calculated results are in good agreement with the experimental values.

Criticality calculation method for mixer-settlers

International Nuclear Information System (INIS)

Gonda, Kozo; Aoyagi, Haruki; Nakano, Ko; Kamikawa, Hiroshi.

1980-01-01

A new criticality calculation code MACPEX has been developed to evaluate and manage the criticality of the process in the extractor of mixer-settler type. MACPEX can perform the combined calculation with the PUREX process calculation code MIXSET, to get the neutron flux and the effective multiplication constant in the mixer-settlers. MACPEX solves one-dimensional diffusion equation by the explicit difference method and the standard source-iteration technique. The characteristics of MACPEX are as follows. 1) Group constants of 4 energy groups for the 239 Pu-H 2 O solution, water, polyethylene and SUS 28 are provided. 2) The group constants of the 239 Pu-H 2 O solution are given by the functional formulae of the plutonium concentration, which is less than 50 g/l. 3) Two boundary conditions of the vacuum condition and the reflective condition are available in this code. 4) The geometrical bucklings can be calculated for a certain energy group and/or region by using the three dimentional neutron flux profiles obtained by CITATION. 5) The buckling correction search can be carried out in order to get a desired k sub(eff). (author)

SU-E-T-397: Include Organ Deformation Into Dose Calculation of Prostate Brachytherapy

International Nuclear Information System (INIS)

Shao, Y; Shen, D; Chen, R; Wang, A; Lian, J

2014-01-01

Purpose: Prostate brachytherapy is an important curative treatment for patients with localized prostate cancer. In brachytherapy, rectal balloon is generally needed to adjust for unfavorable prostate position for seed placement. However, rectal balloon causes prostate deformation, which is not accounted for in dosimetric planning. Therefore, it is possible that brachytherapy dosimetry deviates significantly from initial plan when prostate returns to its non-deformed state (after procedure). The goal of this study is to develop a method to include prostate deformation into the treatment planning of brachytherapy dosimetry. Methods: We prospectively collected ultrasound images of prostate pre- and post- rectal balloon inflation from thirty five consecutive patients undergoing I-125 brachytherapy. Based on the cylinder coordinate systems, we learned the initial coordinate transformation parameters between the manual segmentations of both deformed and non-deformed prostates of each patient in training set. With the nearest-neighbor interpolation, we searched the best transformation between two coordinate systems to maximum the mutual information of deformed and non-deformed images. We then mapped the implanted seeds of five selected patients from the deformed prostate into non-deformed prostate. The seed position is marked on original pre-inflation US image and it is imported into VariSeed software for dose calculation. Results: The accuracy of image registration is 87.5% as quantified by Dice Index. The prostate coverage V100% dropped from 96.5±0.5% of prostate deformed plan to 91.9±2.6% (p<0.05) of non-deformed plan. The rectum V100% decreased from 0.44±0.26 cc to 0.10±0.18 cc (p<0.05). The dosimetry of the urethra showed mild change but not significant: V150% changed from 0.05±0.10 cc to 0.14±0.15 cc (p>0.05) and D1% changed from 212.9±37.3 Gy to 248.4±42.8 Gy (p>0.05). Conclusion: We have developed a deformable image registration method that allows

RCS Leak Rate Calculation with High Order Least Squares Method

International Nuclear Information System (INIS)

Lee, Jeong Hun; Kang, Young Kyu; Kim, Yang Ki

2010-01-01

As a part of action items for Application of Leak before Break(LBB), RCS Leak Rate Calculation Program is upgraded in Kori unit 3 and 4. For real time monitoring of operators, periodic calculation is needed and corresponding noise reduction scheme is used. This kind of study was issued in Korea, so there have upgraded and used real time RCS Leak Rate Calculation Program in UCN unit 3 and 4 and YGN unit 1 and 2. For reduction of the noise in signals, Linear Regression Method was used in those programs. Linear Regression Method is powerful method for noise reduction. But the system is not static with some alternative flow paths and this makes mixed trend patterns of input signal values. In this condition, the trend of signal and average of Linear Regression are not entirely same pattern. In this study, high order Least squares Method is used to follow the trend of signal and the order of calculation is rearranged. The result of calculation makes reasonable trend and the procedure is physically consistence

Comparison of Two-Block Decomposition Method and Chebyshev Rational Approximation Method for Depletion Calculation

International Nuclear Information System (INIS)

Lee, Yoon Hee; Cho, Nam Zin

2016-01-01

The code gives inaccurate results of nuclides for evaluation of source term analysis, e.g., Sr- 90, Ba-137m, Cs-137, etc. A Krylov Subspace method was suggested by Yamamoto et al. The method is based on the projection of solution space of Bateman equation to a lower dimension of Krylov subspace. It showed good accuracy in the detailed burnup chain calculation if dimension of the Krylov subspace is high enough. In this paper, we will compare the two methods in terms of accuracy and computing time. In this paper, two-block decomposition (TBD) method and Chebyshev rational approximation method (CRAM) are compared in the depletion calculations. In the two-block decomposition method, according to the magnitude of effective decay constant, the system of Bateman equation is decomposed into short- and longlived blocks. The short-lived block is calculated by the general Bateman solution and the importance concept. Matrix exponential with smaller norm is used in the long-lived block. In the Chebyshev rational approximation, there is no decomposition of the Bateman equation system, and the accuracy of the calculation is determined by the order of expansion in the partial fraction decomposition of the rational form. The coefficients in the partial fraction decomposition are determined by a Remez-type algorithm.

Comparison of Two-Block Decomposition Method and Chebyshev Rational Approximation Method for Depletion Calculation

Energy Technology Data Exchange (ETDEWEB)

Lee, Yoon Hee; Cho, Nam Zin [KAERI, Daejeon (Korea, Republic of)

2016-05-15

The code gives inaccurate results of nuclides for evaluation of source term analysis, e.g., Sr- 90, Ba-137m, Cs-137, etc. A Krylov Subspace method was suggested by Yamamoto et al. The method is based on the projection of solution space of Bateman equation to a lower dimension of Krylov subspace. It showed good accuracy in the detailed burnup chain calculation if dimension of the Krylov subspace is high enough. In this paper, we will compare the two methods in terms of accuracy and computing time. In this paper, two-block decomposition (TBD) method and Chebyshev rational approximation method (CRAM) are compared in the depletion calculations. In the two-block decomposition method, according to the magnitude of effective decay constant, the system of Bateman equation is decomposed into short- and longlived blocks. The short-lived block is calculated by the general Bateman solution and the importance concept. Matrix exponential with smaller norm is used in the long-lived block. In the Chebyshev rational approximation, there is no decomposition of the Bateman equation system, and the accuracy of the calculation is determined by the order of expansion in the partial fraction decomposition of the rational form. The coefficients in the partial fraction decomposition are determined by a Remez-type algorithm.

The adaptation of methods in multilayer optics for the calculation of specular neutron reflection

International Nuclear Information System (INIS)

Penfold, J.

1988-10-01

The adaptation of standard methods in multilayer optics to the calculation of specular neutron reflection is described. Their application is illustrated with examples which include a glass optical flat and a deuterated Langmuir-Blodgett film. (author)

Method of characteristics - Based sensitivity calculations for international PWR benchmark

International Nuclear Information System (INIS)

Suslov, I. R.; Tormyshev, I. V.; Komlev, O. G.

2013-01-01

Method to calculate sensitivity of fractional-linear neutron flux functionals to transport equation coefficients is proposed. Implementation of the method on the basis of MOC code MCCG3D is developed. Sensitivity calculations for fission intensity for international PWR benchmark are performed. (authors)

Research on neutron noise analysis stochastic simulation method for α calculation

International Nuclear Information System (INIS)

Zhong Bin; Shen Huayun; She Ruogu; Zhu Shengdong; Xiao Gang

2014-01-01

The prompt decay constant α has significant application on the physical design and safety analysis in nuclear facilities. To overcome the difficulty of a value calculation with Monte-Carlo method, and improve the precision, a new method based on the neutron noise analysis technology was presented. This method employs the stochastic simulation and the theory of neutron noise analysis technology. Firstly, the evolution of stochastic neutron was simulated by discrete-events Monte-Carlo method based on the theory of generalized Semi-Markov process, then the neutron noise in detectors was solved from neutron signal. Secondly, the neutron noise analysis methods such as Rossia method, Feynman-α method, zero-probability method, and cross-correlation method were used to calculate a value. All of the parameters used in neutron noise analysis method were calculated based on auto-adaptive arithmetic. The a value from these methods accords with each other, the largest relative deviation is 7.9%, which proves the feasibility of a calculation method based on neutron noise analysis stochastic simulation. (authors)

Comparison of methods for calculating decay lifetimes

International Nuclear Information System (INIS)

Tobocman, W.

1978-01-01

A simple scattering model is used to test alternative methods for calculating decay lifetimes, or equivalently, resonance widths. We consider the scattering of s-wave particles by a square well with a square barrier. Exact values for resonance energies and resonance widths are compared with values calculated from Wigner-Weisskopf perturbation theory and from the Garside-MacDonald projection operator formalism. The Garside-MacDonald formalism gives essentially exact results while the predictions of the Wigner-Weisskopf formalism are fairly poor

Structural reliability calculation method based on the dual neural network and direct integration method.

Science.gov (United States)

Li, Haibin; He, Yun; Nie, Xiaobo

2018-01-01

Structural reliability analysis under uncertainty is paid wide attention by engineers and scholars due to reflecting the structural characteristics and the bearing actual situation. The direct integration method, started from the definition of reliability theory, is easy to be understood, but there are still mathematics difficulties in the calculation of multiple integrals. Therefore, a dual neural network method is proposed for calculating multiple integrals in this paper. Dual neural network consists of two neural networks. The neural network A is used to learn the integrand function, and the neural network B is used to simulate the original function. According to the derivative relationships between the network output and the network input, the neural network B is derived from the neural network A. On this basis, the performance function of normalization is employed in the proposed method to overcome the difficulty of multiple integrations and to improve the accuracy for reliability calculations. The comparisons between the proposed method and Monte Carlo simulation method, Hasofer-Lind method, the mean value first-order second moment method have demonstrated that the proposed method is an efficient and accurate reliability method for structural reliability problems.

New method of ionization energy calculation for two-electron ions

International Nuclear Information System (INIS)

Ershov, D.K.

1997-01-01

A new method for calculation of the ionization energy of two-electron ions is proposed. The method is based on the calculation of the energy of second electron interaction with the field of an one-electron ion the potential of which is well known

Methods for tornado frequency calculation of nuclear power plant

International Nuclear Information System (INIS)

Liu Haibin; Li Lin

2012-01-01

In order to take probabilistic safety assessment of nuclear power plant tornado attack event, a method to calculate tornado frequency of nuclear power plant is introduced based on HAD 101/10 and NUREG/CR-4839 references. This method can consider history tornado frequency of the plant area, construction dimension, intensity various along with tornado path and area distribution and so on and calculate the frequency of different scale tornado. (authors)

Virial-statistic method for calculation of atom and molecule energies

International Nuclear Information System (INIS)

Borisov, Yu.A.

1977-01-01

A virial-statistical method has been applied to the calculation of the atomization energies of the following molecules: Mo(CO) 6 , Cr(CO) 6 , Fe(CO) 5 , MnH(CO) 5 , CoH(CO) 4 , Ni(CO) 4 . The principles of this method are briefly presented. Calculation results are given for the individual contributions to the atomization energies together with the calculated and experimental atomization energies (D). For the Mo(CO) 6 complex Dsub(calc) = 1759 and Dsub(exp) = 1763 kcal/mole. Calculated and experimental combination heat values for carbonyl complexes are presented. These values are shown to be adequately consistent [ru

3D electric field calculation with surface charge method

International Nuclear Information System (INIS)

Yamada, S.

1992-01-01

This paper describes an outline and some examples of three dimensional electric field calculations with a computer code developed at NIRS. In the code, a surface charge method is adopted because of it's simplicity in the mesh establishing procedure. The charge density in a triangular mesh is assumed to distribute with a linear function of the position. The electric field distribution is calculated for a pair of drift tubes with the focusing fingers on the opposing surfaces. The field distribution in an acceleration gap is analyzed with a Fourier-Bessel series expansion method. The calculated results excellently reproduces the measured data with a magnetic model. (author)

Validation of calculational methods for nuclear criticality safety - approved 1975

International Nuclear Information System (INIS)

Anon.

1977-01-01

The American National Standard for Nuclear Criticality Safety in Operations with Fissionable Materials Outside Reactors, N16.1-1975, states in 4.2.5: In the absence of directly applicable experimental measurements, the limits may be derived from calculations made by a method shown to be valid by comparison with experimental data, provided sufficient allowances are made for uncertainties in the data and in the calculations. There are many methods of calculation which vary widely in basis and form. Each has its place in the broad spectrum of problems encountered in the nuclear criticality safety field; however, the general procedure to be followed in establishing validity is common to all. The standard states the requirements for establishing the validity and area(s) of applicability of any calculational method used in assessing nuclear criticality safety

Determination of water pH using absorption-based optical sensors: evaluation of different calculation methods

Science.gov (United States)

Wang, Hongliang; Liu, Baohua; Ding, Zhongjun; Wang, Xiangxin

2017-02-01

Absorption-based optical sensors have been developed for the determination of water pH. In this paper, based on the preparation of a transparent sol-gel thin film with a phenol red (PR) indicator, several calculation methods, including simple linear regression analysis, quadratic regression analysis and dual-wavelength absorbance ratio analysis, were used to calculate water pH. Results of MSSRR show that dual-wavelength absorbance ratio analysis can improve the calculation accuracy of water pH in long-term measurement.

Calculation of neutron importance function in fissionable assemblies using Monte Carlo method

International Nuclear Information System (INIS)

Feghhi, S. A. H.; Afarideh, H.; Shahriari, M.

2007-01-01

The purpose of the present work is to develop an efficient solution method to calculate neutron importance function in fissionable assemblies for all criticality conditions, using Monte Carlo Method. The neutron importance function has a well important role in perturbation theory and reactor dynamic calculations. Usually this function can be determined by calculating adjoint flux through out solving the Adjoint weighted transport equation with deterministic methods. However, in complex geometries these calculations are very difficult. In this article, considering the capabilities of MCNP code in solving problems with complex geometries and its closeness to physical concepts, a comprehensive method based on physical concept of neutron importance has been introduced for calculating neutron importance function in sub-critical, critical and supercritical conditions. For this means a computer program has been developed. The results of the method has been benchmarked with ANISN code calculations in 1 and 2 group modes for simple geometries and their correctness has been approved for all three criticality conditions. Ultimately, the efficiency of the method for complex geometries has been shown by calculation of neutron importance in MNSR research reactor

Efficient pseudospectral methods for density functional calculations

International Nuclear Information System (INIS)

Murphy, R. B.; Cao, Y.; Beachy, M. D.; Ringnalda, M. N.; Friesner, R. A.

2000-01-01

Novel improvements of the pseudospectral method for assembling the Coulomb operator are discussed. These improvements consist of a fast atom centered multipole method and a variation of the Head-Gordan J-engine analytic integral evaluation. The details of the methodology are discussed and performance evaluations presented for larger molecules within the context of DFT energy and gradient calculations. (c) 2000 American Institute of Physics

Use of the Streaming Matrix Hybrid Method for discrete-ordinates fusion reactor calculations

International Nuclear Information System (INIS)

Battat, M.E.; Davidson, J.W.; Dudziak, D.J.; Thayer, G.R.

1984-01-01

The use of the discrete-ordinates method for solving two-dimensional, neutral-particle transport in fusion reactor blankets and shields is often limited by inherent inaccuracies due to the ray-effect. This effect presents a particular problem in the case of neutron streaming in the large internal void regions of a fusion reactor. A deterministic streaming technique called the Streaming Matrix Hybrid Method (SMHM) has been incorporated in the two-dimensional discrete-ordinates code TRIDENT-CTR. Calculations have been performed for an actual inertial-confinement fusion (ICF) reactor design using TRIDENT-CTR both with and without the SMHM. Comparisons of the calculated fluxes indicate that substantial mitigation of the ray effect can be achieved with the SMHM. Calculations were performed for the Los Alamos FIRST STEP hybrid ICF reactor designed for tritium production. Conventional 238 U fuel rod assemblies surround the spherical steel target chamber to form an annular cylindrical blanket. An axial fuel region is included to complete the blanket

Biostatistics with emphasis on life table survival rate calculations (including Kaplan Meier) and the logrank test

International Nuclear Information System (INIS)

Mould, Richard F.

1995-01-01

Purpose/Objective: To explain some of the most useful statistical calculation procedures which are relevant to radiation oncologists and to provide insights on what tests and procedures should be used in various situations such as when survival rates and their associated standard errors have to be determined. To describe some of the problems and pitfalls in clinical trial designs which have to be overcome if a trial is to have the possibility of reaching a successful conclusion. To review methods of computing criteria to quantitatively describe criteria of success (eg. quality of life, long-term survival, cure) of radiation oncology and to suggest possible future statistical improvements in this area. Chi-Squared Test: The chi-squared test is probably the most useful of the tests of statistical significance for the radiation oncologist. Applications will be described, including goodness of fit tests and 2x2 contingency tables which are the simplest of the generalized nxm contingency tables. Degrees of Freedom and P<0.05 for Significance Testing: An Introduction will be given to the meaning of P<0.05 in relation to significance testing and the use of tables of critical values of a test statistic (eg. chi-squared) which are given as a function of degrees of freedom and P-values. Survival Rate Calculations for Grouped and Ungrouped Data: The life-table method (sometimes termed the actuarial method) will be explained for both grouped data (eg. survival times grouped in annual intervals for patients who have died and for those who are still alive or lost to follow-up) and for ungrouped data (when individual survival times are used). The method for ungrouped data is variously termed the Kaplan-Meier or Product Limit method. Logrank Test: This is the most useful test for comparison of the survival experience of two groups of patients and its use will be explained. In part the computation is similar to that for the Kaplan-Meier/Product Limit method

Biostatistics with emphasis on life table survival rate calculations (including Kaplan Meier) and the logrank test

Energy Technology Data Exchange (ETDEWEB)

Mould, Richard F

1995-07-01

Purpose/Objective: To explain some of the most useful statistical calculation procedures which are relevant to radiation oncologists and to provide insights on what tests and procedures should be used in various situations such as when survival rates and their associated standard errors have to be determined. To describe some of the problems and pitfalls in clinical trial designs which have to be overcome if a trial is to have the possibility of reaching a successful conclusion. To review methods of computing criteria to quantitatively describe criteria of success (eg. quality of life, long-term survival, cure) of radiation oncology and to suggest possible future statistical improvements in this area. Chi-Squared Test: The chi-squared test is probably the most useful of the tests of statistical significance for the radiation oncologist. Applications will be described, including goodness of fit tests and 2x2 contingency tables which are the simplest of the generalized nxm contingency tables. Degrees of Freedom and P<0.05 for Significance Testing: An Introduction will be given to the meaning of P<0.05 in relation to significance testing and the use of tables of critical values of a test statistic (eg. chi-squared) which are given as a function of degrees of freedom and P-values. Survival Rate Calculations for Grouped and Ungrouped Data: The life-table method (sometimes termed the actuarial method) will be explained for both grouped data (eg. survival times grouped in annual intervals for patients who have died and for those who are still alive or lost to follow-up) and for ungrouped data (when individual survival times are used). The method for ungrouped data is variously termed the Kaplan-Meier or Product Limit method. Logrank Test: This is the most useful test for comparison of the survival experience of two groups of patients and its use will be explained. In part the computation is similar to that for the Kaplan-Meier/Product Limit method.

High order aberrations calculation of a hexapole corrector using a differential algebra method

Energy Technology Data Exchange (ETDEWEB)

Kang, Yongfeng, E-mail: yfkang@mail.xjtu.edu.cn [Key Laboratory for Physical Electronics and Devices of the Ministry of Education, Xi' an Jiaotong University, Xi' an 710049 (China); Liu, Xing [Key Laboratory for Physical Electronics and Devices of the Ministry of Education, Xi' an Jiaotong University, Xi' an 710049 (China); Zhao, Jingyi, E-mail: jingyi.zhao@foxmail.com [School of Science, Chang’an University, Xi’an 710064 (China); Tang, Tiantong [Key Laboratory for Physical Electronics and Devices of the Ministry of Education, Xi' an Jiaotong University, Xi' an 710049 (China)

2017-02-21

A differential algebraic (DA) method is proved as an unusual and effective tool in numerical analysis. It implements conveniently differentiation up to arbitrary high order, based on the nonstandard analysis. In this paper, the differential algebra (DA) method has been employed to compute the high order aberrations up to the fifth order of a practical hexapole corrector including round lenses and hexapole lenses. The program has been developed and tested as well. The electro-magnetic fields of arbitrary point are obtained by local analytic expressions, then field potentials are transformed into new forms which can be operated in the DA calculation. In this paper, the geometric and chromatic aberrations up to fifth order of a practical hexapole corrector system are calculated by the developed program.

Resampling Approach for Determination of the Method for Reference Interval Calculation in Clinical Laboratory Practice▿

Science.gov (United States)

Pavlov, Igor Y.; Wilson, Andrew R.; Delgado, Julio C.

2010-01-01

Reference intervals (RI) play a key role in clinical interpretation of laboratory test results. Numerous articles are devoted to analyzing and discussing various methods of RI determination. The two most widely used approaches are the parametric method, which assumes data normality, and a nonparametric, rank-based procedure. The decision about which method to use is usually made arbitrarily. The goal of this study was to demonstrate that using a resampling approach for the comparison of RI determination techniques could help researchers select the right procedure. Three methods of RI calculation—parametric, transformed parametric, and quantile-based bootstrapping—were applied to multiple random samples drawn from 81 values of complement factor B observations and from a computer-simulated normally distributed population. It was shown that differences in RI between legitimate methods could be up to 20% and even more. The transformed parametric method was found to be the best method for the calculation of RI of non-normally distributed factor B estimations, producing an unbiased RI and the lowest confidence limits and interquartile ranges. For a simulated Gaussian population, parametric calculations, as expected, were the best; quantile-based bootstrapping produced biased results at low sample sizes, and the transformed parametric method generated heavily biased RI. The resampling approach could help compare different RI calculation methods. An algorithm showing a resampling procedure for choosing the appropriate method for RI calculations is included. PMID:20554803

Some experience of shielding calculations by combinatorial method

International Nuclear Information System (INIS)

Korobejnikov, V.V.; Oussanov, V.I.

1996-01-01

Some aspects of the compound systems shielding calculations by a combinatorial approach are discussed. The effectiveness of such an approach is based on the fundamental characteristic of a compound system: if some element of the system have in itself mathematical or physical properties favorable for calculation, these properties may be used in a combinatorial approach and are lost when the system is being calculated in the whole by a direct approach. The combinatorial technique applied is well known. A compound system are being splitting for two or more auxiliary subsystems (so that calculation each of them is a more simple problem than calculation of the original problem (or at last is a soluble problem if original one is not). Calculation of every subsystem are carried out by suitable method and code, the coupling being made through boundary conditions or boundary source. The special consideration in the paper is given to a fast reactor shielding combinatorial analysis and to the testing of the results received. (author)

A method for the computation of turbulent polymeric liquids including hydrodynamic interactions and chain entanglements

Energy Technology Data Exchange (ETDEWEB)

Kivotides, Demosthenes, E-mail: demosthenes.kivotides@strath.ac.uk

2017-02-12

An asymptotically exact method for the direct computation of turbulent polymeric liquids that includes (a) fully resolved, creeping microflow fields due to hydrodynamic interactions between chains, (b) exact account of (subfilter) residual stresses, (c) polymer Brownian motion, and (d) direct calculation of chain entanglements, is formulated. Although developed in the context of polymeric fluids, the method is equally applicable to turbulent colloidal dispersions and aerosols. - Highlights: • An asymptotically exact method for the computation of polymer and colloidal fluids is developed. • The method is valid for all flow inertia and all polymer volume fractions. • The method models entanglements and hydrodynamic interactions between polymer chains.

A New Power Calculation Method for Single-Phase Grid-Connected Systems

DEFF Research Database (Denmark)

Yang, Yongheng; Blaabjerg, Frede

2013-01-01

A new method to calculate average active power and reactive power for single-phase systems is proposed in this paper. It can be used in different applications where the output active power and reactive power need to be calculated accurately and fast. For example, a grid-connected photovoltaic...... system in low voltage ride through operation mode requires a power feedback for the power control loop. Commonly, a Discrete Fourier Transform (DFT) based power calculation method can be adopted in such systems. However, the DFT method introduces at least a one-cycle time delay. The new power calculation...... method, which is based on the adaptive filtering technique, can achieve a faster response. The performance of the proposed method is verified by experiments and demonstrated in a 1 kW single-phase grid-connected system operating under different conditions.Experimental results show the effectiveness...

Earthquake analysis of structures including structure-soil interaction by a substructure method

International Nuclear Information System (INIS)

Chopra, A.K.; Guttierrez, J.A.

1977-01-01

A general substructure method for analysis of response of nuclear power plant structures to earthquake ground motion, including the effects of structure-soil interaction, is summarized. The method is applicable to complex structures idealized as finite element systems and the soil region treated as either a continuum, for example as a viscoelastic halfspace, or idealized as a finite element system. The halfspace idealization permits reliable analysis for sites where essentially similar soils extend to large depths and there is no rigid boundary such as soil-rock interface. For sites where layers of soft soil are underlain by rock at shallow depth, finite element idealization of the soil region is appropriate; in this case, the direct and substructure methods would lead to equivalent results but the latter provides the better alternative. Treating the free field motion directly as the earthquake input in the substructure eliminates the deconvolution calculations and the related assumption-regarding type and direction of earthquake waves-required in the direct method. (Auth.)

Perturbation method for calculation of narrow-band impedance and trapped modes

International Nuclear Information System (INIS)

Heifets, S.A.

1987-01-01

An iterative method for calculation of the narrow-band impedance is described for a system with a small variation in boundary conditions, so that the variation can be considered as a perturbation. The results are compared with numeric calculations. The method is used to relate the origin of the trapped modes with the degeneracy of the spectrum of an unperturbed system. The method also can be applied to transverse impedance calculations. 6 refs., 6 figs., 1 tab

Comparison of matrix exponential methods for fuel burnup calculations

International Nuclear Information System (INIS)

Oh, Hyung Suk; Yang, Won Sik

1999-01-01

Series expansion methods to compute the exponential of a matrix have been compared by applying them to fuel depletion calculations. Specifically, Taylor, Pade, Chebyshev, and rational Chebyshev approximations have been investigated by approximating the exponentials of bum matrices by truncated series of each method with the scaling and squaring algorithm. The accuracy and efficiency of these methods have been tested by performing various numerical tests using one thermal reactor and two fast reactor depletion problems. The results indicate that all the four series methods are accurate enough to be used for fuel depletion calculations although the rational Chebyshev approximation is relatively less accurate. They also show that the rational approximations are more efficient than the polynomial approximations. Considering the computational accuracy and efficiency, the Pade approximation appears to be better than the other methods. Its accuracy is better than the rational Chebyshev approximation, while being comparable to the polynomial approximations. On the other hand, its efficiency is better than the polynomial approximations and is similar to the rational Chebyshev approximation. In particular, for fast reactor depletion calculations, it is faster than the polynomial approximations by a factor of ∼ 1.7. (author). 11 refs., 4 figs., 2 tabs

Correlation expansion: a powerful alternative multiple scattering calculation method

International Nuclear Information System (INIS)

Zhao Haifeng; Wu Ziyu; Sebilleau, Didier

2008-01-01

We introduce a powerful alternative expansion method to perform multiple scattering calculations. In contrast to standard MS series expansion, where the scattering contributions are grouped in terms of scattering order and may diverge in the low energy region, this expansion, called correlation expansion, partitions the scattering process into contributions from different small atom groups and converges at all energies. It converges faster than MS series expansion when the latter is convergent. Furthermore, it takes less memory than the full MS method so it can be used in the near edge region without any divergence problem, even for large clusters. The correlation expansion framework we derive here is very general and can serve to calculate all the elements of the scattering path operator matrix. Photoelectron diffraction calculations in a cluster containing 23 atoms are presented to test the method and compare it to full MS and standard MS series expansion

A numerical method for resonance integral calculations

International Nuclear Information System (INIS)

Tanbay, Tayfun; Ozgener, Bilge

2013-01-01

A numerical method has been proposed for resonance integral calculations and a cubic fit based on least squares approximation to compute the optimum Bell factor is given. The numerical method is based on the discretization of the neutron slowing down equation. The scattering integral is approximated by taking into account the location of the upper limit in energy domain. The accuracy of the method has been tested by performing computations of resonance integrals for uranium dioxide isolated rods and comparing the results with empirical values. (orig.)

Simplified dose calculation method for mantle technique

International Nuclear Information System (INIS)

Scaff, L.A.M.

1984-01-01

A simplified dose calculation method for mantle technique is described. In the routine treatment of lymphom as using this technique, the daily doses at the midpoints at five anatomical regions are different because the thicknesses are not equal. (Author) [pt

A fast dose calculation method based on table lookup for IMRT optimization

International Nuclear Information System (INIS)

Wu Qiuwen; Djajaputra, David; Lauterbach, Marc; Wu Yan; Mohan, Radhe

2003-01-01

This note describes a fast dose calculation method that can be used to speed up the optimization process in intensity-modulated radiotherapy (IMRT). Most iterative optimization algorithms in IMRT require a large number of dose calculations to achieve convergence and therefore the total amount of time needed for the IMRT planning can be substantially reduced by using a faster dose calculation method. The method that is described in this note relies on an accurate dose calculation engine that is used to calculate an approximate dose kernel for each beam used in the treatment plan. Once the kernel is computed and saved, subsequent dose calculations can be done rapidly by looking up this kernel. Inaccuracies due to the approximate nature of the kernel in this method can be reduced by performing scheduled kernel updates. This fast dose calculation method can be performed more than two orders of magnitude faster than the typical superposition/convolution methods and therefore is suitable for applications in which speed is critical, e.g., in an IMRT optimization that requires a simulated annealing optimization algorithm or in a practical IMRT beam-angle optimization system. (note)

Point kernels and superposition methods for scatter dose calculations in brachytherapy

International Nuclear Information System (INIS)

Carlsson, A.K.

2000-01-01

Point kernels have been generated and applied for calculation of scatter dose distributions around monoenergetic point sources for photon energies ranging from 28 to 662 keV. Three different approaches for dose calculations have been compared: a single-kernel superposition method, a single-kernel superposition method where the point kernels are approximated as isotropic and a novel 'successive-scattering' superposition method for improved modelling of the dose from multiply scattered photons. An extended version of the EGS4 Monte Carlo code was used for generating the kernels and for benchmarking the absorbed dose distributions calculated with the superposition methods. It is shown that dose calculation by superposition at and below 100 keV can be simplified by using isotropic point kernels. Compared to the assumption of full in-scattering made by algorithms currently in clinical use, the single-kernel superposition method improves dose calculations in a half-phantom consisting of air and water. Further improvements are obtained using the successive-scattering superposition method, which reduces the overestimates of dose close to the phantom surface usually associated with kernel superposition methods at brachytherapy photon energies. It is also shown that scatter dose point kernels can be parametrized to biexponential functions, making them suitable for use with an effective implementation of the collapsed cone superposition algorithm. (author)

Calculation of neutron importance function in fissionable assemblies using Monte Carlo method

International Nuclear Information System (INIS)

Feghhi, S.A.H.; Shahriari, M.; Afarideh, H.

2007-01-01

The purpose of the present work is to develop an efficient solution method for the calculation of neutron importance function in fissionable assemblies for all criticality conditions, based on Monte Carlo calculations. The neutron importance function has an important role in perturbation theory and reactor dynamic calculations. Usually this function can be determined by calculating the adjoint flux while solving the adjoint weighted transport equation based on deterministic methods. However, in complex geometries these calculations are very complicated. In this article, considering the capabilities of MCNP code in solving problems with complex geometries and its closeness to physical concepts, a comprehensive method based on the physical concept of neutron importance has been introduced for calculating the neutron importance function in sub-critical, critical and super-critical conditions. For this propose a computer program has been developed. The results of the method have been benchmarked with ANISN code calculations in 1 and 2 group modes for simple geometries. The correctness of these results has been confirmed for all three criticality conditions. Finally, the efficiency of the method for complex geometries has been shown by the calculation of neutron importance in Miniature Neutron Source Reactor (MNSR) research reactor

Weather data for simplified energy calculation methods. Volume IV. United States: WYEC data

Energy Technology Data Exchange (ETDEWEB)

Olsen, A.R.; Moreno, S.; Deringer, J.; Watson, C.R.

1984-08-01

The objective of this report is to provide a source of weather data for direct use with a number of simplified energy calculation methods available today. Complete weather data for a number of cities in the United States are provided for use in the following methods: degree hour, modified degree hour, bin, modified bin, and variable degree day. This report contains sets of weather data for 23 cities using Weather Year for Energy Calculations (WYEC) source weather data. Considerable overlap is present in cities (21) covered by both the TRY and WYEC data. The weather data at each city has been summarized in a number of ways to provide differing levels of detail necessary for alternative simplified energy calculation methods. Weather variables summarized include dry bulb and wet bulb temperature, percent relative humidity, humidity ratio, wind speed, percent possible sunshine, percent diffuse solar radiation, total solar radiation on horizontal and vertical surfaces, and solar heat gain through standard DSA glass. Monthly and annual summaries, in some cases by time of day, are available. These summaries are produced in a series of nine computer generated tables.

Influence of Individual Differences on the Calculation Method for FBG-Type Blood Pressure Sensors.

Science.gov (United States)

Koyama, Shouhei; Ishizawa, Hiroaki; Fujimoto, Keisaku; Chino, Shun; Kobayashi, Yuka

2016-12-28

In this paper, we propose a blood pressure calculation and associated measurement method that by using a fiber Bragg grating (FBG) sensor. There are several points at which the pulse can be measured on the surface of the human body, and when a FBG sensor located at any of these points, the pulse wave signal can be measured. The measured waveform is similar to the acceleration pulse wave. The pulse wave signal changes depending on several factors, including whether or not the individual is healthy and/or elderly. The measured pulse wave signal can be used to calculate the blood pressure using a calibration curve, which is constructed by a partial least squares (PLS) regression analysis using a reference blood pressure and the pulse wave signal. In this paper, we focus on the influence of individual differences from calculated blood pressure based on each calibration curve. In our study, the calculated blood pressure from both the individual and overall calibration curves were compared, and our results show that the calculated blood pressure based on the overall calibration curve had a lower measurement accuracy than that based on an individual calibration curve. We also found that the influence of the individual differences on the calculated blood pressure when using the FBG sensor method were very low. Therefore, the FBG sensor method that we developed for measuring the blood pressure was found to be suitable for use by many people.

Research of coincidence method for calculation model of the specific detector

Energy Technology Data Exchange (ETDEWEB)

Guangchun, Hu; Suping, Liu; Jian, Gong [China Academy of Engineering Physics, Mianyang (China). Inst. of Nuclear Physics and Chemistry

2003-07-01

The physical size of specific detector is known normally, but production business is classified for some sizes that is concerned with the property of detector, such as the well diameter, well depth of detector and dead region. The surface source of even distribution and the sampling method of source particle isotropy sport have been established with the method of Monte Carlo, and gamma ray respond spectral with the {sup 152}Eu surface source been calculated. The experiment have been performed under the same conditions. Calculation and experiment results are compared with relative efficiency coincidence method and spectral similar degree coincidence method. According to comparison as a result, detector model is revised repeatedly to determine the calculation model of detector and to calculate efficiency of detector and spectra. (authors)

Reliable method for fission source convergence of Monte Carlo criticality calculation with Wielandt's method

International Nuclear Information System (INIS)

Yamamoto, Toshihiro; Miyoshi, Yoshinori

2004-01-01

A new algorithm of Monte Carlo criticality calculations for implementing Wielandt's method, which is one of acceleration techniques for deterministic source iteration methods, is developed, and the algorithm can be successfully implemented into MCNP code. In this algorithm, part of fission neutrons emitted during random walk processes are tracked within the current cycle, and thus a fission source distribution used in the next cycle spread more widely. Applying this method intensifies a neutron interaction effect even in a loosely-coupled array where conventional Monte Carlo criticality methods have difficulties, and a converged fission source distribution can be obtained with fewer cycles. Computing time spent for one cycle, however, increases because of tracking fission neutrons within the current cycle, which eventually results in an increase of total computing time up to convergence. In addition, statistical fluctuations of a fission source distribution in a cycle are worsened by applying Wielandt's method to Monte Carlo criticality calculations. However, since a fission source convergence is attained with fewer source iterations, a reliable determination of convergence can easily be made even in a system with a slow convergence. This acceleration method is expected to contribute to prevention of incorrect Monte Carlo criticality calculations. (author)

Overview of multifluid-flow-calculation methods

International Nuclear Information System (INIS)

Stewart, H.B.

1981-01-01

Two categories of numerical methods which may be useful in multiphase flow research are discussed. The first category includes methods which are specifically intended for accurate computation of discontinuities, such as the method of characteristics, particle-in-cell method, flux-corrected transport, and random choice methods. Methods in this category could be applied to research on rocket exhaust plumes and interior ballistics. The second category includes methods for smooth, subsonic flows, such as fractional step methods, semi-implicit method, and methods which treat convection implicitly. The subsonic flow methods could be of interest for ice flows

Modifications of the Fourier approach for magnetic field calculations to include axial shields in superconducting magnets

International Nuclear Information System (INIS)

Caldwell, J.

1984-01-01

Martinelli and Morini have used an analytical method for calculating values and distribution of the magnetic field in superconducting magnets. Using Fourier series the magnetic field is determined by carrying out a series expansion of the current density distribution of the system of coils. This Fourier method can be modified to include axial iron to a far greater accuracy (for finite permeability) by incorporating the image series approach of Caldwell and Zisserman. Also an exact solution can be obtained for the case of infinite permeability. A comparison of the results derived from the expansion of Martinelli and Morini with the exact solution of Caldwell and Zisserman shows excellent agreement for the iron-free case but the accuracy deteriorates as the permeability μ/sub z/ increases. The exact solution should be used for infinite permeability and also gives satisfactory results for permeability μ/sub z/ >100. A symmetric geometry is used throughout the communication for simplicity of presentation

Proposed method to calculate FRMAC intervention levels for the assessment of radiologically contaminated food and comparison of the proposed method to the U.S. FDA's method to calculate derived intervention levels

Energy Technology Data Exchange (ETDEWEB)

Kraus, Terrence D.; Hunt, Brian D.

2014-02-01

This report reviews the method recommended by the U.S. Food and Drug Administration for calculating Derived Intervention Levels (DILs) and identifies potential improvements to the DIL calculation method to support more accurate ingestion pathway analyses and protective action decisions. Further, this report proposes an alternate method for use by the Federal Emergency Radiological Assessment Center (FRMAC) to calculate FRMAC Intervention Levels (FILs). The default approach of the FRMAC during an emergency response is to use the FDA recommended methods. However, FRMAC recommends implementing the FIL method because we believe it to be more technically accurate. FRMAC will only implement the FIL method when approved by the FDA representative on the Federal Advisory Team for Environment, Food, and Health.

An Efficient Method for Electron-Atom Scattering Using Ab-initio Calculations

Energy Technology Data Exchange (ETDEWEB)

Xu, Yuan; Yang, Yonggang; Xiao, Liantuan; Jia, Suotang [Shanxi University, Taiyuan (China)

2017-02-15

We present an efficient method based on ab-initio calculations to investigate electron-atom scatterings. Those calculations profit from methods implemented in standard quantum chemistry programs. The new approach is applied to electron-helium scattering. The results are compared with experimental and other theoretical references to demonstrate the efficiency of our method.

Application of nonparametric statistic method for DNBR limit calculation

International Nuclear Information System (INIS)

Dong Bo; Kuang Bo; Zhu Xuenong

2013-01-01

Background: Nonparametric statistical method is a kind of statistical inference method not depending on a certain distribution; it calculates the tolerance limits under certain probability level and confidence through sampling methods. The DNBR margin is one important parameter of NPP design, which presents the safety level of NPP. Purpose and Methods: This paper uses nonparametric statistical method basing on Wilks formula and VIPER-01 subchannel analysis code to calculate the DNBR design limits (DL) of 300 MW NPP (Nuclear Power Plant) during the complete loss of flow accident, simultaneously compared with the DL of DNBR through means of ITDP to get certain DNBR margin. Results: The results indicate that this method can gain 2.96% DNBR margin more than that obtained by ITDP methodology. Conclusions: Because of the reduction of the conservation during analysis process, the nonparametric statistical method can provide greater DNBR margin and the increase of DNBR margin is benefited for the upgrading of core refuel scheme. (authors)

Methods for calculating the speed-up characteristics of steam-water turbines

International Nuclear Information System (INIS)

Golovach, E.A.

1981-01-01

The methods of approximate and specified calculations of speed- up characteristics of steam-water turbines are considered. The specified non-linear method takes into account change of thermal efficiency, heat drop and losses in the turbine as well as vacuum break-up the condenser. Speed-up characteristics of the K-1000-60-1500 turbine are presented. The calculational results obtained by the non-linear method are compared with the calculations conducted by the approximate linearized method. Differences in the frequency speed up of the turbine rotor rotation calculated by the two methods constitute only 0.5-2.0%. That is why it is necessary to take into account in the specified calculations first of all the most important factors following the rotor speed- up in the following consequence: valve shift of the high pressure cylinder (HPC); steam volume in front of the HPC; shift of the valves behind the separator-steam superheater (SSS); steam volumes and moisture boiling in the SSS; steam consumption for regenerating heating of feed water, steam volumes at the intermediate elements of the turbine, losses in the turbine, heat drop and thermal efficiency [ru

Standard Guide for Selection and Use of Mathematical Methods for Calculating Absorbed Dose in Radiation Processing Applications

CERN Document Server

American Society for Testing and Materials. Philadelphia

2010-01-01

1.1 This guide describes different mathematical methods that may be used to calculate absorbed dose and criteria for their selection. Absorbed-dose calculations can determine the effectiveness of the radiation process, estimate the absorbed-dose distribution in product, or supplement or complement, or both, the measurement of absorbed dose. 1.2 Radiation processing is an evolving field and annotated examples are provided in Annex A6 to illustrate the applications where mathematical methods have been successfully applied. While not limited by the applications cited in these examples, applications specific to neutron transport, radiation therapy and shielding design are not addressed in this document. 1.3 This guide covers the calculation of radiation transport of electrons and photons with energies up to 25 MeV. 1.4 The mathematical methods described include Monte Carlo, point kernel, discrete ordinate, semi-empirical and empirical methods. 1.5 General purpose software packages are available for the calcul...

PKI, Gamma Radiation Reactor Shielding Calculation by Point-Kernel Method

International Nuclear Information System (INIS)

Li Chunhuai; Zhang Liwu; Zhang Yuqin; Zhang Chuanxu; Niu Xihua

1990-01-01

1 - Description of program or function: This code calculates radiation shielding problem of gamma-ray in geometric space. 2 - Method of solution: PKI uses a point kernel integration technique, describes radiation shielding geometric space by using geometric space configuration method and coordinate conversion, and makes use of calculation result of reactor primary shielding and flow regularity in loop system for coolant

On the resolvents methods in quantum perturbation calculations

International Nuclear Information System (INIS)

Burzynski, A.

1979-01-01

This paper gives a systematic review of resolvent methods in quantum perturbation calculations. The case of discrete spectrum of hamiltonian is considered specially (in the literature this is the fewest considered case). The topics of calculations of quantum transitions by using of the resolvent formalism, quantum transitions between states from particular subspaces, the shifts of energy levels, are shown. The main ideas of stationary perturbation theory developed by Lippmann and Schwinger are considered too. (author)

Ca(AlH4)2, CaAlH5, and CaH2+6LiBH4: Calculated dehydrogenation enthalpy, including zero point energy, and the structure of the phonon spectra.

Science.gov (United States)

Marashdeh, Ali; Frankcombe, Terry J

2008-06-21

The dehydrogenation enthalpies of Ca(AlH(4))(2), CaAlH(5), and CaH(2)+6LiBH(4) have been calculated using density functional theory calculations at the generalized gradient approximation level. Harmonic phonon zero point energy (ZPE) corrections have been included using Parlinski's direct method. The dehydrogenation of Ca(AlH(4))(2) is exothermic, indicating a metastable hydride. Calculations for CaAlH(5) including ZPE effects indicate that it is not stable enough for a hydrogen storage system operating near ambient conditions. The destabilized combination of LiBH(4) with CaH(2) is a promising system after ZPE-corrected enthalpy calculations. The calculations confirm that including ZPE effects in the harmonic approximation for the dehydrogenation of Ca(AlH(4))(2), CaAlH(5), and CaH(2)+6LiBH(4) has a significant effect on the calculated reaction enthalpy. The contribution of ZPE to the dehydrogenation enthalpies of Ca(AlH(4))(2) and CaAlH(5) calculated by the direct method phonon analysis was compared to that calculated by the frozen-phonon method. The crystal structure of CaAlH(5) is presented in the more useful standard setting of P2(1)c symmetry and the phonon density of states of CaAlH(5), significantly different to other common complex metal hydrides, is rationalized.

Practical Calculational Scheme Implementing the Wilsonian RG Results for Nuclear Effective Field Theory Including Pions

International Nuclear Information System (INIS)

Kubo, H.; Harada, K.; Sakaeda, T.; Yamamoto, Y.

2013-01-01

On the basis of the Wilsonian renormalization group (WRG) analysis of nuclear effective field theory (NEFT) including pions, we propose a practical calculational scheme in which the short-distance part of one-pion exchange (S-OPE) is removed and represented as contact terms. The long-distance part of one-pion exchange (L-OPE) is treated as perturbation. The use of dimensional regularization (DR) for diagrams consisting only of contact interactions considerably simplifies the calculation of scattering amplitude and the renormalization group equations. NLO results for nucleon-nucleon elastic scattering in the S-waves are obtained and compared with experiments. A brief comment on NNLO calculations is given. (author)

Large-scale atomic calculations using variational methods

International Nuclear Information System (INIS)

Joensson, Per.

1995-01-01

Atomic properties, such as radiative lifetimes, hyperfine structures and isotope shift, have been studied both theoretically and experimentally. Computer programs which calculate these properties from multiconfiguration Hartree-Fock (MCHF) and configuration interaction (CI) wave functions have been developed and tested. To study relativistic effects, a program which calculates hyperfine structures from multiconfiguration Dirac-Fock (MCDF) wave functions has also been written. A new method of dealing with radial non-orthogonalities in transition matrix elements has been investigated. This method allows two separate orbital sets to be used for the initial and final states, respectively. It is shown that, once the usual orthogonality restrictions have been overcome, systematic MCHF calculations are able to predict oscillator strengths in light atoms with high accuracy. In connection with recent high-power laser experiments, time-dependent calculations of the atomic response to intense laser fields have been performed. Using the frozen-core approximation, where the atom is modeled as an active electron moving in the average field of the core electrons and the nucleus, the active electron has been propagated in time under the influence of the laser field. Radiative lifetimes and hyperfine structures of excited states in sodium and silver have been experimentally determined using time-resolved laser spectroscopy. By recording the fluorescence light decay following laser excitation in the vacuum ultraviolet spectral region, the radiative lifetimes and hyperfine structures of the 7p 2 P states in silver have been measured. The delayed-coincidence technique has been used to make very accurate measurements of the radiative lifetimes and hyperfine structures of the lowest 2P states in sodium and silver. 77 refs, 2 figs, 14 tabs

Large-scale atomic calculations using variational methods

Energy Technology Data Exchange (ETDEWEB)

Joensson, Per

1995-01-01

Atomic properties, such as radiative lifetimes, hyperfine structures and isotope shift, have been studied both theoretically and experimentally. Computer programs which calculate these properties from multiconfiguration Hartree-Fock (MCHF) and configuration interaction (CI) wave functions have been developed and tested. To study relativistic effects, a program which calculates hyperfine structures from multiconfiguration Dirac-Fock (MCDF) wave functions has also been written. A new method of dealing with radial non-orthogonalities in transition matrix elements has been investigated. This method allows two separate orbital sets to be used for the initial and final states, respectively. It is shown that, once the usual orthogonality restrictions have been overcome, systematic MCHF calculations are able to predict oscillator strengths in light atoms with high accuracy. In connection with recent high-power laser experiments, time-dependent calculations of the atomic response to intense laser fields have been performed. Using the frozen-core approximation, where the atom is modeled as an active electron moving in the average field of the core electrons and the nucleus, the active electron has been propagated in time under the influence of the laser field. Radiative lifetimes and hyperfine structures of excited states in sodium and silver have been experimentally determined using time-resolved laser spectroscopy. By recording the fluorescence light decay following laser excitation in the vacuum ultraviolet spectral region, the radiative lifetimes and hyperfine structures of the 7p{sup 2}P states in silver have been measured. The delayed-coincidence technique has been used to make very accurate measurements of the radiative lifetimes and hyperfine structures of the lowest 2P states in sodium and silver. 77 refs, 2 figs, 14 tabs.

Cluster monte carlo method for nuclear criticality safety calculation

International Nuclear Information System (INIS)

Pei Lucheng

1984-01-01

One of the most important applications of the Monte Carlo method is the calculation of the nuclear criticality safety. The fair source game problem was presented at almost the same time as the Monte Carlo method was applied to calculating the nuclear criticality safety. The source iteration cost may be reduced as much as possible or no need for any source iteration. This kind of problems all belongs to the fair source game prolems, among which, the optimal source game is without any source iteration. Although the single neutron Monte Carlo method solved the problem without the source iteration, there is still quite an apparent shortcoming in it, that is, it solves the problem without the source iteration only in the asymptotic sense. In this work, a new Monte Carlo method called the cluster Monte Carlo method is given to solve the problem further

Polarizable Embedded RI-CC2 Method for Two-Photon Absorption Calculations

DEFF Research Database (Denmark)

Hršak, Dalibor; Khah, Alireza Marefat; Christiansen, Ove

2015-01-01

We present a novel polarizable embedded resolution-of-identity coupled cluster singles and approximate doubles (PERI-CC2) method for calculation of two-photon absorption (TPA) spectra of large molecular systems. The method was benchmarked for three types of systems: a water-solvated molecule...... of formamide, a uracil molecule in aqueous solution, and a set of mutants of the channelrhodopsin (ChR) protein. The first test case shows that the PERI-CC2 method is in excellent agreement with the PE-CC2 method and in good agreement with the PE-CCSD method. The uracil test case indicates that the effects...... of hydrogen bonding on the TPA of a chromophore with the nearest environment is well-described with the PERI-CC2 method. Finally, the ChR calculation shows that the PERI-CC2 method is well-suited and efficient for calculations on proteins with medium-sized chromophores....

An accelerated hologram calculation using the wavefront recording plane method and wavelet transform

Science.gov (United States)

Arai, Daisuke; Shimobaba, Tomoyoshi; Nishitsuji, Takashi; Kakue, Takashi; Masuda, Nobuyuki; Ito, Tomoyoshi

2017-06-01

Fast hologram calculation methods are critical in real-time holography applications such as three-dimensional (3D) displays. We recently proposed a wavelet transform-based hologram calculation called WASABI. Even though WASABI can decrease the calculation time of a hologram from a point cloud, it increases the calculation time with increasing propagation distance. We also proposed a wavefront recoding plane (WRP) method. This is a two-step fast hologram calculation in which the first step calculates the superposition of light waves emitted from a point cloud in a virtual plane, and the second step performs a diffraction calculation from the virtual plane to the hologram plane. A drawback of the WRP method is in the first step when the point cloud has a large number of object points and/or a long distribution in the depth direction. In this paper, we propose a method combining WASABI and the WRP method in which the drawbacks of each can be complementarily solved. Using a consumer CPU, the proposed method succeeded in performing a hologram calculation with 2048 × 2048 pixels from a 3D object with one million points in approximately 0.4 s.

Biexponential analysis of diffusion-weighted imaging: comparison of three different calculation methods in transplanted kidneys.

Science.gov (United States)

Heusch, Philipp; Wittsack, Hans-Jörg; Pentang, Gael; Buchbender, Christian; Miese, Falk; Schek, Julia; Kröpil, Patric; Antoch, Gerald; Lanzman, Rotem S

2013-12-01

Biexponential analysis has been used increasingly to obtain contributions of both diffusion and microperfusion to the signal decay in diffusion-weighted imaging DWI of different parts of the body. To compare biexponential diffusion parameters of transplanted kidneys obtained with three different calculation methods. DWI was acquired in 15 renal allograft recipients (eight men, seven women; mean age, 52.4 ± 14.3 years) using a paracoronal EPI sequence with 16 b-values (b = 0-750 s/mm(2)) and six averages at 1.5T. No respiratory gating was used. Three different calculation methods were used for the calculation of biexponential diffusion parameters: Fp, ADCP, and ADCD were calculated without fixing any parameter a priori (calculation method 1); ADCP was fixed to 12.0 µm(2)/ms, whereas Fp and ADCD were calculated using the biexponential model (calculation method 2); multistep approach with monoexponential fitting of the high b-value portion (b ≥ 250 s/mm(2)) for determination of ADCD and assessment of the low b intercept for determination of Fp (calculation method 3). For quantitative analysis, ROI measurements were performed on the according parameter maps. Mean ADCD values of the renal cortex using calculation method 1 were significantly lower than using calculation methods 2 and 3 (P < 0.001). There was a significant correlation between calculation methods 1 and 2 (r = 0.69 (P < 0.005) and calculation methods 1 and 3 (r = 0.59; P < 0.05) as well as calculation methods 2 and 3 (r = 0.98; P < 0.001). Mean Fp values of the renal cortex were higher with calculation method 1 than with calculation methods 2 and 3 (P < 0.001). For Fp, only the correlation between calculation methods 2 and 3 was significant (r = 0.98; P < 0.001). Biexponential diffusion parameters differ significantly depending on the calculation methods used for their calculation.

A method for calculating effective lifetime risk of radiation-induced cancer from screening mammography

International Nuclear Information System (INIS)

Ali, R.M.; England, A.; McEntee, M.F.; Hogg, P.

2015-01-01

Purpose: To propose a method for evaluating the effective lifetime risk of radiation-induced cancer from screening mammography and to present initial data for the UK National Breast Screening Programme. Material and methods: The imaging was undertaken using a Hologic Selenia full field digital mammographic unit. The proposed method utilises an ATOM phantom containing thermoluminescent dosimeters and a perspex-polyethylene breast phantom to measure organ doses during a standard four view screening mammogram. Effective dose was calculated and effective risk was modelled for a range of client ages. The total lifetime effective risk was then calculated for the UK national screening programme. Calculation of effective risk includes the radiation dose to examined and contralateral breasts in addition to other body organs; this is an advantage over the mean glandular dose. Results: The contralateral breast, thyroid, thymus, brain, lung, salivary glands, and bone marrow all receive more than 1 μGy radiation dose during screening mammography. A major difference exists for total effective lifetime risk of radiation-induced cancer between clients with average and high breast cancer risk. Differences are attributed to the commencement age of screening and time interval between screens. Conclusion: This study proposes a method to evaluate effective lifetime risk of radiation-induced cancer from screening mammography in order to compare different mammography screening programmes. - Highlights: • We proposed a method for the calculation of radiation-induced cancer from screening mammography. • We measured the radiation absorbed dose of different organs during screening mammography. • There are major differences between mammography screening programme categories with regard to radiation effective risk.

Method for calculating required shielding in medical x-ray rooms

International Nuclear Information System (INIS)

Karppinen, J.

1997-10-01

The new annual radiation dose limits - 20 mSv (previously 50 mSv) for radiation workers and 1 mSv (previously 5 mSv) for other persons - implies that the adequacy of existing radiation shielding must be re-evaluated. In principle, one could assume that the thicknesses of old radiation shields should be increased by about one or two half-value layers in order to comply with the new dose limits. However, the assumptions made in the earlier shielding calculations are highly conservative; the required shielding was often determined by applying the maximum high-voltage of the x-ray tube for the whole workload. A more realistic calculation shows that increased shielding is typically not necessary if more practical x-ray tube voltages are used in the evaluation. We have developed a PC-based calculation method for calculating the x-ray shielding which is more realistic than the highly conservative method formerly used. The method may be used to evaluate an existing shield for compliance with new regulations. As examples of these calculations, typical x-ray rooms are considered. The lead and concrete thickness requirements as a function of x-ray tube voltage and workload are also given in tables. (author)

Methods for Calculating Empires in Quasicrystals

Directory of Open Access Journals (Sweden)

Fang Fang

2017-10-01

Full Text Available This paper reviews the empire problem for quasiperiodic tilings and the existing methods for generating the empires of the vertex configurations in quasicrystals, while introducing a new and more efficient method based on the cut-and-project technique. Using Penrose tiling as an example, this method finds the forced tiles with the restrictions in the high dimensional lattice (the mother lattice that can be cut-and-projected into the lower dimensional quasicrystal. We compare our method to the two existing methods, namely one method that uses the algorithm of the Fibonacci chain to force the Ammann bars in order to find the forced tiles of an empire and the method that follows the work of N.G. de Bruijn on constructing a Penrose tiling as the dual to a pentagrid. This new method is not only conceptually simple and clear, but it also allows us to calculate the empires of the vertex configurations in a defected quasicrystal by reversing the configuration of the quasicrystal to its higher dimensional lattice, where we then apply the restrictions. These advantages may provide a key guiding principle for phason dynamics and an important tool for self error-correction in quasicrystal growth.

Methods for calculating radiation attenuation in shields

Energy Technology Data Exchange (ETDEWEB)

Butler, J; Bueneman, D; Etemad, A; Lafore, P; Moncassoli, A M; Penkuhn, H; Shindo, M; Stoces, B

1964-10-01

In recent years the development of high-speed digital computers of large capacity has revolutionized the field of reactor shield design. For compact special-purpose reactor shields, Monte-Carlo codes in two- and three dimensional geometries are now available for the proper treatment of both the neutron and gamma- ray problems. Furthermore, techniques are being developed for the theoretical optimization of minimum-weight shield configurations for this type of reactor system. In the design of land-based power reactors, on the other hand, there is a strong incentive to reduce the capital cost of the plant, and economic considerations are also relevant to reactors designed for merchant ship propulsion. In this context simple methods are needed which are economic in their data input and computing time requirements and which, at the same time, are sufficiently accurate for design work. In general the computing time required for Monte-Carlo calculations in complex geometry is excessive for routine design calculations and the capacity of the present codes is inadequate for the proper treatment of large reactor shield systems in three dimensions. In these circumstances a wide range of simpler techniques are currently being employed for design calculations. The methods of calculation for neutrons in reactor shields fall naturally into four categories: Multigroup diffusion theory; Multigroup diffusion with removal sources; Transport codes; and Monte Carlo methods. In spite of the numerous Monte- Carlo techniques which are available for penetration and back scattering, serious problems are still encountered in practice with the scattering of gamma rays from walls of buildings which contain critical facilities and also concrete-lined discharge shafts containing irradiated fuel elements. The considerable volume of data in the unclassified literature on the solution of problems of this type in civil defence work appears not to have been evaluated for reactor shield design. In

Uneconomical top calculation method

International Nuclear Information System (INIS)

De Noord, M.; Vanm Sambeek, E.J.W.

2003-08-01

The methodology used to calculate the financial gap of renewable electricity sources and technologies is described. This methodology is used for calculating the production subsidy levels (MEP subsidies) for new renewable electricity projects in 2004 and 2005 in the Netherlands [nl

A Gaussian quadrature method for total energy analysis in electronic state calculations

Science.gov (United States)

Fukushima, Kimichika

This article reports studies by Fukushima and coworkers since 1980 concerning their highly accurate numerical integral method using Gaussian quadratures to evaluate the total energy in electronic state calculations. Gauss-Legendre and Gauss-Laguerre quadratures were used for integrals in the finite and infinite regions, respectively. Our previous article showed that, for diatomic molecules such as CO and FeO, elliptic coordinates efficiently achieved high numerical integral accuracy even with a numerical basis set including transition metal atomic orbitals. This article will generalize straightforward details for multiatomic systems with direct integrals in each decomposed elliptic coordinate determined from the nuclear positions of picked-up atom pairs. Sample calculations were performed for the molecules O3 and H2O. This article will also try to present, in another coordinate, a numerical integral by partially using the Becke's decomposition published in 1988, but without the Becke's fuzzy cell generated by the polynomials of internuclear distance between the pair atoms. Instead, simple nuclear weights comprising exponential functions around nuclei are used. The one-center integral is performed with a Gaussian quadrature pack in a spherical coordinate, included in the author's original program in around 1980. As for this decomposition into one-center integrals, sample calculations are carried out for Li2.

Fast calculation method for computer-generated cylindrical holograms.

Science.gov (United States)

Yamaguchi, Takeshi; Fujii, Tomohiko; Yoshikawa, Hiroshi

2008-07-01

Since a general flat hologram has a limited viewable area, we usually cannot see the other side of a reconstructed object. There are some holograms that can solve this problem. A cylindrical hologram is well known to be viewable in 360 deg. Most cylindrical holograms are optical holograms, but there are few reports of computer-generated cylindrical holograms. The lack of computer-generated cylindrical holograms is because the spatial resolution of output devices is not great enough; therefore, we have to make a large hologram or use a small object to fulfill the sampling theorem. In addition, in calculating the large fringe, the calculation amount increases in proportion to the hologram size. Therefore, we propose what we believe to be a new calculation method for fast calculation. Then, we print these fringes with our prototype fringe printer. As a result, we obtain a good reconstructed image from a computer-generated cylindrical hologram.

Methods for reactor physics calculations for control rods in fast reactors

International Nuclear Information System (INIS)

Grimstone, M.J.; Rowlands, J.L.

1988-12-01

The IAEA Specialists' Meeting on ''Methods for Reactor Physics Calculations for Control Rods in Fast Reactors'' was held in Winfrith, United Kingdom, on 6-8 December, 1988. The meeting was attended by 23 participants from nine countries. The purpose of the meeting was to review the current calculational methods and their accuracy as assessed by theoretical studies and comparisons with measurements, and then to identify the requirements for improved methods or additional studies and comparisons. The control rod properties or effects to be considered were their reactivity worths, their effect on the power distribution through the core, and the reaction rates and energy deposition both within and adjacent to the rods. The meeting was divided into five sessions, in the first of which each national delegation presented a brief overview of their programme of work on calculational methods for fast reactor control rods. In the next three sessions a total of seventeen papers were presented describing calculational methods and assessments of their accuracy. The final session was a discussion to draw conclusions regarding the current status of methods and the further developments and validation work required. A separate abstract was prepared for each of the 23 papers presented at the meeting. Refs, figs and tabs

A comparison of the calculation methods of the maze shielding dose

International Nuclear Information System (INIS)

Li Wenqian; Li Junli; Li Pengyu; Tao Yinghua

2009-01-01

This paper gives a theoretical calculating method for the dose rate of the maze of the low-energy accelerators or high-energy accelerators, based on the NCRP report Nos.49, 51 and 151. The multi-legged maze of the Miyun CT workshop of the NUCTECH Company Limited and the arc maze of the radiation laboratory of the Academy of Military Medical Sciences were calculated using this method. The calculating results were compared with the MCNP simulating results and the measured results. For the commonly estimation of the maze dose rate, as long as the parameters chosen properly, this method can give a conservative result, and save time from simulation. It's hoped that this work could offer a reference for the maze design and the dose estimation method in the aftertime. (authors)

Improvement of methods for calculation of sound insulation in buildings

OpenAIRE

Mašović, Draško B.

2015-01-01

The main object of this work are the methods for calculation of sound insulation based on the classical model of sound propagation in buildings and single-number rating of sound insulation. The aim of the work is inspection of the possibilities for improvement of standard methods for quantification and calculation of sound insulation, in order to achieve higher accuracy of the obtained numerical values and their correlation with subjective impression of the acoustic comfort in buildings. Proc...

Application of the variational method for calculation of neutron spectra and group constants - Master thesis

International Nuclear Information System (INIS)

Milosevic, M.

1979-01-01

One-dimensional variational method for cylindrical configuration was applied for calculating group constants, together with effects of elastic slowing down, anisotropic elastic scattering, inelastic scattering, heterogeneous resonance absorption with the aim to include the presence of a number of different isotopes and effects of neutron leakage from the reactor core. Neutron flux shape P 3 and adjoint function are proposed in order to enable calculation of smaller size reactors and inclusion of heterogeneity effects by cell calculations. Microscopic multigroup constants were prepared based on the UKNDL data library. Analytical-numerical approach was applied for solving the equations of the P 3 approximation to obtain neutron flux moments and adjoint functions

Fourier Bessel transform method for efficiently calculating the magnetic field of solenoids

International Nuclear Information System (INIS)

Nachamkin, J.; Maggiore, C.J.

1980-01-01

A numerical procedure for calculating the magnetic field of a selenoid is derived. Based on the properties of Bessel functions, the procedure is shown to be convergent everywhere, including within the windings of the solenoid. The most critical part of the procedure is detailed in the main text. A simple method is used to ensure numerical significance while allowing economical computational times. In the appendix the procedure is generalized to universal convergence by appropriate partitioning of the solenoid windings

Temperature Calculation of Annular Fuel Pellet by Finite Difference Method

Energy Technology Data Exchange (ETDEWEB)

Yang, Yong Sik; Bang, Je Geon; Kim, Dae Ho; Kim, Sun Ki; Lim, Ik Sung; Song, Kun Woo [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)

2009-10-15

KAERI has started an innovative fuel development project for applying dual-cooled annular fuel to existing PWR reactor. In fuel design, fuel temperature is the most important factor which can affect nuclear fuel integrity and safety. Many models and methodologies, which can calculate temperature distribution in a fuel pellet have been proposed. However, due to the geometrical characteristics and cooling condition differences between existing solid type fuel and dual-cooled annular fuel, current fuel temperature calculation models can not be applied directly. Therefore, the new heat conduction model of fuel pellet was established. In general, fuel pellet temperature is calculated by FDM(Finite Difference Method) or FEM(Finite Element Method), because, temperature dependency of fuel thermal conductivity and spatial dependency heat generation in the pellet due to the self-shielding should be considered. In our study, FDM is adopted due to high exactness and short calculation time.

A simple method for calculation of Glauber's amplitude

International Nuclear Information System (INIS)

Omboo, Z.

1983-01-01

A method of calculating the terms of Glauber series expansions for elastic scattering of composed systems are presented. The inclusion of general scattering diagram simplifies essentially the calculation procedure. In this case the complicated combinatorical problem of reduction of similar terms in Glauber series is solved easily and determinant corresponding to various terms of the series decreases at least by a factor of two, if numbers of constituents of scattered systems are equal. If these numbers are not equal, the determinant order is equal to the smallest one

A coupling method for a cardiovascular simulation model which includes the Kalman filter.

Science.gov (United States)

Hasegawa, Yuki; Shimayoshi, Takao; Amano, Akira; Matsuda, Tetsuya

2012-01-01

Multi-scale models of the cardiovascular system provide new insight that was unavailable with in vivo and in vitro experiments. For the cardiovascular system, multi-scale simulations provide a valuable perspective in analyzing the interaction of three phenomenons occurring at different spatial scales: circulatory hemodynamics, ventricular structural dynamics, and myocardial excitation-contraction. In order to simulate these interactions, multiscale cardiovascular simulation systems couple models that simulate different phenomena. However, coupling methods require a significant amount of calculation, since a system of non-linear equations must be solved for each timestep. Therefore, we proposed a coupling method which decreases the amount of calculation by using the Kalman filter. In our method, the Kalman filter calculates approximations for the solution to the system of non-linear equations at each timestep. The approximations are then used as initial values for solving the system of non-linear equations. The proposed method decreases the number of iterations required by 94.0% compared to the conventional strong coupling method. When compared with a smoothing spline predictor, the proposed method required 49.4% fewer iterations.

Novel and Efficient Methods for Calculating Pressure in Polymer Lattice Models

Science.gov (United States)

Zhang, Pengfei; Wang, Qiang

2014-03-01

Pressure calculation in polymer lattice models is an important but nontrivial subject. The three existing methods - thermodynamic integration, repulsive wall, and sedimentation equilibrium methods - all have their limitations and cannot be used to accurately calculate the pressure at all polymer volume fractions φ. Here we propose two novel methods. In the first method, we combine Monte Carlo simulation in an expanded grand-canonical ensemble with the Wang-Landau - Optimized Ensemble (WL-OE) simulation to calculate the pressure as a function of polymer volume fraction, which is very efficient at low to intermediate φ and exhibits negligible finite-size effects. In the second method, we introduce a repulsive plane with bridging bonds, which is similar to the repulsive wall method but eliminates its confinement effects, and estimate the two-dimensional density of states (in terms of the number of bridging bonds and the contact number) using the 1/ t version of Wang-Landau algorithm. This works well at all φ, especially at high φ where all the methods involving chain insertion trial moves fail.

KaKs_Calculator 2.0: A Toolkit Incorporating Gamma-Series Methods and Sliding Window Strategies

KAUST Repository

Wang, Dapeng

2010-05-05

We present an integrated stand-alone software package named KaKs_Calculator 2.0 as an updated version. It incorporates 17 methods for the calculation of nonsynonymous and synonymous substitution rates; among them, we added our modified versions of several widely used methods as the gamma series including γ-NG, γ-LWL, γ-MLWL, γ-LPB, γ-MLPB, γ-YN and γ-MYN, which have been demonstrated to perform better under certain conditions than their original forms and are not implemented in the previous version. The package is readily used for the identification of positively selected sites based on a sliding window across the sequences of interests in 5\\' to 3\\' direction of protein-coding sequences, and have improved the overall performance on sequence analysis for evolution studies. A toolbox, including C++ and Java source code and executable files on both Windows and Linux platforms together with a user instruction, is downloadable from the website for academic purpose at https://sourceforge.net/projects/kakscalculator2/.

KaKs_Calculator 2.0: A Toolkit Incorporating Gamma-Series Methods and Sliding Window Strategies

KAUST Repository

Wang, Dapeng; Zhang, Yubin; Zhang, Zhang; Zhu, Jiang; Yu, Jun

2010-01-01

We present an integrated stand-alone software package named KaKs_Calculator 2.0 as an updated version. It incorporates 17 methods for the calculation of nonsynonymous and synonymous substitution rates; among them, we added our modified versions of several widely used methods as the gamma series including γ-NG, γ-LWL, γ-MLWL, γ-LPB, γ-MLPB, γ-YN and γ-MYN, which have been demonstrated to perform better under certain conditions than their original forms and are not implemented in the previous version. The package is readily used for the identification of positively selected sites based on a sliding window across the sequences of interests in 5' to 3' direction of protein-coding sequences, and have improved the overall performance on sequence analysis for evolution studies. A toolbox, including C++ and Java source code and executable files on both Windows and Linux platforms together with a user instruction, is downloadable from the website for academic purpose at https://sourceforge.net/projects/kakscalculator2/.

Method for consequence calculations for severe accidents

International Nuclear Information System (INIS)

Nielsen, F.; Thykier-Nielsen, S.; Walmod-Larsen, O.

1986-08-01

This report was commissioned by the Swedish State Power Board, who wanted a method for calculation of radiation doses in the surroundings of nuclear power plants caused by severe accidents. The PLUCON4 code were used for the calculations. A TC-SV-accident at Ringhals 1 wer chosen as example. A transient without shutdown leads to core meltdown through the reactor vessel. The pressure peak at the moment of vessel failure opens a safety valve in the dry well. Meteorolgical data for two years from the Ringhals meteorological tower were analysed to find representative weather situations. As typical weather were chosen Pasquill D with wind speed 8 m/s, and as extreme weather were chosen Pasquill F with wind speed 4.8 m/s. (author)

Internal quality control of RIA with Tonks error calculation method

International Nuclear Information System (INIS)

Chen Xiaodong

1996-01-01

According to the methodology feature of RIA, an internal quality control chart with Tonks error calculation method which is suitable for RIA is designed. The quality control chart defines the value of the allowance error with normal reference range. The method has the simplicity of its performance and directly perceived through the senses. Taking the example of determining T 3 and T 4 , the calculation of allowance error, drawing of quality control chart and the analysis of result are introduced

A parallel orbital-updating based plane-wave basis method for electronic structure calculations

International Nuclear Information System (INIS)

Pan, Yan; Dai, Xiaoying; Gironcoli, Stefano de; Gong, Xin-Gao; Rignanese, Gian-Marco; Zhou, Aihui

2017-01-01

Highlights: • Propose three parallel orbital-updating based plane-wave basis methods for electronic structure calculations. • These new methods can avoid the generating of large scale eigenvalue problems and then reduce the computational cost. • These new methods allow for two-level parallelization which is particularly interesting for large scale parallelization. • Numerical experiments show that these new methods are reliable and efficient for large scale calculations on modern supercomputers. - Abstract: Motivated by the recently proposed parallel orbital-updating approach in real space method , we propose a parallel orbital-updating based plane-wave basis method for electronic structure calculations, for solving the corresponding eigenvalue problems. In addition, we propose two new modified parallel orbital-updating methods. Compared to the traditional plane-wave methods, our methods allow for two-level parallelization, which is particularly interesting for large scale parallelization. Numerical experiments show that these new methods are more reliable and efficient for large scale calculations on modern supercomputers.

Calculation of mixed mode stress intensity factors using an alternating method

International Nuclear Information System (INIS)

Sakai, Takayuki

1999-01-01

In this study, mixed mode stress intensity factors (K I and K II ) of a square plate with a notch were calculated using a finite element alternating method. The obtained results were compared with the ones by a finite element method, and it was shown that the finite element alternating method can accurately estimate mixed mode stress intensity factors. Then, using this finite element alternating method, mixed mode stress intensity factors were calculated as changing the size and position of the notch, and its simplified equations were proposed. (author)

Pellet by pellet neutron flux calculations coupled with nodal expansion method

International Nuclear Information System (INIS)

Aldo, Dall'Osso

2003-01-01

We present a technique whose aim is to replace 2-dimensional pin by pin de-homogenization, currently done in core reactor calculations with the nodal expansion method (NEM), by a 3-dimensional finite difference diffusion calculation. This fine calculation is performed as a zoom in each node taking as boundary conditions the results of the NEM calculations. The size of fine mesh is of the order of a fuel pellet. The coupling between fine and NEM calculations is realised by an albedo like boundary condition. Some examples are presented showing fine neutron flux shape near control rods or assembly grids. Other fine flux behaviour as the thermal flux rise in the fuel near the reflector is emphasised. In general the results show the interest of the method in conditions where the separability of radial and axial directions is not granted. (author)

Numerical calculation of acoustic radiation from band-vibrating structures via FEM/FAQP method

Directory of Open Access Journals (Sweden)

GAO Honglin

2017-08-01

Full Text Available The Finite Element Method (FEM combined with the Frequency Averaged Quadratic Pressure method (FAQP are used to calculate the acoustic radiation of structures excited in the frequency band. The surface particle velocity of stiffened cylindrical shells under frequency band excitation is calculated using finite element software, the normal vibration velocity is converted from the surface particle velocity to calculate the average energy source (frequency averaged across intensity, frequency averaged across pressure and frequency averaged across velocity, and the FAQP method is used to calculate the average sound pressure level within the bandwidth. The average sound pressure levels are then compared with the bandwidth using finite element and boundary element software, and the results show that FEM combined with FAQP is more suitable for high frequencies and can be used to calculate the average sound pressure level in the 1/3 octave band with good stability, presenting an alternative to applying frequency-by-frequency calculation and the average frequency process. The FEM/FAQP method can be used as a prediction method for calculating acoustic radiation while taking the randomness of vibration at medium and high frequencies into consideration.

Thick-Restart Lanczos Method for Electronic Structure Calculations

International Nuclear Information System (INIS)

Simon, Horst D.; Wang, L.-W.; Wu, Kesheng

1999-01-01

This paper describes two recent innovations related to the classic Lanczos method for eigenvalue problems, namely the thick-restart technique and dynamic restarting schemes. Combining these two new techniques we are able to implement an efficient eigenvalue problem solver. This paper will demonstrate its effectiveness on one particular class of problems for which this method is well suited: linear eigenvalue problems generated from non-self-consistent electronic structure calculations

Faster method for the calculation of the chattering signal at the ct-detector by monte-carlo-method

International Nuclear Information System (INIS)

Schmidt, B.; Kalender, W.A.

2003-01-01

Multislice spiral CT scanners allow to acquire multiple slices simultaneously. With increasing numbers of slices, not only the total extent of slice collimation increases, but also the contribution of scatter radiation to the detector signal. A fast method for calculating the scatter signal would offer the possibility to correct the measured detector signal. Monte Carlo methods allow to simulate the paths of photons through a 3D volume, both in a patient- and scanner-specific fashion. If a scatter photon leaves the volume, its path can be followed and its interaction with an element of the detector be checked. This conventional way of calculating the scatter signal is time-consuming. In order to reduce the calculation time, a more efficient method was developed (Method of Weights). Every time an interaction occurs inside of the 3D volume, the probability of a detector hit due to photon scattering is calculated for each detector channel. The respective value is added to the scatter signal per detector with the corresponding weight. Simulated values of scatter-to-primary-signal ratios were confirmed by data available in the literature. Both the conventional and fast methods for the calculation of scatter signals yielded identical values within the range of statistical accuracy. Assuming the same computing time, the standard deviation for the conventional method was 5 times higher than for the fast one. The presented method allows to significantly reduce the computation time. It may therefore provide a basis for ''real time'' methods to correct for the scatter signal, especially in case of increasing numbers of slices. (orig.) [de

LEGO-Method--New Strategy for Chemistry Calculation

Science.gov (United States)

Molnar, Jozsef; Molnar-Hamvas, Livia

2011-01-01

The presented strategy of chemistry calculation is based on mole-concept, but it uses only one fundamental relationship of the amounts of substance as a basic panel. The name of LEGO-method comes from the famous toy of LEGO[R] because solving equations by grouping formulas is similar to that. The relations of mole and the molar amounts, as small…

The analytic method for calculating the control rod worth

International Nuclear Information System (INIS)

Kim, Han Gon; Lee, Byeong Ho; Chang, Soon Heung

1989-01-01

We calculated the control rod worth in this paper. To avoid complexity, we did not consider burnable poisons and soluble boron. The system was localized within one assembly. The control rod was treated as not an absorber but an another boundary. Thus all of the group constants were unchanged before and after control rod insertion. And we discussed the method for calculation of the reactivity of the whole core

Emergy Algebra: Improving Matrix Methods for Calculating Tranformities

Science.gov (United States)

Transformity is one of the core concepts in Energy Systems Theory and it is fundamental to the calculation of emergy. Accurate evaluation of transformities and other emergy per unit values is essential for the broad acceptance, application and further development of emergy method...

THE PARTICULARITIES OF THE COST CALCULATION METHOD ON COMMANDS IN FURNITURE INDUSTRY

Directory of Open Access Journals (Sweden)

Felicia Sabou

2014-10-01

Full Text Available The paper present the importance of the method on commands in cost calculation and the particularities of the cost calculation method on commands in the furniture industry. This paper presents a hypotetical study on the method on commands, considering the observations made during 2013-2014, on how it is organized and managed accounts management using method on commands.By presenting this hypothetical model about the accounting in management accounting using the method on commands, the paper contributes to the correct application of this method in practice, specifically in management accounting in companies from the furniture industry. In my opinion the method on commands is an appropriate method for achieving management accounting for companies that have as main activity the production of furniture. When applying the method on commands in cost calculation and in management accounting, the companies must to consider the particularities of the cost calculation, in the furniture industry, like: technical and economic factors from this sector, the technical details of each command, the codification of the commands, planning materials and labor costs for each command, monitoring and recording production costs, registration of the direct costs, distribution of the indirect costs on commands, registration of the indirect costs and registration in management accounting.

A Novel TRM Calculation Method by Probabilistic Concept

Science.gov (United States)

Audomvongseree, Kulyos; Yokoyama, Akihiko; Verma, Suresh Chand; Nakachi, Yoshiki

In a new competitive environment, it becomes possible for the third party to access a transmission facility. From this structure, to efficiently manage the utilization of the transmission network, a new definition about Available Transfer Capability (ATC) has been proposed. According to the North American ElectricReliability Council (NERC)’s definition, ATC depends on several parameters, i. e. Total Transfer Capability (TTC), Transmission Reliability Margin (TRM), and Capacity Benefit Margin (CBM). This paper is focused on the calculation of TRM which is one of the security margin reserved for any uncertainty of system conditions. The TRM calculation by probabilistic method is proposed in this paper. Based on the modeling of load forecast error and error in transmission line limitation, various cases of transmission transfer capability and its related probabilistic nature can be calculated. By consideration of the proposed concept of risk analysis, the appropriate required amount of TRM can be obtained. The objective of this research is to provide realistic information on the actual ability of the network which may be an alternative choice for system operators to make an appropriate decision in the competitive market. The advantages of the proposed method are illustrated by application to the IEEJ-WEST10 model system.

Simple method to calculate percolation, Ising and Potts clusters

International Nuclear Information System (INIS)

Tsallis, C.

1981-01-01

A procedure ('break-collapse method') is introduced which considerably simplifies the calculation of two - or multirooted clusters like those commonly appearing in real space renormalization group (RG) treatments of bond-percolation, and pure and random Ising and Potts problems. The method is illustrated through two applications for the q-state Potts ferromagnet. The first of them concerns a RG calculation of the critical exponent ν for the isotropic square lattice: numerical consistence is obtained (particularly for q→0) with den Nijs conjecture. The second application is a compact reformulation of the standard star-triangle and duality transformations which provide the exact critical temperature for the anisotropic triangular and honeycomb lattices. (Author) [pt

Efficient Calculation of Near Fields in the FDTD Method

DEFF Research Database (Denmark)

Franek, Ondrej

2011-01-01

When calculating frequency-domain near fields by the FDTD method, almost 50 % reduction in memory and CPU operations can be achieved if only E-fields are stored during the main time-stepping loop and H-fields computed later. An improved method of obtaining the H-fields from Faraday's Law is prese...

Hamming method for solving the delayed neutron precursor concentration for reactivity calculation

International Nuclear Information System (INIS)

Díaz, Daniel Suescún; Ospina, Juan Felipe Flórez; Sarasty, Jesús Andrés Rodríguez

2012-01-01

Highlights: ► We present a new formulation to calculate the reactivity using the Hamming method. ► This method shows better accuracy than existing methods for reactivity calculation. ► The reactivity is calculated without limitation of the nuclear power form. ► The method can be implemented in reactivity meters with time step of up to 0.1 s. - Abstract: We propose a new method for numerically solving the inverse point kinetic equation for a nuclear reactor using the Hamming method, without requiring the nuclear power history and without using the Laplace transform. This new method converges with accuracy of order h 5 , where h is the step in the computation time. The procedure is validated for different forms of the nuclear power and with different time steps. The results indicate that this method has a better accuracy and lower computational effort compared with other conventional methods that use the nuclear power history.

Numerical calculation of electromagnetic properties including chirality parameters for uniaxial bianisotropic media

International Nuclear Information System (INIS)

Amirkhizi, Alireza V; Nemat-Nasser, Sia

2008-01-01

Through the use of conductive straight wires or coils the electromagnetic properties of a composite material can be modified. The asymmetric geometry of the coils creates an overall chiral response. The polarization vectors rotate as an electromagnetic wave travels through such a medium. To calculate the chirality of a medium prior to its manufacturing, we developed a method to extract all four electromagnetic material parameter tensors for a general uniaxial bianisotropic composite based on the numerical simulation of the electromagnetic fields. Our method uses appropriate line and surface field averages in a single unit cell of the periodic structure of the composite material. These overall field quantities have physical meaning only when the microscopic variation of the electromagnetic fields in the scale of the unit cell is not important, that is when the wavelength of interest is significantly larger than the maximum linear dimension of the unit cell. The overall constitutive relations of the periodic structure can then be obtained from the relations among the average quantities

Development of approximate shielding calculation method for high energy cosmic radiation on LEO satellites

International Nuclear Information System (INIS)

Sin, M. W.; Kim, M. H.

2002-01-01

To calculate total dose effect on semi-conductor devices in satellite for a period of space mission effectively, two approximate calculation models for a comic radiation shielding were proposed. They are a sectoring method and a chord-length distribution method. When an approximate method was applied in this study, complex structure of satellite was described into multiple 1-dimensional slabs, structural materials were converted to reference material(aluminum), and the pre-calculated dose-depth conversion function was introduced to simplify the calculation process. Verification calculation was performed for orbit location and structure geometry of KITSAT-1 and compared with detailed 3-dimensional calculation results and experimental values. The calculation results from approximate method were estimated conservatively with acceptable error. However, results for satellite mission simulation were underestimated in total dose rate compared with experimental values

Development of approximate shielding calculation method for high energy cosmic radiation on LEO satellites

Energy Technology Data Exchange (ETDEWEB)

Sin, M. W.; Kim, M. H. [Kyunghee Univ., Yongin (Korea, Republic of)

2002-10-01

To calculate total dose effect on semi-conductor devices in satellite for a period of space mission effectively, two approximate calculation models for a comic radiation shielding were proposed. They are a sectoring method and a chord-length distribution method. When an approximate method was applied in this study, complex structure of satellite was described into multiple 1-dimensional slabs, structural materials were converted to reference material(aluminum), and the pre-calculated dose-depth conversion function was introduced to simplify the calculation process. Verification calculation was performed for orbit location and structure geometry of KITSAT-1 and compared with detailed 3-dimensional calculation results and experimental values. The calculation results from approximate method were estimated conservatively with acceptable error. However, results for satellite mission simulation were underestimated in total dose rate compared with experimental values.

Cluster-cell calculation using the method of generalized homogenization

International Nuclear Information System (INIS)

Laletin, N.I.; Boyarinov, V.F.

1988-01-01

The generalized-homogenization method (GHM), used for solving the neutron transfer equation, was applied to calculating the neutron distribution in the cluster cell with a series of cylindrical cells with cylindrically coaxial zones. Single-group calculations of the technological channel of the cell of an RBMK reactor were performed using GHM. The technological channel was understood to be the reactor channel, comprised of the zirconium rod, the water or steam-water mixture, the uranium dioxide fuel element, and the zirconium tube, together with the adjacent graphite layer. Calculations were performed for channels with no internal sources and with unit incoming current at the external boundary as well as for channels with internal sources and zero current at the external boundary. The PRAKTINETs program was used to calculate the symmetric neutron distributions in the microcell and in channels with homogenized annular zones. The ORAR-TsM program was used to calculate the antisymmetric distribution in the microcell. The accuracy of the calculations were compared for the two channel versions

Using the SAPFIR_95&RC_WWER combined method for calculation of pin-wise power distribution in the WWER reactor during power manoeuvres

International Nuclear Information System (INIS)

Artemov, V.G.; Ivanov, A.S.; Kuznetsov, A.N.; Shemaev, Yu.P.

2015-01-01

The method of calculating the pin-wise power distribution by a combined method has been tested using the example of calculating the maneuvering regime. The combined method makes it possible to clarify the behavior of the linear fuel element load in comparison with the results of the superposition method in fuel assemblies with the regulator of the control and protection system in the boundaries of the control groups. The use of prepared on the basis of the results of macro calculation of the pin-wise characteristics, including 135Xe and 135I concentrations, allows to perform an adequate fine-grid calculation of the pin-wise power distribution without expensive calculation of feedbacks on a fine grid, which limits the application of a direct method for solving practical problems [ru

Realistic methods for calculating the releases and consequences of a large LOCA

International Nuclear Information System (INIS)

Stephenson, W.; Dutton, L.M.C.; Handy, B.J.; Smedley, C.

1992-01-01

This report describes a calculational route to predict realistic radiological consequences for a successfully terminated large-loss-of-coolant accident (LOCA) at a pressurized-water reactor (PWR). All steps in the calculational route are considered. For each one, a brief comment is made on the significant differences between the methods of calculation that were identified in the benchmark studies and recommendations are made for the methods and data for carrying out realistic calculations. These are based on the best supportable methods and data and the technical basis for each recommendation is given. Where the lack of well-validated methods or data means that the most realistic method that can be justified is considered to be very conservative, the need for further research is identified. The behaviour of inorganic iodine and the removal of aerosols from the atmosphere of the reactor building are identified as areas of particular importance. Where the retention of radioactivity is sensitive to design features, these are identified and, for the most importance features, the impact of different designs on the release of activity is indicated. The predictions of the proposed model are calculated for each stage and compared with the releases of activity predicted by the licensing methods that were used in the earlier benchmark studies. The conservative nature of the latter is confirmed. Methods and data are also presented for calculating the resulting doses to members of the public of the National Radiological Protection Boards as a result of work carried out by several national bodies in the UK. Other, equally acceptable, models are used in other countries of the Community and some examples are given

Comparison of calculational methods for liquid metal reactor shields

International Nuclear Information System (INIS)

Carter, L.L.; Moore, F.S.; Morford, R.J.; Mann, F.M.

1985-09-01

A one-dimensional comparison is made between Monte Carlo (MCNP), discrete ordinances (ANISN), and diffusion theory (MlDX) calculations of neutron flux and radiation damage from the core of the Fast Flux Test Facility (FFTF) out to the reactor vessel. Diffusion theory was found to be reasonably accurate for the calculation of both total flux and radiation damage. However, for large distances from the core, the calculated flux at very high energies is low by an order of magnitude or more when the diffusion theory is used. Particular emphasis was placed in this study on the generation of multitable cross sections for use in discrete ordinates codes that are self-shielded, consistent with the self-shielding employed in the generation of cross sections for use with diffusion theory. The Monte Carlo calculation, with a pointwise representation of the cross sections, was used as the benchmark for determining the limitations of the other two calculational methods. 12 refs., 33 figs

Assessment of chemical exposures: calculation methods for environmental professionals

National Research Council Canada - National Science Library

Daugherty, Jack E

1997-01-01

... on by scientists, businessmen, and policymakers. Assessment of Chemical Exposures: Calculation Methods for Environmental Professionals addresses the expanding scope of exposure assessments in both the workplace and environment...

The calculations of small molecular conformation energy differences by density functional method

Science.gov (United States)

Topol, I. A.; Burt, S. K.

1993-03-01

The differences in the conformational energies for the gauche (G) and trans(T) conformers of 1,2-difluoroethane and for myo-and scyllo-conformer of inositol have been calculated by local density functional method (LDF approximation) with geometry optimization using different sets of calculation parameters. It is shown that in the contrast to Hartree—Fock methods, density functional calculations reproduce the correct sign and value of the gauche effect for 1,2-difluoroethane and energy difference for both conformers of inositol. The results of normal vibrational analysis for1,2-difluoroethane showed that harmonic frequencies calculated in LDF approximation agree with experimental data with the accuracy typical for scaled large basis set Hartree—Fock calculations.

A combination between the differential and the perturbation theory methods for calculating sensitivity coefficients

International Nuclear Information System (INIS)

Borges, Antonio Andrade

1998-01-01

A new method for the calculation of sensitivity coefficients is developed. The new method is a combination of two methodologies used for calculating theses coefficients, which are the differential and the generalized perturbation theory methods. The method utilizes as integral parameter the average flux in an arbitrary region of the system. Thus, the sensitivity coefficient contains only the component corresponding to the neutron flux. To obtain the new sensitivity coefficient, the derivatives of the integral parameter, Φ, with respect to σ are calculated using the perturbation method and the functional derivatives of this generic integral parameter with respect to σ and Φ are calculated using the differential method. (author)

Simple Calculation Programs for Biology Methods in Molecular ...

Indian Academy of Sciences (India)

First page Back Continue Last page Overview Graphics. Simple Calculation Programs for Biology Methods in Molecular Biology. GMAP: A program for mapping potential restriction sites. RE sites in ambiguous and non-ambiguous DNA sequence; Minimum number of silent mutations required for introducing a RE sites; Set ...

METHOD OF CALCULATING THE OPTIMAL HEAT EMISSION GEOTHERMAL WELLS

Directory of Open Access Journals (Sweden)

A. I. Akaev

2015-01-01

Full Text Available This paper presents a simplified method of calculating the optimal regimes of the fountain and the pumping exploitation of geothermal wells, reducing scaling and corrosion during operation. Comparative characteristics to quantify the heat of formation for these methods of operation under the same pressure at the wellhead. The problem is solved graphic-analytical method based on a balance of pressure in the well with the heat pump. 

On Calculation Methods and Results for Straight Cylindrical Roller Bearing Deflection, Stiffness, and Stress

Science.gov (United States)

Krantz, Timothy L.

2011-01-01

The purpose of this study was to assess some calculation methods for quantifying the relationships of bearing geometry, material properties, load, deflection, stiffness, and stress. The scope of the work was limited to two-dimensional modeling of straight cylindrical roller bearings. Preparations for studies of dynamic response of bearings with damaged surfaces motivated this work. Studies were selected to exercise and build confidence in the numerical tools. Three calculation methods were used in this work. Two of the methods were numerical solutions of the Hertz contact approach. The third method used was a combined finite element surface integral method. Example calculations were done for a single roller loaded between an inner and outer raceway for code verification. Next, a bearing with 13 rollers and all-steel construction was used as an example to do additional code verification, including an assessment of the leading order of accuracy of the finite element and surface integral method. Results from that study show that the method is at least first-order accurate. Those results also show that the contact grid refinement has a more significant influence on precision as compared to the finite element grid refinement. To explore the influence of material properties, the 13-roller bearing was modeled as made from Nitinol 60, a material with very different properties from steel and showing some potential for bearing applications. The codes were exercised to compare contact areas and stress levels for steel and Nitinol 60 bearings operating at equivalent power density. As a step toward modeling the dynamic response of bearings having surface damage, static analyses were completed to simulate a bearing with a spall or similar damage.

Bending Moment Calculations for Piles Based on the Finite Element Method

Directory of Open Access Journals (Sweden)

Yu-xin Jie

2013-01-01

Full Text Available Using the finite element analysis program ABAQUS, a series of calculations on a cantilever beam, pile, and sheet pile wall were made to investigate the bending moment computational methods. The analyses demonstrated that the shear locking is not significant for the passive pile embedded in soil. Therefore, higher-order elements are not always necessary in the computation. The number of grids across the pile section is important for bending moment calculated with stress and less significant for that calculated with displacement. Although computing bending moment with displacement requires fewer grid numbers across the pile section, it sometimes results in variation of the results. For displacement calculation, a pile row can be suitably represented by an equivalent sheet pile wall, whereas the resulting bending moments may be different. Calculated results of bending moment may differ greatly with different grid partitions and computational methods. Therefore, a comparison of results is necessary when performing the analysis.

Calculations of hydrogen detonations in nuclear containments by the random choice method

International Nuclear Information System (INIS)

Delichatsios, M.A.; Genadry, M.B.

1983-01-01

Computer codes were developed for the prediction of pressure histories at different points of a nuclear containment wall due to postulated internal hydrogen detonations. These pressure histories are required to assess the structural response of a nuclear containment to hydrogen detonations. The compressible flow equations including detonation, which was treated as a sharp fluid discontinuity, were solved by the random choice method which reproduces maximum pressures and discontinuities sharply. The computer codes were validated by calculating pressure profiles and maximum wall pressures for plane and spherical geometries and comparing the results with exact analytic solutions. The two-dimensional axisymmetric program was used to calculate wall pressure histories in an actual nuclear containment. The numerical results for wall pressures are presented in a dimensionless form, which allows their use for different combinations of hydrogen concentration, and initial conditions. (orig.)

Calculations of Neutron Flux Distributions by Means of Integral Transport Methods

Energy Technology Data Exchange (ETDEWEB)

Carlvik, I

1967-05-15

Flux distributions have been calculated mainly in one energy group, for a number of systems representing geometries interesting for reactor calculations. Integral transport methods of two kinds were utilised, collision probabilities (CP) and the discrete method (DIT). The geometries considered comprise the three one-dimensional geometries, planes, sphericals and annular, and further a square cell with a circular fuel rod and a rod cluster cell with a circular outer boundary. For the annular cells both methods (CP and DIT) were used and the results were compared. The purpose of the work is twofold, firstly to demonstrate the versatility and efficacy of integral transport methods and secondly to serve as a guide for anybody who wants to use the methods.

The method of calculation of pipelines laid on supports

OpenAIRE

Benin D.M.

2017-01-01

this article focuses on the issue of laying pipelines on supports and the method of calculation of vertical and horizontal loads acting on the support. As pipelines can be water piping systems, heat networks, oil and mazout lines, condensate lines, steam lines, etc. this article describes the calculations of supports for pipelines laid above ground, in crowded channels, premises, on racks, in impassable channels, hanging supports, etc. The paper explores recommendations for placement of the s...

Nuclear data and multigroup methods in fast reactor calculations

International Nuclear Information System (INIS)

Gur, Y.

1975-03-01

The work deals with fast reactor multigroup calculations, and the efficient treatment of basic nuclear data, which serves as raw material for the calculations. Its purpose is twofold: to build a computer code system that handles a large, detailed library of basic neutron cross section data, (such as ENDF/B-III) and yields a compact set of multigroup cross sections for reactor calculations; to use the code system for comparative analysis of different libraries, in order to discover basic uncertainties that still exist in the measurement of neutron cross sections, and to determine their influence upon uncertainties in nuclear calculations. A program named NANICK which was written in two versions is presented. The first handles the American basic data library, ENDF/B-III, while the second handles the German basic data library, KEDAK. The mathematical algorithm is identical in both versions, and only the file management is different. This program calculates infinitely diluted multigroup cross sections and scattering matrices. It is complemented by the program NASIF that calculates shielding factors from resonance parameters. Different versions of NASIF were written to handle ENDF/B-III or KEDAK. New methods for evaluating in reactor calculations the long term behavior of the neutron flux as well as its fine structure are described and an efficient calculation of the shielding factors from resonance parameters is offered. (B.G.)

A Method of Calculating Motion Error in a Linear Motion Bearing Stage

Directory of Open Access Journals (Sweden)

Gyungho Khim

2015-01-01

Full Text Available We report a method of calculating the motion error of a linear motion bearing stage. The transfer function method, which exploits reaction forces of individual bearings, is effective for estimating motion errors; however, it requires the rail-form errors. This is not suitable for a linear motion bearing stage because obtaining the rail-form errors is not straightforward. In the method described here, we use the straightness errors of a bearing block to calculate the reaction forces on the bearing block. The reaction forces were compared with those of the transfer function method. Parallelism errors between two rails were considered, and the motion errors of the linear motion bearing stage were measured and compared with the results of the calculations, revealing good agreement.

A Method of Calculating Motion Error in a Linear Motion Bearing Stage

Science.gov (United States)

Khim, Gyungho; Park, Chun Hong; Oh, Jeong Seok

2015-01-01

We report a method of calculating the motion error of a linear motion bearing stage. The transfer function method, which exploits reaction forces of individual bearings, is effective for estimating motion errors; however, it requires the rail-form errors. This is not suitable for a linear motion bearing stage because obtaining the rail-form errors is not straightforward. In the method described here, we use the straightness errors of a bearing block to calculate the reaction forces on the bearing block. The reaction forces were compared with those of the transfer function method. Parallelism errors between two rails were considered, and the motion errors of the linear motion bearing stage were measured and compared with the results of the calculations, revealing good agreement. PMID:25705715

Cell homogenization methods for pin-by-pin core calculations tested in slab geometry

International Nuclear Information System (INIS)

Yamamoto, Akio; Kitamura, Yasunori; Yamane, Yoshihiro

2004-01-01

In this paper, performances of spatial homogenization methods for fuel or non-fuel cells are compared in slab geometry in order to facilitate pin-by-pin core calculations. Since the spatial homogenization methods were mainly developed for fuel assemblies, systematic study of their performance for the cell-level homogenization has not been carried out. Importance of cell-level homogenization is recently increasing since the pin-by-pin mesh core calculation in actual three-dimensional geometry, which is less approximate approach than current advanced nodal method, is getting feasible. Four homogenization methods were investigated in this paper; the flux-volume weighting, the generalized equivalence theory, the superhomogenization (SPH) method and the nonlinear iteration method. The last one, the nonlinear iteration method, was tested as the homogenization method for the first time. The calculations were carried out in simplified colorset assembly configurations of PWR, which are simulated by slab geometries, and homogenization performances were evaluated through comparison with the reference cell-heterogeneous calculations. The calculation results revealed that the generalized equivalence theory showed best performance. Though the nonlinear iteration method can significantly reduce homogenization error, its performance was not as good as that of the generalized equivalence theory. Through comparison of the results obtained by the generalized equivalence theory and the superhomogenization method, important byproduct was obtained; deficiency of the current superhomogenization method, which could be improved by incorporating the 'cell-level discontinuity factor between assemblies', was clarified

A combination of differential method and perturbation theory for the calculation of sensitivity coefficients

International Nuclear Information System (INIS)

Santos, Adimir dos; Borges, A.A.

2000-01-01

A new method for the calculation of sensitivity coefficients is developed. The new method is a combination of two methodologies used for calculating these coefficients, which are the differential and the generalized perturbation theory methods. The proposed method utilizes as integral parameter the average flux in an arbitrary region of the system. Thus, the sensitivity coefficient contains only the component corresponding to the neutron flux. To obtain the new sensitivity coefficient, the derivates of the integral parameter, φ(ξ), with respect to σ are calculated using the perturbation method and the functional derivates of this generic integral parameter with respect to σ and φ are calculated using the differential method. The new method merges the advantages of the differential and generalized perturbation theory methods and eliminates their disadvantages. (author)

New sampling method in continuous energy Monte Carlo calculation for pebble bed reactors

International Nuclear Information System (INIS)

Murata, Isao; Takahashi, Akito; Mori, Takamasa; Nakagawa, Masayuki.

1997-01-01

A pebble bed reactor generally has double heterogeneity consisting of two kinds of spherical fuel element. In the core, there exist many fuel balls piled up randomly in a high packing fraction. And each fuel ball contains a lot of small fuel particles which are also distributed randomly. In this study, to realize precise neutron transport calculation of such reactors with the continuous energy Monte Carlo method, a new sampling method has been developed. The new method has been implemented in the general purpose Monte Carlo code MCNP to develop a modified version MCNP-BALL. This method was validated by calculating inventory of spherical fuel elements arranged successively by sampling during transport calculation and also by performing criticality calculations in ordered packing models. From the results, it was confirmed that the inventory of spherical fuel elements could be reproduced using MCNP-BALL within a sufficient accuracy of 0.2%. And the comparison of criticality calculations in ordered packing models between MCNP-BALL and the reference method shows excellent agreement in neutron spectrum as well as multiplication factor. MCNP-BALL enables us to analyze pebble bed type cores such as PROTEUS precisely with the continuous energy Monte Carlo method. (author)

A method of paralleling computer calculation for two-dimensional kinetic plasma model

International Nuclear Information System (INIS)

Brazhnik, V.A.; Demchenko, V.V.; Dem'yanov, V.G.; D'yakov, V.E.; Ol'shanskij, V.V.; Panchenko, V.I.

1987-01-01

A method for parallel computer calculation and OSIRIS program complex realizing it and designed for numerical plasma simulation by the macroparticle method are described. The calculation can be carried out either with one or simultaneously with two computers BESM-6, that is provided by some package of interacting programs functioning in every computer. Program interaction in every computer is based on event techniques realized in OS DISPAK. Parallel computer calculation with two BESM-6 computers allows to accelerate the computation 1.5 times

Advances in methods of commercial FBR core characteristics analyses. Investigations of a treatment of the double-heterogeneity and a method to calculate homogenized control rod cross sections

Energy Technology Data Exchange (ETDEWEB)

Sugino, Kazuteru [Power Reactor and Nuclear Fuel Development Corp., Oarai, Ibaraki (Japan). Oarai Engineering Center; Iwai, Takehiko

1998-07-01

A standard data base for FBR core nuclear design is under development in order to improve the accuracy of FBR design calculation. As a part of the development, we investigated an improved treatment of double-heterogeneity and a method to calculate homogenized control rod cross sections in a commercial reactor geometry, for the betterment of the analytical accuracy of commercial FBR core characteristics. As an improvement in the treatment of double-heterogeneity, we derived a new method (the direct method) and compared both this and conventional methods with continuous energy Monte-Carlo calculations. In addition, we investigated the applicability of the reaction rate ratio preservation method as a advanced method to calculate homogenized control rod cross sections. The present studies gave the following information: (1) An improved treatment of double-heterogeneity: for criticality the conventional method showed good agreement with Monte-Carlo result within one sigma standard deviation; the direct method was consistent with conventional one. Preliminary evaluation of effects in core characteristics other than criticality showed that the effect of sodium void reactivity (coolant reactivity) due to the double-heterogeneity was large. (2) An advanced method to calculate homogenize control rod cross sections: for control rod worths the reaction rate ratio preservation method agreed with those produced by the calculations with the control rod heterogeneity included in the core geometry; in Monju control rod worth analysis, the present method overestimated control rod worths by 1 to 2% compared with the conventional method, but these differences were caused by more accurate model in the present method and it is considered that this method is more reliable than the conventional one. These two methods investigated in this study can be directly applied to core characteristics other than criticality or control rod worth. Thus it is concluded that these methods will

A new method for calculating gas saturation of low-resistivity shale gas reservoirs

Directory of Open Access Journals (Sweden)

Jinyan Zhang

2017-09-01

Full Text Available The Jiaoshiba shale gas field is located in the Fuling area of the Sichuan Basin, with the Upper Ordovician Wufeng–Lower Silurian Longmaxi Fm as the pay zone. At the bottom of the pay zone, a high-quality shale gas reservoir about 20 m thick is generally developed with high organic contents and gas abundance, but its resistivity is relatively low. Accordingly, the gas saturation calculated by formulas (e.g. Archie using electric logging data is often much lower than the experiment-derived value. In this paper, a new method was presented for calculating gas saturation more accurately based on non-electric logging data. Firstly, the causes for the low resistivity of shale gas reservoirs in this area were analyzed. Then, the limitation of traditional methods for calculating gas saturation based on electric logging data was diagnosed, and the feasibility of the neutron–density porosity overlay method was illustrated. According to the response characteristics of neutron, density and other porosity logging in shale gas reservoirs, a model for calculating gas saturation of shale gas was established by core experimental calibration based on the density logging value, the density porosity and the difference between density porosity and neutron porosity, by means of multiple methods (e.g. the dual-porosity overlay method by optimizing the best overlay coefficient. This new method avoids the effect of low resistivity, and thus can provide normal calculated gas saturation of high-quality shale gas reservoirs. It works well in practical application. This new method provides a technical support for the calculation of shale gas reserves in this area. Keywords: Shale gas, Gas saturation, Low resistivity, Non-electric logging, Volume density, Compensated neutron, Overlay method, Reserves calculation, Sichuan Basin, Jiaoshiba shale gas field

Three-dimensional static and dynamic reactor calculations by the nodal expansion method

International Nuclear Information System (INIS)

Christensen, B.

1985-05-01

This report reviews various method for the calculation of the neutron-flux- and power distribution in an nuclear reactor. The nodal expansion method (NEM) is especially described in much detail. The nodal expansion method solves the diffusion equation. In this method the reactor core is divided into nodes, typically 10 to 20 cm in each direction, and the average flux in each node is calculated. To obtain the coupling between the nodes the local flux inside each node is expressed by use of a polynomial expansion. The expansion is one-dimensional, so inside each node such three expansions occur. To calculate the expansion coefficients it is necessary that the polynomial expansion is a solution to the one-dimensional diffusion equation. When the one-dimensional diffusion equation is established a term with the transversal leakage occur, and this term is expanded after the same polynomials. The resulting equation system with the expansion coefficients as the unknowns is solved with weigthed residual technique. The nodal expansion method is built into a computer program (also called NEM), which is divided into two parts, one part for steady-state calculations and one part for dynamic calculations. It is possible to take advantage of symmetry properties of the reactor core. The program is very flexible with regard to the number of energy groups, the node size, the flux expansion order and the transverse leakage expansion order. The boundary of the core is described by albedos. The program and input to it are described. The program is tested on a number of examples extending from small theoretical one up to realistic reactor cores. Many calculations are done on the wellknown IAEA benchmark case. The calculations have tested the accuracy and the computing time for various node sizes and polynomial expansions. In the dynamic examples various strategies for variation of the time step-length have been tested. (author)

Closure and Sealing Design Calculation

International Nuclear Information System (INIS)

T. Lahnalampi; J. Case

2005-01-01

The purpose of the ''Closure and Sealing Design Calculation'' is to illustrate closure and sealing methods for sealing shafts, ramps, and identify boreholes that require sealing in order to limit the potential of water infiltration. In addition, this calculation will provide a description of the magma that can reduce the consequences of an igneous event intersecting the repository. This calculation will also include a listing of the project requirements related to closure and sealing. The scope of this calculation is to: summarize applicable project requirements and codes relating to backfilling nonemplacement openings, removal of uncommitted materials from the subsurface, installation of drip shields, and erecting monuments; compile an inventory of boreholes that are found in the area of the subsurface repository; describe the magma bulkhead feature and location; and include figures for the proposed shaft and ramp seals. The objective of this calculation is to: categorize the boreholes for sealing by depth and proximity to the subsurface repository; develop drawing figures which show the location and geometry for the magma bulkhead; include the shaft seal figures and a proposed construction sequence; and include the ramp seal figure and a proposed construction sequence. The intent of this closure and sealing calculation is to support the License Application by providing a description of the closure and sealing methods for the Safety Analysis Report. The closure and sealing calculation will also provide input for Post Closure Activities by describing the location of the magma bulkhead. This calculation is limited to describing the final configuration of the sealing and backfill systems for the underground area. The methods and procedures used to place the backfill and remove uncommitted materials (such as concrete) from the repository and detailed design of the magma bulkhead will be the subject of separate analyses or calculations. Post-closure monitoring will not

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Some Airfoils Using a Discrete Vortex Method

Science.gov (United States)

Riley, Donald R.

2015-01-01

This paper contains a collection of some results of four individual studies presenting calculated numerical values for airfoil aerodynamic stability derivatives in unseparated inviscid incompressible flow due separately to angle-of-attack, pitch rate, flap deflection, and airfoil camber using a discrete vortex method. Both steady conditions and oscillatory motion were considered. Variables include the number of vortices representing the airfoil, the pitch axis / moment center chordwise location, flap chord to airfoil chord ratio, and circular or parabolic arc camber. Comparisons with some experimental and other theoretical information are included. The calculated aerodynamic numerical results obtained using a limited number of vortices provided in each study compared favorably with thin airfoil theory predictions. Of particular interest are those aerodynamic results calculated herein (such as induced drag) that are not readily available elsewhere.

Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method

Science.gov (United States)

Gilbreth, C. N.; Alhassid, Y.

2015-03-01

Quantum Monte Carlo methods are powerful techniques for studying strongly interacting Fermi systems. However, implementing these methods on computers with finite-precision arithmetic requires careful attention to numerical stability. In the auxiliary-field Monte Carlo (AFMC) method, low-temperature or large-model-space calculations require numerically stabilized matrix multiplication. When adapting methods used in the grand-canonical ensemble to the canonical ensemble of fixed particle number, the numerical stabilization increases the number of required floating-point operations for computing observables by a factor of the size of the single-particle model space, and thus can greatly limit the systems that can be studied. We describe an improved method for stabilizing canonical-ensemble calculations in AFMC that exhibits better scaling, and present numerical tests that demonstrate the accuracy and improved performance of the method.

Weather data for simplified energy calculation methods. Volume II. Middle United States: TRY data

Energy Technology Data Exchange (ETDEWEB)

Olsen, A.R.; Moreno, S.; Deringer, J.; Watson, C.R.

1984-08-01

The objective of this report is to provide a source of weather data for direct use with a number of simplified energy calculation methods available today. Complete weather data for a number of cities in the United States are provided for use in the following methods: degree hour, modified degree hour, bin, modified bin, and variable degree day. This report contains sets of weather data for 22 cities in the continental United States using Test Reference Year (TRY) source weather data. The weather data at each city has been summarized in a number of ways to provide differing levels of detail necessary for alternative simplified energy calculation methods. Weather variables summarized include dry bulb and wet bulb temperature, percent relative humidity, humidity ratio, wind speed, percent possible sunshine, percent diffuse solar radiation, total solar radiation on horizontal and vertical surfaces, and solar heat gain through standard DSA glass. Monthly and annual summaries, in some cases by time of day, are available. These summaries are produced in a series of nine computer generated tables.

A Non-parametric Method for Calculating Conditional Stressed Value at Risk

Directory of Open Access Journals (Sweden)

Kohei Marumo

2017-01-01

Full Text Available We consider the Value at Risk (VaR of a portfolio under stressed conditions. In practice, the stressed VaR (sVaR is commonly calculated using the data set that includes the stressed period. It tells us how much the risk amount increases if we use the stressed data set. In this paper, we consider the VaR under stress scenarios. Technically, this can be done by deriving the distribution of profit or loss conditioned on the value of risk factors. We use two methods; the one that uses the linear model and the one that uses the Hermite expansion discussed by Marumo and Wolff (2013, 2016. Numerical examples shows that the method using the Hermite expansion is capable of capturing the non-linear effects such as correlation collapse and volatility clustering, which are often observed in the markets.

Moment methods with effective nuclear Hamiltonians; calculations of radial moments

International Nuclear Information System (INIS)

Belehrad, R.H.

1981-02-01

A truncated orthogonal polynomial expansion is used to evaluate the expectation value of the radial moments of the one-body density of nuclei. The expansion contains the configuration moments, , , and 2 >, where R/sup (k)/ is the operator for the k-th power of the radial coordinate r, and H is the effective nuclear Hamiltonian which is the sum of the relative kinetic energy operator and the Bruckner G matrix. Configuration moments are calculated using trace reduction formulae where the proton and neutron orbitals are treated separately in order to find expectation values of good total isospin. The operator averages are taken over many-body shell model states in the harmonic oscillator basis where all particles are active and single-particle orbitals through six major shells are included. The radial moment expectation values are calculated for the nuclei 16 O, 40 Ca, and 58 Ni and find that is usually the largest term in the expansion giving a large model space dependence to the results. For each of the 3 nuclei, a model space is found which gives the desired rms radius and then we find that the other 5 lowest moments compare favorably with other theoretical predictions. Finally, we use a method of Gordon (5) to employ the lowest 6 radial moment expectation values in the calculation of elastic electron scattering from these nuclei. For low to moderate momentum transfer, the results compare favorably with the experimental data

A new shielding calculation method for X-ray computed tomography regarding scattered radiation.

Science.gov (United States)

Watanabe, Hiroshi; Noto, Kimiya; Shohji, Tomokazu; Ogawa, Yasuyoshi; Fujibuchi, Toshioh; Yamaguchi, Ichiro; Hiraki, Hitoshi; Kida, Tetsuo; Sasanuma, Kazutoshi; Katsunuma, Yasushi; Nakano, Takurou; Horitsugi, Genki; Hosono, Makoto

2017-06-01

The goal of this study is to develop a more appropriate shielding calculation method for computed tomography (CT) in comparison with the Japanese conventional (JC) method and the National Council on Radiation Protection and Measurements (NCRP)-dose length product (DLP) method. Scattered dose distributions were measured in a CT room with 18 scanners (16 scanners in the case of the JC method) for one week during routine clinical use. The radiation doses were calculated for the same period using the JC and NCRP-DLP methods. The mean (NCRP-DLP-calculated dose)/(measured dose) ratios in each direction ranged from 1.7 ± 0.6 to 55 ± 24 (mean ± standard deviation). The NCRP-DLP method underestimated the dose at 3.4% in fewer shielding directions without the gantry and a subject, and the minimum (NCRP-DLP-calculated dose)/(measured dose) ratio was 0.6. The reduction factors were 0.036 ± 0.014 and 0.24 ± 0.061 for the gantry and couch directions, respectively. The (JC-calculated dose)/(measured dose) ratios ranged from 11 ± 8.7 to 404 ± 340. The air kerma scatter factor κ is expected to be twice as high as that calculated with the NCRP-DLP method and the reduction factors are expected to be 0.1 and 0.4 for the gantry and couch directions, respectively. We, therefore, propose a more appropriate method, the Japanese-DLP method, which resolves the issues of possible underestimation of the scattered radiation and overestimation of the reduction factors in the gantry and couch directions.

Complex Correlation Kohn-T Method of Calculating Total and Elastic Cross Sections. Part 1; Electron-Hydrogen Elastic Scattering

Science.gov (United States)

Bhatia, A. K.; Temkin, A.; Fisher, Richard R. (Technical Monitor)

2001-01-01

We report on the first part of a study of electron-hydrogen scattering, using a method which allows for the ab initio calculation of total and elastic cross sections at higher energies. In its general form the method uses complex 'radial' correlation functions, in a (Kohn) T-matrix formalism. The titled method, abbreviated Complex Correlation Kohn T (CCKT) method, is reviewed, in the context of electron-hydrogen scattering, including the derivation of the equation for the (complex) scattering function, and the extraction of the scattering information from the latter. The calculation reported here is restricted to S-waves in the elastic region, where the correlation functions can be taken, without loss of generality, to be real. Phase shifts are calculated using Hylleraas-type correlation functions with up to 95 terms. Results are rigorous lower bounds; they are in general agreement with those of Schwartz, but they are more accurate and outside his error bounds at a couple of energies,

Method for dose calculation in intracavitary irradiation of endometrical carcinoma

International Nuclear Information System (INIS)

Zevrieva, I.F.; Ivashchenko, N.T.; Musapirova, N.A.; Fel'dman, S.Z.; Sajbekov, T.S.

1979-01-01

A method for dose calculation for the conditions of intracavitary gamma therapy of endometrial carcinoma using spherical and linear 60 Co sources was elaborated. Calculations of dose rates for different amount and orientation of spherical radiation sources and for different planes were made with the aid of BEhSM-4M computer. Dosimet were made with the aid of BEhSM-4M computer. Dosimetric study of dose fields was made using a phantom imitating the real conditions of irradiation. Discrepancies between experimental and calculated values are within the limits of the experiment accuracy

Calculation of the collective mass-parameter including RPA corrections

International Nuclear Information System (INIS)

Pal, M.K.; Zawischa, D.; Speth, J.

1975-01-01

A derivation of the vibrational mass-parameter B is given which makes the consistency with RPA calculations explicit. The expected enhancement by the residual particle-hole and particle-particle interaction is demonstrated by solving the quasiparticle-RPA for deformed nuclei in the rare earth region. (orig.) [de

Development of concept and neutronic calculation method for large LMFBR core

International Nuclear Information System (INIS)

Shirakata, K.; Ishikawa, M.; Ikegami, T.; Sanda, T.; Kaneto, K.; Kawashima, M.; Kaise, Y.; Shirakawa, M.; Hibi, K.

1991-01-01

Presented in this paper is the state of the art of reactor physics R and Ds for the development of concept and neutronic calculation method for large Liquid Metal Fast Breeder Reactor (LMFBR) core. Physics characteristics of concepts for mixed oxide (MOX) fueled large FBR core were investigated by a series of benchmark critical experiments. Next, an adequacy and accuracy of the current neutronic calculation method was assessed by the experiments analyses, and then neutronic prediction accuracies by the method were evaluated for physics characteristics of the large core. Concerns on core development were discussed in terms of neutronics. (author)

Accurate pKa Calculation of the Conjugate Acids of Alkanolamines, Alkaloids and Nucleotide Bases by Quantum Chemical Methods

NARCIS (Netherlands)

Gangarapu, S.; Marcelis, A.T.M.; Zuilhof, H.

2013-01-01

The pKa of the conjugate acids of alkanolamines, neurotransmitters, alkaloid drugs and nucleotide bases are calculated with density functional methods (B3LYP, M08-HX and M11-L) and ab initio methods (SCS-MP2, G3). Implicit solvent effects are included with a conductor-like polarizable continuum

A finite element method for a time dependence soil-structure interactions calculations

International Nuclear Information System (INIS)

Ni, X.M.; Gantenbein, F.; Petit, M.

1989-01-01

The method which is proposed is based on a finite element modelisation for the soil and the structure and a time history calculation. It has been developed for plane and axisymmetric geometries. The principle of this method will be presented, then applications will be given, first to a linear calculation for which results will be compared to those obtained by standard methods. Then results for a non linear behavior will be described [fr

Heterogeneous Calculation of {epsilon}

Energy Technology Data Exchange (ETDEWEB)

Jonsson, Alf

1961-02-15

A heterogeneous method of calculating the fast fission factor given by Naudet has been applied to the Carlvik - Pershagen definition of {epsilon}. An exact calculation of the collision probabilities is included in the programme developed for the Ferranti - Mercury computer.

Heterogeneous Calculation of ε

International Nuclear Information System (INIS)

Jonsson, Alf

1961-02-01

A heterogeneous method of calculating the fast fission factor given by Naudet has been applied to the Carlvik - Pershagen definition of ε. An exact calculation of the collision probabilities is included in the programme developed for the Ferranti - Mercury computer

Development of High Precision Tsunami Runup Calculation Method Coupled with Structure Analysis

Science.gov (United States)

Arikawa, Taro; Seki, Katsumi; Chida, Yu; Takagawa, Tomohiro; Shimosako, Kenichiro

2017-04-01

The 2011 Great East Japan Earthquake (GEJE) has shown that tsunami disasters are not limited to inundation damage in a specified region, but may destroy a wide area, causing a major disaster. Evaluating standing land structures and damage to them requires highly precise evaluation of three-dimensional fluid motion - an expensive process. Our research goals were thus to develop a coupling STOC-CADMAS (Arikawa and Tomita, 2016) coupling with the structure analysis (Arikawa et. al., 2009) to efficiently calculate all stages from tsunami source to runup including the deformation of structures and to verify their applicability. We also investigated the stability of breakwaters at Kamaishi Bay. Fig. 1 shows the whole of this calculation system. The STOC-ML simulator approximates pressure by hydrostatic pressure and calculates the wave profiles based on an equation of continuity, thereby lowering calculation cost, primarily calculating from a e epi center to the shallow region. As a simulator, STOC-IC solves pressure based on a Poisson equation to account for a shallower, more complex topography, but reduces computation cost slightly to calculate the area near a port by setting the water surface based on an equation of continuity. CS3D also solves a Navier-Stokes equation and sets the water surface by VOF to deal with the runup area, with its complex surfaces of overflows and bores. STR solves the structure analysis including the geo analysis based on the Biot's formula. By coupling these, it efficiently calculates the tsunami profile from the propagation to the inundation. The numerical results compared with the physical experiments done by Arikawa et. al.,2012. It was good agreement with the experimental ones. Finally, the system applied to the local situation at Kamaishi bay. The almost breakwaters were washed away, whose situation was similar to the damage at Kamaishi bay. REFERENCES T. Arikawa and T. Tomita (2016): "Development of High Precision Tsunami Runup

Application of the Activity-Based Costing Method for Unit-Cost Calculation in a Hospital

Science.gov (United States)

Javid, Mahdi; Hadian, Mohammad; Ghaderi, Hossein; Ghaffari, Shahram; Salehi, Masoud

2016-01-01

Background: Choosing an appropriate accounting system for hospital has always been a challenge for hospital managers. Traditional cost system (TCS) causes cost distortions in hospital. Activity-based costing (ABC) method is a new and more effective cost system. Objective: This study aimed to compare ABC with TCS method in calculating the unit cost of medical services and to assess its applicability in Kashani Hospital, Shahrekord City, Iran. Methods: This cross-sectional study was performed on accounting data of Kashani Hospital in 2013. Data on accounting reports of 2012 and other relevant sources at the end of 2012 were included. To apply ABC method, the hospital was divided into several cost centers and five cost categories were defined: wage, equipment, space, material, and overhead costs. Then activity centers were defined. ABC method was performed into two phases. First, the total costs of cost centers were assigned to activities by using related cost factors. Then the costs of activities were divided to cost objects by using cost drivers. After determining the cost of objects, the cost price of medical services was calculated and compared with those obtained from TCS. Results: The Kashani Hospital had 81 physicians, 306 nurses, and 328 beds with the mean occupancy rate of 67.4% during 2012. Unit cost of medical services, cost price of occupancy bed per day, and cost per outpatient service were calculated. The total unit costs by ABC and TCS were respectively 187.95 and 137.70 USD, showing 50.34 USD more unit cost by ABC method. ABC method represented more accurate information on the major cost components. Conclusion: By utilizing ABC, hospital managers have a valuable accounting system that provides a true insight into the organizational costs of their department. PMID:26234974

Radioactive cloud dose calculations

International Nuclear Information System (INIS)

Healy, J.W.

1984-01-01

Radiological dosage principles, as well as methods for calculating external and internal dose rates, following dispersion and deposition of radioactive materials in the atmosphere are described. Emphasis has been placed on analytical solutions that are appropriate for hand calculations. In addition, the methods for calculating dose rates from ingestion are discussed. A brief description of several computer programs are included for information on radionuclides. There has been no attempt to be comprehensive, and only a sampling of programs has been selected to illustrate the variety available

Charged-particle calculations using Boltzmann transport methods

International Nuclear Information System (INIS)

Hoffman, T.J.; Dodds, H.L. Jr.; Robinson, M.T.; Holmes, D.K.

1981-01-01

Several aspects of radiation damage effects in fusion reactor neutron and ion irradiation environments are amenable to treatment by transport theory methods. In this paper, multigroup transport techniques are developed for the calculation of charged particle range distributions, reflection coefficients, and sputtering yields. The Boltzmann transport approach can be implemented, with minor changes, in standard neutral particle computer codes. With the multigroup discrete ordinates code, ANISN, determination of ion and target atom distributions as functions of position, energy, and direction can be obtained without the stochastic error associated with atomistic computer codes such as MARLOWE and TRIM. With the multigroup Monte Carlo code, MORSE, charged particle effects can be obtained for problems associated with very complex geometries. Results are presented for several charged particle problems. Good agreement is obtained between quantities calculated with the multigroup approach and those obtained experimentally or by atomistic computer codes

The method of calculation of pipelines laid on supports

Directory of Open Access Journals (Sweden)

Benin D.M.

2017-08-01

Full Text Available this article focuses on the issue of laying pipelines on supports and the method of calculation of vertical and horizontal loads acting on the support. As pipelines can be water piping systems, heat networks, oil and mazout lines, condensate lines, steam lines, etc. this article describes the calculations of supports for pipelines laid above ground, in crowded channels, premises, on racks, in impassable channels, hanging supports, etc. The paper explores recommendations for placement of the supports on the route of the pipelines, calculation of loads on rotating and stationary supports of pipelines; inspection of stresses in the metal pipe, resulting from elongation of the piping from the temperature from the thermal expansion of the metal during operation.

Conjugate-gradient optimization method for orbital-free density functional calculations.

Science.gov (United States)

Jiang, Hong; Yang, Weitao

2004-08-01

Orbital-free density functional theory as an extension of traditional Thomas-Fermi theory has attracted a lot of interest in the past decade because of developments in both more accurate kinetic energy functionals and highly efficient numerical methodology. In this paper, we developed a conjugate-gradient method for the numerical solution of spin-dependent extended Thomas-Fermi equation by incorporating techniques previously used in Kohn-Sham calculations. The key ingredient of the method is an approximate line-search scheme and a collective treatment of two spin densities in the case of spin-dependent extended Thomas-Fermi problem. Test calculations for a quartic two-dimensional quantum dot system and a three-dimensional sodium cluster Na216 with a local pseudopotential demonstrate that the method is accurate and efficient. (c) 2004 American Institute of Physics.

Calculation of degenerated Eigenmodes with modified power method

Energy Technology Data Exchange (ETDEWEB)

Zhang, Peng; Lee, Hyun Suk; Lee, Deok Jung [School of Mechanical and Nuclear Engineering, Ulsan National Institute of Science and Technology, Ulsan (Korea, Republic of)

2017-02-15

The modified power method has been studied by many researchers to calculate the higher Eigenmodes and accelerate the convergence of the fundamental mode. Its application to multidimensional problems may be unstable due to degenerated or near-degenerated Eigenmodes. Complex Eigenmode solutions are occasionally encountered in such cases, and the shapes of the corresponding eigenvectors may change during the simulation. These issues must be addressed for the successful implementation of the modified power method. Complex components are examined and an approximation method to eliminate the usage of the complex numbers is provided. A technique to fix the eigenvector shapes is also provided. The performance of the methods for dealing with those aforementioned problems is demonstrated with two dimensional one group and three dimensional one group homogeneous diffusion problems.

Method for calculating carbon footprint of cattle feeds – including contribution from soil carbon changes and use of cattle manure

DEFF Research Database (Denmark)

Mogensen, Lisbeth; Kristensen, Troels; Nguyen, T Lan T

2014-01-01

fodder crop, an individual production scheme was set up as the basis for calculating the carbon footprint (CF). In the calculations, all fodder crops were fertilized by artificial fertilizer based on the assumption that the environmental burden of using manure is related to the livestock production......Greenhouse gas emissions (GHG) related to feed production is one of the hotspots in livestock production. The aim of this paper was to estimate the carbon footprint of different feedstuffs for dairy cattle using life cycle assessment (LCA). The functional unit was ‘1 kg dry matter (DM) of feed...

The Impact of Harmonics Calculation Methods on Power Quality Assessment in Wind Farms

DEFF Research Database (Denmark)

Kocewiak, Lukasz Hubert; Hjerrild, Jesper; Bak, Claus Leth

2010-01-01

Different methods of calculating harmonics in measurements obtained from offshore wind farms are shown in this paper. Appropriate data processing methods are suggested for harmonics with different origin and nature. Enhancements of discrete Fourier transform application in order to reduce...... measurement data processing errors are proposed and compared with classical methods. Comparison of signal processing methods for harmonic studies is presented and application dependent on harmonics origin and nature recommended. Certain aspects related to magnitude and phase calculation in stationary...... measurement data are analysed and described. Qualitative indices of measurement data harmonic analysis in order to assess the calculation accuracy are suggested and used....

Calculation of seismic response of a flexible rotor by complex modal method, 1

International Nuclear Information System (INIS)

Azuma, Takao; Saito, Shinobu

1984-01-01

In rotary machines, at the time of earthquakes, whether the rotating part and stationary part touch or whether the bearings and seals are damaged or not are problems. In order to examine these problems, it is necessary to analyze the seismic response of a rotary shaft or sometimes a casing system. But the conventional analysis methods are unsatisfactory. Accordingly, in the case of a general shaft system supported with slide bearings and on which gyro effect acts, complex modal method must be used. This calculation method is explained in detail in the book of Lancaster, however, when this method is applied to the seismic response of rotary shafts, the calculation time is considerably different according to the method of final integration. In this study, good results were obtained when the method which did not depend on numerical integration was attempted. The equation of motion and its solution, the displacement vector of a foundation, the verification of the calculation program and the example of calculating the seismic response of two coupled rotor shafts are reported. (Kako, I.)

Nucleon-deuteron breakup quantities calculated with separable interactions including tensor forces and P-wave interactions

International Nuclear Information System (INIS)

Bruinsma, J.; Wageningen, R. van

1977-01-01

Nucleon-deuteron breakup calculations at a nucleon bombarding energy of 22.7 MeV have been performed with separable interactions including a tensor force and P-wave interactions. Differential cross sections and a selection of polarization quantities have been computed for special regions of the phase space. The influence of a tensor force and P-wave interactions on the differential cross section is of the order of 20%. Large discrepancies between theory and experiment occur for the vector analyzing powers, both for the kinematically complete and for the incomplete situation. The calculations show that there are kinematical situations in which the differential cross sections and the tensor analyzing powers are sufficiently large to make measurements feasible. (Auth.)

Dose calculation methods in photon beam therapy using energy deposition kernels

International Nuclear Information System (INIS)

Ahnesjoe, A.

1991-01-01

The problem of calculating accurate dose distributions in treatment planning of megavoltage photon radiation therapy has been studied. New dose calculation algorithms using energy deposition kernels have been developed. The kernels describe the transfer of energy by secondary particles from a primary photon interaction site to its surroundings. Monte Carlo simulations of particle transport have been used for derivation of kernels for primary photon energies form 0.1 MeV to 50 MeV. The trade off between accuracy and calculational speed has been addressed by the development of two algorithms; one point oriented with low computional overhead for interactive use and one for fast and accurate calculation of dose distributions in a 3-dimensional lattice. The latter algorithm models secondary particle transport in heterogeneous tissue by scaling energy deposition kernels with the electron density of the tissue. The accuracy of the methods has been tested using full Monte Carlo simulations for different geometries, and found to be superior to conventional algorithms based on scaling of broad beam dose distributions. Methods have also been developed for characterization of clinical photon beams in entities appropriate for kernel based calculation models. By approximating the spectrum as laterally invariant, an effective spectrum and dose distribution for contaminating charge particles are derived form depth dose distributions measured in water, using analytical constraints. The spectrum is used to calculate kernels by superposition of monoenergetic kernels. The lateral energy fluence distribution is determined by deconvolving measured lateral dose distributions by a corresponding pencil beam kernel. Dose distributions for contaminating photons are described using two different methods, one for estimation of the dose outside of the collimated beam, and the other for calibration of output factors derived from kernel based dose calculations. (au)

Higher order methods for burnup calculations with Bateman solutions

International Nuclear Information System (INIS)

Isotalo, A.E.; Aarnio, P.A.

2011-01-01

Highlights: → Average microscopic reaction rates need to be estimated at each step. → Traditional predictor-corrector methods use zeroth and first order predictions. → Increasing predictor order greatly improves results. → Increasing corrector order does not improve results. - Abstract: A group of methods for burnup calculations solves the changes in material compositions by evaluating an explicit solution to the Bateman equations with constant microscopic reaction rates. This requires predicting representative averages for the one-group cross-sections and flux during each step, which is usually done using zeroth and first order predictions for their time development in a predictor-corrector calculation. In this paper we present the results of using linear, rather than constant, extrapolation on the predictor and quadratic, rather than linear, interpolation on the corrector. Both of these are done by using data from the previous step, and thus do not affect the stepwise running time. The methods were tested by implementing them into the reactor physics code Serpent and comparing the results from four test cases to accurate reference results obtained with very short steps. Linear extrapolation greatly improved results for thermal spectra and should be preferred over the constant one currently used in all Bateman solution based burnup calculations. The effects of using quadratic interpolation on the corrector were, on the other hand, predominantly negative, although not enough so to conclusively decide between the linear and quadratic variants.

A quick method to calculate QTL confidence interval

Indian Academy of Sciences (India)

2011-08-19

Aug 19, 2011 ... experimental design and analysis to reveal the real molecular nature of the ... strap sample form the bootstrap distribution of QTL location. The 2.5 and ..... ative probability to harbour a true QTL, hence x-LOD rule is not stable ... Darvasi A. and Soller M. 1997 A simple method to calculate resolv- ing power ...

Comparison of the accuracy of three angiographic methods for calculating left ventricular volume measurement

International Nuclear Information System (INIS)

Hu Lin; Cui Wei; Shi Hanwen; Tian Yingping; Wang Weigang; Feng Yanguang; Huang Xueyan; Liu Zhisheng

2003-01-01

Objective: To compare the relative accuracy of three methods measuring left ventricular volume by X-ray ventriculography: single plane area-length method, biplane area-length method, and single-plane Simpson's method. Methods: Left ventricular casts were obtained within 24 hours after death from 12 persons who died from non-cardiac causes. The true left ventricular cast volume was measured by water displacement. The calculated volume of the casts was obtained with 3 angiographic methods, i.e., single-plane area-length method, biplane area-length method, and single-plane Simpson's method. Results: The actual average volume of left ventricular casts was (61.17±26.49) ml. The left ventricular volume was averagely (97.50±35.56) ml with single plane area-length method, (90.51±36.33) ml with biplane area-length method, and (65.00± 23.63) ml with single-plane Simpson's method. The left ventricular volumes calculated with single-plane and biplane area-length method were significantly larger than that the actual volumes (P 0.05). The left ventricular volumes calculated with single-plane and biplane area-length method were significantly larger than those calculated with single-plane Simpson's method (P 0.05). The over-estimation of left ventricular volume by single plane area-length method (36.34±17.98) ml and biplane area-length method (29.34±15.59) ml was more obvious than that calculated by single-plane Simpson's method (3.83±8.48) ml. Linear regression analysis showed that there was close correlations between left ventricular volumes calculated with single plane area-length method, biplane area-length method, Simpson's method and the true volume (all r>0.98). Conclusion: Single-plane Simpson's method is more accurate than single plane area-length method and biplane area-length method for left ventricular volume measurement; however, both the single-plane and biplane area-length methods could be used in clinical practice, especially in those imaging modality

New approach to invariant-embedding methods in reactor physics calculations

International Nuclear Information System (INIS)

Forsbacka, M.J.; Rydin, R.A.

1997-01-01

Invariant-embedding methods offer an alternative approach to modeling physical phenomena and solving mathematical problems. Invariant embedding allows one to express traditional boundary-value problems as initial-value problems. In doing this, one effectively reformulates a problem to be solved in terms of an embedding parameter. In this paper, a hybrid method that consists of Monte Carlo-generated response functions that describe the neutronic properties of local spatial cells are coupled together in a global reactor model using the invariant embedding methodology, where the system multiplication factor k eff is used as the embedding parameter. Thus, k eff is computed directly rather than as the result of a secondary eigenvalue calculation. Because the response functions can represent any arbitrary material distribution within a local cell, this method shows promise to accurately assess the change in reactivity due to core disruptive accidents and other changes in system configuration such as changing control rod positions. This paper reports a series of proof-of-concept calculations that assess this method

Wave resistance calculation method combining Green functions based on Rankine and Kelvin source

Directory of Open Access Journals (Sweden)

LI Jingyu

2017-12-01

Full Text Available [Ojectives] At present, the Boundary Element Method(BEM of wave-making resistance mostly uses a model in which the velocity distribution near the hull is solved first, and the pressure integral is then calculated using the Bernoulli equation. However,the process of this model of wave-making resistance is complex and has low accuracy.[Methods] To address this problem, the present paper deduces a compound method for the quick calculation of ship wave resistance using the Rankine source Green function to solve the hull surface's source density, and combining the Lagally theorem concerning source point force calculation based on the Kelvin source Green function so as to solve the wave resistance. A case for the Wigley model is given.[Results] The results show that in contrast to the thin ship method of the linear wave resistance theorem, this method has higher precision, and in contrast to the method which completely uses the Kelvin source Green function, this method has better computational efficiency.[Conclusions] In general, the algorithm in this paper provides a compromise between precision and efficiency in wave-making resistance calculation.

Approximate calculation method for integral of mean square value of nonstationary response

International Nuclear Information System (INIS)

Aoki, Shigeru; Fukano, Azusa

2010-01-01

The response of the structure subjected to nonstationary random vibration such as earthquake excitation is nonstationary random vibration. Calculating method for statistical characteristics of such a response is complicated. Mean square value of the response is usually used to evaluate random response. Integral of mean square value of the response corresponds to total energy of the response. In this paper, a simplified calculation method to obtain integral of mean square value of the response is proposed. As input excitation, nonstationary white noise and nonstationary filtered white noise are used. Integrals of mean square value of the response are calculated for various values of parameters. It is found that the proposed method gives exact value of integral of mean square value of the response.

An improved correlated sampling method for calculating correction factor of detector

International Nuclear Information System (INIS)

Wu Zhen; Li Junli; Cheng Jianping

2006-01-01

In the case of a small size detector lying inside a bulk of medium, there are two problems in the correction factors calculation of the detectors. One is that the detector is too small for the particles to arrive at and collide in; the other is that the ratio of two quantities is not accurate enough. The method discussed in this paper, which combines correlated sampling with modified particle collision auto-importance sampling, and has been realized on the MCNP-4C platform, can solve these two problems. Besides, other 3 variance reduction techniques are also combined with correlated sampling respectively to calculate a simple calculating model of the correction factors of detectors. The results prove that, although all the variance reduction techniques combined with correlated sampling can improve the calculating efficiency, the method combining the modified particle collision auto-importance sampling with the correlated sampling is the most efficient one. (authors)

A New Method to Calculate Internal Rate of Return

Directory of Open Access Journals (Sweden)

azadeh zandi

2015-09-01

Full Text Available A number of methods have been developed to choose the best capital investment projects such as net present value, internal rate of return and etc. Internal rate of return method is probably the most popular method among managers and investors. But despite the popularity there are serious drawbacks and limitations in this method. After decades of efforts made by economists and experts to improve the method and its shortcomings, Magni in 2010 has revealed a new approach that can solves the most of internal rate of return method problems. This paper present a new method which is originated from Magni’s approach but has much more simple calculations and can resolve all the drawbacks of internal rate of return method.

An approximate method to calculate ionization of LTE and non-LTE plasma

International Nuclear Information System (INIS)

Zhang Jun; Gu Peijun

1987-01-01

When matter, especially high Z element, is heated to high temperature, it will be ionized many times. The degree of ionization has a strong effect on many plasma properties. So an approximate method to calculate the mean ionization degree is needed for solving many practical problems. An analytical expression which is convenient for the approximate numerical calculation is given by fitting it to the scaling law and numerical results of the ionization potential of Thomas-Fermi statistical model. In LTE case, the ionization degree of Au calculated by using the approximate method is in agreement with that of the average ion model. By extending the approximate method to non-LTE case, the ionization degree of Au is similarly calculated according to Corona model and Collision-Radiatoin model(C-R). The results of Corona model agree with the published data quite well, while the results of C-R approach those of Corona model as the density is reduced and approach those of LTE as the density is increased. Finally, all approximately calculated results of ionization degree of Au and the comparision of them are given in figures and tables

Application of the Activity-Based Costing Method for Unit-Cost Calculation in a Hospital.

Science.gov (United States)

Javid, Mahdi; Hadian, Mohammad; Ghaderi, Hossein; Ghaffari, Shahram; Salehi, Masoud

2015-05-17

Choosing an appropriate accounting system for hospital has always been a challenge for hospital managers. Traditional cost system (TCS) causes cost distortions in hospital. Activity-based costing (ABC) method is a new and more effective cost system. This study aimed to compare ABC with TCS method in calculating the unit cost of medical services and to assess its applicability in Kashani Hospital, Shahrekord City, Iran.‎ This cross-sectional study was performed on accounting data of Kashani Hospital in 2013. Data on accounting reports of 2012 and other relevant sources at the end of 2012 were included. To apply ABC method, the hospital was divided into several cost centers and five cost categories were defined: wage, equipment, space, material, and overhead costs. Then activity centers were defined. ABC method was performed into two phases. First, the total costs of cost centers were assigned to activities by using related cost factors. Then the costs of activities were divided to cost objects by using cost drivers. After determining the cost of objects, the cost price of medical services was calculated and compared with those obtained from TCS.‎ The Kashani Hospital had 81 physicians, 306 nurses, and 328 beds with the mean occupancy rate of 67.4% during 2012. Unit cost of medical services, cost price of occupancy bed per day, and cost per outpatient service were calculated. The total unit costs by ABC and TCS were respectively 187.95 and 137.70 USD, showing 50.34 USD more unit cost by ABC method. ABC method represented more accurate information on the major cost components. By utilizing ABC, hospital managers have a valuable accounting system that provides a true insight into the organizational costs of their department.

Further Stable methods for the calculation of partition functions

International Nuclear Information System (INIS)

Wilson, B G; Gilleron, F; Pain, J

2007-01-01

The extension to recursion over holes of the Gilleron and Pain method for calculating partition functions of a canonical ensemble of non-interacting bound electrons is presented as well as a generalization for the efficient computation of collisional line broadening

Use of the Local Variation Methods for Nuclear Design Calculations

International Nuclear Information System (INIS)

Zhukov, A.I.

2006-01-01

A new problem-solving method for steady-state equations, which describe neutron diffusion, is presented. The method bases on a variation principal for steady-state diffusion equations and direct search the minimum of a corresponding functional. Benchmark problem calculation for power of fuel assemblies show ∼ 2% relative accuracy

The GMD Method for Inductance Calculation Applied to Conductors with Skin Effect

Directory of Open Access Journals (Sweden)

H. A. Aebischer

2017-09-01

Full Text Available The GMD method (geometric mean distance to calculate inductance offers undoubted advantages over other methods. But so far it seemed to be limited to the case where the current is uniformly distributed over the cross section of the conductor, i.e. to DC (direct current. In this paper, the definition of the GMD is extended to include cases of nonuniform distribution observed at higher frequencies as the result of skin effect. An exact relation between the GMD and the internal inductance per unit length for infinitely long conductors of circularly symmetric cross section is derived. It enables much simpler derivations of Maxwell’s analytical expressions for the GMD of circular and annular disks than were known before. Its salient application, however, is the derivation of exact expressions for the GMD of infinitely long round wires and tubular conductors with skin effect. These expressions are then used to verify the consistency of the extended definition of the GMD. Further, approximate formulae for the GMD of round wires with skin effect based on elementary functions are discussed. Total inductances calculated with the help of the derived formulae for the GMD with and without skin effect are compared to measurement results from the literature. For conductors of square cross section, an analytical approximation for the GMD with skin effect based on elementary functions is presented. It is shown that it allows to calculate the total inductance of such conductors for frequencies from DC up to 25 GHz to a precision of better than 1 %.

The multigrid method for reactor calculations

International Nuclear Information System (INIS)

Douglas, S.R.

1991-07-01

Iterative solutions to linear systems of equations are discussed. The emphasis is on the concepts that affect convergence rates of these solution methods. The multigrid method is described, including the smoothing property, restriction, and prolongation. A simple example is used to illustrate the ideas

Advanced numerical methods for three dimensional two-phase flow calculations in PWR

International Nuclear Information System (INIS)

Toumi, I.; Gallo, D.; Royer, E.

1997-01-01

This paper is devoted to new numerical methods developed for three dimensional two-phase flow calculations. These methods are finite volume numerical methods. They are based on an extension of Roe's approximate Riemann solver to define convective fluxes versus mean cell quantities. To go forward in time, a linearized conservative implicit integrating step is used, together with a Newton iterative method. We also present here some improvements performed to obtain a fully implicit solution method that provides fast running steady state calculations. This kind of numerical method, which is widely used for fluid dynamic calculations, is proved to be very efficient for the numerical solution to two-phase flow problems. This numerical method has been implemented for the three dimensional thermal-hydraulic code FLICA-4 which is mainly dedicated to core thermal-hydraulic transient and steady-state analysis. Hereafter, we will also find some results obtained for the EPR reactor running in a steady-state at 60% of nominal power with 3 pumps out of 4, and a thermal-hydraulic core analysis for a 1300 MW PWR at low flow steam-line-break conditions. (author)

Methods for U.S. shielding calculations: applications to FFTF and CRBR designs

International Nuclear Information System (INIS)

Engle, W.W. Jr.; Mynatt, F.R.; Disney, R.K.

1978-01-01

The primary components of the U.S. reactor shielding methodology consist of: (1) computer code systems based on discrete ordinates or Monte Carlo radiation transport calculational methods; (2) a data base of neutron and gamma-ray interaction and gamma-ray-production cross sections used as input in the codes; (3) a capability for processing the cross sections into multigroup or point energy formats as required by the codes; (4) large-scale integral shielding experiments designed to test cross-section data or techniques utilized in the calculations; and (5) a ''sensitivity'' analysis capability that can identify the most important interactions in a transport calculation and assign uncertainties to the calculated result that are based on uncertainties in all of the input data. The required accuracy for the methodology is to within 5 to 10% for responses at locations near the core to within a factor of 2 for responses at distant locations. Under these criteria, the methodology has proved to be adequate for in-vessel LMFBR calculations of neutron transport through deep sodium and thick iron and stainless steel shields, of neutron streaming through lower axial coolant channels and primary pipe chaseways, and of the effects of fuel stored within the reactor vessel. For ex-vessel LMFBR problems, the methodology requires considerable improvement, the areas of concern including neutron streaming through heating and ventilation ducts, through the cavity surrounding the reactor vessel, and through gaps around rotating plugs in the reactor heat, as well as gamma-ray streaming through plant shield penetrations

Dose calculation method with 60-cobalt gamma rays in total body irradiation

International Nuclear Information System (INIS)

Scaff, Luiz Alberto Malaguti

2001-01-01

Physical factors associated to total body irradiation using 60 Co gamma rays beams, were studied in order to develop a calculation method of the dose distribution that could be reproduced in any radiotherapy center with good precision. The method is based on considering total body irradiation as a large and irregular field with heterogeneities. To calculate doses, or doses rates, of each area of interest (head, thorax, thigh, etc.), scattered radiation is determined. It was observed that if dismagnified fields were considered to calculate the scattered radiation, the resulting values could be applied on a projection to the real size to obtain the values for dose rate calculations. In a parallel work it was determined the variation of the dose rate in the air, for the distance of treatment, and for points out of the central axis. This confirm that the use of the inverse square law is not valid. An attenuation curve for a broad beam was also determined in order to allow the use of absorbers. In this work all the adapted formulas for dose rate calculations in several areas of the body are described, as well time/dose templates sheets for total body irradiation. The in vivo dosimetry, proved that either experimental or calculated dose rate values (achieved by the proposed method), did not have significant discrepancies. (author)

Monte Carlo methods to calculate impact probabilities

Science.gov (United States)

Rickman, H.; Wiśniowski, T.; Wajer, P.; Gabryszewski, R.; Valsecchi, G. B.

2014-09-01

Context. Unraveling the events that took place in the solar system during the period known as the late heavy bombardment requires the interpretation of the cratered surfaces of the Moon and terrestrial planets. This, in turn, requires good estimates of the statistical impact probabilities for different source populations of projectiles, a subject that has received relatively little attention, since the works of Öpik (1951, Proc. R. Irish Acad. Sect. A, 54, 165) and Wetherill (1967, J. Geophys. Res., 72, 2429). Aims: We aim to work around the limitations of the Öpik and Wetherill formulae, which are caused by singularities due to zero denominators under special circumstances. Using modern computers, it is possible to make good estimates of impact probabilities by means of Monte Carlo simulations, and in this work, we explore the available options. Methods: We describe three basic methods to derive the average impact probability for a projectile with a given semi-major axis, eccentricity, and inclination with respect to a target planet on an elliptic orbit. One is a numerical averaging of the Wetherill formula; the next is a Monte Carlo super-sizing method using the target's Hill sphere. The third uses extensive minimum orbit intersection distance (MOID) calculations for a Monte Carlo sampling of potentially impacting orbits, along with calculations of the relevant interval for the timing of the encounter allowing collision. Numerical experiments are carried out for an intercomparison of the methods and to scrutinize their behavior near the singularities (zero relative inclination and equal perihelion distances). Results: We find an excellent agreement between all methods in the general case, while there appear large differences in the immediate vicinity of the singularities. With respect to the MOID method, which is the only one that does not involve simplifying assumptions and approximations, the Wetherill averaging impact probability departs by diverging toward

Reactor calculation in coarse mesh by finite element method applied to matrix response method

International Nuclear Information System (INIS)

Nakata, H.

1982-01-01

The finite element method is applied to the solution of the modified formulation of the matrix-response method aiming to do reactor calculations in coarse mesh. Good results are obtained with a short running time. The method is applicable to problems where the heterogeneity is predominant and to problems of evolution in coarse meshes where the burnup is variable in one same coarse mesh, making the cross section vary spatially with the evolution. (E.G.) [pt

Coarse-mesh method for multidimensional, mixed-lattice diffusion calculations

International Nuclear Information System (INIS)

Dodds, H.L. Jr.; Honeck, H.C.; Hostetler, D.E.

1977-01-01

A coarse-mesh finite difference method has been developed for multidimensional, mixed-lattice reactor diffusion calculations, both statics and kinetics, in hexagonal geometry. Results obtained with the coarse-mesh (CM) method have been compared with a conventional mesh-centered finite difference method and with experiment. The results of this comparison indicate that the accuracy of the CM method for highly heterogeneous (mixed) lattices using one point per hexagonal mesh element (''hex'') is about the same as the conventional method with six points per hex. Furthermore, the computing costs (i.e., central processor unit time and core storage requirements) of the CM method with one point per hex are about the same as the conventional method with one point per hex

A comparison of published methods of calculation of defect significance

International Nuclear Information System (INIS)

Ingham, T.; Harrison, R.P.

1982-01-01

This paper presents some of the results obtained in a round-robin calculational exercise organised by the OECD Committee on the Safety of Nuclear Installations (CSNI). The exercise was initiated to examine practical aspects of using documented elastic-plastic fracture mechanics methods to calculate defect significance. The extent to which the objectives of the exercise were met is illustrated using solutions to 'standard' problems produced by UKAEA and CEGB using the methods given in ASME XI, Appendix A, BSI PD6493, and the CEGB R/H/R6 Document. Differences in critical or tolerable defect size defined using these procedures are examined in terms of their different treatments and reasons for discrepancies are discussed. (author)

Thermal disadvantage factor calculation by the multiregion collision probability method

International Nuclear Information System (INIS)

Ozgener, B.; Ozgener, H.A.

2004-01-01

A multi-region collision probability formulation that is capable of applying white boundary condition directly is presented and applied to thermal neutron transport problems. The disadvantage factors computed are compared with their counterparts calculated by S N methods with both direct and indirect application of white boundary condition. The results of the ABH and collision probability method with indirect application of white boundary condition are also considered and comparisons with benchmark Monte Carlo results are carried out. The studies show that the proposed formulation is capable of calculating thermal disadvantage factor with sufficient accuracy without resorting to the fictitious scattering outer shell approximation associated with the indirect application of the white boundary condition in collision probability solutions

Kernel polynomial method for a nonorthogonal electronic-structure calculation of amorphous diamond

International Nuclear Information System (INIS)

Roeder, H.; Silver, R.N.; Drabold, D.A.; Dong, J.J.

1997-01-01

The Kernel polynomial method (KPM) has been successfully applied to tight-binding electronic-structure calculations as an O(N) method. Here we extend this method to nonorthogonal basis sets with a sparse overlap matrix S and a sparse Hamiltonian H. Since the KPM method utilizes matrix vector multiplications it is necessary to apply S -1 H onto a vector. The multiplication of S -1 is performed using a preconditioned conjugate-gradient method and does not involve the explicit inversion of S. Hence the method scales the same way as the original KPM method, i.e., O(N), although there is an overhead due to the additional conjugate-gradient part. We apply this method to a large scale electronic-structure calculation of amorphous diamond. copyright 1997 The American Physical Society

Numerical method for three dimensional steady-state two-phase flow calculations

International Nuclear Information System (INIS)

Raymond, P.; Toumi, I.

1992-01-01

This paper presents the numerical scheme which was developed for the FLICA-4 computer code to calculate three dimensional steady state two phase flows. This computer code is devoted to steady state and transient thermal hydraulics analysis of nuclear reactor cores 1,3 . The first section briefly describes the FLICA-4 flow modelling. Then in order to introduce the numerical method for steady state computations, some details are given about the implicit numerical scheme based upon an approximate Riemann solver which was developed for calculation of flow transients. The third section deals with the numerical method for steady state computations, which is derived from this previous general scheme and its optimization. We give some numerical results for steady state calculations and comparisons on required CPU time and memory for various meshing and linear system solvers

A computer program for unilateral renal clearance calculation by a modified Oberhausen method

International Nuclear Information System (INIS)

Brueggemann, G.

1980-01-01

A FORTAN program is presented which, on the basis of data obtained with NUKLEOPAN M, calculates the glomerular filtration rate with sup(99m)Tc-DTPA, the unilateral effective renal plasma flow with 131 I-hippuran, and the parameters for describing the isotope rephrogram (ING) with 131 I-hippuran. The results are calculated fully automatically upon entry of the data, and the results are processed and printed out. The theoretical fundamentals of ING and whole-body clearance calculation are presented as well as the methods available for unilateral clearance calculation, and the FORTAN program is described in detail. The standard values of the method are documented, as well as a comparative gamma camera study of 48 patients in order to determine the accuracy of unilateral imaging with the NUKLEOPAN M instrument, a comparison of unilateral clearances by the Oberhausen and Taplin methods, and a comparison between 7/17' plasma clearance and whole-body clearance. Problems and findings of the method are discussed. (orig./MG) [de

Earthquake analysis of structures including structure-soil interaction by a substructure method

International Nuclear Information System (INIS)

Chopra, A.K.; Guttierrez, J.A.

1977-01-01

A general substructure method for analysis of response of nuclear power plant structures to earthquake ground motion, including the effects of structure-soil interaction, is summarized. The method is applicable to complex structures idealized as finite element systems and the soil region treated as either a continuum, for example as a viscoelastic halfspace, or idealized as a finite element system. The halfspace idealization permits reliable analysis for sites where essentially similar soils extend to large depths and there is no rigid boundary such as soil-rock interface. For sites where layers of soft soil are underlain by rock at shallow depth, finite element idealization of the soil region is appropriate; in this case, the direct and substructure methods would lead to equivalent results but the latter provides the better alternative. Treating the free field motion directly as the earthquake input in the substructure method eliminates the deconvolution calculations and the related assumption -regarding type and direction of earthquake waves- required in the direct method. The substructure method is computationally efficient because the two substructures-the structure and the soil region- are analyzed separately; and, more important, it permits taking advantage of the important feature that response to earthquake ground motion is essentially contained in the lower few natural modes of vibration of the structure on fixed base. For sites where essentially similar soils extend to large depths and there is no obvious rigid boundary such as a soil-rock interface, numerical results for earthquake response of a nuclear reactor structure are presented to demonstrate that the commonly used finite element method may lead to unacceptable errors; but the substructure method leads to reliable results

Self-consistent DFT +U method for real-space time-dependent density functional theory calculations

Science.gov (United States)

Tancogne-Dejean, Nicolas; Oliveira, Micael J. T.; Rubio, Angel

2017-12-01

We implemented various DFT+U schemes, including the Agapito, Curtarolo, and Buongiorno Nardelli functional (ACBN0) self-consistent density-functional version of the DFT +U method [Phys. Rev. X 5, 011006 (2015), 10.1103/PhysRevX.5.011006] within the massively parallel real-space time-dependent density functional theory (TDDFT) code octopus. We further extended the method to the case of the calculation of response functions with real-time TDDFT+U and to the description of noncollinear spin systems. The implementation is tested by investigating the ground-state and optical properties of various transition-metal oxides, bulk topological insulators, and molecules. Our results are found to be in good agreement with previously published results for both the electronic band structure and structural properties. The self-consistent calculated values of U and J are also in good agreement with the values commonly used in the literature. We found that the time-dependent extension of the self-consistent DFT+U method yields improved optical properties when compared to the empirical TDDFT+U scheme. This work thus opens a different theoretical framework to address the nonequilibrium properties of correlated systems.

Ca(AlH4)2, CaAlH5, and CaH2+6LiBH4 : Calculated dehydrogenation enthalpy, including zero point energy, and the structure of the phonon spectra

NARCIS (Netherlands)

Marashdeh, A.; Frankcombe, T.J.

2008-01-01

The dehydrogenation enthalpies of Ca(AlH4)2, CaAlH5, and CaH2+6LiBH4 have been calculated using density functional theory calculations at the generalized gradient approximation level. Harmonic phonon zero point energy (ZPE) corrections have been included using Parlinski’s direct method. The

Transport calculation of medium-energy protons and neutrons by Monte Carlo method

International Nuclear Information System (INIS)

Ban, Syuuichi; Hirayama, Hideo; Katoh, Kazuaki.

1978-09-01

A Monte Carlo transport code, ARIES, has been developed for protons and neutrons at medium energy (25 -- 500 MeV). Nuclear data provided by R.G. Alsmiller, Jr. were used for the calculation. To simulate the cascade development in the medium, each generation was represented by a single weighted particle and an average number of emitted particles was used as the weight. Neutron fluxes were stored by the collisions density method. The cutoff energy was set to 25 MeV. Neutrons below the cutoff were stored to be used as the source for the low energy neutron transport calculation upon the discrete ordinates method. Then transport calculations were performed for both low energy neutrons (thermal -- 25 MeV) and secondary gamma-rays. Energy spectra of emitted neutrons were calculated and compared with those of published experimental and calculated results. The agreement was good for the incident particles of energy between 100 and 500 MeV. (author)

Optimization method for quantitative calculation of clay minerals in soil

Indian Academy of Sciences (India)

However, no reliable method for quantitative analysis of clay minerals has been established so far. In this study, an attempt was made to propose an optimization method for the quantitative ... 2. Basic principles. The mineralogical constitution of soil is rather complex. ... K2O, MgO, and TFe as variables for the calculation.

Calculation of the fast multiplication factor by the fission matrix method

International Nuclear Information System (INIS)

Naumov, V.A.; Rozin, S.G.; Ehl'perin, T.I.

1976-01-01

A variation of the Monte Carlo method to calculate an effective breeding factor of a nuclear reactor is described. The evaluation procedure of reactivity perturbations by the Monte Carlo method in the first order perturbation theory is considered. The method consists in reducing an integral neutron transport equation to a set of linear algebraic equations. The coefficients of this set are elements of a fission matrix. The fission matrix being a Grin function of the neutron transport equation, is evaluated by the Monte Carlo method. In the program realizing the suggested algorithm, the game for initial neutron energy of a fission spectrum and then for the region of neutron birth, ΔVsub(f)sup(i)has been played in proportion to the product of Σsub(f)sup(i)ΔVsub(f)sup(i), where Σsub(f)sup(i) is a macroscopic cross section in the region numbered at the birth energy. Further iterations of a space distribution of neutrons in the system are performed by the generation method. In the adopted scheme of simulation of neutron histories the emission of secondary neutrons is controlled by weights; it occurs at every collision and not only in the end on the history. The breeding factor is calculated simultaneously with the space distribution of neutron worth in the system relative to the fission process and neutron flux. Efficiency of the described procedure has been tested on the calculation of the breeding factor for the Godiva assembly, simulating a fast reactor with a hard spectrum. A high accuracy of calculations at moderate number of zones in the core and reasonable statistics has been stated

A new method for the automatic calculation of prosody

International Nuclear Information System (INIS)

GUIDINI, Annie

1981-01-01

An algorithm is presented for the calculation of the prosodic parameters for speech synthesis. It uses the melodic patterns, composed of rising and falling slopes, suggested by G. CAELEN, and rests on: 1. An analysis into units of meaning to determine a melodic pattern 2. the calculation of the numeric values for the prosodic variations of each syllable; 3. The use of a table of vocalic values for the three parameters for each vowel according to the consonantal environment and of a table of standard duration for consonants. This method was applied in the 'SARA' program of synthesis with satisfactory results. (author) [fr

4He binding energy calculation including full tensor-force effects

Science.gov (United States)

Fonseca, A. C.

1989-09-01

The four-body equations of Alt, Grassberger, and Sandhas are solved in the version where the (2)+(2) subamplitudes are treated exactly by convolution, using one-term separable Yamaguchy nucleon-nucleon potentials in the 1S0 and 3S1-3D1 channels. The resulting jp=1/2+ and (3/2+ three-body subamplitudes are represented in a separable form using the energy-dependent pole expansion. Converged bound-state results are calculated for the first time using the full interaction, and are compared with those obtained from a simplified treatment of the tensor force. The Tjon line that correlates three-nucleon and four-nucleon binding energies is shown using different nucleon-nucleon potentials. In all calculations the Coulomb force has been neglected.

Calculations of pair production by Monte Carlo methods

International Nuclear Information System (INIS)

Bottcher, C.; Strayer, M.R.

1991-01-01

We describe some of the technical design issues associated with the production of particle-antiparticle pairs in very large accelerators. To answer these questions requires extensive calculation of Feynman diagrams, in effect multi-dimensional integrals, which we evaluate by Monte Carlo methods on a variety of supercomputers. We present some portable algorithms for generating random numbers on vector and parallel architecture machines. 12 refs., 14 figs

Bulk Electric Load Cost Calculation Methods: Iraqi Network Comparative Study

Directory of Open Access Journals (Sweden)

Qais M. Alias

2016-09-01

Full Text Available It is vital in any industry to regain the spent capitals plus running costs and a margin of profits for the industry to flourish. The electricity industry is an everyday life touching industry which follows the same finance-economic strategy. Cost allocation is a major issue in all sectors of the electric industry, viz, generation, transmission and distribution. Generation and distribution service costing’s well documented in the literature, while the transmission share is still of need for research. In this work, the cost of supplying a bulk electric load connected to the EHV system is calculated. A sample basic lump-average method is used to provide a rough costing guide. Also, two transmission pricing methods are employed, namely, the postage-stamp and the load-flow based MW-distance methods to calculate transmission share in the total cost of each individual bulk load. The three costing methods results are then analyzed and compared for the 400kV Iraqi power grid considered for a case study.

Monte-Carlo Method Python Library for dose distribution Calculation in Brachytherapy

Energy Technology Data Exchange (ETDEWEB)

Randriantsizafy, R D; Ramanandraibe, M J [Madagascar Institut National des Sciences et Techniques Nucleaires, Antananarivo (Madagascar); Raboanary, R [Institut of astro and High-Energy Physics Madagascar, University of Antananarivo, Antananarivo (Madagascar)

2007-07-01

The Cs-137 Brachytherapy treatment is performed in Madagascar since 2005. Time treatment calculation for prescribed dose is made manually. Monte-Carlo Method Python library written at Madagascar INSTN is experimentally used to calculate the dose distribution on the tumour and around it. The first validation of the code was done by comparing the library curves with the Nucletron company curves. To reduce the duration of the calculation, a Grid of PC's is set up with listner patch run on each PC. The library will be used to modelize the dose distribution in the CT scan patient picture for individual and better accuracy time calculation for a prescribed dose.

Monte-Carlo Method Python Library for dose distribution Calculation in Brachytherapy

International Nuclear Information System (INIS)

Randriantsizafy, R.D.; Ramanandraibe, M.J.; Raboanary, R.

2007-01-01

The Cs-137 Brachytherapy treatment is performed in Madagascar since 2005. Time treatment calculation for prescribed dose is made manually. Monte-Carlo Method Python library written at Madagascar INSTN is experimentally used to calculate the dose distribution on the tumour and around it. The first validation of the code was done by comparing the library curves with the Nucletron company curves. To reduce the duration of the calculation, a Grid of PC's is set up with listner patch run on each PC. The library will be used to modelize the dose distribution in the CT scan patient picture for individual and better accuracy time calculation for a prescribed dose.

HIRFL-SSC trim coil currents calculation by conjugate gradients method

International Nuclear Information System (INIS)

Liu, W.

2005-01-01

For accelerating different kinds of ions to various energies, the HIRFL-SSC should form the corresponding isochronous magnetic field by its main coil and trim coils. Previously, there were errors in fitting the theoretical isochronous magnetic field in the small radius region, which led to some operation difficulties for ion acceleration in the inject region. After further investigation of the restrictive condition of the maximum current limitation, the trim coil currents for fitting the theoretical isochronous magnetic field were recalculated by the conjugate gradients method. Better results were obtained in the operation of HIRFL-SSC. This article introduces the procedure to calculate the trim coil currents. The calculation method of conjugate gradients is introduced and the fitting error is analysed. (author)

A modified method of calculating the lateral build-up ratio for small electron fields

International Nuclear Information System (INIS)

Tyner, E; McCavana, P; McClean, B

2006-01-01

This note outlines an improved method of calculating dose per monitor unit values for small electron fields using Khan's lateral build-up ratio (LBR). This modified method obtains the LBR directly from the ratio of measured, surface normalized, electron beam percentage depth dose curves. The LBR calculated using this modified method more accurately accounts for the change in lateral scatter with decreasing field size. The LBR is used along with Khan's dose per monitor unit formula to calculate dose per monitor unit values for a set of small fields. These calculated dose per monitor unit values are compared to measured values to within 3.5% for all circular fields and electron energies examined. The modified method was further tested using a small triangular field. A maximum difference of 4.8% was found. (note)

Calculation of Pressure Distribution at Rotary Body Surface with the Vortex Element Method

Directory of Open Access Journals (Sweden)

S. A. Dergachev

2014-01-01

Full Text Available Vortex element method allows to simulate unsteady hydrodynamic processes in incompressible environment, taking into account the evolution of the vortex sheet, including taking into account the deformation or moving of the body or part of construction.For the calculation of the hydrodynamic characteristics of the method based on vortex element software package was developed MVE3D. Vortex element (VE in program is symmetrical Vorton-cut. For satisfying the boundary conditions at the surface used closed frame of vortons.With this software system modeled incompressible flow around a cylindrical body protection elongation L / D = 13 with a front spherical blunt with the angle of attack of 10 °. We analyzed the distribution of the pressure coefficient on the body surface of the top and bottom forming.The calculate results were compared with known Results of experiment.Considered design schemes with different number of Vorton framework. Also varied radius of VE. Calculation make possible to establish the degree of sampling surface needed to produce close to experiment results. It has been shown that an adequate reproducing the pressure distribution in the transition region spherical cylindrical surface, on the windward side requires a high degree of sampling.Based on these results Can be possible need to improve on the design scheme of body's surface, allowing more accurate to describe the flow vorticity in areas with abrupt changes of geometry streamlined body.

CALCULATION METHODS OF OPTIMAL ADJUSTMENT OF CONTROL SYSTEM THROUGH DISTURBANCE CHANNEL

Directory of Open Access Journals (Sweden)

I. M. Golinko

2014-01-01

Full Text Available In the process of automatic control system debugging the great attention is paid to determining formulas’ parameters of optimal dynamic adjustment of regulators, taking into account the dynamics of Objects control. In most cases the known formulas are oriented on design of automatic control system through channel “input-output definition”. But practically in all continuous processes the main task of all regulators is stabilization of output parameters. The Methods of parameters calculation for dynamic adjustment of regulations were developed. These methods allow to optimize the analog and digital regulators, taking into account minimization of regulated influences. There were suggested to use the fact of detuning and maximum value of regulated influence. As the automatic control system optimization with proportional plus reset controllers on disturbance channel is an unimodal task, the main algorithm of optimization is realized by Hooke – Jeeves method. For controllers optimization through channel external disturbance there were obtained functional dependences of parameters calculations of dynamic proportional plus reset controllers from dynamic characteristics of Object control. The obtained dependences allow to improve the work of controllers (regulators of automatic control on external disturbance channel and so it allows to improve the quality of regulation of transient processes. Calculation formulas provide high accuracy and convenience in usage. In suggested method there are no nomographs and this fact expels subjectivity of investigation in determination of parameters of dynamic adjustment of proportional plus reset controllers. Functional dependences can be used for calculation of adjustment of PR controllers in a great range of change of dynamic characteristics of Objects control.

Comparison of the methods for calculating the interfacial heat transfer coefficient in hot stamping

International Nuclear Information System (INIS)

Zhao, Kunmin; Wang, Bin; Chang, Ying; Tang, Xinghui; Yan, Jianwen

2015-01-01

This paper presents a hot stamping experimentation and three methods for calculating the Interfacial Heat Transfer Coefficient (IHTC) of 22MnB5 boron steel. Comparison of the calculation results shows an average error of 7.5% for the heat balance method, 3.7% for the Beck's nonlinear inverse estimation method (the Beck's method), and 10.3% for the finite-element-analysis-based optimization method (the FEA method). The Beck's method is a robust and accurate method for identifying the IHTC in hot stamping applications. The numerical simulation using the IHTC identified by the Beck's method can predict the temperature field with a high accuracy. - Highlights: • A theoretical formula was derived for direct calculation of IHTC. • The Beck's method is a robust and accurate method for identifying IHTC. • Finite element method can be used to identify an overall equivalent IHTC

Applying Activity Based Costing (ABC) Method to Calculate Cost Price in Hospital and Remedy Services.

Science.gov (United States)

Rajabi, A; Dabiri, A

2012-01-01

Activity Based Costing (ABC) is one of the new methods began appearing as a costing methodology in the 1990's. It calculates cost price by determining the usage of resources. In this study, ABC method was used for calculating cost price of remedial services in hospitals. To apply ABC method, Shahid Faghihi Hospital was selected. First, hospital units were divided into three main departments: administrative, diagnostic, and hospitalized. Second, activity centers were defined by the activity analysis method. Third, costs of administrative activity centers were allocated into diagnostic and operational departments based on the cost driver. Finally, with regard to the usage of cost objectives from services of activity centers, the cost price of medical services was calculated. The cost price from ABC method significantly differs from tariff method. In addition, high amount of indirect costs in the hospital indicates that capacities of resources are not used properly. Cost price of remedial services with tariff method is not properly calculated when compared with ABC method. ABC calculates cost price by applying suitable mechanisms but tariff method is based on the fixed price. In addition, ABC represents useful information about the amount and combination of cost price services.

Use of results from microscopic methods in optical model calculations

International Nuclear Information System (INIS)

Lagrange, C.

1985-11-01

A concept of vectorization for coupled-channel programs based upon conventional methods is first presented. This has been implanted in our program for its use on the CRAY-1 computer. In a second part we investigate the capabilities of a semi-microscopic optical model involving fewer adjustable parameters than phenomenological ones. The two main ingredients of our calculations are, for spherical or well-deformed nuclei, the microscopic optical-model calculations of Jeukenne, Lejeune and Mahaux and nuclear densities from Hartree-Fock-Bogoliubov calculations using the density-dependent force D1. For transitional nuclei deformation-dependent nuclear structure wave functions are employed to weigh the scattering potentials for different shapes and channels [fr

Effectiveness of the current method of calculating member states' contributions

CERN Document Server

2002-01-01

At its Two-hundred and eighty-sixth Meeting of 19 September 2001, the Finance Committee requested the Management to re-assess the effectiveness of the current method of forecasting Net National Income (NNI) for the purposes of calculating the Member States' contributions by comparing the results of the current weighted average method with a method based on a simple arithmetic average. The Finance Committee is invited to take note of this information.

Advances in computational methods for Quantum Field Theory calculations

NARCIS (Netherlands)

Ruijl, B.J.G.

2017-01-01

In this work we describe three methods to improve the performance of Quantum Field Theory calculations. First, we simplify large expressions to speed up numerical integrations. Second, we design Forcer, a program for the reduction of four-loop massless propagator integrals. Third, we extend the R*

Comparison of optimization methods for electronic-structure calculations

International Nuclear Information System (INIS)

Garner, J.; Das, S.G.; Min, B.I.; Woodward, C.; Benedek, R.

1989-01-01

The performance of several local-optimization methods for calculating electronic structure is compared. The fictitious first-order equation of motion proposed by Williams and Soler is integrated numerically by three procedures: simple finite-difference integration, approximate analytical integration (the Williams-Soler algorithm), and the Born perturbation series. These techniques are applied to a model problem for which exact solutions are known, the Mathieu equation. The Williams-Soler algorithm and the second Born approximation converge equally rapidly, but the former involves considerably less computational effort and gives a more accurate converged solution. Application of the method of conjugate gradients to the Mathieu equation is discussed

Nuclear calculation methods for light water moderated reactors

International Nuclear Information System (INIS)

Hicks, D.

1961-02-01

This report is intended as an introductory review. After a brief discussion of problems encountered in the nuclear design of water moderated reactors a comprehensive scheme of calculations is described. This scheme is based largely on theoretical methods and computer codes developed in the U.S.A. but some previously unreported developments made in this country are also described. It is shown that the effective reproduction factor of simple water moderated lattices may be estimated to an accuracy of approximately 1%. Methods for treating water gap flux peaking and control absorbers are presented in some detail, together with a brief discussion of temperature coefficients, void coefficients and burn-up problems. (author)

Method of the reduced-added Green function in the calculation of atomic polarizabilities

International Nuclear Information System (INIS)

Chernov, V.E.; Dorofeev, D.L.; Kretinin, I.Yu.; Zon, B.A.

2005-01-01

The Green function in the quantum defect theory provides an exact account for high-excited and continuum electronic states. We modify it by taking into account the ground and low-excited states using their wave functions calculated ab initio. As an application, we present a simple and efficient semianalytical method for the calculation of atomic electric frequency-dependent scalar dipole polarizability, for both real and imaginary frequencies. The polarizabilities calculated for some atoms (Li, Na, K, Be, Mg, Ca, Si, P, S, O, Al, Ge, C, N, F, He, Ne, Ar, Kr, and Xe) are compared with existing methods of computational quantum chemistry and with experiments; good accuracy of the proposed method is demonstrated

A method of calculation on the airloading of vertical axis wind turbine

Science.gov (United States)

Azuma, A.; Kimura, S.

A new method of analyzing the aerodynamic characteristics of the Darrieus Vertical-Axis Wind Turbine (VAWT) by applying the local circulation method is described. The validity of this method is confirmed by analyzing the air load acting on a curved blade. The azimuthwise variation of spanwise airloading, torque, and longitudinal forces are accurately calculated for a variety of operational conditions. The results are found to be in good agreement with experimental ones obtained elsewhere. It is concluded that the present approach can calculate the aerodynamic characteristics of the VAWT with much less computational time than that used by the free vortex model.

A method for calculation of dose per unit concentration values for aquatic biota

International Nuclear Information System (INIS)

Batlle, J Vives i; Jones, S R; Gomez-Ros, J M

2004-01-01

A dose per unit concentration database has been generated for application to ecosystem assessments within the FASSET framework. Organisms are represented by ellipsoids of appropriate dimensions, and the proportion of radiation absorbed within the organisms is calculated using a numerical method implemented in a series of spreadsheet-based programs. Energy-dependent absorbed fraction functions have been derived for calculating the total dose per unit concentration of radionuclides present in biota or in the media they inhabit. All radionuclides and reference organism dimensions defined within FASSET for marine and freshwater ecosystems are included. The methodology has been validated against more complex dosimetric models and compared with human dosimetry based on ICRP 72. Ecosystem assessments for aquatic biota within the FASSET framework can now be performed simply, once radionuclide concentrations in target organisms are known, either directly or indirectly by deduction from radionuclide concentrations in the surrounding medium

The calculation of viscosity of liquid n-decane and n-hexadecane by the Green-Kubo method

Science.gov (United States)

Cui, S. T.; Cummings, P. T.; Cochran, H. D.

This short commentary presents the result of long molecular dynamics simulation calculations of the shear viscosity of liquid n-decane and n-hexadecane using the Green-Kubo integration method. The relaxation time of the stress-stress correlation function is compared with those of rotation and diffusion. The rotational and diffusional relaxation times, which are easy to calculate, provide useful guides for the required simulation time in viscosity calculations. Also, the computational time required for viscosity calculations of these systems by the Green-Kubo method is compared with the time required for previous non-equilibrium molecular dynamics calculations of the same systems. The method of choice for a particular calculation is determined largely by the properties of interest, since the efficiencies of the two methods are comparable for calculation of the zero strain rate viscosity.

Nested element method in multidimensional neutron diffusion calculations

International Nuclear Information System (INIS)

Altiparmakov, D.V.

1983-01-01

A new numerical method is developed that is particularly efficient in solving the multidimensional neutron diffusion equation in geometrically complex systems. The needs for a generally applicable and fast running computer code have stimulated the inroad of a nonclassical (R-function) numerical method into the nuclear field. By using the R-functions, the geometrical components of the diffusion problem are a priori analytically implemented into the approximate solution. The class of functions, to which the approximate solution belongs, is chosen as close to the exact solution class as practically acceptable from the time consumption point of view. That implies a drastic reduction of the number of degrees of freedom, compared to the other methods. Furthermore, the reduced number of degrees of freedom enables calculation of large multidimensional problems on small computers

New nonlinear methods for linear transport calculations

International Nuclear Information System (INIS)

Adams, M.L.

1993-01-01

We present a new family of methods for the numerical solution of the linear transport equation. With these methods an iteration consists of an 'S N sweep' followed by an 'S 2 -like' calculation. We show, by analysis as well as numerical results, that iterative convergence is always rapid. We show that this rapid convergence does not depend on a consistent discretization of the S 2 -like equations - they can be discretized independently from the S N equations. We show further that independent discretizations can offer significant advantages over consistent ones. In particular, we find that in a wide range of problems, an accurate discretization of the S 2 -like equation can be combined with a crude discretization of the S N equations to produce an accurate S N answer. We demonstrate this by analysis as well as numerical results. (orig.)

The application of advanced rotor (performance) methods for design calculations

Energy Technology Data Exchange (ETDEWEB)

Bussel, G.J.W. van [Delft Univ. of Technology, Inst. for Wind Energy, Delft (Netherlands)

1997-08-01

The calculation of loads and performance of wind turbine rotors has been a topic for research over the last century. The principles for the calculation of loads on rotor blades with a given specific geometry, as well as the development of optimal shaped rotor blades have been published in the decades that significant aircraft development took place. Nowadays advanced computer codes are used for specific problems regarding modern aircraft, and application to wind turbine rotors has also been performed occasionally. The engineers designing rotor blades for wind turbines still use methods based upon global principles developed in the beginning of the century. The question what to expect in terms of the type of methods to be applied in a design environment for the near future is addressed here. (EG) 14 refs.

Comparing Methods of Calculating Expected Annual Damage in Urban Pluvial Flood Risk Assessments

DEFF Research Database (Denmark)

Skovgård Olsen, Anders; Zhou, Qianqian; Linde, Jens Jørgen

2015-01-01

Estimating the expected annual damage (EAD) due to flooding in an urban area is of great interest for urban water managers and other stakeholders. It is a strong indicator for a given area showing how vulnerable it is to flood risk and how much can be gained by implementing e.g., climate change...... adaptation measures. This study identifies and compares three different methods for estimating the EAD based on unit costs of flooding of urban assets. One of these methods was used in previous studies and calculates the EAD based on a few extreme events by assuming a log-linear relationship between cost...... of an event and the corresponding return period. This method is compared to methods that are either more complicated or require more calculations. The choice of method by which the EAD is calculated appears to be of minor importance. At all three case study areas it seems more important that there is a shift...

Methods of producing adsorption media including a metal oxide

Science.gov (United States)

Mann, Nicholas R; Tranter, Troy J

2014-03-04

Methods of producing a metal oxide are disclosed. The method comprises dissolving a metal salt in a reaction solvent to form a metal salt/reaction solvent solution. The metal salt is converted to a metal oxide and a caustic solution is added to the metal oxide/reaction solvent solution to adjust the pH of the metal oxide/reaction solvent solution to less than approximately 7.0. The metal oxide is precipitated and recovered. A method of producing adsorption media including the metal oxide is also disclosed, as is a precursor of an active component including particles of a metal oxide.

Simplified method of calculating residual stress in circumferential welding of piping

International Nuclear Information System (INIS)

Umemoto, Tadahiro

1984-01-01

Many circumferential joints of piping are used in as-welded state, but in these welded joints, the residual stress as high as the yield stress of materials arises, and causes to accelerate stress corrosion cracking and corrosion fatigue. The experiment or the finite element method to clarify welding residual stress requires much time and labor, and is expensive, therefore, the author proposed the simplified method of calculation. The heating and cooling process of welding is very complex, and cannot be modeled as it is, therefore, it was assumed that in multiple layer welding, the welding condition of the last layer determines the residual stress, that material constants are invariable regardless of temperature, that the temperature distribution and residual stress are axisymmetric, and that there is repeated stress-strain relation in the vicinity of welded parts. The temperature distribution at the time of welding, thermal stress and welding residual stress are analyzed, and the material constants used for the calculation of residual stress are given. As the example of calculation, the effect of welding heat input and materials is shown. The extension of the method to a thick-walled pipe is discussed. (Kako, I.)

Detent Force Calculations of a PMLSM Using the Finite Element Method

Science.gov (United States)

Remy, Ghislain; Krebs, Guillaume; Tounzi, Abdelmounaïm; Barre, Pierre-Jean

This paper presents a Finite Element Analysis of a Permanent Magnet Linear Synchronous Motor. The aim is to obtain an accurate estimation of the detent force without oversize computation. First, some usual techniques dedicated to the calculation of the forces in electromagnetic devices, such as the Virtual Work Method and the Maxwell Stress Tensor, are described. Some keypoints of the meshing method using a commercial FEM software are presented and used in order to improve the thrust computations. After that, the topology and features of the studied motor are described to highlight specific problems of the modelling process. In the 2D FEM case, new meshing techniques are proposed, according to the force calculations. The FEM results obtained from the different methods are analysed and compared with the experimental ones. Second, using FEM results, a study of the independence of the cogging and the end-effect forces is presented. Particularly, an original approach is suggested in order to compute the cogging force only, using the same mesh for each motion step. Then, the PMLSM geometry is adapted to calculate the end-effect forces only.

Methodes de calcul des forces aerodynamiques pour les etudes des interactions aeroservoelastiques

Science.gov (United States)

Biskri, Djallel Eddine

L'aeroservoelasticite est un domaine ou interagissent la structure flexible d'un avion, l'aerodynamique et la commande de vol. De son cote, la commande du vol considere l'avion comme une structure rigide et etudie l'influence du systeme de commande sur la dynamique de vol. Dans cette these, nous avons code trois nouvelles methodes d'approximation de forces aerodynamiques: Moindres carres corriges, Etat minimal corrige et Etats combines. Dans les deux premieres methodes, les erreurs d'approximation entre les forces aerodynamiques approximees par les methodes classiques et celles obtenues par les nouvelles methodes ont les memes formes analytiques que celles des forces aerodynamiques calculees par LS ou MS. Quant a la troisieme methode, celle-ci combine les formulations des forces approximees avec les methodes standards LS et MS. Les vitesses et frequences de battement et les temps d'executions calcules par les nouvelles methodes versus ceux calcules par les methodes classiques ont ete analyses.

Efficient parallel implicit methods for rotary-wing aerodynamics calculations

Science.gov (United States)

Wissink, Andrew M.

Euler/Navier-Stokes Computational Fluid Dynamics (CFD) methods are commonly used for prediction of the aerodynamics and aeroacoustics of modern rotary-wing aircraft. However, their widespread application to large complex problems is limited lack of adequate computing power. Parallel processing offers the potential for dramatic increases in computing power, but most conventional implicit solution methods are inefficient in parallel and new techniques must be adopted to realize its potential. This work proposes alternative implicit schemes for Euler/Navier-Stokes rotary-wing calculations which are robust and efficient in parallel. The first part of this work proposes an efficient parallelizable modification of the Lower Upper-Symmetric Gauss Seidel (LU-SGS) implicit operator used in the well-known Transonic Unsteady Rotor Navier Stokes (TURNS) code. The new hybrid LU-SGS scheme couples a point-relaxation approach of the Data Parallel-Lower Upper Relaxation (DP-LUR) algorithm for inter-processor communication with the Symmetric Gauss Seidel algorithm of LU-SGS for on-processor computations. With the modified operator, TURNS is implemented in parallel using Message Passing Interface (MPI) for communication. Numerical performance and parallel efficiency are evaluated on the IBM SP2 and Thinking Machines CM-5 multi-processors for a variety of steady-state and unsteady test cases. The hybrid LU-SGS scheme maintains the numerical performance of the original LU-SGS algorithm in all cases and shows a good degree of parallel efficiency. It experiences a higher degree of robustness than DP-LUR for third-order upwind solutions. The second part of this work examines use of Krylov subspace iterative solvers for the nonlinear CFD solutions. The hybrid LU-SGS scheme is used as a parallelizable preconditioner. Two iterative methods are tested, Generalized Minimum Residual (GMRES) and Orthogonal s-Step Generalized Conjugate Residual (OSGCR). The Newton method demonstrates good

Methods and codes for neutronic calculations of the MARIA research reactor

International Nuclear Information System (INIS)

Andrzejewski, K.; Kulikowska, T.; Bretscher, M.M.; Hanan, N.A.; Matos, J.E.

1998-01-01

The core of the MARIA high flux multipurpose research reactor is highly heterogeneous. It consists of beryllium blocks arranged in 6x8 matrix, tubular fuel assemblies, control rods and irradiation channels. The reflector is also heterogeneous and consists of graphite blocks clad with aluminium. Its structure is perturbed by the experimental beam tubes. This paper presents methods and codes used to calculate the MARIA reactor neutronics characteristics and experience gained thus far at IAE and ANL. At ANL the methods of MARIA calculations were developed in connection with RERTR program. At IAE the package of programs was developed to help its operator in optimization of fuel utilization. (author)

Methods of calculating the post-closure performance of high-level waste repositories

Energy Technology Data Exchange (ETDEWEB)

Ross, B. (ed.)

1989-02-01

This report is intended as an overview of post-closure performance assessment methods for high-level radioactive waste repositories and is designed to give the reader a broad sense of the state of the art of this technology. As described here, ''the state of the art'' includes only what has been reported in report, journal, and conference proceedings literature through August 1987. There is a very large literature on the performance of high-level waste repositories. In order to make a review of this breadth manageable, its scope must be carefully defined. The essential principle followed is that only methods of calculating the long-term performance of waste repositories are described. The report is organized to reflect, in a generalized way, the logical order to steps that would be taken in a typical performance assessment. Chapter 2 describes ways of identifying scenarios and estimating their probabilities. Chapter 3 presents models used to determine the physical and chemical environment of a repository, including models of heat transfer, radiation, geochemistry, rock mechanics, brine migration, radiation effects on chemistry, and coupled processes. The next two chapters address the performance of specific barriers to release of radioactivity. Chapter 4 treats engineered barriers, including containers, waste forms, backfills around waste packages, shaft and borehole seals, and repository design features. Chapter 5 discusses natural barriers, including ground water systems and stability of salt formations. The final chapters address optics of general applicability to performance assessment models. Methods of sensitivity and uncertainty analysis are described in Chapter 6, and natural analogues of repositories are treated in Chapter 7. 473 refs., 19 figs., 2 tabs.

Methods of calculating the post-closure performance of high-level waste repositories

International Nuclear Information System (INIS)

Ross, B.

1989-02-01

This report is intended as an overview of post-closure performance assessment methods for high-level radioactive waste repositories and is designed to give the reader a broad sense of the state of the art of this technology. As described here, ''the state of the art'' includes only what has been reported in report, journal, and conference proceedings literature through August 1987. There is a very large literature on the performance of high-level waste repositories. In order to make a review of this breadth manageable, its scope must be carefully defined. The essential principle followed is that only methods of calculating the long-term performance of waste repositories are described. The report is organized to reflect, in a generalized way, the logical order to steps that would be taken in a typical performance assessment. Chapter 2 describes ways of identifying scenarios and estimating their probabilities. Chapter 3 presents models used to determine the physical and chemical environment of a repository, including models of heat transfer, radiation, geochemistry, rock mechanics, brine migration, radiation effects on chemistry, and coupled processes. The next two chapters address the performance of specific barriers to release of radioactivity. Chapter 4 treats engineered barriers, including containers, waste forms, backfills around waste packages, shaft and borehole seals, and repository design features. Chapter 5 discusses natural barriers, including ground water systems and stability of salt formations. The final chapters address optics of general applicability to performance assessment models. Methods of sensitivity and uncertainty analysis are described in Chapter 6, and natural analogues of repositories are treated in Chapter 7. 473 refs., 19 figs., 2 tabs

A keff calculation method by Monte Carlo

International Nuclear Information System (INIS)

Shen, H; Wang, K.

2008-01-01

The effective multiplication factor (k eff ) is defined as the ratio between the number of neutrons in successive generations, which definition is adopted by most Monte Carlo codes (e.g. MCNP). Also, it can be thought of as the ratio of the generation rate of neutrons by the sum of the leakage rate and the absorption rate, which should exclude the effect of the neutron reaction such as (n, 2n) and (n, 3n). This article discusses the Monte Carlo method for k eff calculation based on the second definition. A new code has been developed and the results are presented. (author)

Study on application of green's function method in thermal stress rapid calculation

International Nuclear Information System (INIS)

Zhang Guihe; Duan Yuangang; Xu Xiao; Chen Rong

2013-01-01

This paper presents a quick and accuracy thermal stress calculation method, the Green's Function Method, which is a combination of finite element method and numerical algorithm method. Thermal stress calculation of Safe Injection Nozzle of Reactor Coolant Line of PWR plant is performed with Green's function method for heatup and cooldown thermal transients as a demonstration example, and the result is compared with finite element method to verify the rationality and accuracy of this method. The advantage and disadvantage of the Green's function method and the finite element method are also compared. (authors)

Neutron flux calculation by means of Monte Carlo methods

International Nuclear Information System (INIS)

Barz, H.U.; Eichhorn, M.

1988-01-01

In this report a survey of modern neutron flux calculation procedures by means of Monte Carlo methods is given. Due to the progress in the development of variance reduction techniques and the improvements of computational techniques this method is of increasing importance. The basic ideas in application of Monte Carlo methods are briefly outlined. In more detail various possibilities of non-analog games and estimation procedures are presented, problems in the field of optimizing the variance reduction techniques are discussed. In the last part some important international Monte Carlo codes and own codes of the authors are listed and special applications are described. (author)

An Efficient numerical method to calculate the conductivity tensor for disordered topological matter

Science.gov (United States)

Garcia, Jose H.; Covaci, Lucian; Rappoport, Tatiana G.

2015-03-01

We propose a new efficient numerical approach to calculate the conductivity tensor in solids. We use a real-space implementation of the Kubo formalism where both diagonal and off-diagonal conductivities are treated in the same footing. We adopt a formulation of the Kubo theory that is known as Bastin formula and expand the Green's functions involved in terms of Chebyshev polynomials using the kernel polynomial method. Within this method, all the computational effort is on the calculation of the expansion coefficients. It also has the advantage of obtaining both conductivities in a single calculation step and for various values of temperature and chemical potential, capturing the topology of the band-structure. Our numerical technique is very general and is suitable for the calculation of transport properties of disordered systems. We analyze how the method's accuracy varies with the number of moments used in the expansion and illustrate our approach by calculating the transverse conductivity of different topological systems. T.G.R, J.H.G and L.C. acknowledge Brazilian agencies CNPq, FAPERJ and INCT de Nanoestruturas de Carbono, Flemish Science Foundation for financial support.

Memory-efficient calculations of adjoint-weighted tallies by the Monte Carlo Wielandt method

International Nuclear Information System (INIS)

Choi, Sung Hoon; Shim, Hyung Jin

2016-01-01

Highlights: • The MC Wielandt method is applied to reduce memory for the adjoint estimation. • The adjoint-weighted kinetics parameters are estimated in the MC Wielandt calculations. • The MC S/U analyses are conducted in the MC Wielandt calculations. - Abstract: The current Monte Carlo (MC) adjoint-weighted tally techniques based on the iterated fission probability (IFP) concept require a memory amount which is proportional to the numbers of the adjoint-weighted tallies and histories per cycle to store history-wise tally estimates during the convergence of the adjoint flux. Especially the conventional MC adjoint-weighted perturbation (AWP) calculations for the nuclear data sensitivity and uncertainty (S/U) analysis suffer from the huge memory consumption to realize the IFP concept. In order to reduce the memory requirement drastically, we present a new adjoint estimation method in which the memory usage is irrelevant to the numbers of histories per cycle by applying the IFP concept for the MC Wielandt calculations. The new algorithms for the adjoint-weighted kinetics parameter estimation and the AWP calculations in the MC Wielandt method are implemented in a Seoul National University MC code, McCARD and its validity is demonstrated in critical facility problems. From the comparison of the nuclear data S/U analyses, it is demonstrated that the memory amounts to store the sensitivity estimates in the proposed method become negligibly small.

An efficient method for hybrid density functional calculation with spin-orbit coupling

Science.gov (United States)

Wang, Maoyuan; Liu, Gui-Bin; Guo, Hong; Yao, Yugui

2018-03-01

In first-principles calculations, hybrid functional is often used to improve accuracy from local exchange correlation functionals. A drawback is that evaluating the hybrid functional needs significantly more computing effort. When spin-orbit coupling (SOC) is taken into account, the non-collinear spin structure increases computing effort by at least eight times. As a result, hybrid functional calculations with SOC are intractable in most cases. In this paper, we present an approximate solution to this problem by developing an efficient method based on a mixed linear combination of atomic orbital (LCAO) scheme. We demonstrate the power of this method using several examples and we show that the results compare very well with those of direct hybrid functional calculations with SOC, yet the method only requires a computing effort similar to that without SOC. The presented technique provides a good balance between computing efficiency and accuracy, and it can be extended to magnetic materials.

Study on calculation of influence function with fracture mechanics analysis method for circumferential through-crack pipe

International Nuclear Information System (INIS)

Zheng Bin; Lu Yuechuan; Zang Fenggang; Sun Yingxue

2009-01-01

In order to widen the application of the engineering method of EPRI, with a series of analysis on the 3D elastic and elastic-plastic fracture mechanics finite element, the crack open displacements (COD) of cracked pipe were calculated and a key influence function h 2 in EPRI engineering method was studied against the COD results of FEM. A calculation method of h2 under the condition of tension and bending combined load was introduced in detail. In order to validate this method, the calculated h 2 results were compared with that of EPRI, and the calculated COD results based on the h 2 results were compared with that of PICEP. The compared results indicated that the calculated h 2 results as well as the COD results and the corresponding reference values were respectively accordant, and the calculation method in this paper was validated accordingly. (authors)

Calculation of passive earth pressure of cohesive soil based on Culmann's method

Directory of Open Access Journals (Sweden)

Hai-feng Lu

2011-03-01

Full Text Available Based on the sliding plane hypothesis of Coulumb earth pressure theory, a new method for calculation of the passive earth pressure of cohesive soil was constructed with Culmann's graphical construction. The influences of the cohesive force, adhesive force, and the fill surface form were considered in this method. In order to obtain the passive earth pressure and sliding plane angle, a program based on the sliding surface assumption was developed with the VB.NET programming language. The calculated results from this method were basically the same as those from the Rankine theory and Coulumb theory formulas. This method is conceptually clear, and the corresponding formulas given in this paper are simple and convenient for application when the fill surface form is complex.

TORT/MCNP coupling method for the calculation of neutron flux around a core of BWR

International Nuclear Information System (INIS)

Kurosawa, M.

2005-01-01

For the analysis of BWR neutronics performance, accurate data are required for neutron flux distribution over the In-Reactor Pressure Vessel equipments taking into account the detailed geometrical arrangement. The TORT code can calculate neutron flux around a core of BWR in a three-dimensional geometry model, but has difficulties in fine geometrical modelling and lacks huge computer resource. On the other hand, the MCNP code enables the calculation of the neutron flux with a detailed geometry model, but requires very long sampling time to give enough number of particles. Therefore, a TORT/MCNP coupling method has been developed to eliminate the two problems mentioned above in each code. In this method, the TORT code calculates angular flux distribution on the core surface and the MCNP code calculates neutron spectrum at the points of interest using the flux distribution. The coupling method will be used as the DOT-DOMINO-MORSE code system. This TORT/MCNP coupling method was applied to calculate the neutron flux at points where induced radioactivity data were measured for 54 Mn and 60 Co and the radioactivity calculations based on the neutron flux obtained from the above method were compared with the measured data. (authors)

TORT/MCNP coupling method for the calculation of neutron flux around a core of BWR.

Science.gov (United States)

Kurosawa, Masahiko

2005-01-01

For the analysis of BWR neutronics performance, accurate data are required for neutron flux distribution over the In-Reactor Pressure Vessel equipments taking into account the detailed geometrical arrangement. The TORT code can calculate neutron flux around a core of BWR in a three-dimensional geometry model, but has difficulties in fine geometrical modelling and lacks huge computer resource. On the other hand, the MCNP code enables the calculation of the neutron flux with a detailed geometry model, but requires very long sampling time to give enough number of particles. Therefore, a TORT/MCNP coupling method has been developed to eliminate the two problems mentioned above in each code. In this method, the TORT code calculates angular flux distribution on the core surface and the MCNP code calculates neutron spectrum at the points of interest using the flux distribution. The coupling method will be used as the DOT-DOMINO-MORSE code system. This TORT/MCNP coupling method was applied to calculate the neutron flux at points where induced radioactivity data were measured for 54Mn and 60Co and the radioactivity calculations based on the neutron flux obtained from the above method were compared with the measured data.

Perturbation method for calculating impurity binding energy in an ...

Indian Academy of Sciences (India)

Nilanjan Sil

2017-12-18

Dec 18, 2017 ... Abstract. In the present paper, we have studied the binding energy of the shallow donor hydrogenic impurity, which is confined in an inhomogeneous cylindrical quantum dot (CQD) of GaAs-AlxGa1−xAs. Perturbation method is used to calculate the binding energy within the framework of effective mass ...

A grey diffusion acceleration method for time-dependent radiative transfer calculations: analysis and application

International Nuclear Information System (INIS)

Nowak, P.F.

1993-01-01

A grey diffusion acceleration method is presented and is shown by Fourier analysis and test calculations to be effective in accelerating radiative transfer calculations. The spectral radius is bounded by 0.9 for the continuous equations, but is significantly smaller for the discretized equations, especially in the optically thick regimes characteristic to radiation transport problems. The GDA method is more efficient than the multigroup DSA method because its slightly higher iteration count is more than offset by the much lower cost per iteration. A wide range of test calculations confirm the efficiency of GDA compared to multifrequency DSA. (orig.)

An axial calculation method for accurate two-dimensional PWR core simulation

International Nuclear Information System (INIS)

Grimm, P.

1985-02-01

An axial calculation method, which improves the agreement of the multiplication factors determined by two- and three-dimensional PWR neutronic calculations, is presented. The axial buckling is determined at each time point so as to reproduce the increase of the leakage due to the flattening of the axial power distribution and the effect of the axial variation of the group constants of the fuel on the reactivity is taken into account. The results of a test example show that the differences of k-eff and cycle length between two- and three-dimensional calculations, which are unsatisfactorily large if a constant buckling is used, become negligible if the results of the axial calculation are used in the two-dimensional core simulation. (Auth.)

Microfluidic devices and methods including porous polymer monoliths

Science.gov (United States)

Hatch, Anson V; Sommer, Gregory J; Singh, Anup K; Wang, Ying-Chih; Abhyankar, Vinay V

2014-04-22

Microfluidic devices and methods including porous polymer monoliths are described. Polymerization techniques may be used to generate porous polymer monoliths having pores defined by a liquid component of a fluid mixture. The fluid mixture may contain iniferters and the resulting porous polymer monolith may include surfaces terminated with iniferter species. Capture molecules may then be grafted to the monolith pores.

Integrated Power Flow and Short Circuit Calculation Method for Distribution Network with Inverter Based Distributed Generation

Directory of Open Access Journals (Sweden)

Shan Yang

2016-01-01

Full Text Available Power flow calculation and short circuit calculation are the basis of theoretical research for distribution network with inverter based distributed generation. The similarity of equivalent model for inverter based distributed generation during normal and fault conditions of distribution network and the differences between power flow and short circuit calculation are analyzed in this paper. Then an integrated power flow and short circuit calculation method for distribution network with inverter based distributed generation is proposed. The proposed method let the inverter based distributed generation be equivalent to Iθ bus, which makes it suitable to calculate the power flow of distribution network with a current limited inverter based distributed generation. And the low voltage ride through capability of inverter based distributed generation can be considered as well in this paper. Finally, some tests of power flow and short circuit current calculation are performed on a 33-bus distribution network. The calculated results from the proposed method in this paper are contrasted with those by the traditional method and the simulation method, whose results have verified the effectiveness of the integrated method suggested in this paper.

Generalized Coarse-Mesh Rebalance Method for Acceleration of Neutron Transport Calculations

International Nuclear Information System (INIS)

Yamamoto, Akio

2005-01-01

This paper proposes a new acceleration method for neutron transport calculations: the generalized coarse-mesh rebalance (GCMR) method. The GCMR method is a unified scheme of the traditional coarse-mesh rebalance (CMR) and the coarse-mesh finite difference (CMFD) acceleration methods. Namely, by using an appropriate acceleration factor, formulation of the GCMR method becomes identical to that of the CMR or CMFD method. This also indicates that the convergence property of the GCMR method can be controlled by the acceleration factor since the convergence properties of the CMR and CMFD methods are generally different. In order to evaluate the convergence property of the GCMR method, a linearized Fourier analysis was carried out for a one-group homogeneous medium, and the results clarified the relationship between the acceleration factor and the spectral radius. It was also shown that the spectral radius of the GCMR method is smaller than those of the CMR and CMFD methods. Furthermore, the Fourier analysis showed that when an appropriate acceleration factor was used, the spectral radius of the GCMR method did not exceed unity in this study, which was in contrast to the results of the CMR or the CMFD method. Application of the GCMR method to practical calculations will be easy when the CMFD acceleration is already adopted in a transport code. By multiplying a suitable acceleration factor to a coefficient (D FD ) of a finite difference formulation, one can improve the numerical instability of the CMFD acceleration method

Physics methods for calculating light water reactor increased performances

International Nuclear Information System (INIS)

Vandenberg, C.; Charlier, A.

1988-01-01

The intensive use of light water reactors (LWRs) has induced modification of their characteristics and performances in order to improve fissile material utilization and to increase their availability and flexibility under operation. From the conceptual point of view, adequate methods must be used to calculate core characteristics, taking into account present design requirements, e.g., use of burnable poison, plutonium recycling, etc. From the operational point of view, nuclear plants that have been producing a large percentage of electricity in some countries must adapt their planning to the need of the electrical network and operate on a load-follow basis. Consequently, plant behavior must be predicted and accurately followed in order to improve the plant's capability within safety limits. The Belgonucleaire code system has been developed and extensively validated. It is an accurate, flexible, easily usable, fast-running tool for solving the problems related to LWR technology development. The methods and validation of the two computer codes LWR-WIMS and MICROLUX, which are the main components of the physics calculation system, are explained

Calculation method for the seasonal performance of heat pump compact units and validation. Final report

Energy Technology Data Exchange (ETDEWEB)

Wemhoener, C.; Dott, R.; Afjei, Th. [University of Applied Sciences Northwestern Switzerland, Institute of Energy in Buildings, Muttenz (Switzerland); Huber, H.; Helfenfinger, D.; Keller, P.; Furter, R. [University of Applied Sciences Lucerne (HTA), Test center HLKS, Horw (Switzerland)

2007-02-15

This comprehensive final report for the Swiss Federal Office of Energy (SFOE) takes a look at compact heat pump units that have been developed for the heating of low energy consumption houses built to MINERGIE or MINERGIE-P standards. These units, which combine the functions of space heating, domestic hot water preparation and ventilation in one unit are described. A testing procedure developed at the University of Applied Science in Lucerne, Switzerland, using a test rig for the measurement of the seasonal performance factor (SPF) is described. A calculation method based on temperature classes for the calculation of the SPF of combined heat pump systems for space heating and domestic hot water preparation that was developed by the Institute of Energy in Buildings at the University of Applied Sciences Northwestern Switzerland is examined. Two pilot plants allowing detailed field monitoring of two compact units are described. One pilot plant installed in a single-family house built to MINERGIE standard in Gelterkinden, Switzerland, provided data on a compact unit. These results of measurements made on this and a further installation in a MINERGIE-P ultra-low energy consumption house in Zeiningen, Switzerland, are presented and discussed. Calculation methods, including exergy considerations are reviewed and their validation is discussed.

Application of Indenting Method for Calculation of Activation Energy

International Nuclear Information System (INIS)

Kim, Jong-Seog; Kim, Tae-Ryong

2006-01-01

For the calculation of activation energy of cable materials, we used to apply the break-elongation test in accordance with ASTM D412(Stand Test Methods for Rubber Properties in Tension). For the cable jacket and insulation which have regular thickness, break-elongation test had been preferred since it showed linear character in the activation energy curve. But, for the cable which has irregular thickness or rugged surface of cable inside, break-elongation test show scattered data which can not be used for the calculation of activation energy. It is not easy to prepare break-elongation specimen for the cable smaller than 13mm diameter in accordance with ASTM D412. In the cases of above, we sometime use TGA method which heat the specimen from 50 .deg. C to 700 .deg. C at heating rates of 10, 15, 20 .deg. C/min. But, TGA is suspected for the representative of natural aging in the plant since it measure the weight decreasing rate during burning which may have different aging mechanism with that of natural aging. To solve above problems, we investigated alternatives such as indenter test. Indenter test is very convenient since it does not ask for a special test specimen as the break-elongation test does. Regular surface of cable outside is the only requirement of indenter test. Experience of activation energy calculation by using the indenter test is described herein

Magnetism of hexagonal close-packed nickel calculated by full-potential linearized augmented plane wave method

International Nuclear Information System (INIS)

Tian, F.; Tian, H.; Whitmore, L.; Ye, L.Y.

2015-01-01

The energy dependent on volume of hexagonal close-packed (hcp) nickel with different magnetism is calculated by full-potential linearized augmented plane wave method. Based on the calculation ferromagnetic state is found to be the most stable state. The magnetic moment of hcp Ni is calculated and compared to those calculated by different pseudo-potential methods. Furthermore, it is also compared to that of face-centered cubic (fcc) one with the reason discussed

Application of the particle-in-cell method in propagation calculations

International Nuclear Information System (INIS)

Voelz, E.

1979-01-01

The Particle-in-Cell-Method that is capable of calculating the spreading of a plume in the atmosphere under instationary and inhomogeneous conditions, has a systematical advantage over the steady state Gaussian plume model usually used. Especially the fixed-point concentration time integral is calculated realistically instead of the locally integrated concentration at a constant time as is done in the plume model. Inaccuracies due to the computational techniques may be avoided in this way. On the other hand, at first the turbulent diffusion coefficients that describe the diffusion in the particle-in-cell method, must be prepared for all diffusion types. Thereby the diffusion coefficients can be seen to be mainly deduced in the steady state. This is one reason why they cannot be used in an optimal sense in a model that actually works instationary. (orig.) [de

Calculation methods for single-sided natural ventilation - simplified or detailed?

DEFF Research Database (Denmark)

Larsen, Tine Steen; Plesner, Christoffer; Leprince, Valérie

2016-01-01

A great energy saving potential lies within increased use of natural ventilation, not only during summer and midseason periods, where it is mainly used today, but also during winter periods, where the outdoor air holds a great cooling potential for ventilative cooling if draft problems can...... be handled. This paper presents a newly developed simplified calculation method for single-sided natural ventilation, which is proposed for the revised standard FprEN 16798-7 (earlier EN 15242:2007) for design of ventilative cooling. The aim for predicting ventilative cooling is to find the most suitable......, while maintaining an acceptable correlation with measurements on average and the authors consider the simplified calculation method well suited for the use in standards such as FprEN 16798-7 for the ventilative cooling effects from single-sided natural ventilation The comparison of different design...

Assessment of four calculation methods proposed by the EC for waste hazardous property HP 14 'Ecotoxic'.

Science.gov (United States)

Hennebert, Pierre; Humez, Nicolas; Conche, Isabelle; Bishop, Ian; Rebischung, Flore

2016-02-01

Legislation published in December 2014 revised both the List of Waste (LoW) and amended Appendix III of the revised Waste Framework Directive 2008/98/EC; the latter redefined hazardous properties HP 1 to HP 13 and HP 15 but left the assessment of HP 14 unchanged to allow time for the Directorate General of the Environment of the European Commission to complete a study that is examining the impacts of four different calculation methods for the assessment of HP 14. This paper is a contribution to the assessment of the four calculation methods. It also includes the results of a fifth calculation method; referred to as "Method 2 with extended M-factors". Two sets of data were utilised in the assessment; the first (Data Set #1) comprised analytical data for 32 different waste streams (16 hazardous (H), 9 non-hazardous (NH) and 7 mirror entries, as classified by the LoW) while the second data set (Data Set #2), supplied by the eco industries, comprised analytical data for 88 waste streams, all classified as hazardous (H) by the LoW. Two approaches were used to assess the five calculation methods. The first approach assessed the relative ranking of the five calculation methods by the frequency of their classification of waste streams as H. The relative ranking of the five methods (from most severe to less severe) is: Method 3>Method 1>Method 2 with extended M-factors>Method 2>Method 4. This reflects the arithmetic ranking of the concentration limits of each method when assuming M=10, and is independent of the waste streams, or the H/NH/Mirror status of the waste streams. A second approach is the absolute matching or concordance with the LoW. The LoW is taken as a reference method and the H wastes are all supposed to be HP 14. This point is discussed in the paper. The concordance for one calculation method is established by the number of wastes with identical classification by the considered calculation method and the LoW (i.e. H to H, NH to NH). The discordance is

Benchmark calculations for evaluation methods of gas volumetric leakage rate

International Nuclear Information System (INIS)

Asano, R.; Aritomi, M.; Matsuzaki, M.

1998-01-01

A containment function of radioactive materials transport casks is essential for safe transportation to prevent the radioactive materials from being released into environment. Regulations such as IAEA standard determined the limit of radioactivity to be released. Since is not practical for the leakage tests to measure directly the radioactivity release from a package, as gas volumetric leakages rates are proposed in ANSI N14.5 and ISO standards. In our previous works, gas volumetric leakage rates for several kinds of gas from various leaks were measured and two evaluation methods, 'a simple evaluation method' and 'a strict evaluation method', were proposed based on the results. The simple evaluation method considers the friction loss of laminar flow with expansion effect. The strict evaluating method considers an exit loss in addition to the friction loss. In this study, four worked examples were completed for on assumed large spent fuel transport cask (Type B Package) with wet or dry capacity and at three transport conditions; normal transport with intact fuels or failed fuels, and an accident in transport. The standard leakage rates and criteria for two kinds of leak test were calculated for each example by each evaluation method. The following observations are made based upon the calculations and evaluations: the choked flow model of ANSI method greatly overestimates the criteria for tests ; the laminar flow models of both ANSI and ISO methods slightly overestimate the criteria for tests; the above two results are within the design margin for ordinary transport condition and all methods are useful for the evaluation; for severe condition such as failed fuel transportation, it should pay attention to apply a choked flow model of ANSI method. (authors)

Multi-scale calculation based on dual domain material point method combined with molecular dynamics

Energy Technology Data Exchange (ETDEWEB)

Dhakal, Tilak Raj [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2017-02-27

This dissertation combines the dual domain material point method (DDMP) with molecular dynamics (MD) in an attempt to create a multi-scale numerical method to simulate materials undergoing large deformations with high strain rates. In these types of problems, the material is often in a thermodynamically non-equilibrium state, and conventional constitutive relations are often not available. In this method, the closure quantities, such as stress, at each material point are calculated from a MD simulation of a group of atoms surrounding the material point. Rather than restricting the multi-scale simulation in a small spatial region, such as phase interfaces, or crack tips, this multi-scale method can be used to consider non-equilibrium thermodynamic e ects in a macroscopic domain. This method takes advantage that the material points only communicate with mesh nodes, not among themselves; therefore MD simulations for material points can be performed independently in parallel. First, using a one-dimensional shock problem as an example, the numerical properties of the original material point method (MPM), the generalized interpolation material point (GIMP) method, the convected particle domain interpolation (CPDI) method, and the DDMP method are investigated. Among these methods, only the DDMP method converges as the number of particles increases, but the large number of particles needed for convergence makes the method very expensive especially in our multi-scale method where we calculate stress in each material point using MD simulation. To improve DDMP, the sub-point method is introduced in this dissertation, which provides high quality numerical solutions with a very small number of particles. The multi-scale method based on DDMP with sub-points is successfully implemented for a one dimensional problem of shock wave propagation in a cerium crystal. The MD simulation to calculate stress in each material point is performed in GPU using CUDA to accelerate the

An Effective Method to Accurately Calculate the Phase Space Factors for β"-β"- Decay

International Nuclear Information System (INIS)

Horoi, Mihai; Neacsu, Andrei

2016-01-01

Accurate calculations of the electron phase space factors are necessary for reliable predictions of double-beta decay rates and for the analysis of the associated electron angular and energy distributions. We present an effective method to calculate these phase space factors that takes into account the distorted Coulomb field of the daughter nucleus, yet it allows one to easily calculate the phase space factors with good accuracy relative to the most exact methods available in the recent literature.

Calculation of the shock-wave in the region close to an underground nuclear explosion (method Cades); Calcul de l'onde de choc en zone proche d'une explosion nucleaire souterraine (methode cades)

Energy Technology Data Exchange (ETDEWEB)

Supiot, F; Brugies, J [Commissariat a l' Energie Atomique, Bruyeres-le-Chatel (France). Centre d' Etudes

1969-07-01

The outline of a method is presented for calculating the characteristics of a shock wave produced by an underground nuclear explosion (pressure, wave velocity, velocity of the medium, energy left in the medium by the shock, etc.). By means of an application to a granitic medium and of a comparison with results obtained during French nuclear explosions, it has been possible to show the good agreement existing between the calculations and the experimental results. The advantages of such a method for studying the industrial applications of underground nuclear explosions are stressed. (authors) [French] On expose les grandes lignes d'une methode de calcul des caracteristiques de l'onde de choc issue d'une explosion nucleaire souterraine (pression, vitesse de l'onde, vitesse du milieu, energie deposee par le choc dans le milieu...). Une application a un milieu granitique et une comparaison aux resultats obtenus au cours d'explosions nucleaires francaises permet de montrer la bonne concordance entre le calcul et les resultats experimentaux. On souligne l'interet d'une telle, methode pour l'etude d'applications industrielles des explosions nucleaires souterraines. (auteurs)

Process control and optimization with simple interval calculation method

DEFF Research Database (Denmark)

Pomerantsev, A.; Rodionova, O.; Høskuldsson, Agnar

2006-01-01

for the quality improvement in the course of production. The latter is an active quality optimization, which takes into account the actual history of the process. The advocate approach is allied to the conventional method of multivariate statistical process control (MSPC) as it also employs the historical process......Methods of process control and optimization are presented and illustrated with a real world example. The optimization methods are based on the PLS block modeling as well as on the simple interval calculation methods of interval prediction and object status classification. It is proposed to employ...... the series of expanding PLS/SIC models in order to support the on-line process improvements. This method helps to predict the effect of planned actions on the product quality and thus enables passive quality control. We have also considered an optimization approach that proposes the correcting actions...

Calculation methods of reactivity using derivatives of nuclear power and Filter fir

International Nuclear Information System (INIS)

Diaz, Daniel Suescun

2007-01-01

This work presents two new methods for the solution of the inverse point kinetics equation. The first method is based on the integration by parts of the integral of the inverse point kinetics equation, which results in a power series in terms of the nuclear power in time dependence. Applying some conditions to the nuclear power, the reactivity is represented as first and second derivatives of this nuclear power. This new calculation method for reactivity has special characteristics, amongst which the possibility of using different sampling periods, and the possibility of restarting the calculation, after its interruption associated it with a possible equipment malfunction, allowing the calculation of reactivity in a non-continuous way. Apart from this reactivity can be obtained with or without dependency on the nuclear power memory. The second method is based on the Laplace transform of the point kinetics equations, resulting in an expression equivalent to the inverse kinetics equation as a function of the power history. The reactivity can be written in terms of the summation of convolution with response to impulse, characteristic of a linear system. For its digital form the Z-transform is used, which is the discrete version of the Laplace transform. In this method it can be pointed out that the linear part is equivalent to a filter named Finite Impulse Response (Fir). The Fir filter will always be, stable and non-varying in time, and, apart from this, it can be implemented in the non-recursive way. This type of implementation does not require feedback, allowing the calculation of reactivity in a continuous way. The proposed methods were validated using signals with random noise and showing the relationship between the reactivity difference and the degree of the random noise. (author)

Nonlinear optimization method of ship floating condition calculation in wave based on vector

Science.gov (United States)

Ding, Ning; Yu, Jian-xing

2014-08-01

Ship floating condition in regular waves is calculated. New equations controlling any ship's floating condition are proposed by use of the vector operation. This form is a nonlinear optimization problem which can be solved using the penalty function method with constant coefficients. And the solving process is accelerated by dichotomy. During the solving process, the ship's displacement and buoyant centre have been calculated by the integration of the ship surface according to the waterline. The ship surface is described using an accumulative chord length theory in order to determine the displacement, the buoyancy center and the waterline. The draught forming the waterline at each station can be found out by calculating the intersection of the ship surface and the wave surface. The results of an example indicate that this method is exact and efficient. It can calculate the ship floating condition in regular waves as well as simplify the calculation and improve the computational efficiency and the precision of results.

Calculation of proton-deuteron phase parameters including the Coulomb force

International Nuclear Information System (INIS)

Alt, E.O.; Sandhas, W.; Ziegelmann, H.

1985-04-01

A previously proposed exact method for including the Coulomb force in three-body collisions is applied to proton-deuteron scattering. We present phase shifts for angular momenta up to L=9, from elastic threshold to 50 MeV proton laboratory energy. Separable rank-one potentials are taken for the nuclear interactions. A charge-independent and a charge-symmetric choise, while leading to different neutron-deuteron and proton-deuteron phase parameters, nevertheless yields practically the same Coulomb corrections. We, moreover, investigate the question of P-wave resonances.A critical comparison of our results with those obtained in a co-ordinate space formalism is performed. Furthermore, proposals for an approximate inclusion of the Coulomb potential are tested, and found unsatisfactory. (orig.)

A method for calculating Bayesian uncertainties on internal doses resulting from complex occupational exposures

International Nuclear Information System (INIS)

Puncher, M.; Birchall, A.; Bull, R. K.

2012-01-01

Estimating uncertainties on doses from bioassay data is of interest in epidemiology studies that estimate cancer risk from occupational exposures to radionuclides. Bayesian methods provide a logical framework to calculate these uncertainties. However, occupational exposures often consist of many intakes, and this can make the Bayesian calculation computationally intractable. This paper describes a novel strategy for increasing the computational speed of the calculation by simplifying the intake pattern to a single composite intake, termed as complex intake regime (CIR). In order to assess whether this approximation is accurate and fast enough for practical purposes, the method is implemented by the Weighted Likelihood Monte Carlo Sampling (WeLMoS) method and evaluated by comparing its performance with a Markov Chain Monte Carlo (MCMC) method. The MCMC method gives the full solution (all intakes are independent), but is very computationally intensive to apply routinely. Posterior distributions of model parameter values, intakes and doses are calculated for a representative sample of plutonium workers from the United Kingdom Atomic Energy cohort using the WeLMoS method with the CIR and the MCMC method. The distributions are in good agreement: posterior means and Q 0.025 and Q 0.975 quantiles are typically within 20 %. Furthermore, the WeLMoS method using the CIR converges quickly: a typical case history takes around 10-20 min on a fast workstation, whereas the MCMC method took around 12-hr. The advantages and disadvantages of the method are discussed. (authors)

An algorithm of α-and γ-mode eigenvalue calculations by Monte Carlo method

International Nuclear Information System (INIS)

Yamamoto, Toshihiro; Miyoshi, Yoshinori

2003-01-01

A new algorithm for Monte Carlo calculation was developed to obtain α- and γ-mode eigenvalues. The α is a prompt neutron time decay constant measured in subcritical experiments, and the γ is a spatial decay constant measured in an exponential method for determining the subcriticality. This algorithm can be implemented into existing Monte Carlo eigenvalue calculation codes with minimum modifications. The algorithm was implemented into MCNP code and the performance of calculating the both mode eigenvalues were verified through comparison of the calculated eigenvalues with the ones obtained by fixed source calculations. (author)

An integral nodal variational method for multigroup criticality calculations

International Nuclear Information System (INIS)

Lewis, E.E.; Tsoulfanidis, N.

2003-01-01

An integral formulation of the variational nodal method is presented and applied to a series of benchmark critically problems. The method combines an integral transport treatment of the even-parity flux within the spatial node with an odd-parity spherical harmonics expansion of the Lagrange multipliers at the node interfaces. The response matrices that result from this formulation are compatible with those in the VARIANT code at Argonne National Laboratory. Either homogeneous or heterogeneous nodes may be employed. In general, for calculations requiring higher-order angular approximations, the integral method yields solutions with comparable accuracy while requiring substantially less CPU time and memory than the standard spherical harmonics expansion using the same spatial approximations. (author)

Some observations concerning blade-element-momentum (BEM) methods and vortex wake methods, including numerical experiments with a simple vortex model

Energy Technology Data Exchange (ETDEWEB)

Snel, H. [Netherlands Energy Research Foundation ECN, Renewable Energy, Wind Energy (Netherlands)

1997-08-01

Recently the Blade Element Momentum (BEM) method has been made more versatile. Inclusion of rotational effects on time averaged profile coefficients have improved its achievements for performance calculations in stalled flow. Time dependence as a result of turbulent inflow, pitching actions and yawed operation is now treated more correctly (although more improvement is needed) than before. It is of interest to note that adaptations in modelling of unsteady or periodic induction stem from qualitative and quantitative insights obtained from free vortex models. Free vortex methods and further into the future Navier Stokes (NS) calculations, together with wind tunnel and field experiments, can be very useful in enhancing the potential of BEM for aero-elastic response calculations. It must be kept in mind however that extreme caution must be used with free vortex methods, as will be discussed in the following chapters. A discussion of the shortcomings and the strength of BEM and of vortex wake models is given. Some ideas are presented on how BEM might be improved without too much loss of efficiency. (EG)

Development of 2-D/1-D fusion method for three-dimensional whole-core heterogeneous neutron transport calculations

International Nuclear Information System (INIS)

Lee, Gil Soo

2006-02-01

To describe power distribution and multiplication factor of a reactor core accurately, it is necessary to perform calculations based on neutron transport equation considering heterogeneous geometry and scattering angles. These calculations require very heavy calculations and were nearly impossible with computers of old days. From the limitation of computing power, traditional approach of reactor core design consists of heterogeneous transport calculation in fuel assembly level and whole core diffusion nodal calculation with assembly homogenized properties, resulting from fuel assembly transport calculation. This approach may be effective in computation time, but it gives less accurate results for highly heterogeneous problems. As potential for whole core heterogeneous transport calculation became more feasible owing to rapid development of computing power during last several years, the interests in two and three dimensional whole core heterogeneous transport calculations by deterministic method are increased. For two dimensional calculation, there were several successful approaches using even parity transport equation with triangular meshes, S N method with refined rectangular meshes, the method of characteristics (MOC) with unstructured meshes, and so on. The work in this thesis originally started from the two dimensional whole core heterogeneous transport calculation by using MOC. After successful achievement in two dimensional calculation, there were efforts in three-dimensional whole-core heterogeneous transport calculation using MOC. Since direct extension to three dimensional calculation of MOC requires too much computing power, indirect approach to three dimensional calculation was considered.Thus, 2D/1D fusion method for three dimensional heterogeneous transport calculation was developed and successfully implemented in a computer code. The 2D/1D fusion method is synergistic combination of the MOC for radial 2-D calculation and S N -like methods for axial 1

Multiband method for resonance self-shielding calculation of fuel assembly in arbitrary geometries

International Nuclear Information System (INIS)

Huang Weibin; Wu Hongchun; Cao Liangzhi; Yang Weiyan

2009-01-01

A formula to calculate the multiband parameters is derived based on the multiband method. Adopting the method combining two-band and three-band, and based on the WIMSD4-69 library, a code named RESCAL is developed. The validation shows that the results of RESCAL code are well in accordance with MCNP's, and the numerical errors meet the practical requirement. Due to the limitation of WIMSD4 69-group library and the method adopted to calculate multiband parameters, the precision of RESCAL code is highly affected by the ratio of water to uranium. (authors)

Method for calculating internal radiation and ventilation with the ADINAT heat-flow code

International Nuclear Information System (INIS)

Butkovich, T.R.; Montan, D.N.

1980-01-01

One objective of the spent fuel test in Climax Stock granite (SFTC) is to correctly model the thermal transport, and the changes in the stress field and accompanying displacements from the application of the thermal loads. We have chosen the ADINA and ADINAT finite element codes to do these calculations. ADINAT is a heat transfer code compatible to the ADINA displacement and stress analysis code. The heat flow problem encountered at SFTC requires a code with conduction, radiation, and ventilation capabilities, which the present version of ADINAT does not have. We have devised a method for calculating internal radiation and ventilation with the ADINAT code. This method effectively reproduces the results from the TRUMP multi-dimensional finite difference code, which correctly models radiative heat transport between drift surfaces, conductive and convective thermal transport to and through air in the drifts, and mass flow of air in the drifts. The temperature histories for each node in the finite element mesh calculated with ADINAT using this method can be used directly in the ADINA thermal-mechanical calculation

Accuracy of cell calculation methods used for analysis of high conversion light water reactor lattice

International Nuclear Information System (INIS)

Jeong, Chang-Joon; Okumura, Keisuke; Ishiguro, Yukio; Tanaka, Ken-ichi

1990-01-01

Validation tests were made for the accuracy of cell calculation methods used in analyses of tight lattices of a mixed-oxide (MOX) fuel core in a high conversion light water reactor (HCLWR). A series of cell calculations was carried out for the lattices referred from an international HCLWR benchmark comparison, with emphasis placed on the resonance calculation methods; the NR, IR approximations, the collision probability method with ultra-fine energy group. Verification was also performed for the geometrical modelling; a hexagonal/cylindrical cell, and the boundary condition; mirror/white reflection. In the calculations, important reactor physics parameters, such as the neutron multiplication factor, the conversion ratio and the void coefficient, were evaluated using the above methods for various HCLWR lattices with different moderator to fuel volume ratios, fuel materials and fissile plutonium enrichments. The calculated results were compared with each other, and the accuracy and applicability of each method were clarified by comparison with continuous energy Monte Carlo calculations. It was verified that the accuracy of the IR approximation became worse when the neutron spectrum became harder. It was also concluded that the cylindrical cell model with the white boundary condition was not so suitable for MOX fuelled lattices, as for UO 2 fuelled lattices. (author)

Accurate pKa calculation of the conjugate acids of alkanolamines, alkaloids and nucleotide bases by quantum chemical methods.

Science.gov (United States)

Gangarapu, Satesh; Marcelis, Antonius T M; Zuilhof, Han

2013-04-02

The pKa of the conjugate acids of alkanolamines, neurotransmitters, alkaloid drugs and nucleotide bases are calculated with density functional methods (B3LYP, M08-HX and M11-L) and ab initio methods (SCS-MP2, G3). Implicit solvent effects are included with a conductor-like polarizable continuum model (CPCM) and universal solvation models (SMD, SM8). G3, SCS-MP2 and M11-L methods coupled with SMD and SM8 solvation models perform well for alkanolamines with mean unsigned errors below 0.20 pKa units, in all cases. Extending this method to the pKa calculation of 35 nitrogen-containing compounds spanning 12 pKa units showed an excellent correlation between experimental and computational pKa values of these 35 amines with the computationally low-cost SM8/M11-L density functional approach. Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Methods for calculating population dose from atmospheric dispersion of radioactivity

Energy Technology Data Exchange (ETDEWEB)

Cohen, B L; Jow, H N; Lee, I S [Pittsburgh Univ., PA (USA)

1978-06-01

Curves are computed from which population dose (man-rem) due to dispersal of radioactivity from a point source can be calculated in the gaussian plume model by simple multiplication, and methods of using them and their limitations are considered. Illustrative examples are presented.

SU-E-T-02: 90Y Microspheres Dosimetry Calculation with Voxel-S-Value Method: A Simple Use in the Clinic

International Nuclear Information System (INIS)

Maneru, F; Gracia, M; Gallardo, N; Olasolo, J; Fuentemilla, N; Bragado, L; Martin-Albina, M; Lozares, S; Pellejero, S; Miquelez, S; Rubio, A; Otal, A

2015-01-01

Purpose: To present a simple and feasible method of voxel-S-value (VSV) dosimetry calculation for daily clinical use in radioembolization (RE) with 90 Y microspheres. Dose distributions are obtained and visualized over CT images. Methods: Spatial dose distributions and dose in liver and tumor are calculated for RE patients treated with Sirtex Medical miscrospheres at our center. Data obtained from the previous simulation of treatment were the basis for calculations: Tc-99m maggregated albumin SPECT-CT study in a gammacamera (Infinia, General Electric Healthcare.). Attenuation correction and ordered-subsets expectation maximization (OSEM) algorithm were applied.For VSV calculations, both SPECT and CT were exported from the gammacamera workstation and registered with the radiotherapy treatment planning system (Eclipse, Varian Medical systems). Convolution of activity matrix and local dose deposition kernel (S values) was implemented with an in-house developed software based on Python code. The kernel was downloaded from www.medphys.it. Final dose distribution was evaluated with the free software Dicompyler. Results: Liver mean dose is consistent with Partition method calculations (accepted as a good standard). Tumor dose has not been evaluated due to the high dependence on its contouring. Small lesion size, hot spots in health tissue and blurred limits can affect a lot the dose distribution in tumors. Extra work includes: export and import of images and other dicom files, create and calculate a dummy plan of external radiotherapy, convolution calculation and evaluation of the dose distribution with dicompyler. Total time spent is less than 2 hours. Conclusion: VSV calculations do not require any extra appointment or any uncomfortable process for patient. The total process is short enough to carry it out the same day of simulation and to contribute to prescription decisions prior to treatment. Three-dimensional dose knowledge provides much more information than other

SU-E-T-02: 90Y Microspheres Dosimetry Calculation with Voxel-S-Value Method: A Simple Use in the Clinic

Energy Technology Data Exchange (ETDEWEB)

Maneru, F; Gracia, M; Gallardo, N; Olasolo, J; Fuentemilla, N; Bragado, L; Martin-Albina, M; Lozares, S; Pellejero, S; Miquelez, S; Rubio, A [Complejo Hospitalario de Navarra, Pamplona, Navarra (Spain); Otal, A [Hospital Clinica Benidorm, Benidorm, Alicante (Spain)

2015-06-15

Purpose: To present a simple and feasible method of voxel-S-value (VSV) dosimetry calculation for daily clinical use in radioembolization (RE) with {sup 90}Y microspheres. Dose distributions are obtained and visualized over CT images. Methods: Spatial dose distributions and dose in liver and tumor are calculated for RE patients treated with Sirtex Medical miscrospheres at our center. Data obtained from the previous simulation of treatment were the basis for calculations: Tc-99m maggregated albumin SPECT-CT study in a gammacamera (Infinia, General Electric Healthcare.). Attenuation correction and ordered-subsets expectation maximization (OSEM) algorithm were applied.For VSV calculations, both SPECT and CT were exported from the gammacamera workstation and registered with the radiotherapy treatment planning system (Eclipse, Varian Medical systems). Convolution of activity matrix and local dose deposition kernel (S values) was implemented with an in-house developed software based on Python code. The kernel was downloaded from www.medphys.it. Final dose distribution was evaluated with the free software Dicompyler. Results: Liver mean dose is consistent with Partition method calculations (accepted as a good standard). Tumor dose has not been evaluated due to the high dependence on its contouring. Small lesion size, hot spots in health tissue and blurred limits can affect a lot the dose distribution in tumors. Extra work includes: export and import of images and other dicom files, create and calculate a dummy plan of external radiotherapy, convolution calculation and evaluation of the dose distribution with dicompyler. Total time spent is less than 2 hours. Conclusion: VSV calculations do not require any extra appointment or any uncomfortable process for patient. The total process is short enough to carry it out the same day of simulation and to contribute to prescription decisions prior to treatment. Three-dimensional dose knowledge provides much more information than

Improvement of calculation method for temperature coefficient of HTTR by neutronics calculation code based on diffusion theory. Analysis for temperature coefficient by SRAC code system

International Nuclear Information System (INIS)

Goto, Minoru; Takamatsu, Kuniyoshi

2007-03-01

The HTTR temperature coefficients required for the core dynamics calculations had been calculated from the HTTR core calculation results by the diffusion code with which the corrections had been performed using the core calculation results by the Monte-Carlo code MVP. This calculation method for the temperature coefficients was considered to have some issues to be improved. Then, the calculation method was improved to obtain the temperature coefficients in which the corrections by the Monte-Carlo code were not required. Specifically, from the point of view of neutron spectrum calculated by lattice calculations, the lattice model was revised which had been used for the calculations of the temperature coefficients. The HTTR core calculations were performed by the diffusion code with the group constants which were generated by the lattice calculations with the improved lattice model. The core calculations and the lattice calculations were performed by the SRAC code system. The HTTR core dynamics calculation was performed with the temperature coefficient obtained from the core calculation results. In consequence, the core dynamics calculation result showed good agreement with the experimental data and the valid temperature coefficient could be calculated only by the diffusion code without the corrections by Monte-Carlo code. (author)

Hybrid Monte-Carlo method for ICF calculations

International Nuclear Information System (INIS)

Clouet, J.F.; Samba, G.

2003-01-01

) conduction and ray-tracing for laser description. Radiation transport is usually solved by a Monte-Carlo method. In coupling diffusion approximation and transport description, the difficult part comes from the need for an implicit discretization of the emission-absorption terms: this problem was solved by using the symbolic Monte-Carlo method. This means that at each step of the simulation a matrix is computed by a Monte-Carlo method which accounts for the radiation energy exchange between the cells. Because of time step limitation by hydrodynamic motion, energy exchange is limited to a small number of cells and the matrix remains sparse. This matrix is added to usual diffusion matrix for thermal and radiative conductions: finally we arrive at a non-symmetric linear system to invert. A generalized Marshak condition describe the coupling between transport and diffusion. In this paper we will present the principles of the method and numerical simulation of an ICF hohlraum. We shall illustrate the benefits of the method by comparing the results with full implicit Monte-Carlo calculations. In particular we shall show how the spectral cut-off evolves during the propagation of the radiative front in the gold wall. Several issues are still to be addressed (robust algorithm for spectral cut- off calculation, coupling with ALE capabilities): we shall briefly discuss these problems. (authors)

Continuous energy Monte Carlo method based homogenization multi-group constants calculation

International Nuclear Information System (INIS)

Li Mancang; Wang Kan; Yao Dong

2012-01-01

The efficiency of the standard two-step reactor physics calculation relies on the accuracy of multi-group constants from the assembly-level homogenization process. In contrast to the traditional deterministic methods, generating the homogenization cross sections via Monte Carlo method overcomes the difficulties in geometry and treats energy in continuum, thus provides more accuracy parameters. Besides, the same code and data bank can be used for a wide range of applications, resulting in the versatility using Monte Carlo codes for homogenization. As the first stage to realize Monte Carlo based lattice homogenization, the track length scheme is used as the foundation of cross section generation, which is straight forward. The scattering matrix and Legendre components, however, require special techniques. The Scattering Event method was proposed to solve the problem. There are no continuous energy counterparts in the Monte Carlo calculation for neutron diffusion coefficients. P 1 cross sections were used to calculate the diffusion coefficients for diffusion reactor simulator codes. B N theory is applied to take the leakage effect into account when the infinite lattice of identical symmetric motives is assumed. The MCMC code was developed and the code was applied in four assembly configurations to assess the accuracy and the applicability. At core-level, A PWR prototype core is examined. The results show that the Monte Carlo based multi-group constants behave well in average. The method could be applied to complicated configuration nuclear reactor core to gain higher accuracy. (authors)

Ground-state inversion method applied to calculation of molecular photoionization cross-sections by atomic extrapolation: Interference effects at low energies

International Nuclear Information System (INIS)

Hilton, P.R.; Nordholm, S.; Hush, N.S.

1980-01-01

The ground-state inversion method, which we have previously developed for the calculation of atomic cross-sections, is applied to the calculation of molecular photoionization cross-sections. These are obtained as a weighted sum of atomic subshell cross-sections plus multi-centre interference terms. The atomic cross-sections are calculated directly for the atomic functions which when summed over centre and symmetry yield the molecular orbital wave function. The use of the ground-state inversion method for this allows the effect of the molecular environment on the atomic cross-sections to be calculated. Multi-centre terms are estimated on the basis of an effective plane-wave expression for this contribution to the total cross-section. Finally the method is applied to the range of photon energies from 0 to 44 eV where atomic extrapolation procedures have not previously been tested. Results obtained for H 2 , N 2 and CO show good agreement with experiment, particularly when interference effects and effects of the molecular environment on the atomic cross-sections are included. The accuracy is very much better than that of previous plane-wave and orthogonalized plane-wave methods, and can stand comparison with that of recent more sophisticated approaches. It is a feature of the method that calculation of cross-sections either of atoms or of large molecules requires very little computer time, provided that good quality wave functions are available, and it is then of considerable potential practical interest for photoelectorn spectroscopy. (orig.)

Calculation of pole profiles for the design of particle accelerators magnets. Equivalent current and conformal representation methods; Calcul de profils polaires destines a des aimants d'accelerateurs de particules. Methode des courants et methode des transformations conformes

Energy Technology Data Exchange (ETDEWEB)

Jaidane, S [Commissariat a l' Energie Atomique, Saclay (France). Centre d' Etudes Nucleaires

1968-04-01