Calculation method of Tesla coil
Коломієць, Роман Олександрович
2015-01-01
Tesla coil, despite the simplicity of its design may be called one of the least studied electronic devices. The article is an attempt to bring in various experimental results of general theoretical framework, which is the basis of exact calculation method of Tesla coils. Such calculation should be the starting point to create devices based on it. In order to develop such methods were considered the general principles of designing Tesla coil, reviewed the most famous mathematical models of its...
Methods of core neutronic calculation
Core neutronic calculations lead to the determination of geometry, composition, controls systems and to the core exploitation limits in agreement with the expected performances, with safety rules, technological choices and fuel management methods. Neutronic calculations object are described with physics justifications of hypothesis and approximations. A description and a definition of reactivity and power distribution are also given. A panorama of calculation methods used in the conception of fast breeder and pressure water reactors, are described with numerical aspects and general interest considerations related to the field of these methods and to the industrial options chosen. A complete industrial uses panorama of methods derived from the classical or generalized perturbation theory is followed by the qualification and the definition of the validity field of numerical codes.(A.B.). 88 refs., 6 figs
Friction and wear calculation methods
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
Calculation Methods for Wallenius’ Noncentral Hypergeometric Distribution
Fog, Agner
2008-01-01
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...
Methods in nuclear reactors calculations
Studies are made of the neutron transport equation corresponding to the the real and virtual reactors, as well as the starting hypotheses. Methods are developed to solve the transport equation in slab geometry, and Pl; Bl; Ml; Sn and discrete ordinates approximations. (Author)
Pile Load Capacity – Calculation Methods
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.
Assessment of seismic margin calculation methods
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.
Assessment of seismic margin calculation methods
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
Comparison of methods for calculating water erosion
Svobodová, Pavlína
2011-01-01
Bachelor thesis presents a comparison of methods for calculating water erosion. The aim is to summarize available evidence concerning the problems of water erosion. There are presented some methods how to calculate average annual erosion of soils, and selected models for calculating the erosion immediately. There are also listed possible erosion control measures through which we can at least slow the effects of erosion, rather than stop completely.
A novel algorithmic method for piezoresistance calculation
A novel algorithmic method, based on the different stress distribution on the surface of thin film in an SOI microstructure, is put forward to calculate the value of the silicon piezoresistance on the sensitive film. In the proposed method, we take the Ritz method as an initial theoretical model to calculate the rate of piezoresistance ΔR/R through an integral (the closed area Ω where the surface piezoresistance of the film lies as the integral area and the product of stress σ and piezoresistive coefficient π as the integral object) and compare the theoretical values with the experimental results. Compared with the traditional method, this novel calculation method is more accurate when applied to calculating the value of the silicon piezoresistance on the sensitive film of an SOI pieoresistive pressure sensor. (semiconductor devices)
Methods of bone marrow dose calculation
Several methods of bone marrow dose calculation for photon irradiation were analised. After a critical analysis, the author proposes the adoption, by the Instituto de Radioprotecao e Dosimetria/CNEN, of Rosenstein's method for dose calculations in Radiodiagnostic examinations and Kramer's method in case of occupational irradiation. It was verified by Eckerman and Simpson that for monoenergetic gamma emitters uniformly distributed within the bone mineral of the skeleton the dose in the bone surface can be several times higher than dose in skeleton. In this way, is also proposed the Calculation of tissue-air ratios for bone surfaces in some irradiation geometries and photon energies to be included in the Rosenstein's method for organ dose calculation in Radiodiagnostic examinations. (Author)
Relaxation Method For Calculating Quantum Entanglement
Tucci, R R
2001-01-01
In a previous paper, we showed how entanglement of formation can be defined as a minimum of the quantum conditional mutual information (a.k.a. quantum conditional information transmission). In classical information theory, the Arimoto-Blahut method is one of the preferred methods for calculating extrema of mutual information. We present a new method akin to the Arimoto-Blahut method for calculating entanglement of formation. We also present several examples computed with a computer program called Causa Comun that implements the ideas of this paper.
Reactor perturbation calculations by Monte Carlo methods
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)
Willow growing - Methods of calculation and profitability
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
DNBR limit calculation by sampling statistical method
The parametric uncertainties of DNBR and exit quality were calculated using sampling statistical method based on Wilks formula and VIPRE-W code. Then the DNBR design limit and exit quality limit were got by combining with the uncertainties of models and DNB correlation. This method can gain the more DNBR margin than RTDP methodology which is developed by Westinghouse by comparison of these two methods. (authors)
A method for tokamak neutronics calculations
This paper presents a new method for neutron transport calculation in tokamak fusion reactors. The computational procedure is based on the solution of the even-parity transport equation in a toroidal geometry. The angular neutron distribution is treated by even-parity spherical harmonic expansion, while the spatial dependence is approximated by using R-function finite elements that are defined for regions of arbitrary geometric shape. In order to test the method, calculation of a simplified tokamak model is carried out. The results are compared with the results from the literature and for the same order of accuracy a reduction of the number of spatial unknowns is shown. (author)
Multigrid Methods in Electronic Structure Calculations
Briggs, E L; Bernholc, J
1996-01-01
We describe a set of techniques for performing large scale ab initio calculations using multigrid accelerations and a real-space grid as a basis. The multigrid methods provide effective convergence acceleration and preconditioning on all length scales, thereby permitting efficient calculations for ill-conditioned systems with long length scales or high energy cut-offs. We discuss specific implementations of multigrid and real-space algorithms for electronic structure calculations, including an efficient multigrid-accelerated solver for Kohn-Sham equations, compact yet accurate discretization schemes for the Kohn-Sham and Poisson equations, optimized pseudo\\-potentials for real-space calculations, efficacious computation of ionic forces, and a complex-wavefunction implementation for arbitrary sampling of the Brillioun zone. A particular strength of a real-space multigrid approach is its ready adaptability to massively parallel computer architectures, and we present an implementation for the Cray-T3D with essen...
A New Iterative Method to Calculate [pi
Dion, Peter; Ho, Anthony
2012-01-01
For at least 2000 years people have been trying to calculate the value of [pi], the ratio of the circumference to the diameter of a circle. People know that [pi] is an irrational number; its decimal representation goes on forever. Early methods were geometric, involving the use of inscribed and circumscribed polygons of a circle. However, real…
Covariant method for calculating helicity amplitudes
We present an alternative approach for calculating helicity amplitudes for processes involving both massless and massive fermions. With this method one can easily obtain covariant expressions for the helicity amplitudes. The final expressions involve only four-vector products and are independent of the basis for γ matrices or specific form of the spinors. We use the method to obtain the helicity amplitudes for several processes involving top quark production. copyright 1996 The American Physical Society
Hemielectron method generalization for complex compound calculations
A hemielectron method is considered as applied to the systems having one open shell (n electrons on m-fold degenerated level) or several open shells. General expressions for corrections to a total energy of the system and to energies of monoelectron levels, corresponding to an unfilled shell, are obtained. Using the [RuCl6]3- complex as an example, expendiency of the approach application to the calculation of electron structure of coordination compounds is shown
Nodal expansion method for reactor core calculations
To perform realistic space dependent reactor dynamics analyses in large power reactor with all asymmetric material, control and shutdown devices, a full three dimensional calculation model is essential. A code FEMINA (Flux Expansion Method In Nodal Analysis) implementing a higher order nodal scheme employing a nodal flux expansion method in 3D is being developed. In this report the first part of this code viz., the theory of the static version and its validation with well known benchmark problems are described. The code has been found to be quite accurate as well as fast. It is available on DEC 10'', CYBER 170/730 and ND 540 computers. (author)
Acceleration methods and models in Sn calculations
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)
Methods for calculating radiation attenuation in shields
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
Discontinuous finite element methods for reactor calculations
Variational principles which employ discontinuous shape functions for the angular and/or the spatial component of the neutron flux are established to obtain numerical solutions for neutron diffusion and transport equations. Implementing discontinuous finite element methods reduces the total nodal unknowns and hence the over all computational efforts. This reduction varies from one problem to another. In this paper one group neutron transport problems are solved by varying only the order of spherical harmonic expansion for the angular component of the flux. A comparison of the solutions obtained from the discontinuous approach with either a published solutions or a conventional finite element solutions shows that the method is a very effective tool for reactor calculations
Methods of core neutronic calculation; Methodes de calcul neutronique de coeur
Bruna, G.B.; Guesdon, B. [Societe Franco-Americaine de Constructions Atomiques (FRAMATOME), 92 - Paris-La-Defense (France)
1996-02-01
Core neutronic calculations lead to the determination of geometry, composition, controls systems and to the core exploitation limits in agreement with the expected performances, with safety rules, technological choices and fuel management methods. Neutronic calculations object are described with physics justifications of hypothesis and approximations. A description and a definition of reactivity and power distribution are also given. A panorama of calculation methods used in the conception of fast breeder and pressure water reactors, are described with numerical aspects and general interest considerations related to the field of these methods and to the industrial options chosen. A complete industrial uses panorama of methods derived from the classical or generalized perturbation theory is followed by the qualification and the definition of the validity field of numerical codes.(A.B.). 88 refs., 6 figs.
ANALYTICAL METHODS FOR CALCULATING FAN AERODYNAMICS
Jan Dostal
2015-12-01
Full Text Available This paper presents results obtained between 2010 and 2014 in the field of fan aerodynamics at the Department of Composite Technology at the VZLÚ aerospace research and experimental institute in Prague – Letnany. The need for rapid and accurate methods for the preliminary design of blade machinery led to the creation of a mathematical model based on the basic laws of turbomachine aerodynamics. The mathematical model, the derivation of which is briefly described below, has been encoded in a computer programme, which enables the theoretical characteristics of a fan of the designed geometry to be determined rapidly. The validity of the mathematical model is assessed continuously by measuring model fans in the measuring unit, which was developed and manufactured specifically for this purpose. The paper also presents a comparison between measured characteristics and characteristics determined by the mathematical model as the basis for a discussion on possible causes of measured deviations and calculation deviations.
The matrix method to calculate page rank
H. Barboucha, M. Nasri
2014-06-01
Full Text Available Choosing the right keywords is relatively easy, whereas getting a high PageRank is more complicated. The index Page Rank is what defines the position in the result pages of search engines (for Google of course, but the other engines are now using more or less the same kind of algorithm. It is therefore very important to understand how this type of algorithm functions to hope to appear on the first page of results (the only page read in 95 % of cases or at least be among the first. We propose in this paper to clarify the operation of this algorithm using a matrix method and a JavaScript program enabling to experience this type of analysis. It is of course a simplified version, but it can add value to the website and achieve a high ranking in the search results and reach a larger customer base. The interest is to disclose an algorithm to calculate the relevance of each page. This is in fact a mathematical algorithm based on a web graph. This graph is formed of all the web pages that are modeled by nodes, and hyperlinks that are modeled by arcs.
METHODS OF CALCULATING THE ELECTRONIC AND ATOMIC STRUCTURES OF INTERFACES
Sutton, A
1985-01-01
Methods of calculating the electronic and atomic structures of interfaces are described. An introduction to pseudopotentials and LCAO methods is given. Methods of calculating the electronic structure of an interface with a given atomic structure are considered. The feasibility of total energy calculations, in which the atomic and electronic structures are calculated simultaneously, is discussed.
Overview of multifluid-flow-calculation methods
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
Method of integral transforms for calculating few-body reactions
Efros, V. D.; Leidemann, W.; Orlandini, G.
1998-01-01
A non-conventional approach to calculating reactions in quantum mechanics is presented. Reaction observables are obtained with bound state calculation techniques. The accuracy of the method to calculate few-nucleon response functions is discussed.
Method for consequence calculations for severe accidents
This report was commissioned by the Swedish State Power Board. The report contains a calculation of radiation doses in the surroundings caused by a theoretical core meltdown accident at Forsmark reactor No 3. The assumption used for the calculations were a 0.06% release of iodine and cesium corresponding to a 0.1% release through the FILTRA plant at Barsebaeck. The calculations were made by means of the PLUCON4 code. Meteorological data for two years from the Forsmark meteorological tower were analysed to find representative weather situations. As typical weather pasquill D was chosen with wind speed 5 m/s, and as extreme weather, Pasquill F with wind speed 2 m/s. 23 tabs., 36 ills., 21 refs. (author)
A method to calculate Fresnel lenses
In solar engineering, in contrast to image optics, Fresnel lenses are intended for securing the required concentrations of solar radiation and its distribution over a receiver's surface. It is also important to secure a high use coefficient of the concentrated flux. In particular, this defines the features of calculation of Fresnel lenses: it is necessary to take into account inaccuracies in fabrication of Fresnel lenses and solar radiation redistribution by means of selecting the respective parameters of Fresnel lens belts. In the present work, we examine the procedure for the calculating geometrical parameters of Fresnel lenses on a flat base by considering the mentioned requirements. A corresponding software for calculating the geometrical parameters and concentrating characteristics of the Fresnel lenses is developed, and examples of calculation are given. For a constant refractive index of Fresnel lens material, it is shown that the Fresnel lens can secure a concentration of about 1000, but in this case the optical efficiency of the Fresnel lens will not be higher than 70%. The procedure that has been developed may be the basic one for determining the parameters and concentrating characteristics of Fresnel lenses by considering refractive index variance. (author)
Analytic methods for calculating coupling impedances
These lecture notes describe a variety of analytic techniques to calculate the longitudinal and transverse impedances of obstacles in a beam pipe. They also treat the effort to shield these impedances from the beam by appropriate use of thin conducting layers. (orig.)
Computational methods for probability of instability calculations
Wu, Y.-T.; Burnside, O. H.
1990-01-01
This paper summarizes the development of the methods and a computer program to compute the probability of instability of a dynamic system than can be represented by a system of second-order ordinary linear differential equations. Two instability criteria based upon the roots of the characteristics equation or Routh-Hurwitz test functions are investigated. Computational methods based on system reliability analysis methods and importance sampling concepts are proposed to perform efficient probabilistic analysis. Numerical examples are provided to demonstrate the methods.
Method for consequence calculations for severe accidents
This report was commissioned by the Swedish State Power Board. The report contains a calculation of radiation doses in the surroundings caused by a theoretical core meltdown accident at Ringhals reactor No 3/4. The accident sequence chosen for the calcualtions was a release caused by total power failure. The calculations were made by means of the PLUCON4 code. A decontamination factor of 500 is used to account for the scrubber effect. Meteorological data for two years from the Ringhals meteorological tower were analysed to find representative weather situations. As typical weather, Pasquill D, was chosen with a wind speed of 10 m/s, and as extreme weather, Pasquill E, with a wind speed of 2 m/s. 19 refs. (author)
COSTS CALCULATION OF TARGET COSTING METHOD
UNGUREANU Sebastian
2014-01-01
Cost information system plays an important role in every organization in the decision making process. An important task of management is ensuring control of the operations, processes, sectors, and not ultimately on costs. Although in achieving the objectives of an organization compete more control systems (production control, quality control, etc.), the cost information system is important because monitors results of the other. Detailed analysis of costs, production cost calculation, quantifi...
Methods for calculating anisotropic transfer cross sections
The Legendre moments of the group transfer cross section, which are widely used in the numerical solution of the transport calculation can be efficiently and accurately constructed from low-order (K = 1--2) successive partial range moments. This is convenient for the generation of group constants. In addition, a technique to obtain group-angle correlation transfer cross section without Legendre expansion is presented. (author)
Monte Carlo method application to shielding calculations
CANDU spent fuel discharged from the reactor core contains Pu, so it must be stressed in two directions: tracing for the fuel reactivity in order to prevent critical mass formation and personnel protection during the spent fuel manipulation. The basic tasks accomplished by the shielding calculations in a nuclear safety analysis consist in dose rates calculations in order to prevent any risks both for personnel protection and impact on the environment during the spent fuel manipulation, transport and storage. To perform photon dose rates calculations the Monte Carlo MORSE-SGC code incorporated in SAS4 sequence from SCALE system was used. The paper objective was to obtain the photon dose rates to the spent fuel transport cask wall, both in radial and axial directions. As source of radiation one spent CANDU fuel bundle was used. All the geometrical and material data related to the transport cask were considered according to the shipping cask type B model, whose prototype has been realized and tested in the Institute for Nuclear Research Pitesti. (authors)
Calculation of radon concentration in water by toluene extraction method
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.)
A review of calculation methods for fast and intermediate reactors
This paper discusses the development of methods for calculating intermediate and fast reactors. It deals with various approaches to the problems of physical calculation. The calculation of resonance effects is discussed. Consideration is given to multi-group systems of fundamental and conjugate equations, various applications of perturbation theory to the problems of physical reactor calculation, and numerical methods of solving fundamental and conjugate reactor equations, which approximate the method of spherical harmonics. The paper describes an application of the response method to the solution of critical-mass problems, and methods of calculating reactors with hydrogeneous moderators. The fundamental features of an effective one-group reactor model are described. (author)
COSTS CALCULATION OF TARGET COSTING METHOD
Sebastian UNGUREANU
2014-06-01
Full Text Available Cost information system plays an important role in every organization in the decision making process. An important task of management is ensuring control of the operations, processes, sectors, and not ultimately on costs. Although in achieving the objectives of an organization compete more control systems (production control, quality control, etc., the cost information system is important because monitors results of the other. Detailed analysis of costs, production cost calculation, quantification of losses, estimate the work efficiency provides a solid basis for financial control. Knowledge of the costs is a decisive factor in taking decisions and planning future activities. Managers are concerned about the costs that will appear in the future, their level underpinning the supply and production decisions as well as price policy. An important factor is the efficiency of cost information system in such a way that the information provided by it may be useful for decisions and planning of the work.
Investigation of Calculation Techniques of Finite Difference Method
Audrius Krukonis
2011-03-01
Full Text Available Finite difference method used for microstrip transmission line analysis is considered in this article. Paper mainly deals with iterative and bound matrix calculation techniques of finite difference method. Mathematical model for microstrip transmission line electrical potential calculations using both techniques is described. Results of characteristic impedance calculation using iterative and bound matrix techniques are presented and analyzed.Article in Lithuanian
Approximate methods in gamma-ray skyshine calculations
An approximate computational method for gamma-ray skyshine calculations is described. The method is suitable for a source collimated uniformly about the vertical and accounts for uniform overhead concrete shielding above the source. Results of calculations are compared to measurements as well as results of other calculations
Transportation channels calculation method in MATLAB
Output devices and charged particles transport channels are necessary components of any modern particle accelerator. They differ both in sizes and in terms of focusing elements depending on particle accelerator type and its destination. A package of transport line designing codes for magnet optical channels in MATLAB environment is presented in this report. Charged particles dynamics in a focusing channel can be studied easily by means of the matrix technique. MATLAB usage is convenient because its information objects are matrixes. MATLAB allows the use the modular principle to build the software package. Program blocks are small in size and easy to use. They can be executed separately or commonly. A set of codes has a user-friendly interface. Transport channel construction consists of focusing lenses (doublets and triplets). The main of the magneto-optical channel parameters are total length and lens position and parameters of the output beam in the phase space (channel acceptance, beam emittance - beam transverse dimensions, particles divergence and image stigmaticity). Choice of the channel operation parameters is based on the conditions for satisfying mutually competing demands. And therefore the channel parameters calculation is carried out by using the search engine optimization techniques.
Calculation methods for neutron radiography spatial resolution
Spatial resolution is an important parameter for neutron radiography facility. In this paper, different methods to define the spatial resolution,such as point spread function (PSF), line spread function (LSF), edge spread function (ESF) and modulation transfer function (MTF), are analyzed and compared. MTF turns out to be the best, as it is derived from the linear system theory in a given frequency domain, and gives the maximum amount of useful information on system signal modulation. (authors)
Progress and prospects of calculation methods for radiation shielding
Progress in calculation methods for radiation shielding are reviewed based on the activities of research committees related to radiation shielding fields established in the Atomic Energy Society of Japan. A technological roadmap for the field of radiation shielding; progress and prospects for specific shielding calculation methods such as the Monte Carlo, discrete ordinate Sn transport, and simplified methods; and shielding experiments used to validate calculation methods are presented in this paper. (author)
Nodal methods in numerical reactor calculations
The present work describes the antecedents, developments and applications started in 1972 with Prof. Hennart who was invited to be part of the staff of the Nuclear Engineering Department at the School of Physics and Mathematics of the National Polytechnic Institute. Since that time and up to 1981, several master theses based on classical finite element methods were developed with applications in point kinetics and in the steady state as well as the time dependent multigroup diffusion equations. After this period the emphasis moved to nodal finite elements in 1, 2 and 3D cartesian geometries. All the thesis were devoted to the numerical solution of the neutron multigroup diffusion and transport equations, few of them including the time dependence, most of them related with steady state diffusion equations. The main contributions were as follows: high order nodal schemes for the primal and mixed forms of the diffusion equations, block-centered finite-differences methods, post-processing, composite nodal finite elements for hexagons, and weakly and strongly discontinuous schemes for the transport equation. Some of these are now being used by several researchers involved in nuclear fuel management. (Author)
Sectoring method for cosmic radiation shielding calculation for LEO satellite
One of an approximate calculation model (sectoring method) is developed for a cosmic radiation shielding in satellite. Shielding calculation is performed for KITSAT-1 at the assumed SAA (South Atlatic Anomaly) location with AP-8 model radiation spectrum. When sectoring method is applied, calculation error is expected compared with 3-D detailed geometry calculation because of straight knock-on assumption neglecting the deflection of incident proton. However, sectoring method shows good agreements with 3-dimensional detailed Monte Carlo calculation in two TID detector locations
Quantum Monte Carlo diagonalization method as a variational calculation
Mizusaki, Takahiro; Otsuka, Takaharu [Tokyo Univ. (Japan). Dept. of Physics; Honma, Michio
1997-05-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)
The TAB method for numerical calculation of spray droplet breakup
Orourke, P. J.; Amsden, A. A.
A short history is given of the major milestones in the development of the stochastic particle method for calculating liquid fuel sprays. The most recent advance has been the discovery of the importance of drop breakup in engine sprays. A new method, called TAB, for calculating drop breakup is presented. Some theoretical properties of the method are derived; its numerical implementation in the computer program KIVA is described; and comparisons are presented between TAB-method calculations and experiments and calculations using another breakup model.
Comments on Simplified Calculation Method for Fire Exposed Concrete Columns
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 exposed...... 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....
A Damping Characteristics Calculation Method of Metal Dry Friction Isolators
JIANG Hong-yuan; HAO De-gang; XIA Yu-hong; ULANOV A M; PONOMAREV Yu K
2008-01-01
The dry friction ring-type vibration isolator is considered as an isotropic continuous medium. A method of dry friction hysteresis loop calculation is proposed based on friction force analysis of contact beam. The friction force is modeled as an equivalent distributed moment to use the finite element method (FEM) to calculate the dry friction vibration isolator hysteresis loop, so the damping characteristics can be obtained. A comparison of the hysteresis loop calculation results and the experimental results shows the average relative error is 2.7%, it proves the calculation method is feasible.
Some methods for calculation of perturbations in nuclear reactors
Abramov, B. D., E-mail: abramov@ippe.ru [Leypunsky Institute of Physics and Power Engineering (Russian Federation)
2015-12-15
Some methods for calculation of local perturbations of neutron fields and reactivity effects accompanying them are considered. Existence, uniqueness, properties and methods for finding solutions to the considered problems are discussed.
Evolution of calculation methods taking into account severe accidents
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
A New Method for Calculating the Thermoelectric Efficiency
吴一东; 王志敏; 何元金
2004-01-01
We present an approximate method for calculating the thermoelectric effciency. The method has a high precision and is applicable to almost all of the thermoelectric devices. The expression for the thermoelectric efficiency we obtained does not involve the position variable, so the calculations are simplified greatly.
Calculation of transonic flows using an extended integral equation method
Nixon, D.
1976-01-01
An extended integral equation method for transonic flows is developed. In the extended integral equation method velocities in the flow field are calculated in addition to values on the aerofoil surface, in contrast with the less accurate 'standard' integral equation method in which only surface velocities are calculated. The results obtained for aerofoils in subcritical flow and in supercritical flow when shock waves are present compare satisfactorily with the results of recent finite difference methods.
Iterative acceleration methods for Monte Carlo and deterministic criticality calculations
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
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.
Approximate design calculation methods for radiation streaming in shield irregularities
Investigation and assessment are made for approximate design calculation methods of radiation streaming in shield irregularities. Investigation is made for (1) source, (2) definition of streaming radiation components, (3) calculation methods of streaming radiation, (4) streaming formulas for each irregularity, (5) difficulties in application of streaming formulas, etc. Furthermore, investigation is made for simple calculation codes and albedo data. As a result, it is clarified that streaming calculation formulas are not enough to cover various irregularities and their accuracy or application limit is not sufficiently clear. Accurate treatment is not made in the formulas with respect to the radiation behavior for slant incidence, bend part, offset etc., that results in too much safety factors in the design calculation and distrust of the streaming calculation. To overcome the state and improve the accuracy of the design calculation for shield irregularities, it is emphasized to assess existing formulas and develop better formulas based on systematic experimental studies. (author)
无
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.
A finite element method for SSI time history calculation
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
Ore reserve calculation methods used by Eldorado Nuclear Limited
Uranium-bearing pitchblende deposits of the Beaverlodge area in northern Saskatchewan are highly complex. The ideas and concepts of ore reserve calculation methods employed by Eldorado Nuclear Limited in assessing and planning the mining of these deposits are described. A manual block-system of ore reserve calculation was used before the adoption of the current computerized system. Four classifications are used for ore reserves calculated by the system, which provides two main program jobs for calculating ore reserves and several additional ones that involve calculations and graphical presentation of ore reserve information for use in mine planning. A comparison of production statistics and ore reserve calculations illustrates the accuracy of the method. (author)
Comparison of calculational methods for EBT reactor nucleonics
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
New method for calculation of integral characteristics of thermal plumes
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...... of a sitting occupant. The improvement in calculation of the characteristics of the thermal plume achieved with the developed method, in comparison with methods used and reported in the literature, is demonstrated....... 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...
Elongation method for electronic structure calculations of random DNA sequences.
Orimoto, Yuuichi; Liu, Kai; Aoki, Yuriko
2015-10-30
We applied ab initio order-N elongation (ELG) method to calculate electronic structures of various deoxyribonucleic acid (DNA) models. We aim to test potential application of the method for building a database of DNA electronic structures. The ELG method mimics polymerization reactions on a computer and meets the requirements for linear scaling computational efficiency and high accuracy, even for huge systems. As a benchmark test, we applied the method for calculations of various types of random sequenced A- and B-type DNA models with and without counterions. In each case, the ELG method maintained high accuracy with small errors in energy on the order of 10(-8) hartree/atom compared with conventional calculations. We demonstrate that the ELG method can provide valuable information such as stabilization energies and local densities of states for each DNA sequence. In addition, we discuss the "restarting" feature of the ELG method for constructing a database that exhaustively covers DNA species. PMID:26337429
Fluid-structure interaction calculations using a linear perturbation method
Aim of the work is to present and validate FSI (Fluid-Structure Interaction) calculations by using a linear perturbation method and commercial Computational Fluid Dynamics (CFD) and structural analysis codes. Star-CD is used for CFD calculations and ABAQUS for structural analysis. The external MpCCI code is used for coupling the CFD and structural analysis codes
Protein-protein binding affinities calculated using the LIE method
Andberg, Tor Arne Heim
2011-01-01
Absolute binding free energies for the third domain of the turkey ovomucoid inhibitor in complex with Streptomyces griseus proteinase B and porcine pancreatic elastase has been calculated using the linear interaction energy method.
A new method to calculate hole subbands of semiconductors hetrostructures
A new method, ordinary boundary method (OBM), is suggested to calculate the hole subbands and wavefunctions of semiconductor hetrostructures. Based on it, a transfer matrix is established. By means of it not only the hole subbands and wavefunctions of hetrostructures can be calculated rigorously and quickly but also the asymptotic transfer method can be generalized to calculate the hole subbands of complex systems in which the valence bandedges are oblique lines or curves. As examples we calculated the hole subbands of GaAs/AlxGa1-xAs quantum well and GaAs/AlxGa1-xAs superlattice in terms of both OBM and transfer matrix method. (author). 7 refs, 4 figs
A valid method of calculating virtual scene depth
QUAN Hong-yan; ZHANG Tian-wen; LIN Xiang-hong
2005-01-01
A valid method of virtual scene depth calculating is put forward. In this method cameras rotate in three different viewpoints in the plane and we calculate the depth of panorama using three stitching cylinder panoramas. In the investigation, the column of panorama is regarded as a slot image. Using the conic intersected by the epipolar plane and the cylinder, we can obtain the perpendicularity disparity. In order to obtain dense correspondence fast and accurately, a new method of obtaining horizontal disparity using depth continuity is also put forward. It converts the problem of panorama dense correspondence to the problem of searching points in the conic. The occlusion problem is dealt with using three cylinders in the depth calculation. It is verified that this method is convenient, useful and efficient in calculating the depth of a virtual scene.
Comparison study on cell calculation method of fast reactor
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)
Seismic response based on transient calculations. Spectral and stochastic methods
Further to the recent development in the ASTER code of functionalities enabling random dynamic responses to be calculated, notably a stochastic type seismic analysis, we propose a combination of three calculation methods to estimate a probabilistic seismic response on an N4 reactor building stick-model. Transient calculations involves time-and cost-consuming repetition. The conventional oscillator response spectrum calculation yields only the maximum response expectation. On the other hand, the stochastic approach in this context gives the response corresponding to selected probabilities. (authors). 12 figs., 3 tabs., 4 refs
Methods for tornado frequency calculation of nuclear power plant
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)
Calculating expected deaths: a comparison of two methods.
Chovil, A C
1987-01-01
A comparison is presented between the traditional "person-years" and more recently described "prospective model" methods for calculating mortality expectations. Problems arising from the fact that expectations under the person-years method are calculated on the basis that a null hypothesis is true, which results in artificial figures that, at least theoretically, are meaningless if the hypothesis is rejected, are discussed. Data are presented from two studies in which expectations have been c...
CALCULATION OF COMPANY COSTS THROUGH THE DIRECT-COSTING CALCULATION METHOD
Florin-Constantin DIMA
2013-06-01
Full Text Available The cost of production has as its starting point the purchase cost of raw materials and consumables, as well as their processing cost and the calculation of the production cost involves complex aspects. This article is based on the two major concepts of costs calculation, namely the concept of full costs and the concept of partial costs, and it analyses the direct-costing calculation method. Necessity of the Development of calculation methods to ensure rapid determination of the cost of production, and the establishment of indicators broad spectrum of information necessary for making decisions to streamline a business activity conducted by direct-costing method. Direct-costing method appeared in the U.S. for the first time in 1934 (applied by Jonathan Harris and G. Charter Harrison. Subsequently, this method was applied to European countries (England, France, Germany etc.. We stopped on this method because it is considered a modern method of costing. Therefore, we analyzed both advantages and limitations of the method in question
Simplified hourly method to calculate summer temperatures in dwellings
Mortensen, Lone Hedegaard; Aggerholm, Søren
2012-01-01
with further simplifications. The method is used for calculating room temperatures for all hours of a reference year. It is essential that the simplified method is able to predict the temperature in the room with the highest heat load. The heat load is influenced by the solar load, internal load...... and solar load. The developed method can calculate the number of hours above a given temperature limit. The limits are a prerequisite for the development of the simplified method, and a supplementary maximum temperature limit is suggested to ensure robustness. The setting of the ventilation rate is...
Development of 3-D detailed FBR core calculation method based on method of characteristics
A new detailed 3-D transport calculation method taking into account the heterogeneity of fuel assemblies has been developed in hexagonal-z geometry by combining the method of characteristics (MOC) and the nodal transport method. From the nodal transport calculation which uses assembly homogenized cross sections, the axial leakage is calculated, and it is used for the MOC calculation which treats the heterogeneity of fuel assemblies. Series of homogeneous MOC calculations which use assembly homogeneous cross sections are carried out to obtain effective cross sections, which preserve assembly reaction rates. This effective cross sections are again used in the 3-dimensional nodal transport calculation. The numerical calculations have been performed to verify 3-dimensional radial calculations of FBR (fast breeder reactor) 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. (authors)
Comparison of MCNPX and Albedo method in criticality calculation
This study aims to conduct a computer simulation that will calculate the reactivity of a homogeneous reactor and compare the results with the calculations made by the albedo method. The simulation will be developed using the MCNPX. The study compared the results calculated for a hypothetical reactor by the albedo method for four groups of energy with those obtained by the MCNPX simulation. The design of the reactor is spherical and homogeneous with a reflector of finite thickness. The value obtained for the neutron effective multiplication factor - keff will be compared. Different situations were simulated in order to obtain results closer to the compared method and reality. The was Good consistency could be noticed between the calculated results. (author)
Calculation method for gamma dose rates from Gaussian puffs
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 Rxy. 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 spherical puffs
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.)
According to the acceptance of ICRP Publication 60 (1990), the dose equivalent limit for the boarder of controlled area will be defined as 1.3 mSv/3 months in the Regulation for the Enforcement of the Medical Service Law which is scheduled to be revised. The calculating methods of radiation shielding to be considered are as follows: The first method is calculating the dose equivalent for each nuclide using 3-month maximum estimated use dose. The second method is calculating the dose equivalent using 3-month maximum estimated use dose after the conversion of all nuclide dose into that of 131I. The third method is calculating the dose equivalent using 1 day maximum estimated use dose after the conversion of all nuclide dose into that of 131I. We've investigated which of methods can meet the new regulation value (1.3 mSv/3 months). In modeled facility, we've tried to calculate the dose by the first method to confirm if we can perform the reasonable control in safe. Total dose equivalent for the boarder of controlled area (B) was 883 μSv/3 months by the first method, and its value turned out to be about 1/4 of that of the third method. Only the result by the first method was found to be within the confines of new dose equivalent limit of 1.3 mSv/3 months. The results of both method the second and the third were found to be within the confines of existing dose equivalent limit. The method as to calculate the shielding for each nuclide by using 3-month maximum estimated use dose has been accepted in the Law Concerning Prevention from Radiation Hazards due to Radioisotopes, etc. As the method is practically in accordance with the current use of radioisotope in nuclear medicine facility, the possibility of it coping with the new dose equivalent limit was indicated. (author)
Application of nonparametric statistic method for DNBR limit calculation
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)
A simple method for calculating Clebsch-Gordan coefficients
Klink, W H; Wickramasekara, S, E-mail: william-klink@uiowa.ed, E-mail: wickrama@grinnell.ed, E-mail: s-wickram@uiowa.ed [Department of Physics and Astronomy, University of Iowa, Iowa City, IA 52242 (United States)
2010-09-15
This paper presents a simple method for calculating Clebsch-Gordan coefficients for the tensor product of two unitary irreducible representations (UIRs) of the rotation group. The method also works for multiplicity-free irreducible representations appearing in the tensor product of any number of UIRs of the rotation group. The generalization to representations with multiplicity is straightforward and briefly discussed.
Efficient Calculation of Near Fields in the FDTD Method
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...
New methods for neutron response calculations with MCNP
MCNP4B was released for international distribution in February, 1997. The author summarized the new MCNP4B features since the release of MCNP4A over three years earlier and compare some results. Then he describes new methods being developed for future code releases. The focus is methods and applications of ex-core neutron response calculations
Experimental verification of a recursive method to calculate evapotranspiration
Recently, a recursive combination method (RCM) to calculate potential and crop evapotranspiration (ET) was given by Lascano and Van Bavel (Agron. J. 2007, 99:585–590). The RCM differs from the Penman-Monteith (PM) method, the main difference being that the assumptions made regarding the temperature ...
Calculation of Solar Radiation by Using Regression Methods
Kızıltan, Ö.; Şahin, M.
2016-04-01
In this study, solar radiation was estimated at 53 location over Turkey with varying climatic conditions using the Linear, Ridge, Lasso, Smoother, Partial least, KNN and Gaussian process regression methods. The data of 2002 and 2003 years were used to obtain regression coefficients of relevant methods. The coefficients were obtained based on the input parameters. Input parameters were month, altitude, latitude, longitude and landsurface temperature (LST).The values for LST were obtained from the data of the National Oceanic and Atmospheric Administration Advanced Very High Resolution Radiometer (NOAA-AVHRR) satellite. Solar radiation was calculated using obtained coefficients in regression methods for 2004 year. The results were compared statistically. The most successful method was Gaussian process regression method. The most unsuccessful method was lasso regression method. While means bias error (MBE) value of Gaussian process regression method was 0,274 MJ/m2, root mean square error (RMSE) value of method was calculated as 2,260 MJ/m2. The correlation coefficient of related method was calculated as 0,941. Statistical results are consistent with the literature. Used the Gaussian process regression method is recommended for other studies.
On-line reactivity calculation using Lagrange method
Highlights: • Lagrange method is proposed for on-line reactivity calculation in nuclear reactors. • The need for nuclear power history or the Laplace transform is vanished. • The three- and five-point formulas are presented and examined in different benchmark cases. • Computational time-steps of up to 1 s lead to highly reliable reactivity calculations. • The main advantage of the proposed approach is its stability and convergence in large time-step calculations. - Abstract: In this paper, a novel multi-step method is proposed for solving the inverse point kinetics problem using Lagrange polynomial method. By use of this approach, the need for nuclear power history or the Laplace transform is vanished. Furthermore, the accuracy of the method is of order hn for the (n + 1)-point formula, where h is the computational time-step. The three- and five-point formulas of the Lagrange method are used for on-line reactivity calculations and results are benchmarked against reference solutions for different nuclear power forms. Moreover, results for different computational time-steps are compared in each case. The results show the accuracy of the proposed method in all benchmarking cases. For slow transients (large reactor periods), it is shown that time-steps of up to 1 s lead to highly reliable reactivity calculations. However, the optimal time-step in almost all cases is shown to be 0.1 s. The main advantage of the proposed approach, in contrast with previous numerical methods, is its stability and convergence in large time-step calculations. The proposed method can be used as real time reactivity meter in all nuclear reactors without limitation of nuclear power form
Calculating atomic and molecular properties using variational Monte Carlo methods
The authors compute a number of properties for the 1 1S, 21S, and 23S states of helium as well as the ground states of H2 and H/+3 using Variational Monte Carlo. These are in good agreement with previous calculations (where available). Electric-response constants for the ground states of helium, H2 and H+3 are computed as derivatives of the total energy. The method used to calculate these quantities is discussed in detail
Refinement of thermal imager minimum resolvable temperature difference calculating method
Kolobrodov, V. G.; Mykytenko, V. I.
2015-11-01
Calculating methods, which accurately predict minimum resolvable temperature difference (MRTD), are of significant interest for many years. The article deals with improvement the accuracy of determining the thermal imaging system MRTD by elaboration the visual perception model. We suggest MRTD calculating algorithm, which is based on a reliable approximation of the human visual system modulation transfer function (MTF) proposed by N. Nill. There was obtained a new expression for the bandwidth evaluation, which is independent of angular size of the Foucault bar target.
The analytic method for calculating the control rod worth
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
Neutron transport calculations using Quasi-Monte Carlo methods
Moskowitz, B.S.
1997-07-01
This paper examines the use of quasirandom sequences of points in place of pseudorandom points in Monte Carlo neutron transport calculations. For two simple demonstration problems, the root mean square error, computed over a set of repeated runs, is found to be significantly less when quasirandom sequences are used ({open_quotes}Quasi-Monte Carlo Method{close_quotes}) than when a standard Monte Carlo calculation is performed using only pseudorandom points.
Nuclear data and multigroup methods in fast reactor calculations
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.)
Shape integral method for magnetospheric shapes. [boundary layer calculations
Michel, F. C.
1979-01-01
A method is developed for calculating the shape of any magnetopause to arbitrarily high precision. The method uses an integral equation which is evaluated for a trial shape. The resulting values of the integral equation as a function of auxiliary variables indicate how close one is to the desired solution. A variational method can then be used to improve the trial shape. Some potential applications are briefly mentioned.
Improvement of fitting method of multiband parameters for cell calculations
To accurately perform cell calculations of nuclear reactors, a new fitting procedure has been developed for calculating multiband parameters, which are necessary for effective cross section calculations. By using the new fitting procedure, the error of multiband parameters becomes always zero. Reactor cell calculations have been performed to compare the effective cross sections and the infinite multiplication factors etc. calculated using the multiband parameters obtained by the new and the conventional fitting procedures by using the cross section set based on the JENDL-3.1 library with 107 energy groups. It is found that there is a small difference of the calculational results between the two fitting procedures and it is found from burnup calculations that the difference of the infinite multiplication factors is not dependent on the burnup period up to about 30 GWd/t. The onion skin effect can be exactly treated by dividing a fuel pellet to multiple regions and by using the multiband method. Thus the difference of burnup properties between two fitting procedures are investigated for the divided and the undivided fueled cells. The total inventory of Pu, Am etc. at the divided case is almost the same to the undivided case at the end of the burnup period. However it is found that the radial distribution of atomic density is slightly different between the two fitting procedures. (author)
Composite electron propagator methods for calculating ionization energies
Díaz-Tinoco, Manuel; Dolgounitcheva, O.; Zakrzewski, V. G.; Ortiz, J. V.
2016-06-01
Accurate ionization energies of molecules may be determined efficiently with composite electron-propagator (CEP) techniques. These methods estimate the results of a calculation with an advanced correlation method and a large basis set by performing a series of more tractable calculations in which large basis sets are used with simpler approximations and small basis sets are paired with more demanding correlation techniques. The performance of several CEP methods, in which diagonal, second-order electron propagator results with large basis sets are combined with higher-order results obtained with smaller basis sets, has been tested for the ionization energies of closed-shell molecules from the G2 set. Useful compromises of accuracy and computational efficiency employ complete-basis-set extrapolation for second-order results and small basis sets in third-order, partial third-order, renormalized partial-third order, or outer valence Green's function calculations. Analysis of results for vertical as well as adiabatic ionization energies leads to specific recommendations on the best use of regular and composite methods. Results for 22 organic molecules of interest in the design of photovoltaic devices, benzo[a]pyrene, Mg-octaethylporphyrin, and C60 illustrate the capabilities of CEP methods for calculations on large molecules.
Composite electron propagator methods for calculating ionization energies.
Díaz-Tinoco, Manuel; Dolgounitcheva, O; Zakrzewski, V G; Ortiz, J V
2016-06-14
Accurate ionization energies of molecules may be determined efficiently with composite electron-propagator (CEP) techniques. These methods estimate the results of a calculation with an advanced correlation method and a large basis set by performing a series of more tractable calculations in which large basis sets are used with simpler approximations and small basis sets are paired with more demanding correlation techniques. The performance of several CEP methods, in which diagonal, second-order electron propagator results with large basis sets are combined with higher-order results obtained with smaller basis sets, has been tested for the ionization energies of closed-shell molecules from the G2 set. Useful compromises of accuracy and computational efficiency employ complete-basis-set extrapolation for second-order results and small basis sets in third-order, partial third-order, renormalized partial-third order, or outer valence Green's function calculations. Analysis of results for vertical as well as adiabatic ionization energies leads to specific recommendations on the best use of regular and composite methods. Results for 22 organic molecules of interest in the design of photovoltaic devices, benzo[a]pyrene, Mg-octaethylporphyrin, and C60 illustrate the capabilities of CEP methods for calculations on large molecules. PMID:27305999
Solar opacity calculations using the super-transition-array method
Krief, M; Gazit, D
2016-01-01
An opacity model based on the Super-Transition-Array (STA) method for the calculation of monochromatic opacities of local thermodynamic equilibrium plasmas was developed. The model is described and used to calculate spectral opacities for a solar model implementing the recent AGSS09 composition. Calculations are carried throughout the solar radiative zone. The relative contributions of different chemical elements and photon-matter processes to the total Rosseland mean opacity are analyzed in detail. Monochromatic opacities and charge state distributions were compared with the widely used Opacity-Project (OP) code, for several elements near the radiation-convection interface. STA 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 of the full AGSS09 photospheric mixture, including all heavy metals was performed. It was shown that due to their extremely low abundance, and despite being very go...
Calculations of EURACOS iron benchmark experiment using the HYBRID method
In this paper, the HYBRID method is used in the calculations of the iron benchmark experiment at the EURACOS-II device. The saturation activities of the 32S(n,p)32P reaction at different depths in an iron block are computed with ENDF/B-IV data to compare with the measurements. At the outer layers of the iron block, the HYBRID calculation gives increasingly higher results than the VITAMIN-C multigroup calculation. With the adjustment of the two- to one-dimensional ratios, the HYBRID results agree with the measurements to within 10% at most penetration depths, a considerable improvement over the VITAMIN-C multigroup results. The development of a collapsing method for the HYBRID cross sections provides a more direct and practical way of using the HYBRID method in the two-dimensional calculations. It is observed that half of the window effect is smeared in the collapsing treatment, but it still provides a better cross-section set than the VITAMIN-C cross sections for the deep-penetration calculations
A method to calculate displacement factors using SVM
Li Peixian; Tan Zhixiang; Yan Lili; Deng Kazhong
2011-01-01
In order to improve the precision of mining subsidence prediction,a mathematical model using Support Vector Machine (SVM) was established to calculate the displacement factor.The study is based on a comprehensive analysis of factors affecting the displacement factor,such as mechanical properties of the cover rock,the ratio of mining depth to seam thickness,dip angle of the coal seam and the thickness of loose layer.Data of 63 typical observation stations were used as a training and testing sample set.A SVM regression model of the displacement factor and the factors affecting it was established with a kernel function,an insensitive loss factor and a properly selected penalty factor.Given an accurate calculation algorithm for testing and analysis,the results show that an SVM regression model can calculate displacement factor precisely and reliable precision can be obtained which meets engineering requirements.The experimental results show that the method to calculation of the displacement factor,based on the SVM method,is feasible.The many factors affecting the displacement factor can be considered with this method.The research provides an efficient and accurate approach for the calculation of displacement in mining subsidence prediction.
A New Method to Calculate Internal Rate of Return
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.
Optimization method for quantitative calculation of clay minerals in soil
Libo Hao; Qiaoqiao Wei; Yuyan Zhao; Zilong Lu; Xinyun Zhao
2015-04-01
Determination of types and amounts for clay minerals in soil are important in environmental, agricultural, and geological investigations. Many reliable methods have been established to identify clay mineral types. 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 determination of clay minerals in soil based on bulk chemical composition data. The fundamental principles and processes of the calculation are elucidated. Some samples were used for reliability verification of the method and the results prove the simplicity and efficacy of the approach.
Investigation of methods used in calculations of solar cell parameters
Shvets, E. Ya.; Khrypko, S. L.; Zubko, E. I.
2009-01-01
Analytical expressions have been obtained for extracting the electrical parameters and characteristics of solar cells, including series and shunt resistances, and the saturation current. The method of Lagrange multipliers was used for computing the shape factor of the current–voltage characteristic (CVC) of solar cell. The calculation results demonstrated a satisfactory agreement with experimental data.
Emergy Algebra: Improving Matrix Methods for Calculating Tranformities
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...
Thick-Restart Lanczos Method for Electronic Structure Calculations
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
Theories and calculation methods for regional objective ET (evapotranspiration): Applications
LIU diaHong; QIN DaYong; WANG MingNa; L(U) JinYan; SANG XueFeng; ZHANG RuiMei
2009-01-01
The regional objective ET (evapotranspiration) is defined as the quantity of water that could be con-sumed in a particular region. It varies with the water conditions and economic development stages in the region. It is also constrained by the requirement of benign environment cycle. At the same time, it must meet the demands of sustainable economic growth and the construction of harmony society.Objective ET based water resources distribution will replace the conventional method, which empha-sizes the balance between the water demand and the water supply. It puts focus on the reasonable water consumption instead of the forecasted water demand, which is usually greater than the actual one. In this paper, we calculated the objective ET of 2010 year level in Tianjin by an analysis-integra-tion-assessment method. Objective ET can be classified into two parts: controllable ET and uncontrol-lable ET. Controllable ET includes the ET from irrigation land and the ET from resident land, among which the former can be calculated with soil moisture model and evapotranspiration model, while the latter can be calculated by water use ration and water consumption rate. The uncontrollable ET can be calculated with the distributed hydrological model and the remote sensing monitoring model. The two models can be mutually calibrated. In this paper, eight schemes are put forward based on different portfolios of water resources. The objective ET of each scheme was calculated and the results were assessed and analyzed. Finally, an optimal scheme was recommended.
Studies on the calculation method of regional solar radiation
Studies on the Calculation Method of Regional Solar Radiation 1. The significance and question of regional solar radiation The significance of regional solar radiation in agriculture is clear. To estimate regional agricultural producing potential, we need to know the regional solar radiation. In the field of hydrology, regional solar radiation is also important to estimate evapotranspiration of the region. There are so many slopes with different slope angles and slope directions in a region. So, we have to know how we can calculate slope radiation. The conversion
Load calculation methods for offshore wind turbine foundations
Passon, Patrik; Branner, Kim
2014-01-01
turbine manufacturer provides the FD with dynamic responses obtained from aeroelastic simulations at a predefined interface. These responses are subsequently expanded to the corresponding dynamic responses in all structural parts of the foundation. In this article, a novel load calculation method, for the......Calculation of design loads for offshore wind turbine (OWT) foundations is typically performed in a joint effort between wind turbine manufactures and foundation designers (FDs). Ideally, both parties would apply the same fully integrated design tool and model for that purpose. However, such...
Slab Storage Calculation Method for Continuous Casting-Hot Rolling
PENG Qi-chun; LIU Qing; TIAN Nai-yuan
2004-01-01
Based on load-oriented manufacturing control theory, different combining modes and slab storage calculation method for continuous casting and hot rolling were discussed. The buffer capacity index of continuous casting-rolling was introduced, and the reasonable slab storage under different combining modes was calculated with buffer capacity index of 120.00 h for CCR, 79.20 h for HCR, 19.68 h for DHCR and 3.84 h for DR. Thin slab is 1.20 h, and the strip is zero. Theory gist was provided for steel enterprise to decrease storage.
A new method for the automatic calculation of prosody
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)
Short circuit currents calculation by using the impedance correction method
Schaefer, W.
1985-01-01
By introduction of correction factors for the impedances of generators and power station unit transformers into the known guidings according to VDE 0102 (basis principle of short circuit current calculation) it is possible to limit the deviation to less then +-5% of the value of the complete superposition method apart from some exceptions. The correction factors are valid for calculation of the maximum short circuit currents considering the practically admissible operating conditions. The separated correction of the impedances of generators and power station unit transformers can be used without restriction for all short circuit points of interest.
Simplified method for trace anomaly calculations in d≤6
We discuss a simplified method for computing trace anomalies in dimensions d≤6. It is known that in the quantum mechanical approach trace anomalies in d dimensions are given by a (d/2+1)-loop computation in an auxiliary 1D sigma model with arbitrary geometry. We show how one can obtain the same information using a simpler (d/2)-loop calculation on an arbitrary geometry supplemented by a (d/2+1)-loop calculation on the simplified geometry of a maximally symmetric space
Refinement of the substructure method for integral transport calculations
A new generalization of the interface-current method for coupling two-dimensional heterogeneous assemblies, called substructures, has been developed. The method has been designed for fine-structure burnup calculations in large, very heterogeneous media. For the calculations, the medium is divided into rectangular substructures, which can have internal symmetries, containing rectangular and/or cylindrical structure elements, divided into homogeneous zones. A zonewise flat or linear expansion is used to formulate a direct collision-probability problem within each substructure. The substructures are coupled by making a piecewise uniform or linear expansion for the partial currents entering and leaving the substructures. The method has also been used to implement an approximate piecewise isotropic reflection for two-dimensional x-y collision probabilities calculations. The accuracies and computing times achieved are illustrated by one-group fixed-source numerical calculations for a typical 7 x 7 pin pressurized water reactor assembly as well as for a set of fuel slabs imbedded in a water moderator
Calculation of interface curvature with the level-set method
Lervåg, Karl Yngve
2014-01-01
The level-set method is a popular method for interface capturing. One of the advantages of the level-set method is that the curvature and the normal vector of the interface can be readily calculated from the level-set function. However, in cases where the level-set method is used to capture topological changes, the standard discretization techniques for the curvature and the normal vector do not work properly. This is because they are affected by the discontinuities of the signed-distance function half-way between two interfaces. This article addresses the calculation of normal vectors and curvatures with the level-set method for such cases. It presents a discretization scheme that is relatively easy to implement in to an existing code. The improved discretization scheme is compared with a standard discretization scheme, first for a case with no flow, then for a case where two drops collide in a shear flow. The results show that the improved discretization yields more robust calculations in areas where topolo...
Accelerating molecular property calculations with nonorthonormal Krylov space methods
Furche, Filipp; Krull, Brandon T.; Nguyen, Brian D.; Kwon, Jake
2016-05-01
We formulate Krylov space methods for large eigenvalue problems and linear equation systems that take advantage of decreasing residual norms to reduce the cost of matrix-vector multiplication. The residuals are used as subspace basis without prior orthonormalization, which leads to generalized eigenvalue problems or linear equation systems on the Krylov space. These nonorthonormal Krylov space (nKs) algorithms are favorable for large matrices with irregular sparsity patterns whose elements are computed on the fly, because fewer operations are necessary as the residual norm decreases as compared to the conventional method, while errors in the desired eigenpairs and solution vectors remain small. We consider real symmetric and symplectic eigenvalue problems as well as linear equation systems and Sylvester equations as they appear in configuration interaction and response theory. The nKs method can be implemented in existing electronic structure codes with minor modifications and yields speed-ups of 1.2-1.8 in typical time-dependent Hartree-Fock and density functional applications without accuracy loss. The algorithm can compute entire linear subspaces simultaneously which benefits electronic spectra and force constant calculations requiring many eigenpairs or solution vectors. The nKs approach is related to difference density methods in electronic ground state calculations and particularly efficient for integral direct computations of exchange-type contractions. By combination with resolution-of-the-identity methods for Coulomb contractions, three- to fivefold speed-ups of hybrid time-dependent density functional excited state and response calculations are achieved.
An economical method to calculate eigenvalues of the Schroedinger equation
A method is presented which is an extension to negative energies of a spectral integral equation method to solve the Schroedinger equation, developed previously for scattering applications. One important innovation is a re-scaling procedure in order to compensate for the exponential behaviour of the negative energy Green's function. Another is the need to find approximate energy eigenvalues, to serve as starting values for a subsequent iteration procedure. In order to illustrate the new method, the binding energy of the He-He dimer is calculated, using the He-He TTY potential. In view of the small value of the binding energy, the wavefunction has to be calculated to a distance of 3000 au. Two hundred and twenty mesh points were sufficient to obtain an accuracy of three significant figures for the binding energy, and with 320 mesh points the accuracy increased to six significant figures. An application to a potential with two wells, separated by a barrier, is also made
Calculation of Turbulent Boundary Layers Using the Dissipation Integral Method
MatthiasBuschmann
1999-01-01
This paper gives an introduction into the dissipation integral method.The general integral equations for the three-dimensional case are derved.It is found that for a practical calculation algorithm the integral monentum equation and the integral energy equation are msot useful.Using Two different sets of mean velocity profiles the hyperbolical character of a dissipation integral method is shown.Test cases for two-and three-dimensional boundary layers are analysed and discussed.The paper concludes with a discussion of the advantages and limits of dissipation integral methods.
Analytical methods to calculate correlation functions in quantum statistical physics
In the work there is presented a brief but clear and quite reserved account of two analytical methods to calculate correlation functions in quantum statistical physics, proceeding from the first principles, i.e., the most broadly used at present two-time temperature Green's functions method and a new, so-called 'direct algebraic' method (DAM). The aim of this work is to show on the concrete examples of live the most broadly used models of quantum statistical physics, mathematical and technical clarity and simplicity of DAM and hence its practical value
Comparison between calculation methods of dose rates in gynecologic brachytherapy
In treatments with radiations for gynecologic tumors is necessary to evaluate the quality of the results obtained by different calculation methods for the dose rates on the points of clinical interest (A, rectal, vesicle). The present work compares the results obtained by two methods. The Manual Calibration Method (MCM) tri dimensional (Vianello E., et.al. 1998), using orthogonal radiographs for each patient in treatment, and the Theraplan/T P-11 planning system (Thratonics International Limited 1990) this last one verified experimentally (Vianello et.al. 1996). The results show that MCM can be used in the physical-clinical practice with a percentile difference comparable at the computerized programs. (Author)
Artificial Boundary Method for Calculating Ship Wave Resistance
文新; 韩厚德
2003-01-01
The calculation of wave resistance for a ship moving at constant speed near a free surface is considered. This wave resistance is calculated with a linearized steady potential model. To deal with the unboundedness of the physical domain in the potential flow problem, we introduce one vertical side as an artificial upstream boundary and two vertical sides as the artificial downstream boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the potential flow problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. The solution of the variational problem by the finite element method gives the numerical approximation of the potential flow around the ship, which was used to calculate the wave resistance. The numerical examples show the accuracy and efficiency of the proposed numerical scheme.
The application of advanced rotor (performance) methods for design calculations
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.
Implicit calculations of transonic flows using monotone methods
Goorjian, P. M.; van Buskirk, R.
1981-01-01
Implicit approximate-factorization algorithms have been developed that use monotone methods for the calculation of steady and unsteady transonic flows governed by the small-disturbance-potential equation. These algorithms use the new Engquist-Osher switch in the type-dependent differencing in place of the standard Murman-Cole switch. The resulting algorithms are more stable; hence, calculations can be done more efficiently. For steady flows, the convergence rate is about 35% faster, and for unsteady flows the allowable time step is about 10 times larger. These improvements are achieved with no increase in computer storage and with only minor modifications in codes that use the Murman-Cole switch. Also an implicit algorithm has been developed for the steady full-potential equation in one-dimension, which uses monotone methods.
Methods and computer codes for nuclear systems calculations
B P Kochurov; A P Knyazev; A Yu Kwaretzkheli
2007-02-01
Some numerical methods for reactor cell, sub-critical systems and 3D models of nuclear reactors are presented. The methods are developed for steady states and space–time calculations. Computer code TRIFON solves space-energy problem in (, ) systems of finite height and calculates heterogeneous few-group matrix parameters of reactor cells. These parameters are used as input data in the computer code SHERHAN solving the 3D heterogeneous reactor equation for steady states and 3D space–time neutron processes simulation. Modification of TRIFON was developed for the simulation of space–time processes in sub-critical systems with external sources. An option of SHERHAN code for the system with external sources is under development.
Problems in radiation shielding calculations with Monte Carlo methods
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)
Helicity methods in LO and NLO QCD calculations
The goal of this thesis is the acceleration of numerical calculations of QCD observables, both at leading order and next-to-leading order in the coupling constant. In particular, the optimization of helicity and spin summation in the context of VEGAS Monte Carlo algorithms is investigated. In the literature, two such methods are mentioned but without detailed analyses. Only one of these methods can be used at next-to-leading order. This work presents a total of five different methods that replace the helicity sums with a Monte Carlo integration. This integration can be combined with the existing phase space integral, in the hope that this causes less overhead than the complete summation. For three of these methods, an extension to existing subtraction terms is developed which is required to enable next-to-leading order calculations. All methods are analyzed with respect to efficiency, accuracy, and ease of implementation before they are compared with each other. In this process, one method shows clear advantages in relation to all others.
Benchmark calculations for evaluation methods of gas volumetric leakage rate
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)
Testing the QA Method for Calculating Jet v_{2}
Mueller, Jason
2014-01-01
For the summer, I was assigned to work on the ALICE experiment with Alice Ohlson. I wrote several programs throughout the summer that were used to calculate jet v 2 using a non-standard method described by my supervisor in her Ph.D. thesis. Though the project is not yet complete, significant progress has been made, and the results so far seem promising.
Calculations of pair production by Monte Carlo methods
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
Calculations of pair production by Monte Carlo methods
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.
The New Performance Calculation Method of Fouled Axial Flow Compressor
Huadong Yang; Hong Xu
2014-01-01
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 in...
Fast Stiffness Matrix Calculation for Nonlinear Finite Element Method
Emir Gülümser
2014-01-01
Full Text Available We propose a fast stiffness matrix calculation technique for nonlinear finite element method (FEM. Nonlinear stiffness matrices are constructed using Green-Lagrange strains, which are derived from infinitesimal strains by adding the nonlinear terms discarded from small deformations. We implemented a linear and a nonlinear finite element method with the same material properties to examine the differences between them. We verified our nonlinear formulation with different applications and achieved considerable speedups in solving the system of equations using our nonlinear FEM compared to a state-of-the-art nonlinear FEM.
Methods for calculating conjugate problems of heat transfer
Kalinin, E. K.; Dreitser, G. A.; Kostiuk, V. V.; Berlin, I. I.
Methods are examined for calculating various conjugate problems of heat transfer in channels and closed vessels in cases of single-phase and two-phase flow in steady and unsteady conditions. The single-phase-flow studies involve the investigation of gaseous and liquid heat-carriers in pipes, annular and plane channels, and pipe bundles in cases of cooling and heating. General relationships are presented for heat transfer in cases of film, transition, and nucleate boiling, as well as for boiling crises. Attention is given to methods for analyzing the filling and cooling of conduits and tanks by cryogenic liquids; and ways to intensify heat transfer in these conditions are examined.
Unfolding method for first-principles LCAO electronic structure calculations
Lee, Chi-Cheng; Yamada-Takamura, Yukiko; Ozaki, Taisuke
2013-08-01
Unfolding the band structure of a supercell to a normal cell enables us to investigate how symmetry breakers such as surfaces and impurities perturb the band structure of the normal cell. We generalize the unfolding method, originally developed based on Wannier functions, to the linear combination of atomic orbitals (LCAO) method, and present a general formula to calculate the unfolded spectral weight. The LCAO basis set is ideal for the unfolding method because the basis functions allocated to each atomic species are invariant regardless of the existence of surface and impurity. The unfolded spectral weight is well defined by the property of the LCAO basis functions. In exchange for the property, the non-orthogonality of the LCAO basis functions has to be taken into account. We show how the non-orthogonality can be properly incorporated in the general formula. As an illustration of the method, we calculate the dispersive quantized spectral weight of a ZrB2 slab and show strong spectral broadening in the out-of-plane direction, demonstrating the usefulness of the unfolding method.
Unfolding method for first-principles LCAO electronic structure calculations
Unfolding the band structure of a supercell to a normal cell enables us to investigate how symmetry breakers such as surfaces and impurities perturb the band structure of the normal cell. We generalize the unfolding method, originally developed based on Wannier functions, to the linear combination of atomic orbitals (LCAO) method, and present a general formula to calculate the unfolded spectral weight. The LCAO basis set is ideal for the unfolding method because the basis functions allocated to each atomic species are invariant regardless of the existence of surface and impurity. The unfolded spectral weight is well defined by the property of the LCAO basis functions. In exchange for the property, the non-orthogonality of the LCAO basis functions has to be taken into account. We show how the non-orthogonality can be properly incorporated in the general formula. As an illustration of the method, we calculate the dispersive quantized spectral weight of a ZrB2 slab and show strong spectral broadening in the out-of-plane direction, demonstrating the usefulness of the unfolding method. (paper)
Comparison of dose calculation methods for brachytherapy of intraocular tumors
Rivard, Mark J.; Chiu-Tsao, Sou-Tung; Finger, Paul T.; Meigooni, Ali S.; Melhus, Christopher S.; Mourtada, Firas; Napolitano, Mary E.; Rogers, D. W. O.; Thomson, Rowan M.; Nath, Ravinder [Department of Radiation Oncology, Tufts University School of Medicine, Boston, Massachusetts 02111 (United States); Quality MediPhys LLC, Denville, New Jersey 07834 (United States); New York Eye Cancer Center, New York, New York 10065 (United States); Department of Radiation Oncology, Comprehensive Cancer Center of Nevada, Las Vegas, Nevada 89169 (United States); Department of Radiation Oncology, Tufts University School of Medicine, Boston, Massachusetts 02111 (United States); Department of Radiation Physics, University of Texas, M.D. Anderson Cancer Center, Houston, Texas 77030 (United States) and Department of Experimental Diagnostic Imaging, University of Texas, M.D. Anderson Cancer Center, Houston, Texas 77030 (United States); Physics, Elekta Inc., Norcross, Georgia 30092 (United States); Department of Physics, Carleton University, Ottawa, Ontario K1S 5B6 (Canada); Department of Therapeutic Radiology, Yale University School of Medicine, New Haven, Connecticut 06520 (United States)
2011-01-15
Purpose: To investigate dosimetric differences among several clinical treatment planning systems (TPS) and Monte Carlo (MC) codes for brachytherapy of intraocular tumors using {sup 125}I or {sup 103}Pd plaques, and to evaluate the impact on the prescription dose of the adoption of MC codes and certain versions of a TPS (Plaque Simulator with optional modules). Methods: Three clinical brachytherapy TPS capable of intraocular brachytherapy treatment planning and two MC codes were compared. The TPS investigated were Pinnacle v8.0dp1, BrachyVision v8.1, and Plaque Simulator v5.3.9, all of which use the AAPM TG-43 formalism in water. The Plaque Simulator software can also handle some correction factors from MC simulations. The MC codes used are MCNP5 v1.40 and BrachyDose/EGSnrc. Using these TPS and MC codes, three types of calculations were performed: homogeneous medium with point sources (for the TPS only, using the 1D TG-43 dose calculation formalism); homogeneous medium with line sources (TPS with 2D TG-43 dose calculation formalism and MC codes); and plaque heterogeneity-corrected line sources (Plaque Simulator with modified 2D TG-43 dose calculation formalism and MC codes). Comparisons were made of doses calculated at points-of-interest on the plaque central-axis and at off-axis points of clinical interest within a standardized model of the right eye. Results: For the homogeneous water medium case, agreement was within {approx}2% for the point- and line-source models when comparing between TPS and between TPS and MC codes, respectively. For the heterogeneous medium case, dose differences (as calculated using the MC codes and Plaque Simulator) differ by up to 37% on the central-axis in comparison to the homogeneous water calculations. A prescription dose of 85 Gy at 5 mm depth based on calculations in a homogeneous medium delivers 76 Gy and 67 Gy for specific {sup 125}I and {sup 103}Pd sources, respectively, when accounting for COMS-plaque heterogeneities. For off
Solar Opacity Calculations Using the Super-transition-array Method
Krief, M.; Feigel, A.; Gazit, D.
2016-04-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.
Investigations on the Migration Mode Method (MMM) for reactor calculations
Current calculation codes for reactor analysis are based on the multi-group method to evaluate energy distribution of neutron flux. Usually a two energy group diffusion equation is adopted. This choice is adequate for PWRs associated to cross sections libraries tabulated versus fuel exposure and other state parameters as moderator density, fuel temperature, boron concentration. An improvement of this approach is represented by the Migration Mode Method (MMM) by which the neutron spectrum is expanded in terms of base functions corresponding to the different modes of migration of the neutrons in the energy dimension. For a thermal reactor, three such functions may be easily identified: the Maxwellian distribution of the neutrons at thermal equilibrium with the moderator, the 1/E slowing down distribution (corrected to take into account resonance absorption effects) and the fission neutron spectrum. The (space-dependent) coefficients of the expansion are calculated by solving a differential equation which results having a structure similar to the one relevant to multi-group theory. The method can therefore be easily implemented adopting existing diffusion theory codes. With the present work, some investigations on the MMM are described relevant to UO2 fuelled PWR systems. Demonstrative results are given to validate the potentiality of the method. (authors)
Investigations on the migration mode method (MMM) for reactor calculations
Current calculation codes for reactor analysis are based on the multi-group method to evaluate energy distribution of neutron flux. Usually a two energy group diffusion equation is adopted. This choice is adequate for PWRs associated to cross sections libraries tabulated versus fuel exposure and other state parameters as moderator density, fuel temperature, boron concentration. An improvement of this approach is represented by the migration mode method (MMM) by which the neutron spectrum is expanded in terms of base functions corresponding to the different modes of migration of the neutrons in the energy dimension. For a thermal reactor, three such functions may be easily identified: the Maxwellian distribution of the neutrons at thermal equilibrium with the moderator, the 1/E slowing down distribution (corrected to take into account resonance absorption effects) and the fission neutron spectrum. The (space-dependent) coefficients of the expansion are calculated by solving a differential equation which results having a structure similar to the one relevant to multi-group theory. The method can therefore be easily implemented adopting existing diffusion theory codes. With the present work, some investigations on the MMM are described relevant to UO2 fuelled PWR systems. Demonstrative results are given to validate the potentiality of the method
The new performance calculation method of fouled axial flow compressor.
Yang, Huadong; Xu, Hong
2014-01-01
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. PMID:25197717
The New Performance Calculation Method of Fouled Axial Flow Compressor
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.
Investigations on the migration mode method (MMM) for reactor calculations
Dall' Osso, Aldo [AREVA NP, Tour Areva, 92084 Paris La Defense Cedex (France); Gandini, Augusto [DINCE, University of Rome ' La Sapienza' , Piazzale Aldo Moro 5, 00185 Rome (Italy)], E-mail: augusto.gandini@uniroma1.it; Rotella, Rossella [S.R.S. GROUP S.r.l., Via dei Prefetti 26, 00186 Rome (Italy)
2008-07-15
Current calculation codes for reactor analysis are based on the multi-group method to evaluate energy distribution of neutron flux. Usually a two energy group diffusion equation is adopted. This choice is adequate for PWRs associated to cross sections libraries tabulated versus fuel exposure and other state parameters as moderator density, fuel temperature, boron concentration. An improvement of this approach is represented by the migration mode method (MMM) by which the neutron spectrum is expanded in terms of base functions corresponding to the different modes of migration of the neutrons in the energy dimension. For a thermal reactor, three such functions may be easily identified: the Maxwellian distribution of the neutrons at thermal equilibrium with the moderator, the 1/E slowing down distribution (corrected to take into account resonance absorption effects) and the fission neutron spectrum. The (space-dependent) coefficients of the expansion are calculated by solving a differential equation which results having a structure similar to the one relevant to multi-group theory. The method can therefore be easily implemented adopting existing diffusion theory codes. With the present work, some investigations on the MMM are described relevant to UO{sub 2} fuelled PWR systems. Demonstrative results are given to validate the potentiality of the method.
Convergence and accuracy of numerical methods for trajectory calculations
Computation of trajectories by a kinematic method requires the numerical solution of the differential equation by which the trajectory is defined. A widely used method is the iterative scheme of Petterssen which has second-order accuracy. The convergence and accuracy of this scheme is investigated for some simple flow types where analytical solutions are available. The results are discussed in comparison to other methods, especially a method similar to the Patterssen scheme, which has been recommended for use in semi-Lagrangian advection schemes. The truncation error in trajectory calculations should be kept about one order of magnitude smaller than the total uncertainty, which is mainly due to analysis errors and limited resolution of the wind data. It is shown that for trajectory calculations based on the typical output of current numerical weather prediction models or comparable data, this requires a time step 15% of the time step necessary to achieve convergence. If a fixed time step is used, it should not exceed about 0.5 h. Flexible time steps, based on the estimate of the truncation error which is provided by the difference between the first and the second iteration, are an attractive alternative. 26 refs., 8 figs
Calculation of VPP basing on functional analyzing method
Bai Kaixiang; Wang Dexun; Han Jiurui
2007-01-01
The establishment and realization of the VPP calucation's model for the functional analytic theory are discussed in this paper. Functional analyzing method is a theoretical model of the VPP calculation which can eliminate the influence of the sail and board's size skillfully, so it can be regarded as a brief standard of the sailboard's VPP results. As a brief watery dynamical model, resistance on board can be regarded as a direct proportion to the square of the boat-velocity. The boat-velocities at the state of six wind-velocities (3 m/s-8 m/s) with angles of 25°-180° are obtained by calculating, which provides an important gist of the sailing-route's selection in upwind-sailing.
Calculation method of solar radiation incident upon slopes considering topography
When radiation in a basin is calculated, slope inclination, slope orientation and topography of surroundings have to be taken into account. The method of approximation to topography by triangles proposed by Miura et al. is employed to take slope characteristics and topography of surroundings into account. Authors prepared 360 directions' shades altitudes, i.e. every degree of angle, for each triangle in advance, and used these shades' altitudes to calculate both direct radiation on a slope diffuse radiation taking topography of surroundings into account. And authors show how to estimate hourly direct and diffuse solar radiation from hourly horizontal global radiation and synthesize hourly slope global radiation on slopes
Large-scale atomic calculations using variational methods
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.
Large-scale atomic calculations using variational methods
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 7p2P 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
A unique manual method for emergency offsite dose calculations
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
Nested element method in multidimensional neutron diffusion calculations
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, here presented, 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 priory 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
Nested element method in multidimensional neutron diffusion calculations
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
Large Scale Electronic Structure Calculations using Quantum Chemistry Methods
Scuseria, Gustavo E.
1998-03-01
This talk will address our recent efforts in developing fast, linear scaling electronic structure methods for large scale applications. Of special importance is our fast multipole method( M. C. Strain, G. E. Scuseria, and M. J. Frisch, Science 271), 51 (1996). (FMM) for achieving linear scaling for the quantum Coulomb problem (GvFMM), the traditional bottleneck in quantum chemistry calculations based on Gaussian orbitals. Fast quadratures(R. E. Stratmann, G. E. Scuseria, and M. J. Frisch, Chem. Phys. Lett. 257), 213 (1996). combined with methods that avoid the Hamiltonian diagonalization( J. M. Millam and G. E. Scuseria, J. Chem. Phys. 106), 5569 (1997) have resulted in density functional theory (DFT) programs that can be applied to systems containing many hundreds of atoms and ---depending on computational resources or level of theory-- to many thousands of atoms.( A. D. Daniels, J. M. Millam and G. E. Scuseria, J. Chem. Phys. 107), 425 (1997). Three solutions for the diagonalization bottleneck will be analyzed and compared: a conjugate gradient density matrix search (CGDMS), a Hamiltonian polynomial expansion of the density matrix, and a pseudo-diagonalization method. Besides DFT, our near-field exchange method( J. C. Burant, G. E. Scuseria, and M. J. Frisch, J. Chem. Phys. 105), 8969 (1996). for linear scaling Hartree-Fock calculations will be discussed. Based on these improved capabilities, we have also developed programs to obtain vibrational frequencies (via analytic energy second derivatives) and excitation energies (through time-dependent DFT) of large molecules like porphyn or C_70. Our GvFMM has been extended to periodic systems( K. N. Kudin and G. E. Scuseria, Chem. Phys. Lett., in press.) and progress towards a Gaussian-based DFT and HF program for polymers and solids will be reported. Last, we will discuss our progress on a Laplace-transformed \\cal O(N^2) second-order pertubation theory (MP2) method.
Score Calculation in Informatics Contests Using Multiple Criteria Decision Methods
Jurate SKUPIENE
2011-04-01
Full Text Available The Lithuanian Informatics Olympiad is a problem solving contest for high school students. The work of each contestant is evaluated in terms of several criteria, where each criterion is measured according to its own scale (but the same scale for each contestant. Several jury members are involved in the evaluation. This paper analyses the problem how to calculate the aggregated score for whole submission in the above mentioned situation. The chosen methodology for solving this problem is Multiple Criteria Decision Analysis (MCDA. The outcome of this paper is the score aggregation method proposed to be applied in LitIO developed using MCDA approaches.
Dose calculation of 6 MV Truebeam using Monte Carlo method
The purpose of this work is to simulate 6 MV Varian Truebeam linac dosimeter characteristics using Monte Carlo method and to investigate the availability of phase space file and the accuracy of the simulation. With the phase space file at linac window supplied by Varian to be a source, the patient-dependent part was simulated. Dose distributions in a water phantom with a 10 cm × 10 cm field were calculated and compared with measured data for validation. Evident time reduction was obtained from 4-5 h which a whole simulation cost on the same computer to around 48 minutes. Good agreement between simulations and measurements in water was observed. Dose differences are less than 3% for depth doses in build-up region and also for dose profiles inside the 80% field size, and the effect in penumbra is good. It demonstrate that the simulation using existing phase space file as the EGSnrc source is efficient. Dose differences between calculated data and measured data could meet the requirements for dose calculation. (authors)
A CNS calculation line based on a Monte Carlo method
Full text: The design of the moderator cell of a Cold Neutron Source (CNS) involves many different considerations regarding geometry, location, and materials. Decisions taken in this sense affect not only the neutron flux in the source neighborhood, which can be evaluated by a standard empirical method, but also the neutron flux values in experimental positions far away of the neutron source. At long distances from the neutron source, very time consuming 3D deterministic methods or Monte Carlo transport methods are necessary in order to get accurate figures. Standard and typical terminology such as average neutron flux, neutron current, angular flux, luminosity, are magnitudes very difficult to evaluate in positions located several meters away from the neutron source. The Monte Carlo method is a unique and powerful tool to transport neutrons. Its use in a bootstrap scheme appears to be an appropriate solution for this type of systems. The proper use of MCNP as the main tool leads to a fast and reliable method to perform calculations in a relatively short time with low statistical errors. The design goal is to evaluate the performance of the neutron sources, their beam tubes and neutron guides at specific experimental locations in the reactor hall as well as in the neutron or experimental hall. In this work, the calculation methodology used to design Cold, Thermal and Hot Neutron Sources and their associated Neutron Beam Transport Systems, based on the use of the MCNP code, is presented. This work also presents some changes made to the cross section libraries in order to cope with cryogenic moderators such as liquid hydrogen and liquid deuterium. (author)
Study on calculation methods for the effective delayed neutron fraction
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)
Method for calculating annual energy efficiency improvement of TV sets
Varman, M. [Department of Mechanical Engineering, University of Malaya, Lembah Pantai, 50603 Kuala Lumpur (Malaysia); Mahlia, T.M.I. [Department of Mechanical Engineering, University of Malaya, Lembah Pantai, 50603 Kuala Lumpur (Malaysia)]. E-mail: indra@um.edu.my; Masjuki, H.H. [Department of Mechanical Engineering, University of Malaya, Lembah Pantai, 50603 Kuala Lumpur (Malaysia)
2006-10-15
The popularization of 24 h pay-TV, interactive video games, web-TV, VCD and DVD are poised to have a large impact on overall TV electricity consumption in the Malaysia. Following this increased consumption, energy efficiency standard present a highly effective measure for decreasing electricity consumption in the residential sector. The main problem in setting energy efficiency standard is identifying annual efficiency improvement, due to the lack of time series statistical data available in developing countries. This study attempts to present a method of calculating annual energy efficiency improvement for TV set, which can be used for implementing energy efficiency standard for TV sets in Malaysia and other developing countries. Although the presented result is only an approximation, definitely it is one of the ways of accomplishing energy standard. Furthermore, the method can be used for other appliances without any major modification.
Method for calculating annual energy efficiency improvement of TV sets
The popularization of 24 h pay-TV, interactive video games, web-TV, VCD and DVD are poised to have a large impact on overall TV electricity consumption in the Malaysia. Following this increased consumption, energy efficiency standard present a highly effective measure for decreasing electricity consumption in the residential sector. The main problem in setting energy efficiency standard is identifying annual efficiency improvement, due to the lack of time series statistical data available in developing countries. This study attempts to present a method of calculating annual energy efficiency improvement for TV set, which can be used for implementing energy efficiency standard for TV sets in Malaysia and other developing countries. Although the presented result is only an approximation, definitely it is one of the ways of accomplishing energy standard. Furthermore, the method can be used for other appliances without any major modification
A new method for calculation of an air quality index
Ilvessalo, P. [Finnish Meteorological Inst., Helsinki (Finland). Air Quality Dept.
1995-12-31
Air quality measurement programs in Finnish towns have expanded during the last few years. As a result of this it is more and more difficult to make use of all the measured concentration data. Citizens of Finnish towns are nowadays taking more of an interest in the air quality of their surroundings. The need to describe air quality in a simplified form has increased. Air quality indices permit the presentation of air quality data in such a way that prevailing conditions are more easily understandable than when using concentration data as such. Using an air quality index always means that some of the information about concentrations of contaminants in the air will be lost. How much information is possible to extract from a single index number depends on the calculation method. A new method for the calculation of an air quality index has been developed. This index always indicates the overstepping of an air quality guideline level. The calculation of this air quality index is performed using the concentrations of all the contaminants measured. The index gives information both about the prevailing air quality and also the short-term trend. It can also warn about the expected exceeding of guidelines due to one or several contaminants. The new index is especially suitable for the real-time monitoring and notification of air quality values. The behaviour of the index was studied using material from a measurement period in the spring of 1994 in Kaepylae, Helsinki. Material from a pre-operational period in the town of Oulu was also available. (author)
An accurate δf method for neoclassical transport calculation
A δf method, solving drift kinetic equation, for neoclassical transport calculation is presented in detail. It is demonstrated that valid results essentially rely on the correct evaluation of marker density g in weight calculation. A general and accurate weighting scheme is developed without using some assumed g in weight equation for advancing particle weights, unlike the previous schemes. This scheme employs an additional weight function to directly solve g from its kinetic equation using the idea of δf method. Therefore the severe constraint that the real marker distribution must be consistent with the initially assumed g during a simulation is relaxed. An improved like-particle collision scheme is presented. By performing compensation for momentum, energy and particle losses arising from numerical errors, the conservations of all the three quantities are greatly improved during collisions. Ion neoclassical transport due to self-collisions is examined under finite banana case as well as zero banana limit. A solution with zero particle and zero energy flux (in case of no temperature gradient) over whole poloidal section is obtained. With the improvement in both like-particle collision scheme and weighting scheme, the δf simulation shows a significantly upgraded performance for neoclassical transport study. (author)
Simulation of FEL pulse length calculation with THz streaking method.
Gorgisyan, I; Ischebeck, R; Prat, E; Reiche, S; Rivkin, L; Juranić, P
2016-05-01
Having accurate and comprehensive photon diagnostics for the X-ray pulses delivered by free-electron laser (FEL) facilities is of utmost importance. Along with various parameters of the photon beam (such as photon energy, beam intensity, etc.), the pulse length measurements are particularly useful both for the machine operators to measure the beam parameters and monitor the stability of the machine performance, and for the users carrying out pump-probe experiments at such facilities to better understand their measurement results. One of the most promising pulse length measurement techniques used for photon diagnostics is the THz streak camera which is capable of simultaneously measuring the lengths of the photon pulses and their arrival times with respect to the pump laser. This work presents simulations of a THz streak camera performance. The simulation procedure utilizes FEL pulses with two different photon energies in hard and soft X-ray regions, respectively. It recreates the energy spectra of the photoelectrons produced by the photon pulses and streaks them by a single-cycle THz pulse. Following the pulse-retrieval procedure of the THz streak camera, the lengths were calculated from the streaked spectra. To validate the pulse length calculation procedure, the precision and the accuracy of the method were estimated for streaking configuration corresponding to previously performed experiments. The obtained results show that for the discussed setup the method is capable of measuring FEL pulses with about a femtosecond accuracy and precision. PMID:27140142
An accurate {delta}f method for neoclassical transport calculation
Wang, W.X.; Nakajima, N.; Murakami, S.; Okamoto, M. [National Inst. for Fusion Science, Toki, Gifu (Japan)
1999-03-01
A {delta}f method, solving drift kinetic equation, for neoclassical transport calculation is presented in detail. It is demonstrated that valid results essentially rely on the correct evaluation of marker density g in weight calculation. A general and accurate weighting scheme is developed without using some assumed g in weight equation for advancing particle weights, unlike the previous schemes. This scheme employs an additional weight function to directly solve g from its kinetic equation using the idea of {delta}f method. Therefore the severe constraint that the real marker distribution must be consistent with the initially assumed g during a simulation is relaxed. An improved like-particle collision scheme is presented. By performing compensation for momentum, energy and particle losses arising from numerical errors, the conservations of all the three quantities are greatly improved during collisions. Ion neoclassical transport due to self-collisions is examined under finite banana case as well as zero banana limit. A solution with zero particle and zero energy flux (in case of no temperature gradient) over whole poloidal section is obtained. With the improvement in both like-particle collision scheme and weighting scheme, the {delta}f simulation shows a significantly upgraded performance for neoclassical transport study. (author)
An improved method for calculation of interface pressure force in PLIC-VOF methods
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)
Present HTR physics calculational methods at the Paul Scherrer Institute
In this paper a general description of HTR related calculational methods and data available at the Paul Scherrer Institute (PSI) is given. The cell codes used are MICROX-2, WIMS-D, and TRAMIX. MICROX-2 is an integral transport theory spectrum code which solves the neutron slowing down and thermalization equations on a detailed energy grid for a two region lattice cell. A second level of heterogeneity can be treated, i.e. the inner region may include two different types of grains (particles). WIMS-D contains tabulations of temperature dependent resonance integrals accurately evaluated for homogeneous mixtures of moderator and absorber at many energy points. Equivalence theorems are utilized to obtain few-group effective cross sections in heterogeneous problems. The intermediate resonance absorption shielding (IR) method of resolved resonances is employed. The spherical geometry of the pebble can be converted into an equivalent cylindrical geometry, and the double heterogeneity can be treated using the available cluster option with a method developed by Segev. TRAMIX is a flexible computer (cell) code which reads fine group nuclear data from a library in the Los Alamos National Laboratory (LANL) format MATXS and produces fine group cross sections using the IR methods. Self-shielded multigroup cross sections produced with these cell codes can be used in connection with the one- and two-dimensional discrete-ordinates finite difference transport codes ONEDANT and TWODANT from LANL for full reactor calculations. Additional interface modules which are available include a BN module for leakage calculations and some library management codes plus the PERT-V diffusion perturbation theory code. Most of the nuclear data libraries for the cell codes are based on the JEF-1 (Joint European File) evaluation. They were processed using the NJOY nuclear data processing system from LANL with additional PSI developments. These libraries include the fast, thermal and resonance data tapes
About possibilities using of theoretical calculation methods in radioecology
Full text: Increasing the radiation level into environment is accompanied by accumulation of radioactive compounds into organism and/or their migration into biosphere. Radiotoxins are accumulated into irradiated plants and animals in result of violation of exchanging processes. The are play an important role at the pathogenesis of irradiation. To date, there is well known that even small quantity of the pesticides capable intensified the radiation effect. To understand the mechanism of radiation effect on physiologically active compounds and their complexes, the knowledge of such molecules three-dimensional organization and electron structure is essential. This work is devoted to study the pesticides of carbamate range, i.e. 'sevin' and its derivatives the physiological activity of which has been connected with cholinesterase degradation. Spatial organization and conformational possibilities of the pesticides has been studied using a method of the theoretical conformational analysis on the base of computational program worked out in laboratory of Molecular Biophysics at the Baku State University. Quantum-chemical methods CNDO/2, AM1 and PM3 and complex programs 'LEV' were used in studies of electronic structures of 'sevin' and number of its analogues. Charge distribution on the atoms, optimization of geometrical electrooptic parameters, as well as molecular electrostatic potentials, electron density and nuclear forces were calculated. Visual maps and surface of valence electron density distribution in the given plane and surface of electron-nuclear forces distribution projection were constructed. The geometrical and energetic characteristics, charges on the atoms of investigated pesticides, as well as the maps and relief of the valence electron density distribution on the atoms have been received. According to calculation results, the changing of charge distribution in naphthalene ring is observed. The conclusion was made that the carbonyl group is essential for
METHOD FOR CALCULATION OF STRESSED STATE SUBSTANTIATED BY DYNAMIC MICROTWIN
V. V. Vlashevich
2015-01-01
Full Text Available Method for calculation of the stressed state in a dynamic twin has been developed on the basis of a non-thin non-coherent micro-twin model with continuous distribution of twinning dislocations at twin boundaries. In this case there is no additional generation with the help of twinning dislocation source. The model takes into account that the twin has coherent and noncoherent boundary sections. The developed model has made it possible to take into consideration a form of non-coherent sections of twinning boundaries in calculations of stressed and deformed state at dynamic twins. It has been established that localized stresses are migrating together with non-coherent sections of the twin. Normal stresses σxx change their sign in relation to direction of the twin development. Shear stresses σxy are alternating in signs in relation to an axis which is perpendicular to the direction of the twin development and which is passing through a mid-point of non-coherent twin section. Distribution of stresses σyy и σyz has similar configuration. Stresses σzx in the second and fourth quarters of XOY plane are negative and the stresses in the first and third quarters are positive. Distribution of stresses σzz practically does not differ from distribution of stresses σyy according to configuration but numerical values of stress tensor component data are different.The results have been obtained without thin twin model that permits to consider only elastic stage of the twinning process. The executed stress calculations at dynamic twin are important for forecasting at the accumulation stage of damage origination which is caused by twinning destruction and permit to improve forecasting accuracy of technical system resources on the basis of twinning materials such as alloys based on iron, copper, zinc, aluminium, titanium.
Gradient type optimization methods for electronic structure calculations
Zhang, Xin; Wen, Zaiwen; Zhou, Aihui
2013-01-01
The density functional theory (DFT) in electronic structure calculations can be formulated as either a nonlinear eigenvalue or direct minimization problem. The most widely used approach for solving the former is the so-called self-consistent field (SCF) iteration. A common observation is that the convergence of SCF is not clear theoretically while approaches with convergence guarantee for solving the latter are often not competitive to SCF numerically. In this paper, we study gradient type methods for solving the direct minimization problem by constructing new iterations along the gradient on the Stiefel manifold. Global convergence (i.e., convergence to a stationary point from any initial solution) as well as local convergence rate follows from the standard theory for optimization on manifold directly. A major computational advantage is that the computation of linear eigenvalue problems is no longer needed. The main costs of our approaches arise from the assembling of the total energy functional and its grad...
NEW METHOD FOR CALCULATION OF STATISTIC MISTAKE IN MARKETING INVESTIGATIONS
V. A. Koldachiov
2015-01-01
Full Text Available An idea of a new method is that while breaking-down analysis sample in some sub-samples there is a probability that an actual value for general body will be inside the interval between the highest and lowest average meaning of sub-sample is much higher of the probability that the given value will be beyond the limits of the indicated interval. In this case a size of the interval appears to be less than analogous parameter while making calculation with the help of the Stewdent formula.Thus, it is possible to reach high accuracy in results of marketing investigations while preserving analysis sample size or reducing the necessary size of analysis sample while preserving level of statistical mistake.
Hamiltonian lattice field theory: Computer calculations using variational methods
A variational method is developed for systematic numerical computation of physical quantities-bound state energies and scattering amplitudes-in quantum field theory. An infinite-volume, continuum theory is approximated by a theory on a finite spatial lattice, which is amenable to numerical computation. An algorithm is presented for computing approximate energy eigenvalues and eigenstates in the lattice theory and for bounding the resulting errors. It is shown how to select basis states and choose variational parameters in order to minimize errors. The algorithm is based on the Rayleigh-Ritz principle and Kato's generalizations of Temple's formula. The algorithm could be adapted to systems such as atoms and molecules. It is shown how to compute Green's functions from energy eigenvalues and eigenstates in the lattice theory, and relate these to physical (renormalized) coupling constants, bound state energies and Green's functions. Thus one can compute approximate physical quantities in a lattice theory that approximates a quantum field theory with specified physical coupling constants. The author discusses the errors in both approximations. In principle, the errors can be made arbitrarily small by increasing the size of the lattice, decreasing the lattice spacing and computing sufficiently long. Unfortunately, the author does not understand the infinite-volume and continuum limits well enough to quantify errors due to the lattice approximation. Thus the method is currently incomplete. The method is applied to real scalar field theories using a Fock basis of free particle states. All needed quantities can be calculated efficiently with this basis. The generalization to more complicated theories is straightforward. The author describes a computer implementation of the method and present numerical results for simple quantum mechanical systems
Hamiltonian lattice field theory: Computer calculations using variational methods
I develop a variational method for systematic numerical computation of physical quantities -- bound state energies and scattering amplitudes -- in quantum field theory. An infinite-volume, continuum theory is approximated by a theory on a finite spatial lattice, which is amenable to numerical computation. I present an algorithm for computing approximate energy eigenvalues and eigenstates in the lattice theory and for bounding the resulting errors. I also show how to select basis states and choose variational parameters in order to minimize errors. The algorithm is based on the Rayleigh-Ritz principle and Kato's generalizations of Temple's formula. The algorithm could be adapted to systems such as atoms and molecules. I show how to compute Green's functions from energy eigenvalues and eigenstates in the lattice theory, and relate these to physical (renormalized) coupling constants, bound state energies and Green's functions. Thus one can compute approximate physical quantities in a lattice theory that approximates a quantum field theory with specified physical coupling constants. I discuss the errors in both approximations. In principle, the errors can be made arbitrarily small by increasing the size of the lattice, decreasing the lattice spacing and computing sufficiently long. Unfortunately, I do not understand the infinite-volume and continuum limits well enough to quantify errors due to the lattice approximation. Thus the method is currently incomplete. I apply the method to real scalar field theories using a Fock basis of free particle states. All needed quantities can be calculated efficiently with this basis. The generalization to more complicated theories is straightforward. I describe a computer implementation of the method and present numerical results for simple quantum mechanical systems
Improved Numerical Method for Calculation of 4-Body Transition Amplitudes
Harris, A. L.
2013-01-01
In order to study 4-body atomic collisions such as excitation-ionization, transfer with target excitation, and double electron capture, the calculation of a nine-dimensional numerical integral is often required. This calculation can become computationally expensive, especially when calculating fully differential cross sections (FDCS), where the positions and momenta of all the particles are known. We have developed a new technique for calculating FDCS using fewer computing hours, but more mem...
Reactor calculation in coarse mesh by finite element method applied to matrix response method
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.)
CALCULATION OF MILL RIGIDITY BY THREE DIMENSION CONTACT BOUNDARY ELEMENT METHOD
无
2001-01-01
Vertical rigidity of the space self-adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three-dimension elastic contact problem,which can update the existed deforming separation calculating theory and corresponding methods of material mechanics,elastic mechanics and finite element method.The method has less hypotheses and stronger synthesis in contact-type calculating model.The advantages of the method are high calculating rate,high calculating accuracy,etc..
Dynamic Calculation Method of Beam System Under Low Velocity Impact
LI Wen-pei; WANG De-rong; SONG Chun-ming; WANG Ming-yang
2008-01-01
The beating beams and the supporting beams under low velocity impact may be in four different strain stages of deformation depending on the impact intensity and beam structure strength. Based on the different judging conditions of deformation stages, the corresponding calculation models are proposed, the calculation formulae for the determination of the impact force and the beam's lateral displacement are obtained. Calculation shows that the beam's total deflection is small when the flexibility of the supporting component is high and the effect of diminishing deflection disappears almost when the stiffness of the supporting component is high.
Measurement and Calculation Methods of a Stem Image Information
Yang Hua; Meng Xianyu; Liu Yan; Cheng Jun
2006-01-01
The paper shows a study on 2-D stem image information collected by a digital camera.Information on a single stem is obtained through calculations after the application of the direct linear transformation model of close-range photogrammetry and binocular stereo vision technology,so that the calculating problem between stem image information and its 2-D coordinate can be solved.Furthermore,the 2-D processing methodology for measuring tree image information simplifies calculating equations and increases calculating speed.Although computer stereo vision techniques for collecting parameters of a single stem shape are comparative,complicated,and expensive,research indicates the efficiently and feasibility of closerange photogrammetry for stem image information.
Estonian oil shale resources calculated by GIS method
A digital map of Estonian oil shale mining was created for joining the data about technological, environmental, and social limitations in the deposit. For evaluating potential resource of oil shale, based on borehole database, its amount, tonnage and energy were calculated. Thereafter the quantity of economical oil shale for power plants and shale oil resource were calculated. Energy rating is the most important factor for determining oil shale reserves in the case of using it for electricity generation. In the case of oil production, data on oil yield and potential resources in oil shale are the most important figures to determine the value of the deposit. Basing on the models, oil resource has been calculated. Resource data can be used for composing master plans for the deposit considering both power generation and oil production. The data can be also used for composing development plans of mines and for logistics calculations. (author)
Comparison of Methods for Calculating Radiative Heat Transfer
Schock, Alfred; Abbate, M J
2012-01-19
Various approximations for calculating radioactive heat transfer between parallel surfaces are evaluated. This is done by applying the approximations based on total emissivities to a special case of known spectral emissivities, for which exact heat transfer calculations are possible. Comparison of results indicates that the best approximation is obtained by basing the emissivity of the receiving surface primarily on the temperature of the emitter. A specific model is shown to give excellent agreement over a very wide range of values.
An approximate method for calculating the deformation of rotating nuclei
The author presents as a collective model where the potential surface at spin I=0 is calculated in the Nilsson-Strutinsky model, an analytical expression for the moment of inertia is used which depends on the deformation and the pairing gaps for protons and neutrons, and the energy is minimized with respect to these gaps. Calculations in this model are performed for 16Oyb. (HSI)
A mathematical method to calculate efficiency of BF3 detectors
SI Fenni; HU Qingyuan; PENG Taiping
2009-01-01
In order to calculate absolute efficiency of the BF3 detector, MCNP/4C code is applied to calculate rela-tive efficiency of the BF3 detector first, and then absolute efficiency is figured out through mathematical techniques. Finally an energy response curve of the BF3 detector for 1~20 MeV neutrons is derived. It turns out that efficiency of BF3 detector are relatively uniform for 2~16 MeV neutrons.
Computer development has a bearing on the choice of methods and their possible uses. The authors discuss the possible uses of the diffusion and transport theories and their limitations. Most of the problems encountered in regard to criticality involve fissile materials in simple or multiple assemblies. These entail the use of methods of calculation based on different principles. There are approximate methods of calculation, but very often, for economic reasons or with a view to practical application, a high degree of accuracy is required in determining the reactivity of the assemblies in question, and the methods based on the Monte Carlo principle are then the most valid. When these methods are used, accuracy is linked with the calculation time, so that the usefulness of the codes derives from their speed. With a view to carrying out the work in the best conditions, depending on the geometry and the nature of the materials involved, various codes must be used. Four principal codes are described, as are their variants; some typical possibilities and certain fundamental results are presented. Finally the accuracies of the various methods are compared. (author)
Benchmark calculations for hexagonal lattices with different methods
Necessity to increase the safety conditions of exploitation of recently designed core of modern nuclear reactors causes stronger requirements to the precision of neutron-physical analysis. To get more precise characteristics of nuclear reactor cells and assembly one can increase the accuracy of neutron-physical calculation analysis by taking account the spectral effects. This paper deals with the analysis of the ZR-6 series of experiments using some components of the KARATE code system. The goal of our investigations is the comparison of measured and calculated parameters of perturbed hexagonal lattices containing Gd2O3 in Al2O3 matrix or water holes/ The quoted results include: the critical y parameters Hcr, dρ/dh and the absorber rod efficiency: Δρ. The experiments are based on doubling time measurements. The calculations have been compared not only to the measured data but to the Monte Carlo code results, too (Authors)
Methods for calculating SEU rates for bipolar and NMOS circuits
McNulty, P. J.; Abdel-Kader, W. G.; Bisgrove, J. M.
1985-12-01
Computer codes developed at Clarkson for simulating charge generation by proton-induced nuclear reactions in well-defined silicon microstructures can be used to calculate SEU rates for specific devices when the critical charge and the dimensions of all SEU sensitive junctions on the device are known, provided one can estimate the contribution from externally-generated charge which enters the sensitive junction by drift and diffusion. Calculations for two important bipolar devices, the AMD 2901B bit slice and the Fairchild 93L422 RAM, for which the dimensions of the sensitive volumes were estimated from available heavy-ion test data, have been found to be in agreement with experimental data. Circuit data for the Intel 2164A, an alpha sensitive dRAM, was provided by the manufacturer. Calculations based on crude assumptions regarding which nuclear recoils and which alphas trigger upsets in the 2164A were found to agree with experimental data.
Calculation of VPP basing on functional analyzing method
无
2007-01-01
VPP can be used to deter mine the maxi mumvelocities of a sailboard at various sailing-routes,byestablishing the forces and moments balance-equa-tions on the sail and board in accordance with theprinciple of the maxi mal drive-force.Selectingroute is the most i mportant issue in upwind-sailing,and VPP calculations could provide the basis for de-ter mining the opti mal routes.VPP calculation of the sailboard perfor mance isa complex and difficult research task,and there arefew projects in this research-field...
Comparison of hardenability calculation methods of the heat-treatable constructional steels
Evaluation has been made of the consistency of calculation of the hardenability curves of the selected heat-treatable alloyed constructional steels with the experimental data. The study has been conducted basing on the analysis of present state of knowledge on hardenability calculation employing the neural network methods. Several calculation examples and comparison of the consistency of calculation methods employed are included. (author)
The correlated-k method and related methods for broadband radiation calculations
The fundamentals of radiation calculations in planetary atmospheres have been known for a long time, but the practical problems of calculating radiation over a broad spectral range both efficiently and accurately remain an area of active research. The Goody et al. paper of 1989 established a milestone in that long effort. Goody et al. highlighted attempts to formulate the problem in terms of a few absorption coefficients that could represent typically tens of thousands of coefficients in a rigorous line-by-line calculation. Here we provide a brief background, point out that the correlated-k method is a special case of a broader spectral mapping concept, and mention some new ideas that have emerged recently.
Comparison between ASHRAE and ISO thermal transmittance calculation methods
Blanusa, Petar; Goss, William P.; Roth, Hartwig;
2007-01-01
The intent of this paper is to describe and compare the two different two-dimensional frame/spacer heat transfer calculation methodologies used in North America (FRAME [EEL. The FRAMEplus Toolkit for Heat Transfer Assessment of Building Components, Version 3.0, Enermodal Engineering, Kichener, On...
Classification of methods for annual energy harvesting calculations of photovoltaic generators
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
Analysis of some splitting and roulette algorithms in shield calculations by the Monte Carlo method
Different schemes of using the splitting and roulette methods in calculation of radiation transport in nuclear facility shields by the Monte Carlo method are considered. Efficiency of the considered schemes is estimated on the example of test calculations
Calculations of NMR chemical shifts with APW-based methods
Laskowski, Robert; Blaha, Peter
2012-01-01
We present a full potential, all electron augmented plane wave (APW) implementation of first-principles calculations of NMR chemical shifts. In order to obtain the induced current we follow a perturbation approach [Pickard and Mauri, Phys. Rev. BPRBMDO1098-012110.1103/PhysRevB.63.245101 63, 245101 (2001)] and extended the common APW + local orbital (LO) basis by several LOs at higher energies. The calculated all-electron current is represented in traditional APW manner as Fourier series in the interstitial region and with a spherical harmonics representation inside the nonoverlapping atomic spheres. The current is integrated using a “pseudocharge” technique. The implementation is validated by comparison of the computed chemical shifts with some “exact” results for spherical atoms and for a set of solids and molecules with available published data.
Calculation of gamma-ray families by Monte Carlo method
Extensive Monte Carlo calculation on gamma-ray families was carried out under appropriate model parameters which are currently used in high energy cosmic ray phenomenology. Characteristics of gamma-ray families are systematically investigated by the comparison of calculated results with experimental data obtained at mountain altitudes. The main point of discussion is devoted to examine the validity of Feynman scaling in the fragmentation region of the multiple meson production. It is concluded that experimental data cannot be reproduced under the assumption of scaling law when primary cosmic rays are dominated by protons. Other possibilities on primary composition and increase of interaction cross section are also examined. These assumptions are consistent with experimental data only when we introduce intense dominance of heavy primaries in E0>1015 eV region and very strong increase of interaction cross section (say sigma varies as Esub(0)sup(0.06)) simultaneously
Score Calculation in Informatics Contests Using Multiple Criteria Decision Methods
Jurate SKUPIENE
2011-01-01
The Lithuanian Informatics Olympiad is a problem solving contest for high school students. The work of each contestant is evaluated in terms of several criteria, where each criterion is measured according to its own scale (but the same scale for each contestant). Several jury members are involved in the evaluation. This paper analyses the problem how to calculate the aggregated score for whole submission in the above mentioned situation. The chosen methodology for solving this problem is Mult...
Multi CPU clusters and calculations by molecular dynamics method
The technical characteristics of multi CPU (Central Processor Unit) clusters in Institute of Ion-Plasma and Laser Technologies AS RUz and Institute of Mathematics and Information Technologies AS RUz are described. There is detail information about cluster s architecture, installed programs and their productivity for decision of molecular dynamics tasks. Molecular dynamics program packages GROMACS, OPENMX and AutoDock-4.2.3 are described. The results of calculations using these program packages are presented. (author)
Hargreaves and other reduced-set methods for calculating evapotranspiration
Shahidian, Shakib; Serralheiro, Ricardo; Serrano, João; Teixeira, José; Haie, Naim; Santos, Francisco
2012-01-01
Globally, irrigation is the main user of fresh water, and with the growing scarcity of this essential natural resource, it is becoming increasingly important to maximize efficiency of water usage. This implies proper management of irrigation and control of application depths in order to apply water effectively according to crop needs. Daily calculation of the Reference Potential Evapotranspiration (ETo) is an important tool in determining the water needs of different crops. The United Nations...
Experiences with leak rate calculations methods for LBB application
Grebner, H.; Kastner, W.; Hoefler, A.; Maussner, G. [and others
1997-04-01
In this paper, three leak rate computer programs for the application of leak before break analysis are described and compared. The programs are compared to each other and to results of an HDR Reactor experiment and two real crack cases. The programs analyzed are PIPELEAK, FLORA, and PICEP. Generally, the different leak rate models are in agreement. To obtain reasonable agreement between measured and calculated leak rates, it was necessary to also use data from detailed crack investigations.
Cost Accounting Methods and Calculation Agricultural Products` Cost
Saule B. Spatayeva
2015-01-01
In the condition of the current market the effective manage of expenses and calculation accountancy of cost production in agriculture must be aimed to control for resources usage at any level of technology process and getting the accountancy database needed for gaining the management targets.The improving the technologies and set up aspects of business entity activity, taken place for the last decades, which caused a significant influence on condition and structure expenses but could not prov...
7 CFR 51.308 - Methods of sampling and calculation of percentages.
2010-01-01
... 7 Agriculture 2 2010-01-01 2010-01-01 false Methods of sampling and calculation of percentages. 51..., 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...
Recently developed methods in neutral-particle transport calculations: overview
It has become increasingly apparent that successful, general methods for the solution of the neutral particle transport equation involve a close connection between the spatial-discretization method used and the source-acceleration method chosen. The first form of the transport equation, angular discretization which is discrete ordinates is considered as well as spatial discretization based upon a mesh arrangement. Characteristic methods are considered briefly in the context of future, desirable developments. The ideal spatial-discretization method is described as having the following attributes: (1) positive-positive boundary data yields a positive angular flux within the mesh including its boundaries; (2) satisfies the particle balance equation over the mesh, that is, the method is conservative; (3) possesses the diffusion limit independent of spatial mesh size, that is, for a linearly isotropic flux assumption, the transport differencing reduces to a suitable diffusion equation differencing; (4) the method is unconditionally acceleratable, i.e., for each mesh size, the method is unconditionally convergent with a source iteration acceleration. It is doubtful that a single method possesses all these attributes for a general problem. Some commonly used methods are outlined and their computational performance and usefulness are compared; recommendations for future development are detailed, which include practical computational considerations
Single-Sweep Methods for Free Energy Calculations
Maragliano, Luca; Vanden-Eijnden, Eric
2007-01-01
A simple, efficient, and accurate method is proposed to map multi-dimensional free energy landscapes. The method combines the temperature-accelerated molecular dynamics (TAMD) proposed in [Maragliano & Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006)] with a variational reconstruction method using radial-basis functions for the representation of the free energy. TAMD is used to rapidly sweep through the important regions of the free energy landscape and compute the gradient of the free energ...
Methods of parameters calculation for wastes hightemperature plasma processing
Kolesnyk, Vasyl; Orlyk, Volodymyr
2015-01-01
The article presents the methodology of plasma parameters calculation for different ratios of its constituents and taking into account that plasma jet is both a source of thermal energy and a part of initial reaction mixture for various physical and chemical transformations, in particular, those associated with processing of different solid wastes. As the methodology application example carbon conversion by steam-and-air plasma jet is investigated. У роботі запропоновано методику розрах...
Cost Accounting Methods and Calculation Agricultural Products` Cost
Saule B. Spatayeva
2015-04-01
Full Text Available In the condition of the current market the effective manage of expenses and calculation accountancy of cost production in agriculture must be aimed to control for resources usage at any level of technology process and getting the accountancy database needed for gaining the management targets.The improving the technologies and set up aspects of business entity activity, taken place for the last decades, which caused a significant influence on condition and structure expenses but could not provide the increase of economic effectiveness in agriculture.
Projection methods for the calculation of incompressible or dilatable flows
This thesis treats of time resolution methods for the Navier-Stokes equations. Based on the well-known projection method of Chorin and Temam, an original pressure correction method, named 'projection-penalty' is developed. Its specificity concerns the addition of a penalty term in the prediction step, which constrains the predicted velocity to fit with the mass balance. The precision improvements added by this method are demonstrated by some analysis results and by some numerical experiments of incompressible or dilatable flows. Finally, the potentialities offered by the use of the joint finite elements method in this type of fractionary step scheme is studied. Two applications are presented, one for local refinement purpose, the other for the resolution of a multi-physics problem. (J.S.)
A method for the calculation of annual limits on intake
The German Radiation Protection Ordinance issued in 1989 - as well as ICRP Publication 30 - included tables of annual limits on intake (ALIs). These are secondary limits derived to indicate those intake values which are not allowed to be exceeded within a year to meet the primary annual dose limits. The latest ICRP Publications dealing with doses for workers give no ALIs any more and also the IAEA Basic Safety Standards and the Council Directive of the European Commission do not list such values. And also the new Radiation Protection Ordinance as well as the federal gazette with dose coefficients do not include such values any more. However, it may be that they will be included as guiding levels in the calculation principles for the assessment of internal doses. This paper will give reasons why there are no ALIs given any more, it will show how ALIs could be calculated, it will show the influence of additional dose restrictions for female workers of child-bearing age, and it will discuss if radiation protection for the unborn child would be fulfilled if the dose limits of the radiation protection ordinance are met. (orig.)
Godunov Method for Calculating Multicomponent Heterogeneous Medium Flows
Surov, V. S.
2014-03-01
The modified Godunov method intended for integrating the nondivergent systems that describe a multivelocity heterogeneous mixture flow is presented. The linearized Riemann solver has been used in solving the Riemann problems.
A nonlinear analytic function expansion nodal method for transient calculations
Joo, Han Gyn; Park, Sang Yoon; Cho, Byung Oh; Zee, Sung Quun [Korea Atomic Energy Research Institute, Taejon (Korea, Republic of)
1998-12-31
The nonlinear analytic function expansion nodal (AFEN) method is applied to the solution of the time-dependent neutron diffusion equation. Since the AFEN method requires both the particular solution and the homogeneous solution to the transient fixed source problem, the derivation of the solution method is focused on finding the particular solution efficiently. To avoid complicated particular solutions, the source distribution is approximated by quadratic polynomials and the transient source is constructed such that the error due to the quadratic approximation is minimized, In addition, this paper presents a new two-node solution scheme that is derived by imposing the constraint of current continuity at the interface corner points. The method is verified through a series of application to the NEACRP PWR rod ejection benchmark problems. 6 refs., 2 figs., 1 tab. (Author)
A Green's function nodal expansion method for LWR diffusion calculation
A Green's Function Nodal Expansion Method (GNEM) has been developed for the efficient numerical solution of the LWR multi-dimensional neutron diffusion equation. It is an improved version of Nodal Expansion Method (NEM) and Nodal Green's Function Method (NGFM). The code interior fluxes are approximated by a high order polynomial expansion as in NEM. The nodal surface fluxes are coupled with the net currents by using the Green's function method to improve the accuracy. A computer code GNEM has been developed and tested. The numerical results demonstrate that GNEM has the same accuracy as NGFM, while it is twice as fast as NGFM. Especially, the numerical results of TMI-1 core depletion cycles 1 and 6 demonstrate that GNEM is about two times faster than ADMARC and possesses better accuracy
Hourly Calculation Method of Air Source Heat Pump Behavior
Ludovico Danza; Lorenzo Belussi; Italo Meroni; Michele Mililli; Francesco Salamone
2016-01-01
The paper describes an hourly simplified model for the evaluation of the energy performance of heat pumps in cooling mode maintaining a high accuracy and low computational cost. This approach differs from the methods used for the assessment of the overall energy consumption of the building, normally placed in the so-called white or black box models, where the transient conduction equation is deterministically and stochastically solved, respectively. The present method wants to be the expressi...
Uranyl minerals form by oxidation and alteration of uraninite, UO2+x, and the UO2 in used nuclear fuels. The thermodynamic database for these phases is extremely limited. However, the Gibbs free energies and enthalpies for uranyl phases may be estimated based on a method that sums polyhedral contributions. The molar contributions of the structural components to Δf Gm0 and Δf Hm0 are derived by multiple regression using the thermodynamic data of phases for which the crystal structures are known. In comparison with experimentally determined values, the average residuals associated with the predicted Δf Gm0 and Δf Hm0 for the uranyl phases used in the model are 0.08 and 0.10%, respectively. There is also good agreement between the predicted mineral stability relations and field occurrences, thus providing confidence in this method for the estimation of Δf Gm0 and Δf Hm0 of the U(VI) phases. This approach provides a means of generating estimated thermodynamic data for performance assessment calcination and a basic for making bounding calcination of phase stabilities and solubilities. (author)
A New Pseudospectral Method for Calculations of Hydrogen Atom in Arbitrary External Fields
QIAO Hao-Xue; LI Bai-Wen1
2002-01-01
A new pseudospectral method was introduced to calculate wavefunctions and energy levels of hydrogen atom in arbitrary potential. Some results of hydrogen atom in uniform magnetic fields were presented, high accuracy of results was obtained with simple calculations, and our calculations show very fast convergence. It suggests a new methodfor calculations of hydrogen atom in external fields.
[Calculating method for crop water requirement based on air temperature].
Tao, Guo-Tong; Wang, Jing-Lei; Nan, Ji-Qin; Gao, Yang; Chen, Zhi-Fang; Song, Ni
2014-07-01
The importance of accurately estimating crop water requirement for irrigation forecast and agricultural water management has been widely recognized. Although it has been broadly adopted to determine crop evapotranspiration (ETc) via meteorological data and crop coefficient, most of the data in whether forecast are qualitative rather than quantitative except air temperature. Therefore, in this study, how to estimate ETc precisely only using air temperature data in forecast was explored, the accuracy of estimation based on different time scales was also investigated, which was believed to be beneficial to local irrigation forecast as well as optimal management of water and soil resources. Three parameters of Hargreaves equation and two parameters of McClound equation were corrected by using meteorological data of Xinxiang from 1970 to 2010, and Hargreaves equation was selected to calculate reference evapotranspiration (ET0) during the growth period of winter wheat. A model of calculating crop water requirement was developed to predict ETc at time scales of 1, 3, and 7 d intervals through combining Hargreaves equation and crop coefficient model based on air temperature. Results showed that the correlation coefficients between measured and predicted values of ETc reached 0.883 (1 d), 0.933 (3 d), and 0.959 (7 d), respectively. The consistency indexes were 0.94, 0.95 and 0.97, respectively, which showed that forecast error decreased with the increasing time scales. Forecasted accuracy with an error less than 1 mm x d(-1) was more than 80%, and that less than 2 mm x d(-1) was greater than 90%. This study provided sound basis for irrigation forecast and agricultural management in irrigated areas since the forecasted accuracy at each time scale was relatively high. PMID:25345053
The paper aims to present the main physical principles for selection of design characteristics of the fast reactor control rods (CR) system. The brief analysis of problems of CR physical calculations is given. Four components are described for the correction to the control rod worth calculated by the routine method based on the few - group three - dimensional diffusion code (TRIGEX) in hexagonal geometry. Principle considerations are given for the choice of the original task discretization methods implemented in this code to minimize the total error. Brief information is given about methods and codes used for the evaluation of error components of control rod worths calculated in a standard way. The results of experimental and calculational investigations of control rod physical characteristics are presented. These results were obtained at BFS critical assemblies simulating LMFBR cores. The investigations have been carried out for different types of core configurations. The experimental and calculated values are given on the distortion of power distribution due to the control rod insertion in the core. (author). 51 refs, 9 figs, 5 tabs
Many methods are used to calculate the positron lifetime, these methods could be divided into two main types. The first method is atomic superposition approximation method and the second one is the so called energy band calculation method. They are also known as the non-self-consistent field method and self-consistent field method respectively. In this paper, we first introduce the two basic methods and then, we take Si as an example and give our calculation results, these results coincide with our latest experimental results, finally, we discuss the advantages and disadvantages of the two methods
Transport survey calculations using the spectral collocation method
A novel transport survey code has been developed and is being used to study the sensitivity of stellarator reactor performance to various transport assumptions. Instead of following one of the usual approaches, the steady-state transport equation are solved in integral form using the spectral collocation method. This approach effectively combine the computational efficiency of global models with the general nature of 1-D solutions. A compact torsatron reactor test case was used to study the convergence properties and flexibility of the new method. The heat transport model combined Shaing's model for ripple-induced neoclassical transport, the Chang-Hinton model for axisymmetric neoclassical transport, and neoalcator scaling for anomalous electron heat flux. Alpha particle heating, radiation losses, classical electron-ion heat flow, and external heating were included. For the test problem, the method exhibited some remarkable convergence properties. As the number of basis functions was increased, the maximum, pointwise error in the integrated power balance decayed exponentially until the numerical noise level as reached. Better than 10% accuracy in the globally-averaged quantities was achieved with only 5 basis functions; better than 1% accuracy was achieved with 10 basis functions. The numerical method was also found to be very general. Extreme temperature gradients at the plasma edge which sometimes arise from the neoclassical models and are difficult to resolve with finite-difference methods were easily resolved. 8 refs., 6 figs
Comparison of methods for numerical calculation of continuum damping
Bowden, George; Hole, Matthew; Gorelenkov, Nikolai; Dennis, Graham
2014-01-01
Continuum resonance damping is an important factor in determining the stability of certain global modes in fusion plasmas. A number of analytic and numerical approaches have been developed to compute this damping, particularly in the case of the toroidicity-induced shear Alfv\\'en eigenmode. This paper compares results obtained using an analytical perturbative approach with those found using resistive and complex contour numerical approaches. It is found that the perturbative method does not provide accurate agreement with reliable numerical methods for the range of parameters examined. This discrepancy exists even in the limit where damping approaches zero. When the perturbative technique is implemented using a standard finite element method, the damping estimate fails to converge with radial grid resolution. The finite elements used cannot accurately represent the eigenmode in the region of the continuum resonance, regardless of the number of radial grid points used.
Analytical method of spectra calculations in the Bargmann representation
We formulate a universal method for solving an arbitrary quantum system which, in the Bargmann representation, is described by a system of linear equations with one independent variable, such as one- and multi-photon Rabi models, or N level systems interacting with a single mode of the electromagnetic field and their various generalizations. We explain three types of conditions that determine the spectrum and show their usage for two deformations of the Rabi model. We prove that the spectra of both models are just zeros of transcendental functions, which in one case are given explicitly in terms of confluent Heun functions. - Highlights: • Analytical method of spectrum determination in Bargmann representation is proposed. • Three types of conditions determining spectrum are identified. • Method to two generalizations of the Rabi system is applied
Three-dimensional calculations using the quiet implicit PIC method
Solution of the time-implicit moment equations for electron and ion species, combined with Maxwell's equations, is the kernel of the moment-implicit particle method and of the quiet implicit PIC method, a generalized δf method for electromagnetic plasma simulation. Previous implementations have used a variety of direct and iterative approaches for obtaining solutions to large sparse linear systems. These methods suffer from both excessive computational cost, sometimes negating the advantage of implicit time differencing; and from lack of convergence in some regimes of interest, rendering the method inapplicable. We describe a new formulation of the coupled problem, leading to a symmetric, positive definite system. We also show that this symmetric problem may be efficiently and reliably solved by a conjugate gradient method. A three-dimensional algorithm has been constructed, using a pseudospectral Fourier treatment of the poloidal and toroidal directions, and a finite difference treatment of the radial direction. The radially-dependent, poloidal and toroidal averaged operator is used as a preconditioner. Convergence is rapid, with a typical iteration count of 10 for 10-5 convergence. New results of the two-fluid form of this code to an internal m = 1 internal kink mode will be presented. The algorithm's ability to reproduce kinetic properties of plasmas is being tested with a one-dimensional code, which has recently been modified to permit periodic boundary conditions. When a nonuniform temperature distribution is imposed as an initial condition, the gradients axe reduced by long mean-free-path particles that stream parallel to the magnetic field. Tests of collisionless wave damping have also been conducted for the ion-cyclotron range of frequencies, and results axe compared with analytic predictions
Ab initio calculations of mechanical properties: Methods and applications
Pokluda, J.; Černý, Miroslav; Šob, Mojmír; Umeno, Y.
2015-01-01
Roč. 73, AUG (2015), s. 127-158. ISSN 0079-6425 R&D Projects: GA ČR(CZ) GAP108/12/0311 Institutional support: RVO:68081723 Keywords : Ab initio methods * Elastic moduli * Intrinsic hardness * Stability analysis * Theoretical strength * Intrinsic brittleness/ductility Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 27.417, year: 2014
Hourly Calculation Method of Air Source Heat Pump Behavior
Ludovico Danza
2016-04-01
Full Text Available The paper describes an hourly simplified model for the evaluation of the energy performance of heat pumps in cooling mode maintaining a high accuracy and low computational cost. This approach differs from the methods used for the assessment of the overall energy consumption of the building, normally placed in the so-called white or black box models, where the transient conduction equation is deterministically and stochastically solved, respectively. The present method wants to be the expression of the grey box model, taking place between the previous approaches. The building envelope is defined using a building thermal model realized with a 3 Resistance 1 Capacitance (3R1C thermal network based on the solution of the lumped capacitance method. The simplified model evaluates the energy efficiency ratio (EER of a heat pump through the determination of the hourly second law efficiency of a reversed Carnot cycle. The results of the simplified method were finally compared with those provided by EnergyPlus, a dynamic building energy simulation program, and those collected from an outdoor test cell in real working conditions. The results are presented in temperatures and energy consumptions profiles and are validated using the Bland-Altman test.
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)
A Direct Iteration Method using Resonance Integral Table for the Self-Shielding Calculations
In this paper, a direct iteration method using the resonance integral table is introduced for the self-shielding calculations. The basic purpose of this paper is to show the possibility that the HELIOS subgroup method can be replaced with this method. This method doesn't use the subgroup data but only the resonance integral tables given in library. The basic idea of this method is to use the Bondarenko's iteration in order to obtain the self-shielded effective cross sections with the background cross sections which are calculated by the heterogeneous transport calculation. This method is implemented in the KARMA lattice calculation code and tested
A Direct Iteration Method using Resonance Integral Table for the Self-Shielding Calculations
Hong, Ser Gi; Kim, Kang Seog; Song, Jae Seung [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)
2009-10-15
In this paper, a direct iteration method using the resonance integral table is introduced for the self-shielding calculations. The basic purpose of this paper is to show the possibility that the HELIOS subgroup method can be replaced with this method. This method doesn't use the subgroup data but only the resonance integral tables given in library. The basic idea of this method is to use the Bondarenko's iteration in order to obtain the self-shielded effective cross sections with the background cross sections which are calculated by the heterogeneous transport calculation. This method is implemented in the KARMA lattice calculation code and tested.
Comparison of methods for numerical calculation of continuum damping
Bowden, George; Könies, Axel; Hole, Matthew; Gorelenkov, Nikolai; Dennis, Graham
2014-01-01
Continuum resonance damping is an important factor in determining the stability of certain global modes in fusion plasmas. A number of analytic and numerical approaches have been developed to compute this damping, particularly in the case of the toroidicity-induced shear Alfv\\'en eigenmode. This paper compares results obtained using an analytical perturbative approach with those found using resistive and complex contour numerical approaches. It is found that the perturbative method does not p...
Relief valve discharge piping calculational methods and results
The analysis of a safety relief valve and its connecting discharge pipe is an important parameter in the design of a light water reactor plant. This paper addresses the computer modelling techniques and methods used for this analysis and presents guidelines, cautions and improvements which should be used for design. Specific discussions include considerations for loop seal design, heat transfer effects, pipe submergence, reflood effects and vacuum breaker design. (orig.)
Free-energy calculation methods for collective phenomena in membranes
Smirnova, Yuliya G.; Fuhrmans, Marc; Barragan Vidal, Israel A.; Müller, Marcus
2015-09-01
Collective phenomena in membranes are those which involve the co-operative reorganization of many molecules. Examples of these are membrane fusion, pore formation, bending, adhesion or fission. The time and length scales, on which these processes occur, pose a challenge for atomistic simulations. Therefore, in order to solve the length scale problem it is popular to introduce a coarse-grained representation. To facilitate sampling of the relevant states additional computational techniques, which encourage the system to explore the free-energy landscape far from equilibrium and visit transition states, are needed. These computational techniques provide insights about the free-energy changes involved in collective transformations of membranes, yielding information about the rate limiting states, the transformation mechanism and the influence of architectural, compositional and interaction parameters. A common approach is to identify an order parameter (or reaction coordinate), which characterizes the pathway of membrane reorganization. However, no general strategy exists to define such an order parameter that can properly describe cooperative reorganizations in membranes. Recently developed methods can overcome this problem of the order-parameter choice and allow us to study collective phenomena in membranes. We will discuss such methods as thermodynamic integration, umbrella sampling, and the string method and results provided by their applications to particle-based simulations, particularly focusing on membrane fusion and pore formation.
Free-energy calculation methods for collective phenomena in membranes
Collective phenomena in membranes are those which involve the co-operative reorganization of many molecules. Examples of these are membrane fusion, pore formation, bending, adhesion or fission. The time and length scales, on which these processes occur, pose a challenge for atomistic simulations. Therefore, in order to solve the length scale problem it is popular to introduce a coarse-grained representation. To facilitate sampling of the relevant states additional computational techniques, which encourage the system to explore the free-energy landscape far from equilibrium and visit transition states, are needed. These computational techniques provide insights about the free-energy changes involved in collective transformations of membranes, yielding information about the rate limiting states, the transformation mechanism and the influence of architectural, compositional and interaction parameters. A common approach is to identify an order parameter (or reaction coordinate), which characterizes the pathway of membrane reorganization. However, no general strategy exists to define such an order parameter that can properly describe cooperative reorganizations in membranes. Recently developed methods can overcome this problem of the order-parameter choice and allow us to study collective phenomena in membranes. We will discuss such methods as thermodynamic integration, umbrella sampling, and the string method and results provided by their applications to particle-based simulations, particularly focusing on membrane fusion and pore formation. (topical review)
Rate Constant Calculation for Thermal Reactions Methods and Applications
DaCosta, Herbert
2011-01-01
Providing an overview of the latest computational approaches to estimate rate constants for thermal reactions, this book addresses the theories behind various first-principle and approximation methods that have emerged in the last twenty years with validation examples. It presents in-depth applications of those theories to a wide range of basic and applied research areas. When doing modeling and simulation of chemical reactions (as in many other cases), one often has to compromise between higher-accuracy/higher-precision approaches (which are usually time-consuming) and approximate/lower-preci
Visual Method for Spectral Energy Distribution Calculation of Blazars
Y. Huang; J. H. Fan
2014-09-01
In this work, we propose to use `The Geometer’s Sketchpad’ to the fitting of a spectral energy distribution of blazar based on three effective spectral indices, RO, OX, and RX and the flux density in the radio band. It can make us to see the fitting in detail with both the peak frequency and peak luminosity given immediately. We used our method to those sources whose peak frequency and peak luminosity are given and found that our results are consistent with those given in the work of Sambruna et al. (1996).
New nuclear medicine method of calculating left ventricular stroke volume
Gieschke, R.; Luig, H.; Reuter, R.; Figulla, H.R.
1983-12-01
A new non-invasive nuclear medicine procedure for determining the left ventricular stroke volume is described. The procedure exhibits the following features: 1. individual calibration of scintigraphic counts in activity by first-pass evaluation; 2. no need for a delta-shaped bolus injection; and 3. determination of different stroke volumes, e.g. during different grades of exercise, by only one injection and by only one blood sample. 36 results obtained at rest and during exercise are compared with corresponding results of the thermodilution method (r = 0.86).
Study on increasing calculation precision and convergence speed of streamline strip element method
彭艳; 刘宏民
2004-01-01
The calculation precision and convergence speed of streamline strip element method are increased by using the method whose initial value of the exit lateral displacement is determined with strip element variation method, and the accurate tension lateral distribution model is adopted based on the original third power spline function streamline strip element method. The basic theory of the strip element method is developed. The calculated results by the improved streamline strip element method and the original streamline strip element method are compared with the measured results, showing that the calculated results of the improved method are in good agreement with the measured results.
Comparison of hardenability calculation methods of the heat-treatable constructional steels
Dobrzanski, L.A.; Sitek, W. [Division of Tool Materials and Computer Techniques in Metal Science, Silesian Technical University, Gliwice (Poland)
1995-12-31
Evaluation has been made of the consistency of calculation of the hardenability curves of the selected heat-treatable alloyed constructional steels with the experimental data. The study has been conducted basing on the analysis of present state of knowledge on hardenability calculation employing the neural network methods. Several calculation examples and comparison of the consistency of calculation methods employed are included. (author). 35 refs, 2 figs, 3 tabs.
Finite element method in density functional theory electronic structure calculations
Vackář, Jiří; Čertík, Ondřej; Cimrman, R.; Novák, M.; Šipr, Ondřej; Plešek, Jiří
Berlin : Springer, 2012 - (Hoggan, P.; Brändas, E.; Maruani, J.; Piecuch, P.; Delgado- Barrio , G.), s. 199-217 ISBN 978-94-007-2075-6. - (Progress in Theoretical Chemistry and Physics. vol. 12) R&D Projects: GA ČR GA101/09/1630; GA ČR(CZ) GAP108/11/0853; GA MŠk(CZ) LC06040 Institutional research plan: CEZ:AV0Z10100520; CEZ:AV0Z10100521; CEZ:AV0Z20760514 Keywords : finite-element method * pseudopotentials * density functional theory Subject RIV: BM - Solid Matter Physics ; Magnetism http://www.springerlink.com/content/p7k3j7047720270r/
A numerical method for calculating resonant-state wave functions
An initial-value method of numerical solving of Sturm-Liouville problems is applied to find the solution to the Schroedinger equation which corresponds to a resonance situation. The depth of the nuclear potential is regarded as an eigenvalue, which is obtained by iteration. Having established the nuclear potential, the resonant wavefunction is generated by integrating numerically the Schroedinger differential equation inwards from larger radii using the initial conditions of G(r), where G is the irregular Coulomb function. Because the solution is exactly on resonance, nosearching for the phase shift is required. Consequently, the suggested procedure may be employed even if the resonance widths are extremely narrow (e.g., 10-16 MeV)
Study on the Processing Method for Resonance Self-shielding Calculations
无
2011-01-01
We investigate a new approach for resonance self-shielding calculations, based on a straightforward subgroup method, used in association with characteristics method. Subgroup method is actually the subdivision of cross section range for resonance energy range.
A high-precision calculation method for interface normal and curvature on an unstructured grid
Ito, Kei; Kunugi, Tomoaki; Ohno, Shuji; Kamide, Hideki; Ohshima, Hiroyuki
2014-09-01
In the volume-of-fluid algorithm, the calculations of the interface normal and curvature are crucially important for accurately simulating interfacial flows. However, few methods have been proposed for the high-precision interface calculation on an unstructured grid. In this paper, the authors develop a height function method that works appropriately on an unstructured grid. In the process, the definition of the height function is discussed, and the high-precision calculation method of the interface normal is developed to meet the necessary condition for a second-order method. This new method has highly reduced computational cost compared with a conventional high-precision method because the interface normal calculation is completed by solving relatively simple algebraic equations. The curvature calculation method is also discussed and the approximated quadric curve of an interface is employed to calculate the curvature. Following a basic verification, the developed height function method is shown to successfully provide superior calculation accuracy and highly reduced computational cost compared with conventional calculation methods in terms of the interface normal and curvature. In addition, the height function method succeeds in calculating accurately the slotted-disk revolution problem and the oscillating drop on unstructured grids. Therefore, the developed height function method is confirmed to be an efficient technique for the high-precision numerical simulation of interfacial flows on an unstructured grid.
Theories and calculation methods for regional objective ET
QIN DaYong; LO JinYan; LIU JiaHong; WANG MingNa
2009-01-01
The regional objective ET (Evapotranspiration) is a new concept in water resources research, which refers to the total amount of water that could be exhausted from a region in the form of vapor per year. The objective-ET based water resources management allocates water to different regions in terms of ET. It controls the water exhausted from a region to meet the objective ET. The regional objective ET must be adapted to fit the region's local available water resources. By improving the water utilization effi-ciency and reducing the unrecoverable water in the social water circle, it is saved so that water related production is maintained or even increased under the same water consumption conditions. Regional water balance is realized by rationally deploying the available water among different industries, adjusting industrial structures, and adopting new water-saving technologies, therefore to meeting the requirements for groundwater conservation, agricultural income stability, and avoiding environmental damages. Furthermore, water competition among various departments and industries (including envi-ronmental and ecological water use) may be avoided. This paper proposes an innovative definition of objective ET, and its principles, sub-index systems. Besides, a computational method for regional ob-jective ET is developed by combining the distributed hydrological model and the soil moisture model.
Comparison of organ dosimetry methods and effective dose calculation methods for paediatric CT
Computed tomography (CT) is the single biggest ionising radiation risk from anthropogenic exposure. Reducing unnecessary carcinogenic risks from this source requires the determination of organ and tissue absorbed doses to estimate detrimental stochastic effects. In addition, effective dose can be used to assess comparative risk between exposure situations and facilitate dose reduction through optimisation. Children are at the highest risk from radiation induced carcinogenesis and therefore dosimetry for paediatric CT recipients is essential in addressing the ionising radiation health risks of CT scanning. However, there is no well-defined method in the clinical environment for routinely and reliably performing paediatric CT organ dosimetry and there are numerous methods utilised for estimating paediatric CT effective dose. Therefore, in this study, eleven computational methods for organ dosimetry and/or effective dose calculation were investigated and compared with absorbed doses measured using thermoluminescent dosemeters placed in a physical anthropomorphic phantom representing a 10 year old child. Three common clinical paediatric CT protocols including brain, chest and abdomen/pelvis examinations were evaluated. Overall, computed absorbed doses to organs and tissues fully and directly irradiated demonstrated better agreement (within approximately 50 %) with the measured absorbed doses than absorbed doses to distributed organs or to those located on the periphery of the scan volume, which showed up to a 15-fold dose variation. The disparities predominantly arose from differences in the phantoms used. While the ability to estimate CT dose is essential for risk assessment and radiation protection, identifying a simple, practical dosimetry method remains challenging.
Preconditioned Conjugate Gradient methods for low speed flow calculations
Ajmani, Kumud; Ng, Wing-Fai; Liou, Meng-Sing
1993-01-01
An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations are integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and the convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the lower-upper (L-U)-successive symmetric over-relaxation iterative scheme is more efficient than a preconditioner based on incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional line Gauss-Seidel relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.
Odru P.
2006-11-01
Full Text Available Cet article présente une méthode de calcul des structures composites fibres-résine appliquée aux cas des tubes minces. Outre l'établissement des relations contraintes - déformations généralisées des tubes à partir des caractéristiques des matériaux de base et de leur orientation, on pose les relations permettant de calculer leur comportement et leur dimensionnement sous des charges axisymétriques combinées de traction, pression et flexion. Une méthode simplifiée applicable au cas des composites microfissurés est aussi présentée. On montre ensuite, à travers quelques exemples concrets d'applications, les propriétés intéressantes ou inhabituelles que le matériau permet de conférer aux structures. This article presents a method of calculation of composite structures applied to thin-walled tubes. Starting from the characteristics and orientation of the basic materials, the generalized stress-strain equations of the tubes are determined ; then the relationship allowing the calculation of their design and behavior under combined axisymmetrical loads of tension, pressure and bending are established. A simplified method applicable to microcracked composite materials is also described. Several complete examples of applications illustrate the interesting or unusual properties that this material can impart to structures
Lesniak, Joseph; Behrman, Elizabeth; Zandler, Melvin; Kumar, Preethika
2008-03-01
Very few quantum algorithms are currently useable today. When calculating molecular energies, using a quantum algorithm takes advantage of the quantum nature of the algorithm and calculation. A few small molecules have been used to show that this method is possible. This method will be applied to larger molecules and compared to classical computer methods.
Zhang, Wen-juan; Huang, Shou-dao; Chen, Zhe
2013-01-01
An analytic electromagnetic calculation method for doubly fed induction generator (DFIG) in wind turbine system was presented. Based on the operation principles, steady state equivalent circuit and basic equations of DFIG, the modeling for electromagnetic calculation of DFIG was proposed. The...... electromagnetic calculation of DFIG was divided into three steps: the magnetic flux calculation, parameters derivation and performance checks. For each step, the detailed numeric calculation formulas were all derived. Combining the calculation formulas, the whole electromagnetic calculation procedure was...... established, which consisted of three iterative calculation loops, including magnetic saturation coefficient, electromotive force and total output power. All of the electromagnetic and performance data of DIFG can be calculated conveniently by the established calculation procedure, which can be used to...
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
Perturbation method for calculation of narrow-band impedance and trapped modes
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
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.)
Perfetti, Christopher M [ORNL; Rearden, Bradley T [ORNL
2014-01-01
This work introduces a new approach for calculating sensitivity coefficients for generalized neutronic responses to nuclear data uncertainties using continuous-energy Monte Carlo methods. The approach presented in this paper, known as the GEAR-MC method, allows for the calculation of generalized sensitivity coefficients for multiple responses in a single Monte Carlo calculation with no nuclear data perturbations or knowledge of nuclear covariance data. The theory behind the GEAR-MC method is presented here, and proof of principle is demonstrated by using the GEAR-MC method to calculate sensitivity coefficients for responses in several 3D, continuous-energy Monte Carlo applications.
THE PARTICULARITIES OF THE COST CALCULATION METHOD ON COMMANDS IN FURNITURE INDUSTRY
Felicia Sabou
2014-01-01
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 pape...
Perfetti, Christopher M [ORNL; Martin, William R [University of Michigan; Rearden, Bradley T [ORNL; Williams, Mark L [ORNL
2012-01-01
Three methods for calculating continuous-energy eigenvalue sensitivity coefficients were developed and implemented into the SHIFT Monte Carlo code within the Scale code package. The methods were used for several simple test problems and were evaluated in terms of speed, accuracy, efficiency, and memory requirements. A promising new method for calculating eigenvalue sensitivity coefficients, known as the CLUTCH method, was developed and produced accurate sensitivity coefficients with figures of merit that were several orders of magnitude larger than those from existing methods.
A new method for calculating the activity of stable compound from binary phase diagram
CHEN Dengfu; DONG Lingyan; BAI Chenguang; LIU Qingcai; WANG Chuanjun
2006-01-01
A new method to calculate the activity of a stable compound in a binary phase diagram was presented and dis cussed. According to the formula for calculating activity from the binary phase diagram, the equilibrium constant can be calculated through the mass action principle after the activities of two pure components were computed respectively. Based on that, the activity of a stable compound can be easily obtained at last. The activity of the stable compound InSb is calculated in the In-Sb binary system by using this method. The result is well consistent with another calculation value.
Shan Yang; Xiangqian Tong
2016-01-01
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 inverte...
Comparison of the methods for calculating the interfacial heat transfer coefficient in hot stamping
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
JM-ECS: A hybrid method combining the $J$-matrix and ECS methods for scattering calculations
Bidasyuk, Y; Broeckhove, J; Arickx, F; Vasilevsky, V
2010-01-01
The paper proposes a hybrid method for calculating scattering processes that combines the $J$-matrix method with exterior complex scaling as an absorbing boundary condition. It represents the wave function as a finite sum of oscillator eigenstates in the inner region and on grid in the outer region. The method is validated for one and two dimensional model partial wave equation equations. Finally, the method calculates nuclear $p$-shell scattering.
The shell model Monte Carlo (SMMC) method enables calculations in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods, and is particularly suitable for the calculation of level densities in the presence of correlations. We review recent advances and applications of SMMC for the microscopic calculation of level densities. Recent developments include (1) a method to calculate accurately the ground-state energy of an odd-mass nucleus, circumventing a sign problem that originates in the projection on an odd number of particles, and (2) a method to calculate directly level densities, which, unlike state densities, do not include the spin degeneracy of the levels. We calculated the level densities of a family of nickel isotopes 59-64Ni and of a heavy deformed rare-earth nucleus 162Dy and found them to be in close agreement with various experimental data sets. (author)
Alhassid, Y; Liu, S; Mukherjee, A; Nakada, H
2014-01-01
The shell model Monte Carlo (SMMC) method enables calculations in model spaces that are many orders of magnitude larger than those that can be treated by conventional methods, and is particularly suitable for the calculation of level densities in the presence of correlations. We review recent advances and applications of SMMC for the microscopic calculation of level densities. Recent developments include (i) a method to calculate accurately the ground-state energy of an odd-mass nucleus, circumventing a sign problem that originates in the projection on an odd number of particles, and (ii) a method to calculate directly level densities, which, unlike state densities, do not include the spin degeneracy of the levels. We calculated the level densities of a family of nickel isotopes $^{59-64}$Ni and of a heavy deformed rare-earth nucleus $^{162}$Dy and found them to be in close agreement with various experimental data sets.
Efficient method for the calculation of dissipative quantum transport in quantum cascade lasers.
Greck, Peter; Birner, Stefan; Huber, Bernhard; Vogl, Peter
2015-03-01
We present a novel and very efficient method for calculating quantum transport in quantum cascade lasers (QCLs). It follows the nonequilibrium Green's function (NEGF) framework but sidesteps the calculation of lesser self-energies by replacing them by a quasi-equilibrium expression. This method generalizes the phenomenological Büttiker probe model by taking into account individual scattering mechanisms. It is orders of magnitude more efficient than a fully self-consistent NEGF calculation for realistic devices. We apply this method to a new THz QCL design which works up to 250 K - according to our calculations. PMID:25836876
This work introduces a new approach for calculating the sensitivity of generalized neutronic responses to nuclear data uncertainties using continuous-energy Monte Carlo methods. The GEneralized Adjoint Responses in Monte Carlo (GEAR-MC) method has enabled the calculation of high resolution sensitivity coefficients for multiple, generalized neutronic responses in a single Monte Carlo calculation with no nuclear data perturbations or knowledge of nuclear covariance data. The theory behind the GEAR-MC method is presented here and proof of principle is demonstrated by calculating sensitivity coefficients for responses in several 3D, continuous-energy Monte Carlo applications. (author)