Search NASA⌕ Search

SEARCH · Search NASA

Results for “FEYNMAN METHOD”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

COUPLING SJÖSTRAND AND FEYNMAN METHODS IN PROMPT NEUTRON DECAY CONSTANT ANALYSES

The Sjöstrand and Feynman methods have been widely applied to subcritical assemblies analyses. Within the point kinetics framework, the Sjöstrand method determines the system reactivity by the ratio between the prompt and delayed areas from pulsed neutron source experiments. In addition, the slope of the prompt area gives the prompt neutron decay constant. The latter parameter can be also obtained from fitting the Feynman curve. Consequently, the goal of this work is to combine the Sjöstrand and Feynman methods in order to find the best fitting domain for calculating the prompt neutron decay constant. The experiments and MCNP simulations performed to accomplish this goal have analyzed the KUKA subcritical assembly facility of Japan driven by 100 MeV protons.

Feynman↗

Kinetic parameters of the GUINEVERE reference configuration in VENUS-F reactor obtained from a pile noise experiment using Rossi and Feynman methods

A pile noise measurement campaign has been conducted by the CEA in the VENUS-F reactor (SCK-CEN, Mol Belgium) in April 2011 in the reference critical configuration of the GUINEVERE experimental program. The experimental setup made it possible to estimate the core kinetic parameters: the prompt neutron decay constant, the delayed neutron fraction and the generation time. A precise assessment of these constants is of prime importance. In particular, the effective delayed neutron fraction is used to normalize and compare calculated reactivities of different subcritical configurations, obtained by modifying either the core layout or the control rods position, with experimental ones deduced from the analysis of measurements. This paper presents results obtained with a CEA-developed time stamping acquisition system. Data were analyzed using Rossi-α and Feynman-α methods. Results were normalized to reactor power using a calibrated fission chamber with a deposit of Np-237. Calculated factors were necessary to the analysis: the Diven factor was computed by the ENEA (Italy) and the power calibration factor by the CNRS/IN2P3/LPC Caen. Results deduced with both methods are consistent with respect to calculated quantities. Recommended values are given by the Rossi-α estimator, that was found to be the most robust. The neutron generation time was found equal to 0.438 ± 0.009 μs and the effective delayed neutron fraction is 765 ± 8 pcm. Discrepancies with the calculated value (722 pcm, calculation from ENEA) are satisfactory: -5.6% for the Rossi-α estimate and -2.7% for the Feynman-α estimate. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Improved methods for Feynman path integral calculations and their application to calculate converged vibrational–rotational partition functions, free energies, enthalpies, entropies, and heat capacities for methane

We present an improved version of our “path-by-path” enhanced same path extrapolation scheme for Feynman path integral (FPI) calculations that permits rapid convergence with discretization errors ranging from O(P{sup −6}) to O(P{sup −12}), where P is the number of path discretization points. We also present two extensions of our importance sampling and stratified sampling schemes for calculating vibrational–rotational partition functions by the FPI method. The first is the use of importance functions for dihedral angles between sets of generalized Jacobi coordinate vectors. The second is an extension of our stratification scheme to allow some strata to be defined based only on coordinate information while other strata are defined based on both the geometry and the energy of the centroid of the Feynman path. These enhanced methods are applied to calculate converged partition functions by FPI methods, and these results are compared to ones obtained earlier by vibrational configuration interaction (VCI) calculations, both calculations being for the Jordan–Gilbert potential energy surface. The earlier VCI calculations are found to agree well (within ∼1.5%) with the new benchmarks. The FPI partition functions presented here are estimated to be converged to within a 2σ statistical uncertainty of between 0.04% and 0.07% for the given potential energy surface for temperatures in the range 300–3000 K and are the most accurately converged partition functions for a given potential energy surface for any molecule with five or more atoms. We also tabulate free energies, enthalpies, entropies, and heat capacities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the Feynman-alpha method for reflected fissile assemblies

The Feynman-alpha method is a neutron noise technique that is used to estimate the prompt neutron period of fissile assemblies. The method and quantity are of widespread interest including in applications such as nuclear criticality safety, safeguards and nonproliferation, and stockpile stewardship; the prompt neutron period may also be used to infer the k eff multiplication factor. The Feynman-alpha method is predicated on time-correlated neutron detections that deviate from a Poisson random variable due to multiplication. Traditionally, such measurements are diagnosed with one-region point kinetics, but two-region models are required when the fissile assembly is reflected. This paper presents a derivation of the two-region point kinetics Feynman-alpha equations based on a double integration of the Rossi-alpha equations, develops novel propagation of measurement uncertainty, and validates the theory. Validation is achieved with organic scintillator measurements of weapons-grade plutonium reflected by various amounts of copper to achieve k eff values of 0.83–0.94 and prompt periods of 5–75 ns. The results demonstrate that Feynman-alpha measurements should use the two-region model instead of the one-region model. The simplified one-region model deviates from the validated two-region models by as much as 10% in the estimate of the prompt neutron period, and the two-region model reduces to the one-region model for small amounts of reflector. The Feynman-alpha estimates of the prompt neutron period are compared to those of the Rossi-alpha approach. The comparative results demonstrate that the Feynman-alpha method is more precise than the Rossi-alpha method and more accurate for k eff < 0.92, whereas the Rossi-alpha method is generally more accurate for higher multiplications. Here, the uncertainty propagation developed in this work should be used for all Feynman-alpha measurements and will therein improve fitting accuracy and appropriate precision estimates.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

MCNP6 Modification for the Calculation of Rossi and Feynman Distributions

The Rossi distribution is given by the number of captured neutron pairs as a function of the time interval between the two capture events. This distribution contains a random part, which is given by neutron captures caused by neutrons of different fission chains, and a correlated part, which is given by neutron captures caused by neutrons from the same fission chain. The correlated part exponentially vanishes with the {Delta}t time interval following the decay of prompt neutrons. The Rossi distribution peaks when the time interval between two capture events is small since this condition is satisfied for neutrons generated from the same fission chain. In the Poisson probability distribution, the variance is equal to the mean and the next event probability is independent of the elapsed time since the previous event. Neutron captures events do not follow the Poisson probability distribution because they are correlated events via fission chains. Consequently, the Feynman distribution Y is defined as the variance to mean ratio minus one, of the neutron counts array as a function of a fixed time interval. The MCNP6 source code has been modified to allow the calculation of the Rossi and Feynman distributions. The source code modification takes into account the time when a neutron is emitted by the external neutron source and gives a precision in the timestamp lower than a fraction of nanosecond. The modified MCNP6 version has been used to model KUKA facility driven by a californium neutron source. The prompt neutron decay constants obtained by the Rossi and Feynman methods are in agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of neutron and gamma-ray noise in large subcritical condition for criticality safety - 14485

Various experimental methods have been proposed for subcriticality measurement and neutron noise analysis method is one of important experimental techniques for this purpose. The variance to mean value method, called Feynman-alpha method, is one of well-known noise analysis methods. In this method, variation of neutron count rate distribution from the Poisson distribution caused by fission chain reactions, which is called Y-value, is measured, however, it is difficult to measure the Y-value in large subcritical state such as nuclear fuel facility, because number of emitted neutrons by fission chain reaction is small, which causes large statistical error. We developed new experimental technique based on Feynman-alpha method combined with pulsed neutron source, called Delayed Neutron Noise Analysis (DNNA) method, to monitor the change of subcriticality. In this new method, neutron and gamma-ray counts measured only appeared in the delayed neutron region after injection of pulsed neutron are analyzed, and the prompt neutron region is completely ignored. The advantage of DNNA compared with conventional method is that spatial dependency of detector positions on measured results can be greatly decreased, which result in increase of accuracy of measured subcriticality. The validity of the DNNA method was experimentally investigated in the Kyoto University Critical Assembly (KUCA) at Kyoto University Research Reactor Institute. In the experiments, pulsed neutrons by D-T reactions are injected into the core and we used {sup 3}He detectors to measure neutron noise signals and BGO scintillation detectors to measure gamma-ray noise signals, and all data in the delayed neutron region were analyzed based on Feynman-alpha method. From the experimental results, it was found that the Y-values are clearly larger than zero in the case of existence of fuel assemblies and this method can be applied to detection of hidden special nuclear materials for example in a cargo container. (authors)

42 ENGINEERING↗

Advances in the computation of the Sjöstrand, Rossi, and Feynman distributions

This study illustrates recent computational advances in the application of the Sjöstrand (area), Rossi, and Feynman methods to estimate the effective multiplication factor of a subcritical system driven by an external neutron source. The methodologies introduced in this study have been validated with the experimental results from the KUKA facility of Japan by Monte Carlo (MCNP6 and MCNPX) and deterministic (ERANOS, VARIANT, and PARTISN) codes. When the assembly is driven by a pulsed neutron source generated by a particle accelerator and delayed neutrons are at equilibrium, the Sjöstrand method becomes extremely fast if the integral of the reaction rate from a single pulse is split into two parts. These two integrals distinguish between the neutron counts during and after the pulse period. To conclude, when the facility is driven by a spontaneous fission neutron source, the timestamps of the detector neutron counts can be obtained up to the nanosecond precision using MCNP6, which allows obtaining the Rossi and Feynman distributions.

97 MATHEMATICS AND COMPUTING↗

Exchange and spin states in quantum dots under strong spatial correlations. Computer simulation by the Feynman path integral method

The fundamental laws in the behavior of electrons in model quantum dots that are caused by exchange and strong Coulomb correlations are studied. The ab initio path integral method is used to numerically simulate systems of two, three, four, and six interacting identical electrons confined in a three-dimensional spherical potential well with a parabolic confining potential against the background of thermal fluctuations. The temperature dependences of spin and collective spin magnetic susceptibility are calculated for model quantum dots of various spatial sizes. A basically exact procedure is proposed for taking into account the permutation symmetry and the spin state of electrons, which makes it possible to perform numerical calculations using modern computer facilities. The conditions of applicability of a virial energy estimator and its optimum form in exchange systems are determined. A correlation estimator of kinetic energy, which is an alternative to a basic estimator, is suggested. A fundamental relation between the kinetic energy of a quantum particle and the character of its virtual diffusion in imaginary time is demonstrated. The process of natural 'pairing' of electron spins during the compression of a quantum dot and cooling of a system is numerically reproduced in terms of path integrals. The temperature dependences of the spin magnetic susceptibility of electron pairs with a characteristic maximum caused by spin pairing are obtained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. In this work, we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Easy-plane anisotropic-exchange magnets on a honeycomb lattice: Quantum effects and dealing with them

We provide analytical and numerical insights into the phase diagram and other properties of the extended Kitaev-Heisenberg model on the honeycomb lattice in the easy-plane limit, in which interactions are only between spin components that belong to the plane of magnetic ions. This parameter subspace allows for a much-needed systematic quantitative investigation of spin excitations in the ordered phases and of their generic features. Specifically, we demonstrate that in this limit one can consistently take into account magnon interactions in both zero-field zigzag and field-polarized phases. For the nominally polarized phase, we propose a regularization of the unphysical divergences that occur at the critical field and are plaguing the 1/S approximation in this class of models. For the explored parameter subspace, all symmetry-allowed terms of the standard parametrization of the extended Kitaev-Heisenberg model, such as K, J, and Γ, are significant, making the offered consideration relevant to a much wider parameter space. Furthermore, the dynamical structure factor near the paramagnetic critical point illustrates this relevance by showing features that are reminiscent of the ones observed in α–RuCl 3 , underscoring that they are not unique and should be common to a wide range of parameters of the model and, by extension, to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reactivity measurements using the Zolotukhin-Mogilner Method

The zero count probability method (also called Zolotukhin-Mogilner Method, after its inventors) is a method for measuring the reactivity in nuclear reactors, being a competitor of the [more well-known in the West] Feynman-alpha Method and the Rossi-alpha Method. The modus operandi of this method is using the detector channels where there was no counting at all in order to calculate reactor parameters. In deep subcriticality few models have been tested and this work tries out the Zolotukhin-Mogilner Method in one of such scenarios: measurements will be made in environments below -3,5 k pcm in the zero-potency water-moderated reactor IPEN/MB-01 which is fuelled by UO{sub 2} enriched by 4.3%. These extremely low reactivity environments are required because the chance of no counts on the detector must be significant: otherwise, the method would demand a large time of acquisition of data. Besides that, the method is very simple and straightforward. One of the advantages of this method is that it needs very little data reduction, since the essential data is directly given by the measuring apparatus. The detection will be in charge of modern BF{sub 3} detectors. It will be assumed that, in these deep subcritical systems, the function K{sub eff} = f(N-Nkp) has a linear portion in its first part, resulting of the decomposition of it in a Taylor series. The value of alpha is related to the reactivity ρ with linear dependency. The results will be compared with recent studies of the two main methods described above. Presented for the first time in the 60's, this tool has seen little use in the west hemisphere. This work shows its use in the measurements of the nuclear reactor IPEN/MB-01, as well as the code developed for its employment. It will be the first time this method is used in the south hemisphere. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SU-E-E-08: Applications of the Quantization of Coupled Circuits in Radiation Physics (design of Klystron, Bremsstrahlung, Synchrotron)

Purpose: During the past decade the quantization of coupled/forced electromagnetic circuits with or without Ohm’s resistance has gained the subject of some fundamental studies, since even problems of quantum electrodynamics can be solved in an elegant manner, e.g. the creation of quantized electromagnetic fields. In this communication, we shall use these principles to describe optimization procedures in the design of klystrons, synchrotron irradiation and high energy bremsstrahlung. Methods: The base is the Hamiltonian of an electromagnetic circuit and the extension to coupled circuits, which allow the study of symmetries and perturbed symmetries in a very apparent way (SU2, SU3, SU4). The introduction resistance and forced oscillators for the emission and absorption in such coupled systems provides characteristic resonance conditions, and atomic orbitals can be described by that. The extension to virtual orbitals leads to creation of bremsstrahlung, if the incident electron (velocity v nearly c) is described by a current, which is associated with its inductivitance and the virtual orbital to the charge distribution (capacitance). Coupled systems with forced oscillators can be used to amplify drastically the resonance frequencies to describe klystrons and synchrotron radiation. Results: The cross-section formula for bremsstrahlung given by the propagator method of Feynman can readily be derived. The design of klystrons and synchrotrons inclusive the radiation outcome can be described and optimized by the determination of the mutual magnetic couplings between the oscillators induced by the currents. Conclusions: The presented methods of quantization of circuits inclusive resistance provide rather a straightforward way to understand complex technical processes such as creation of bremsstrahlung or creation of radiation by klystrons and synchrotrons. They can either be used for optimization procedures and, last but not least, for pedagogical purposes with regard to a qualified understanding of radiation physics for students.

43 PARTICLE ACCELERATORS↗

Optimization of Binning Parameters in the Feynman Variance-to-Mean Method for Delayed Neutron Reinterrogation [Poster]

The purpose of this work is to determine the optimal input parameters for analyzing active interrogation data that utilizes the physics of the delayed neutron re-interrogation technique. This work compares the signature from an active neutron measurement of a highly enriched uranium (HEU) system to an active background measurement (i.e., without the HEU item in place). In addition, varying analyses of the NoMAD list-mode detector data are performed to determine whether it is preferential to use a particular analysis bin size for producing Feynman histograms, selection of pulses to prevent discrepancies from buildup, or window of time during the pulsed neutron interrogation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimization of Binning Parameters in the Feynman Variance-to-Mean Method for Delayed Neutron Reinterrogation

The focus of this work is to determine the optimal input parameters for analyzing active interrogation data that utilizes the physics of the delayed neutron re-interrogation technique. This work compares the signature from an active neutron measurement of a highly enriched uranium (HEU) system to an active background measurement (i.e., without the HEU item in place). The item being measured represents a 21.8 kg metal sphere of HEU (93 weight percent) assembled from a subset of the Rocky Flats (RF) uranium hemishells. The detection system consisted of two NoMAD neutron multiplicity detectors, a SNAP III neutron detector, and a Thermo Fisher P211 pulsed 14 MeV neutron source. In addition, varying analyses of the NoMAD list-mode detector data are performed to determine whether it is preferential to use a particular analysis bin size for producing Feynman histograms, selection of pulses to prevent discrepancies from buildup, or window of time during the pulsed neutron interrogation. This work is still currently ongoing.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Feynman’s proof of the commutativity of the Calogero integrals of motion

During his last year of life Feynman became interested in integrable models. In the course of his study of one dimensional particles with inverse square potentials (the Calogero model) he came up with a field theory-inspired proof of the commutativity of the integrals of motion of the system, which remained unpublished. I present an exact transcript of the proof together with a scan of the first page of Feynman’s original manuscript and some historical background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Information-based physics, influence, and forces

In recent works, Knuth and Bahreyni have demonstrated that the concepts of space and time are emergent in a coarse-grained model of direct particle-particle influence. In addition, Knuth demonstrated that observer-made inferences regarding the free particle, which is defined as a particle that influences others, but is not itself influenced, result in a situation identical to the Feynman checkerboard model of the Dirac equation. This suggests that the same theoretical framework that gives rise to an emergent spacetime is consistent with quantum mechanics. In this paper, we begin to explore the effect of influence on the emergent properties of a particle. This initial study suggests that when a particle is influenced, it is interpreted as accelerating in a manner consistent with special relativity implying that, at least in this situation, influence can be conceived of as a force.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Feynman-α analysis using BGO gamma-ray detector in UTR-KINKI reactor

In conventional reactor noise measurement based on neutron detection, a neutron detector must be placed as close to a core as possible, to observe a largest possible correlation amplitude. For most of reactors other than critical assemblies, however, the placement is hardly possible. We focused on an alternative reactor noise measurement using prompt gamma-rays having a long flight range, to overcome the above constraint. A Feynman-α analysis using a NaI(Tl) scintillation detector was previously performed in UTR-KINKI, however, white noise originating from radioactivation of the scintillator significantly degraded the correlation amplitude. In this study, another Feynman-α analysis using a BGO gamma-ray detector, which never suffers from any radioactivation, is performed. A series of Feynman-α analyses for a critical state and several subcritical states of UTR-KINKI reactor has been carried out using two BGO detectors and a BF{sub 3} neutron detector to determine prompt-neutron decay constant. In the present analysis, the dead-time effect of these detectors and the delayed-neutron effect are considered. The prompt-neutron decay constant determined from Feynman-α analysis based on gamma-ray detection agrees with that from the analysis based on neutron detection. The Feynman-α analysis based on detecting gamma rays appears to be able to determine prompt-neutron decay constant at least as efficiently as the conventional analysis detecting neutrons.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗