Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 235 records · Page 13

Experimental investigation of helicopter vibration reduction using rotor blade aeroelastic tailoring

A wind tunnel investigation has been conducted to parametrically investigate the effect of blade nonstructural mass on helicopter fixed- and rotating-system vibratory loads. The data were obtained using Mach- and aeroelastically-scaled model rotor blades which allowed for the addition of concentrated nonstructural masses at multiple points along the blade radius. Testing was conducted for advance ratios ranging from 0.10 to 0.35 for ten blade mass configurations. Three thrust levels were obtained at representative full-scale shaft angles for each configuration. Results indicate that proper placement of blade nonstructural mass can provide reductions in fixed-system vibratory loads, but that correct mass placement and the loads reduction realized are dependent upon flight condition. The data base obtained with this investigation provides a comprehensive set of fixed-system shears and moments, blade moment, and blade flap and lag response. The data set is well suited for use in the correlation and development of advanced rotorcraft analyses.

Wilbur, Matthew L.↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

The Influence of Dynamic Contact Angle on Wetting Dynamics

When surface tension forces dominate, and regardless of whether the situation is static or dynamic, the contact angle (the angle the interface between two immiscible fluids makes when it contacts a solid) is the key parameter that determines the shape of a fluid-fluid interface. The static contact angle is easy to measure and implement in models predicting static capillary surface shapes and such associated quantities as pressure drops. By contrast, when the interface moves relative to the solid (as in dynamic wetting processes) the dynamic contact angle is not identified unambiguously because it depends on the geometry of the system Consequently, its determination becomes problematic and measurements in one geometry cannot be applied in another for prediction purposes. However, knowing how to measure and use the dynamic contact angle is crucial to determine such dynamics as a microsystem throughput reliably. In this talk we will present experimental and analytical efforts aimed at resolving modeling issues present in dynamic wetting. We will review experiments that show the inadequacy of the usual hydrodynamic model when a fluid-fluid meniscus moves over a solid surface such as the wall of a small tube or duct. We will then present analytical results that show how to parametrize these problems in a predictive manner. We will illustrate these ideas by showing how to implement the method in numerical fluid mechanical calculations.

Rame, Enrique↗

Kinetic Theory and Fast Wind Observations of the Electron Strahl

We develop a model for the strahl population in the solar wind - a narrow, low-density and high-energy electron beam centred on the magnetic field direction. Our model is based on the solution of the electron drift-kinetic equation at heliospheric distances where the plasma density, temperature and the magnetic field strength decline as power laws of the distance along a magnetic flux tube. Our solution for the strahl depends on a number of parameters that, in the absence of the analytic solution for the full electron velocity distribution function (eVDF), cannot be derived from the theory. We however demonstrate that these parameters can be efficiently found from matching our solution with observations of the eVDF made by the Wind satellite's SWE strahl detector. The model is successful at predicting the angular width (FWHM) of the strahl for the Wind data at 1 au, in particular by predicting how this width scales with particle energy and background density. We find that the strahl distribution is largely determined by the local temperature Knudsen number γ ∼ |T dT/dx|/n, which parametrizes solar wind collisionality. We compute averaged strahl distributions for typical Knudsen numbers observed in the solar wind, and fit our model to these data. The model can be matched quite closely to the eVDFs at 1 au; however, it then overestimates the strahl amplitude at larger heliocentric distances. This indicates that our model may be improved through the inclusion of additional physics, possibly through the introduction of 'anomalous diffusion' of the strahl electrons.

Horaites, Konstantinos↗

Mach number scaling of helicopter rotor blade/vortex interaction noise

A parametric study of model helicopter rotor blade slap due to blade vortex interaction (BVI) was conducted in a 5 by 7.5-foot anechoic wind tunnel using model helicopter rotors with two, three, and four blades. The results were compared with a previously developed Mach number scaling theory. Three- and four-bladed rotor configurations were found to show very good agreement with the Mach number to the sixth power law for all conditions tested. A reduction of conditions for which BVI blade slap is detected was observed for three-bladed rotors when compared to the two-bladed baseline. The advance ratio boundaries of the four-bladed rotor exhibited an angular dependence not present for the two-bladed configuration. The upper limits for the advance ratio boundaries of the four-bladed rotors increased with increasing rotational speed.

Leighton, Kenneth P.↗

Software for Estimating Costs of Testing Rocket Engines

A high-level parametric mathematical model for estimating the costs of testing rocket engines and components at Stennis Space Center has been implemented as a Microsoft Excel program that generates multiple spreadsheets. The model and the program are both denoted, simply, the Cost Estimating Model (CEM). The inputs to the CEM are the parameters that describe particular tests, including test types (component or engine test), numbers and duration of tests, thrust levels, and other parameters. The CEM estimates anticipated total project costs for a specific test. Estimates are broken down into testing categories based on a work-breakdown structure and a cost-element structure. A notable historical assumption incorporated into the CEM is that total labor times depend mainly on thrust levels. As a result of a recent modification of the CEM to increase the accuracy of predicted labor times, the dependence of labor time on thrust level is now embodied in third- and fourth-order polynomials.

Hines, Merion M.↗

Software for Estimating Costs of Testing Rocket Engines

A high-level parametric mathematical model for estimating the costs of testing rocket engines and components at Stennis Space Center has been implemented as a Microsoft Excel program that generates multiple spreadsheets. The model and the program are both denoted, simply, the Cost Estimating Model (CEM). The inputs to the CEM are the parameters that describe particular tests, including test types (component or engine test), numbers and duration of tests, thrust levels, and other parameters. The CEM estimates anticipated total project costs for a specific test. Estimates are broken down into testing categories based on a work-breakdown structure and a cost-element structure. A notable historical assumption incorporated into the CEM is that total labor times depend mainly on thrust levels. As a result of a recent modification of the CEM to increase the accuracy of predicted labor times, the dependence of labor time on thrust level is now embodied in third- and fourth-order polynomials.

Hines, Merlon M.↗

Software for Estimating Costs of Testing Rocket Engines

A high-level parametric mathematical model for estimating the costs of testing rocket engines and components at Stennis Space Center has been implemented as a Microsoft Excel program that generates multiple spreadsheets. The model and the program are both denoted, simply, the Cost Estimating Model (CEM). The inputs to the CEM are the parameters that describe particular tests, including test types (component or engine test), numbers and duration of tests, thrust levels, and other parameters. The CEM estimates anticipated total project costs for a specific test. Estimates are broken down into testing categories based on a work-breakdown structure and a cost-element structure. A notable historical assumption incorporated into the CEM is that total labor times depend mainly on thrust levels. As a result of a recent modification of the CEM to increase the accuracy of predicted labor times, the dependence of labor time on thrust level is now embodied in third- and fourth-order polynomials.

Hines, Merlon M.↗

Bayesian Inference and the Effects of Varying Uncertainty Models in Charring Ablator Calibration and Uncertainty Quantification Problems

The Mars Science Laboratory (MSL) vehicle utilized a heat shield constructed from NASA’s Phenolic-Impregnated Carbon Ablator (PICA) material to protect the main structure from the high enthalpy environment encountered during hypersonic atmospheric entry. During the vehicle’s descent through Martian atmosphere, multiple thermocouples embedded within the heat shield captured in-depth material temperature data that allow for studies to be conducted on current material response reconstruction tools. In the present work, material temperature data obtained from thermocouples within the MISP-4 plug (MEDLI (Mars Science Laboratory Entry, Descent, and Landing Instrument) Integrated Sensor Plug) are utilized in the calibration of Theoretical Ablative Composite for Open Testing (TACOT) model parameters in conjunction with NASA’s Porous material Analysis Toolbox (PATO) through Bayesian inference where uncertainty due to parametric, modeling, and experimental sources is simultaneously quantified. Prior to the study, a sensitivity analysis is performed through computation of the robust Sobol indices in an effort to study the relationship between input space and model response and to reduce the dimensionality of the statistical inverse problem. The Bayesian inference methodology necessitates an a-priori choice to be made for the uncertainty model for which numerous possibilities are available. Across most works, however, only basic additive or multiplicative models are utilized with pre-defined magnitudes of uncertainty based on a-priori knowledge or to-be-calibrated multipliers of static covariance matrix structures. The present effort explores the effects of informed uncertainty models, ones with temporal dependence that are simultaneously calibrated through Bayesian inference, on calibrated results for parameters that make up the uncertain input space.

Sensitivity Analysis↗

Investigation of CAST-10-2/DOA 2 airfoil in NAE high Reynolds number two-dimensional test facility

A common airfoil model with the CAST 10-2/DOA-2 profile and 228 mm (9 inches) chord length was tested. The tests performed in NAE covered the Mach numbers from 0.3 to 0.8 and Reynolds numbers from 10 to 30 million. The model was tested with transition free and with transition fixed at 5 percent chord for both the upper and the lower surfaces. The data obtained were analyzed for the effects of Reynolds number, transition fixing and Mach number. The role of the boundary layer on the displacement effect, the interaction with the shock wave and the trailing edge separation are examined. The results are summarized as follows: (1) the airfoil performance depends strongly on Reynolds number and transition fixing; (2) with transition fixed, the aerodynamic quantities such as lift, pitching moment and drag show a monotonic variation with Reynolds number; (3) with transition free, the aerodynamic quantities vary less regularly with Reynolds number and a slight parametric dependency is shown. The weak dependency is due to the compensatory effect of the forward shift of the transition position and the thinning of the turbulent boundary layer as Reynolds number increases; (4) the shock Mach number and the shock position are weakly dependent on Reynolds number; and (5) the long extent of the laminar boundary layer at transonic speeds reduces the drag appreciably at low Reynolds numbers. The drag bucket around the design Mach number can be observed below Reynolds number 15 million.

Chan, Y. Y.↗

Towards Simultaneous Multi-Parameter Measurements Using Burst-Mode Frequency-Scanning Planar Laser-Induced Fluorescence

While planar laser-induced fluorescence (PLIF) can provide valuable flow visualization, there are many challenges in extracting quantitative measurements from the fluorescence intensity. To overcome the multi-parameter dependence of the signal intensity, frequency scanning planar laser-induced fluorescence (FS-PLIF) provides the potential for 2D temperature, pressure, and velocity measurements by spectrally resolving the excitation spectrum of one or more fluorescence transitions to extract the Doppler shift, Doppler- and collisional-broadening of absorption lines, and the relative signal intensities of multiple transitions. This work aims to characterize a burst-mode FS-PLIF system consisting of a Spectral Energies, LLC burst-mode laser coupled with a Spectral Energies, LLC optical parametric oscillator using nitric oxide (NO) static pressure cell measurements. The effect of the OPO cavity on intensity fluctuations during frequency-scanning is investigated. A NO PLIF model is implemented to investigate the feasibility of property measurements using single- and two-NO-line scans.

burst-mode laser↗

Hypervelocity Aeroheating Measurements in Wake of Mars Mission Entry Vehicle

Detailed measurements of aerodynamic heating rates in the wake of a Mars-Pathfinder configuration model have been made. Heating data were obtained in a conventional wind tunnel, the NASA LaRC 31" Mach 10 Air Tunnel, and in a high-enthalpy impulse facility, the NASA HYPULSE expansion tube, in which air and CO2 were employed as test gases. The enthalpy levels were 0.7 MJ/kg in the Mach 10 Tunnel, 12 MJ/kg at Mach 9.8 for HYPULSE CO2 tests and 14 MJ/kg at Mach 7.9 for HYPULSE air tests. Wake heating rates were also measured on three similar parametric configurations, and forebody heating measurements were made in order to facilitate CFD comparisons. The ratio of peak wake heating to forebody stagnation point heating in the Mach 10 Tunnel varied from 7% to 15% depending on the freestream Reynolds number. In HYPULSE, the ratio was ~5% for both air and CO 2. It was observed that an increase in the ratio of forebody corner radius to nose radius resulted in a decrease in peak wake heating, and moved the peak closer to the base of the forebody. The wake flow establishment process in HYPULSE was studied, and a method was developed to determine when the wake has become fully established.

Hollis, Brian R.↗

Orbiter Capability for Providing Water to the International Space Station according to the Most Probable Flight Attitudes

Water to be generated by, delivered to, and processed by the International Space Station (ISS) is a critical Environmental Control and Life Support (ECLS) element, especially for the early ISS missions. A significant portion of the water required by the ISS shall be provided by the Shuttle Transportation System (STS) Orbiter. The balance of water generated by the Orbiter Fuel Cells (FC), minus that water consumed by the Orbiter Flash Evaporator System (FES) and crew, is available for transfer to the ISS. During later missions, crew respired and perspired water, as well as effluent water from the Orbiter LiOH canisters, will be collected as condensate and available for transfer to the ISS. Orbiter radiator performance provides the most variance in determining the amount of net Orbiter water available for transfer to the ISS. As radiator performance decreases, the dependence upon the FES (and FC water) increases for rejecting Orbiter waste heat. Generally, radiator performance decreases as the ISS assembly size increases (especially as solar arrays are added), and also as beta angle increases. ISS solar array deployment necessitates the use of models with articulating solar arrays (for Earth local-vertical attitudes), as array position dramatically affects Orbiter radiator performance. Recent developments in the relaxation of beta angle limitations have also increased the complexity and difficulty of providing water to the ISS. Other factors that may hinder the ability to transfer water are the number of empty Contingency Water Containers (CWCs) available, duration of open-hatch time, crew activity timeline, and full CWC storage capability. A parametric study has been accomplished that provides a quick-reference table for determining expected water generation rates for ISS missions 2A.2 through 7A.1. An hourly Orbiter water generation rate is reported according to a matrix that consists of: (1) (six) significant changes in ISS assembly configuration; (2) (four) beta angles (0 deg. , +37 deg., +53 deg. , and +75 deg.); (3) the (three) most representative ISS attitudes (XPOP-O, XPOP-180 and +XVV); (4) (four) Orbiter radiator configurations (both stowed, starboard deployed, port deployed, and both deployed) and (5) the (two) conditions (radiator inlet temperatures and fuel cell power) most consistent with sleep and wake periods. Those permutations of higher probability of occurrence than others have been identified. Another parametric study has been accomplished that provides a quick-reference table for determining expected water generation rates for ISS assembly complete missions. An hourly Orbiter water generation rate is reported according to a matrix that consists of: (1) (seven) beta angles (-75 deg., -60 deg., -30 deg., 0 deg., +30 deg., +60 deg., and +75 deg.); (2) the (nine) PYR angles that define the corners of the envelope; (3) (four) Orbiter radiator configurations (both stowed, starboard deployed, port deployed, and both deployed) and (4) the (two) conditions (radiator inlet temperatures and fuel cell power) most consistent with sleep and wake periods.

Dunaway, Brian↗

Uranian methane abundance, rotational temperature, and effective pressure from the 6800 A band

Measurements of the methane absorption features at 6800 A in two spectrograms of Uranus have been analyzed on the working hypothesis that the features are the R branch of the 5 mu(3) rotation-vibration overtone. Internal and external consistencies in temperature calculations from the measurements give increased confidence that the band identification and rotational-quantum-number assignments to band members are correct. Parametric fitting of synthetic reflecting-layer and homogeneous scattering-layer spectra to the spectrograms is attempted. The best-fit methane abundances, taken with recent estimates of the molecular hydrogen abundance, correspond to an average C:H ratio between 1.5 and 25 times larger than the solar ratio. Even greater enrichment of the C:H ratio may occur below the visible cloud boundary. High spectroscopic resolution of the present data permits direct measurement of the methane line halfwidth in Uranus's spectrum. The halfwidth corresponds to an effective pressure greater than 3 atm, depending on the mixture of atmospheric gases. Pressure in this range, at the methane rotational temperature of between 93 and 100 K, is inconsistent with current Uranus atmosphere models having solar elemental-abundance ratios. It is consistent with models which are considerably enriched in 'heavier' constituents.

Bergstralh, J. T.↗

Computational Analysis for Rocket-Based Combined-Cycle Systems During Rocket-Only Operation

A series of Reynolds-averaged Navier-Stokes calculations were employed to study the performance of rocket-based combined-cycle systems operating in an all-rocket mode. This parametric series of calculations were executed within a statistical framework, commonly known as design of experiments. The parametric design space included four geometric and two flowfield variables set at three levels each, for a total of 729 possible combinations. A D-optimal design strategy was selected. It required that only 36 separate computational fluid dynamics (CFD) solutions be performed to develop a full response surface model, which quantified the linear, bilinear, and curvilinear effects of the six experimental variables. The axisymmetric, Reynolds-averaged Navier-Stokes simulations were executed with the NPARC v3.0 code. The response used in the statistical analysis was created from Isp efficiency data integrated from the 36 CFD simulations. The influence of turbulence modeling was analyzed by using both one- and two-equation models. Careful attention was also given to quantify the influence of mesh dependence, iterative convergence, and artificial viscosity upon the resulting statistical model. Thirteen statistically significant effects were observed to have an influence on rocket-based combined-cycle nozzle performance. It was apparent that the free-expansion process, directly downstream of the rocket nozzle, can influence the Isp efficiency. Numerical schlieren images and particle traces have been used to further understand the physical phenomena behind several of the statistically significant results.

Steffen, C. J., Jr.↗

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.

Venkateswaran,S.↗

Contamination of grazing incidence EUV mirrors - An assessment

Contamination assessment for space optical systems requires an understanding of the sensitivity of component performance, e.g. mirror reflectance, to materials deposited on the mirror surface. In a previous study, the sensitivity of typical normal incidence mirror coatings to surface deposits of generic hydrocarbons was reported. Recent activity in the development of grazing incidence telescopes for extreme ultraviolet space astronomy has stimulated the need for a similar assessment in the spectral region extending from approximately 100 A to 1000 A. The model used for analysis treats the contamination layer as a continuous thin film deposited on the mirror surface. The mirror surfaces selected for this study are opaque vacuum deposited gold and the uncoated and polished Zerodur. Scatter caused by film irregularities or particulates are not included in this assessment. Parametric evaluations at 100, 500, and 1000 A determine the sensitivity of mirror reflectance to a range of optical constants selected for the generic contaminants. This sensitivity analysis combined with the limited amount of optical data in the EUV for hydrocarbons, is used to select representative optical constants for the three wavelength regions. Reflectance versus contamination layer thickness curves are then calculated and used to determine critical thickness limits based on allowable reflectance change. Initial observations indicate that thickness limits will be highly dependent on the real part of the complex index of refraction of the contaminant film being less than 1.0. Preliminary laboratory measurements of samples contaminated with some commonly encountered hydrocarbons confirm trends indicated in the analytical studies.

Osantowski, John F.↗

An Acoustic Charge Transport Imager for High Definition Television Applications: Reliability Modeling and Parametric Yield Prediction of GaAs Multiple Quantum Well Avalanche Photodiodes

Reliability modeling and parametric yield prediction of GaAs/AlGaAs multiple quantum well (MQW) avalanche photodiodes (APDs), which are of interest as an ultra-low noise image capture mechanism for high definition systems, have been investigated. First, the effect of various doping methods on the reliability of GaAs/AlGaAs multiple quantum well (MQW) avalanche photodiode (APD) structures fabricated by molecular beam epitaxy is investigated. Reliability is examined by accelerated life tests by monitoring dark current and breakdown voltage. Median device lifetime and the activation energy of the degradation mechanism are computed for undoped, doped-barrier, and doped-well APD structures. Lifetimes for each device structure are examined via a statistically designed experiment. Analysis of variance shows that dark-current is affected primarily by device diameter, temperature and stressing time, and breakdown voltage depends on the diameter, stressing time and APD type. It is concluded that the undoped APD has the highest reliability, followed by the doped well and doped barrier devices, respectively. To determine the source of the degradation mechanism for each device structure, failure analysis using the electron-beam induced current method is performed. This analysis reveals some degree of device degradation caused by ionic impurities in the passivation layer, and energy-dispersive spectrometry subsequently verified the presence of ionic sodium as the primary contaminant. However, since all device structures are similarly passivated, sodium contamination alone does not account for the observed variation between the differently doped APDs. This effect is explained by the dopant migration during stressing, which is verified by free carrier concentration measurements using the capacitance-voltage technique.

Hunt, W. D.↗