Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

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 613 records · Page 34

Implementation of a Blowing Boundary Condition in the LAURA Code

Preliminary steps toward modeling a coupled ablation problem using a finite-volume Navier-Stokes code (LAURA) are presented in this paper. Implementation of a surface boundary condition with mass transfer (blowing) is described followed by verification and validation through comparisons with analytic results and experimental data. Application of the code to a carbon-nosetip ablation problem is demonstrated and the results are compared with previously published data. It is concluded that the code and coupled procedure are suitable to support further ablation analyses and studies.

Thompson, Richard a.↗

Initial comparison of SSPM ground test results and flight data to NASCAP simulations

The satellite surface potential monitor (SSPM) has been developed for the P78-2 SCATHA (Spacecraft Charging At The High Altitudes) satellite to determine the response of selected spacecraft materials to charged-particle environmental fluxes. Since the monitor infers total surface voltages from a single point on the interior side of insulators, a ground simulation program was undertaken to develop analytical techniques to model the monitors and to obtain an experimental calibration of the relationship between the flight measurement techniques and actual measurements. The experimental testing was conducted using monoenergetic electron beams irradiating samples in the dark. An analytical computer model was developed in the NASCAP (NASA Charging Analyzer Program) code. The analytical model material properties for Kapton that controlled backscatter and secondary yield were adjusted to obtain a single set of values that produced reasonable fits for both voltages and currents. The analytical techniques developed in the ground technology investigation have been applied to space flight conditions. Predictions were compared with limited flight data. The agreement is very good indicating that the technique, and NASCAP, can be used to predict spacecraft material charging behavior. Details of the testing, the analytical modelling technique and flight data comparisons are presented.

Stevens, N. J.↗

Studies in the parameterization of cloudiness in climate models and the analysis of radiation fields in general circulation models

Broad-band parameterizations for atmospheric radiative transfer were developed for clear and cloudy skies. These were in the shortwave and longwave regions of the spectrum. These models were compared with other models in an international effort called ICRCCM (Intercomparison of Radiation Codes for Climate Models). The radiation package developed was used for simulations of a General Circulation Model (GCM). A synopsis is provided of the research accomplishments in the two areas separately. Details are available in the published literature.

HARSHVARDHAN↗

The role of global cloud climatologies in validating numerical models

Reliable estimates of the components of the surface radiation budget are important in studies of ocean-atmosphere interaction, land-atmosphere interaction, ocean circulation and in the validation of radiation schemes used in climate models. The methods currently under consideration must necessarily make certain assumptions regarding both the presence of clouds and their vertical extent. Because of the uncertainties in assumed cloudiness, all these methods involve perhaps unacceptable uncertainties. Here, a theoretical framework that avoids the explicit computation of cloud fraction and the location of cloud base in estimating the surface longwave radiation is presented. Estimates of the global surface downward fluxes and the oceanic surface net upward fluxes were made for four months (April, July, October and January) in 1985 to 1986. These estimates are based on a relationship between cloud radiative forcing at the top of the atmosphere and the surface obtained from a general circulation model. The radiation code is the version used in the UCLA/GLA general circulation model (GCM). The longwave cloud radiative forcing at the top of the atmosphere as obtained from Earth Radiation Budget Experiment (ERBE) measurements is used to compute the forcing at the surface by means of the GCM-derived relationship. This, along with clear-sky fluxes from the computations, yield maps of the downward longwave fluxes and net upward longwave fluxes at the surface. The calculated results are discussed and analyzed. The results are consistent with current meteorological knowledge and explainable on the basis of previous theoretical and observational works; therefore, it can be concluded that this method is applicable as one of the ways to obtain the surface longwave radiation fields from currently available satellite data.

HARSHVARDHAN↗

Analytical combustion/emissions research related to the NASA high-speed research program

Increasing the pressure and temperature of the engines of new generation supersonic airliners increases the emissions of nitrogen oxides to a level that would have an adverse impact on the Earth's protective ozone layer. In the process of implementing low emissions combustor technologies, NASA Lewis Research Center has pursued a combustion analysis program to guide combustor design processes, to identify potential concepts of greatest promise, and to optimize them at low cost, with short turn-around time. The approach is to upgrade and apply advanced computer programs for gas turbine applications. Efforts have been made to improve the code capabilities of modeling the physics. Test cases and experiments are used for code validation. To provide insight into the combustion process and combustor design, two-dimensional and three-dimensional codes such as KIVA-II and LeRC 3D have been used. These codes are operational and calculations have been performed to guide low emissions combustion experiments.

Nguyen, H. Lee↗

Comparison of analysis and experiment for gearbox noise

Low contact ratio spur gears were tested in the NASA gear-noise rig to study the noise radiated from the top of the gearbox. Experimental results were compared with a NASA acoustics code to validate the code for predicting transmission noise. The analytical code is based on the boundary element method (BEM) which models the gearbox top as a plate in an infinite baffle. Narrow band vibration spectra measured at 63 nodes on the gearbox top were used to produce input data for the BEM model. The BEM code predicted the total sound power based on the measured vibration. The measured sound power was obtained from an acoustic intensity scan taken near the surface of the gearbox at the same 63 nodes used for vibration measurement. Analytical and experimental results were compared at four different speeds for sound power at each of the narrow band frequencies over the range of 400 to 3200 Hz. Results are also compared for the sound power level at meshing frequency plus three sideband pairs and at selected gearbox resonant frequencies. The difference between predicted and measure sound power is typically less than 3 dB with the predicted value generally less than the measured value.

Oswald, Fred B.↗

Simscape Modeling of a Custom Closed-Volume Tank

The library for Mathworks Simscape does not currently contain a model for a closed volume fluid tank where the ullage pressure is variable. In order to model a closed-volume variable ullage pressure tank, it was necessary to consider at least two separate cases: a vertical cylinder, and a sphere. Using library components, it was possible to construct a rough model for the cylindrical tank. It was not possible to construct a model for a spherical tank, using library components, due to the variable area. It was decided that, for these cases, it would be preferable to create a custom library component to represent each case, using the Simscape language. Once completed, the components were added to models, where filling and draining the tanks could be simulated. When the models were performing as expected, it was necessary to generate code from the models and run them in Trick (a real-time simulation program). The data output from Trick was then compared to the output from Simscape and found to be within acceptable limits.

Simscape↗

A Global Magnetohydrodynamic Model of Jovian Magnetosphere

The goal of this project was to develop a new global magnetohydrodynamic model of the interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 11, Voyager I and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, complicated system. The Jovian aurora also has been monitored for several years. Like auroral observations at Earth, these measurements provide us with a global picture of magnetospheric dynamics. Despite this wide range of observations, we have limited quantitative understanding of the Jovian magnetosphere and how it interacts with the solar wind. For the past several years we have been working toward a quantitative understanding of the Jovian magnetosphere and its interaction with the solar wind by employing global magnetohydrodynamic simulations to model the magnetosphere. Our model has been an explicit MHD code (previously used to model the Earth's magnetosphere) to study Jupiter's magnetosphere. We continue to obtain important insights with this code, but it suffers from some severe limitations. In particular with this code we are limited to considering the region outside of 15RJ, with cell sizes of about 1.5R(sub J). The problem arises because of the presence of widely separated time scales throughout the magnetosphere. The numerical stability criterion for explicit MHD codes is the CFL limit and is given by C(sub max)(Delta)t/(Delta)x less than 1 where C(sub max) is the maximum group velocity in a given cell, (Delta)x is the grid spacing and (Delta)t is the time step. If the maximum wave velocity is C(sub w) and the flow speed is C(sub f), C(sub max) = C(sub w) + C(sub f). Near Jupiter the Alfven wave speed becomes very large (it approaches the speed of light at one Jovian radius). Operating with this time step makes the calculation essentially intractable. Therefore under this funding we have been designing a new MHD model that will be able to compute solutions in the wide parameter regime of the Jovian magnetosphere.

Walker, Raymond J.↗

DDFRG: Double Differential FRaGmentation Models for Proton and Light Ion Production in High Energy Nuclear Collisions: Closed Form, Analytic Formulas for Transport Codes and other Applications

New models for Double-Differential FRaGmentation (DDFRG) cross sections for proton and light ion production from high energy nucleus-nucleus collisions are developed. The proton model employs thermal production from the projectile, central fireball and target sources, and also quasi-elastic direct knockout production. Light ion production cross sections employ a hybrid coalescence model. The models are able to describe a wide range of experimental data with only a limited set of model parameters. Closed form analytic formulas for double-differential cross sections as well as single-differential energy cross sections are developed. The analytic formulas enable highly efficient computation for space radiation transport codes and other applications.

John W Norbury↗

Experimental study and improved modeling of high-temperature air radiation

Spectral emission measurements were conducted in an atmospheric pressure air plasma. The thermodynamic state of the plasma was determined to be close to local thermodynamic equilibrium (LTE), and the temperature profile was measured. The spectrum emitted by the plasma over the range 2,000 to 8,000 A was recorded and calibrated. This spectrum comprises the major radiating molecular bands and atomic lines in air, and is therefore proposed as a benchmark to test radiative calculations. A comparison of these results with the predictions of the NEQAIR code induced several modeling improvements in the code. In particular, radiative transition probabilities and spectroscopic constants were updated, and additional band systems of NO (NO Delta, Epsilon, Beta prime and Gamma prime) were included. Since the C state from which the NO Delta transition originates is predissociated, a simplified collision-predissociation model for this state was added to the code. These changes are presented, and their effect discussed.

Laux, Christophe O.↗

Spacecraft Solar Particle Event (SPE) Shielding: Shielding Effectiveness as a Function of SPE model as Determined with the FLUKA Radiation Transport Code

Analysis of both satellite and surface neutron monitor data demonstrate that the widely utilized Exponential model of solar particle event (SPE) proton kinetic energy spectra can seriously underestimate SPE proton flux, especially at the highest kinetic energies. The more recently developed Band model produces better agreement with neutron monitor data ground level events (GLEs) and is believed to be considerably more accurate at high kinetic energies. Here, we report the results of modeling and simulation studies in which the radiation transport code FLUKA (FLUktuierende KAskade) is used to determine the changes in total ionizing dose (TID) and single-event environments (SEE) behind aluminum, polyethylene, carbon, and titanium shielding masses when the assumed form (i. e., Band or Exponential) of the solar particle event (SPE) kinetic energy spectra is changed. FLUKA simulations have fully three dimensions with an isotropic particle flux incident on a concentric spherical shell shielding mass and detector structure. The effects are reported for both energetic primary protons penetrating the shield mass and secondary particle showers caused by energetic primary protons colliding with shielding mass nuclei. Our results, in agreement with previous studies, show that use of the Exponential form of the event

Koontz, Steve↗

Status Of The UPS Space-Marching Flow Solver

The status of the three-dimensional parabolized Navier-Stokes solver UPS is described. The UPS code, initiated at NASA Ames Research Center in 1986, continues to develop and evolve through application to supersonic and hypersonic flow fields. Hypersonic applications have motivated enhancement of the physical modeling capabilities of the code, specifically real gas modeling, boundary conditions, and turbulence and transition modeling. The UPS code has also been modified to enhance robustness and efficiency in order to be practically used in concert with an optimization code for supersonic transport design. These developments are briefly described along with some relevant results for generic test problems obtained during verification of the enhancements. Included developments and results have previously been published and widely disseminated domestically.

Lawerence, Scott L.↗

Compression behavior of delaminated composite plates

The response of delaminated composite plates to compressive in-plane loads was investigated. The delaminated region may be either circular or elliptical, and may be located between any two plies of the laminate. For elliptical delaminations, the axes of the ellipse may be arbitrarily oriented with respect to the applied loads. A model was developed that describes the stresses, strains, and deformation of the sublaminate created by the delamination. The mathematical model is based on a two dimensional nonlinear plate theory that includes the effects of transverse shear deformation. The model takes into account thermal and moisture induced strains, transverse pressures acting on the sublaminate, and contact between the sublaminate and plate. The solution technique used is the Ritz method. A computationally efficient computer implementation of the model was developed. The code can be used to predict the nonlinear-load-strain behavior of the sublaminate including the buckling load, postbuckling behavior, and the onset of delamination growth. The accuracy of the code was evaluated by comparing the model results to benchmark analytical solutions. A series of experiments was conducted on Fiberite T300/976 graphite/epoxy laminates bonded to an aluminum honeycomb core forming a sandwich panel. Either circles or ellipses made from Teflon film were embedded in the laminates, simulating the presence of a delamination. Each specimen was loaded in compression and the strain history of the sublaminate was recorded far into the postbuckling regime. The extent of delamination growth was evaluated by C-scan examination of each specimen. The experimental data were compared to code predictions. The code was found to describe the data with reasonable accuracy. A sensitivity study examined the relative importance of various material properties, the delamination dimensions, the contact model, the transverse pressure differential, the critical strain energy release rate, and the relative growth direction on the buckling load, the postbuckling behavior, and the growth load of the sublaminate.

Peck, Scott O.↗

Nonlinear static and dynamic finite element analysis of an eccentrically loaded graphite-epoxy beam

The Dynamic Crash Analysis of Structures (DYCAT) and NIKE3D nonlinear finite element codes were used to model the static and implulsive response of an eccentrically loaded graphite-epoxy beam. A 48-ply unidirectional composite beam was tested under an eccentric axial compressive load until failure. This loading configuration was chosen to highlight the capabilities of two finite element codes for modeling a highly nonlinear, large deflection structural problem which has an exact solution. These codes are currently used to perform dynamic analyses of aircraft structures under impact loads to study crashworthiness and energy absorbing capabilities. Both beam and plate element models were developed to compare with the experimental data using the DYCAST and NIKE3D codes.

Fasanella, Edwin L.↗

A method for the modelling of porous and solid wind tunnel walls in computational fluid dynamics codes

Porous wall wind tunnels have been used for several decades and have proven effective in reducing wall interference effects in both low speed and transonic testing. They allow for testing through Mach 1, reduce blockage effects and reduce shock wave reflections in the test section. Their usefulness in developing computational fluid dynamics (CFD) codes has been limited, however, by the difficulties associated with modelling the effect of a porous wall in CFD codes. Previous approaches to modelling porous wall effects have depended either upon a simplified linear boundary condition, which has proven inadequate, or upon detailed measurements of the normal velocity near the wall, which require extensive wind tunnel time. The current work was initiated in an effort to find a simple, accurate method of modelling a porous wall boundary condition in CFD codes. The development of such a method would allow data from porous wall wind tunnels to be used more readily in validating CFD codes. This would be beneficial when transonic validations are desired, or when large models are used to achieve high Reynolds numbers in testing. A computational and experimental study was undertaken to investigate a new method of modelling solid and porous wall boundary conditions in CFD codes. The method utilized experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution was then imposed as a pressure boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the effect of porosity changes on the wall was investigated. Also, the effectiveness of this method when only sparse experimental measurements were available has been investigated. The current work demonstrated this approach for low speed flows and compared the results with experimental data obtained from a heavily instrumented variable porosity test section. The approach developed was simple, computationally inexpensive, and did not require extensive or intrusive measurements of the boundary conditions during the wind tunnel test. It may be applied to both solid and porous wall wind tunnel tests.

Beutner, Thomas John↗

Fatigue analysis of mini-mast space truss

The functional, structural adequacy of a 20 meter long generic space truss (Mini-Mast), subjected to fatigue loading, was examined with respect to the failure modes which are most likely to occur during services. The space truss is made of thin-walled tubes having unidirectional, zero degree layups of Celanese G50 graphite fibers/Narmco 5217 epoxy composites. The approach used to investigate the most probable failure mode of the truss under fatigue loading is to determine the stress level, including the types of stress, in the member first, then followed by failure mode analysis based on the stress level just determined. To begin, an approximate beam-parameter truss (BPT) model is analyzed first, followed by a detailed analysis of the truss using a finite element model (FEM) run with NASTRAN code. The response results of the BPT model are used to compare FEM results and to check any major deviation of trend derived from the FEM. The purpose of the work was to search available fatigue data of the tube material, to conduct approximate dynamical stress analysis of the BPT model, to run detailed dynamical stress analysis of the FEM model using NASTRAN code, and to predict the fatigue life of the truss member based on limited fatigue data.

Hwang, Shoi Y.↗

Computer modeling of pulsed CO2 lasers for lidar applications

The object of this effort is to develop code to enable the accurate prediction of the performance of pulsed transversely excited (TE) CO2 lasers prior to their construction. This is of particular benefit to the NASA Laser Atmospheric Wind Sounder (LAWS) project. A benefit of the completed code is that although developed specifically for the pulsed CO2 laser much of the code can be modified to model other laser systems of interest to the lidar community. A Boltzmann equation solver has been developed which enables the electron excitation rates for the vibrational levels of CO2 and N2, together with the electron ionization and attachment coefficients to be determined for any CO2 laser gas mixture consisting of a combination of CO2, N2, CO, He and CO. The validity of the model has been verified by comparison with published material. The results from the Boltzmann equation solver have been used as input to the laser kinetics code which is currently under development. A numerical code to model the laser induced medium perturbation (LIMP) arising from the relaxation of the lower laser level has been developed and used to determine the effect of LIMP on the frequency spectrum of the LAWS laser output pulse. The enclosed figures show representative results for a laser operating at 0.5 atm. with a discharge cross-section of 4.5 cm to produce a 20 J pulse with aFWHM of 3.1 microns. The first four plots show the temporal evolution of the laser pulse power, energy evolution, LIMP frequency chirp and electric field magnitude. The electric field magnitude is taken by beating the calculated complex electric field and beating it with a local oscillator signal. The remaining two figures show the power spectrum and energy distribution in the pulse as a function of the varying pulse frequency. The LIMP theory has been compared with experimental data from the NOAA Windvan Lidar and has been found to be in good agreement.

Spiers, Gary D.↗

Antenna pattern control using impedance surfaces

The modification of a moment method code is described. Analyses of horn antennas were conducted extensively. Excellent agreements with measured results were observed for pyramidal horn antennas with perfectly conducting surfaces. The predicted results are so accurate that even the finest ripple structures in the far-field radiation patterns are computed correctly. Preliminary results for the radiation patterns of pyramidal horns with impedance surfaces were also obtained. Discussions of using the code in accurate modeling of practical pyramidal horn antennas are also included. After the code for the analysis of the horn antennas with impedance materials is finalized, the synthesis problem will be examined in the coming research period.

Balanis, Constantine A.↗