Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 325 records · Page 18

Parametric instabilities of circularly polarized Alfven waves in high-beta plasmas

CGL relations including the effect of finite ion Larmor radius are used to consider a class of parametric instabilities of finite-amplitude, circularly polarized Alfven waves in high-beta plasmas. The disperison relation governing the instabilities is a sixth-order polynomial which is solved numerically. There are two types of instabilities: a modulational instability at k is less than k(0) and a relatively weak and narrow bandwidth instability at k is less than approximately k(0), where k and k(0) are the wavenumbers of the unstable density fluctuation and the 'pump' wave, respectively. It is shown that these instabilities can occur for left-handed pump waves and that the modulational instability is unstable over a very broad band in k with a maximum growth rate at finite k is not equal to 0.

Hamabata, Hiromitsu↗

Power spectral density estimation by spline smoothing in the frequency domain

An approach, based on a global averaging procedure, is presented for estimating the power spectrum of a second order stationary zero-mean ergodic stochastic process from a finite length record. This estimate is derived by smoothing, with a cubic smoothing spline, the naive estimate of the spectrum obtained by applying FFT techniques to the raw data. By means of digital computer simulated results, a comparison is made between the features of the present approach and those of more classical techniques of spectral estimation.

Defigueiredo, R. J. P.↗

Power spectral density estimation by spline smoothing in the frequency domain.

An approach, based on a global averaging procedure, is presented for estimating the power spectrum of a second order stationary zero-mean ergodic stochastic process from a finite length record. This estimate is derived by smoothing, with a cubic smoothing spline, the naive estimate of the spectrum obtained by applying Fast Fourier Transform techniques to the raw data. By means of digital computer simulated results, a comparison is made between the features of the present approach and those of more classical techniques of spectral estimation.-

De Figueiredo, R. J. P.↗

Wall boundary equations with slip and catalysis for multicomponent, nonequilibrium gas flows

Boundary equations obtained for a low Reynolds number, high enthalpy gas flow in regions of velocity slip and temperature jump are presented. The formulation treats, through a first-order distribution function used to include multicomponent diffusion, a multicomponent gas mixture that may be in nonequilibrium with finite-rate catalytic recombination occurring on the wall. In the boundary equations, which are obtained for use in flow-field calculations applicable to low-density flow regimes, a simplified gas/wall interaction is assumed wherein individual atoms or molecules either reflect specularly off the wall or stick and are fully accommodated. Fluxes in terms of evaluated integrals over the distribution function and integrals necessary for determining fluxes are given.

Scott, C. D.↗

Saturn - UBV photoelectric pinhole scans of the disk. II

During the 1980 Saturn apparition, UBV pinhole scans of the disk were obtained with a photoelectric area-scanning photometer. An analysis of these data reveals that the atmosphere of Saturn can be represented by a finite clear H2 layer overlying a semiinfinite absorbent aerosol haze. The extent of the clear H2 region appears to be latitude-dependent; the H2 column density varies systematically from about 15 km-am over the equatorial and polar regions to about 31 km-am at temperate latitudes. A previous conclusion that the aerosol haze is strongly absorbent in the UV is confirmed; its effective U-band single-scattering albedo is about 0.4. Latitudinal disk structure at visual wavelengths appears to be the result of local variations in the volume density of absorbent particles in the aerosol layer.

Price, M. J.↗

Transfer impedance measurements of the space shuttle Solid Rocket Motor (SRM) joints, wire meshes and a carbon graphite motor case

Using three different current injection methods and a simple voltage probe, transfer impedances for Solid Rocket Motor (SRM) joints, wire meshes, aluminum foil, Thorstrand and a graphite composite motor case were measured. In all cases, the surface current distribution for the particular current injection device was calculated analytically or by finite difference methods. The results of these calculations were used to generate a geometric factor which was the ratio of total injected current to surface current density. The results were validated in several ways. For wire mesh measurements, results showed good agreement with calculated results for a 14 by 18 Al screen. SRM joint impedances were independently verified. The filiment wound case measurement results were validated only to the extent that their curve shape agrees with the expected form of transfer impedance for a homogeneous slab excited by a plane wave source.

Papazian, Peter B.↗

On positivity preserving finite volume schemes for compressible Euler equations

Positivity preserving property of first and higher order finite volume schemes for one and two dimensional compressible Euler equations of gas dynamics is considered. A general framework is established which shows the positivity of density and pressure whenever the underlying one dimensional first order building block based on exact or approximate Riemann solver and the reconstruction are both positivity preserving. Appropriate limitation to achieve high order positivity preserving reconstruction is described.

Perthame, Benoit↗

Experimental Verification of a Progressive Damage Model for IM7/5260 Laminates Subjected to Tension-Tension Fatigue

The durability and damage tolerance of laminated composites are critical design considerations for airframe composite structures. Therefore, the ability to model damage initiation and growth and predict the life of laminated composites is necessary to achieve structurally efficient and economical designs. The purpose of this research is to experimentally verify the application of a continuum damage model to predict progressive damage development in a toughened material system. Damage due to monotonic and tension-tension fatigue was documented for IM7/5260 graphite/bismaleimide laminates. Crack density and delamination surface area were used to calculate matrix cracking and delamination internal state variables to predict stiffness loss in unnotched laminates. A damage dependent finite element code predicted the stiffness loss for notched laminates with good agreement to experimental data. It was concluded that the continuum damage model can adequately predict matrix damage progression in notched and unnotched laminates as a function of loading history and laminate stacking sequence.

Coats, Timothy W.↗

Modal Propellant Gauging: High-resolution and non-invasive gauging of both settled and unsettled liquids in reduced gravity.

The modal response of a liquid-filled tank to external acoustic excitation can be used to infer with high resolutionthe mass of contained liquid, the mass flow rate of liquids into and out of the tank, and changes in tank pressure.Both contained liquid mass and internal ullage pressure affect the modal response of the tank walls through fluidmass-loading of the tank walls and pressure-induced wall stiffening, respectively. Modal Propellant Gauging refersto the technology that exploits these shifts in modal frequencies to infer the mass of propellant in a tank. MPG is anon-invasive gauging technology that has demonstrated gauging resolutions of 1% for settled propellants and 2-3%for unsettled, sloshing propellants. Extensive parabolic flight testing of the MPG system on model tanks has beenconducted to validate the technology in reduced gravity. MPG testing on a qualification tank for the Orion Program'sEuropean Service Module has also been conducted and is reported here. Finite element modeling of the OrionESM "upper" tank is discussed and compared with measurement data. Three computational approaches to massdetermination, Peak Tracking, Point Sensor, and Spectral Density methods, are described here. Use cases are definedand analyzed in the context of the Orion ESM Qualification tank data, and an implementation scheme for continuousmass gauging on the Orion ESM is discussed.Keywords: Propellant Mass Gauging, Modal Analysis, Orion ESM

Propellant Mass Gauging↗

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. The transient finite-element analysis is done using TACS, and the cell-discharge battery model uses OpenMDAO and dymos. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗

A Flexible Method for Producing F.E.M. Analysis of Bone Using Open-Source Software

This project, performed in support of the NASA GRC Space Academy summer program, sought to develop an open-source workflow methodology that segmented medical image data, created a 3D model from the segmented data, and prepared the model for finite-element analysis. In an initial step, a technological survey evaluated the performance of various existing open-source software that claim to perform these tasks. However, the survey concluded that no single software exhibited the wide array of functionality required for the potential NASA application in the area of bone, muscle and bio fluidic studies. As a result, development of a series of Python scripts provided the bridging mechanism to address the shortcomings of the available open source tools. The implementation of the VTK library provided the most quick and effective means of segmenting regions of interest from the medical images; it allowed for the export of a 3D model by using the marching cubes algorithm to build a surface mesh. To facilitate the development of the model domain from this extracted information required a surface mesh to be processed in the open-source software packages Blender and Gmsh. The Preview program of the FEBio suite proved to be sufficient for volume filling the model with an unstructured mesh and preparing boundaries specifications for finite element analysis. To fully allow FEM modeling, an in house developed Python script allowed assignment of material properties on an element by element basis by performing a weighted interpolation of voxel intensity of the parent medical image correlated to published information of image intensity to material properties, such as ash density. A graphical user interface combined the Python scripts and other software into a user friendly interface. The work using Python scripts provides a potential alternative to expensive commercial software and inadequate, limited open-source freeware programs for the creation of 3D computational models. More work will be needed to validate this approach in creating finite-element models.

gravitational physiology↗

Monte Carlo simulation of wave sensing with a short pulse radar

A Monte Carlo simulation is used to study the ocean wave sensing potential of a radar which scatters short pulses at small off-nadir angles. In the simulation, realizations of a random surface are created commensurate with an assigned probability density and power spectrum. Then the signal scattered back to the radar is computed for each realization using a physical optics analysis which takes wavefront curvature and finite radar-to-surface distance into account. In the case of a Pierson-Moskowitz spectrum and a normally distributed surface, reasonable assumptions for a fully developed sea, it has been found that the cumulative distribution of time intervals between peaks in the scattered power provides a measure of surface roughness. This observation is supported by experiments.

Levine, D. M.↗

A magnetohydrodynamical model of the inner coma of Comet Halley

The magnetometer onboard the Giotto spacecraft observed a diamagnetic cavity surrounding the nucleus of Comet Halley. The location of the boundary of this diamagnetic cavity is determined by a balance between an inward magnetic pressure gradient force and an outward ion-neutral drag force, associated with collisions between the outwardly flowing neutrals and the stagnated ions. A one-dimensional time-dependent MHD model has been developed for the inner coma of Comet Halley, and includes ion-neutral collisions, photochemical production and loss of plasma, and finite conductivity. This model is used to investigate the plasma dynamics in the vicinity of the diamagnetic cavity boundary surface. A narrow transition layer with enhanced plasma density is shown to exist just inside the boundary, although a full understanding of this layer will require a two- or three-dimensional MHD model. The flux of cometary ions into this shocklike layer is removed by electron-ion recombination. The thickness of this layer is determined by the Mach number of the incident flow.

Cravens, T. E.↗

The PDF method for turbulent combustion

Probability Density Function (PDF) methods provide a means of calculating the properties of turbulent reacting flows. They have been successfully applied to many turbulent flames, including some with finite rate kinetic effects. Here the methods are reviewed with an emphasis on computational issues and their application to turbulent combustion.

Pope, S. B.↗

Structural Anomalies Detected in Ceramic Matrix Composites Using Combined Nondestructive Evaluation and Finite Element Analysis (NDE and FEA)

Most reverse engineering approaches involve imaging or digitizing an object and then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. The rapid prototyping technique builds high-quality physical prototypes directly from computer-aided design files. This fundamental technique for interpreting and interacting with large data sets is being used here via Velocity2 (an integrated image-processing software, ref. 1) using computed tomography (CT) data to produce a prototype three-dimensional test specimen model for analyses. A study at the NASA Glenn Research Center proposes to use these capabilities to conduct a combined nondestructive evaluation (NDE) and finite element analysis (FEA) to screen pretest and posttest structural anomalies in structural components. A tensile specimen made of silicon nitrite (Si3N4) ceramic matrix composite was considered to evaluate structural durability and deformity. Ceramic matrix composites are being sought as candidate materials to replace nickel-base superalloys for turbine engine applications. They have the unique characteristics of being able to withstand higher operating temperatures and harsh combustion environments. In addition, their low densities relative to metals help reduce component mass (ref. 2). Detailed three-dimensional volume rendering of the tensile test specimen was successfully carried out with Velocity2 (ref. 1) using two-dimensional images that were generated via computed tomography. Subsequent, three-dimensional finite element analyses were performed, and the results obtained were compared with those predicted by NDE-based calculations and experimental tests. It was shown that Velocity2 software can be used to render a three-dimensional object from a series of CT scan images with a minimum level of complexity. The analytical results (ref. 3) show that the high-stress regions correlated well with the damage sites identified by the CT scans and the experimental data. Furthermore, modeling of the voids collected via NDE offered an analytical advantage that resulted in more accurate assessments of the material s structural strength. The top figure shows a CT scan image of the specimen test section illustrating various hidden structural entities in the material and an optical image of the test specimen considered in this study. The bottom figure represents the stress response predicted from the finite element analyses (ref .3 ) for a selected CT slice where it clearly illustrates the correspondence of the high stress risers due to voids in the material with those predicted by the NDE. This study is continuing, and efforts are concentrated on improving the modeling capabilities to imitate the structural anomalies as detected.

Abdul-Aziz, Ali↗

Correlated Electromagnetic Levitation Actuator: A Reaction Sphere Based Attitude Control System

To address problems experienced by current reaction wheels and control moment gyroscopebased attitude control systems (ACS), researchers at NASA’s Marshall Space Flight Center have begun developing a reaction sphere actuator based on correlated electromagnetic levitation that will be immune to destructive bearing friction, momentum saturation, and gimbal lock. The Correlated Electromagnetic Levitation Actuator (CELA) advances the state of the art of reaction sphere ACSs by employing the concept of correlated magnetics. It is a frictionless, direct-drive reaction sphere that harnesses a unique technology with an array of applications across multiple disciplines. Correlated electromagnets function in a manner that is analogous to a matched filter; the convolution of two signals is peaked at the index representing the greatest match. For CELA, the signals are the patterns of magnetic flux density as a function of position. The magnitude of the convolution equates to an attractive or repulsive force, and these forces can be azimuthal or radial. The development of CELA is based in four distinct disciplines: Advanced Manufacturing, Prototype Development, Electromagnetic Modeling, and Controls. We are developing novel manufacturing techniques required to build arrays of permanent and electromagnet dipoles on curved surfaces. To print the permanent magnetic array, we have developed a probe with pyramidal magnets that will reside on a robotic arm to induce localized magnetic fields on a surface. The probe also includes the ability to erase dipole patterns from a permanent magnet by heating the surface to its Curie temperature. A number of test articles and prototypes have been developed using additive manufacturing methods. These prototypes have included hemispherical motors to test the drive algorithm, and a levitation test bed that demonstrates a magnetic bearing method based on attractive magnetic forces and ratiometric Hall effect sensors. We developed an array of electromagnetic dipoles on a printed circuit board (PCB) with individual H-bridges controlling each coil. This device created various flux density patterns and we measured their magnetic fields using a custom Hall effect 3-D probe and a LabVIEW virtual instrument. These data will serve as a benchmark for characterizing the accuracy of future models. Current work is focused on modeling the magnetic fields of our prototype arrays using COMSOL Finite Element Analysis and verifying the model against our test data. Accurate modeling will allow us to quickly test new patterns of electromagnets and their macro behavior. Eventually, the magnetic field models will be implemented in our controls simulations to facilitate precise control of the reaction sphere. Initial model results agree with field measurements to within 1 G (5% of measured flux density). Currently, we are testing different material properties of the electromagnets and their magnetic fields and thermal effects. These results will be used to refine the design of the electromagnetic dipoles. Our control efforts have centered on developing commutation, levitation, and field pattern shaping hardware in the form of breadboards and PCBs with software running on a local microcontroller. In addition, our partners developed MATLAB Simulink models to demonstrate a PID controller thatmitigates disturbance forces resulting from the interaction of drive and levitation magnetics. Finally, we have designed a three-axis test stand that will be used in future work to demonstrate CELA’s orientation control capability.

controls↗

Correlated Electromagnetic Levitation Actuator: A Reaction Sphere-Based Attitude Control System

To address problems experienced by current reaction wheels and control moment gyroscopebased attitude control systems (ACS), researchers at National Aeronautics and Space Administration’s (NASA’s) Marshall Space Flight Center (MSFC) have begun developing a reaction sphere actuator based on correlated electromagnetic levitation that will be immune to destructive bearing friction, momentum saturation, and gimbal lock. The Correlated Electromagnetic Levitation Actuator (CELA) advances the state of the art of reaction sphere ACSs by employing the concept of correlated magnetics. It is a frictionless, direct-drive reaction sphere that harnesses a unique technology with an array of applications across multiple disciplines. Correlated electromagnets function in a manner that is analogous to a matched filter; the convolution of two signals is peaked at the index representing the greatest match. For CELA, the signals are the patterns of magnetic flux density as a function of position. The magnitude of the convolution equates to an attractive or repulsive force, and these forces can be azimuthal or radial. The development of CELA is based in four distinct disciplines: Advanced Manufacturing, Prototype Development, Electromagnetic Modeling, and Controls. We are developing novel manufacturing techniques required to build arrays of permanent and electromagnet dipoles on curved surfaces. To print the permanent magnetic array, we have developed a probe with pyramidal magnets that will reside on a robotic arm to induce localized magnetic fields on a surface. We have also used high temperature ovens to erase dipole patterns from a permanent magnet by heating the surface to its Curie temperature. A number of test articles and prototypes have been developed using additive manufacturing methods. These prototypes have included hemispherical motors to test the drive algorithm, and a levitation test bed that demonstrates a magnetic bearing method based on attractive magnetic forces and ratiometric Hall effect sensors. We developed an array of electromagnetic dipoles on a printed circuit board (PCB) with individual H-bridges controlling each coil. This device created various flux density patterns and we measured their magnetic fields using a custom Hall effect 3-D probe and a LabVIEW virtual instrument. These data will serve as a benchmark for characterizing the accuracy of future models. Current work is focused on modeling the magnetic fields of our prototype arrays using COMSOL finite element analysis (FEA) and verifying the model against our test data. Accurate modeling will allow us to quickly test new patterns of electromagnets and their macro behavior. Eventually, the magnetic field models will be implemented in our controls simulations to facilitate precise control of the reaction sphere. Initial model results agree with field measurements to within 1 G (5% of measured flux density). Currently, we are testing different material properties of the electromagnets and their magnetic fields and thermal effects. These results will be used to refine the design of the electromagnetic dipoles. Our control efforts have centered on developing commutation, levitation, and field pattern shaping hardware in the form of breadboards and PCBs with software running on a local microcontroller. In addition, our partners developed MATLAB Simulink models to demonstrate a PID controller that mitigates disturbance forces resulting from the interaction of drive and levitation magnetics. Finally, we have designed a three-axis test stand that will be used in future work to demonstrate CELA’s orientation control capability.

reaction sphere↗