Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Modeling”

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

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

At least 1,675 records · Page 93

Dynamic responses of railroad car models to vertical and lateral rail inputs

Simplified dynamic models were applied in a study of vibration in a high-speed railroad car. The mathematical models used were a four-degree-of-freedom model for vertical responses to vertical rail inputs and a ten-degree-of-freedom model for lateral response to lateral or rolling (cross-level) inputs from the rails. Elastic properties of the passenger car body were represented by bending and torsion of a uniform beam. Rail-to-car (truck) suspensions were modeled as spring-mass-dashpot oscillators. Lateral spring nonlinearities approximating certain complicated truck mechanisms were introduced. The models were excited by displacement and, in some cases, velocity inputs from the rails by both deterministic (including sinusoidal) and random input functions. Results were obtained both in the frequency and time domains. Solutions in the time domain for the lateral model were obtained for a wide variety of transient and random inputs generated on-line by an analog computer. Variations in one of the damping properties of the lateral car suspension gave large fluctuations in response over a range of car speeds for a given input. This damping coefficient was significant in reducing lateral car responses that were higher for nonlinear springs for three different inputs.

Sewall, J. L.↗

Time dependent deformation and stress in the lithosphere

Efficient computer programs incorporating frontal solution and time stepping procedure were developed for the modelling of geodynamic problems. This scheme allows for investigating the quasi static phenomena including the effects of the rheological structure of a tectonically active region. From three dimensional models of strike slip earthquakes, it was found that lateral variation of viscosity affects the characteristics of surface deformations. The vertical deformation is especially informative about the viscosity structure in a strike slip fault zone. A three dimensional viscoelastic model of a thrust earthquake indicated that the transient disturbance on plate velocity due to a great plate boundary earthquake is significant at intermediate distances, but becomes barely measurable 1000 km away from the source.

Yang, M.↗

Spacecraft loads methodology - Transient vs. shock spectra method

The methodology for the establishment of spacecraft loads is strongly influenced by project constraints which include the cost, schedule and allowable weight. The most rigorous approach is the transient loads analysis which requires a composite mathematical model of the spacecraft and launch vehicle. The structural member loads for the entire composite structure are computed by applying the forcing functions, which represent various dynamic environments during the mission, to the composite model. Although this method ideally leads to a lightweight design, it is costly and time consuming due to complex interfaces involving many organizations. To reduce complexity and cost a shock spectra method has been used to design spacecraft structures. This method utilizes envelopes of shock spectra of launch vehicle accelerations obtained from analysis and/or flight measurements. Since only limited information on the launch vehicle model is involved in this process the design loads iteration cycle can be rapidly performed within the payload organization. In the present paper, these two methods will be evaluated by comparing the loads for several spacecraft. Flight measured loads will also be used in the evaluation.

Chen, J. C.↗

Molecular clouds and galactic spiral structure

Galactic CO line emission at 115 GHz was surveyed in order to study the distribution of molecular clouds in the inner galaxy. Comparison of this survey with similar H1 data reveals a detailed correlation with the most intense 21 cm features. To each of the classical 21 cm H1 spiral arms of the inner galaxy there corresponds a CO molecular arm which is generally more clearly defined and of higher contrast. A simple model is devised for the galactic distribution of molecular clouds. The modeling results suggest that molecular clouds are essentially transient objects, existing for 15 to 40 million years after their formation in a spiral arm, and are largely confined to spiral features about 300 pc wide.

Dame, T. M.↗

Dynamic characteristics of two 300 kW class dual keel space station concepts

Results from an investigation of the dynamic behavior of a 300 kW class solar dynamic powered, dual kell space station are presented. The purpose of the investigation was to determine and assess the influence of space station truss bay size on station controllability during rigid body attitude adjustment and orbit reboost maneuvers. The dual keel space station concept is defined and two finite element models (one which has a truss bay size of 5m and another with a truss bay size of 9 feet) are described. Rigid and flexible body characteristics of the two space station models are also presented. Finally, results from a transient response analysis, where the stations are subjected to an orbit reboost maneuver, are summarized.

Dorsey, J. T.↗

User's manual for a computer program for the emulation/simulation of a space station Environmental Control and Life Support System (ESCM)

This manual describes how to use the Emulation Simulation Computer Model (ESCM). Based on G189A, ESCM computes the transient performance of a Space Station atmospheric revitalization subsystem (ARS) with CO2 removal provided by a solid amine water desorbed subsystem called SAWD. Many performance parameters are computed some of which are cabin CO2 partial pressure, relative humidity, temperature, O2 partial pressure, and dew point. The program allows the user to simulate various possible combinations of man loading, metabolic profiles, cabin volumes and certain hypothesized failures that could occur.

Yanosy, James L.↗

Ocean alkalinity and the Cretaceous/Tertiary boundary

A biogeochemical cycle model resolving ocean carbon and alkalinity content is applied to the Maestrichtian and Danian. The model computes oceanic concentrations and distributions of Ca(2+), Mg(2+), and Sigma-CO2. From these values an atmospheric pCO2 value is calculated, which is used to estimate rates of terrestrial weathering of calcite, dolomite, and calcium and magnesium silicates. Metamorphism of carbonate rocks and the subsequent outgassing of CO2 to the atmosphere are parameterized in terms of carbonate rock reservoir sizes, total land area, and a measure of overall tectonic activity, the sea-floor generation rate. The ocean carbon reservoir computed by the model is used with Deep Sea Drilling Project (DSDP) C-13 data to estimate organic detrital fluxes under a variety of ocean mixing rate assumptions. Using Redfield ratios, the biogenic detrital flux estimate is used to partition the ocean carbon and alkalinity reservoirs between the mixed layer and deep ocean. The calcite flux estimate and carbonate ion concentrations are used to determine the rate of biologically mediated CaCO3 titration. Oceanic productivity was severely limited for approximately 500 kyr following the K/T boundary resulting in significant increases in total ocean alkalinity. As productivity returned to the ocean, excess carbon and alkalinity was removed from the ocean as CaCO3. Model runs indicate that this resulted in a transient imbalance in the other direction. Ocean chemistry returned to near-equilibrium by about 64 mybp.

Caldeira, K. G.↗

Ionospheric chemical releases

Ionospheric plasma density irregularities can be produced by chemical releases into the upper atmosphere. F-region plasma modification occurs by: (1) chemically enhancing the electron number density; (2) chemically reducing the electron population; or (3) physically convecting the plasma from one region to another. The three processes (production, loss, and transport) determine the effectiveness of ionospheric chemical releases in subtle and surprising ways. Initially, a chemical release produces a localized change in plasma density. Subsequent processes, however, can lead to enhanced transport in chemically modified regions. Ionospheric modifications by chemical releases excites artificial enhancements in airglow intensities by exothermic chemical reactions between the newly created plasma species. Numerical models were developed to describe the creation and evolution of large scale density irregularities and airglow clouds generated by artificial means. Experimental data compares favorably with theses models. It was found that chemical releases produce transient, large amplitude perturbations in electron density which can evolve into fine scale irregularities via nonlinear transport properties.

Bernhardt, Paul A.↗

Some further observations on droplet combustion characteristics - NASA LeRC-Princeton results

Experimental and numerical studies are reviewed which are designed to examine the effects of droplet/gas motion, product-intermediate absorption, extinction, and sooting on droplet combustion. The experimental work at the NASA-Lewis Research Center involves a 2.2-s droptower for investigating microgravitational effects of droplet combustion over a relatively extended range. The droplet-gas velocities are very low because the spherosymmetrical nature of the major combustion processes produces a quiescent environment. The refined experimental results are combined with numerical modeling based on a technique that is fully transient, comprehensive, and has few empirical simplifications. The combination of techniques improves the present understanding of convection-induced effects, reducing soot formation, and promoting quiescent droplet combustion.

Choi, Mun Y.↗

Response of laminated composite plates to low-speed impact by airgun-propelled and dropped-weight impactors

An analytical procedure has been developed to determine the transient response of simply supported, rectangular laminated composite plates subjected to impact loads from airgun-propelled or dropped-weight impactors. A first-order shear-deformation theory has been included in the analysis to represent properly any local short-wavelength transient bending response. The impact force has been modeled as a locally distributed load with a cosine-cosine distribution. A double Fourier series expansion and the Timoshenko small increment method have been used to determine the contact force, out-of-plane deflections, and inplane strains and stresses at any plate location due to an impact force at any plate location. The results of experimental and analytical studies are compared for quasi-isotropic laminates. The results indicate the importance of including transverse shear deformation effects in the analysis for predicting the response of laminated plates subjected to both airgun-propelled and dropped-weight impactors. The results also indicate that plate boundary conditions influence the axial strains more significantly than the contact force for a dropped-weight impactor. The results of parametric studies identify a scaling approach based on impactor momentum that suggests an explanation for the differences in the responses of plates impacted by airgun-propelled or dropped-weight impactors.

Prasad, Chunchu B.↗

Response of laminated composite plates to low-speed impact by different impactors

An analytic procedure has been developed to determine the transient response of simply supported, retangular laminated composite plates subjected to impact loads from airgun-propelled or drop-weight impactors. A first-order shear deformation theory has been included in the analysis to represent properly any local short-wavelength transient bending response. The impact force has been modeled as a locally distributed load with cosine-cosine distribution. A double Fourier series expansion and the Timoshenko small increment method have been used to determine the contact force, out-of-plane deflections, and in-plane strains and stresses at any plate location due to an impact force at any plate location. The results of experimental and analytical studies are compared for quasi-isotropic laminates. The results indicate the importance of including transverse shear deformation effects in the analysis for predicting the response of laminated plates subjected to both airgun-propelled and dropped-weight impactors. The results also indicate that plate boundary conditions influence the axial strains more significantly than the contact force for a dropped-weight impactor. The results of parametric studies identify a scaling approach based on impactor momentum that may account for the differences in the responses of plates impacted by airgun-propelled or dropped-weight impactors.

Prasad, Chunchu↗

Interannual Variability of Martian Global Dust Storms: Simulations with a Low-Order Model of the General Circulation

We present simulations of the interannual variability of martian global dust storms (GDSs) with a simplified low-order model (LOM) of the general circulation. The simplified model allows one to conduct computationally fast long-term simulations of the martian climate system. The LOM is constructed by Galerkin projection of a 2D (zonally averaged) general circulation model (GCM) onto a truncated set of basis functions. The resulting LOM consists of 12 coupled nonlinear ordinary differential equations describing atmospheric dynamics and dust transport within the Hadley cell. The forcing of the model is described by simplified physics based on Newtonian cooling and Rayleigh friction. The atmosphere and surface are coupled: atmospheric heating depends on the dustiness of the atmosphere, and the surface dust source depends on the strength of the atmospheric winds. Parameters of the model are tuned to fit the output of the NASA AMES GCM and the fit is generally very good. Interannual variability of GDSs is possible in the IBM, but only when stochastic forcing is added to the model. The stochastic forcing could be provided by transient weather systems or some surface process such as redistribution of the sand particles in storm generating zones on the surface. The results are sensitive to the value of the saltation threshold, which hints at a possible feedback between saltation threshold and dust storm activity. According to this hypothesis, erodable material builds up its a result of a local process, whose effect is to lower the saltation threshold until a GDS occurs. The saltation threshold adjusts its value so that dust storms are barely able to occur.

Pankine, A. A.↗

Heat Transfer in High-Temperature Fibrous Insulation

The combined radiation/conduction heat transfer in high-porosity, high-temperature fibrous insulations was investigated experimentally and numerically. The effective thermal conductivity of fibrous insulation samples was measured over the temperature range of 300-1300 K and environmental pressure range of 1.33 x 10(exp -5)-101.32 kPa. The fibrous insulation samples tested had nominal densities of 24, 48, and 72 kilograms per cubic meter and thicknesses of 13.3, 26.6 and 39.9 millimeters. Seven samples were tested such that the applied heat flux vector was aligned with local gravity vector to eliminate natural convection as a mode of heat transfer. Two samples were tested with reverse orientation to investigate natural convection effects. It was determined that for the fibrous insulation densities and thicknesses investigated no heat transfer takes place through natural convection. A finite volume numerical model was developed to solve the governing combined radiation and conduction heat transfer equations. Various methods of modeling the gas/solid conduction interaction in fibrous insulations were investigated. The radiation heat transfer was modeled using the modified two-flux approximation assuming anisotropic scattering and gray medium. A genetic-algorithm based parameter estimation technique was utilized with this model to determine the relevant radiative properties of the fibrous insulation over the temperature range of 300-1300 K. The parameter estimation was performed by least square minimization of the difference between measured and predicted values of effective thermal conductivity at a density of 24 kilograms per cubic meters and at nominal pressures of 1.33 x 10(exp -4) and 99.98 kPa. The numerical model was validated by comparison with steady-state effective thermal conductivity measurements at other densities and pressures. The numerical model was also validated by comparison with a transient thermal test simulating reentry aerodynamic heating conditions.

Daryabeigi, Kamran↗

Transient Two-Dimensional Analysis of Side Load in Liquid Rocket Engine Nozzles

Two-dimensional planar and axisymmetric numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to develop a computational methodology to identify nozzle side load physics using simplified two-dimensional geometries, in order to come up with a computational strategy to eventually predict the three-dimensional side loads. The computational methodology is based on a multidimensional, finite-volume, viscous, chemically reacting, unstructured-grid, and pressure-based computational fluid dynamics formulation, and a transient inlet condition based on an engine system modeling. The side load physics captured in the low aspect-ratio, two-dimensional planar nozzle include the Coanda effect, afterburning wave, and the associated lip free-shock oscillation. Results of parametric studies indicate that equivalence ratio, combustion and ramp rate affect the side load physics. The side load physics inferred in the high aspect-ratio, axisymmetric nozzle study include the afterburning wave; transition from free-shock to restricted-shock separation, reverting back to free-shock separation, and transforming to restricted-shock separation again; and lip restricted-shock oscillation. The Mach disk loci and wall pressure history studies reconfirm that combustion and the associated thermodynamic properties affect the formation and duration of the asymmetric flow.

Wang, Ten-See↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A High-Fidelity Simulation of a Generic Commercial Aircraft Engine and Controller

A new high-fidelity simulation of a generic 40,000 lb thrust class commercial turbofan engine with a representative controller, known as CMAPSS40k, has been developed. Based on dynamic flight test data of a highly instrumented engine and previous engine simulations developed at NASA Glenn Research Center, this non-proprietary simulation was created especially for use in the development of new engine control strategies. C-MAPSS40k is a highly detailed, component-level engine model written in MATLAB/Simulink (The MathWorks, Inc.). Because the model is built in Simulink, users have the ability to use any of the MATLAB tools for analysis and control system design. The engine components are modeled in C-code, which is then compiled to allow faster-than-real-time execution. The engine controller is based on common industry architecture and techniques to produce realistic closed-loop transient responses while ensuring that no safety or operability limits are violated. A significant feature not found in other non-proprietary models is the inclusion of transient stall margin debits. These debits provide an accurate accounting of the compressor surge margin, which is critical in the design of an engine controller. This paper discusses the development, characteristics, and capabilities of the C-MAPSS40k simulation

May, Ryan D.↗

Dynamic Altitude Simulation System Performance Modeling

The "Dynamic Altitude Simulation Prediction Program," DASSPP, is a program to predict the transient response of an engine, test cell, ejector system under engine shutdown conditions. These transients are important to know so that corrective modifications can be adapted to prevent any damage to the engine or test cell. The "Dynamic Altitude System Simulation Prediction Program," DASSPP, is a major rewrite of the existing program "RL-1000" written in BASIC. The RL-1000 program was written to analyze the transients of the RL-10 system only. The new program is written to run in Excell 97 and utilizes the Visual BASIC language in Excell. The program has many added features not included in the original "Rl-1000" program. The program utilizes the ejector models developed during the summer of 1997. The new program is very user friendly and utilizes a dialog box for data input.

LaFrance, Leo J.↗