Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous variables”

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,387 records · Page 77

Effects of Heat Flux, Oxygen Concentration and Glass Fiber Volume Fraction on Pyrolysate Mass Flux from Composite Solids

Experimental work on the effects of heat flux, oxygen concentration and glass fiber volume fraction on pyrolysate mass flux from samples of polypropylene/glass fiber composite (PP/G) is underway. The research is conducted as part of a larger project to develop a test methodology for flammability of materials, particularly composites, in the microgravity and variable oxygen concentration environment of spacecraft and space structures. Samples of PP/G sized at 30 x 30 x 10 mm are flush mounted in a flow tunnel, which provides a flow of oxidizer over the surface of the samples at a fixed value of 1 m/s and oxygen concentrations varying between 18 and 30%. Each sample is exposed to a constant external radiant heat flux at a given value, which varies between tests from 10 to 24 kW/sq m. Continuous sample mass loss and surface temperature measurements are recorded for each test. Some tests are conducted with an igniter and some are not. In the former case, the research goal is to quantify the critical mass flux at ignition for the various environmental and material conditions described above. The later case generates a wider range of mass flux rates than those seen prior to ignition, providing an opportunity to examine the protective effects of blowing on oxidative pyrolysis and heating of the surface. Graphs of surface temperature and sample mass loss vs. time for samples of 30% PPG at oxygen concentrations of 18 and 21% are presented in the figures below. These figures give a clear indication of the lower pyrolysis rate and extended time to ignition that accompany a lower oxygen concentration. Analysis of the mass flux rate at the time of ignition gives good repeatability but requires further work to provide a clear indication of mass flux trends accompanying changes in environmental and material properties.

Rich, D. B.↗

Uncertainty Analysis for a Jet Flap Airfoil

An analysis of variance (ANOVA) study was performed to quantify the potential uncertainties of lift and pitching moment coefficient calculations from a computational fluid dynamics code, relative to an experiment, for a jet flap airfoil configuration. Uncertainties due to a number of factors including grid density, angle of attack and jet flap blowing coefficient were examined. The ANOVA software produced a numerical model of the input coefficient data, as functions of the selected factors, to a user-specified order (linear, 2-factor interference, quadratic, or cubic). Residuals between the model and actual data were also produced at each of the input conditions, and uncertainty confidence intervals (in the form of Least Significant Differences or LSD) for experimental, computational, and combined experimental / computational data sets were computed. The LSD bars indicate the smallest resolvable differences in the functional values (lift or pitching moment coefficient) attributable solely to changes in independent variable, given just the input data points from selected data sets. The software also provided a collection of diagnostics which evaluate the suitability of the input data set for use within the ANOVA process, and which examine the behavior of the resultant data, possibly suggesting transformations which should be applied to the data to reduce the LSD. The results illustrate some of the key features of, and results from, the uncertainty analysis studies, including the use of both numerical (continuous) and categorical (discrete) factors, the effects of the number and range of the input data points, and the effects of the number of factors considered simultaneously.

Green, Lawrence L.↗

Ultra-Wideband Radar Measurements of Thickness of Snow Over Sea Ice

An accurate knowledge of snow thickness and its variability over sea ice is crucial for determining the overall polar heat and freshwater budget, which influences the global climate. Recently, algorithms have been developed to extract snow thicknesses from passive microwave satellite data. However, validation of these data over the large footprint of the passive microwave sensor has been a challenge. The only method used thus far has been with meter sticks during ship cruises. To address this problem, we developed an ultra wideband frequency-modulated continuous-wave (FM-CW) radar to measure snow thickness over sea ice. We made snow-thickness measurements over Antarctic sea ice by operating the radar from a sled during September and October, 2003. We performed radar measurements over 11 stations with varying snow thickness between 4 and 85 cm. We observed excellent agreement between radar estimates of snow thickness with physical measurements, achieving a correlation coefficient of 0.95 and a vertical resolution of about 3 cm.

Kanagaratnam, P.↗

Trajectory Optimization: OTIS 4

The latest release of the Optimal Trajectories by Implicit Simulation (OTIS4) allows users to simulate and optimize aerospace vehicle trajectories. With OTIS4, one can seamlessly generate optimal trajectories and parametric vehicle designs simultaneously. New features also allow OTIS4 to solve non-aerospace continuous time optimal control problems. The inputs and outputs of OTIS4 have been updated extensively from previous versions. Inputs now make use of objectoriented constructs, including one called a metastring. Metastrings use a greatly improved calculator and common nomenclature to reduce the user s workload. They allow for more flexibility in specifying vehicle physical models, boundary conditions, and path constraints. The OTIS4 calculator supports common mathematical functions, Boolean operations, and conditional statements. This allows users to define their own variables for use as outputs, constraints, or objective functions. The user-defined outputs can directly interface with other programs, such as spreadsheets, plotting packages, and visualization programs. Internally, OTIS4 has more explicit and implicit integration procedures, including high-order collocation methods, the pseudo-spectral method, and several variations of multiple shooting. Users may switch easily between the various methods. Several unique numerical techniques such as automated variable scaling and implicit integration grid refinement, support the integration methods. OTIS4 is also significantly more user friendly than previous versions. The installation process is nearly identical on various platforms, including Microsoft Windows, Apple OS X, and Linux operating systems. Cross-platform scripts also help make the execution of OTIS and post-processing of data easier. OTIS4 is supplied free by NASA and is subject to ITAR (International Traffic in Arms Regulations) restrictions. Users must have a Fortran compiler, and a Python interpreter is highly recommended.

Riehl, John P.↗

Variability of Solar Radiation under Cloud-Free Skies in China: The Role of Aerosols

In this study, we analyzed long-term surface global and diffuse solar radiation, aerosol single scattering albedo (SSA), and relative humidity (RH) from China. Our analysis reveals that much of China experienced significant decreases in global solar radiation (GSR) and increases in diffuse solar radiation under cloud-free skies between the 1960s and 1980s. With RH and aerosol SSA being rather constant during that time period, we suggest that the increasing aerosol loading from emission of pollutants is responsible for the observed reduced GSR and increased diffuse radiation in cloud-free skies. Although pollutant emissions continue to increase after the 1980s, the increment of aerosol SSA since 1980s can partly explain the transition of GSR from a decreasing trend to no apparent trend around that time. Preliminary analysis is also provided on the potential role of RH in affecting the global and diffuse solar radiation reaching the earth surface.

Qian, Yun↗

Relativistic Astrophysics in Black Hole and Low-Mass Neutron Star Binaries LTSA98

My group, in close collaboration with Dr. Zhang's group at University of Alabama-Huntsville, have been systematically analyzing and re-analyzing a substantial amount of archival data from previous and ongoing X-ray missions, in order to study possible relativistic effects around stellar-mass black holes and neutron stars. Our effort has been focused primarily on the data from the Rossi X-ray Timing Explorer. We carefully studied interesting quasi-periodic X-ray variability in newly discovered black hole candidates (XTE J1859+226 and XTE J1550-564), which, as we had proposed earlier, could be caused by general relativistic process (e.g., frame dragging) around the central black hole. We also discovered an intriguing temporal correlation between X-ray photons at different energies that is associated with the quasi-periodic signals of interest. The results provided new insights into the physical origin of the phenomena. Furthermore, we studied the spectral lines of black hole candidates which provide another avenue for studying general relativistic processes around black holes. The lines-may originate in the relativistic jets (which could be powered by the spin of the black hole) or in the disk around the black hole, as in the cases of 4U 1630-47 and GX 339-4 (two well-known black hole candidates), and may thus be distorted or shifted due to relativistic effects. Of course, neutron star systems were not forgotten either. After examining the properties of newly discovered fast quasi-periodic variability (at kiloHertz) associated with such systems, we proposed a relativistic model to explain the origin of the signals. We have also started to use new great observatories in orbit (such as Chandra and XMM-Newton) to observe the sources that are of interest to us. Finally, interesting results were also been obtained from our collaborations with other groups who are interested in some of the same objects. Such collaborative efforts have greatly enhanced the project and will likely continue in the future.

Cui, Wei↗

Preliminary Results from a Coordinated Hisaki/Chandra/XMM-Newton Study of the Jovian Aurora and Io Plasma Torus

We present preliminary results from a coordinated Hisaki/Chandra/XMM-Newton observational campaign of the Jovian aurora and Io plasma torus. The data were taken over a three week period in April, 2014. Jupiter was observed continuously with Hisaki, six times with the Chandra/HRC instrument for roughly 12 hours per observation, and twice by XMM-Newton. The goal of this observational campaign was to understand how energy and matter are exchanged between the Jovian aurora, the IPT, and the Solar wind. X-ray observations provide key diagnostics on highly stripped ions and keV electrons in the Jovian magnetosphere. We use the temporal, spatial, and spectral capabilities of the three instruments to search for correlated variability between the Solar wind, the EUV-emitting plasma of the IPT and UV aurora, and the ions responsible for the X-ray aurora. Preliminary analysis suggests a strong 45 min periodicity in the EUV emission from the electron aurora. There is some evidence for complex variability of the X-ray auroras on scales of tens of minutes. There is also clear morphological changes in the X-ray aurora that do not appear to be correlated with either variations in the IPT or Solar wind.

Kraft, Ralph↗

Torsional vibrations and dynamic loads in a basic planetary gear system

An iterative method has been developed for analyzing dynamic loads in a light weight basic planetary gear system. The effects of fixed, semi-floating, and fully-floating sun gear conditions have been emphasized. The load dependent variable gear mesh stiffness were incorporated into a practical torsional dynamic model of a planetary gear system. The dynamic model consists of input and output units, shafts, and a planetary train. In this model, the sun gear has three degrees of freedom; two transverse and one rotational. The planets, ring gear, and the input and output units have one degree of freedom, (rotation) thus giving a total of nine degrees of freedoms for the basic system. The ring gear has a continuous radial support. The results indicate that the fixed sun gear arrangement with accurate or errorless gearing offers in general better performance than the floating sun gear system.

August, R.↗

The nature of the cold filaments in the California Current system

The velocity fields and water properties associated with cold filaments in the California Current are described on the basis of data from the Coastal Transition Zone experiment. Combined with previous field surveys and satellite imagery, these show seasonal variability with maximum dynamic height ranges and velocities in summer and minimum values in late winter and early spring. North of Point Arena in spring-summer, the flow field on the outer edge of the cold water exhibits the character of a meandering jet, carrying fresh, nutrient-poor water from the farther north on its offshore side and cold, salty, nutrient-rich water on its inshore side. At Point Arena in midsummer, the jet often flows offshore and continues south without meandering back onshore as strongly as it does farther north. At the surface, the jet often separates biological communities and may appear as a barrier to cross-jet transport, especially north of Point Arena in March-May.

Strub, P. T.↗

Assimilation of MLS and OMI Ozone Data

Ozone data from Aura Microwave Limb Sounder (MLS) and Ozone Monitoring Instrument (OMI) were assimilated into the ozone model at NASA's Global Modeling and Assimilation Office (GMAO). This assimilation produces ozone fields that are superior to those from the operational GMAO assimilation of Solar Backscatter Ultraviolet (SBUV/2) instrument data. Assimilation of Aura data improves the representation of the "ozone hole" and the agreement with independent Stratospheric Aerosol and Gas Experiment (SAGE) III and ozone sonde data. Ozone in the lower stratosphere is captured better: mean state, vertical gradients, spatial and temporal variability are all improved. Inclusion of OMI and MLS data together, or separately, in the assimilation system provides a way of checking how consistent OMI and MLS data are with each other, and with the ozone model. We found that differences between OMI total ozone column data and model forecasts decrease after MLS data are assimilated. This indicates that MLS stratospheric ozone profiles are consistent with OMI total ozone columns. The evaluation of error characteristics of OMI and MLS ozone will continue as data from newer versions of retrievals becomes available. We report on the initial step in obtaining global assimilated ozone fields that combine measurements from different Aura instruments, the ozone model at the GMAO, and their respective error characteristics. We plan to use assimilated ozone fields in estimation of tropospheric ozone. We also plan to investigate impacts of assimilated ozone fields on numerical weather prediction through their use in radiative models and in the assimilation of infrared nadir radiance data from NASA's Advanced Infrared Sounder (AIRS).

Stajner, I.↗

Metrics for Linear Kinematic Features in Sea Ice

The treatment of leads as cracks or discontinuities (see Coon et al. presentation) requires some shift in the procedure of evaluation and comparison of lead-resolving models and their validation against observations. Common metrics used to evaluate ice model skills are by and large an adaptation of a least square "metric" adopted from operational numerical weather prediction data assimilation systems and are most appropriate for continuous fields and Eilerian systems where the observations and predictions are commensurate. However, this class of metrics suffers from some flaws in areas of sharp gradients and discontinuities (e.g., leads) and when Lagrangian treatments are more natural. After a brief review of these metrics and their performance in areas of sharp gradients, we present two new metrics specifically designed to measure model accuracy in representing linear features (e.g., leads). The indices developed circumvent the requirement that both the observations and model variables be commensurate (i.e., measured with the same units) by considering the frequencies of the features of interest/importance. We illustrate the metrics by scoring several hypothetical "simulated" discontinuity fields against the lead interpreted from RGPS observations.

Levy, G.↗

Numerical Implementation of a Multiple-ISV Thermodynamically-Based Work Potential Theory for Modeling Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗

A Thermodynamically-Based Mesh Objective Work Potential Theory for Predicting Intralaminar Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗

Intercomparison of Field Methods for Acquiring Ground Reflectance at Railroad Valley Playa for Spectral Calibration of Satellite Data

Ground reflectance was acquired at the Railroad Valley Playa calibration site in Nevada USA using different methods of collection. The data was collected near the time and date of Landsat 8 OLI and Sentinel-2 satellite overpasses so an inter-comparison could be made with the reflectance products to determine which method was more suitable for vicarious calibration. The field spectrometers and reference panels were characterized before the field campaign. A continuous acquisition method was compared to stop and measure collections. Both acquisition methods were collected along an 80 m east-west transect as well as for a series of north-south transects over an 80 x 320 m area, with the stop and measure method being performed at random sampling locations. The measurements were performed using two field spectrometers by three teams of two people to compare the repeatability. The aim of the field campaign was to determine the variability due to the operator and the method of collection.

Lau, Ian C.↗

Intercomparison of Field Methods for Acquiring Ground Reflectance at Railroad Valley Playa for Spectral Calibration of Satellite Data

Ground reflectance was acquired at the Railroad Valley Playa calibration site in Nevada USA using different methods of collection. The data was collected near the time and date of Landsat 8 OLI and Sentinel-2 satellite overpasses so an inter-comparison could be made with the reflectance products to determine which method was more suitable for vicarious calibration. The field spectrometers and reference panels were characterized before the field campaign. A continuous acquisition method was compared to stop and measure collections. Both acquisition methods were collected along an 80 m east-west transect as well as for a series of north-south transects over an 80 x 320 m area, with the stop and measure method being performed at random sampling locations. The measurements were performed using two field spectrometers by three teams of two people to compare the repeatability. The aim of the field campaign was to determine the variability due to the operator and the method of collection.

Lau, Ian C.↗

Winds at the Mars 2020 Landing Site: 1. Near-Surface Wind Patterns at Jezero Crater

This is the first part of a two-part paper. NASA's Mars 2020 Perseverance rover measured winds on the Jezero crater floor close to the delta of an ancient river. A mostly repeatable diurnal cycle was observed and presented two regimes: (a) a convective regime, from dawn to sunset, with average easterly to southeasterly winds, during which maximum wind speeds were measured, and (b) a nighttime regime with westerly-northwesterly winds followed by a relatively calm period with highly variable wind directions as a function of sol and time of night. The timing and magnitude of the observed regimes are consistent with primary control by regional and local slope flows. Data suggest that the surface circulation at Jezero region in northern spring and summer is highly unaffected by large-scale circulation except during particular periods in the diurnal cycle or generally during dust storms, which is supported by MarsWRF model simulations. Consequently, the observed seasonal variability was weak. However, sol-to-sol and seasonal variability were measured, most of it during certain nighttime periods. Traveling waves consistent with baroclinic instability were clearly observed in surface winds at L_s ∼ 75°. The early MY36/2022A regional dust storm at L_s ∼ 153° disturbed the wind patterns with changes suggesting enhanced tidal flows. After sunset, the dust storm also produced detectable gravity wave activity, increasing the mixing in the nighttime planetary boundary layer during storm conditions. Inferred wind directions from dust devil movies strongly suggest that prevailing winds continued to be slope-driven during the late summer, fall and early winter seasons.

D. Viúdez-Moreiras↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

Implementing neural nets with programmable logic

Networks of Boolean programmable logic modules are presented as one purely digital class of artificial neural nets. The approach contrasts with the continuous analog framework usually suggested. Programmable logic networks are capable of handling many neural-net applications. They avoid some of the limitations of threshold logic networks and present distinct opportunities. The network nodes are called dynamically programmable logic modules. They can be implemented with digitally controlled demultiplexers. Each node performs a Boolean function of its inputs which can be dynamically assigned. The overall network is therefore a combinational circuit and its outputs are Boolean global functions of the network's input variables. The approach offers definite advantages for VLSI implementation, namely, a regular architecture with limited connectivity, simplicity of the control machinery, natural modularity, and the support of a mature technology.

Vidal, Jacques J.↗