Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modelling Capabilities”

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 721 records · Page 40

Development of a unified constitutive model for an isotropic nickel base superalloy Rene 80

Accurate analysis of stress-strain behavior is of critical importance in the evaluation of life capabilities of hot section turbine engine components such as turbine blades and vanes. The constitutive equations used in the finite element analysis of such components must be capable of modeling a variety of complex behavior exhibited at high temperatures by cast superalloys. The classical separation of plasticity and creep employed in most of the finite element codes in use today is known to be deficient in modeling elevated temperature time dependent phenomena. Rate dependent, unified constitutive theories can overcome many of these difficulties. A new unified constitutive theory was developed to model the high temperature, time dependent behavior of Rene' 80 which is a cast turbine blade and vane nickel base superalloy. Considerations in model development included the cyclic softening behavior of Rene' 80, rate independence at lower temperatures and the development of a new model for static recovery.

Ramaswamy, V. G.↗

QBO Influence on Polar Stratospheric Variability in the GEOS Chemistry-Climate Model

The quasi-biennial oscillation modulates the strength of both the Arctic and Antarctic stratospheric vortices. Model and observational studies have found that the phase and characteristics of the quasi-biennial oscillation (QBO) contribute to the high degree of variability in the Arctic stratosphere in winter. While the Antarctic stratosphere is less variable, recent work has shown that Southern Hemisphere planetary wave driving increases in response to "warm pool" El Nino events that are coincident with the easterly phase of the QBO. These events hasten the breakup of the Antarctic polar vortex. The Goddard Earth Observing System (GEOS) chemistry-climate model (CCM) is now capable of generating a realistic QBO, due a new parameterization of gravity wave drag. In this presentation, we will use this new model capability to assess the influence of the QBO on polar stratospheric variability. Using simulations of the recent past, we will compare the modeled relationship between QBO phase and mid-winter vortex strength with the observed Holton-Tan relation, in both hemispheres. We will use simulations of the 21 St century to estimate future trends in the relationship between QBO phase and vortex strength. In addition, we will evaluate the combined influence of the QBO and El Nino/Southern Oscillation (ENSO) on the timing of the breakup of the polar stratospheric vortices in the GEOS CCM. We will compare the influence of these two natural phenomena with trends in the vortex breakup associated with ozone recovery and increasing greenhouse gas concentrations.

Hurwitz, M. M.↗

Investigation of Thermospheric and Ionospheric Changes during Ionospheric Storms with Satellite and Ground-Based Data and Modeling

The purpose of this proposed research is to improve our basic understanding of the causes of ionospheric storm behavior in the midlatitude F region ionosphere. This objective will be achieved by detailed comparisons between ground based measurements of the peak electron density (N(sub m)F(sub 2)), Atmosphere Explorer satellite measurements of ion and neutral composition, and output from the Field Line Interhemispheric Plasma (FLIP) model. The primary result will be a better understanding of changes in the neutral densities and ion chemistry during magnetic storms that will improve our capability to model the weather of the ionosphere which will be needed as a basis for ionospheric prediction. Specifically, this study seeks to answer the following questions: (1) To what extent are negative ionospheric storm phases caused by changes in the atomic to molecular ratio? (2) Are the changes in neutral density ratio due to increased N2, or decreased O, or both? (3) Are there other chemical processes (e.g., excited N2) that increase O+ loss rates during negative storms? (4) Do neutral density altitude distributions differed from hydrostatic equilibrium? (5) Why do near normal nighttime densities often follow daytime depletions of electron density; and (6) Can changes in h(sub m)F2 fully account for positive storm phases? To answer these questions, we plan to combine ground-based and space-based measurements with the aid of our ionospheric model which is ideally suited to this purpose. These proposed studies will lead to a better capability to predict long term ionospheric variability, leading to better predictions of ionospheric weather.

Richards, Philip G.↗

Longitudinal splitting in unidirectional composites, analysis and experiments

An experimental study is conducted to determine the fracture behavior of center notched, unidirectional graphite/epoxy laminates when subjected to tensile loading. The actual behavior is compared to the behavior predicted by a mathematical model based on classical shear-lag assumptions. The model allows for damage to occur in the form of longitudinal matrix yielding and splitting with the matrix assumed to fail in pure shear. Results indicate that the model is capable of predicting split initiation stress levels accurately, but does not describe the subsequent split growth adequately. The model predicts rapid split growth following split initiation due to shear failure, while the actual behavior involves a slow split growth region prior to the rapid growth region. It is suggested that transverse matrix normal stresses are responsible for split initiation and the early, slow split growth. The model predicts the actual initiation stress levels reliably, and also appears to be able to predict the point at which the shear failure mode begins to dominate. The shear failure mode does eventually dominate, but at a slower rate than predicted. The nonuniform structure of the graphite/epoxy laminates is thought to be responsible for decreasing the split growth rate due to shear failure.

Goree, J. G.↗

Understanding the Atmosphere of 51 Eri b: Do Photochemical Hazes Cloud the Planets Spectrum?

The first young giant planet to be discovered by the Gemini Planet Imager was the (is) approximately 2MJ planet 51 Eri b. This approximately 20 Myr old young Jupiter is the first directly imaged planet to show unmistakable methane in H band. To constrain the planet's mass, atmospheric temperature, and composition, the GPI J and H band spectra as well as some limited photometric points were compared to the predictions of substellar atmosphere models. The best fitting models reported in the discovery paper (Macintosh et al. 2015) relied upon a combination of clear and cloudy atmospheric columns to reproduce the data. However for an object as cool as 700 K, the origin of the cloud coverage is somewhat puzzling, as the global silicate and iron clouds would be expected to have sunk well below the photosphere by this effective temperature. While strong vertical mixing in these low gravity atmospheres remains a plausible explanation, we have explored whether atmospheric photochemistry, driven by the UV flux from the primary star, may yield hazes that also influence the observed spectrum of the planet. To explore this possibility we have modeled the atmospheric photochemistry of 51 Eri b using two state-of-the-art photochemical models, both capable of predicting yields of complex hydrocarbons under various atmospheric conditions. In our presentation we will summarize the modeling approach employed to characterize 51 Eri b, explaining constraints on the planet's effective temperature, gravity, and atmospheric composition and also present results of our studies of atmospheric photochemistry. We will discuss whether photochemical hazes could indeed be responsible for the particulate opacity that apparently sculpts the spectrum of the planet.

GPI - Gemini Planet Imager↗

Thermal systems analysis for the Space Infrared Telescope Facility dewar

Thermal systems analysis models were used to design SFHe cooled dewar for the Space Infrared Telescope Facility (SIRTF), a 1 m class cryogenically cooled observatory for IR astronomy. The models are capable of computing both the heat leaks into the dewar and the operating temperature of a SFHe tank. The models are aimed at predicting the ability of the SIRTF cryogenic system to satisfy a five-year mission lifetime requirement and maintain the SFHe tank operating temperature of 1.25 K to provide sufficient cooling for science instruments and the optical system. The thermal models are very detailed and very fast with a typical steady state run of about 20 sec on a VAX minicomputer.

Bhandari, Pradeep↗

Theoretical Tools and Software for Modeling, Simulation and Control Design of Rocket Test Facilities

A rocket test stand and associated subsystems are complex devices whose operation requires that certain preparatory calculations be carried out before a test. In addition, real-time control calculations must be performed during the test, and further calculations are carried out after a test is completed. The latter may be required in order to evaluate if a particular test conformed to specifications. These calculations are used to set valve positions, pressure setpoints, control gains and other operating parameters so that a desired system behavior is obtained and the test can be successfully carried out. Currently, calculations are made in an ad-hoc fashion and involve trial-and-error procedures that may involve activating the system with the sole purpose of finding the correct parameter settings. The goals of this project are to develop mathematical models, control methodologies and associated simulation environments to provide a systematic and comprehensive prediction and real-time control capability. The models and controller designs are expected to be useful in two respects: 1) As a design tool, a model is the only way to determine the effects of design choices without building a prototype, which is, in the context of rocket test stands, impracticable; 2) As a prediction and tuning tool, a good model allows to set system parameters off-line, so that the expected system response conforms to specifications. This includes the setting of physical parameters, such as valve positions, and the configuration and tuning of any feedback controllers in the loop.

Richter, Hanz↗

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Majumdar, Alok↗

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗

Computational simulation of liquid fuel rocket injectors

A major component of any liquid propellant rocket is the propellant injection system. Issues of interest include the degree of liquid vaporization and its impact on the combustion process, the pressure and temperature fields in the combustion chamber, and the cooling of the injector face and chamber walls. The Finite Difference Navier-Stokes (FDNS) code is a primary computational tool used in the MSFC Computational Fluid Dynamics Branch. The branch has dedicated a significant amount of resources to development of this code for prediction of both liquid and solid fuel rocket performance. The FDNS code is currently being upgraded to include the capability to model liquid/gas multi-phase flows for fuel injection simulation. An important aspect of this effort is benchmarking the code capabilities to predict existing experimental injection data. The objective of this MSFC/ASEE Summer Faculty Fellowship term was to evaluate the capabilities of the modified FDNS code to predict flow fields with liquid injection. Comparisons were made between code predictions and existing experimental data. A significant portion of the effort included a search for appropriate validation data. Also, code simulation deficiencies were identified.

Landrum, D. Brian↗

DSN Simulator

The DSN Simulator (wherein DSN signifies NASA's Deep Space Network) is an updated version of the software described in DSN Array Simulator (NPO-44506), Software Tech Briefs (Special supplement to NASA Tech Briefs), Vol. 32, No. 9 (September 2008), page 26. To recapitulate: This software is used for computational modeling of proposed DSN facilities comprising arrays of antennas and transmitting and receiving equipment for microwave communication with spacecraft on interplanetary missions. Such modeling is performed to estimate facility performance, evaluate requirements that govern facility design, and evaluate proposed improvements in hardware and/or software. The software includes a Monte Carlo simulation component that enables rapid generation of key mission-set metrics (e.g., numbers of links, data rates, and data volumes), and statistical distributions thereof as functions of time. The prior version of the software could model only one DSN facility at a time and included hard-coded, unconfigurable metrics. The present updated version is capable of modeling the entire DSN and provides for configurable metrics, making it possible to perform loading analyses for alternative future DSN architectures and mission-set scenarios. The present version also features an improved user interface and interfaces for exchange of data with other DSN software and with a DSN mission model database.

Fijany, Amir↗

A Dual-Plane PIV Study of Turbulent Heat Transfer Flows

Thin film cooling is a widely used technique in turbomachinery and rocket propulsion applications, where cool injection air protects a surface from hot combustion gases. The injected air typically has a different velocity and temperature from the free stream combustion flow, yielding a flow field with high turbulence and large temperature differences. These thin film cooling flows provide a good test case for evaluating computational model prediction capabilities. The goal of this work is to provide a database of flow field measurements for validating computational flow prediction models applied to turbulent heat transfer flows. In this work we describe the application of a Dual-Plane Particle Image Velocimetry (PIV) technique in a thin film cooling wind tunnel facility where the injection air stream velocity and temperatures are varied in order to provide benchmark turbulent heat transfer flow field measurements. The Dual-Plane PIV data collected include all three components of velocity and all three components of vorticity, spanning the width of the tunnel at multiple axial measurement planes.

Turbulent heat transfer↗

Coupling a Computational Fluid Dynamics Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop® (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

Heat transfer↗

Coupling a Computational Fluid Dynamics (CFD) Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

heat transfer↗

Coupled electromechanical response of composite beams with embedded piezoelectric sensors and actuators

Unified mechanics are developed with the capability to model both sensory and active composite laminates with embedded piezoelectric layers. A discrete-layer formulation enables analysis of both global and local electromechanical response. The mechanics include the contributions from elastic, piezoelectric, and dielectric components. The incorporation of electric potential into the state variables permits representation of general electromechanical boundary conditions. Approximate finite element solutions for the static and free-vibration analysis of beams are presented. Applications on composite beams demonstrate the capability to represent either sensory or active structures and to model the complicated stress-strain fields, the interactions between passive/active layers, interfacial phenomena between sensors and composite plies, and critical damage modes in the material. The capability to predict the dynamic characteristics under various electrical boundary conditions is also demonstrated.

Saravanos, D. A.↗

Modelling of thick composites using a layerwise laminate theory

The layerwise laminate theory of Reddy (1987) is used to develop a layerwise, two-dimensional, displacement-based, finite element model of laminated composite plates that assumes a piecewise continuous distribution of the tranverse strains through the laminate thickness. The resulting layerwise finite element model is capable of computing interlaminar stresses and other localized effects with the same level of accuracy as a conventional 3D finite element model. Although the total number of degrees of freedom are comparable in both models, the layerwise model maintains a 2D-type data structure that provides several advantages over a conventional 3D finite element model, e.g. simplified input data, ease of mesh alteration, and faster element stiffness matrix formulation. Two sample problems are provided to illustrate the accuracy of the present model in computing interlaminar stresses for laminates in bending and extension.

Robbins, D. H., Jr.↗

Quasi-biennial oscillations of ozone and diabatic circulation in the equatorial stratosphere

The quasi-biennial oscillation (QBO) in ozone in the equatorial stratosphere is obtained by analyzing the Stratospheric Aerosol and Gas Experiment (SAGE) data from 1984 to 1989. The phase of the ozone QBO in the lower stratosphere is found to precede the zonal wind QBO by several months as opposed to the theoretically expected in-phase relationship between the two. A mechanistic model is developed to explore possible reasons for this disagreement. The model is capable of simulating the actual time evolution of the ozone QBO by introducing the observed zonal wind profile as input. The modeled results confirm the conventional view that the ozone QBO is generated by the vertical ozone advection that is driven to maintain the temperature structure against radiative damping. However, a series of experiments emphasizes the importance of the feedback of the ozone QBO to the diabatic heating through the absorption of solar radiation. Due to this effect, the phase of the ozone QBO shifts up to a quarter cycle ahead and approaches that of the temperature QBO. Because of this inphase relationship, the feedback of the ozone QBO to the diabatic heating acts to compensate for the radiative damping of the temperature structure, thus reducing the magnitude of the induced diabatic circulation. Because the reduction of the magnitude of the vertical motion facilitates downward transport of easterly momentum by the mean flow, this feedback process can help to resolve the insufficiency of the easterly momentum in driving the dynamical QBO in general circulation models (GCMs). It should be emphasized that more sophisticated models that allow for full interaction between the chemical species and radiative and dynamical processes should be developed to improve our understanding of both dynamical and ozone QBOs.

Hasebe, Fumio↗

Modeling of LEO Orbital Debris Populations in Centimeter and Millimeter Size Regimes

The building of the NASA Orbital Debris Engineering Model, whether ORDEM2000 or its recently updated version ORDEM2010, uses as its foundation a number of model debris populations, each truncated at a minimum object-size ranging from 10 micron to 1 m. This paper discusses the development of the ORDEM2010 model debris populations in LEO (low Earth orbit), focusing on centimeter (smaller than 10 cm) and millimeter size regimes. Primary data sets used in the statistical derivation of the cm- and mm-size model populations are from the Haystack radar operated in a staring mode. Unlike cataloged objects of sizes greater than approximately 10 cm, ground-based radars monitor smaller-size debris only in a statistical manner instead of tracking every piece. The mono-static Haystack radar can detect debris as small as approximately 5 mm at moderate LEO altitudes. Estimation of millimeter debris populations (for objects smaller than approximately 6 mm) rests largely on Goldstone radar measurements. The bi-static Goldstone radar can detect 2- to 3-mm objects. The modeling of the cm- and mm-debris populations follows the general approach to developing other ORDEM2010-required model populations for various components and types of debris. It relies on appropriate reference populations to provide necessary prior information on the orbital structures and other important characteristics of the debris objects. NASA's LEO-to-GEO Environment Debris (LEGEND) model is capable of furnishing such reference populations in the desired size range. A Bayesian statistical inference process, commonly adopted in ORDEM2010 model-population derivations, changes a priori distribution into a posteriori distribution and thus refines the reference populations in terms of data. This paper describes key elements and major steps in the statistical derivations of the cm- and mm-size debris populations and presents results. Due to lack of data for near 1-mm sizes, the model populations of 1- to 3.16-mm objects are an empirical extension from larger debris. The extension takes into account the results of micro-debris (from 10 micron to 1 mm) population modeling that is based on shuttle impact data, in the hope of making a smooth transition between micron and millimeter size regimes. This paper also includes a brief discussion on issues and potential future work concerning the analysis and interpretation of Goldstone radar data.

Xu, Y.-L.↗