Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,675 records · Page 93

A Multi-Dimensional Heat Transfer Model of a Tie-Tube and Hexagonal Fuel Element for Nuclear Thermal Propulsion

The Space Capable Cryogenic Thermal Engine (SCCTE) effort considers a nuclear thermal rocket design based around a Low-Enriched Uranium (LEU) design fission reactor. The reactor core is comprised of bundled hexagonal fuel elements that directly heat hydrogen for expansion in a thrust chamber and hexagonal tie-tubes that house zirconium hydride moderator mass for the purpose of thermalizing fast neutrons resulting from fission events. Created 3D steady state Hex fuel rod model with 1D flow channels. Hand Calculation were used to set up initial conditions for fluid flow. The Hex Fuel rod uses 1D flow paths to model the channels using empirical correlations for heat transfer in a pipe. Created a 2-D axisymmetric transient to steady state model using the CFD turbulent flow and Heat Transfer module in COMSOL. This model was developed to find and understand the hydrogen flow that might effect the thermal gradients axially and at the end of the tie tube where the flow turns and enters an annulus. The Hex fuel rod and Tie tube models were made based on requirements given to us by CSNR and the SCCTE team. The models helped simplify and understand the physics and assumptions. Using pipe correlations reduced the complexity of the 3-D fuel rod model and is numerically more stable and computationally more time-efficient compared to the CFD approach. The 2-D axisymmetric tie tube model can be used as a reference "Virtual test model" for comparing and improving 3-D Models.

Gomez, C. F.↗

Field Validation of Wake Steering Control with Wind Direction Variability: Preprint

Wake steering is a type of wind farm control strategy wherein upstream turbines are misaligned with the wind direction to redirect their wakes away from downstream turbines, increasing overall wind plant power. Wake steering is often analyzed assuming steady mean wind directions across the wind farm. However, in practice, the wind direction varies considerably over time because of large-scale weather phenomena. Wind direction variability causes wake steering controllers to increase power production less than predicted by steady state models, but more robust wake steering strategies can be designed that account for variable wind conditions. This paper compares the achieved yaw offsets and power gains from two 2-turbine wake steering experiments at a commercial wind farm with model predictions using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented model, assuming both fixed and variable wind directions. The impact of wind direction variability is modeled by including wind direction and yaw uncertainty in the FLORIS calculations. The field results match the trends predicted assuming wind direction variability. Specifically, the yaw offsets achieved in the intended control regions are lower than desired, resulting in less power gain, while a slight loss in power occurs for wind directions outside of the intended control region because of unintentional yaw misalignment. The agreement between the model and field results suggests that the wind direction variability model can be used to design wake steering controllers that are more robust to variable wind conditions present in the field.

17 WIND ENERGY↗

Newton-Okounkov bodies of chemical reaction systems

Despite their noted potential in polynomial-system solving, there are few concrete examples of Newton-Okounkov bodies arising from applications. Accordingly, in this paper, we introduce a new application of Newton-Okounkov body theory to the study of chemical reaction networks and compute several examples. An important invariant of a chemical reaction network is its maximum number of positive steady states Here, we introduce a new upper bound on this number, namely the ‘Newton-Okounkov body bound’ of a chemical reaction network. Further, through explicit examples, we show that the Newton-Okounkov body bound of a network gives a good upper bound on its maximum number of positive steady states. We also compare this Newton-Okounkov body bound to a related upper bound, namely the mixed volume of a chemical reaction network, and find that it often achieves better bounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Testing the DIII-D co/counter off-axis neutral beam injected power and ability to balance injected torque

DIII-D has undergone a major upgrade and successfully injected high power off-axis neutral beams (~4 MW) in both co-current and counter-current directions. This capability of high power co/counter steerable off-axis neutral beams on a major tokamak opens a unique parameter space of broad pressure and current profiles for high beta steady-state advanced tokamak (AT) scenarios, while retaining the ability to balance the injected torque for low rotation studies. This co/counter off-axis neutral beam capability is being used to validate physics-based energetic particle and thermal transport models for designing next-step facilities based on the steady-state AT approach. This paper reports on the critical evaluation of the transmitted power and energetic ion population produced by this heating and current drive system, which is assessed through visible imaging, neutron measurements and rotation profile measurements at balanced torque. Minimal losses of neutral beam power have been achieved by optimizing the strongly focused ion sources required to pass through the aperture. Tilting of the ion source has been guided by fast visible imaging and resulted in neutral beam injection along the design centerline with empirical characterization of each beam's divergence derived from the imaging data and used in the NUBEAM description of the beam injection. Through exclusive power injection of each neutral beam into MHD quiescent plasmas across a range of neutral beam voltage, perveance and plasma current we conclude that a modest reduction (~10%–15%) of transmitted power compared to on-axis, standard focus has been incurred. We report corrections that more accurately represent the injected power. Good ability to balance the neutral beam torque has been demonstrated by injecting the new off-axis counter injecting beam against the existing off-axis co-injected beam in 2.0 T, 1.0 MA, MHD quiescent L-mode plasmas. Furthermore, the torque balance studies verify the ability to operate with balanced injection, which is critical for achieving low torque and low rotation operation for physics studies and in ITER demonstration discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mariner Jupiter/Saturn LCSSE thruster/valve assembly and injection propulsion unit rocket engine assemblies: 0.2-lbf T/VA development and margin limit test report

Thruster valve assemblies (T/VA's) were subjected to the development test program for the combined JPL Low-Cost Standardized Spacecraft Equipment (LCSSE) and Mariner Jupiter/Saturn '77 spacecraft (MJS) programs. The development test program was designed to achieve the following program goals: (1) demonstrate T/VA design compliance with JPL Specifications, (2) to conduct a complete performance Cf map of the T/VA over the full operating range of environment, (3) demonstrate T/VA life capability and characteristics of life margin for steady-state limit cycle and momentum wheel desaturation duty cycles, (4) verification of structural design capability, and (5) generate a computerized performance model capable of predicting T/VA operation over pressures ranging from 420 to 70 psia, propellant temperatures ranging from 140 F to 40 F, pulse widths of 0.008 to steady-state operation with unlimited duty cycle capability, and finally predict the transient performance associated with reactor heatup during any given duty cycle, start temperature, feed pressure, and propellant temperature conditions.

Clark, E. C.↗

The Eddington limit and supercritical accretion. I - Time-independent calculations

Spherically symmetric, steady state accretion of an ionized hydrogen plasma onto a neutron star is considered for accretion rates which exceed a critical rate at which the Eddington limiting luminosity is produced. The coupled hydrodynamic and frequency integrated, radiative transfer equations are solved for accretion rates up to 10 times the nominal limit. Steady state solutions are presented that imply a multiplicity of different luminosity solutions for a single accretion rate in this 'supercritical' regime.

Burger, H. L.↗

A theoretical model for the cross spectra between pressure and temperature downstream of a combustor

A theoretical model developed to calculate pressure-temperature cross spectra, pressure spectra, temperature spectra and pressure cross spectra in a ducted combustion system is presented. The model assumes the presence of a fluctuating-volumetric-heat-release-rate disk source and takes into account the spatial distribution of the steady-state volumetric-heat flux. Using the model, pressure, velocity, and temperature perturbation relationships can be obtained. The theoretical results show that, at a given air mass flow rate, the calculated pressure-temperature cross spectra phase angle at the combustor exit depends on the model selected for the steady-state volumetric-heat flux in the combustor. Using measurements of the phase angle, an appropriate source region model was selected. The model calculations are compared with the data. The comparison shows good agreement and indicates that with the use of this model the pressure-temperature cross spectra measurements provide useful information on the physical mechanisms active at the combustion noise source.

Miles, J. H.↗

Time-marching methods for three-dimensional steady and unsteady viscous imcompressible flows

An implicit algorithm for the solution of three-dimensional, steady and unsteady, viscous, incompressible flows is presented. The algorithm is based on an upwind-relaxation finite-difference method. Steady-state solutions are carried out using a time-marching solution technique in combination with a local time-stepping strategy. To obtain time-accurate solutions, a subiterative procedure is employed at each physical time step using a global time step to ensure the divergence-free condition. Steady-state flows in several straight ducts and in a square duct with a 90-degree bend are computed and compared with analytical and experimental results. The classical problem of starting flow in a circular pipe is chosen to verify the time accuracy of the present scheme. Finally, the three-dimensional bubble-type vortex breakdown of a slender cylindrical vortex in an unbounded flow is investigated.

Hsu, C.-H.↗

Reconnection in Planetary Magnetospheres

Current sheets in planetary magnetospheres that lie between regions of "oppositely-directed" magnetic field are either magnetopause-like, separating plasmas with different properties, or tail-like, separating plasmas of rather similar properties. The magnetopause current sheets generally have a nearly limitless supply of magnetized plasma that can reconnect, possibly setting up steady-state reconnection. In contrast, the plasma on either side of a tail current sheet is stratified so that, as reconnection occurs, the plasma properties, in particular the Alfven velocity, change. If the density drops and the magnetic field increases markedly perpendicular to the sheet, explosive reconnection can occur. Even though steady state reconnection can take place at magnetopause current sheets, the process often appears to be periodic as if a certain low average rate was demanded by the conditions but only a rapid rate was available. Reconnection of sheared fields has been postulated to create magnetic ropes in the solar corona, at the Earth's magnetopause, and in the magnetotail. However, this is not the only way to produce magnetic ropes as the Venus ionosphere shows. The geometry of the reconnecting regions and the plasma conditions both can affect the rate of reconnection. Sorting out the various controlling factors can be assisted through the examination of reconnection in planetary settings. In particular we observe similar small-scale tearing in the magnetopause current layers of the Earth, Saturn. Uranus and Neptune and the magnetodisk current sheet at Jupiter. These sites may be seeds for rapid reconnection if the reconnection site reaches a high Alfven velocity region. In the Jupiter magnetosphere this appears to be achieved with resultant substorm activity. Similar seeds may be present in the Earth's magnetotail with the first one to reach explosive growth dominating the dynamics of the tail.

Russell, C. T.↗

Experimental and Numerical Analysis of Additively Manufactured Inconel 718 Coupons With Lattice Structure

Abstract In the present paper, two lattice geometries suitable for near surface and double wall cooling were developed and tested. The first type of unit cell consisted of six ligaments of 0.5 mm diameter joined at a common vertex near the middle. The second type of unit cell was derived from the first type by adding four mutually perpendicular ligaments in the middle plane. Two lattice configurations, referred to as L1 and L2, respectively, were obtained by repeating the corresponding unit cell in streamwise and spanwise directions in an inline fashion. Test coupons consisting of these lattice geometries embedded inside rectangular cooling channel with dimensions of 2.54 mm height, 38.07 mm width, and 38.1 mm in length were fabricated using Inconel 718 powder and selective laser sintering (SLS) process. The heat transfer and pressure drop performance was then evaluated using steady-state tests with constant wall temperature boundary condition and for channel Reynolds number ranging from 2800 to 15,000. The lattices depicted a higher heat transfer compared with a smooth channel and both the heat transfer and pressure drop increased with a decrease in the porosity from L1 to L2. Steady-state conjugate numerical results revealed formation of prominent vortical structures in the inter-unit cell spaces, which diverted the flow toward the top end wall and created an asymmetric heat transfer between the two end walls. In conclusion, these lattice structures provided an augmented heat transfer while favorably redistributing the coolant within channel.

Engineering↗

Development of a Pre-Combustion CO 2 Capture Process Using High-Temperature PBI Hollow-Fiber Membranes

The overall objective of this project was to evaluate the advantages of transformational polybenzimidazole (PBI) polymer hollow-fiber membrane (HFM)-based, carbon dioxide (CO 2 ) capture and purification technology at bench-scale using an actual coal-derived syngas stream from a coal gasification facility. The project was carried out over two budget periods. The technical objectives in Budget Period 1 (BP1) included preparing HFs and modules and upgrading the available skid for field testing. The technical objectives for BP2 were to field-test the skid unit with actual coal-derived syngas from an oxygen-blown gasifier to obtain performance data, update the Techno-Economic Analysis (TEA) that would assist with future process scale-up, and provide information on the design of a small pilot-scale test unit. The goal was to advance the PBI-HFM CO 2 capture and gas separation system for pre-combustion applications beyond second-generation economic performance predictions and make progress toward meeting overall fossil energy performance goals of CO 2 capture with 95% CO 2 purity at a cost of electricity (COE) 30% less than baseline capture approaches. The research program was designed with progressive technical tasks leading to both dynamic and steady-state testing of the PBI-HFM skid with actual coal-derived syngas. The work plan was to: (1) fabricate sufficient Generation-2 (GEN-2) fibers for module fabrication; (2) upgrade the fiber skid to accommodate large fiber modules for bench-scale field testing; (3) conduct dynamic and steady-state testing with coal-derived syngas from an oxygen-blown gasifier and obtain system performance data; (4) perform a TEA and environmental, health, and safety (EH&S) assessment; (5) update the State-Point Data Table, Technology Gap Analysis (TGA), and Technology Maturation Plan (TMP); (6) uninstall and return the test skid to the Recipient’s facilities; and (7) submit a Final Report that describes the results and analysis of the project research effort.

03 NATURAL GAS↗

FAST-1.2.2: Integral Assessment

This integral assessment quantifies the predictive capabilities of Fuel Analysis under Steady-state and Transients (FAST), a thermal-mechanical nuclear fuel performance code designed to analyze fuel behavior from beginning of life to burnup levels allowed by the U.S. Nuclear Regulatory Commission. FAST code calculations are shown to compare satisfactorily to a pre-selected set of experimental data with both steady-state and anticipated operating occurrence conditions and design basis accident transient operating conditions. This document describes the assessment of FAST-1.2.2, the latest version of FAST, released April 2026.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Generalized dynamic engine simulation techniques for the digital computers

Recently advanced simulation techniques have been developed for the digital computer and used as the basis for development of a generalized dynamic engine simulation computer program, called DYNGEN. This computer program can analyze the steady state and dynamic performance of many kinds of aircraft gas turbine engines. Without changes to the basic program, DYNGEN can analyze one- or two-spool turbofan engines. The user must supply appropriate component performance maps and design point information. Examples are presented to illustrate the capabilities of DYNGEN in the steady state and dynamic modes of operation. The analytical techniques used in DYNGEN are briefly discussed, and its accuracy is compared with a comparable simulation using the hybrid computer. The impact of DYNGEN and similar digital programs on future engine simulation philosophy is also discussed.

Sellers, J.↗

Parametric dependence of ion temperature and relative density in the NASA Lewis SUMMA facility

Further hot-ion plasma experiments were conducted in the SUMMA superconducting magnetic mirror facility. A steady-state ExB plasma was formed by applying a strong radially inward dc electric field between cylindrical anodes and hollow cathodes located near the magnetic mirror maxima. Extending the use of water cooling to the hollow cathodes, in addition to the anodes, resulted in higher maximum power input to the plasma. Steady-state hydrogen plasmas with ion kinetic temperatures as high as 830 eV were produced. Functional relations were obtained empirically among the plasma current, voltage, magnetic flux density, ion temperature, and relative ion density. The functional relations were deduced by use of a multiple correlation analysis. Data were obtained for midplane magnetic fields from 0.5 to 3.37 tesla and input power up to 45 kW. Also, initial absolute electron density measurements are reported from a 90 deg Thomson scattering laser system.

Snyder, A.↗

A Report from the Thermal Science Research Center (TSRC)

A vertical flow loop was designed and assembled to determine the local (circumferential and axial) and mean wall temperature distributions for single-phase and two-phase (subcooled and saturated) downward flow in both uniformly-heated and single-side heated vertical channels. Freon-11 was used as the working fluid in order to directly relate and compare the results with a previous experimental campaign which employed this same working fluid. For a given steady-state experiment, the following parameters were held constant: (1) exit pressure, (2) inlet temperature, and (3) mass velocity. For a given configuration of the 2.2 m long cylindrical channel test section, which had a 1.2 m long heated section, the applied heat rate was varied from zero through successive quasi-steady states to a level which corresponded to localized film boiling in the test section. The measurements showed that the boiling curve changes significantly at higher mass velocities with respect to both the circumferential and axial directions. The slope of the boiling curve changes in a non-monotonic fashion with respect to the circumferential directions. The slope of the boiling curve changes in a non-monotonic fashion with respect to the circumferential direction. The measurements point to the existence of a dry-out phenomenon occurring at multiple levels of the applied heat for the single-side heated channel. In comparing the heat transfer for horizontal channel flow with a vertically downward flow, the results show that significantly lower heat transfer occurs in the horizontal flow. However, this trend reverses as both the Reynolds number and the applied heat rate increase. Both the Liu-Winterton and Shah correlations were compared with the experimental data. The Shah correlation predicted the uniformly heated tube data better. When a thermal hydraulic diameter approach was used for the single-side heated case, the data at upstream locations for Z/L less than 0.5 was bounded above by the Liu-Winterton correlation and below by the Shah correlation. At Z/L = 0.5, the Shah correlation bounded the data; and for Z/L greater than 0.5, both correlations overpredicted the data with the Shah correlation being closest to the data. The present results indicate that additional correlational development is needed. In addressing some of the advanced space thermal management objectives concerning accommodating high heat fluxes in non-uniformly heated systems, a large battery of experiments 88 have been completed where local two-dimensional wall temperature variations were measured for both single-phase and two-phase flow in a single-side heated circular tube. As noted above, the results show significant axial and circumferential variations. Accurately accounting for such variations can result in optimized future advanced space, enhanced (high heat flux) thermal management systems.

Boyd, Ronald D.↗

Higher Order Time Integration Schemes for the Unsteady Navier-Stokes Equations on Unstructured Meshes

The rapid increase in available computational power over the last decade has enabled higher resolution flow simulations and more widespread use of unstructured grid methods for complex geometries. While much of this effort has been focused on steady-state calculations in the aerodynamics community, the need to accurately predict off-design conditions, which may involve substantial amounts of flow separation, points to the need to efficiently simulate unsteady flow fields. Accurate unsteady flow simulations can easily require several orders of magnitude more computational effort than a corresponding steady-state simulation. For this reason, techniques for improving the efficiency of unsteady flow simulations are required in order to make such calculations feasible in the foreseeable future. The purpose of this work is to investigate possible reductions in computer time due to the choice of an efficient time-integration scheme from a series of schemes differing in the order of time-accuracy, and by the use of more efficient techniques to solve the nonlinear equations which arise while using implicit time-integration schemes. This investigation is carried out in the context of a two-dimensional unstructured mesh laminar Navier-Stokes solver.

Jothiprasad, Giridhar↗