Engineering Papers⌕ Search

SEARCH · Engineering Papers

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,405 records · Page 78

Estimating normal mixture parameters from the distribution of a reduced feature vector

A FORTRAN computer program was written and tested. The measurements consisted of 1000 randomly chosen vectors representing 1, 2, 3, 7, and 10 subclasses in equal portions. In the first experiment, the vectors are computed from the input means and covariances. In the second experiment, the vectors are 16 channel measurements. The starting covariances were constructed as if there were no correlation between separate passes. The biases obtained from each run are listed.

Guseman, L. F.↗

Branch-and-Bound algorithm applied to uncertainty quantification of a Boiling Water Reactor Station Blackout

Evaluation of the impacts of uncertainty and sensitivity in modeling presents a significant set of challenges in particular to high fidelity modeling. Computational costs and validation of models creates a need for cost effective decision making with regards to experiment design. Experiments designed to validate computation models can be used to reduce uncertainty in the physical model. In some cases, large uncertainty in a particular aspect of the model may or may not have a large impact on the final results. For example, modeling of a relief valve may result in large uncertainty, however, the actual effects on final peak clad temperature in a reactor transient may be small and the large uncertainty with respect to valve modeling may be considered acceptable. Additionally, the ability to determine the adequacy of a model and the validation supporting it should be considered within a risk informed framework. Low fidelity modeling with large uncertainty may be considered adequate if the uncertainty is considered acceptable with respect to risk. In other words, models that are used to evaluate the probability of failure should be evaluated more rigorously with the intent of increasing safety margin. Probabilistic risk assessment (PRA) techniques have traditionally been used to identify accident conditions and transients. Traditional classical event tree methods utilize analysts’ knowledge and experience to identify the important timing of events in coordination with thermal-hydraulic modeling. These methods lack the capability to evaluate complex dynamic systems. In these systems, time and energy scales associated with transient events may vary as a function of transition times and energies to arrive at a different physical state. Dynamic PRA (DPRA) methods provide a more rigorous analysis of complex dynamic systems. Unfortunately DPRA methods introduce issues associated with combinatorial explosion of states. This study presents a methodology to address combinatorial explosion using a Branch-and-Bound algorithm applied to Dynamic Event Trees (DET), which utilize LENDIT (L – Length, E – Energy, N – Number, D – Distribution, I – Information, and T – Time) as well as a set theory to describe system, state, resource, and response (S2R2) sets to create bounding functions for the DET. The optimization of the DET in identifying high probability failure branches is extended to create a Phenomenological Identification and Ranking Table (PIRT) methodology to evaluate modeling parameters important to safety of those failure branches that have a high probability of failure. The PIRT can then be used as a tool to identify and evaluate the need for experimental validation of models that have the potential to reduce risk. Finally, in order to demonstrate this methodology, a Boiling Water Reactor (BWR) Station Blackout (SBO) case study is presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantum advantage in learning from experiments

Quantum technology promises to revolutionize how we learn about the physical world. An experiment that processes quantum data with a quantum computer could have substantial advantages over conventional experiments in which quantum states are measured and outcomes are processed with a classical computer. We proved that quantum machines could learn from exponentially fewer experiments than the number required by conventional experiments. This exponential advantage is shown for predicting properties of physical systems, performing quantum principal component analysis, and learning about physical dynamics. Furthermore, the quantum resources needed for achieving an exponential advantage are quite modest in some cases. Conducting experiments with 40 superconducting qubits and 1300 quantum gates, we demonstrated that a substantial quantum advantage is possible with today’s quantum processors.

Science & Technology - Other Topics↗

Scout trajectory error propagation computer program

Since 1969, flight experience has been used as the basis for predicting Scout orbital accuracy. The data used for calculating the accuracy consists of errors in the trajectory parameters (altitude, velocity, etc.) at stage burnout as observed on Scout flights. Approximately 50 sets of errors are used in Monte Carlo analysis to generate error statistics in the trajectory parameters. A covariance matrix is formed which may be propagated in time. The mechanization of this process resulted in computer program Scout Trajectory Error Propagation (STEP) and is described herein. Computer program STEP may be used in conjunction with the Statistical Orbital Analysis Routine to generate accuracy in the orbit parameters (apogee, perigee, inclination, etc.) based upon flight experience.

Myler, T. R.↗

Evaluation of the concrete shield compositions from the 2010 criticality accident alarm system Benchmark experiments at the CEA VALDUC silene facility - 14454

In October 2010, a series of benchmark experiments were conducted at the French Commissariat a l'energie Atomique et aux energies Alternatives (CEA) Valduc SILENE facility. These experiments were a joint effort between the United States Department of Energy Nuclear Criticality Safety Program and the CEA. The purpose of these experiments was to create three benchmarks for the verification and validation of radiation transport codes and evaluated nuclear data used in the analysis of criticality accident alarm systems. This series of experiments consisted of three single-pulsed experiments with the SILENE reactor. For the first experiment, the reactor was bare (unshielded), whereas in the second and third experiments, it was shielded by lead and polyethylene, respectively. The polyethylene shield of the third experiment had a cadmium liner on its internal and external surfaces, which vertically was located near the fuel region of SILENE. During each experiment, several neutron activation foils and thermoluminescent dosimeters (TLDs) were placed around the reactor. Nearly half of the foils and TLDs had additional high-density magnetite concrete, high-density barite concrete, standard concrete, and/or BoroBond shields. CEA Saclay provided all the concrete, and the US Y-12 National Security Complex provided the BoroBond. Measurement data from the experiments were published at the 2011 International Conference on Nuclear Criticality (ICNC 2011) and the 2013 Nuclear Criticality Safety Division (NCSD 2013) topical meeting. Preliminary computational results for the first experiment were presented in the ICNC 2011 paper, which showed poor agreement between the computational results and the measured values of the foils shielded by concrete. Recently the hydrogen content, boron content, and density of these concrete shields were further investigated within the constraints of the previously available data. New computational results for the first experiment are now available that show much better agreement with the measured values. (authors)

42 ENGINEERING↗

Computer program to minimize prediction error in models from experiments with 16 hypercube points and 0 to 6 center points

A previous report described a backward deletion procedure of model selection that was optimized for minimum prediction error and which used a multiparameter combination of the F - distribution and an order statistics distribution of Cochran's. A computer program is described that applies the previously optimized procedure to real data. The use of the program is illustrated by examples.

Holms, A. G.↗

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: 1) a k-epsilon turbulence model computation on a 6.8 million point grid using wall functions, 2) a k-epsilon turbulence model computation on a 14 million point grid integrating to the wall, and 3) a k-omega turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-epsilon and k-omega computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-epsilon computations which extended from 70% to 100% span at the exit rating plane, whereas the k-omega computation had reversed flow from 95% to 100% span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-epsilon case resulted in an under-prediction in adiabatic efficiency of 8.3 points, whereas the k-omega case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: (1) a k-ε turbulence model computation on a 6.8 million point grid using wall functions, (2) a k-ε turbulence model computation on a 14 million point grid integrating to the wall, and (3) a k-ω turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-ε and k-ω computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-ε computations which extended from 70 to 100 percent span at the exit rating plane, whereas the k-ω computation had reversed flow from 95 to 100 percent span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-ε case resulted in an underprediction in adiabatic efficiency of 8.3 points, whereas the k-ω case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

Constraining computational modeling of indirect drive double shell capsule implosions using experiments

Double shell capsule implosions are an alternative approach to achieving alpha heating on the National Ignition Facility. Current machining techniques construct the outer shell as two hemispheres that are glued together, and the deuterium and tritium (DT) liquid inside the inner shell will be injected by a fill tube. These features introduce asymmetries and jetting that may disrupt the confinement of the DT fuel if not carefully controlled. Simulations indicate that in order to achieve high yields in the laboratory, these features as well as susceptibility to the Rayleigh–Taylor instability (RTI) must be mitigated. Furthermore, due to uncertainties in computational models and the expense of using the best physics models at adequate resolution in three dimensions, our computational modeling must be constrained by experiments. We report on the results of recent hydrogrowth radiography and dual-axis keyhole experiments with double shell targets that have been used to evaluate our modeling of the outer shell joint as well as the impacts of high-energy x-ray preheat that strongly impacts RTI growth. Our simulations show good agreement with the experimental data and inform several important modeling choices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled EM-PIC/Radiation Transport Simulations of HERMES Courtyard Experiments

A suite of coupled computational models for simulating the radiation, plasma, and electromagnetic (EM) environment in the High-Energy Radiation Megavolt Electron Source (HERMES) courtyard has been developed. In principle, this provides a predictive forward-simulation capability based solely on measured upstream anode and cathode current waveforms in the Magnetically Insulated Transmission Line (MITL). First, 2D R-Z ElectroMagnetic Particle-in-Cell (EM-PIC) simulations model the MITL and diode to compute a history of all electrons incident on the converter. Next, radiation transport simulations use these electrons as a source to compute the time-dependent dose rate and volumetric electron production in the courtyard. Finally, the radiation transport output is used as sources for EM-PIC simulations of the courtyard to com- pute electromagnetic responses. This suite has been applied to the November 2016 trials, shots 10268-10313. Modeling and experiment differ in significant ways. This is just the first iteration of a long process to improve the agreement, as outlined in the summary.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Collaboration with NSTX in Calculations of Radiofrequency and Neutral Beam Heating and Current Drive

The overall objective of this project was to advance the understanding of heating and current drive necessary for toroidal magnetic confinement fusion energy reactors, and to provide accurate models for these sources to aid in understanding magneto-hydrodynamic, transport and boundary physics issues. Increased understanding is provided by computer simulation using the underlying physics laws, and will lead to greater reliability and increased efficiency of a power producing fusion reactor. The proposed work entailed applications of existing CompX Fokker-Planck, ray tracing, full-wave microwave, and particle diffusion calculator codes to the National Spherical Torus Experiment at Princeton Plasma Physics Laboratory. Special features of NSTX include that large ion orbits are emphasized in its geometry and thus offer a major opportunity within fusion research for focused exploration of their ubiquitous effects, including comparison with simulation codes. The CompX CQL3D Fokker-Planck code was be further adapted to the finite-orbit-width physics. Comparison with experiment ensured that the computer codes have sufficient accuracy to accurately model the experiments. Nonthermal distributions of energetic ions are produced in experiments-- and modeled with the codes-- exhibiting the signatures of auxiliary heating, transport, and radiation processes within the plasma. Using physics-based modeling, in conjunction with experimental data, nonthermal distributions can provide keys to validation and quantification of auxiliary heating physics. Nonthermal distributions also hold particular information on radial transport processes, not obscured by the smoothing Maxwellization effects of Coulomb collisions. Generally, radial transport is not well understood but is a key to efficient fusion energy production. By comparison with Fokker-Planck modeling, the velocity dependence of radial transport coefficients becomes experimentally accessible, particularly in cases of well-modeled sources such as neutral beams, fast wave power, or the toroidal electric field. It is clear from modeling of fast-wave interactions with ions that finite-orbit-width effects play a significant role. For example, a different set of ions is computed to intersect a given cyclotron resonance depending on whether or not radial orbit drifts are taken into account. Also, radial transport, associated bootstrap current, and alpha channeling depend on radial drifts. Ion cyclotron heating above the first harmonic interacts preferentially with higher velocity ions, giving large radial drifts. The CQL3D Fokker-Planck code, maintained by the proposers, is a physics based model used extensively within the US fusion community to compare with and interpret experimental observations. At the beginning of this grant, bounce-averages of the rf quasilinear and other operators have been performed neglecting the radial drifts which give finite orbit widths to the particle orbits. The proposed work included a major upgrade to CQL3D to include finite-orbit-width effects be accounted for by performing the bounce-averages of rf and neutral beam fast ions along guiding center orbits. This CQL3D enhancement provided microwave and neutral beam contributions to neoclassical radial transport coefficients for nonthermal distributions, and substantially increased the comprehensiveness of the ion model. Additional work validated Fokker-Planck and ray tracing modeling of neutral beam and high harmonic fast wave absorption and current drive against the NSTX experiment, as had been very successful for electron cyclotron current drive in the DIII-D tokamak at General Atomics and lower hybrid waves on C-Mod at MIT. The validation work was in conjunction with staff and graduate students at Princeton Plasma Physics Laboratory. This DOE grant supported fusion energy research, a potential long-term solution to the world's energy needs. The research grant under consideration, DE-SC0006614, provided research support at the 0.3 Full-Time-Equivalent level for application of computer models to aid in understanding and projecting efficacy of heating and current drive sources in the National Spherical Torus Experiment, a tokamak variant, at the Princeton Plasma Physics Laboratory. The NSTX experiment explores the physics of very tight aspect ratio, almost spherical tokamaks, aiming at producing steady-state fusion plasmas. The current drive is an integral part of the steady-state concept, maintaining the magnetic geometry in the steady-state tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental Data Sets for CFD Calculations of Flow and Heat Transfer in Deformed Fuel Assemblies

Liquid-metal fast-breeder reactors (LMFBRs) have traditionally used fuel assemblies consisting of a hexagonal array of pins enclosed in a hexagonally shaped duct. Each fuel pin is surrounded by a wire that is wrapped in a helical configuration, which is used to separate the pins sufficiently to allow sodium coolant to flow through the assembly. The helical wire wrap also promotes mixing of the coolant. Experience with these assemblies has shown that they begin to deform with time during normal operation leading to significant changes in core thermal-hydraulics. For example, deformations of assemblies in the Experimental Breeder Reactor-2 and the Fast Flux Test Facility led to pin overheating resulting from pin-to-duct and pin-to-pin contact. This problem is well known and was investigated extensively in the 1970's and early 1980's. This was a period, however, before the development of modern Computational Fluid Dynamics (CFD) computer codes. Therefore, the experiments performed at the time did not record their results in sufficient detail to be useful for compiling validation datasets for CFD codes. Most of the experiments focused only on bulk measurements, such as the pressure drop across the bundle. When the effects of deformation of the fuel assembly were investigated, the studies performed to quantify these effects considered only bare pins without the wire wrap. With the resurgence of interest in LMFBRs in recent years, particularly in designs that leverage long-life performance of the fuel, the need has arisen for data that are suitable for the validation and benchmarking of CFD codes to predict the thermal-hydraulic behavior of flow in wire-wrapped fuel bundles. Such data must include high-resolution, transient, full-field measurements of the fluid velocity and temperature and an accurate record of the test conditions, including a precise description of the geometry. This paper describes a project to produce a validation-quality set of data that can be used to benchmark CFD codes for modeling flow and thermal performance in wire-wrapped fuel assemblies of the kind used by LMFBR designs. In addition, this effort investigates the changes to the flow and other phenomena that occur in the deformed geometries that can form later in a fuel assembly's life. This work is a joint effort of four organizations. The experimental work is being performed by AREVA and Texas A and M University. The computational work, both to support the experiments and to perform blind benchmark calculations, is being done by TerraPower, LLC, and Argonne National Laboratory (ANL). The experiments use a geometry that consists of a complete assembly of 61 wire-wrapped pins, which is surrounded by a hexagonal duct. This configuration is considered the minimum size necessary to adequately capture the important results of interest, such as the thermal gradients across the assembly and the velocity field effects in the assembly resulting from the surrounding duct. These types of effects cannot be captured by a geometry consisting of only a small group of fuel pins. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Tackling the Reproducibility Problem in Systems Research with Declarative Experiment Specifications

Validating experimental results in the field of computer systems is a challenging task, mainly due to the many changes in software and hardware that computational environments go through. Determining if an experiment is reproducible entails two separate tasks: re-executing the experiment and validating the results. Existing reproducibility efforts have focused on the former, envisioning techniques and infrastructures that make it easier to re-execute an experiment. In this work we focus on the latter by analyzing the validation workflow that an experiment re-executioner goes through. We notice that validating results is done on the basis of experiment design and high-level goals, rather than exact quantitative metrics. Based on this insight, we introduce a declarative format for specifying the high-level components of an experiment as well as describing generic, testable conditions that serve as the basis for validation. We present a use case in the area of storage systems to illustrate the usefulness of this approach. We also discuss limitations and potential benefits of using this approach in other areas of experimental systems research.

97 MATHEMATICS AND COMPUTING↗

Mixing of multiple jets with a confined subsonic crossflow - Summary of NASA-supported experiments and modeling

Experimental and computational results on the mixing of single, double, and opposed rows of jets with an isothermal or variable temperature mainstream in a confined subsonic crossflow are summarized. The studies were performed to investigate flow and geometric variations typical of the complex 3D flowfield in the dilution zone of combustion chambers in gas turbine engines. The principal observations from the experiments were that the momentum-flux ratio was the most significant flow variable, and that temperature distributions were similar (independent of orifice diameter) when the orifice spacing and the square-root of the momentum-flux ratio were inversely proportional. The experiments and empirical model for the mixing of a single row of jets from round holes were extended to include several variations typical of gas turbine combustors.

Holdeman, James D.↗

Computers vs. wind tunnels for aerodynamic flow simulations

It is pointed out that in other fields of computational physics, such as ballistics, celestial mechanics, and neutronics, computations have already displaced experiments as the principal means of obtaining dynamic simulations. In the case of aerodynamic investigations, the complexity of the computational work involved in solving the Navier-Stokes equations is the reason that such investigations rely currently mainly on wind-tunnel testing. However, because of inherent limitations of the wind-tunnel approach and economic considerations, it appears that at some time in the future aerodynamic studies will chiefly rely on computational flow data provided by the computer. Taking into account projected development trends, it is estimated that computers with the required capabilities for a solution of the complete viscous, time-dependent Navier-Stokes equations will be available in the mid-1980s.

Chapman, D. R.↗

Mixing of multiple jets with a confined subsonic crossflow. Summary of NASA-supported experiments and modeling

Experimental and computational results on the mixing of single, double, and opposed rows of jets with an isothermal or variable temperature mainstream in a confined subsonic crossflow are summarized. The studies were performed to investigate flow and geometric variations typical of the complex 3-D flowfield in the dilution zone of combustion chambers in gas turbine engines. The principal observations from the experiments were that the momentum-flux ratio was the most significant flow variable, and that temperature distributions were similar (independent of orifice diameter) when the orifice spacing and the square-root of the momentum-flux ratio were inversely proportional. The experiments and empirical model for the mixing of a single row of jets from round holes were extended to include several variations typical of gas turbine combustors. Combinations of flow and geometry that gave optimum mixing were identified from the experimental results. Based on results of calculations made with a 3-D numerical model, the empirical model was further extended to model the effects of curvature and convergence. The principle conclusions from this study were that the orifice spacing and momentum-flux relationships were the same as observed previously in a straight duct, but the jet structure was significantly different for jets injected from the inner wall wall of a turn than for those injected from the outer wall. Also, curvature in the axial direction caused a drift of the jet trajectories toward the inner wall, but the mixing in a turning and converging channel did not seem to be inhibited by the convergence, independent of whether the convergence was radial or circumferential. The calculated jet penetration and mixing in an annulus were similar to those in a rectangular duct when the orifice spacing was specified at the radius dividing the annulus into equal areas.

Holdeman, James D.↗

Testing and Simulation of an Updated Cartridge Loop Vehicle

The Versatile Test Reactor (VTR) is a sodium-cooled, fast-spectrum test reactor that is being developed in the United States and will support a variety of irradiation test vehicle configurations, including cartridge loops. This work includes out-of-pile experimental results from a single-phase, natural circulation cartridge loop vehicle with geometry relevant to VTR irradiation sites, as well as comparisons between the experimental results and results predicted using the TRAC/RELAP Advanced Computational Engine (TRACE) modeling tool. The experiments were conducted in the thermosyphon test loop (TSTL) facility at Oak Ridge National Laboratory. Comparisons are also made between the current experimental data and results from natural circulation experiments previously conducted in the TSTL in a cartridge vehicle that is similar in design but has smaller flow areas. This cartridge vehicle and the experimental program were developed to add to the single-phase, natural circulation data collected in the previous iteration of the cartridge loop design, which supports future irradiation experiments and adds to a database that is useful for validating computer models. Comparisons of experimental results to TRACE model predictions is a pertinent step in validating the computational tool for supporting future irradiation experiment design and safety calculations, and comparisons to previous cartridge loop results highlight the impact of the design changes made to the test vehicle. The experiments conducted include several steady state tests and transients, including power ramp, loss of offsite power, and loss of external flow scenarios. This work shows that TRACE can accurately predict temperatures and flow conditions in the cartridge loop and the updated vehicle design achieves higher mass flow rates at the same steady state power levels.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗