Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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 829 records · Page 46

The Impact of Firn Models on Ultrawideband Brightness Temperatures in the Partially Coherent Model

The density profile of a polar ice sheet is an important parameter for the estimation of ice mass balance. Wave reflections caused by density variations are also a key uncertainty in the retrieval of ice sheet temperature profiles in Ultra-Wide band radiometry. In this paper, we examine different firn density profile models and analyze the subsurface reflections they cause using an analytical partially coherent approach. We also examine firn density profiles obtained from borehole measurements, from past UWBRAD modeling studies, from a community firn model, and from snow radar echo measurements. In previous studies, the ice sheet has been model as a 1D random medium with density variations in depth. However, horizontal density variations also exist, so that the ice sheet is a 3D random medium. Analyses using the partially coherent model show that in the presence of horizontal fluctuations, contributions from short scale variations vanish as the horizontal correlation length decreases due to the diffraction of waves.

Haokui Xu↗

Effect of Trailing Edge Flow Injection on Fan Noise and Aerodynamic Performance

An experimental investigation using trailing edge blowing for reducing fan rotor/guide vane wake interaction noise was completed in the NASA Glenn 9- by 15-foot Low Speed Wind Tunnel. Data were acquired to measure noise, aerodynamic performance, and flow features for a 22" tip diameter fan representative of modern turbofan technology. The fan was designed to use trailing edge blowing to reduce the fan blade wake momentum deficit. The test objective was to quantify noise reductions, measure impacts on fan aerodynamic performance, and document the flow field using hot-film anemometry. Measurements concentrated on approach, cutback, and takeoff rotational speeds as those are the primary conditions of acoustic interest. Data are presented for a 2% (relative to overall fan flow) trailing edge injection rate and show a 2 dB reduction in Overall Sound Power Level (OAPWL) at all fan test speeds. The reduction in broadband noise is nearly constant and is approximately 1.5 dB up to 20 kHz at all fan speeds. Measurements of tone noise show significant variation, as evidenced by reductions of up to 6 dB in the 2 BPF tone at 6700 rpm.: and increases of nearly 2 dB for the 4 BPF tone at approach speed. Aerodynamic performance measurements show the fan with 2 % injection has an overall efficiency that is comparable to the baseline fan and operates, as intended, with nearly the same pressure ratio and mass flow parameters. Hot-film measurements obtained at the approach operating condition indicate that mean blade wake filling in the tip region was not as significant as expected. This suggests that additional acoustic benefits could be realized if the trailing edge blowing could be modified to provide better filling of the wake momentum deficit. Nevertheless, the hot-film measurements indicate that the trailing edge blowing provided significant reductions in blade wake turbulence. Overall, these results indicate that further work may be required to fully understand the proper implementation of injecting flow at/near the trailing edge as a wake filling strategy. However, data do support the notion that noise reductions can be realized not only for tones but perhaps more importantly, also for broadband. Furthermore, the technique can be implemented without adversely effecting overall fan aerodynamic performance.

Fite, E. Brian↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

New Global Electron Density Observations from GPS-RO in the D- and E-Region Ionosphere

A novel retrieval technique is developed for electron density (N(sub e)) in the D- and E-region (80-120 km) using the high-quality 50-Hz GPS radio occultation (GPS-RO) phase measurements. The new algorithm assumes a slow, linear variation in the F-region background when the GPS-RO passes through the D- and E-region, and extracts the N(sub e) profiles at 80-130 km from the phase advance signal caused by N(sub e). Unlike the conventional Abel function, the new approach produces a sharp N(sub e) weighting function in the lower ionosphere, and the N(sub e) retrievals are in good agreement with the IRI (International Reference Ionosphere) model in terms of monthly maps, zonal means and diurnal variations. The daytime GPS-RO N(sub e) profiles can be well characterized by the alpha-Chapman function of three parameters (N(sub mE), h(sub mE) and H), showing that the bottom of E-region is deepening and sharpening towards the summer pole. At high latitudes the monthly GPS-RO N(sub e) maps at 80-120 km reveal clear enhancement in the auroral zones, more prominent at night, as a result of energetic electron precipitation (EEP) from the outer radiation belt. The D-/E-region auroral N(sub e) is strongly correlated with K(sub p) on a daily basis. The new N(sub e) data allow further comprehensive analyses of the sporadic E (E(sub s)) phenomena in connection with the background N(sub e) in the E-region. The layered (2-10 km) and fluctuated (less than 2 km) E(sub s) components, namely N(sub e_Layer) than N(sub e_Pert), are extracted with respect to the background N( sub e_Region) on a profile-by-profile basis. The N(sub e_Layer) component has a strong but highly-refined peak at approximately 105 km, with an amplitude smaller than N(sub e_Region) approximately by an order of magnitude. The N(sub e_Pert) component, which was studied extensively in the past, is approximately 2 orders of magnitude weaker than N(sub e_Layer). Both N(sub e_Layer) and N(sub e_Pert) are subject to significant diurnal and semidiurnal variations, showing downward progression with local time in amplitude. The 11-year solar cycle dominates the N(sub e) interannual variations, showing larger N(sub e_Region) and N(sub e_Layer) but smaller N(sub e_Pert) amplitudes in the solar maximum years. Enhanced Ne profiles are often observed in the polar winter, showing good correlation with solar proton events (SPEs) and geomagnetic activity. The new methodology offers great potential for retrieving low N(sub e) in the D-region, where radio propagation and communication blackouts can occur due to enhanced ionization. For space weather applications it is recommended for GPSRO operations to raise the top of high-rate data acquisition to approximately 140 km in the future.

GSP-R↗

Computing Highly Accurate Spectroscopic Line Lists for Characterization of Planetary Atmospheres: CO2 and SO2 Line Lists Needed for Modeling Venus

Over the last decade, it has become apparent that the most effective approach for determining highly accurate rotational and rovibrational line lists for molecules of interest in planetary atmospheres and other astrophysical environments is through a combination of highresolution laboratory experiments coupled with state-of-the art ab initio quantum chemistry methods. The approach involves computing the most accurate potential energy surface (PES) possible using state-of-the art electronic structure methods, followed by computing rotational and rovibrational energy levels using an exact variational method to solve the nuclear Schrödinger equation. Then, reliable experimental data from high-resolution experiments is used to refine the ab initio PES in order to improve the accuracy of the computed energy levels and transition energies. From the refinement step, we have been able to achieve an accuracy of approximately 0.015 cm-1 for rovibrational transition energies, and even better for purely rotational transitions. This combined "experiment / theory" approach allows for determination of essentially a complete line list, with hundreds of millions of transitions, and having the transition energies and intensities be highly accurate. Our group has successfully applied this approach to determine highly accurate line lists for NH3, CO2 and isotopologues, and SO2 and isotopologues. Here I will report our latest results for CO2 and SO2 including all isotopologues. Comparisons to the available data in HITRAN2012 and other available databases will be shown, though we note that our line lists for SO2 are significantly more complete than any other databases. Since it is important to span a large temperature range in order to model the spectral signature of Venus as well as exoplanets, we will demonstrate how the spectra change on going from low temperatures (100 K) to higher temperatures (500 K to 1500 K).

potential energy surface (PES)5405 PLANETARY SCIEN↗

Multi Model Monte Carlo with Python (MXMCPy)

Multi Model Monte Carlo with Python (\mxmc {}) is a software package developed as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Motivated by uncertainty propagation problems where classical Monte Carlo (MC) simulation is computationally intractable, various multi-model MC approaches have recently emerged that yield unbiased estimators with significantly reduced variance relative to MC for the same cost. These existing methods include multi-level Monte Carlo (MLMC), multi-fidelity Monte Carlo (MFMC), and approximate control variates (ACV). Given a fixed computational budget and a collection of models with varying cost/accuracy, each method seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. \mxmc {} is a versatile tool that enables convenient access to many existing multi-model MC approaches within one modular and extensible package. With \mxmc {}, users can easily compare existing methods to determine the best choice for their particular problem, while developers have a basis for implementing and sharing new variance reduction approaches. This report introduces the \mxmc {} software, providing a summary of the problem-solving workflow for users as well as a brief overview of the code layout for developers.

Geoffrey F Bomarito↗

A topological approach to computer-aided sensitivity analysis

Sensitivities of any arbitrary system are calculated using general purpose digital computer with available software packages for transfer function analysis. Sensitivity shows how element variation within system affects system performance. Signal flow graph illustrates topological system behavior and relationship among parameters in system.

Chan, S. P.↗

Assimilation of Geosat Altimetric Data in a Nonlinear Shallow-Water Model of the Indian Ocean by Adjoint Approach. Part II: Some Validation and Interpretation of the Assimilated Results: Some Validation and Interpretation of the Assimilated Results - Part 2

This paper examines the results of assimilating Geosat sea level variations relative to the November 1986-November 1988 mean reference, in a nonlinear reduced-gravity model of the Indian Ocean, Data have been assimilated during one year starting in November 1986 with the objective of optimizing the initial conditions and the yearly averaged reference surface. The thermocline slope simulated by the model with or without assimilation is validated by comparison with the signal, which can be derived from expandable bathythermograph measurements performed in the Indian Ocean at that time. The topography simulated with assimilation on November 1986 is in very good agreement with the hydrographic data. The slopes corresponding to the South Equatorial Current and to the South Equatorial Countercurrent are better reproduced with assimilation than without during the first nine months. The whole circulation of the cyclonic gyre south of the equator is then strongly intensified by assimilation. Another assimilation experiment is run over the following year starting in November 1987. The difference between the two yearly mean surfaces simulated with assimilation is in excellent agreement with Geosat. In the southeastern Indian Ocean, the correction to the yearly mean dynamic topography due to assimilation over the second year is negatively correlated to the one the year before. This correction is also in agreement with hydrographic data. It is likely that the signal corrected by assimilation is not only due to wind error, because simulations driven by various wind forcings present the same features over the two years. Model simulations run with a prescribed throughflow transport anomaly indicate that assimilation is rather correcting in the interior of the model domain for inadequate boundary conditions with the Pacific.

Greiner, Eric↗

Use of Smoothed Measured Winds to Predict and Assess Launch Environments

This paper describes an approach which was used to predict Ares I-X launch vehicle structural airloads and controllability prior to launch. This prediction along with the proper application of protection for variations in key environment and trajectory parameters was used to predict the potential for the launch trajectory to stay within the integrated vehicle certification envelopes. This data was applied by the launch team as a key element of the launch day go/no-go recommendation. NASA s Space Shuttle and the heavy lift versions of the Evolved Expendable Launch Vehicles assessments require the use of a high fidelity wind measurement which limits the number of wind measurement sources which can be used for a pre-launch go/no-go assessment. The approach described in this paper has the potential to allow the pre-launch assessment team to use wind measurements from a number of additional sources which will enhance launch availability and also has the potential to improve mission assurance and enhance safety

Cordova, Henry S.↗

Usefulness of heart measures in flight simulation

The results of three studies performed at the NASA Langley Research Center are presented to indicate the areas in which heart measures are useful for detecting differences in the workload state of subjects. Tasks that involve the arousal of the sympathetic nervous system, such as landing approaches, were excellent candidates for the use of average heart-rate and/or the increase in heart-rate during a task. The latter of these two measures was the better parameter because it removed the effects of diurnal variations in heart-rate and some of the intersubject variability. Tasks which differ in the amount of mental resources required are excellent candidates for heart-rate variability measures. Heart-rate variability measures based upon power spectral density techniques were responsive to the changing task demands of landing approach tasks, approach guidance options, and 2 versus 20 second interstimulus-intervals of a monitoring task. Heart-rate variability measures were especially sensitive to time-on-task when the task was characterized by minimal novelty, complexity, and uncertainty (i.e., heart-rate variability increases as a function of the subjects boredom).

Harris, Randall L., Sr.↗

Space shuttle orbit maneuvering engine reusable thrust chamber program

Tests were conducted on the regenerative cooled thrust chamber of the space shuttle orbit maneuvering engine. The conditions for the tests and the durations obtained are presented. The tests demonstrated thrust chamber operation over the nominal ranges of chamber pressure mixture ratio. Variations in auxiliary film coolant flowrate were also demonstrated. High pressure tests were conducted to demonstrate the thrust chamber operation at conditions approaching the design chamber pressure for the derivative space tug application.

Pauckert, R. P.↗

Computation of three-dimensional nozzle-exhaust flow fields with the GIM code

A methodology is introduced for constructing numerical analogs of the partial differential equations of continuum mechanics. A general formulation is provided which permits classical finite element and many of the finite difference methods to be derived directly. The approach, termed the General Interpolants Method (GIM), can combined the best features of finite element and finite difference methods. A quasi-variational procedure is used to formulate the element equations, to introduce boundary conditions into the method and to provide a natural assembly sequence. A derivation is given in terms of general interpolation functions from this procedure. Example computations for transonic and supersonic flows in two and three dimensions are given to illustrate the utility of GIM. A three-dimensional nozzle-exhaust flow field is solved including interaction with the freestream and a coupled treatment of the shear layer. Potential applications of the GIM code to a variety of computational fluid dynamics problems is then discussed in terms of existing capability or by extension of the methodology.

Spradley, L. W.↗

Cometary nucleus and active regions

On the basis of the icy conglomerate model of cometary nuclei, various observations demonstrate the spotted nature of many or most nuclei, i.e., regions of unusual activity, either high or low. Rotation periods, spin axes and even precession of the axes are determined. The observational evidence for variations in activity over the surfaces of cometary nuclei are listed and discussed. On June 11 the comet IRAS-ARAKI-ALCOCK approached the Earth to a distance of 0.031 AU, the nearest since C/Lexell, 1770 I, providing a unique opportunity for near-nucleus observations. Preliminary analysis of these images establishes the spin axis of the nucleus, with an oblioquity to the orbit plane of approximately 50 deg, and a lag angle of sublimation approximately 35 deg from the solar meridian on the nucleus. Asymmetries of the inner coma suggests a crazy-quilt distribution of ices with differing volatility over the surface of the nucleus. The observations of Comet P/Homes 1892 III, exhibiting two 8-10 magnitude bursts, are carefully analyzed. The grazing encounter produced, besides the first great burst, an active area on the nucleus, which was rotating retrograde with a period of 16.3hr and inclination nearly 180 deg. After the first burst the total magnitude fell less than two magnitudes from November 7 to November 30 (barely naked eye) while the nuclear region remained diffuse or complex, rarely if ever showing a stellar appearance. The fading was much more rapid after the second burst. The grazing encounter distributed a volume of large chunks in the neighborhood of the nucleus, maintaining activity for weeks.

Whipple, F. L.↗

Aboveground Biomass Estimation in a Tidal Brackish Marsh Using Simulated Thematic Mapper Spectral Data

Spectral radiance data were collected from the ground and from a low altitude aircraft in an attempt to gain some insight into the potential utility of actual Thematic Mapper data for biomass estimation in wetland plant communities. No attempt was made to distinguish individual plant species within brackish marsh plant associations. Rather, it was decided to lump plant species with similar canopy morphologies and then estimate from spectral radiance data the biomass of the group. The rationale for such an approach is that plants with a similar morphology will produce a similar reflecting or absorping surface (i.e., canopy) for incoming electromagnetic radiation. Variations in observed reflectance from different plant communities with a similar canopy morphology are more likely to be a result of biomass differences than a result of differences in canopy architecture. If the hypothesis that plants with a similar morphology exhibit similar reflectance characteristics is true, then biomass can be estimated based on a model for the dominant plant morphology within a plant association and the need for species discrimination has effectively been eliminated.

Hardisky, M.↗

Satellite passive microwave rain measurement techniques for land and ocean

Multiseasonal rainfall was found to be measurable over land with satellite passive microwave data, based upon comparisons between Nimbus 7 Scanning Multichannel Microwave Radiometer (SMME) brightness temperatures (T sub B) and operational WSR-57 radar rain rates. All of the SMMR channels (bipolarized 37, 21, 18, 10.7, and 6.6. GHz T sub B) were compared to radar reflectivities for 25 SMMR passes and 234 radar scans over the U.S. during the spring, summer, and fall of 1979. It was found that the radar rain rates were closely related to the difference between 37 and 21 GHz T sub B. This result is due to the volume scattering effects of precipitation which cause emissivity decreases with frequency, as opposed to emissive surfaces (e.g., water) whose emissivities increase with frequency. Two frequencies also act to reduce the effects of thermometric temperature variations on T sub B to a miminum. During summer and fall, multiple correlation coefficients of 0.80 and 0.75 were obtained. These approach the limit of correlation that can be expected to exist between two very different data sources, especially in light of the errors attributable to manual digitization of PPI photographs of variable quality from various operational weather radar not calibrated for research purposes. During the spring, a significantly lower (0.63) correlation was found. This poorer performance was traced to cases of wet, unvegetated soil being sensed at the lower frequencies through light rain, partly negating the rain scattering signal.

Spencer, R. W.↗

Extension of a semi-implicit shock-capturing algorithm for 3-D fully coupled, chemically reacting flows in generalized coordinates

A semiimplicit high-resolution shock-capturing method for multidimensional systems of hyperbolic conservation laws with stiff source terms has been developed by Yee and Shinn (1987). The goal of this work is to extend this method to solve the three-dimensional fully coupled Navier-Stokes equations for a hypersonic chemically reacting flow in generalized coordinates. In this formulation, the global continuity equation was replaced by all the species continuity equations. The shock-capturing technique is a second-order-accurate, symmetric total-variation-diminishing method which accounts fully and directly for the coupling among the fluid and all the species. To verify the current approach, it was implemented into an existing computer code which contained the MacCormack method. Test results for a five-species reacting flow are shown to be oscillation-free around the shock, and the time spent per iteration only doubles when compared to the result using classical way of supplying numerical dissipation. The extra computation is more than justified by the elimination of spurious oscillation and nonlinear instability associated with the classical shock-capturing schemes in computing hypersonic reacting flows.

Shinn, J. L.↗

Stall Characteristics Obtained from Flight 10 of Northrop X-4 No. 2 Airplane (USAF No. 46-677)

NACA instrumentation has been installed ii the X-J4 airplanes to obtain stability and control data during the acceptance tests conducted by the Northrop Aircraft Corporation. This report presents data obtained on the stalling characteristics of the airplane in the clean and gear- down configurations. The center of gravity was located at approximately 18 percent of the mean aerodynamic chord during the tests. The results indicated that the airplane was not completely stalled when stall was gradually approached during nominally U accelerated flight but that it was completely stalled during a more abruptly approached stall in accelerated flight. The stall in accelerated flight was relatively mild, and this was attributed to the nature of the variation of lift with angle of attack for the 001-614 airfoil section, the plan form of the wing, and to the fact that the initial sideslip at the stall produced (as shown by wind-tunnel tests of a model of the airplane) a more symmetrical stall pattern.

Sadoff, Melvin↗

Electro-Optically Active Monomers: Synthesis and Characterization of Thin Films of Liquid Crystalline Substituted Polyacetylenes

The overall objective of this study was the description of the behavior of mesogen substituted acetylene monomers and polymers in monolayer films at the air/water interface and as multilayer films including the formation of such films. Fundamental knowledge to be gained would include the effect of balancing hydrophilic and hydrophobic tendencies in a molecule more complex than the classical fatty acids or lipids. The effect of molecular shape on the packing and thus the ultimate stability of monolayers formed from these new molecules was explored. The work takes on the challenge of preorienting monomers in well-ordered arrays prior to attempting polymerization with the hope that order would be preserved in any resulting polymer. New knowledge gained with regard to the acetylenic monomers includes processing of the acetylene monomer into multi-layer films, followed by the design and synthesis of a second generation of improved monomer structure for superior LBK film transfer properties. A third generation of acetylenic monomer was synthesized which approaches more closely the goal of solid state polymerization of these materials. A parallel study took a different approach. The materials are pre-formed poly(phenylene-acetylene) polymers so questions about reactivity are mute. The materials are a variation on the well-known hairy-rod polymers with regard to their Langmuir film-forming properties. Overall, the goal was to demonstrate that these polymers could be processed into NLO materials with novel polar order.

Duran, R. S.↗