Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code development and validation”

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 343 records · Page 19

Connect microstructure evolution to property degradation with validated simulation

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗

Proposed New STELLA-2 Tests

In 2025, KAERI and Argonne initiated a collaboration based on software validation using data from tests performed at KAERI’s STELLA-2 facility. STELLA-2 is a large-scale sodium thermal-hydraulic test facility that was originally developed to support KAERI’s development of the PGSFR design concept. KAERI has performed a large number of tests at STELLA-2, providing valuable data to support sodium fast reactor code validation. KAERI will be conducting new tests in 2026, targeting system conditions that were not achieved during previous test campaigns. KAERI has requested from Argonne proposed tests that could be performed during the upcoming test campaign. This report documents Argonne’s proposed tests in fulfillment of KAERI’s request.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Investigation of aerothermodynamics and optical radiation in the AFE hypersonic flow field

Research for the radiation experiments on the Aeroassist Flight Experiment (AFE) was performed to obtain a data base for development of engineering requirements for aerobrakes. Due to funding restrictions, the necessary Aerobrake design data were obtained from ground based experiments, specifically arc-jet wind tunnels. Except for the instrument windows, final development of the AFE radiometers was completed. Window definition included several designs to be flight validated in arc-jets. This work was completed, and successful designs have evolved which will yield full scientific return from the flight experiment. The theoretical work includes final code development to describe the spacecraft environment to support instrument definition, optical radiation codes that operate at significantly improved speeds, and calculation of radiation forebody loads and effects on experiment objectives for varying AFE weights and trajectories. Furthermore, radiant flux vectors at the AFE base were predicted to be used for afterbody instrument definition. The Vacuum Ultraviolet Experiment for the Arc-Jet was completed, and arc-jet experiment conditions and experiment objectives were refined to incorporate the most current thinking. Ballistic range experiments were completed. Several shadowgraphs indicating corner turning angles and waviness structures were obtained, and the results were used to estimate the compression region location for the AFE flight. The waviness measured has helped to model the temporal variations of wake radiation for purposes of defining radiometers to measure the energy content of wake dynamics.

Palumbo, Giusseppe↗

Pressurization of cryogens: A review of current technology and its applicability to low-gravity conditions

A review of technology, history, and current status for pressurized expulsion of cryogenic tankage is presented. Use of tank pressurization to expel cryogenic fluids will continue to be studied for future spacecraft applications over a range of operating conditions in the low-gravity environment. The review examines experimental test results and analytical model development for quiescent and agitated conditions in normal-gravity, followed by a discussion of pressurization and expulsion in low-gravity. Validated, 1-D, finite difference codes exist for the prediction of pressurant mass requirements within the range of quiescent normal-gravity test data. To date, the effects of liquid sloshing have been characterized by tests in normal-gravity, but analytical models capable of predicting pressurant gas requirements remain unavailable. Efforts to develop multidimensional modeling capabilities in both normal and low-gravity have recently occurred. Low-gravity cryogenic fluid transfer experiments are needed to obtain low-gravity pressurized expulsion data. This data is required to guide analytical model development and to verify code performance.

Vandresar, N. T.↗

Pressurization of cryogens - A review of current technology and its applicability to low-gravity conditions

A review of technology, history, and current status for pressurized expulsion of cryogenic tankage is presented. Use of tank pressurization to expel cryogenic fluid will continue to be studied for future spacecraft applications over a range of operating conditions in the low-gravity environment. The review examines experimental test results and analytical model development for quiescent and agitated conditions in normal-gravity followed by a discussion of pressurization and expulsion in low-gravity. Validated, 1-D, finite difference codes exist for the prediction of pressurant mass requirements within the range of quiescent normal-gravity test data. To date, the effects of liquid sloshing have been characterized by tests in normal-gravity, but analytical models capable of predicting pressurant gas requirements remain unavailable. Efforts to develop multidimensional modeling capabilities in both normal and low-gravity have recently occurred. Low-gravity cryogenic fluid transfer experiments are needed to obtain low-gravity pressurized expulsion data. This data is required to guide analytical model development and to verify code performance.

Van Dresar, N. T.↗

Improved Spacecraft Materials for Radiation Protection

Methods by which radiation shielding is optimized need to be developed and materials of improved shielding characteristics identified and validated. The galactic cosmic rays (GCR) are very penetrating and the energy absorbed by the astronaut behind the shield is nearly independent of shield composition and even the shield thickness. However, the mix of particles in the transmitted beam changes rapidly with shield material composition and thickness. This results in part from the breakup of the high-energy heavy ions of the GCR which make contributions to biological effects out of proportion to their deposited energy. So the mixture of particles in the radiation field changes with shielding and the control of risk contributions from dominant particle types is critical to reducing the hazard to the astronaut. The risk of biological injury for a given particle type depends on the type of biological effect and is specific to cell or tissue type. Thus, one is faced with choosing materials which may protect a given tissue against a given effect but leave unchanged or even increase the risk of other effects in the same tissue or increase the risks to other adjacent tissues of a different type in the same individual. The optimization of shield composition will then be tied to a specific tissue and risk to that tissue. Such peculiarities arise from the complicated mixture of particles, the nature of their biological response, and the details of their interaction with material constituents. Aside from the understanding of the biological response to specific components, one also needs an accurate understanding of the radiation emerging from the shield material. This latter subject has been a principal element of this project. In the past ten years our understanding of space radiation interactions with materials has changed radically, with a large impact on shield design. For example, the NCRP estimated that only 2 g/sq cm. of aluminum would be required to meet the annual 500 mSv limit for the exposure of the blood forming organs (this limit is strictly for LEO but can be used as a guideline for the Mars mission analysis). The current estimates require aluminum shield thicknesses above 50 g/sq cm., which is impractical. In such a heavily shielded vehicle, the neutrons produced throughout the vehicle also contribute significantly to the exposure and this demands greater care in describing the angular dependence of secondary particle production processes. As such the continued testing of databases and transport procedures in laboratory and spaceflight experiments has continued. This has been the focus of much of the last year's activity and has resulted in improved neutron prediction capability. These new methods have also improved our understanding of the surface environment of Mars. The Mars 2003 NRA HEDS related surface science requirements were driven by the need to validate predictions on the upward flux of neutrons produced in the Martian regolith and bedrock made by the codes developed under this project. The codes used in the surface environment definition are also being used to look at in situ resources for the development of construction material for Martian surface facilities. For example, synthesis of polyimides and polyethylene as binders of regolith for developing basic structural elements has been studied and targets built for accelerator beam testing of radiation shielding properties. Preliminary mechanical tests have also been promising. Improved spacecraft materials have been identified (using the criteria reported by this project at the last conference) as potentially important for future shielding materials. These are liquid hydrogen, hydrogenated nanofibers, liquid methane, LiH, Polyethylene, Polysulfone, and Polyetherimide (in order of decreasing shield performance). Some of the materials are multifunctional and are required for other onboard systems. We are currently preparing software for trade studies with these materials relative to the Mars Reference Mission as required in the project's final year.

Wilson, John W.↗

New free-wake analysis of rotorcraft hover performance using influence coefficients

Free-wake analyses of helicopter rotor wakes in hover using time stepping have been shown to encounter instabilities which preclude convergence to valid free-vortex solutions for rotor-wake geometries. Previous work has demonstrated that these convergence difficulties can be overcome by implementing a new free-wake analysis method based on the use of influence coefficients. The present paper reviews this approach and documents its incorporation into a hover performance analysis called Evaluation of Hover Performance using Influence Coefficients (EHPIC). The technical principles underlying the EHPIC code are described with emphasis on steps taken to develop the single-filament wake models used in previous work into a multifilament wake valid for realistic hover performance predictions. The coupling of the wake model to a lifting surface loads analysis is described, and sample problems are solved that illustrate the robustness of the method. Performance calculations are also undertaken for hover to illustrate the utility of EHPIC in the analysis of rotorcraft performance.

Quackenbush, Todd R.↗

Tone and Broadband Noise Separation from Acoustic Data of a Scale-Model Counter-Rotating Open Rotor

Renewed interest in contra-rotating open rotor technology for aircraft propulsion application has prompted the development of advanced diagnostic tools for better design and improved acoustical performance. In particular, the determination of tonal and broadband components of open rotor acoustic spectra is essential for properly assessing the noise control parameters and also for validating the open rotor noise simulation codes. The technique of phase averaging has been employed to separate the tone and broadband components from a single rotor, but this method does not work for the two-shaft contra-rotating open rotor. A new signal processing technique was recently developed to process the contra-rotating open rotor acoustic data. The technique was first tested using acoustic data taken of a hobby aircraft open rotor propeller, and reported previously. The intent of the present work is to verify and validate the applicability of the new technique to a realistic one-fifth scale open rotor model which has 12 forward and 10 aft contra-rotating blades operating at realistic forward flight Mach numbers and tip speeds. The results and discussions of that study are presented in this paper.

Acoustic Data↗

Tone and Broadband Noise Separation from Acoustic Data of a Scale-Model Contra-Rotating Open Rotor

Renewed interest in contra-rotating open rotor technology for aircraft propulsion application has prompted the development of advanced diagnostic tools for better design and improved acoustical performance. In particular, the determination of tonal and broadband components of open rotor acoustic spectra is essential for properly assessing the noise control parameters and also for validating the open rotor noise simulation codes. The technique of phase averaging has been employed to separate the tone and broadband components from a single rotor, but this method does not work for the two-shaft contra-rotating open rotor. A new signal processing technique was recently developed to process the contra-rotating open rotor acoustic data. The technique was first tested using acoustic data taken of a hobby aircraft open rotor propeller, and reported previously. The intent of the present work is to verify and validate the applicability of the new technique to a realistic one-fifth scale open rotor model which has 12 forward and 10 aft contra-rotating blades operating at realistic forward flight Mach numbers and tip speeds. The results and discussions of that study are presented in this paper.

Noise Separation↗

Development and Implementation of a New AI-Based Tool to Support Fast Reactor Software Model Generation and Validation

This report summarizes FY26 work to develop Maggie, an artificial intelligence-based assistant designed to support software model generation and validation activities for fast reactor analysis codes. The project established a modular, code-agnostic software architecture that separates reusable agent capabilities from code-specific knowledge and tools, with initial implementation focused on the FRP-supported fast reactor safety analysis code SAS4A/SASSYS1 (SAS). A curated SAS-specific knowledge base was assembled from the code manual, training materials, historical analysis reports, and representative input files, and was integrated through retrieval-augmented generation to ground Maggie’s responses in authoritative sources. Maggie was deployed on the internal Argonne network, where it demonstrated practical user-facing capability as a chatbot for answering natural language questions about SAS and retrieving relevant technical information. Demonstration cases also showed that Maggie can generate useful snippets of SAS input for selected modeling tasks, while highlighting current limitations in reliability and consistency for more complex input generation tasks. Overall, the FY26 effort established the technical foundation for an AI-assisted capability intended to improve the efficiency, consistency, and accessibility of fast reactor software model development at Argonne and, with further improvements, to support eventual use by the broader fast reactor community, including industry users of FRP-supported analysis tools.

Thomas, Rachel [Argonne National Laboratory (ANL),↗

Technical Assessment of the Application of Digital Twin and Prognostic Tools for Condition Monitoring

This report was prepared for the U.S. Nuclear Regulatory Commission (NRC) to present use cases of the application of advanced technologies toward meeting the current and future regulatory requirements for maintenance and condition monitoring of structures, systems, and components (SSCs). The advanced technologies considered in this work, collectively referred to as digital twin (DT) technologies, are advanced sensors and instrumentation, data analytics, machine learning and artificial intelligence (ML/AI), and physics-based models. The report presents two use cases of reactor coolant pumps (RCPs) and heat pipes in nuclear power plants (NPPs) with technical and regulatory considerations and opportunities in using advanced technologies for conditional monitoring. Key findings from the exploration of these considerations are as follows: - Uncertainties in sensor data and model predictions must be rigorously addressed through validation and verification processes - Regulatory compliance is paramount, necessitating data driven models to be developed in line with existing codes and standards, as well as considering potential future guidelines for advanced reactors - Explainability and transparency in ML/AI models are essential for developing operator trust and regulatory review, including methods that enhance the interpretability of complex data-driven predictions - Condition monitoring programs must be evaluated for their effectiveness in reducing maintenance-preventable function failures (MPFF) and aligning with plant performance criteria - The deployment of advanced technologies for condition monitoring could lead to a transition from periodic to continuous monitoring, thereby optimizing maintenance schedules - Collaborative efforts between industry stakeholders, regulatory bodies, and technology developers are crucial for the successful adoption of advanced technologies for condition monitoring systems in nuclear facilities In summary, the introduction of advanced technologies into condition monitoring programs represents a significant leap forward in the domain of NPP maintenance. By harnessing the capabilities of advanced sensors, data analytics, and ML/AI, NPP operators can transition from a time-based to a condition-based maintenance approach. This shift can potentially enhance the reliability and safety of critical plant components while optimizing maintenance efforts and minimizing unnecessary outages. The NRC is continuing to explore the regulatory aspects of advanced technologies as part of inservice inspection and inservice testing (ISI and IST) programs by pursuing additional research in this technical area.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Theoretical investigation of gas-surface interactions

The goal of this project was to develop computational tools for the calculation of the electronic structure of molecules containing heavy atoms, and to use these tools in the study of catalytic processes, with the overall objective of gaining an understanding of the catalytic process which could be used to design more efficient catalysts. The main catalytic system of interest was the combustion of hydrogen on platinum surfaces. Under this project, a flexible Dirac-Hartree-Fock (DHF) program has been developed, a code to calculate DHF correlation energies at the second-order Moeller-Plesset perturbation (PP2) level is almost complete, and code to include correlation at the MCSCF and MCSCF/MP2 level is planned. The tools so far developed have been validated and used to calibrate some more approximate methods, and applied to investigate the importance of relativistic effects in the bonding of hydrogen to platinum.

Dyall, Kenneth G.↗

Description, Usage, and Validation of the MVL-15 Modified Vortex Lattice Analysis Capability

MVL-15 is the most recent version of the Modified Vortex-Lattice (MVL) code developed within the Aerodynamics Systems Analysis Branch (ASAB) at NASA LaRC. The term "modified" refers to the primary modification of the core vortex-lattice methodology: inclusion of viscous aerodynamics tables that are linked to the linear solution via iterative processes. The inclusion of the viscous aerodynamics inherently converts the MVL-15 from a purely analytic linearized method to a semi-empirical blend which retains the rapid execution speed of the linearized method while empirically characterizing the section aerodynamics at all spanwise lattice points. The modification provides a means to assess non-linear effects on lift that occur at angles of attack near stall, and provides a means to determine the drag associated with the application of design strategies for lift augmentation such as the use of flaps or blowing. The MVL-15 code is applicable to the analyses of aircraft aerodynamics during cruise, but it is most advantageously applied to the analysis of aircraft operating in various high-lift configurations. The MVL methodology has been previously conceived and implemented; the initial concept version was delivered to the ASAB in 2001 (van Dam, C.), subsequently revised (Gelhausen, P. and Ozoroski, T. 2002 / AVID Inc., Gelhausen, P., and Roberts, M. 2004), and then overhauled (Ozoroski, T., Hahn, A. 2008). The latest version, MVL-15 has been refined to provide analysis transparency and enhanced to meet the analysis requirements of the Environmentally Responsible Aviation (ERA) Project. Each revision has been implemented with reasonable success. Separate applications of the methodology are in use, including a similar in-house capability, developed by Olson, E. that is tailored for structural and acoustics analyses. A central premise of the methodology is that viscous aerodynamic data can be associated with analytic inviscid aerodynamic results at each spanwise wing section, thereby providing a pathway to map viscous data to the inviscid results. However, a number of factors can sidetrack the analysis consistency during various stages of this process. For example, it should be expected that the final airplane lift curve and drag polar results depend strongly on the geometry and aerodynamics of the airfoil section; however, flap deflections and flap chord extensions change the local reference geometry of the input airfoil, the airplane wing, the tabulated non-dimensional viscous aerodynamics, and the spanwise links between the linear and the viscous aerodynamics. These changes also affect the bound circulation and therefore, calculation and integration of the induced angle of attack and induced drag. MVL-15 is configured to ensure these types of challenges are properly addressed. This report is a comprehensive manual describing the theory, use, and validation of the MVL-15 analysis tool. Section 3 summarizes theoretical, procedural, and characteristic features of MVL-15, and includes a list of the files required to setup, execute, and summarize an analysis. Section 4, Section 5, Section 6, and Section 7 combine to comprise the User's Guide portions of this report. The MVL-15 input and output files are described in Section 4 and Section 5, respectively; the descriptions are supplemented with example files and information about the file formats, parameter definitions, and typical parameter values. Section 6 describes the Wing Geometry Setup Utility and the 2d-Variants Utility files that simplify and assist setting up a consistent set of MVL-15 geometry and aerodynamics input parameters and input files. Section 7 describes the use of the 3d-Results Presentation Utility file that can be used to automatically create summary tables and charts from the MVL-15 output files. Section 8 documents the Validation Results of an extensive and varied validation test matrix, including results of an airplane analysis representative of the ERA Program. A start-to-finish example of the airplane analysis procedure is described in Section 7.

Ozoroski, Thomas A.↗

A wall interference assessment/correction system

A Wall Signature method, the Hackett method, has been selected to be adapted for the 12-ft Wind Tunnel wall interference assessment/correction (WIAC) system in the present phase. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating wall interferences at the model location. The Wall Signature method will be formulated for application to the unique geometry of the 12-ft Tunnel. The development and implementation of a working prototype will be completed, delivered and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method-Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free-air simulation. Each set of wind tunnel/test article numerical simulations provides data to validate the WIAC method. A numerical wind tunnel test simulation is initiated to validate the WIAC methods developed in the project. In the present reported period, the blockage correction has been developed and implemented for a rectangular tunnel as well as the 12-ft Pressure Tunnel. An improved wall interference assessment and correction method for three-dimensional wind tunnel testing is presented in the appendix.

Lo, Ching F.↗

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗

Development of a generic combustion stability code for liquid propellant rocket engines

The mathematical framework for a combustion stability analysis code is outlined. The goal for the code is to be general enough in problem treatment so that its validity and accuracy extend over a wide range of problem applications and that it lends the convenience for any future model improvement if necessary. An approach for modeling the combustion dynamics is devised to meet both requirements. An open-loop numerical procedure is also developed to mechanistically model various combustion processes for determining the stability parameters.

Fang, J. J.↗

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Proteus Navier-Stokes code

An effort is currently underway at NASA Lewis to develop two and three dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. Proteus solves the Reynolds-averaged, unsteady, compressible Navier-Stokes equations in strong conservation law form. Turbulence is modeled using a Baldwin-Lomax based algebraic eddy viscosity model. In addition, options are available to solve thin layer or Euler equations, and to eliminate the energy equation by assuming constant stagnation enthalpy. An extensive series of validation cases have been run, primarily using the two dimensional planar/axisymmetric version of the code. Several flows were computed that have exact solution such as: fully developed channel and pipe flow; Couette flow with and without pressure gradients; unsteady Couette flow formation; flow near a suddenly accelerated flat plate; flow between concentric rotating cylinders; and flow near a rotating disk. The two dimensional version of the Proteus code has been released, and the three dimensional code is scheduled for release in late 1991.

Towne, Charles E.↗