Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical Physics”

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 829 records · Page 46

MCNP6 Results for the Phase III Sensitivity Benchmark of the OCED/NEA Expert Group on Uncertainty Analysis for Criticality Safety Assessment

Within the last decade, there has been increasing interest in the calculation of cross section sensitivity coefficients of k eff for integral experiment design and uncertainty analysis. The OECD/NEA has an Expert Group devoted to Sensitivity and Uncertainty Analysis within the Working Party for Nuclear Criticality Safety. This expert group has developed benchmarks to assess code capabilities and performance for doing sensitivity and uncertainty analysis. Phase III of a set of sensitivity benchmarks evaluates capabilities for computing sensitivity coefficients. MCNP6 has the capability to compute cross section sensitivities for k eff using continuous-energy physics. To help verify this capability, results for the Phase III benchmark cases are generated and submitted to the Expert Group for comparison. The Phase III benchmark has three cases: III.1, an array of MOX fuel pins, III.2, a series of infinite lattices of MOX fuel pins with varying pitches, and III.3 two spheres with homogeneous mixtures of UF 4 and polyethylene with different enrichments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Systematic and statistical uncertainties in simulated r -process abundances due to uncertain nuclear masses

Unknown nuclear masses are a major source of nuclear physics uncertainty for r -process nucleosynthesis calculations. Here we examine the systematic and statistical uncertainties that arise in r -process abundance predictions due to uncertainties in the masses of nuclear species on the neutron-rich side of stability. There is a long history of examining systematic uncertainties by the application of a variety of different mass models to r-process calculations. Here we expand upon such efforts by examining six DFT mass models, where we capture the full impact of each mass model by updating the other nuclear properties — including neutron capture rates, β-decay lifetimes, and β-delayed neutron emission probabilities — that depend on the masses. Unlike systematic effects, statistical uncertainties in the r-process pattern have just begun to be explored. Here we apply a global Monte Carlo approach, starting from the latest FRDM masses and considering random mass variations within the FRDM r ms error. Here, we find in each approach that uncertain nuclear masses produce dramatic uncertainties in calculated r -process yields, which can be reduced in upcoming experimental campaigns.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The 2018 Edition of the IRPhEP Handbook

The International Reactor Physics Experiment Evaluation Project (IRPhEP) was established to preserve and evaluate integral reactor physics experiment data, including separate or special effects measurements for nuclear energy and technology applications. Numerous international experiments have been performed over the past several decades to support reactor operations, safety, measurements, and design. Those experiments represent a significant investment of infrastructure, cost, and expertise that, in many instances, cannot be as easily reproduced. They represent valuable resources of data for present and future research, and provide the basis for recording, development, and validation of methods. Contributors worldwide collaborate to evaluate and review selected benchmark experiments for inclusion in the International Handbook of Evaluated Reactor Physics Benchmark Experiments (IRPhEP Handbook). Several new evaluations, and revisions to existing benchmarks have been prepared for inclusion in the 2018 edition of the IRPhEP Handbook. Additionally, an IRPhEP Uncertainty Guide was discussed to provide further guidance regarding the evaluation of reactor physics measurements besides criticality. Recent changes to the IRPhEP Handbook include five minor revisions, the revision of a previous benchmark evaluation to include new data, publication of a new draft benchmark evaluation, final approval for the formal publication of two formerly draft benchmark evaluations, and the inclusion of ten new benchmark evaluations. The IRPhEP Handbook now contains a total of 159 experimental measurement series from 54 unique reactor facilities. International participants are actively engaged in ongoing evaluation activities to prepare content for future editions of the IRPhEP Handbook. There are numerous experiment and reactor facilities worldwide from which benchmark data can be evaluated. Those interested in serving as evaluators and/or reviewers for the IRPhEP are encouraged to participate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An evidential reasoning extension to quantitative model-based failure diagnosis

The detection and diagnosis of failures in physical systems characterized by continuous-time operation are studied. A quantitative diagnostic methodology has been developed that utilizes the mathematical model of the physical system. On the basis of the latter, diagnostic models are derived each of which comprises a set of orthogonal parity equations. To improve the robustness of the algorithm, several models may be used in parallel, providing potentially incomplete and/or conflicting inferences. Dempster's rule of combination is used to integrate evidence from the different models. The basic probability measures are assigned utilizing quantitative information extracted from the mathematical model and from online computation performed therewith.

Gertler, Janos J.↗

Internal Material Melting Development Progress [Slides]

Slides for an AFC workshop giving a brief overview of the current status of Bison code development for correct energy balance of melting material. The slides provide background information of why the development is needed. The planned approach of casting the energy balance as a change in enthalpy is described. Status of current code additions is outlined. Extra slides are provided with more mathematical details of the approach derivation and verification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Some Problems in Density Functional Theory

Though calculations based on density functional theory (DFT) are used remarkably widely in chemistry, physics, materials science, and biomolecular research and though the modern form of DFT has been studied for almost 60 years, some mathematical problems remain. From a physical science perspective, it is far from clear whether those problems are of major import. For context, we provide an outline of the basic structure of DFT as it is presented and used conventionally in physical sciences, note some unresolved mathematical difficulties with those conventional demonstrations, then pose several questions regarding both the time-independent and time-dependent forms of DFT that could benefit from attention in applied mathematics. Furthermore, progress on any of these would aid in development of better approximate functionals and in interpretation of DFT.

36 MATERIALS SCIENCE↗

Beyond Archimedean Space-Time Structure

It took two millennia after Euclid and until in the early 1880s, when we went beyond the ancient axiom of parallels, and inaugurated geometries of curved spaces. In less than one more century, General Relativity followed. At present, physical thinking is still beheld by the yet deeper and equally ancient Archimedean assumption. In view of that, it is argued with some rather easily accessible mathematical support that Theoretical Physics may at last venture into the non-Archimedean realms. In this introductory paper we stress two fundamental consequences of the non-Archimedean approach to Theoretical Physics: one of them for quantum theory and another for relativity theory. From the non-Archimedean viewpoint, the assumption of the existence of minimal quanta of light (of the fixed frequency) is an artifact of the present Archimedean mathematical basis of quantum mechanics. In the same way the assumption of the existence of the maximal velocity, the velocity of light, is a feature of the real space-time structure which is fundamentally Archimedean. Both these assumptions are not justified in corresponding non-Archimedean models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Group action in topos quantum physics

Topos theory has been suggested first by Isham and Butterfield, and then by Isham and Doering, as an alternative mathematical structure within which to formulate physical theories. In particular, it has been used to reformulate standard quantum mechanics in such a way that a novel type of logic is used to represent propositions. In this paper, we extend this formulation to include the notion of a group and group transformation in such a way that we overcome the problem of twisted presheaves. In order to implement this we need to change the type of topos involved, so as to render the notion of continuity of the group action meaningful.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NLOX, a one-loop provider for Standard Model processes

NLOX is a computer program for calculations in high-energy particle physics. It provides fully renormalized scattering matrix elements in the Standard Model of particle physics, up to one-loop accuracy for all possible coupling-power combinations in the strong and electroweak couplings, and for processes with up to six external particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmark of the KGMf with a coupled Boltzmann equation solver

The Kinetic Global Model framework (KGMf) is an open-source general-purpose global model (spatially averaged) simulation code developed to explore the reaction kinetics and pathways in plasma discharge systems. It contains species continuity and electron energy balance equations, with a time-dependent evaluated electron energy distribution function (EEDF) for electron impact reactions. The EEDF is utilized to determine the rate coefficients for electron impact reactions, which can have profound impact on the temporal evolutions of plasma parameters. Previously, the EEDF was commonly assumed as Maxwellian or an analytical function of the effective electron temperature. In this work, the KGMf is coupled with a Boltzmann equation (BE) solver to self-consistently compute the EEDF. The EEDF evolution frequency is determined based on relative changes of the reduced electric field. The KGMf is benchmarked with the ZDPlasKin code based on high-pressure low-temperature argon plasma discharge cases. Additionally, the temporal evolutions of reduced electric field, electron temperature, EEDF, reaction rates, and species densities, are obtained and compared under different discharge conditions, showing good agreement between the KGMf and the ZDPlasKin simulations. The application of the KGMf for predicting breakdown times in high power microwave discharges is also presented, which shows qualitative agreement with particle-in-cell simulations. The KGMf can be further applied for more complicated plasma discharge systems, where the reaction kinetics are more intricate (e.g., plasma-assisted combustion systems).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT

A consistent “2D/1D” neutron transport method is derived from the 3D Boltzmann transport equation, to calculate fuel-pin-resolved neutron fluxes for realistic full-core Pressurized Water Reactor (PWR) problems. The 2D/1D method employs the Method of Characteristics to discretize the radial variables and a lower order transport solution to discretize the axial variable. This paper describes the theory of the 2D/1D method and its implementation in the MPACT code, which has become the whole-core deterministic neutron transport solver for the Consortium for Advanced Simulations of Light Water Reactors (CASL) core simulator VERA-CS. Several applications have been performed on both leadership-class and industry-class computing clusters. Results are presented for whole-core solutions of the Watts Bar Nuclear Power Station Unit 1 and compared to both continuous-energy Monte Carlo results and plant data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An implicit, sharp numerical treatment of viscous terms at arbitrarily shaped liquid-gas interfaces in evaporative flows

We report a novel methodology for accurate treatment of viscous terms in evaporation problems. The proposed scheme is an extension of the sharp viscous treatment of Kang et al. (2000) [7] to 3D phase change problems. To ensure accuracy and grid converging solutions, a new implicit approach to computing viscous fluxes across the phase interfaces is proposed, which was previously unavailable in fixed grid numerical schemes. Analytical relations were derived for the jump in velocity gradients across the 2D and 3D phase interfaces, an important constituent of the proposed scheme. The relations show a non-vanishing jump in the tangential gradients across the phase interface that are associated with evaporative flux and interfacial curvature. The proposed methodology demonstrated first order accuracy in canonical test cases. It is general and applicable to arbitrarily oriented interfaces, and can be readily implemented in existing evaporation flow solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

When and why PINNs fail to train: A neural tangent kernel perspective

Physics-informed neural networks (PINNs) have lately received great attention thanks to their flexibility in tackling a wide range of forward and inverse problems involving partial differential equations. However, despite their noticeable empirical success, little is known about how such constrained neural networks behave during their training via gradient descent. More importantly, even less is known about why such models sometimes fail to train at all. Here in this work, we aim to investigate these questions through the lens of the Neural Tangent Kernel (NTK); a kernel that captures the behavior of fully-connected neural networks in the infinite width limit during training via gradient descent. Specifically, we derive the NTK of PINNs and prove that, under appropriate conditions, it converges to a deterministic kernel that stays constant during training in the infinite-width limit. This allows us to analyze the training dynamics of PINNs through the lens of their limiting NTK and find a remarkable discrepancy in the convergence rate of the different loss components contributing to the total training error. To address this fundamental pathology, we propose a novel gradient descent algorithm that utilizes the eigenvalues of the NTK to adaptively calibrate the convergence rate of the total training error. Finally, we perform a series of numerical experiments to verify the correctness of our theory and the practical effectiveness of the proposed algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SympNets: Intrinsic structure-preserving symplectic networks for identifying Hamiltonian systems

Here, we propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Stabilizer Framework for the Contextual Subspace Variational Quantum Eigensolver and the Noncontextual Projection Ansatz

Quantum chemistry is a promising application for noisy intermediate-scale quantum (NISQ) devices. However, quantum computers have thus far not succeeded in providing solutions to problems of real scientific significance, with algorithmic advances being necessary to fully utilize even the modest NISQ machines available today. We discuss a method of ground state energy estimation predicated on a partitioning of the molecular Hamiltonian into two parts: one that is noncontextual and can be solved classically, supplemented by a contextual component that yields quantum corrections obtained via a Variational Quantum Eigensolver (VQE) routine. This approach has been termed Contextual Subspace VQE (CS-VQE); however, there are obstacles to overcome before it can be deployed on NISQ devices. The problem we address here is that of the ansatz, a parametrized quantum state over which we optimize during VQE; it is not initially clear how a splitting of the Hamiltonian should be reflected in the CS-VQE ansätze. We propose a “noncontextual projection” approach that is illuminated by a reformulation of CS-VQE in the stabilizer formalism. This defines an ansatz restriction from the full electronic structure problem to the contextual subspace and facilitates an implementation of CS-VQE that may be deployed on NISQ devices. We validate the noncontextual projection ansatz using a quantum simulator and demonstrate chemically precise ground state energy calculations for a suite of small molecules at a significant reduction in the required qubit count and circuit depth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Applying Bayesian inference and deterministic anisotropy to retrieve the molecular structure ∣Ψ(R)∣2 distribution from gas-phase diffraction experiments

Abstract Currently, our general approach to retrieving molecular structures from ultrafast gas-phase diffraction heavily relies on complex ab initio electronic or vibrational excited state simulations to make conclusive interpretations. Without such simulations, inverting this measurement for the structural probability distribution is typically intractable. This creates a so-called inverse problem. Here we address this inverse problem by developing a broadly applicable method that approximates the molecular frame structure ∣Ψ( R , t )∣ 2 distribution independent of these complex simulations. We retrieve the vibronic ground state ∣Ψ( R )∣ 2 for both simulated stretched NO 2 and measured N 2 O. From measured N 2 O, we observe 40 mÅ coordinate-space resolution from 3.75 Å −1 reciprocal space range and poor signal-to-noise, a 50X improvement over traditional Fourier transform methods. In simulated NO 2 diffraction experiments, typical to high signal-to-noise levels predict 100–1000X resolution improvements, down to 0.1 mÅ. By directly measuring the width of ∣Ψ( R )∣ 2 , we open ultrafast gas-phase diffraction capabilities to measurements beyond current analysis approaches. This method has the potential to effectively turn gas-phase ultrafast diffraction into a discovery-oriented technique to probe systems that are prohibitively difficult to simulate.

74 ATOMIC AND MOLECULAR PHYSICS↗

Transition from gas to plasma kinetic equilibria in gravitating axisymmetric structures

The problem of the transition from gas to plasma in gravitating axisymmetric structures is addressed under the assumption of having initial and final states realized by kinetic Maxwellian-like equilibria. In astrophysics, the theory applies to accretion-disc scenarios around compact objects. A formulation based on non-relativistic kinetic theory for collisionless systems is adopted. Equilibrium solutions for the kinetic distribution functions describing the initial neutral matter and the resulting plasma state are constructed in terms of single-particle invariants and expressed by generalized Maxwellian distributions. The final plasma configuration is related to the initial gas distribution by the introduction of appropriate functional constraints. Qualitative aspects of the solution are investigated and physical properties of the system are pointed out. In particular, the admitted functional dependences of the fluid fields carried by the corresponding equilibrium distributions are determined. Then, the plasma is proved to violate the condition of quasi-neutrality, implying a net charge separation between ions and electrons. This result is shown to be independent of the precise realization of the plasma distribution function, while a physical mechanism able to support a non-neutral equilibrium state is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Skyrmion-antiskyrmion pairs in ultracold atomic gases

We study theoretically the dynamics of two-component Bose-Einstein condensates in two dimensions, which admit topological excitations related to the skyrmions of nuclear physics. We investigate a branch of uniformly propagating solitary waves, which, at high momentum, can be viewed as skyrmion-antiskyrmion pairs. We study these solitary waves for a range of interaction regimes and show that, for experimentally relevant cases, there is a transition to spatially extended spin-wave states at low momentum. We explain how this can be understood by analogy to the two-dimensional ferromagnet. Finally, we discuss how such solitary waves can be generated and studied in experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗