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.
An explicit but simple lumped parameter nonlinear multivariable model of a LN2-GN2-operated closed circuit cryogenic wind tunnel has been developed and its basic features have been experimentally validated. The model describes the mass-energy interaction involved in the cryogenic tunnel process and includes the real gas properties of nitrogen gas.
A simple lumped parameter based inverse design tool has been developed which provides flow path geometry and entrainment estimates subject to operational, acoustic, and design constraints. These constraints are manifested through specification of primary mass flow rate or ejector thrust, fully-mixed exit velocity, and static pressure matching. Fundamentally, integral forms of the conservation equations coupled with the specified design constraints are combined to yield an easily invertible linear system in terms of the flow path cross-sectional areas. Entrainment is computed by back substitution. Initial comparison with experimental and analogous one-dimensional methods show good agreement. Thus, this simple inverse design code provides an analytically based, preliminary design tool with direct application to High Speed Civil Transport (HSCT) design studies.
The equations of motion for a flexible vehicle capable of arbitrary translational and rotational motions in inertial space accompanied by small elastic deformations are derived in an unabridged form. The vehicle is idealized as consisting of a single rigid body with an ensemble of mass particles interconnected by massless elastic structure. The internal elastic restoring forces are quantified in terms of a stiffness matrix. A transformation and truncation of elastic degrees of freedom is made in the interest of numerical integration efficiency. Deformation dependent terms are partitioned into a hierarchy of significance. The final set of motion equations are brought to a fully assembled first order form suitable for direct digital implementation. A FORTRAN program implementing the equations is given and its salient features described.
Nontopological solitons are stable field configurations which may be formed in a primordial phase transition. Their cosmic evolution is studied, and the possibility that such objects could contribute significantly to the energy density of the Universe is examined. As the Universe cools, initially all but the largest lumps evaporate into free particles; those which survive may subsequently enter a brief accretion phase before they freeze out at a final size. Although the minimum critical charges which survive depend on particle masses and couplings, researchers develop an analysis which applies to a wide class of models. In most cases, solitons of moderate size survive the evaporation process only if there is a significant charge asymmetry or if they form at a temperature well below their binding energy per charge.
Bragg scattering in periodic media generates bandgaps, frequency bands where waves attenuate rather than propagate. Yet, a finite periodic structure may exhibit resonance frequencies within these bandgaps. This is caused by boundary effects introduced by the truncation of the nominal infinite medium. Previous studies of discrete systems determined existence conditions for bandgap resonances, although the focus has been limited to mainly periodic chains with free–free boundaries. In this paper, we present closed-form existence conditions for bandgap resonances in discrete diatomic chains with general boundary conditions (free–free, free–fixed, fixed–free, or fixed–fixed), odd or even chain parity (contrasting or identical masses at the ends), and the possibility of attaching a unique component (mass and/or spring) at one or both ends. The derived conditions are consistent with those theoretically presented or experimentally observed in prior studies of structures that can be modeled as linear discrete diatomic chains with free–free boundary conditions. An intriguing case is a free–free chain with even parity and an arbitrary additional mass at one end of the chain. Introducing such an arbitrary mass underscores a transition among a set of distinct existence conditions, depending on the type of chain boundaries and parity. The proposed analysis is applicable to linear periodic chains in the form of lumped-parameter models, examined across the frequency spectrum, as well as continuous granular media models, or similar configurations, examined in the low-frequency regime.
This paper presents the salient features of a preliminary thermal-structural design and analysis study of a hydrogen-fueled, regeneratively cooled, airframe-integrated scramjet. This three-dimensional fixed geometry scramjet concept which was developed at NASA Langley Research Center is designed to operate over a flight Mach number range from 4 to 10. The thermal-structural study was focused on a scramjet application to one concept for a hypersonic research vehicle and was based on technology developed under the NASA Hypersonic Research Engine Project. State-of-the-art analytical methods consisting of lumped system and finite difference steady-state thermal analyses and a finite element structural analysis were used. The results of the study indicated that this scramjet concept is viable from both a structural mass and cooling requirement standpoint. However, advances in material and fabrication technology for hydrogen-cooled structures appear necessary for acceptable engine life for commercial application.
NASA design teams have been investigating options for "detuning" Ares I to prevent oscillations originating in the vehicle solid-rocket main stage from synching up with the natural resonance of the rest of the vehicle. An experimental work started at NASA MSFC center in 2008 using a damping device showed great promise in damping the vibration level of an 8 resonant tank. However, the mechanisms of the vibration damping were not well understood and there were many unknowns such as the physics, scalability, technology readiness level (TRL), and applicability for the Ares I vehicle. The objectives of this study are to understand the physics of intriguing slosh damping observed in the experiments, to further validate a Computational Fluid Dynamics (CFD) software in propellant sloshing against experiments with water, and to study the applicability and efficiency of the slosh damper to a full scale propellant tank and to cryogenic fluids. First a 2D fluid-structure interaction model is built to model the system resonance of liquid sloshing and structure vibration. A damper is then added into the above model to simulate experimentally observed system damping phenomena. Qualitative agreement is found. An analytical solution is then derived from the Newtonian dynamics for the thrust oscillation damper frequency, and a slave mass concept is introduced in deriving the damper and tank interaction dynamics. The paper will elucidate the fundamental physics behind the LOX damper success from the derivation of the above analytical equation of the lumped Newtonian dynamics. Discussion of simulation results using high fidelity multi-phase, multi-physics, fully coupled CFD structure interaction model will show why the LOX damper is unique and superior compared to other proposed mitigation techniques.
MOTIVATION: Spaceflight countermeasures mitigate the harmful effects of the space environment on astronaut health and performance. Exercise has historically been used as a countermeasure to physical deconditioning, and additional countermeasures including lower body negative pressure, blood flow occlusion and artificial gravity are being researched as countermeasures to spaceflight-induced fluid shifts. The NASA Digital Astronaut Project uses computational models of physiological systems to inform countermeasure design and to predict countermeasure efficacy.OVERVIEW: Computational modeling supports the development of the exercise devices that will be flown on NASAs new exploration crew vehicles. Biomechanical modeling is used to inform design requirements to ensure that exercises can be properly performed within the volume allocated for exercise and to determine whether the limited mass, volume and power requirements of the devices will affect biomechanical outcomes. Models of muscle atrophy and bone remodeling can predict device efficacy for protecting musculoskeletal health during long-duration missions. A lumped-parameter whole-body model of the fluids within the body, which includes the blood within the cardiovascular system, the cerebral spinal fluid, interstitial fluid and lymphatic system fluid, estimates compartmental changes in pressure and volume due to gravitational changes. These models simulate fluid shift countermeasure effects and predict the associated changes in tissue strain in areas of physiological interest to aid in predicting countermeasure effectiveness. SIGNIFICANCE: Development and testing of spaceflight countermeasure prototypes are resource-intensive efforts. Computational modeling can supplement this process by performing simulations that reduce the amount of necessary experimental testing. Outcomes of the simulations are often important for the definition of design requirements and the identification of factors essential in ensuring countermeasure efficacy.
Axions are a promising cold dark matter candidate. Haloscopes, which use the conversion of axions to photons in the presence of a magnetic field to detect axions, are the basis of microwave cavity searches such as the Axion Dark Matter eXperiment (ADMX). To search for lighter, low frequency axions in the sub-2 × 10 –7 eV (50 MHz) range, a tunable lumped-element $\textit{LC}$ circuit has been proposed. For the first time, through ADMX SLIC (Superconducting $\textit{LC}$ Circuit Investigating Cold Axions), a resonant $\textit{LC}$ circuit was used to probe this region of axion mass-coupling space. The detector used a superconducting $\textit{LC}$ circuit with piezoelectric driven capacitive tuning. The axion mass and corresponding frequency ranges 1.7498–1.7519 × 10 –7 eV (42.31–42.36 MHz), 1.7734–1.7738 × 10 –7 eV (42.88–42.89 MHz), and 1.8007–1.8015 × 10 –7 eV (43.54–43.56 MHz) were covered at magnetic fields of 4.5 T, 5.0 T, and 7.0 T, respectively. Exclusion results from the search data, for coupling below 10 –12 GeV –1 , are presented.
Parameters for a lumped linear model approximating a distributed elastic structure are determined from dynamic test data comprised of several mode shapes and frequencies. Measurement errors, nonlinear response, and nonmeasurable quantities such as mode slope components are accommodated. Some mass and stiffness parameters may be known accurately, whereas the remainder are to be estimated. The method entails minimizing a quadratic function of the difference between corresponding modes and frequencies of the theoretical model and the test specimen. This technique was applied to some actual vibration test data, and the special techniques that are required to overcome convergence problems are described.
OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.
We developed an CFD aided ML-based tool for rapid assessment and optimization of different compositions of biomass in a fluidized bed reactor. First, we use CFD to simulate fluidized bed reactors with known inlet biomass mixtures and obtain corresponding syngas yields. A lumped kinetic mechanisms represents the conversion of cellulose, hemicellulose, and lignin, as well as subsequent cracking of tars into non-condensable gases (H 2 , CO, CO 2 , CH 4 ). We use Bayesian analysis/optimization to obtain the ideal operational and mass flow conditions for hydrogen production.
One-dimensional oscillating flow models are very useful for designing pulse tubes. They are simple to use, not computationally intensive, and the physical relationship between temperature, pressure and mass flow are easy to understand when used in conjunction with phasor diagrams. They do not possess, however, the ability to directly calculate thermal and momentum diffusion in the direction transverse to the oscillating flow. To account for transverse effects, lumped parameter corrections, which are obtained though experiment, must be used. Or two-dimensional solutions of the differential fluid equations must be obtained. A linear two-dimensional solution to the fluid equations has been obtained. The solution provides lumped parameter corrections for one-dimensional models. The model accounts for heat transfer and shear flow between the gas and the tube. The complex Nusselt number and complex shear wall are useful in describing these corrections, with phase relations and amplitudes scaled with the Prandtl and Valensi numbers. The calculated ratio, a, between a two-dimensional solution of the oscillating temperature and velocity and a one-dimensional solution for the same shows a scales linearly with Va for Va less than 30. In this region alpha less than 0.5, that is, the enthalpy flow calculated with a two-dimensional model is 50% of a calculation using a one-dimensional model. For Va greater than 250, alpha = 0.8, showing that diffusion is still important even when it is confined to a thing layer near the tube wall.
The Balloon Experiment for Galactic INfrared Science (BEGINS) is a concept for a sub-orbital observatory that will characterize dust in the vicinity of high-mass stars. It requires detectors with NEPs from 2x10 -16 W Hz -1/2 to 6x10 -17 W Hz -1/2 from 25-250 microns, respectively. The mission’s sensitivity requirements will be met by utilizing arrays of 1,840 lens-coupled, lumped-element kinetic inductance detectors (KIDs) operating at 300 mK. Each KID will consist of a titanium nitride (TiN) parallel strip absorbing inductive section and parallel plate capacitor deposited on a Silicon (Si) substrate. The parallel plate capacitor geometry allows for reduction of the pixel spacing. We present the optical performance of a prototype BEGINS KID array at 25 microns when coupled to Fresnel zone plate lenses. We present sensitivity, optical efficiency and quasiparticle lifetime measurements taken in a cryogenic test bed with a cryogenic blackbody.
The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers. It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages. A common method for topology optimization is that of compliance minimization which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system. Sherrnan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this. SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.
Spatially distributed prediction of streamflow and nitrogen export dynamics is essential for precision management of agricultural watersheds. While temporal deep learning models such as Long Short-Term Memory (LSTM) have shown strong performance at basin scales, their ability to generalize spatially is limited by insufficient representation of spatial dependencies and flow paths, particularly under data-scarce conditions. To address this gap, we propose HydroGraphNet, a knowledge-guided graph machine learning framework that integrates process-based knowledge and explicit spatial learning into temporal modeling. This framework incorporates directed graph topology to encode watershed connectivity and upstream inflows, with mass balance constraints to improve physical consistency. To enhance generalization in sparsely monitored regions, HydroGraphNet is pretrained on synthetic data generated by the SWAT+ (Soil and Water Assessment Tool Plus) model. We evaluated HydroGraphNet in the Upper Sangamon River Basin (44 HUC-12 subwatersheds, 2001–2020) against two LSTM baselines: a lumped basin-level model and a distributed variant. When benchmarked on SWAT+ simulations in pretraining, HydroGraphNet improved test NSEs by 8.9% (discharge) and 13.7% (NO₃–N load) in temporal extrapolation, and by 27.1% and 34.7% in spatial extrapolation, relative to the Lumped LSTM baseline. After fine-tuning with USGS monitoring data, the model achieved mean test NSE (KGE) scores of 0.768 (0.861) for discharge and 0.626 (0.664) for NO₃–N load, substantially outperforming baselines. Attribution analysis further highlighted the importance of upstream inflow representation and graph-based spatial learning in capturing cross-subwatershed dependencies. The model also reproduced seasonal hydrological and biogeochemical patterns consistent with known processes, demonstrating its robustness and process fidelity for spatially distributed prediction. Altogether, HydroGraphNet advances the integration of physical knowledge and spatially explicit learning in hydrological modeling, offering a generalizable framework for distributed modeling to support spatially targeted water quality management in data-scarce watersheds.
The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system.. Shennan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this (Akgun et al., 1998b). SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.
The author has previously shown that a considerable amount of clumping of interstellar grains is likely to take place during the free-fall collapse phase of an interstellar cloud which is forming the primitive solar nebula, with the assumption of sonic turbulence in the gas. The original estimate involved the crude assumption of hierarchal amalgamation of the grains upon collision. A Monte Carlo simulation of this process confirmed the general features of the results, but it was further found that the introduction of a low sticking probability reduced the size of the lumps quite significantly. A more realistic calculation was therefore carried out in which it was assumed that clumps of grains would tend to stick together if their collisions were approximately head-on, but that they would tend to fragment into smaller pieces if the collisions were more tangential. For typical values of the amalgamation parameter, this tends to spread the mass of the interstellar grains over a wide range of clump sizes, ranging from individual grains to objects in the millimeter or centimeter size.