Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer experiments”

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 1,135 records · Page 63

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Flow physics

A discussion is presented on the coupling of computational analysis and experiment. It is believed that this coupling is critical in developing new aerodynamic insights. Additionally, new methods for analyzing and interpreting data are discussed. These methods need to be developed in small-scale research studies and then applied to large-scale technology programs. The specific objectives of this program are threefold: (1) provide definitive data sets for the assessment of numerical simulations to the Navier-Stokes equations; (2) incorporate advanced instrumentation to measure the spatial and temporal structure of fluid flows; and (3) develop true parallelism between computational and experimental research using the 'scientific workstation' concept. The discussion is presented in viewgraph form.

Davis, Sanford↗

An experimental and analytical study of flow through a supersonic open channel with contoured floor

A wind tunnel experiment was performed to study the characteristics of supersonic airflow (M(infinity) = 2.5-3.86) through an open channel with a contoured floor. The measured static pressures along the centerline of the channel floor exhibited an unexpected rise at the end of the channel. Complex three-dimensional interactions of compression and expansion waves within the channel coupled with external flow perturbations caused by model/tunnel wall interference were the suspected sources of this flow behavior. Three-dimensional inviscid flow analysis procedures were used to investigate and explain this phenomenon. The results of the computations and the experiment are presented and discussed.

Saheli, F. P.↗

Investigating the Dynamic Compaction Behavior of Brittle Powders

The compression of porous materials to full density is a complex physical and mechanical process. Modeling efforts to capture the compaction behavior of distended solids (e.g., porous materials or powders) have produced a variety of compaction models, each functionalized to address specific compaction mechanisms. Of the available models, those suited for high-strain-rate deformation associated with shock loading have generally been developed for ductile materials. Compaction models for brittle materials have only been developed for quasistatic loading conditions. No predictive models currently exist that successfully capture the physical processes and mechanics for shock compaction of brittle distended materials, such as sand and oxide powders. The primary hindrance in producing a dynamic compaction model for brittle distended materials is the lack of availability of high-fidelity, high-precision shock compaction data for validation. The proposed work is aimed at alleviating this issue by developing a high-precision testing method in combination with high-fidelity analysis methods to accurately and precisely determine the shocked state of an ideal test material, ceria powder, under uniaxial strain dynamic compaction. The measured shock compaction data from these brittle powder experiments is then used in initial computations to design a cylindrically converging shock compaction experiment within situ time-resolved shocked density measurements, employing proton radiography. The preliminary experimental and computational results obtained from these uniaxial strain dynamic compaction and cylindrically converging shock compaction experiments are presented. An experimental and computational path forwards to predictive model development is then proposed.

36 MATERIALS SCIENCE↗

Experimental, Theoretical, and Computational Investigation of Separated Nozzle Flows

A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.

Hunter, Craig A.↗

The impact of distributed computing on education

In this paper, developments in digital computer technology since the early Fifties are reviewed briefly, and the parallelism which exists between these developments and developments in analysis and design procedures of structural engineering is identified. The recent trends in digital computer technology are examined in order to establish the fact that distributed processing is now an accepted philosophy for further developments. The impact of this on the analysis and design practices of structural engineering is assessed by first examining these practices from a data processing standpoint to identify the key operations and data bases, and then fitting them to the characteristics of distributed processing. The merits and drawbacks of the present philosophy in educating structural engineers are discussed and projections are made for the industry-academia relations in the distributed processing environment of structural analysis and design. An ongoing experiment of distributed computing in a university environment is described.

Utku, S.↗

Reaction Rate Ratios for Recent Fast Metal Experiments with Large Plutonium Masses

Reaction rate ratios are integral responses that are used within the criticality experiments field because they contain spectral information. While these types of measurements have been utilized for nuclear data validation with historic experiments, few experiments of this type have been utilized for recent experiments, as few exist. This work focuses on measured reaction rate ratios for two nearly bare plutonium critical assemblies with different geometries: one that is cube like (with a Pu mass of 40 kg) and one that is slab like (with a Pu mass of 109 kg). Irradiations were performed with both configurations in which foils were placed near the center of the assembly. Plutonium, highly enriched uranium, depleted uranium, and Au foils were included in the irradiation and counted via high-purity germanium detectors. From these measurements, reaction rate ratios were calculated. Measured and simulated values and uncertainties are presented for the reaction rate ratios. Ratios utilizing the following reactions are given in this work: 197 Au(n, γ), 197 Au(n,2n), 235 U(n,fission), 238 U(n, fission), 238 U(n,2n), 238 U(n,γ), and 239 Pu(n,fission). Uncertainties for the measured reaction rate ratios ranged from 4% to 7%, and the contribution of various parameters to this uncertainty was investigated. The results are compared to historical experiments and should be used for nuclear data validation for future nuclear data library releases. These measurements are part of the EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) project, which utilizes measurement responses in addition to k eff (such as these reaction rate ratios) to help reduce uncertainties in 239 Pu nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of High-Fidelity and Efficient Modeling Capabilities for Enabling Co-Optimization of Fuels and Multi-Mode Engines

The objective of this closely coordinated research project is to develop improved physical sub-models and innovative numerical algorithms to enable accurate and efficient predictions of multi-mode combustion to support EERE’s Co-Optima program. To this end, a research pro-gram is proposed that seeks • to develop accurate submodels for predicting multi-mode combustion regimes, wall-heat transfer, non-equilibrium plasma ignition, and combustion-mode transition, • to develop dynamic load-rebalancing algorithms and efficient semi-implicit time-integration schemes for enabling exascale computing, • to validate the computational submodels against experiments and direct numerical simulation (DNS) data in conjunction with sensitivity analysis, and • to support the direct transition of the resulting sub-models into other engine-simulation tools through the development of self-contained modules. All sub-models will be integrated into the open-source exascale engine-simulation platform Nek5000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dynamics of Nuclear Regions of Galaxies

Current research carried out with the help of the ASEE-NASA Summer Faculty Program, at NASA-Ames, is concentrated on the dynamics of nuclear regions of galaxies. From a dynamical point of view a galaxy is a collection of around 10(sup 11) stars like our Sun, each of which moves in the summed gravitational field of all the remaining stars. Thus galaxy dynamics becomes a self-consistent n-body problem with forces given by Newtonian gravitation. Strong nonlinearity in the gravitational force and the inherent nonlinearity of self-consistent problems both argue for a numerical approach. The technique of numerical experiments consis of constructing an environment in the computer that is as close as possible to the physical conditions in a real galaxy and then carrying out experiments much like laboratory experiments in physics or engineering, in this environment. Computationally, an experiment is an initial value problem, and a good deal of thought and effort goes into the design of the starting conditions that serve as initial values. Experiments are run at Ames because all the 'equipment' is in place-the programs, the necessary computational power, and good facilities for post-run analysis. Our goal for this research program is to study the nuclear regions in detail and this means replacing most of the galaxy by a suitable boundary condition to allow the full capability of numerical experiments to be brought to bear on a small region perhaps 1/1000 of the linear dimensions of an entire galaxy. This is an extremely delicate numerical problem, one in which some small feature overlook, can easily lead to a collapse or blow-up of the entire system. All particles attract each other in gravitational problems, and the 1/r(sup 2) force is: (1) nonlinear; (2) strong at short range; (3) long-range, and (4) unscreened at any distance.

Miller, Richard H.↗

Unlocking hidden information in sparse small-angle neutron scattering measurements

Hypothesis Small-Angle Neutron Scattering (SANS) is a powerful technique for studying soft matter systems such as colloids, polymers, and lyotropic phases, providing nanoscale structural insights. However, its effectiveness is limited by low neutron flux, leading to long acquisition times and noisy data. Here, we hypothesize that Bayesian statistical inference using Gaussian Process Regression (GPR) can reconstruct high-fidelity scattering data from sparse measurements by leveraging intensity smoothness and continuity. Experiments and Simulations The method was benchmarked computationally and validated through SANS experiments on various soft matter systems, including wormlike micelles, colloidal suspensions, polymeric structures, and lyotropic phases. GPR-based inference was applied to both experimental and synthetic data to evaluate its effectiveness in noise reduction and intensity reconstruction. Findings GPR significantly enhances SANS data quality and therefore reducing measurement times by up to two orders of magnitude. This cost-effective approach maximizes experimental efficiency, enabling high-throughput studies and real-time monitoring of dynamic systems. It is particularly beneficial for weakly scattering and time-sensitive studies. Beyond SANS, this framework applies to other low-SNR techniques, including laboratory-based small-angle X-ray scattering and various dynamical scattering methods. Furthermore, it offers transformative potential for compact neutron sources, enhancing their viability for structural analysis in resource-limited settings.

Small angle neutron scattering↗

Reweighting simulated events using machine-learning techniques in the CMS experiment

Data analyses in particle physics rely on an accurate simulation of particle collisions and a detailed simulation of detector effects to extract physics knowledge from the recorded data. Event generators together with a GEANT -based simulation of the detectors are used to produce large samples of simulated events for analysis by the LHC experiments. These simulations come at a high computational cost, where the detector simulation and reconstruction algorithms have the largest CPU demands. This article describes how machine-learning (ML) techniques are used to reweight simulated samples obtained with a given set of parameters to samples with different parameters or samples obtained from entirely different simulation programs. The ML reweighting method avoids the need for simulating the detector response multiple times by incorporating the relevant information in a single sample through event weights. Results are presented for reweighting to model variations and higher-order calculations in simulated top quark pair production at the LHC. This ML-based reweighting is an important element of the future computing model of the CMS experiment and will facilitate precision measurements at the High-Luminosity LHC.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Towards New Metrics for High-Performance Computing Resilience

Ensuring the reliability of applications is becoming an increasingly important challenge as high-performance computing (HPC) systems experience an ever-growing number of faults, errors and failures. While the HPC community has made substantial progress in developing various resilience solutions, it continues to rely on platform-based metrics to quantify application resiliency improvements. The resilience of an HPC application is concerned with the reliability of the application outcome as well as the fault handling efficiency. To understand the scope of impact, effective coverage and performance efficiency of existing and emerging resilience solutions, there is a need for new metrics. In this paper, we develop new ways to quantify resilience that consider both the reliability and the performance characteristics of the solutions from the perspective of HPC applications. As HPC systems continue to evolve in terms of scale and complexity, it is expected that applications will experience various types of faults, errors and failures, which will require applications to apply multiple resilience solutions across the system stack. The proposed metrics are intended to be useful for understanding the combined impact of these solutions on an application's ability to produce correct results and to evaluate their overall impact on an application's performance in the presence of various modes of faults.

Fault Tolerance↗

A Parametric Geometry Computational Fluid Dynamics (CFD) Study Utilizing Design of Experiments (DOE)

Design of Experiments (DOE) was applied to the LAS geometric parameter study to efficiently identify and rank primary contributors to integrated drag over the vehicles ascent trajectory in an order of magnitude fewer CFD configurations thereby reducing computational resources and solution time. SME s were able to gain a better understanding on the underlying flowphysics of different geometric parameter configurations through the identification of interaction effects. An interaction effect, which describes how the effect of one factor changes with respect to the levels of other factors, is often the key to product optimization. A DOE approach emphasizes a sequential approach to learning through successive experimentation to continuously build on previous knowledge. These studies represent a starting point for expanded experimental activities that will eventually cover the entire design space of the vehicle and flight trajectory.

Rhew, Ray D.↗

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Investigation of Improved Wall Heat Flux Behavior Observed When Computing Hypersonic Flow Using a Face-averaged Nodal-gradient Approach On Tetrahedral Grids

We investigate the origin of a significant reduction in computed numerical “noise” in the spatial variation of computed wall heat flux observed when using a node-centered weighted least squares(WLSQ) gradient approach when solving the Reynolds-averaged Navier-Stokes (RANS) equations for hypersonic boundary layer flows on highly stretched tetrahedral grids. This investigation is conducted using a 2nd-order, cell-centered, finite-volume, discretization which we employ to contrast results obtained using two WLSQ gradient approaches: a conventional 3-D, cell-centered, node-neighbor, WLSQ gradient (NN-CCG) approach and a 3-D, face-averaged, node-centered gradient (F-ANG)approach. The F-ANG approach, is shown to produce significantly less “noise” in the spatial variation of wall heat flux, than the NN-CCG approach and we demonstrate the cause of this “noise” via a series of numerical experiments where we compute hypersonic laminar and turbulent boundary layer flows.We further investigate the level of grid resolution required for a tetrahedral grid to approach the fidelity of a grid resolved hexahedral grid solution obtained using an equivalent structured grid method.

Jeffery A. White↗

Convergent energies and anharmonic vibrational spectra of Ca2H2 and Ca2H4 constitutional isomers

Three constitutional isomers of both Ca 2 H 2 and Ca 2 H 4 have been characterized with molecular electronic structure theory. Correlation methods as complete as CCSDT(Q) and basis sets as large as cc-pwCV5Z have been used to converge the relative energies within chemical accuracy (≤1 kcal mol -1 ). Anharmonic vibrational frequencies were computed using second-order vibrational perturbation theory employing CCSD(T)/cc-pwCVTZ cubic and quartic force-fields and a CCSD(T)/cc-pwCVQZ quadratic force field. The monobridged [Ca(μ 2 -H)CaH] and dibridged [Ca(μ 2 -H) 2 Ca] isomers of Ca 2 H 2 were predicted to lie 6.5 and 12.9 kcal mol -1 below the energy of the classical HCaCaH linear isomer, respectively. Despite the energetic favorability of the bridged Ca 2 H 2 isomers, we conclude (surprisingly) that only the higher energy linear structure has been observed in the laboratory. At 0 K, the tribridged [Ca(μ 2 -H) 3 CaH] isomer of Ca 2 H 4 is predicted to be enthalpically favored by 0.9 kcal mol -1 in comparison to the enthalpy of the dibridged [HCa(μ 2 -H) 2 CaH] structure. Comparison of experiment with our computed frequencies suggests that the observed vibrational features arise from both the dibridged and the tribridged Ca 2 H 4 structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclide inventory characterization of EBR-II MOX fuel test pins

Characterization of 3 mixed oxide (MOX) fuel pins from the Experimental Breeder Reactor II (EBR-II) SPA-2/-2B experiment program was performed in support of the Advanced Reactor Experiments for Sodium Fast Reactor Fuels (ARES) joint project between Idaho National Laboratory (INL) and the Japan Atomic Energy Agency (JAEA). The characterization efforts were performed to enable transient testing in the Transient Reactor Test (TREAT) facility at INL and will include post irradiation examination (PIE) and other analyses. MOX fuel pins UW02009, 02011, and 06024 were historically irradiated from 1989 to 1994 for a total of 822.42 effective full power days (EFPDs) and peak burnup of approximately 14.3 at.% (atomic percent of heavy metal atoms, ∼134.4 GWd/t). They had since been placed in storage until a future use was identified. Characterization of the fission products and activation of each fuel pin was necessary to enable removal from storage and subsequent irradiation testing or PIE activities. The initial fuel pin geometry and composition were established using historical documentation. Irradiation history details were utilized to compute irradiated and decayed fuel pin masses, atomic densities, and activities. A pellet-wise axial neutron flux was also computed to assist in evaluation of the axial distribution of {sup 235}U and {sup 239}Pu, per pellet, within each of the test fuel pins. The results computed herein supported experiment design for advanced high-burnup MOX fuel designs under slow transient overpower (MOXTOP) conditions, which was Phase I of the ARES project. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

PETSc Users Manual (Rev. 3.13)

This manual describes the use of PETSc for the numerical solution of partial differential equations and related problems on high-performance computers. The Portable, Extensible Toolkit for Scientific Computation (PETSc) is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all message-passing communication. PETSc includes an expanding suite of parallel linear solvers, nonlinear solvers, and time integrators that may be used in application codes written in Fortran, C, C++, and Python. PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or dbx, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than \rolling them" yourself.

97 MATHEMATICS AND COMPUTING↗