Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer model”

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,117 records · Page 62

Modeling the Effects of Liquid-to-Gas Density Ratio on Slosh Dynamics

Mechanical models are commonly used in Guidance, Navigation and Controls (GN&C) system-level models to represent propellant slosh forces on launch vehicles and spacecraft. The slosh model parameters are typically predicted using semi-empirical analytical methods that are based on experimental data for high liquid-to-gas density ratios (such as water and air at standard sea level conditions). Model parameter calculations typically neglect any contribution to slosh dynamics from the gas phase. However, some cryogenic propellant systems operate in conditions where the density ratio can be orders of magnitude smaller. Analytical and computational modeling was used in this effort to investigate the effects of liquid-to-gas density ratio on slosh dynamics. No experimental data nor previous studies are available at this time that quantify the effect of liquid-to-gas density ratio on slosh dynamics. However, in a recent CFD study, results showed that density ratio can have a significant effect on slosh model parameters. In this study, a dual pendulum slosh model was created that distinctly represents both liquid and gas phase dynamics. Solution of the model, which invokes the Euler-Lagrange equations of motion, demonstrates that slosh frequency is a combination of liquid slosh frequency and gas phase slosh frequency. Additionally, it shows that slosh mass is reduced due to the opposing motion of the gas phase with respect to that of the liquid phase. The slosh model parameter trends were verified using Computational Fluid Dynamics (CFD) analysis results for various liquid-to-gas density ratios, tanks, and fill levels.

Christopher D Moore↗

Numerical Model for Isothermal and Non-Isothermal Crystallization of Liquids and Glasses

A new numerical model of isothermal and non-isothermal first order phase transformations, such as the crystallization of liquids and glasses, is presented. This model computes directly the volume fraction transformed, taking into account time-dependent nucleation rates and cluster-size-dependent growth velocities. The model is applied to the crystallization of lithium disilicate glass, using the appropriate kinetic and thermodynamic parameters. The model is used (1) to determine the validity of common methods for computing the volume fraction transformed as a function of time in isothermal experiments when a time-dependent nucleation rate is expected, (2) to simulate non-isothermal differential scanning calorimetry (DSC) studies of glass devitrification as a function of scan rate, and (3) to compute the effect of preannealing on the DSC peak parameters. A novel behavior of the nucleation rate with scan rate is predicted, arising because the relaxation of the cluster distribution cannot be described by a single relaxation time. Comparisons of the calculations with experimental data on this glass demonstrate the validity of the model.

Kelton, K. F.↗

RL10A-3-3A Rocket Engine Modeling Project

Two RL10A-3-3A rocket engines comprise the main propulsion system for the Centaur upper stage vehicle. Centaur is used with bod Titan and Atlas launch vehicles, carrying military and civilian payloads from high altitudes into orbit and beyond. The RL10 has delivered highly reliable service for the past 30 years. Recently, however, there have been two in-flight failures which have refocused attention on the RL10. This heightened interest has sparked a desire for an independent RL10 modeling capability within NASA and th Air Force. Pratt & Whitney, which presently has the most detailed model of the RL10, also sees merit in having an independent model which could be used as a cross-check with their own simulations. The Space Propulsion Technology Division (SPTD) at the NASA Lewis Research Center has developed a computer model of the RL10A-3-3A. A project team was formed, consisting of experts in the areas of turbomachinery, combustion, and heat transfer. The overall goal of the project was to provide a model of the entire RL10 rocket engine for government use. In the course of the project, the major engine components have been modeled using a combination of simple correlations and detailed component analysis tools (computer codes). The results of these component analyses were verified with data provided by Pratt & Whitney. Select modeling results and test data curves were then integrated to form the RL10 engine system model The purpose of this report is to introduce the reader to the RL10 rocket engine and to describe the engine system model. The RL10 engine and its application to U.S. launch vehicles are described first, followed by a summary of the SPTD project organization, goals, and accomplishments. Simulated output from the system model are shown in comparison with test and flight data for start transient, steady state, and shut-down transient operations. Detailed descriptions of all component analyses, including those not selected for integration with the system model, are included as appendices.

Binder, Michael↗

Prediction of the structure of fuel sprays in gas turbine combustors

The structure of fuel sprays in a combustion chamber is theoretically investigated using computer models of current interest. Three representative spray models are considered: (1) a locally homogeneous flow (LHF) model, which assumes infinitely fast interphase transport rates; (2) a deterministic separated flow (DSF) model, which considers finite rates of interphase transport but ignores effects of droplet/turbulence interactions; and (3) a stochastic separated flow (SSF) model, which considers droplet/turbulence interactions using random sampling for turbulence properties in conjunction with random-walk computations for droplet motion and transport. Two flow conditions are studied to investigate the influence of swirl on droplet life histories and the effects of droplet/turbulence interactions on flow properties. Comparison of computed results with the experimental data show that general features of the flow structure can be predicted with reasonable accuracy using the two separated flow models. In contrast, the LHF model overpredicts the rate of development of the flow. While the SSF model provides better agreement with measurements than the DSF model, definitive evaluation of the significance of droplet/turbulence interaction is not achieved due to uncertainties in the spray initial conditions.

Shuen, J. S.↗

Prediction of the structure of fuel sprays in gas turbine combustors

The structure of fuel sprays in a combustion chamber is theoretically investigated using computer models of current interest. Three representative spray models are considered: (1) a locally homogeneous flow (LHF) model, which assumes infinitely fast interphase transport rates; (2) a deterministic separated flow (DSF) model, which considers finite rates of interphase transport but ignores effects of droplet/turbulence interactions; and (3) a stochastic separated flow (SSF) model, which considers droplet/turbulence interactions using random sampling for turbulence properties in conjunction with random-walk computations for droplet motion and transport. Two flow conditions are studied to investigate the influence of swirl on droplet life histories and the effects of droplet/turbulence interactions on flow properties. Comparison of computed results with the experimental data show that general features of the flow structure can be predicted with reasonable accuracy using the two separated flow models. In contrast, the LHF model overpredicts the rate of development of the flow. While the SSF model provides better agreement with measurements than the DSF model, definitive evaluation of the significance of droplet/turbulence interaction is not achieved due to uncertainties in the spray initial conditions.

Shuen, J.-S.↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Preconditioning electromyographic data for an upper extremity model using neural networks

A back propagation neural network has been employed to precondition the electromyographic signal (EMG) that drives a computational model of the human upper extremity. This model is used to determine the complex relationship between EMG and muscle activation, and generates an optimal muscle activation scheme that simulates the actual activation. While the experimental and model predicted results of the ballistic muscle movement are very similar, the activation function between the start and the finish is not. This neural network preconditions the signal in an attempt to more closely model the actual activation function over the entire course of the muscle movement.

Roberson, D. J.↗

Comparison of Spacecraft Contamination Models with Well-Defined Flight Experiment

During this reporting period, the phase 1 interim report was written and submitted. Detailed computer models were used to predict exposure to atomic oxygen and solar ultraviolet radiation. ESCA survey measurements were started on the tray wall section cut from LDEF tray location EIO-3. At this time there are no technical issues impeding the progress of the required tasks. Cost are in line withn the percentage of completion of tasks. The expended budget is 18% of the total and the project is approximately 18% complete. During the next reporting period we will continue to analyze points from tray locations E10-3 and A4-1 using ESCA and auger surface analysis techniques and a detailed location map will be produced.

Pippin, Gary↗

SPACE-Gateway: Modeling the Electrical Performance of the Gateway Power & Propulsion Element (PPE) Solar Arrays

The System Power Analysis for Capability Evaluation (SPACE) computer model was originally developed in 1988 at NASA Glenn Research Center to model the International Space Station (ISS) electrical power system (EPS). Since then, SPACE development has continued with the rise of new NASA human spaceflight missions. Today, SPACE’s uniquely complex and robust solar array model is being leveraged to model the solar arrays on the Gateway Power & Propulsion Element (PPE) for NASA's Artemis Program. This presentation will discuss the model’s integrated features, including its solar cell current-voltage (IV) curve, degradation, temperature, and on-orbit shadowing components. The presentation will also discuss recent SPACE advancements to meet the needs of NASA's latest spaceflight missions.

SPACE↗

Artemis 1 Radiation Modeling and Analysis Using Operational Methods

The successful completion of the Artemis 1 mission has provided NASA with a significant volume of space radiation dosimetry data that can be used to verify the accuracy of current operational SRAG dose computation models and methods for Beyond LEO missions. High-Fidelity CAD models of the Artemis 1 MPCV Orion were raytraced in 10,000 directions per sensor location, including the location of Artemis HERA sensors, Artemis RAMs (TLDs), and the crew active dosimetry (CADs). The full Artemis 1 mission was modelled using the as-flown trajectory, Ap9-Ae9 IRENE providing the VAB trapped proton and electron environment, and the Badhwar-O’Neill 2020 model providing the freespace GCR environment. IGRF 12 was used to compute cutoff rigidities for the GCR environment in LEO. 1DHZETRN was used to estimate point doses at each sensor position. Computed RAM doses were within 5% of measured values. The computed VAB transit and freespace GCR dose rates showed good qualitative agreement with measured HERA values. An estimate of normalized mean effective dose for a hypothetical male crew member on Artemis 1 using current standard operational methods was computed to be 22.3 mSv.

Artemis 1↗

Interactions of CO2, temperature and management practices: simulations with a modified version of CERES-Wheat

A new growth subroutine was developed for CERES-Wheat, a computer model of wheat (Triticum aestivum) growth and development. The new subroutine simulates canopy photosynthetic response to CO2 concentrations and light levels, and includes the effects of temperature on canopy light-use efficiency. Its performance was compared to the original CERES-Wheat V-2 10 in 30 different cases. Biomass and yield predictions of the two models were well correlated (correlation coefficient r > 0.95). As an application, summer growth of spring wheat was simulated at one site. Modeled crop responses to higher mean temperatures, different amounts of minimum and maximum warming, and doubled CO2 concentrations were compared to observations. The importance of irrigation and nitrogen fertilization in modulating the wheat crop climatic responses were also analyzed. Specifically, in agreement with observations, rainfed crops were found to be more sensitive to CO2 increases than irrigated ones. On the other hand, low nitrogen applications depressed the ability of the wheat crop to respond positively to CO2 increases. In general, the positive effects of high CO2 on grain yield were found to be almost completely counterbalanced by the negative effects of high temperatures. Depending on how temperature minima and maxima were increased, yield changes averaged across management practices ranged from -4% to 8%.

NASA Discipline Life Support Systems↗

Overview of High-Fidelity Modeling Activities in the Numerical Propulsion System Simulations (NPSS) Project

A high-fidelity simulation of a commercial turbofan engine has been created as part of the Numerical Propulsion System Simulation Project. The high-fidelity computer simulation utilizes computer models that were developed at NASA Glenn Research Center in cooperation with turbofan engine manufacturers. The average-passage (APNASA) Navier-Stokes based viscous flow computer code is used to simulate the 3D flow in the compressors and turbines of the advanced commercial turbofan engine. The 3D National Combustion Code (NCC) is used to simulate the flow and chemistry in the advanced aircraft combustor. The APNASA turbomachinery code and the NCC combustor code exchange boundary conditions at the interface planes at the combustor inlet and exit. This computer simulation technique can evaluate engine performance at steady operating conditions. The 3D flow models provide detailed knowledge of the airflow within the fan and compressor, the high and low pressure turbines, and the flow and chemistry within the combustor. The models simulate the performance of the engine at operating conditions that include sea level takeoff and the altitude cruise condition.

Veres, Joseph P.↗

Human transfer functions used to predict system performance parameters

Automatic, parameter-tracking, model-matching technique compares the responses of a human operator with those of an analog computer model of a human operator to predict and analyze the performance of mechanical or electromechanical systems prior to construction. Transfer functions represent the input-output relation of an operator controlling a closed-loop system.

Source record↗

Simulation and evaluation of latent heat thermal energy storage

The relative value of thermal energy storage (TES) for heat pump storage (heating and cooling) as a function of storage temperature, mode of storage (hotside or coldside), geographic locations, and utility time of use rate structures were derived. Computer models used to simulate the performance of a number of TES/heat pump configurations are described. The models are based on existing performance data of heat pump components, available building thermal load computational procedures, and generalized TES subsystem design. Life cycle costs computed for each site, configuration, and rate structure are discussed.

Sigmon, T. W.↗

Indium phosphide solar cell research in the US: Comparison with nonphotovoltaic sources

Highlights of the InP solar cell research program are presented. Homojunction cells with AMO efficiences approaching 19 percent were demonstrated while 17 percent was achieved for indium tin oxide (ITO)/InP cells. The superior radiation resistance of these latter two cell configurations over both Si and GaAs were demonstrated. InP cells on board the LIPS III satellite show no degradation after more than a year in orbit. Computer modeling calculations were directed toward radiation damage predictions and the specification of concentrator cell parameters. Computed array specific powers, for a specific orbit, are used to compare the performance of an InP solar cell array to solar dynamic and nuclear systems.

Weinberg, I.↗

Assessment of higher order turbulence models for complex two- and three-dimensional flowfields

A numerical method is presented to solve the three-dimensional Navier-Stokes equations in combination with a full Reynolds-stress turbulence model. Computations will be shown for three complex flowfields. The results of the Reynolds-stress model will be compared with those predicted by two different versions of the k-omega model. It will be shown that an improved version of the k-omega model gives as accurate results as the Reynolds-stress model.

Menter, Florian R.↗

Bubble Formation on a Wall in Cross-Flowing Liquid and Surrounding Fluid Motion,With and Without Heating

Application of gas-liquid two-phase flow systems for space-based thermal management and for the HEDS program demands a precise control of bubble size distribution in liquid. The necessity of bulk liquid motion for controlling bubble size and frequency in the space environment has been suggested by recent studies on pool, forced convection boiling and bubble formation in flowing liquid. The present work, consisting of two parts, explores bubble generation at wall in a cross-flowing liquid, i.e., in a forced convection boiling configuration. A schematic is shown. The first part looks into the bubble formation process under isothermal conditions in a reduced gravity environment, by injecting gas through a hole in the wall of a flowing liquid channel. In the latter part with channel wall heating, flow and temperature fields near a single bubble are studied under normal (1-g) and micro-gravity (mu-g) conditions. The objective of the isothermal experiments is to experimentally investigate the effects of liquid cross-flow velocity, gas flow rate, and orifice diameter on bubble formation. Data were taken mainly under reduced gravity conditions but some data were taken in normal gravity for comparison. The reduced gravity experiment was conducted aboard the NASA DC-9 Reduced Gravity Aircraft. The results show that the process of bubble formation and detachment depends on gravity, the orifice diameter (D(sub N)), the gas flow rate (Q(sub g)), and the liquid cross-flow velocity (U(sub L)). The reduced gravity data are shown. The data are analyzed based on a force balance, and two different detachment mechanisms are identified. When the gas momentum is large, the bubble detaches from the injection orifice as the gas momentum overcomes the attaching effects of liquid drag and inertia. The surface tension force is much reduced because a large part of the bubble pinning edge at the orifice is lost as the bubble axis is tilted by the liquid flow. When the gas momentum is small, the force balance in the liquid flow direction is important, and the bubble detaches when the bubble axis inclination exceeds a certain angle. With wall heating, liquid motion around an air bubble in cross-flowing 2cs silicone oil is experimentally investigated in 1-g. A spectral element based steady 2D numerical model is also developed. The traces of particles from experimental flow visualization and the corresponding computed streamlines are shown. At the upstream side of the bubble facing the cross-flow, thermocapillary and forced convection create liquid motion away from the wall, up along the surface. At the downstream side, a competing interaction between the two creates a recirculation cell, causing the bulk liquid to stagnate on the surface and separate thereafter. The important dimensionless parameters are - Surface tension and local cross-flow Reynolds numbers R(sub sigma)(U(sub ref)a/v) and Re(sub loc)(U(sub L)a/v), respectively based on reference thermocapillary U(sub ref)(sigma(sub T)Delta.T/mu, Delta T=T(sub wall)-T(sub liquid)) and local cross-flow velocity U(sub L), Prandtl number P(sub r) and Grashoff number Gr(rho.g.beta.DeltaTa(sup 3)/mu.v). Variation of the stagnation point with R(sub sigma) and Re(sub loc) is shown. Also shown are good agreement between experimental and numerical results in 1-g. The computational model is extended to mu-g condition to investigate temperature and velocity on the bubble surface, stagnation and reattachment points of the recirculation cell and wall heat transfer. It is observed that wall heating significantly alters the flow field around the bubble and thus the forces acting on the bubble, which control its detachment. Thus a combination of heating and liquid cross-flow can be utilized to precisely control bubble formation in a mu-g environment.

Bhunia, Avijit↗