Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOUNDARY LAYER TRANSITION - EFFECT OF PRESSURE GRADIENT”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

A magnetohydrodynamical model of the inner coma of Comet Halley

The magnetometer onboard the Giotto spacecraft observed a diamagnetic cavity surrounding the nucleus of Comet Halley. The location of the boundary of this diamagnetic cavity is determined by a balance between an inward magnetic pressure gradient force and an outward ion-neutral drag force, associated with collisions between the outwardly flowing neutrals and the stagnated ions. A one-dimensional time-dependent MHD model has been developed for the inner coma of Comet Halley, and includes ion-neutral collisions, photochemical production and loss of plasma, and finite conductivity. This model is used to investigate the plasma dynamics in the vicinity of the diamagnetic cavity boundary surface. A narrow transition layer with enhanced plasma density is shown to exist just inside the boundary, although a full understanding of this layer will require a two- or three-dimensional MHD model. The flux of cometary ions into this shocklike layer is removed by electron-ion recombination. The thickness of this layer is determined by the Mach number of the incident flow.

Cravens, T. E.↗
DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1984

Convective heat transfer studies at high temperatures with pressure gradient for inlet flow Mach number of 0.45

Heat transfer measurements were determined for a flat plate with and without pressure gradient for various free stream temperatures, wall temperature ratios, and Reynolds numbers for an inlet flow Mach number of 0.45, which is a representative inlet Mach number for gas turbine rotor blades. A shock tube generated the high temperature and pressure air flow, and a variable geometry test section was used to produce inlet flow Mach number of 0.45 and accelerate the flow over the plate to sonic velocity. Thin-film platinum heat gages recorded the local heat flux for laminar, transition, and turbulent boundary layers. The free stream temperatures varied from 611 R (339 K) to 3840 R (2133 K) for a T(w)/T(r,g) temperature ratio of 0.87 to 0.14. The Reynolds number over the heat gages varied from 3000 to 690,000. The experimental heat transfer data were correlated with laminar and turbulent boundary layer theories for the range of temperatures and Reynolds numbers and the transition phenomenon was examined.

Pedrosa, A. C. F.↗
NASA NTRS1998

Transition in Turbine Flows

We have further developed our capabilities to analyze transition in turbine boundary layers from first principles by integrating the nonlinear parabolized stability equations (PSE) with improved initial and boundary conditions. With modified iteration schemes, we are able to proceed deeper into the transition region where skin friction coefficient and heat transfer coefficient significantly increase. Initial and boundary conditions at elevated turbulence levels can be derived by receptivity analysis. Test runs for ERCOFTAC test case T3A at 2.4\% turbulence level provide results in good agreement with the experimental data. The sharper minimum of the skin coefficient also shown by DNS results is likely due to the missing intermittency. The method has been applied to various experimentally studied turbine blades (UTRC, VKI, Zierke, Langston, Hippensteele, and others). The PSE results, though physically reasonable, do not agree as well as expected with the experimental findings. We have, therefore, performed an extensive search for the reasons of the seemingly systematic deviations. A first source of uncertainty has been found in the often insufficient documentation of the experiments (e.g. on blockage by end-wall boundary layers). However, variation of the relevant parameters does not lead to more satisfactory agreement. A second reason has been found in the "standard procedure" which considers a 2D flow at midspan and uses a panel code and subsequent boundary-layer code to obtain the laminar basic flow for the transition analysis. Comparison with the pressure distribution obtained with a 3D design code (RVC3D) shows significant three-dimensionality of the flow (e.g. in the UTRC experiments). The spanwise variation has been neglected in our original PSE code. To overcome this problem, we have developed the PSE/3D for fully 3D boundary layers to account for streamwise and spanwise variations. Since the design code does not provide the boundary-layer flow with sufficient resolution, we have generated the Euler solution and employed a 3D boundary-layer code to obtain the viscous basic flow. Although only the linear stability level of PSE/3D has been implemented so far, the discrepancies with the experiments change but do not disappear. We still find deviations between the computed and experimental variations of C(sub f), and St along the blade for laminar flow. The main reason can be seen by comparing the solution of the boundary-layer code with the viscous results of the design code. The conventional boundary-layer solution exhibits an asymptotic behavior appropriate in external aerodynamics but does not match the steep gradients of the inviscid flow through the passage and consequently provides biased results for C(sub f), and St. An attempt is currently being made to correct this deficiency. Before attempting to perform the transition analysis for the viscous flow provided by the design code, we have analyzed the implementation and "best possible" results. Code and results exhibit flaws that may negatively affect the design and are intolerable for transition analysis. Therefore, we have decided to develop a new code to obtain a reliable basis for stability and transition studies. We expect to report improved results by the time of the meeting.

Herbert, Thorwald↗
NASA NTRS2003

Turbulence in Supercritical O2/H2 and C7H16/N2 Mixing Layers

This report presents a study of numerical simulations of mixing layers developing between opposing flows of paired fluids under supercritical conditions, the purpose of the study being to elucidate chemical-species-specific aspects of turbulence. The simulations were performed for two different fluid pairs O2/H2 and C7H16/N2 at similar reduced initial pressures (reduced pressure is defined as pressure divided by critical pressure). Thermodynamically, O2/H2 behaves more nearly like an ideal mixture and has greater solubility, relative to C7H16/N2, which departs strongly from ideality. Because of a specified smaller initial density stratification, the C7H16/N2 layers exhibited greater levels of growth, global molecular mixing, and turbulence. However, smaller density gradients at the transitional state for the O2/H2 system were interpreted as indicating that locally, this system exhibits enhanced mixing as a consequence of its greater solubility and closer approach to ideality. These thermodynamic features were shown to affect entropy dissipation, which was found to be larger for O2/H2 and concentrated in high-density-gradient-magnitude regions that are distortions of the initial density-stratification boundary. In C7H16/N2, the regions of largest dissipation were found to lie in high-density-gradient-magnitude regions that result from mixing of the two fluids.

Bellan, Josette↗
NASA NTRS1980

Variable-grid-size transformation for solving nonsimilar laminar and turbulent boundary-layers

A numerical code has been developed to calculate nonsimilar boundary-layers (either laminar, transitional, or turbulent). The formulation of the equations of motion employs two transformations for the y-coordinate which provide a dense spacing of points near the wall and, therefore, can provide engineering solutions using relatively few points in the y-direction. The turbulent viscosity is correlated in terms of a length scale, which is in turn a function of the physical y-coordinate and of the mean-flow velocity gradient at that point. The boundary-layer solutions computed using the present code are compared with theoretical solutions presented by other investigators and with data from the literature. Comparisons with the data indicate that the theoretical solutions provide suitable values for the pitot-pressure profiles and the displacement thickness. Significant differences exist in the one comparison between the theoretical and the experimental heat-transfer rates for hypersonic flow past a flat plate.

Bertin, J. J.↗
NASA NTRSDate not supplied

Investigation of a Smooth Local Correlation-based Transition Model in a Discrete-Adjoint Aerodynamic Shape Optimization Algorithm

A smooth local correlation-based transition model is fully coupled to a RANS-based Newton-Krylov flow solver and discrete-adjoint gradient-based optimization algorithm. The free-transition optimization framework is evaluated using lift-constrained drag minimizations of airfoils at design conditions ranging from light to single-aisle aircraft and an infinite swept wing at design conditions representative of a transonic strut-braced wing aircraft. The impact of the streamwise grid resolution on the ability of the optimization algorithm to delay boundary-layer transition is investigated, with the results demonstrating that streamwise grid resolution requirements increase as the transition length decreases with increasing Reynolds number. The optimization problem at the light aircraft design conditions is demonstrated to be multi-modal, with the optimization algorithm producing two distinct designs: one with a thin, reflexed trailing edge and steep pressure recovery regions, the other with increased aft loading, with the latter design outperforming the former. A drag minimization of an airfoil at transonic design conditions demonstrates that the optimization algorithm successfully trades a decrease in viscous drag by delaying boundary-layer transition with an increase in wave drag, while the drag minimization of an infinite swept wing demonstrates the capability of the optimizational gorithm to delay both Tollmien-Schlichting and stationary crossflow instabilities.

AATT↗
NASA NTRS2002

Aspects of Supercritical Turbulence: Direct Numerical Simulation of O2/H2 and C7H16/N2 Temporal Mixing Layers

Results from Direct Numerical Simulations of temporal, supercritical mixing layers for two species systems are analyzed to elucidate species-specific turbulence aspects. The two species systems, O2/H2 and C7HG16/N2, have different thermodynamic characteristics; thus, although the simulations are performed at similar reduced pressure (ratio of the pressure to the critical pressure), the former system is dose to mixture ideality and has a relatively high solubility with respect to the latter, which exhibits strong departures from mixture ideality Due to the specified, smaller initial density stratification, the C7H16/N2 layers display higher growth and increased global molecular mixing as well as larger turbulence levels. However, smaller density gradients at the transitional state for the O2/H2 system indicate that on a local basis, the layer exhibits an enhanced mixing, this being attributed to the increased solubility and to mixture ideality. These thermodynamic features are shown to affect the irreversible entropy production (i.e. the dissipation), which is larger for the O2/H2 layer and is primarily concentrated in high density-gradient magnitude regions that are distortions of the initial density stratification boundary. In contrast, the regions of largest dissipation in the C7H16/N2 layer are located in high density-gradient magnitude regions resulting from the mixing of the two fluids.

Bellan, J.↗
NASA NTRS1985

Measurements in the flow field of a cylinder with a laser transit anemometer and a drag rake in the Langley 0.3 m transonic cryogenic tunnel

An experiment was conducted at the 0.3 m Transonic Cryogenic Tunnel using a Laser Transit Anemometer (LTA) to probe the flow field around a 3.05 centimeter-diameter circular cylinder. Measurements were made along the stagnation line and determination of particle size was evaluated by their ability to follow the flow field. The LTA system was also used to scan into the boundary layer near the 45 degree point on the model. Results of these scans are presented in graphic and tabular form. Flow field particle seeding was accomplished by inbleeding dry service air. The residual moisture (5-10 ppm) condensed and formed ice particles which served as Mie scattering centers for the LTA detection system. Comparison of data taken along the stagnation streamline with theory indicated that these particles tracked the velocity gradient of the flow. Tunnel operating conditions for the tests were a Mach number of 0.3, a pressure of 1.93 x 100000 n/m squared, and a temperature of 225 degrees K. Free stream Mach number and pressure were varied for the particle size determination.

Honaker, W. C.↗
NASA NTRS1991

Spatial Direct Numerical Simulation of Boundary-Layer Transition Mechanisms: Validation of PSE Theory

A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to critically evaluate the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave break-down. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the 2-D test problem. This discrepancy is attributed to far-field boundary- condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.

Joslin, R. D.↗
NASA NTRS1988

On the generation of mean flows by the interaction of Goertler vortices and Tollmien-Schlichting waves in curved channel flows

There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete in the nonlinear regime. The interaction of a centrifugal instability mechanism with the viscous mechanism which causes Tollmein-Schlichting waves is discussed. The interaction between these modes can be strong enough to drive the mean state; the interaction is investigated in the context of curved channel flows so as to avoid difficulties associated with boundary layer growth. Essentially it was found that the mean state adjusts itself so that any modes present are neutrally stable even at finite amplitude. In the first instance, the mean state driven by a vortex of short wavelength in the absence of Tollmein-Schlichting waves is considered. It is shown that for a given channel curvature and vortex wavelength, there is an upper limit to the mass flow rate which the channel can support as the pressure gradient is increased. When the Tollmein-Schlichting waves are present, then the nonlinear differential equation to determine the mean state is modified.

Bassom, Andrew P.↗
NASA NTRS2014

Flight Tests of a Supersonic Natural Laminar Flow Airfoil

IR thermography was used to characterize the transition front on a S-NLF test article at chord Reynolds numbers in excess of 30 million Changes in transition due to Mach number, Reynolds number, and surface roughness were investigated - Regions of laminar flow in excess of 80% chord at chord Reynolds numbers greater than 14 million IR thermography clearly showed the transition front and other flow features such as shock waves impinging upon the surface A series of parallel oblique shocks, of yet unknown origin, were found to cause premature transition at higher Reynolds numbers. NASA has a current goal to eliminate barriers to the development of practical supersonic transport aircraft Drag reduction through the use of supersonic natural laminar flow (S-NLF) is currently being explored as a means of increasing aerodynamic efficiency - Tradeoffs work best for business jet class at M<2 Conventional high-speed designs minimize inviscid drag at the expense of viscous drag - Existence of strong spanwise pressure gradient leads to crossflow (CF) while adverse chordwise pressure gradients amplifies and Tollmien-Schlichting (TS) instabilities Aerion Corporation has patented a S-NLF wing design (US Patent No. 5322242) - Low sweep to control CF - dp/dx < 0 on both wing surfaces to stabilize TS - Thin wing with sharp leading edge to minimize wave drag increase due to reduction in sweep NASA and Aerion have partnered to study S-NLF since 1999 Series of S-NLF experiments flown on the NASA F-15B research test bed airplane Infrared (IR) thermography used to characterize transition - Non-intrusive, global, good spatial resolution - Captures significant flow features well

laminar flow↗
NASA NTRS2008

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave transitionally in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗
NASA NTRS2008

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave 'transitionally' in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗
NASA NTRS1947

Some Recent Contributions to the Study of Transition and Turbulent Boundary Layers

The first part of this paper reviews the present state of the problem of the instability of laminar boundary layers which has formed an important part of the general lectures by von Karman at the first and fourth Congresses and by Taylor at the fifth Congress. This problem may now be considered as essentially solved as the result of work completed since 1938. When the velocity fluctuations of the free-stream flow are less than 0.1 percent of the mean speed, instability occurs as described by the well-known Tollmien-Schlichting theory. The Tollmien-Schlichting waves were first observed experimentally by Schubauer and Skramstad in 1940. They devised methods of introducing controlled small disturbances and obtained measured values of frequency, damping, and wave length at various Reynolds numbers which agreed well with the theoretical results. Their experimental results were confirmed by Liepmann. Much theoretical work was done in Germany in extending the Tol1mien-Schlichting theory to other boundary conditions, in particular to flow along a porous wall to which suction is applied for removing part of the boundary layer. The second part of this paper summarizes the present state of knowledge of the mechanics of turbulent boundary layers, and of the methods now being used for fundamental studies of the turbulent fluctuations in turbulent boundary layers. A brief review is given of the semi-empirical method of approach as developed by Buri, Gruschwitz, Fediaevsky, and Kalikhman. In recent years the National Advisory.Commsittee for Aeronautics has sponsored a detailed study at the National Bureau of Standards of the turbulent fluctuations in a turbulent boundary layer under adverse pressure gradient sufficient to produce separation. The aims of this investigation and its present status are described.

Dryden, Hugh L↗
NASA NTRS2020

Streak Instability Analysis for BOLT Configuration

The BOLT configuration designed for an upcoming flight experiment exemplifies the complexities of hypersonic boundary layer transition in the presence of multiple and potentially interacting instability mechanisms. A numerical study of the boundary layer instabilities over this configuration has been performed with the goal of applying state-of-the-art tools in multidimensional stability analysis that complement the other ongoing research related to the BOLT flight campaign. Specifically, the paper outlines modal instability characteristics of the boundary layer streaks adjacent to the minor-axis symmetry plane of a 1/3 rd scale BOLT configuration that was tested in the low-disturbance hypersonic wind tunnel at Purdue University and both quiet and conventional tunnels at Texas A&M University (TAMU). The basic state computations agree well with the heat flux measurements in the Purdue quiet tunnel as well as with the off-body velocity field described in a previous DNS study at the University of Minnesota. The encouraging agreement between the predicted frequency spectra of the streak instability in the region of increased boundary layer thickness and the surface pressure spectra measured during the TAMU experiments suggests that these streak instabilities may play an important role during the transition process near the symmetry plane. The modest value of the predicted peak N-factor (N max ≈ 5) for the streak instabilities provides a potential explanation for why transition is not observed under quiet flow conditions, but may occur in a noisy wind tunnel. Similar instability analyses are performed for the flight configuration at two selected flow conditions from the ascent portion of the planned flight trajectory. The predicted values of the peak N-factor suggest that transition onset in the vicinity of the center plane is possible at the ascent time of t = 26.98 seconds (N max ≈ 12.5 at X/L = 0.77), but the streak instabilities weaken sufficiently at the later time of t = 28.88 seconds (N max ≈ 11.5 at X/L ≈ 1.00) and are less likely to cause transition at that time.

Laminar turbulent transition↗
DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗