Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservative”

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,117 records · Page 62

Icing Analysis of a Swept NACA 0012 Wing Using LEWICE3D Version 3.48

Icing calculations were performed for a NACA 0012 swept wing tip using LEWICE3D Version 3.48 coupled with the ANSYS CFX flow solver. The calculated ice shapes were compared to experimental data generated in the NASA Glenn Icing Research Tunnel (IRT). The IRT tests were designed to test the performance of the LEWICE3D ice void density model which was developed to improve the prediction of swept wing ice shapes. Icing tests were performed for a range of temperatures at two different droplet inertia parameters and two different sweep angles. The predicted mass agreed well with the experiment with an average difference of 12%. The LEWICE3D ice void density model under-predicted void density by an average of 30% for the large inertia parameter cases and by 63% for the small inertia parameter cases. This under-prediction in void density resulted in an over-prediction of ice area by an average of 115%. The LEWICE3D ice void density model produced a larger average area difference with experiment than the standard LEWICE density model, which doesn't account for the voids in the swept wing ice shape, (115% and 75% respectively) but it produced ice shapes which were deemed more appropriate because they were conservative (larger than experiment). Major contributors to the overly conservative ice shape predictions were deficiencies in the leading edge heat transfer and the sensitivity of the void ice density model to the particle inertia parameter. The scallop features present on the ice shapes were thought to generate interstitial flow and horse shoe vortices which enhance the leading edge heat transfer. A set of changes to improve the leading edge heat transfer and the void density model were tested. The changes improved the ice shape predictions considerably. More work needs to be done to evaluate the performance of these modifications for a wider range of geometries and icing conditions.

Icing CFD↗

2nd-Order CESE Results For C1.1: Transonic Ringleb Flow

The Conservation Element and Solution Element (CESE) method was used as implemented in the NASA research code ez4d (an unstructured Navier-Stokes solver coded in C++ with serial and parallel versions available.) The CESE method is a time-accurate formulation with flux-conservation in both space and time. The method treats the discretized derivatives of space and time identically and while the 2nd-order accurate version was used, high-order versions exist.

Numerical Analysis↗

2nd-Order CESE Results For C1.4: Vortex Transport by Uniform Flow

The Conservation Element and Solution Element (CESE) method was used as implemented in the NASA research code ez4d. The CESE method is a time accurate formulation with flux-conservation in both space and time. The method treats the discretized derivatives of space and time identically and while the 2nd-order accurate version was used, high-order versions exist, the 2nd-order accurate version was used. In regards to the ez4d code, it is an unstructured Navier-Stokes solver coded in C++ with serial and parallel versions available. As part of its architecture, ez4d has the capability to utilize multi-thread and Messaging Passage Interface (MPI) for parallel runs.

Space-time CE/SE method↗

Reconstruction of Twist Torque in Main Parachute Risers

The reconstruction of twist torque in the Main Parachute Risers of the Capsule Parachute Assembly System (CPAS) has been successfully used to validate CPAS Model Memo conservative twist torque equations. Reconstruction of basic, one degree of freedom drop tests was used to create a functional process for the evaluation of more complex, rigid body simulation. The roll, pitch, and yaw of the body, the fly-out angles of the parachutes, and the relative location of the parachutes to the body are inputs to the torque simulation. The data collected by the Inertial Measurement Unit (IMU) was used to calculate the true torque. The simulation then used photogrammetric and IMU data as inputs into the Model Memo equations. The results were then compared to the true torque results to validate the Model Memo equations. The Model Memo parameters were based off of steel risers and the parameters will need to be re-evaluated for different materials. Photogrammetric data was found to be more accurate than the inertial data in accounting for the relative rotation between payload and cluster. The Model Memo equations were generally a good match and when not matching were generally conservative.

Day, Joshua D.↗

Cassini Attitude Control Operations - Guidelines Levied on Science to Extend Reaction Wheel Life

The Cassini spacecraft was launched on October 15, 1997 and arrived at Saturn on June 30, 2004. It has performed detailed observations and remote sensing of Saturn, its rings, and its satellites since that time. Cassini deployed the European-built Huygens probe, which descended through the Titan atmosphere (Saturn's largest moon) and landed on its surface on January 14, 2005. The Cassini mission has recently been approved by NASA to continue through September of 2017. This 7-year extension is called the Solstice mission and it presents challenges to the spacecraft operations team and its ability to maintain the health of the spacecraft. To keep the spacecraft healthy for 7 more years, the spacecraft team must carefully manage hydrazine use (about 48% of the 132 kg launch load remains as of January 2011). A vital part of conserving hydrazine is to use the reaction wheel assembly (RWA) control system for precise pointing and slews wherever possible. In any given week, the Cassini spacecraft is commanded to use RWA control about 99% of the time, with about 1% of the time requiring reaction control system (RCS) thruster control (to perform Delta V course corrections or to bias the RWA momentum). Such extensive use of the RWA hardware throughout the mission requires that the RWAs be operated in a way that minimizes degradation in the RWA electronics, DC motor, and spin bearing for each reaction wheel. Three consumables in particular have been identified for the RWAs: (1) Total number of revolutions for each RWA. (2) Time spent at very low wheel speeds. At these low speeds, good elasto-hydrodynamic (EHD) film lubrication may be compromised. (3) Total number of on/off power cycles. The second of these consumables, minimizing the time spent at very low wheel speeds, is especially important to keep the spin bearing healthy and well-lubricated. These consumables are actively managed by the attitude control operations team throughout the mission. One vital management technique is to predict individual RWA momentum (given the pointing and slews that are needed to collect the best science) and to bias the RWA momentum in a way that reduces both the total number of revolutions as well as the time spent below EHD wheel speed. Another strategy to protect RWA health is to alter the planned pointing of the spacecraft (which can affect science collection) so that the RWA consumables are conserved. This paper focuses on why this second technique is needed, and discusses how guidelines have been developed by the attitude control team which affects the planned science pointing, so that science data can be most optimally collected while still minimizing RWA consumable usage.

interplanetary trajectory↗

Comparison of TID Effects in Space-Like Variable Dose Rates and Constant Dose Rates

The degradation of the LM193 dual voltage comparator has been studied at different TID dose rate profiles, including several different constant dose rates and a variable dose rate that simulates the behavior of a solar flare. A comparison of results following constant dose rate vs. variable dose rates is made to explore how well the constant dose rates used for typical part testing predict the performance during a simulated space-like mission. Testing at a constant dose rate equal to the lowest dose rate seen during the simulated flare provides an extremely conservative estimate of the overall amount of degradation. A constant dose rate equal to the average dose rate is also more conservative than the variable rate. It appears that, for this part, weighting the dose rates by the amount of total dose received at each rate (rather than the amount of time at each dose rate) results in an average rate that produces an amount of degradation that is a reasonable approximation to that received by the variable rate.

Solar Flare↗

Cumulative Interference to Aircraft Radios from Multiple Portable Electronic Devices

Cumulative interference effects from portable electronic devices (PEDs) located inside a passenger cabin are conservatively estimated for aircraft radio receivers. PEDs' emission powers in an aircraft radio frequency band are first scaled according to their locations' interference path loss (IPL) values, and the results are summed to determine the total interference power. The multiple-equipment-factor (MEF) is determined by normalizing the result against the worst case contribution from a single device. Conservative assumptions were made and MEF calculations were performed for Boeing 737's Localizer, Glide-slope, Traffic Collision Avoidance System, and Very High Frequency Communication radio systems where full-aircraft IPL data were available. The results show MEF for the systems to vary between 10 and 14 dB. The same process was also used on the more popular window/door IPL data, and the comparison show the multiple-equipment-factor results came within one decibel (dB) of each other.

Nguyen, Truong X.↗

Robust and Simple Non-Reflecting Boundary Conditions for the Euler Equations - A New Approach based on the Space-Time CE/SE Method

This paper reports on a significant advance in the area of nonreflecting boundary conditions (NRBCs) for unsteady flow computations. As a part of t he development of t he space-time conservation element and solution element (CE/SE) method, sets of NRBCs for 1D Euler problems are developed without using any characteristics- based techniques. These conditions are much simpler than those commonly reported in the literature, yet so robust that they are applicable to subsonic, transonic and supersonic flows even in the presence of discontinuities. In addition, the straightforward multidimensional extensions of the present 1D NRBCs have been shown numerically to be equally simple and robust. The paper details the theoretical underpinning of these NRBCs, and explains t heir unique robustness and accuracy in terms of t he conservation of space-time fluxes. Some numerical results for an extended Sod's shock-tube problem, illustrating the effectiveness of the present NRBCs are included, together with an associated simple Fortran computer program. As a preliminary to the present development, a review of the basic CE/SE schemes is also included.

Chang, S.-C.↗

Enhanced Engine Performance During Emergency Operation Using a Model-Based Engine Control Architecture

This paper discusses the design and application of model-based engine control (MBEC) for use during emergency operation of the aircraft. The MBEC methodology is applied to the Commercial Modular Aero-Propulsion System Simulation 40,000 (CMAPSS40,000) and features an optimal tuner Kalman Filter (OTKF) to estimate unmeasured engine parameters, which can then be used for control. During an emergency scenario, normally-conservative engine operating limits may be relaxed to increase the performance of the engine and overall survivability of the aircraft; this comes at the cost of additional risk of an engine failure. The MBEC architecture offers the advantage of estimating key engine parameters that are not directly measureable. Estimating the unknown parameters allows for tighter control over these parameters, and on the level of risk the engine will operate at. This will allow the engine to achieve better performance than possible when operating to more conservative limits on a related, measurable parameter.

control theory↗

Simplified Methodology to Estimate the Maximum Liquid Helium (LHe) Cryostat Pressure from a Vacuum Jacket Failure

The aircraft-based Stratospheric Observatory for Infrared Astronomy (SOFIA) is a platform for multiple infrared astronomical observation experiments. These experiments carry sensors cooled to liquid helium temperatures. The liquid helium supply is contained in large (i.e., 10 liters or more) vacuum-insulated dewars. Should the dewar vacuum insulation fail, the inrushing air will condense and freeze on the dewar wall, resulting in a large heat flux on the dewar's contents. The heat flux results in a rise in pressure and the actuation of the dewar pressure relief system. A previous NASA Engineering and Safety Center (NESC) assessment provided recommendations for the wall heat flux that would be expected from a loss of vacuum and detailed an appropriate method to use in calculating the maximum pressure that would occur in a loss of vacuum event. This method involved building a detailed supercritical helium compressible flow thermal/fluid model of the vent stack and exercising the model over the appropriate range of parameters. The experimenters designing science instruments for SOFIA are not experts in compressible supercritical flows and do not generally have access to the thermal/fluid modeling packages that are required to build detailed models of the vent stacks. Therefore, the SOFIA Program engaged the NESC to develop a simplified methodology to estimate the maximum pressure in a liquid helium dewar after the loss of vacuum insulation. The method would allow the university-based science instrument development teams to conservatively determine the cryostat's vent neck sizing during preliminary design of new SOFIA Science Instruments. This report details the development of the simplified method, the method itself, and the limits of its applicability. The simplified methodology provides an estimate of the dewar pressure after a loss of vacuum insulation that can be used for the initial design of the liquid helium dewar vent stacks. However, since it is not an exact tool, final verification of the dewar pressure vessel design requires a complete, detailed real fluid compressible flow model of the vent stack. The wall heat flux resulting from a loss of vacuum insulation increases the dewar pressure, which actuates the pressure relief mechanism and results in high-speed flow through the dewar vent stack. At high pressures, the flow can be choked at the vent stack inlet, at the exit, or at an intermediate transition or restriction. During previous SOFIA analyses, it was observed that there was generally a readily identifiable section of the vent stack that would limit the flow – e.g., a small diameter entrance or an orifice. It was also found that when the supercritical helium was approximated as an ideal gas at the dewar condition, the calculated mass flow rate based on choking at the limiting entrance or transition was less than the mass flow rate calculated using the detailed real fluid model2. Using this lower mass flow rate would yield a conservative prediction of the dewar’s wall heat flux capability. The simplified method of the current work was developed by building on this observation.

Ungar, Eugene K.↗

Improving the Timeliness of Winter Wheat Production Forecast in the United States of America, Ukraine and China Using MODIS Data and NCAR Growing Degree Day Information

Wheat is themost important cereal crop traded on international markets and winter wheat constitutes approximately 80% of global wheat production. Thus, accurate and timely production forecasts are critical for making informed agricultural policies and investments, as well as increasing market efficiency and stability. Becker-Reshef et al. (2010) developed an empirical generalized model for forecasting winter wheat production. Their approach combined BRDF-corrected daily surface reflectance from Moderate resolution Imaging Spectroradiometer (MODIS) Climate Modeling Grid (CMG) with detailed official crop statistics and crop typemasks. It is based on the relationship between the Normalized Difference Vegetation Index (NDVI) at the peak of the growing season, percent wheat within the CMG pixel (area within the CMG pixel occupied by wheat), and the final yields. This method predicts the yield approximately one month to six weeks prior to harvest. In this study, we include Growing Degree Day (GDD) information extracted from NCEP/NCAR reanalysis data in order to improve the winter wheat production forecast by increasing the timeliness of the forecasts while conserving the accuracy of the original model. We apply this modified model to three major wheat-producing countries: the Unites States (US), Ukraine and China from 2001 to 2012. We show that a reliable forecast can be made between one month to a month and a half prior to the peak NDVI (meaning two months to two and a half months prior to harvest), while conserving an accuracy of 10% in the production forecast.

Production↗

Enhanced Engine Performance During Emergency Operation Using a Model-Based Engine Control Architecture

This paper discusses the design and application of model-based engine control (MBEC) for use during emergency operation of the aircraft. The MBEC methodology is applied to the Commercial Modular Aero-Propulsion System Simulation 40k (CMAPSS40k) and features an optimal tuner Kalman Filter (OTKF) to estimate unmeasured engine parameters, which can then be used for control. During an emergency scenario, normally-conservative engine operating limits may be relaxed to increase the performance of the engine and overall survivability of the aircraft; this comes at the cost of additional risk of an engine failure. The MBEC architecture offers the advantage of estimating key engine parameters that are not directly measureable. Estimating the unknown parameters allows for tighter control over these parameters, and on the level of risk the engine will operate at. This will allow the engine to achieve better performance than possible when operating to more conservative limits on a related, measurable parameter.

engine control↗

Perfluorinated Alkyl Acids in Plasma of American Alligators (Alligator Mississippiensis) from Florida and South Carolina

This study aimed to quantitate fourteen perfluoroalkyl acids (PFAAs) in 125 adult American alligators at twelve sites across the southeastern US. Of those fourteen PFAAs, nine were detected in 65% - 100% of the samples: PFOA, PFNA, PFDA, PFUnA, PFDoA, PFTriA, PFTA, PFHxS, and PFOS. Males (across all sites) showed significantly higher concentrations of four PFAAs: PFOS (p = 0.01), PFDA (p = 0.0003), PFUnA (p = 0.021), and PFTriA (p = 0.021). Concentrations of PFOS, PFHxS, and PFDA in plasma were significantly different among the sites in each sex. Alligators at Merritt Island National Wildlife Refuge and Kiawah Nature Conservancy both exhibited some of the highest PFOS concentrations (medians 99.5 ng/g and 55.8 ng/g respectively) in plasma measured to date in a crocodilian species. A number of positive correlations between PFAAs and snout-vent length (SVL) were observed in both sexes suggesting PFAA body burdens increase with increasing size. In addition, several significant correlations among PFAAs in alligator plasma may suggest conserved sources of PFAAs at each site throughout the greater study area. This study is the first to report PFAAs in American alligators, reveals potential PFAA hot spots in Florida and South Carolina, and provides and additional contaminant of concern when assessing anthropogenic impacts on ecosystem health.

crocodilians↗

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

Comparison of Fully-Compressible Equation Sets for Atmospheric Dynamics

Traditionally, the equation for the conservation of energy used in atmospheric models is based on potential temperature and is used in place of the total energy conservation. This paper compares the application of the two equations sets for both the Euler and the Navier-Stokes solutions using several benchmark test cases. A high-resolution wave-propagation method which accurately takes into account the source term due to gravity is used for computing the non-hydrostatic atmospheric flows. It is demonstrated that there is little to no difference between the results obtained using the two different equation sets for Euler as well as Navier-Stokes solutions.

Ahmad, Nashat N.↗

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗

A Ratio of Spore to Viable Organisms: A Case Study of the JPL-SAF Cleanroom

Spacecraft surfaces that are destined to land on potential life-harboring celestial bodies are required to be rigorously cleaned and continuously monitored for spore bioburden as a proxy for spacecraft cleanliness. The NASA standard assay (NSA), used for spacecraft bioburden estimates, specifically measures spores that are cultivable, aerobic, resistant to heat shock, and grow at 30 C in a nutrient-rich medium. Since the vast majority of microorganisms cannot be cultivated using the NSA, it is necessary to utilize state-of-the art molecular techniques to better understand the presence of all viable microorganisms, not just those measured with the NSA. In this study, the nutrient-deprived low biomass cleanrooms, where spacecraft are assembled, were used as a surrogate for spacecraft surfaces to measure the ratio of NSA spores in relation to the total viable microorganism population in order to make comparisons with the 2006 Space Studies Board (SSB) estimate of 1 spore per approximately 50,000 viable organisms. Ninety-eight surface wipe samples were collected from the Spacecraft Assembly Facility (SAF) cleanroom at the Jet Propulsion Laboratory (JPL) over a 6-month period. The samples were processed and analyzed using classical microbiology along with molecular methodology. Traditional microbiology plating methods were used to determine the cultivable bacterial, fungal, and spore populations. Molecular assays were used to determine the total organisms (TO, dead and live) and the viable organisms (VO, live). The TO was measured using adenosine triphosphate (ATP) and quantitative polymerase chain reaction (qPCR) assays. The VO was measured using internal ATP, propidium monoazide (PMA)-qPCR, and flow cytometry (after staining for viable microorganisms) assays. Based on the results, it was possible to establish a ratio between spore counts and VO for each viability assay. The ATP-based spore to VO ratio ranged from 149-746, and the bacterial PMA-qPCR assay-based ratio ranged from 314-1,491 VO, per spore. The most conservative estimate came from fluorescent-assisted cell sorting (FACS), which estimated the ratio to be 12,091 VO per 1 NSA spore. Since archaeal (less than 1%) and fungal (approximately 2%) populations were negligible, the spore to VO ratios were based on bacterial population estimates. The most conservative ratio from this study can be used as a replacement for the SSB estimate on nutrient-deprived (oligotrophic) desiccated spacecraft surfaces, to estimate the VO from NSA measurements without utilizing state-of-the art molecular methods that are costly and require more biomass than is typically found on spacecraft surfaces.

Hendrickson, Ryan↗