Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power flow equation”

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 379 records · Page 21

On the Nature of Navier-stokes Turbulence

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, largest Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The solutions, separate exponentially with time, having a positive Liapunov exponent. Thus the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number is reduced. Several qualitatively different flows are noted. These are, fully chaotic, complex period, weakly chaotic, simple periodic, and fixed-point. Of these, only the fully chaotic flows are classified as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, is both time-dependent and aperiodic. Turbulent solutions are obtained in which energy cascades from large to small-scale motions. In general, the spectral energy transfer takes place between wavenumber bands that are considerably separated. The special transfer can occur either as a result of nonlinear turbulence self-interaction or by interaction of turbulence with mean gradients. Turbulent systems are compared with those studied in kinetic theory. The two types of systems are fundamentally different (continuous and dissipative as opposed to discrete and conservative), but there are similarities. For instance, both are nonlinear and show sensitive dependence on initial conditions. Also, the turbulent and molecular stress tensors are identical if the macroscopic velocities for the turbulent stress are replaced by molecular velocities.

Deissler, Robert G.↗

Structural Weight Estimation for Launch Vehicles

This paper describes some of the work in progress to develop automated structural weight estimation procedures within the Vehicle Analysis Branch (VAB) of the NASA Langley Research Center. One task of the VAB is to perform system studies at the conceptual and early preliminary design stages on launch vehicles and in-space transportation systems. Some examples of these studies for Earth to Orbit (ETO) systems are the Future Space Transportation System [1], Orbit On Demand Vehicle [2], Venture Star [3], and the Personnel Rescue Vehicle[4]. Structural weight calculation for launch vehicle studies can exist on several levels of fidelity. Typically historically based weight equations are used in a vehicle sizing program. Many of the studies in the vehicle analysis branch have been enhanced in terms of structural weight fraction prediction by utilizing some level of off-line structural analysis to incorporate material property, load intensity, and configuration effects which may not be captured by the historical weight equations. Modification of Mass Estimating Relationships (MER's) to assess design and technology impacts on vehicle performance are necessary to prioritize design and technology development decisions. Modern CAD/CAE software, ever increasing computational power and platform independent computer programming languages such as JAVA provide new means to create greater depth of analysis tools which can be included into the conceptual design phase of launch vehicle development. Commercial framework computing environments provide easy to program techniques which coordinate and implement the flow of data in a distributed heterogeneous computing environment. It is the intent of this paper to present a process in development at NASA LaRC for enhanced structural weight estimation using this state of the art computational power.

Cerro, Jeff↗

Ground effect for V/STOL aircraft configurations and its simulation in the wind tunnel. Part 2: Experimental studies

Wind tunnel tests on a finite, knee-blowing-flapped wing and on a direct jet lift configuration were performed in a 30 x 42 inch wind tunnel. The objectives of these tests were, in experimental order:(1) to obtain definitive fixed-versus-moving ground comparisons on a powered, high-lift finite wing and on a lifting jet configuration as a function of model height and lift, (2) to understand the flow mechanism which lead to observed differences between moving and fixed-ground results, (3) to discover whether the above effects of a moving ground plane can be simulated by tangential-blowing boundary layer control at the wind tunnel floor, and (4) to check the form and the numerical constants for the wall jet blowing equation derived theoretically in Part I of this report. This relates wall jet blowing quantities to model lift coefficient.

Hackett, J. E.↗

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

Analytical study of laser supported combustion waves in hydrogen

A one-dimensional energy equation, with constant pressure and area, was used to model the LSC wave. This equation balances convection, conduction, laser energy absorption, radiation energy loss and radiation energy transport. Solutions of this energy equation were obtained to give profiles of temperature and other properties, as well as the relation between laser intensity and mass flux through the wave. The flow through the LSC wave was then conducted through a variable pressure, variable area streamtube to accelerate it to high speed, with the propulsion application in mind. A numerical method for coupling the LSC wave model to the streamtube flow was developed, and a sample calculation was performed. The result shows that 42% of the laser power has been radiated away by the time the gas reaches the throat. It was concluded that in the radially confined flows of interest for propulsion applications, transverse velocities would be less important than in the unconfined flows where air experiments have been conducted.

Kemp, N. H.↗

Reduced dimer production in solar-simulator-pumped continuous wave iodine lasers based on model simulations and scaling and pumping studies

A numerical rate equation model for a continuous wave iodine laser with longitudinally flowing gaseous lasant is validated by approximating two experiments that compare the perfluoroalkyl iodine lasants n-C3F7I and t-C4F9I. The salient feature of the simulations is that the production rate of the dimer (C4F9)2 is reduced by one order of magnitude relative to the dimer (C3F7)2. The model is then used to investigate the kinetic effects of this reduced dimer production, especially how it improves output power. Related parametric and scaling studies are also presented. When dimer production is reduced, more monomer radicals (t-C4F9) are available to combine with iodine ions, thus enhancing depletion of the laser lower level and reducing buildup of the principal quencher, molecular iodine. Fewer iodine molecules result in fewer downward transitions from quenching and more transitions from stimulated emission of lasing photons. Enhanced depletion of the lower level reduces the absorption of lasing photons. The combined result is more lasing photons and proportionally increased output power.

Costen, Robert C.↗

Experimental Study of Eutectic Molten Salts NaCl/KCl/ZnCl2 Heat Transfer Inside a Smooth Tube for High-Temperature Application

Abstract Eutectic salts NaCl-KCl-ZnCl2 and NaCl-KCl-MgCl2 are two of the chloride salt systems that are promising for being used as high-temperature heat transfer fluid (HTF) and thermal energy storage (TES) materials in a wide range of temperatures from 230 °C to 750 °C in concentrated solar thermal power systems. To conduct measurement of the heat transfer coefficient of the molten salt at high temperatures, a circulation system and instrumentation of flow and heat transfer was constructed. Experimental measurement of the convective heat transfer coefficients of NaCl-KCl-ZnCl2 (molar fraction: 13.8–41.9–44.3%) inside tubes has been accomplished to find the most suitable heat transfer correlations. Experience of salt processing and operation of the high-temperature heat transfer test system was obtained. Two correlations, Dittus-Boelter equation and Gnielinski’s correlation for Nusselt number against Reynolds number and Prandtl number, are evaluated using the test results, and the latter correlation is recommended due to its better agreement of prediction against tested data.

Energy & Fuels↗

Active Control of Panel Vibrations Induced by a Boundary Layer Flow

In recent years, active and passive control of sound and vibration in aeroelastic structures have received a great deal of attention due to many potential applications to aerospace and other industries. There exists a great deal of research work done in this area. Recent advances in the control of sound and vibration can be found in the several conference proceedings. In this report we will summarize our research findings supported by the NASA grant NAG-1-1175. The problems of active and passive control of sound and vibration has been investigated by many researchers for a number of years. However, few of the articles are concerned with the sound and vibration with flow-structure interaction. Experimental and numerical studies on the coupling between panel vibration and acoustic radiation due to flow excitation have been done by Maestrello and his associates at NASA/Langley Research Center. Since the coupled system of nonlinear partial differential equations is formidable, an analytical solution to the full problem seems impossible. For this reason, we have to simplify the problem to that of the nonlinear panel vibration induced by a uniform flow or a boundary-layer flow with a given wall pressure distribution. Based on this simplified model, we have been able to study the control and stabilization of the nonlinear panel vibration, which have not been treated satisfactorily by other authors. The vibration suppression will clearly reduce the sound radiation power from the panel. The major research findings will be presented in the next three sections. In Section II we shall describe our results on the boundary control of nonlinear panel vibration, with or without flow excitation. Section III is concerned with active control of the vibration and sound radiation from a nonlinear elastic panel. A detailed description of our work on the parametric vibrational control of nonlinear elastic panel will be presented in Section IV. This paper will be submitted to the Journal of Acoustic Society of America for publication.

Chow, Pao-Liu↗

Investigation of parabolic computational techniques for internal high-speed viscous flows

A feasibility study was conducted to assess the applicability of an existing parabolic analysis (ADD-Axisymmetric Diffuser Duct), developed previously for subsonic viscous internal flows, to mixed supersonic/subsonic flows with heat addition simulating a SCRAMJET combustor. A study was conducted with the ADD code modified to include additional convection effects in the normal momentum equation when supersonic expansion and compression waves were present. It is concluded from the present study that for the class of problems where strong viscous/inviscid interactions are present a global iteration procedure is required.

Anderson, O. L.↗

Motorette Thermal Performance Testing and Modeling for an Electric Motor with Additively Manufactured Hollow Conductors with Integrated Heat Pipes

This paper discusses the design, build and test of a motorette to characterize the thermal performance of an additively manufactured coil integrated with heat pipes. The motorette is powered using a variable frequency AC power supply and cooled through two independent cooling systems. The first cooling system cools the stator core using forced convection of air flowing over a finned heat sink mounted on the stator outer diameter. The second cooling system cools the heat pipes using forced convection of 50-50 water ethylene glycol (WEG) mixture flowing through condenser chambers. A 3D thermal FEA model of the setup is built and heat transfer coefficients (HTC) of convective boundaries are computed using coolant flow rates from test data and empirical equations. Temperature at locations on the coil and heat sink are recorded and compared to thermal model predicted values. A maximum temperature error of 22.7% occurs at 180ARMS and 800 Hz operating point for the AM dual coil. The thermal model provides higher temperatures compared to test results and hence it is on the conservative side but in general, there is good correlation between test and model results. Lastly, opportunities for improvement to both test data measurement and 3D thermal FEA modeling are discussed.

additive manufacturing↗

PyAMG: Algebraic Multigrid Solvers in Python

PyAMG is a Python package of algebraic multigrid (AMG) solvers and supporting tools for approximating the solution to large, sparse linear systems of algebraic equations, Ax = b, where A is an n × n sparse matrix. Sparse linear systems arise in a range of problems in science, from fluid flows to solid mechanics to data analysis. While the direct solvers available in SciPy’s sparse linear algebra package (scipy.sparse.linalg) are highly efficient, in many cases iterative methods are preferred due to overall complexity. However, the iterative methods in SciPy, such as CG and GMRES, often require an efficient preconditioner in order to achieve a lower complexity. Preconditioning is a powerful tool whereby the conditioning of the linear system and convergence rate of the iterative method are both dramatically improved. PyAMG constructs multigrid solvers for use as a preconditioner in this setting. A summary of multigrid and algebraic multigrid solvers can be found in Olson (2015a), in Olson (2015b), and in Falgout (2006); a detailed description can be found in Briggs et al. (2000) and Trottenberg et al. (2001).

97 MATHEMATICS AND COMPUTING↗

A CFD-based aerodynamic design procedure for hypersonic wind-tunnel nozzles

A new procedure which unifies the best of current classical design practices, computational fluid dynamics (CFD), and optimization procedures is demonstrated for designing the aerodynamic lines of hypersonic wind-tunnel nozzles. The new procedure can be used to design hypersonic wind tunnel nozzles with thick boundary layers where the classical design procedure has been shown to break down. An efficient CFD code, which solves the parabolized Navier-Stokes (PNS) equations using an explicit upwind algorithm, is coupled to a least-squares (LS) optimization procedure. A LS problem is formulated to minimize the difference between the computed flow field and the objective function, consisting of the centerline Mach number distribution and the exit Mach number and flow angle profiles. The aerodynamic lines of the nozzle are defined using a cubic spline, the slopes of which are optimized with the design procedure. The advantages of the new procedure are that it allows full use of powerful CFD codes in the design process, solves an optimization problem to determine the new contour, can be used to design new nozzles or improve sections of existing nozzles, and automatically compensates the nozzle contour for viscous effects as part of the unified design procedure. The new procedure is demonstrated by designing two Mach 15, a Mach 12, and a Mach 18 helium nozzles. The flexibility of the procedure is demonstrated by designing the two Mach 15 nozzles using different constraints, the first nozzle for a fixed length and exit diameter and the second nozzle for a fixed length and throat diameter. The computed flow field for the Mach 15 least squares parabolized Navier-Stokes (LS/PNS) designed nozzle is compared with the classically designed nozzle and demonstrates a significant improvement in the flow expansion process and uniform core region.

Korte, John J.↗

Transition prediction and control in subsonic flow over a hump

The influence of a surface roughness element in the form of a two-dimensional hump on the transition location in a two-dimensional subsonic flow with a free-stream Mach number up to 0.8 is evaluated. Linear stability theory, coupled with the N-factor transition criterion, is used in the evaluation. The mean flow over the hump is calculated by solving the interacting boundary-layer equations; the viscous-inviscid coupling is taken into consideration, and the flow is solved within the separation bubble. The effects of hump height, length, location, and shape; unit Reynolds number; free-stream Mach number, continuous suction level; location of a suction strip; continuous cooling level; and location of a heating strip on the transition location are evaluated. The N-factor criterion predictions agree well with the experimental correlation of Fage; in addition, the N-factor criterion is more general and powerful than experimental correlations. The theoretically predicted effects of the hump's parameters and flow conditions on transition location are consistent and in agreement with both wind-tunnel and flight observations.

Masad, Jamal A.↗

Experimental and Numerical Investigation of Buoyancy Driven Convection During PDAMNA Thin Film Growth

This paper presents results from numerical simulations as well as laboratory experiments of buoyancy driven convection in an ampoule under varying heating and gravitational acceleration loadings. The modeling effort in this work resolves the large scale natural convective motion that occurs in the fluid during photodeposition of polydiacetelene films which is due to energy absorbed by the growth solution from a UV source. Consequently, the growth kinetics of the film are ignored in the model discussed here, and also a much simplified ampoule geometry is considered. The objective of this work is to validate the numerical prediction on the strength and structure of buoyancy driven convection that could occur under terrestrial conditions during nonlinear optical film growth. The validation is used to enable a reliable predictive capability on the nature and strength of the convective motion under low gravity conditions. The ampoule geometry is in the form of a parallelepiped with rectangular faces. The numerical results obtained from the solution to the Boussinesq equations show that natural convection will occur regardless of the orientation of the UV source with respect to the gravity vector. The least strong convective motion occurred with the UV beam directed at the top face of the parallelepiped. The strength of the convective motion was found to be almost linearly proportional to the total power of the UV source. Also, it was found that the strength of the convective motion decreased linearly with the gravity due to acceleration. The pattern of the convective flow on the other hand, depended on the source location.

Antar, Basil N.↗

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Dugan, E. T.↗

Thermal Spectrum Molten Salt-Fueled Reactor Reference Plant Model

This report details a fully coupled neutronics thermal hydraulics reference plant model for a thermal spectrum molten salt-fueled reactor. The Multiphysics model is developed on the Nuclear Regulatory Commission's Comprehensive Reactor Analysis Bundle available on the Idaho National Laboratory high-performance computer, which natively and seamlessly couples Griffin and Sam codes of the Multiphysics Object-Oriented Simulation Environment (MOOSE) based applications. Griffin provides the reactor physics capabilities, k-eigenvalue, delayed neutron precursor concentrations, and transient solutions for stationary and flowing fuels. Sam provides the solution of mass, momentum, and energy conservation equations of the whole system with a set of system components for one-dimensional single-phase flow. The neutronics feedback model relies primarily on fuel and moderator temperatures, fuel salt density, and delayed neutron precursors distributions. This report presents our results for the coupled steady state core and an unprotected loss of flow event. Although this model is prototypical regarding capabilities in the Comprehensive Reactor Analysis Bundle, its results are consistent with experimental data of the pump transient tests that were performed at zero power of the Molten Salt Reactor Experiment (MSRE).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiphysics Demonstration of Temperature-Driven Assembly Bowing in SFRs using MOOSE-Based Codes

Core bowing is an important passive safety mechanism in liquid metal cooled fast reactors. When the core restraint system is properly designed, temperature and flux gradients influence assemblies in the core to bow into less reactive configurations during accident scenarios, resulting in negative reactivity feedback. Prediction of core bowing involves complex interplay of radiation transport, impacts of fluid flow and heat transfer on duct temperature, and mechanical responses to the induced temperature and flux gradients. Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program [1], an integrated multiphysics approach is being developed to model the core bowing phenomena in liquid metal-cooled fast reactors with the Multiphysics Object Oriented Simulation Environment (MOOSE) [2]. In this methodology, the MOOSE-based reactor physics code Griffin [3] will solve the neutron transport equation and determine the power distribution. With the detailed power distribution from Griffin, the subchannel analysis codes MOOSE-Subchannel [4] and Pronghorn [5] are utilized to calculate the assembly temperature distribution. MOOSE’s Solid Mechanics [6] and Contact [7] Modules are leveraged to calculate the thermal expansion and duct bowing displacement with the duct wall temperature from thermal hydraulics calculation. In this work, an initial one-way coupling demonstration of the integrated multiphysics approach has been performed on a seven-assembly problem based on the sodium-cooled fast reactor ABR-1000 design [8]. The neutronics calculation with Griffin is not yet involved in the current simulation. MOOSE-Subchannel and Pronghorn evaluate fluid and solid temperature based on a fixed power distribution. In addition, one-way coupling is utilized in this coupled calculation, via Pronghorn passing the duct temperature data to the MOOSE Solid Mechanics calculation. An assessment of the Solid Mechanics module was performed in parallel to verify duct bowing behavior with duct-to-duct contact phenomenon [9]. The displacement from MOOSE Solid Mechanics is not yet transferred back and utilized in the Pronghorn and MOOSE-Subchannel calculation. This model will be available on the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) repository [10]. Future stages of this work will involve solving problems of increasing complexity as well as adding more physics (e.g. reactor physics) to the integrated workflow to reach the end goal of modeling the core bowing phenomenon with an integrated multiphysics workflow.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗