SEARCH · Engineering Papers

Results for “simulation methods and programs”

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.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Recent developments of the NESSUS probabilistic structural analysis computer program

The NESSUS probabilistic structural analysis computer program combines state-of-the-art probabilistic algorithms with general purpose structural analysis methods to compute the probabilistic response and the reliability of engineering structures. Uncertainty in loading, material properties, geometry, boundary conditions and initial conditions can be simulated. The structural analysis methods include nonlinear finite element and boundary element methods. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. The scope of the code has recently been expanded to include probabilistic life and fatigue prediction of structures in terms of component and system reliability and risk analysis of structures considering cost of failure. The code is currently being extended to structural reliability considering progressive crack propagation. Several examples are presented to demonstrate the new capabilities.

Millwater, H.↗

NASA NTRS1960

Problems Involved in an Emergency Method of Guiding a Gliding Vehicle from High Altitudes to a High Key Position

An investigation has been conducted to determine the problems involved in an emergency method of guiding a gliding vehicle from high altitudes to a high key position (initial position) above a landing field. A jet airplane in a simulated flameout condition, conventional ground-tracking radar, and a scaled wire for guidance programming on the radar plotting board were used in the tests. Starting test altitudes varied from 30,000 feet to 46,500 feet, and starting positions ranged 8.4 to 67 nautical miles from the high key. Specified altitudes of the high key were 12,000, 10,000 or 4,000 feet. Lift-drag ratios of the aircraft of either 17, 16, or 6 were held constant during any given flight; however, for a few flights the lift-drag ratio was varied from 11 to 6. Indicated airspeeds were held constant at either 160 or 250 knots. Results from these tests indicate that a gliding vehicle having a lift-drag ratio of 16 and an indicated approach speed of 160 knots can be guided to within 800 feet vertically and 2,400 feet laterally of a high key position. When the lift-drag ratio of the vehicle is reduced to 6 and the indicated approach speed is raised to 250 knots, the radar controller was able to guide the vehicle to within 2,400 feet vertically and au feet laterally of the high key. It was also found that radar stations which give only azimuth-distance information could control the glide path of a gliding vehicle as well as stations that receive azimuth-distance-altitude information, provided that altitude information is supplied by the pilot.

Jewel, Joseph W., Jr.↗

NASA NTRS2000

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗

DOE OSTI2020

Comparison of Generation of Higher-Order Neutron Scattering Cross Sections

The generation of high-order neutron scattering cross sections consistent with high-fidelity simulations remains an area of active research. Popular options include generating cross sections from continuous energy Monte Carlo calculations or generating cross sections from a deterministic neutron transport calculation with high-fidelity tabulated cross sections. Both options present challenges. Monte Carlo simulations can naturally process continuous energy cross section data and allow for the general description of anisotropic neutron scattering given a scattering law. However, Monte Carlo simulations are inherently related to particle weighting and it has been suggested that this may be unacceptable for generating high-order neutron scattering cross sections. Deterministic neutron transport calculations can easily calculate high-order moments of the flux to appropriately calculate higher-order neutron scattering cross sections but are generally limited by discretization of space, energy, and angle. In this work, the trade-offs between generation of high-order neutron scattering cross sections via Monte Carlo and deterministic neutron transport methods are investigated. The methods implemented in the Monte Carlo computer program Serpent and the deterministic fast reactor neutron cross section generator MC2-3 are compared. Cross sections resulting from these methods are used in Rattlesnake, a deterministic neutron transport code developed by Idaho National Laboratory (INL), and results are compared to a reference continuous energy Monte Carlo calculation. Whereas previous work investigating the effects of anisotropic neutron scattering has focused on light water reactor (LWR) simulations, this work focuses on high-order neutron scattering cross sections as they relate to fast reactor simulations. To investigate the consequences of the Serpent and MC2-3 methodologies, a test problem is developed. The test problem is a one-dimensional geometry with fast reactor materials designed to demonstrate deep penetration and exacerbate the effects of high-order neutron scattering. Based on the results of the deep penetration test problem, it is concluded that P3 neutron scattering cross sections are sufficient to describe anisotropic scattering in fast reactor materials. Any neutron scattering of order higher than P3 offers negligible change in the eigenvalue of the test problem. Additionally, it is determined that the methodology as implemented in Serpent is applicable for generating high-order neutron scattering cross sections through at least P3 in fast reactor materials.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Application of Quantum Machine Learning to High Energy Physics Analysis at LHC Using Quantum Computer Simulators and Quantum Computer Hardware

Machine learning enjoys widespread success in High Energy Physics (HEP) analyses at LHC. However the ambitious HL-LHC program will require much more computing resources in the next two decades. Quantum computing may offer speed-up for HEP physics analyses at HL-LHC, and can be a new computational paradigm for big data analyses in High Energy Physics.We have successfully employed three methods (1) Variational Quantum Classifier (VQC) method, (2) Quantum Support Vector Machine Kernel (QSVM-kernel) method and (3) Quantum Neural Network (QNN) method for two LHC flagship analyses: ttH (Higgs production in association with two top quarks) and H->mumu (Higgs decay to two muons, the second generation fermions). We shall address the progressive improvements in performance from method (1) to method (3).We will present our experiences and results of a study on LHC High Energy Physics data analyses with IBM Quantum Simulator and Quantum Hardware (using IBM Qiskit framework), Google Quantum Simulator (using Google Cirq framework), and Amazon Quantum Simulator (using Amazon Braket cloud service). The work is in the context of a Qubit platform (a gate-model quantum computer). Taking into account the present limitation of hardware access, different quantum machine learning methods are studied on simulators and the results are compared with classical machine learning methods (BDT, classical Support Vector Machine and classical Neural Network). Furthermore, we do apply quantum machine learning on IBM quantum hardware to compare performance between quantum simulator and quantum hardware. The work is performed by an international and interdisciplinary collaboration with the Department of Physics and Department of Computer Sciences of University of Wisconsin, CERN Quantum Technology Initiative, IBM Research Zurich, IBM T.J. Watson Research Center, Fermilab Quantum Institute, BNL Computational Science Initiative, State University of New York at Stony Brook, and Quantum Computing and AI Research of Amazon Web Services. This work pioneers a close collaboration of academic institutions with industrial corporations in the High Energy Physics analyses effort. Though the size of event samples in future HL-LHC physics and the limited number of qubits pose some challenges to the Quantum Machine learning studies for High Energy Physics, more advanced quantum computers with larger number of qubits, reduced noise and improved running time (as envisioned by IBM and Google) may outperform classical machine learning in both classification power and in speed.Although the era of efficient quantum computing may still be years away, we have made promising progress and obtained preliminary results in applying quantum machine learning to High Energy Physics. A PROOF OF PRINCIPLE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1973

Computer simulation of shock wave boundary layer interaction

Development of a boundary layer correction for an analysis program previously used by Bauer et al. (1972) to study two-dimensional flow past transonic airfoils at off-design conditions where shocks which interact with the boundary layer appear. Using a method which combines conformal mapping with a finite-difference scheme due to Murman and Cole (1971), an analysis is made of the interaction between shock waves on a supercritical wing section and a turbulent boundary layer, assuming that the shocks are weak and that separation is insignificant. It is shown that relatively simple shock wave/boundary layer interactions of the kind envisioned in this study can be treated by determining the displacement thickness iteratively together with the flow and altering the profile accordingly. It is concluded that the proposed boundary layer correction appears to furnish an essential improvement in cases where both pronounced aft loading and a large enclosed supersonic region of flow are present.

Bauer, F.↗

DOE OSTI2025

Data Selection Improvement For MicroBooNE

Data selection is an extremely important part of data analysis for any experiment. Finding a physics result is often the result of sifting through a massive amount of data, keeping data that we believe to be signal, and throwing out data we do not. This process is called data selection. Creating a selection algorithm is an intensive process that must balance keeping enough data to have statistics and maximizing the signal purity of that data. We also need to choose the right reconstruction method, a tool to take raw data from the detector and convert it into physics results. In this study, we used three different reconstruction tools, Pandora, WireCell, and LANTERN, for the MicroBooNE experiment in conjunction to improve the selection algorithm for analysis. For the case of this study, we look into the charged current N proton 0 pions (CCNp0$\pi$) interaction channel. This is the dominant channel for the Short Baseline Neutrino (SBN) program and is expected to be a large contributor to the Deep Underground Neutrino Experiment (DUNE). We first investigated each of the three tools to find out more about their strengths and weaknesses as reconstructions, and compared them to the truth information directly from the MicroBooNE simulation pipeline. We then put together a direct comparison of the three methods to find which method or combination of methods would return the best result for us. While the study is ongoing, we have learned a lot about data selection for the experiment and the differences between the reconstruction tools.

Dillon, Brayden [Fermilab]↗

NASA NTRS2014

HRP's Healthcare Spin-Offs Through Computational Modeling and Simulation Practice Methodologies

Spaceflight missions expose astronauts to novel operational and environmental conditions that pose health risks that are currently not well understood, and perhaps unanticipated. Furthermore, given the limited number of humans that have flown in long duration missions and beyond low Earth-orbit, the amount of research and clinical data necessary to predict and mitigate these health and performance risks are limited. Consequently, NASA's Human Research Program (HRP) conducts research and develops advanced methods and tools to predict, assess, and mitigate potential hazards to the health of astronauts. In this light, NASA has explored the possibility of leveraging computational modeling since the 1970s as a means to elucidate the physiologic risks of spaceflight and develop countermeasures. Since that time, substantial progress has been realized in this arena through a number of HRP funded activates such as the Digital Astronaut Project (DAP) and the Integrated Medical Model (IMM). Much of this success can be attributed to HRP's endeavor to establish rigorous verification, validation, and credibility (VV&C) processes that ensure computational models and simulations (M&S) are sufficiently credible to address issues within their intended scope. This presentation summarizes HRP's activities in credibility of modeling and simulation, in particular through its outreach to the community of modeling and simulation practitioners. METHODS: The HRP requires all M&S that can have moderate to high impact on crew health or mission success must be vetted in accordance to NASA Standard for Models and Simulations, NASA-STD-7009 (7009) [5]. As this standard mostly focuses on engineering systems, the IMM and DAP have invested substantial efforts to adapt the processes established in this standard for their application to biological M&S, which is more prevalent in human health and performance (HHP) and space biomedical research and operations [6,7]. These methods have also generated substantial interest by the broader medical community though institutions like the National Institutes of Health (NIH) and the Food and Drug Administration (FDA) to develop similar standards and guidelines applicable to the larger medical operations and research community. DISCUSSION: Similar to NASA, many leading government agencies, health institutions and medical product developers around the world are recognizing the potential of computational M&S to support clinical research and decision making. In this light, substantial investments are being made in computational medicine and notable discoveries are being realized [8]. However, there is a lack of broadly applicable practice guidance for the development and implementation of M&S in clinical care and research in a manner that instills confidence among medical practitioners and biological researchers [9,10]. In this presentation, we will give an overview on how HRP is working with the NIH's Interagency Modeling and Analysis Group (IMAG), the FDA and the American Society of Mechanical Engineers (ASME) to leverage NASA's biomedical VV&C processes to establish a new regulatory standard for Verification and Validation in Computational Modeling of Medical Devices, and Guidelines for Credible Practice of Computational Modeling and Simulation in Healthcare.

Mulugeta, Lealem↗

DOE OSTI2020

Designing a drone delivery network with automated battery swapping machines

Drones are projected to alter last-mile delivery, but their short travel range is a concern. In this study, we propose a drone delivery network design using automated battery swapping machines (ABSMs) to extend ranges. The design minimizes the long-term delivery costs, including ABSM investment, drone ownership, and cost of the delivery time, and locates ABSMs to serve a set of customers. We build a mixed-integer nonlinear program that captures the nonlinear waiting time of drones at ABSMs. To solve the problem, we create an exact solution algorithm that finds the globally optimal solution using a derivative-supported cutting-plane method. To validate the applicability of our program, we conduct a case study on the Chicago Metropolitan area using cost data from leading ABSM manufacturer and geographical data from the planning and operations language for agent-based regional integrated simulation (more commonly known as POLARIS). A sensitivity analysis identifies that ABSM service times and costs are the key parameters impacting the long-term adoption of drone delivery.

25 ENERGY STORAGE↗

NASA NTRS1991

Probabilistic structural analysis algorithm development for computational efficiency

The PSAM (Probabilistic Structural Analysis Methods) program is developing a probabilistic structural risk assessment capability for the SSME components. An advanced probabilistic structural analysis software system, NESSUS (Numerical Evaluation of Stochastic Structures Under Stress), is being developed as part of the PSAM effort to accurately simulate stochastic structures operating under severe random loading conditions. One of the challenges in developing the NESSUS system is the development of the probabilistic algorithms that provide both efficiency and accuracy. The main probability algorithms developed and implemented in the NESSUS system are efficient, but approximate in nature. In the last six years, the algorithms have improved very significantly.

Wu, Y.-T.↗

NASA NTRS1996

Particle/Continuum Hybrid Simulation in a Parallel Computing Environment

The objective of this study was to modify an existing parallel particle code based on the direct simulation Monte Carlo (DSMC) method to include a Navier-Stokes (NS) calculation so that a hybrid solution could be developed. In carrying out this work, it was determined that the following five issues had to be addressed before extensive program development of a three dimensional capability was pursued: (1) find a set of one-sided kinetic fluxes that are fully compatible with the DSMC method, (2) develop a finite volume scheme to make use of these one-sided kinetic fluxes, (3) make use of the one-sided kinetic fluxes together with DSMC type boundary conditions at a material surface so that velocity slip and temperature slip arise naturally for near-continuum conditions, (4) find a suitable sampling scheme so that the values of the one-sided fluxes predicted by the NS solution at an interface between the two domains can be converted into the correct distribution of particles to be introduced into the DSMC domain, (5) carry out a suitable number of tests to confirm that the developed concepts are valid, individually and in concert for a hybrid scheme.

Baganoff, Donald↗

NASA NTRS1999

Sensor Lead Wires Positioned on SiC-based Monolithic Ceramic and Fiber- reinforced Ceramic Matrix Composite Subcomponents with Flat and Curved Surfaces

There is strong interest in the development of silicon carbide-based monolithic ceramic and composite materials and components for demanding, high-temperature applications. Thorough characterization of material properties, including high-temperature testing under simulated or actual operating conditions, is a high priority for programs involved in developing these silicon carbide- (SiC) based materials and components. Members of the Sensors and Electronics Technology Branch at the NASA Lewis Research Center are developing minimally intrusive methods of measuring the properties (such as the surface temperature, strain, and heat flux characteristics) of components and subelements that are being tested or operated in hostile, high-temperature environments. Their primary goal is to instrument the test article or operating component with durable sensors that have a minimal effect on test conditions such as the gas flow across the surface of the item and the material response (including the through-thickness conduction of heat). Therefore, the main thrust of their work has been the development of thin-film sensors (e.g., thermocouples or strain gauges) for use on various advanced material test articles, including SiC/SiC composite components. There was a need for a better method of securing sensor lead wires on SiC-based components and subelements that would be tested at temperatures to 1000 C (or higher), to enhance the durability of the overall minimally intrusive sensor system. To address this need, Lewis researchers devised an alternative approach for positioning the sensor lead wires (which are connected to the thin-film sensors) on SiC or SiC/SiC components. A reaction-forming method of joining was used to strongly bond hoop-shaped monolithic SiC and SiC/SiC composite attachments of various sizes to both flat and curved surfaces of SiC/SiC composite subelements (see the photos). This approach is based on an affordable, robust ceramic joining technology, named ARCJoinT, which was developed at Lewis for the joining of SiC-based ceramics and fiber-reinforced composites.

Kiser, James D.↗

NASA NTRS2014

Human Performance in Simulated Reduced Gravity Environments

NASA is currently designing a new space suit capable of working in deep space and on Mars. Designing a suit is very difficult and often requires trade‐offs between performance, cost, mass, and system complexity. Our current understanding of human performance in reduced gravity in a planetary environment (the moon or Mars) is limited to lunar observations, studies from the Apollo program, and recent suit tests conducted at JSC using reduced gravity simulators. This study will look at our most recent reduced gravity simulations performed on the new Active Response Gravity Offload System (ARGOS) compared to the C‐9 reduced gravity plane. Methods: Subjects ambulated in reduced gravity analogs to obtain a baseline for human performance. Subjects were tested in lunar gravity (1.6 m/sq s) and Earth gravity (9.8 m/sq s) in shirt‐sleeves. Subjects ambulated over ground at prescribed speeds on the ARGOS, but ambulated at a self‐selected speed on the C‐9 due to time limitations. Subjects on the ARGOS were given over 3 minutes to acclimate to the different conditions before data was collected. Nine healthy subjects were tested in the ARGOS (6 males, 3 females, 79.5 +/- 15.7 kg), while six subjects were tested on the C‐9 (6 males, 78.8 +/- 11.2 kg). Data was collected with an optical motion capture system (Vicon, Oxford, UK) and was analyzed using customized analysis scripts in BodyBuilder (Vicon, Oxford, UK) and MATLAB (MathWorks, Natick, MA, USA). Results: In all offloaded conditions, variation between subjects increased compared to 1‐g. Kinematics in the ARGOS at lunar gravity resembled earth gravity ambulation more closely than the C‐9 ambulation. Toe‐off occurred 10% earlier in both reduced gravity environments compared to earth gravity, shortening the stance phase. Likewise, ankle, knee, and hip angles remained consistently flexed and had reduced peaks compared to earth gravity. Ground reaction forces in lunar gravity (normalized to Earth body weight) were 0.4 +/- 0.2 on the ARGOS, but only 0.2 +/- 0.1 on the C‐9. Discussion: Gait analysis showed differences in joint kinematics and temporal‐spatial parameters between the reduced gravity simulators and with respect to earth gravity. Although most of the subjects chose a somewhat unique ambulation style as a result of learning to ambulate in a new environment, all but two were consistent with keeping an Earth‐like gait. Learning how reduced gravity affects ambulation will help NASA to determine optimal suit designs, influence mission planning, help train crew, and may shed light on the underlying methods the body uses to optimize gait for energetic efficiency. Conclusion: Kinematic and kinetic analysis demonstrated noteworthy differences between an offloaded environment and 1‐g, as would be expected. The analysis showed a trend to change the ambulation style in an offloaded environment to a rolling‐loping walk (resembling crosscountry skiing) with increased swing time. This ambulation modification, particularly in the ARGOS, indicated that the relative kinetic energy of the subject was increased, on average, per the static body weight compared to the 1‐g condition. How much of this was influenced by the active offloading of the ARGOS system is unknown.

Cowley, Matthew↗

NASA NTRSDate not supplied

Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools

The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.

pressurization↗

DOE OSTI2025

Comparing Control Performance Between Simulation and Experiment using the Microreactor Automated Control System Testbed

In the advanced reactor domain, a flexible and scalable software/hardware infrastructure is crucial for integrating and validating various control technologies. This study used the Microreactor Automated Control System (MACS) hardware platform as a testbed. MACS was originally designed to mirror Idaho National Laboratory (INL)'s Microreactor Applications Research Validation and Evaluation (MARVEL), a 85-kW thermal fission microreactor. It features control drums for simulated reactivity control; lights that function as a surrogate reactor core, with the brightness being proportional to the reactor power; and light sensors that emulate neutron detectors. To transform MACS into a physical twin of MARVEL for evaluating control methods, the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software was employed. This software integrated the hardware with two models of the MARVEL core, based on Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The study aimed to demonstrate the gap between control theory and actual practice—a gap that often necessitates empirical adjustments such as control gain retuning, filters, time discretization, and integrator anti-windup measures. Controllers were developed based on increasingly complex simulations without hardware, starting from the base MARVEL model and then introducing actuator saturation constraints and sensor noise. The final control strategy was then tested using MACS, and a comparative performance analysis was conducted.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

NASA NTRS2010

Measurement of Insulation Compaction in the Cryogenic Fuel Tanks at Kennedy Space Center by Fast/Thermal Neutron Techniques

The liquid hydrogen and oxygen cryogenic storage tanks at John F. Kennedy Space Center (KSC) use expanded perlite as thermal insulation. Th ere is evidence that some of the perlite has compacted over time, com promising the thermal performance and possibly also structural integr ity of the tanks. Therefore an Non-destructive Testing (NDT) method for measuring the perlite density or void fraction is urgently needed. Methods based on neutrons are good candidates because they can readil y penetrate through the 1.75 cm outer steel shell and through the ent ire 120 cm thickness of the perlite zone. Neutrons interact with the nuclei of materials to produce characteristic gamma rays which are the n detected. The gamma ray signal strength is proportional to the atom ic number density. Consequently, if the perlite is compacted then the count rates in the individual peaks in the gamma ray spectrum will i ncrease. Perlite is a feldspathic volcanic rock made up of the major elements Si, AI, Na, K and 0 along with some water. With commercially available portable neutron generators it is possible to produce simul taneously fluxes of neutrons in two energy ranges: fast (14 MeV) and thermal (25 meV). Fast neutrons produce gamma rays by inelastic scatt ering which is sensitive to Fe and O. Thermal neutrons produce gamma rays by radiative capture in prompt gamma neutron activation (PGNA) and this is sensitive to Si, AI, Na, Kand H. Thus the two energy ranges produce complementary information. The R&D program has three phases: numerical simulations of neutron and gamma ray transport with MCNP s oftware, evaluation of the system in the laboratory on test articles and finally mapping of the perlite density in the cryogenic tanks at KSC. The preliminary MCNP calculations have shown that the fast/therma l neutron NDT method is capable of distinguishing between expanded an d compacted perlite with excellent statistics.

Livingston, R. A.↗

DOE OSTI2023

Microstructural Simulation of Ion-Irradiated Natural Rocks and Minerals

This report details the contribution of Oak Ridge National Laboratory to the Nuclear Energy University Program (NEUP) project, "Rapid Characterization of Concrete Mineralogy using Multi-Scale Tools" led by the University of Illinois at Urbana-Champaign. To support the subsequent license renewal (SLR) of the US nuclear power plants (NPPs) fleet, the performance of concrete-forming aggregates against neutron irradiation needs to be assessed. Ion irradiation is proposed as a rapid, cost-effective surrogate method to full neutron irradiation testing. To complement the characterization of ion-irradiated rock specimens, ORNL has run finite-elements and fast-Fourier transform (FFT-based) simulations using the codes MARS and Microstructure-Oriented Scientific Analysis of Irradiated Concrete (MOSAIC). The main conclusion of this analysis is that the apparent post–ion-irradiation step-height underestimates the accumulated free radiation-induced volumetric expansion (RIVE) in the ion-implanted depth by about 15% at full amorphization and about 25% toward the beginning of the ion irradiation experiment. The discrepancy is explained by the fact that the step height is proportional to the sum of the RIVE (i.e., one third of the RIVE) and the irradiation-assisted plastic strains in the vertical direction (lower than two thirds of the RIVE). Because of the large lateral compressive stresses caused by the restraining effect of the unirradiated substrate, the stress field in the mineral grains and at the grain boundary (GB) in the ion-implanted layer is much different than the comparable stress field occurring during neutron irradiation. Hence, the mismatch RIVE causing cracks in the rock-forming minerals leads to different cracking patterns. Ion irradiation continues to be used as a rapid technique to assess the RIVE potential of rocks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2019

Scenario Complexity for Unmanned Aircraft System Traffic

This work introduces an approach to estimate the complexity of a low-altitude air traffic scenario involving multiple UASs using mathematical programming. Given a set of multi-point UAS flight trajectories, vehicle dynamics, and a conflict resolution algorithm, an abstract model is developed such that it can be solved quickly using a mathematical programming optimization software without running high-fidelity simulations that can be computationally expensive and may not suit real-time apA quick and accurate assessment of complexity for a given traffic scenario can help plan and schedule flights to alleviate traffic bottleneck and mitigate operation risks, especially for unmanned aerial system traffic management where high traffic density or complexity is expected. This work introduces a traffic scenario complexity metric that was constructed based on the number of potential conflicts weighted by the conflict resolution cost associated. The cost associated with a conflict is calculated based on the corresponding conflict resolution maneuvers. To obtain the conflict resolution maneuvers, a MILP-based optimization was formulated with the vehicle model and conflict management parameters incorporated. To evaluate the complexity metrics, an approach of using measurements from high-fidelity simulations was proposed. The scenario complexity measurements for 920 random-generated scenarios were obtained through high-fidelity simulations and treated as the ground truth. Two statistics methods: Pearson and Alternative Conditional Expectations were applied for analysis. The results showed that the number of flights has low correlation with the scenario complexity according to the correlation coefficients calculated by both methods. The Alternative Conditional Expectations method shows that the proposed scenario complexity metric has better correlation with the ground truth than the number of potential conflicts.plications. In the abstract model, each vehicle is represented by a time-varied vector associated with position, speed, and heading information. The total extra distance that aircraft need to divert from their original routes to avoid collisions is computed and used to setup a quadratic programming formula. The metrics including the number of conflicts and extra distances travelled by all vehicles are then utilized to estimate the complexity of a given UAS flight scenario. Results and verification against high-fidelity simulations will be provided in the final draft.

traffic complexity↗