SEARCH · Search NASA
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.
MAESTROeX is a massively parallel, finite-volume C++/F90 solver for low Mach number astrophysical flows. Our code utilizes a low Mach number equation set allowing for more efficient, long-time integration of highly subsonic flows compared to compressible approaches. The recommended range of applicability is for flows where the Mach number does not exceed ~0.1.
Our objective is to study the collimation of solar jets by nonlinear forces corresponding to torsional Alfvén waves together with external forces. We consider a straight, initially non-rotating, untwisted magnetic cylinder embedded in a plasma with a straight magnetic field, where a shear between the internal and external flows exists. By implementing magnetohydrodynamic theory and taking into account the second-order thin flux tube approximation, the balance between the internal nonlinear forces is visualized. The nonlinear differential equation containing the ponderomotive, magnetic tension, and centrifugal forces in the presence of the shear flow is obtained. The solution presents the scale of influence of the propagating torsional Alfvén wave on compressive perturbations. Explicit expressions for the compressive perturbations caused by the forces connected to the torsional Alfvén wave show that, in the presence of a shear flow, the magnetic tension and centrifugal forces do not cancel each other’s effects as they did in its absence. This shear flow plays in favor of the magnetic tension force, resulting in a more efficient collimation. Regarding the ponderomotive force, the shear flow has no effect. The phase relations highlight the interplay of the shear flow and the plasma- β . As the shear flow and plasma- β increase, compressive perturbation amplitudes emerge. We conclude that the jet collimation due to the torsional Alfvén wave highly depends on the location of the jet. The shear flow tightens the collimation as the jet elevates up to the solar corona.
MAESTROeX is a massively parallel, finite-volume C++/F90 solver for low Mach number astrophysical flows. The code utilizes a low Mach number equation set allowing for more efficient, long-time integration of highly subsonic flows compared to compressible approaches. The recommended range of applicability is for flows where the Mach number does not exceed~0.1.
TetrUSS is a comprehensive suite of computational fluid dynamics (CFD) programs that won the Software of the Year award in 1996 and has found increasing use in government, academia, and industry for solving realistic flow problems (especially in aerodynamics and aeroelastics of aircraft having complex shapes). TetrUSS includes not only programs for solving basic equations of flow but also programs that afford capabilities for efficient generation and utilization of computational grids and for graphical representation of computed flows (see figure). The 2004 version of the Tetrahedral Unstructured Software System (TetrUSS), which is one of two software systems reported in "NASA s 2004 Software of the Year," NASA Tech Briefs, Vol. 28, No. 10 (October 2004), page 18, has been improved greatly since 1996. These improvements include (1) capabilities to simulate viscous flow by solving the Navier-Stokes equations on unstructured grids, (2) portability to personal computers from diverse manufacturers, (3) advanced models of turbulence, (4) a parallel-processing version of one of the unstructured-grid Navier-Stokes-equation-solving programs, and (5) advanced programs for generating unstructured grids.
The present research is aimed at providing a performance model for the Mars Helicopter (MH), to understand the complexity of the flow, and identify future regions of improvement. The low density of the Martian atmosphere and the relatively small MH rotor, result in very low chord-based Reynolds number flows. The low density and Reynolds numbers reduce the lifting force and lifting efficiency, respectively. The high drag coefficients in subcritical flow, especially for thicker sections, are attributed to laminar separation from the rear of the airfoil. In the absence of test data, efforts have been made to explore these effects using prior very low Reynolds number research efforts. The rotor chord-based Reynolds number range is observed to be subcritical, which makes boundary layer transition unlikely to occur. The state of the two-dimensional rotor boundary layer in hover is approximated by calculating the instability point, laminar separation point, and the transition location to provide understanding of the flow state in the high Mach-low Reynolds number regime. The results are used to investigate the need for turbulence modeling in Computational Fluid Dynamics (CFD) calculations afterwards. The goal is to generate a performance model for the MH rotor for a free wake analysis, because the computational budget for a complete Navier-Stokes solution for a rotating body-fitted rotor is substantial. In this study, a Reynolds-Averaged Navier-Stokes (RANS) based approach is used to generate the airfoil deck using C81Gen with stitched experimental data for very high angles of attack. A full Grid Resolution Study is performed and over 4,500 cases are completed to create the full airfoil deck. The laminar separation locations are predicted within the accuracy of the approximate method when compared with the CFD calculations. The model is presented through airfoil data tables (c81 files) that are used by comprehensive rotor analysis codes such as CAMRADII, or the mid-fidelity CFD solver RotCFD. Finally, the rotor performance is compared with experimental data from the 25ft Space Simulator at the NASA Jet Propulsion Laboratory (JPL) and shows good correlation for the rotor Figure of Merit over the available thrust range.
An efficient and user-oriented method has been constructed for calculating flow in and about complex inlet configurations. Efficiency is attained by: the use of a panel method, a technique of superposition for obtaining solutions at any inlet operating condition, and employment of an advanced matrix-iteration technique for solving large full systems of equations, including the nonlinear equations for the Kutta condition. User concerns are addressed by the provision of several novel graphical output options that, taken together, yield a more complete comprehension of the flowfield than had been possible previously. Examples of these features are presented for some complicated configurations, and where possible, comparisons are made between calculation and experiment.
Safe, reliable, and efficient operation of nuclear reactors requires timely and accurate measurement of primary coolant parameters such as temperature, flow, pressure, and level. These measurements may be challenging in small modular reactors (SMRs) with integral vessel design. For SMRs with the integral design, sensors may be located in confined spaces or embedded within pressure vessels. The absence of primary coolant loops is another fundamental design attribute that will affect the location of sensors, sensor interfaces, and types of sensors suited for process parameter measurements. Further, the number and location of pressure vessel penetrations may be constrained compared to traditional nuclear power plants (NPPs). As a result, conventional NPP temperature, flow, pressure, and level sensors may not be appropriate for use in SMRs. This paper addresses the development of new sensor technology for hybrid temperature, flow, and level measurement. The authors demonstrate the alternative use of resistance temperature detectors (RTDs) where cross-correlation of process temperature variations between RTDs is used to calculate process flow velocity. The new sensor and measurement technology will also incorporate RTD response time and self-heating characteristics to provide flow and level measurement in accident conditions. This capability would simplify reactor design and enable optimal placement of reactor pressure vessel penetrations leading to increased operational reliability and safety. The cross-correlation of temperature signals and adaptation for flow measurement in a laboratory setting is presented here. (authors)
Empirical gas–surface interaction models, such as the Maxwell model and the Cercignani–Lampis model, are widely used as the boundary condition in rarefied gas flow simulations. The accuracy of these models in the prediction of macroscopic behavior of rarefied gas flows is less satisfactory in some cases especially the highly non-equilibrium ones. Molecular dynamics simulation can accurately resolve the gas–surface interaction process at atomic scale, and hence can predict accurate macroscopic behavior. They are however too computationally expensive to be applied in real problems. In this work, a statistical physical-based gas–surface interaction model, which complies with the basic relations of boundary condition, is developed based on the framework of the washboard model. In virtue of its physical basis, this new model is capable of capturing some important relations/trends for which the classic empirical models fail to model correctly. As such, the new model is much more accurate than the classic models, and in the meantime is more efficient than MD simulations. Therefore, it can serve as a more accurate and efficient boundary condition for rarefied gas flow simulations. - Highlights: • The new physical-based gas–surface interaction model satisfies the basic relations of boundary condition so that it can be used in rarefied gas flow simulations. • The new physical model possesses several advantages comparing to the empirical gas–surface interaction models, hence it's a more accurate boundary condition for rarefied gas flow study. • The new physical model is a statistical model which can be efficiently implemented by Monte Carlo method. It's much more efficient comparing to the MD simulation.
A numerical method for simultaneously and efficiently coupling an external subsonic potential flow and an interior viscous flow such that the two flows match at an interfacing boundary is discussed. Both a panel method and a simple point compressible vortex model are used for the outer potential field. The interior flow solvers which were used are the Navier-Stokes and Euler codes of T. J. Coakley and the Euler code of A. Verhoff. In order to test compatibility, the panel method is coupled to the less expensive Euler codes since the coupling procedure is identical with the Navier-Stokes code. The results show significant efficiency improvements can be obtained over the uncoupled approach. Results also indicate the outer potential flow is best represented by the simple point compressible vortex model. The panel method couples smoothly to Coakley's implicit code but is numerically incompatible as coupled with the explicit Euler code. An improved Navier-Stokes code is under initial development which extends the Euler code to include the necessary viscous terms. Results are shown for all infinite length channel with one wavy periodic wall with and without laminar separation.
Researchers at the NASA Glenn Research Center have developed a new air-mass-flow sensor to solve the problems of existing mass flow sensor designs. NASA's design consists of thin-film resistors in a Wheatstone bridge arrangement. The resistors are fabricated on a thin, constant-thickness airfoil to minimize disturbance to the airflow being measured. The following photograph shows one of NASA s prototype sensors. In comparison to other air-mass-flow sensor designs, NASA s thin-film sensor is much more robust than hot wires, causes less airflow disturbance than pitot tubes, is more accurate than vane anemometers, and is much simpler to operate than thermocouple rakes. NASA s thin-film air-mass-flow sensor works by converting the temperature difference seen at each leg of the thin-film Wheatstone bridge into a mass-flow rate. The following figure shows a schematic of this sensor with air flowing around it. The sensor operates as follows: current is applied to the bridge, which increases its temperature. If there is no flow, all the arms are heated equally, the bridge remains in balance, and there is no signal. If there is flow, the air passing over the upstream legs of the bridge reduces the temperature of the upstream legs and that leads to reduced electrical resistance for those legs. After the air has picked up heat from the upstream legs, it continues and passes over the downstream legs of the bridge. The heated air raises the temperature of these legs, increasing their electrical resistance. The resistance difference between the upstream and downstream legs unbalances the bridge, causing a voltage difference that can be amplified and calibrated to the airflow rate. Separate sensors mounted on the airfoil measure the temperature of the airflow, which is used to complete the calculation for the mass of air passing by the sensor. A current application for air-mass-flow sensors is as part of the intake system for an internal combustion engine. A mass-flow sensor is used to provide accurate information about the amount of air entering the engine so that the amount of fuel can be adjusted to give the most efficient combustion. The ideal mass-flow sensor would be a rugged design that minimizes the disturbance to the flow stream and provides an accurate reading of both smooth and turbulent flows; NASA's design satisfies these requirements better than any existing design. Most of the mass-flow sensors used today are the hot wire variety. Hot wires can be fragile and cannot accurately measure a turbulent or reversing flow, which is often encountered in an intake manifold. Other types of mass-flow sensors include pitot tubes, vane anemometers, and thermocouple rakes-all of which suffer from some type of performance problem. Because it solves these performance problems while maintaining a simple design that lends itself to low-cost manufacturing techniques, NASA s thin-film resistance temperature detector air-mass-flow sensor should lead to more widespread use of mass-flow sensors.
The world's demand for energy is increasing dramatically, but the best energy conversion systems operate at approximately 30% efficiency. One way to decrease energy loss is in the recovery of waste heat using thermoelectric (TE) generators. A TE generator is device that generates electricity by exploiting heat flow across a thermal gradient. The efficiency of a TE material for power generation and cooling is determined by the dimensionless Figure of Merit (ZT): ZT = S(exp. 2)sigmaT/κ: where S is the Seebeck coefficient, sigma is the electrical conductivity, T is the absolute temperature, and κ is the thermal conductivity. The parameters are not physically independent, but intrinsically coupled since they are a function of the transport properties of electrons. Traditional research on TE materials has focused on synthesizing bulk semiconductor-type materials that have low thermal conductivity and high electrical conductivity affording ZT values of 1. The optimization of the σ/κ ratio is difficult to achieve using current material formats, as these material constants are complementary. Recent areas of research are focusing on using nanostructural artifacts that introduce specific dislocations and boundary conditions that scatter the phonons. This disrupts the physical link between thermal (phonon) and electrical (electron) transport. The result is that κ is decreased without decreasing σ. These material formats give ZT values of up to 2 which represent approximately 18% energy gain from waste heat recovery. The next challenge in developing the next generation of TE materials with superior performance is to tailor the interconnected thermoelectric physical parameters of the material system. In order to approach this problem, the fundamental physics of each parameter S, sigma, and κ need to be physically understood in their context of electron/phonon interaction for the construction of new high ZT thermoelectric devices. Is it possible to overcome the physical limit imposed by of the effect of phonon lattice oscillation and energetic electrons towards thermal conductivity? Is the Seebeck coefficient, based on the difference in voltage over temperature gradient ( deltaV/deltaT), an intrinsic parameter of each material? All these parameters were manipulated using nano-bridge and twin-lattice structural concepts at the NASA Langley Research Center. This talk will review the current trend of TE research to optimize the ZT and discuss about new approaches on increasing ZT within functional limits of each parameter.
Small-scale bed form topographies control hyporheic exchange and biogeochemical processes within aquatic sediments, which ultimately affect water quality and nutrient cycling at the watershed scale. The impact of three-dimensional and small-scale bed form topographies on hyporheic exchange and solute mixing is investigated in the present work. The effect of bed form morphologies on the development of zones of enhanced reaction rates (i.e., hotspots) is also studied. A computational fluid dynamics model to simulate river flow over bed forms is combined with a subsurface flow and multicomponent reactive solute transport model. A wide variety of bed form topographies are generated using geometric models by varying parameters controlling curvature as well as bed form wavelength and amplitude. The results in this research suggest that out-of-phase bed forms generate more complex hyporheic flow patterns which reduce the efficiency of solute transformations. Higher phase shifts in bed form shapes result in overall higher average velocity, larger zones of enhanced pressure and reaction rates, and higher amounts of solute exchange. Moreover, the bed form shapes control the reaction process for a wide range of sediment conductivities. Here this study advances the understanding of the effects of complex and small scale morphological features on hyporheic exchange processes including the rate and spatio-temporal distribution of reaction hotspots.
An advanced system analysis tool is being developed for advanced reactor safety analysis. This paper describes the underlying physics and numerical models used in the code, including the governing equations, the stabilization schemes, the high-order spatial and temporal discretization schemes, and the Jacobian Free Newton Krylov solution method. The effects of the spatial and temporal discretization schemes are investigated. Additionally, a series of verification test problems are presented to confirm the high-order schemes. Furthermore, it is demonstrated that the developed system thermal-hydraulics model can be strictly verified with the theoretical convergence rates, and that it performs very well for a wide range of flow problems with high accuracy, efficiency, and minimal numerical diffusions.
Molten salt thermal energy storage (TES) tanks ensure steady power output of concentrating solar power (CSP) plants; however, recent tank failures have highlighted the need for further analysis. Current studies primarily focus on analyzing the molten salt flow, heat transfer, and thermal efficiency. Additionally, research on the latest tank structures is limited and lacks newest experimental validation. This study measures temperature and molten salt inventory levels in the high-temperature tank at a 50 MW central receiver CSP plant, connected to the power grid in 2019. A multi-physics model was developed to evaluate thermal and mechanical properties of TES tanks by combining computational fluid dynamics and finite element modeling using real plant data. Heat loss, temperature, displacement, and stress distribution of the tank at different inventory levels were investigated. Results show that ambient air velocity near the tank roof reaches 2.14 m/s, much higher than 0.2 m/s near the wall. The temperatures of inventory fluid and tank are close, varying slightly at different levels due to thermal conduction and radiation. Because the heat loss strongly depends on temperature, the total tank loss remains nearly constant across inventory levels. Larger temperature gradients and thermal stresses are primarily localized along the tank floor edge and the air-salt interface. Notably, the maximum thermal stress at the tank edge is three times higher than that at the interface. The magnitude of total stress changes by less than 5 MPa with and without thermal load, indicating that high temperatures exert only a minor impact on tank stress. In contrast, thermal load significantly affects tank deformation, particularly at the roof edge, where values exceed 150 mm. Despite the large variation in molten salt levels, tank wall temperatures and displacements present a minor change, suggesting a weak correlation with inventory levels. In conclusion, the findings obtained in this study provide important insights on the TES tank that could be used to optimize tank design and operation strategies.
Aluminum particles are widely used as a metal fuel in solid propellants. However, poor combustion efficiencies and two-phase flow losses result due in part to particle agglomeration. Engineered composite particles of aluminum (Al) with inclusions of polytetrafluoroethylene (PTFE) or low-density polyethylene (LDPE) have been shown to improve ignition and yield smaller agglomerates in solid propellants, recently. Reductions in agglomeration were attributed to internal pressurization and fragmentation (microexplosions) of the composite particles at the propellant surface. We explore the mechanisms responsible for microexplosions in order to better understand the combustion characteristics of composite fuel particles. Single composite particles of Al/PTFE and Al/LDPE with diameters between 100 and 1200 µm are ignited on a substrate to mimic a burning propellant surface in a controlled environment using a CO 2 laser in the irradiance range of 78–7700 W/cm 2 . Furthermore, the effects of particle size, milling time, and inclusion content on the resulting ignition delay, product particle size distributions, and microexplosion tendencies are reported. For example particles with higher PTFE content (30 wt%) had laser flux ignition thresholds as low as 77 W/cm 2 , exhibiting more burning particle dispersion due to microexplosions compared to the other materials considered. Composite Al/LDPE particles exhibit relatively high ignition thresholds compared to Al/PTFE particles, and microexplosions were observed only with laser fluxes above 5500 W/cm 2 due to low LDPE reactivity with Al resulting in negligible particle self-heating. However, results show that microexplosions can occur for Al containing both low and high reactivity inclusions (LDPE and PTFE, respectively) and that polymer inclusions can be used to tailor the ignition threshold. Furthermore, this class of modified metal particles shows significant promise for application in many different energetic materials that use metal fuels.
Here, this study addresses a large-scale multimodal transit network design problem, with Shared Autonomous Mobility Services (SAMS) as both transit feeders and an origin-to-destination mode. The framework captures spatial demand and modal characteristics, considers intermodal transfers and express services, determines transit infrastructure investment and path flows, and generates transit routes. A system-optimal multimodal transit network is designed with minimum total door-to-door generalized costs of users and operators, satisfying transit origin-destination demand within a pre-set infrastructure budget. Firstly, the geography, demand, and modes in each zone are characterized with continuous approximation. The decisions of network link investment and multimodal path flows in zonal connection optimization are formulated as a minimum-cost multi-commodity network flow (MCNF) problem and solved efficiently with a mixed-integer linear programming (MILP) solver. Subsequently, the route generation problem is solved by expanding the MCNF formulation to minimize intramodal transfers. The model is illustrated through a set of experiments with the Chicago network comprised of 50 zones and seven modes, under three scenarios. The computational results present savings in traveler journey time and operator cost demonstrating the potential benefits of collaboration between multimodal transit systems and SAMS.
Electrochemical synthesis of ammonia via the nitrate reduction reaction (NO3RR) has been intensively researched as an alternative to the traditional Haber–Bosch process. Most research focuses on the low concentration range representative of the nitrate level in wastewater, leaving the high concentration range, which exists in nuclear and fertilizer wastes, unexplored. The use of a concentrated electrolyte (≥1 M) for higher rate production is hampered by poor hydrogen transfer kinetics. Herein, we demonstrate that a cocatalytic system of Ru/Cu 2 O catalyst enables NO3RR at 10.0 A in 1 M nitrate electrolyte in a 16 cm 2 flow electrolyzer, with 100% faradaic efficiency toward ammonia. Detailed mechanistic studies by deuterium labeling and operando Fourier transform infrared (FTIR) spectroscopy allow us to probe the hydrogen transfer rate and intermediate species on Ru/Cu 2 O. Ab initio molecular dynamics (AIMD) simulations reveal that adsorbed hydroxide on Ru nanoparticles increases the density of the hydrogen-bonded water network near the Cu 2 O surface, which promotes the hydrogen transfer rate. Finally, our work highlights the importance of engineering synergistic interactions in cocatalysts for addressing the kinetic bottleneck in electrosynthesis.
This is the second of two large-eddy simulation studies on the mechanisms of mesoscale cellular organization in drizzling (open cells) and nondrizzling marine stratocumulus (closed cells). This study uses a hard nudging approach which maintains fixed horizontal-mean temperature and humidity profiles for a well-mixed boundary layer with a constant boundary layer depth. For the case studied, closed cells develop and broaden by 32 hr to an aspect ratio of 25. Simulations show that the closed-cell mesoscale cellular convection is driven by positive feedback from cloud-induced mesoscale perturbations of longwave radiative cooling. A conceptual model for closed-cell stratocumulus as a mesoscale wavelength hydrodynamic instability in which mesoscale moist and dry anomalies spontaneously grow is presented. In simulations in which long-wavelength sinusoidal moisture anomalies are initially imposed, these anomalies evolve into amplifying closed cells. The cell structure is visualized with a compositing approach based on sorting grid columns by their mesoscale-smoothed total water path. A thermally direct mesoscale circulation pattern develops in the interior of the boundary layer with buoyant mesoscale updrafts, thicker cloud, and a slightly higher capping inversion in the moister columns. There is a mesoscale flow of above-inversion air down the slightly sloping capping inversion from the moist to the dry regions, reinforced by cloud top radiative cooling. This strengthens the mesoscale anomalies by preferentially cooling and drying the already dry regions. The sloping inversion flow is not driven as efficiently if the radiative cooling is artificially horizontally homogenized, partly disrupting this positive feedback and the resulting closed-cell development.