Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational 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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Qutrit and qubit circuits for three-flavor collective neutrino oscillations

We explore the utility of qutrits and qubits for simulating the flavor dynamics of dense neutrino systems. The evolution of such systems impacts some important astrophysical processes, such as core-collapse supernovae and the nucleosynthesis of heavy nuclei. Many-body simulations require classical resources beyond current computing capabilities for physically relevant system sizes. Quantum computers are therefore a promising candidate to efficiently simulate the many-body dynamics of collective neutrino oscillations. Previous quantum simulation efforts have primarily focused on properties of the two-flavor approximation due to their direct mapping to qubits. Furthermore, we present new quantum circuits for simulating three-flavor neutrino systems on qutrit- and qubit-based platforms, and demonstrate their feasibility by simulating systems of two, four, and eight neutrinos on IBM and Quantinuum quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

NSRD-11: Computational Capability to Substantiate DOE-HDBK-3010 Data.

Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Handbook, DOE - HDBK - 3010, Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities, to determine radionuclide source terms. In calculating source terms, analysts tend to use the DOE Handbook's bounding values on airborne release fractions (ARFs) and respirable fractions (RFs) for various categories of insults (representing potential accident release categories). This is typically due to both time constraints and the avoidance of regulatory critique. Unfortunately, these bounding ARFs/RFs represent extremely conservative values. Moreover, they were derived from very limited small-scale bench/laboratory experiments and/or from engineered judgment. Thus, the basis for the data may not be representative of the actual unique accident conditions and configurations being evaluated. The goal of this research is to develop a more accurate and defensible method to determine bounding values for the DOE Handbook using state-of-art multi-physics-based computer codes. This enables us to better understand the fundamental physics and phenomena associated with the types of accidents in the handbook. In this year, this research included improvements of the high-fidelity codes to model particle resuspension and multi-component evaporation for fire scenarios. We also began to model ceramic fragmentation experiments, and to reanalyze the liquid fire and powder release experiments that were done last year. The results show that the added physics better describes the fragmentation phenomena. Thus, this work provides a low-cost method to establish physics-justified safety bounds by taking into account specific geometries and conditions that may not have been previously measured and/or are too costly to perform.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

VAIM-CFF: a variational autoencoder inverse mapper solution to Compton form factor extraction from deeply virtual exclusive reactions

We develop a new methodology for extracting Compton form factors (CFFs) from deeply virtual exclusive reactions such as the unpolarized DVCS cross section using a specialized inverse problem solver, a variational autoencoder inverse mapper (VAIM). The VAIM-CFF framework not only allows us access to a fitted solution set possibly containing multiple solutions in the extraction of all 8 CFFs from a single cross section measurement, but also accesses the lost information contained in the forward mapping from CFFs to cross section. We investigate various assumptions and their effects on the predicted CFFs such as cross section organization, number of extracted CFFs, use of uncertainty quantification technique, and inclusion of prior physics information. We then use dimensionality reduction techniques such as principal component analysis to visualize the missing physics information tracked in the latent space of the VAIM framework. Through re-framing the extraction of CFFs as an inverse problem, we gain access to fundamental properties of the problem not comprehensible in standard fitting methodologies: exploring the limits of the information encoded in deeply virtual exclusive experiments.

Accelerator Physics↗
DOE OSTI2013

Development of a fourth generation predictive capability maturity model.

The Predictive Capability Maturity Model (PCMM) is an expert elicitation tool designed to characterize and communicate completeness of the approaches used for computational model definition, verification, validation, and uncertainty quantification associated for an intended application. The primary application of this tool at Sandia National Laboratories (SNL) has been for physics-based computational simulations in support of nuclear weapons applications. The two main goals of a PCMM evaluation are 1) the communication of computational simulation capability, accurately and transparently, and 2) the development of input for effective planning. As a result of the increasing importance of computational simulation to SNL s mission, the PCMM has evolved through multiple generations with the goal to provide more clarity, rigor, and completeness in its application. This report describes the approach used to develop the fourth generation of the PCMM.

Hills, Richard Guy↗
NASA NTRS1984

CREKID: A computer code for transient, gas-phase combustion of kinetics

A new algorithm was developed for fast, automatic integration of chemical kinetic rate equations describing homogeneous, gas-phase combustion at constant pressure. Particular attention is paid to the distinguishing physical and computational characteristics of the induction, heat-release and equilibration regimes. The two-part predictor-corrector algorithm, based on an exponentially-fitted trapezoidal rule, includes filtering of ill-posed initial conditions, automatic selection of Newton-Jacobi or Newton iteration for convergence to achieve maximum computational efficiency while observing a prescribed error tolerance. The new algorithm was found to compare favorably with LSODE on two representative test problems drawn from combustion kinetics.

Pratt, D. T.↗
DOE OSTI2022

Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2

Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO 2 , it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO 2 . While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO 2 -brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO 2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO 2 obtained from core flow experiments and direct numerical simulation.

Kohanpur, Amir H.↗
DOE OSTI2020

Nonperturbative matching between equal-time and lightcone quantization

We investigate the nonperturbative relation between lightcone (LC) and standard equal-time (ET) quantization in the context of λΦ4 theory in d = 2. We discuss the perturbative matching between bare parameters and the failure of its naive nonperturbative extension. We argue that they are nevertheless the same theory nonperturbatively, and that furthermore the nonperturbative map between bare parameters can be extracted from ET perturbation theory via Borel resummation of the mass gap. We test this map by using it to compare physical quantities computed using numerical Hamiltonian truncation methods in ET and LC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Image Analysis: Comparison Metrics

Despite an increased reliance on computational modeling in engineering and physics, graphically comparing computational and experimental results has often been qualitative, via a so-called ‘viewgraph norm’. Image recognition, computer vision and artificial intelligence (AI) are fields of study in themselves and rapidly progressing, but relying on AI to grade image similarity evokes a notion of asking for an expert judgement, which could be seen as an artificial version of the viewgraph norm. It is, therefore, desirable to use simpler metrics which are more tractable and unambiguous, even though they may not be as ‘intelligent.’ In this document, a few metrics are studied by systematically altering an image and their implications are analyzed

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Digital system for structural dynamics simulation

State-of-the-art digital hardware and software for the simulation of complex structural dynamic interactions, such as those which occur in rotating structures (engine systems). System were incorporated in a designed to use an array of processors in which the computation for each physical subelement or functional subsystem would be assigned to a single specific processor in the simulator. These node processors are microprogrammed bit-slice microcomputers which function autonomously and can communicate with each other and a central control minicomputer over parallel digital lines. Inter-processor nearest neighbor communications busses pass the constants which represent physical constraints and boundary conditions. The node processors are connected to the six nearest neighbor node processors to simulate the actual physical interface of real substructures. Computer generated finite element mesh and force models can be developed with the aid of the central control minicomputer. The control computer also oversees the animation of a graphics display system, disk-based mass storage along with the individual processing elements.

Krauter, A. I.↗
NASA NTRS1984

Computing in the presence of soft bit errors

It is shown that single-event-upsets (SEUs) due to cosmic rays are a significant source of single bit error in spacecraft computers. The physical mechanism of SEU, electron hole generation by means of Linear Energy Transfer (LET), it discussed with reference made to the results of a study of the environmental effects on computer systems of the Galileo spacecraft. Techniques for making software more tolerant of cosmic ray effects are considered, including: reducing the number of registers used by the software; continuity testing of variables; redundant execution of major procedures for error detection; and encoding state variables to detect single-bit changes. Attention is also given to design modifications which may reduce the cosmic ray exposure of on-board hardware. These modifications include: shielding components operating in LEO; removing low-power Schottky parts; and the use of CMOS diodes. The SEU parameters of different electronic components are listed in a table.

Rasmussen, R. D.↗
NASA NTRS1987

Unified viscoplastic constitutive equations and their applications

Unified constitutive equations for time- and temperature-dependent metallic plastic deformation have been applied in FEM simulations of forming processes; increasingly powerful computational tools and physical models are being used to numerically model complex engineering problems. Once confidence has been gained through adequate verification, these numerical models will increasingly replace experimental models. Attention is presently given to the contributions made by physical metallurgy, continuum mechanics, and computational mechanics.

Lindholm, U. S.↗
DOE OSTI2021

Separating Physics and Dynamics Grids for Improved Computational Efficiency in Spectral Element Earth System Models

Previous studies have shown that atmospheric models with a spectral element grid can benefit from putting physics calculations on a relatively coarse finite volume grid. Here we demonstrate an alternative high-order, element-based mapping approach used to implement a quasi-equal-area, finite volume physics grid in E3SM. Unlike similar methods, the new method in E3SM requires topology data purely local to each spectral element, which trivially allows for regional mesh refinement. Simulations with physics grids defined by 2 × 2, 3 × 3, and 4 × 4 divisions of each element are shown to verify that the alternative physics grid does not qualitatively alter the model solution. The model performance is substantially affected by the reduction of physics columns when using the 2 × 2 grid, which can increase the throughput of physics calculations by roughly 60%–120% depending on whether the computational resources are configured to maximize throughput or efficiency. A pair of regionally refined cases are also shown to highlight the refinement capability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Using Python For 3D Visualizations of MCNP output

MCNP is ubiquitously used within industries, universities and national laboratories as a tool for accurately modeling radiation. With so much focus on the details of transport, attention is seldom given to the visualizations of these meticulously constructed outputs. Python is a programing language designed to promote code readability and has steadily grown in popularity since the 1990's. Python is also a heavily crowd-sourced programing language which generally means that the problems that you are trying to solve with computer codes are already solved; The user is often left with the now-simplified task of cobbling together a good and useful compilation of these tools and adapting them to solve the problem at hand. I have taken advantage of this situation and constructed some visualization tools in python that can be used to duplicate some of the 2D radiation fields that MCNP produces and also produce a meaningful 3D representation of a radiation field. (author)

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1992

Domain decomposition: A bridge between nature and parallel computers

Domain decomposition is an intuitive organizing principle for a partial differential equation (PDE) computation, both physically and architecturally. However, its significance extends beyond the readily apparent issues of geometry and discretization, on one hand, and of modular software and distributed hardware, on the other. Engineering and computer science aspects are bridged by an old but recently enriched mathematical theory that offers the subject not only unity, but also tools for analysis and generalization. Domain decomposition induces function-space and operator decompositions with valuable properties. Function-space bases and operator splittings that are not derived from domain decompositions generally lack one or more of these properties. The evolution of domain decomposition methods for elliptically dominated problems has linked two major algorithmic developments of the last 15 years: multilevel and Krylov methods. Domain decomposition methods may be considered descendants of both classes with an inheritance from each: they are nearly optimal and at the same time efficiently parallelizable. Many computationally driven application areas are ripe for these developments. A progression is made from a mathematically informal motivation for domain decomposition methods to a specific focus on fluid dynamics applications. To be introductory rather than comprehensive, simple examples are provided while convergence proofs and algorithmic details are left to the original references; however, an attempt is made to convey their most salient features, especially where this leads to algorithmic insight.

Keyes, David E.↗
DOE OSTI2017

Validation of High-Fidelity Reactor Physics Models for Support of the KJRR Experimental Campaign in the Advanced Test Reactor

The Korea Atomic Energy Research Institute is currently in the process of qualifying a Low-Enriched Uranium fuel element design for the new Ki-Jang Research Reactor (KJRR). As part of this effort, a prototype KJRR fuel element was irradiated for several operating cycles in the Northeast Flux Trap of the Advanced Test Reactor (ATR) at the Idaho National Laboratory. The KJRR fuel element contained a very large quantity of fissile material (618g 235U) in comparison with historical ATR experiment standards (<1g 235U), and its presence in the ATR flux trap was expected to create a neutronic configuration that would be well outside of the approved validation envelope for the reactor physics analysis methods used to support ATR operations. Accordingly it was necessary, prior to high-power irradiation of the KJRR fuel element in the ATR, to conduct an extensive set of new low-power physics measurements with the KJRR fuel element installed in the ATR Critical Facility (ATRC), a companion facility to the ATR that is located in an immediately adjacent building, sharing the same fuel handling and storage canal. The new measurements had the objective of expanding the validation envelope for the computational reactor physics tools used to support ATR operations and safety analysis to include the planned KJRR irradiation in the ATR and similar experiments that are anticipated in the future. The computational and experimental results demonstrated that the neutronic behavior of the KJRR fuel element in the ATRC is well-understood, both in terms of its general effects on core excess reactivity and fission power distributions, its effects on the calibration of the core lobe power measurement system, as well as in terms of its own internal fission rate distribution and total fission power per unit ATRC core power. Taken as a whole, these results have significantly extended the ATR physics validation envelope, thereby enabling an entire new class of irradiation experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Direct observation of dislocation structure evolution in SRF cavity niobium using electron channeling contrast imaging

Two major goals of superconducting radio frequency cavity research and development are achieving higher accelerating gradient and gaining the highest quality factor, a measure of efficiency. However, the consistent improvement of these performance metrics is restricted by many factors, one of which is microstructural defects, such as dislocation substructures within the material. In this work, dislocation evolution is compared in four samples extracted from a 2.8 mm thick large-grain niobium slice, with tensile axes chosen to generate specific dislocation structures during subsequent deformation. The four samples are (1) as-extracted, (2) extracted and annealed, (3) extracted and then deformed to 40% tensile strain, and (4) extracted, annealed at 800 °C 2 h, and then deformed to 40% strain. Electron channeling contrast imaging was performed on all samples to characterize initial dislocation density and dislocation structure evolution due to annealing and deformation, and related to the mechanical behavior observed in stress-strain curves. A fundamental understanding of dislocation evolution in niobium is necessary to develop computational models to simulate cavity forming, which could enable new processing methods for cavity fabrication to be identified that could lead to improved performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mean Field Theory of Self‐Organizing Memristive Connectomes

Abstract Biological neuronal networks are characterized by nonlinear interactions and complex connectivity. Given the growing impetus to build neuromorphic computers, understanding physical devices that exhibit structures and functionalities similar to biological neural networks is an important step toward this goal. Self‐organizing circuits of nanodevices are at the forefront of the research in neuromorphic computing, as their behavior mimics synaptic plasticity features of biological neuronal circuits. However, an effective theory to describe their behavior is lacking. This study provides for the first time an effective mean field theory for the emergent voltage‐induced polymorphism of circuits of a nanowire connectome, showing that the behavior of these circuits can be explained by a low‐dimensional dynamical equation. The equation can be derived from the microscopic dynamics of a single memristive junction in analytical form. The effective model is tested on experiments of nanowire networks and show that it fits both the potentiation and depression of these synapse‐mimicking circuits. It is shown that this theory applies beyond the case of nanowire networks by formulating a general mean‐field theory of conductance transitions in self‐organizing memristive connectomes.

Caravelli, Francesco↗