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 865 records · Page 48

Patched-grid computations of high-speed inlet flows

Two implicit algorithms for the compressible Navier-Stokes equations are applied to the computation of high-speed inlet flows typical of those which are encountered on propulsion systems of advanced hypersonic vehicles. A patched grid formulation is used to treat accurately the blunt cowl interaction with the forebody compression field and efficiently treat variations in inlet geometry downstream of the inlet entrance. Results from the two algorithms, each based on an upwind-biased spatial differencing approach for the convective and pressure terms within a flux-difference-splitting framework, are compared with experimental results for both inlets, and reasonable agreement is found.

Thomas, J. L.↗

Relationship of E-type Apollo asteroid 3103 (1982 BB) to the enstatite achondrite meteorites and the Hungaria asteroids

Visible and near-IR observations are reported which confirm that Apollo asteroid 3103 (1982 BB) is a high-albedo E-type object with a reflectance spectrum consistent with the iron-free silicates, primarily enstatite, present in the aubrites. Visible and thermal IR light curves indicate that 3103 is an elongated object with no evidence of substantial surface albedo variations. It is presently in an orbital resonance with earth and is probably a relatively long-lived member of the earth-approaching population. This is consistent with the relatively long cosmic ray exposure ages of the aubrites compared to other stony meteorites. The fall geometries of the aubrites are consistent with their being derived from 3103. The ubiquitous brecciated nature of the aubrites and the indications of a single parent body for nearly all these meteorite suggest that 3103 is the sole near-earth parent body of this meteorite type.

Gaffey, Michael J.↗

A Generalized Hydrodynamic-Impact Theory for the Loads and Motions of Deeply Immersed Prismatic Bodies

A theory is derived for determining the loads and motions of a deeply immersed prismatic body. The method makes use of a two-dimensional water-mass variation and an aspect-ratio correction for three-dimensional flow. The equations of motion are generalized by using a mean value of the aspect-ratio correction and by assuming a variation of the two-dimensional water mass for the deeply immersed body. These equations lead to impact coefficients that depend on an approach parameter which, in turn, depends upon the initial trim and flight-path angles. Comparison of experiment with theory is shown at maximum load and maximum penetration for the flat-bottom (0 deg dead-rise angle) model with bean-loading coefficients from 36.5 to 133.7 over a wide range of initial conditions. A dead-rise angle correction is applied and maximum-load data are compared with theory for the case of a model with 300 dead-rise angle and beam-loading coefficients from 208 to 530.

Markey, Melvin F.↗

Energy-Dependent Ionization States of Shock-Accelerated Particles in the Solar Corona

We examine the range of possible energy dependence of the ionization states of ions that are shock-accelerated from the ambient plasma of the solar corona. If acceleration begins in a region of moderate density, sufficiently low in the corona, ions above about 0.1 MeV/amu approach an equilibrium charge state that depends primarily upon their speed and only weakly on the plasma temperature. We suggest that the large variations of the charge states with energy for ions such as Si and Fe observed in the 1997 November 6 event are consistent with stripping in moderately dense coronal. plasma during shock acceleration. In the large solar-particle events studied previously, acceleration occurs sufficiently high in the corona that even Fe ions up to 600 MeV/amu are not stripped of electrons.

Reames, Donald V.↗

Estimating Net Primary Productivity Using Satellite and Ancillary Data

The net primary productivity (C) or the annual rate of carbon accumulation per unit ground area by terrestrial plant communities is the difference of gross photosynthesis (A(sub g)) and respiration (R) per unit ground area. Available field observations show that R is a large and variable fraction of A(sub g), although it is generally recognized that there are considerable difficulties in determining these fluxes, and thus pose challenge in assessing the accuracy. Further uncertainties arise in extrapolating field measurements (which are acquired over a hectare or so area) to regional scale. Here, an approach is presented for determining these fluxes using satellite and ancillary data to be representative of regional scale and allow assessment of interannual variation. A, has been expressed as the product of radiation use efficiency for gross photosynthesis by an unstressed canopy and intercepted photosynthetically active radiation, which is then adjusted for stresses due to soil water shortage and temperature away from optimum. R has been calculated as the sum of growth and maintenance components (respectively, R(sub g) and R(sub m)).The R(sub m) has been determined from nitrogen content of plant tissue per unit ground area, while R(sub g) has been obtained as a fraction of the difference of A(sub g) and R(sub m). Results for five consecutive years (1986-1990) are presented for the Amazon-Tocontins, Mississippi, and Ob River basins.

Choudhury, Bhaskar J.↗

GLM Proxy Data Generation: Methods for Stroke/Pulse Level Inter-Comparison of Ground-Based Lightning Reference Networks

In order to produce useful proxy data for the GOES-R Geostationary Lightning Mapper (GLM) in regions not covered by VLF lightning mapping systems, we intend to employ data produced by ground-based (regional or global) VLF/LF lightning detection networks. Before using these data in GLM Risk Reduction tasks, it is necessary to have a quantitative understanding of the performance of these networks, in terms of CG flash/stroke DE, cloud flash/pulse DE, location accuracy, and CLD/CG classification error. This information is being obtained through inter-comparison with LMAs and well-quantified VLF/LF lightning networks. One of our approaches is to compare "bulk" counting statistics on the spatial scale of convective cells, in order to both quantify relative performance and observe variations in cell-based temporal trends provided by each network. In addition, we are using microsecond-level stroke/pulse time correlation to facilitate detailed inter-comparisons at a more-fundamental level. The current development status of our ground-based inter-comparison and evaluation tools will be presented, and performance metrics will be discussed through a comparison of Vaisala s Global Lightning Dataset (GLD360) with the NLDN at locations within and outside the U.S.

Cummins, Kenneth L.↗

Combined Structural and Trajectory Control of Variable-Geometry Planetary Entry Systems

Some of the key challenges of planetary entry are to dissipate the large kinetic energy of the entry vehicle and to land with precision. Past missions to Mars were based on unguided entry, where entry vehicles carried payloads of less than 0.6 T and landed within 100 km of the designated target. The Mars Science Laboratory (MSL) is expected to carry a mass of almost 1 T to within 20 km of the target site. Guided lifting entry is needed to meet these higher deceleration and targeting demands. If the aerodynamic characteristics of the decelerator are variable during flight, more trajectory options are possible, and can be tailored to specific mission requirements. In addition to the entry trajectory modulation, having variable aerodynamic properties will also favor maneuvering of the vehicle prior to descent. For proper supersonic parachute deployment, the vehicle needs to turn to a lower angle of attack. One approach to entry trajectory improvement and angle of attack control is to embed a variable geometry decelerator in the design of the vehicle. Variation in geometry enables the vehicle to adjust its aerodynamic performance continuously without additional fuel cost because only electric power is needed for actuating the mechanisms that control the shape change. Novel structural and control concepts have been developed that enable the decelerator to undergo variation in geometry. Changing the aerodynamic characteristics of a flight vehicle by active means can potentially provide a mechanically simple, affordable, and enabling solution for entry, descent, and landing across a wide range of mission types, sample capture and return, and reentry to Earth, Titan, Venus, or Mars. Unguided ballistic entry is not sufficient to meet this more stringent deceleration, heating, and targeting demands. Two structural concepts for implementing the cone angle variation, a segmented shell, and a corrugated shell, have been presented.

Quadrelli, Marco B.↗

Supersonic Aerodynamic Characteristics of Blunt Body Trim Tab Configurations

Trim tabs are aerodynamic control surfaces that can allow an entry vehicle to meet aerodynamic performance requirements while reducing or eliminating the use of ballast mass and providing a capability to modulate the lift-to-drag ratio during entry. Force and moment data were obtained on 38 unique, blunt body trim tab configurations in the NASA Langley Research Center Unitary Plan Wind Tunnel. The data were used to parametrically assess the supersonic aerodynamic performance of trim tabs and to understand the influence of tab area, cant angle, and aspect ratio. Across the range of conditions tested (Mach numbers of 2.5, 3.5, and 4.5; angles of attack from -4deg to +20deg; angles of sideslip from 0deg to +8deg), the effects of varying tab area and tab cant angle were found to be much more significant than effects from varying tab aspect ratio. Aerodynamic characteristics exhibited variation with Mach number and forebody geometry over the range of conditions tested. Overall, the results demonstrate that trim tabs are a viable approach to satisfy aerodynamic performance requirements of blunt body entry vehicles with minimal ballast mass. For a 70deg sphere-cone, a tab with 3% area of the forebody and canted approximately 35deg with no ballast mass was found to give the same trim aerodynamics as a baseline model with ballast mass that was 5% of the total entry mass.

Korzun, Ashley M.↗

New Mathematical Functions for Vacuum System Analysis

A new bivariate function has been found that provides solutions of integrals having the form u (sup minus eta) e (sup u) du which arise when developing predictions for the behavior of pressure within a rigid volume under high vacuum conditions in the presence of venting as well as sources characterized by power law transient decay over the range [0,1] for eta and for u greater than or equal to 0. A few properties of the new function are explored in this work. For instance the eta equals 1/2 case reproduces the Dawson function. In addition, a slight variation of the solution technique reproduces the exponential integral for eta equals 1. The technique used to generate these functions leads to an approach for solving a more general class of nonlinear ordinary differential equations, with the potential for identifying other new functions that solve other integrals.

Dawson Function↗

Assimilation of Aerosol Observations in the NASA GEOS Model

In the GEOS near real-time system, as well as in MERRA-2 which is the latest reanalysis produced at NASA's Global Modeling Assimilation Office (GMAO), the assimilation of aerosol observations is performed by means of a so-called analysis splitting method. The prognostic model is based on the GEOS model radiatively coupled to GOCART aerosol module and includes assimilation of bias-corrected Aerosol Optical Depth (AOD) at 550 nm from various space-based remote sensing platforms. In line with the transition of the GEOS meteorological data assimilation system to a hybrid Ensemble-Variational formulation, we are updating the aerosol component of our assimilation system to a variational ensemble type of scheme. In this talk we will examine the impact of replacing the current analysis splitting scheme with this new approach. Starting with the assimilation of satellite based single-channel retrievals; we will discuss the impact of this aerosol data assimilation technique on the 3D aerosol distributions by means of innovation statistics and verification against independent datasets such as the Aerosol Robotic Network (AERONET) and surface PM2.5. We will also present preliminary results related to the introduction of new aerosol data types in GEOS, including multi-spectral AOD retrievals.

AERONET↗

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗

Application of a Constant Gain Extended Kalman Filter for In-Flight Estimation of Aircraft Engine Performance Parameters

An approach based on the Constant Gain Extended Kalman Filter (CGEKF) technique is investigated for the in-flight estimation of non-measurable performance parameters of aircraft engines. Performance parameters, such as thrust and stall margins, provide crucial information for operating an aircraft engine in a safe and efficient manner, but they cannot be directly measured during flight. A technique to accurately estimate these parameters is, therefore, essential for further enhancement of engine operation. In this paper, a CGEKF is developed by combining an on-board engine model and a single Kalman gain matrix. In order to make the on-board engine model adaptive to the real engine s performance variations due to degradation or anomalies, the CGEKF is designed with the ability to adjust its performance through the adjustment of artificial parameters called tuning parameters. With this design approach, the CGEKF can maintain accurate estimation performance when it is applied to aircraft engines at offnominal conditions. The performance of the CGEKF is evaluated in a simulation environment using numerous component degradation and fault scenarios at multiple operating conditions.

Kobayashi, Takahisa↗

Data Assimilation and Uncertainties in Early Solar Cycle Predictions

Stochastic nature of solar activity variations together with our limited knowledge of the dynamo mechanism and subsurface dynamics causes uncertainty in predictions of the solar cycle. For improving the physics-­‐based predictions we can take advantage of the mathematical data assimilation approach that allows us to take into account both, observational errors and model uncertainties, and provide estimates of the next solar cycle along with prediction uncertainties. In this study we use the Parker's migratory dynamo model together with the equation of magnetic helicity balance, which reproduces main properties of the sunspot cycles and allow us to minimize discrepancies between the observed global activity variations and the model solution. The test simulation runs show that a reliable prediction can be obtained for two phases of preceding solar cycle: 1) if the polar field reversals shortly after the solar maxima (strong toroidal field and weak poloidal field), and 2) during the solar minima (strongest poloidal and weak toroidal fields). The early estimate of Cycle 25 obtained by this method shows that this cycle will start in 2019 – 2020, reach the maximum in 2023 – 2024, and that the mean sunspot number at the maximum will be about 90 (for the v2.0 sunspot number series).

Data Assimilation↗

Optimal Tuner Selection for Kalman-Filter-Based Aircraft Engine Performance Estimation

An emerging approach in the field of aircraft engine controls and system health management is the inclusion of real-time, onboard models for the inflight estimation of engine performance variations. This technology, typically based on Kalman-filter concepts, enables the estimation of unmeasured engine performance parameters that can be directly utilized by controls, prognostics, and health-management applications. A challenge that complicates this practice is the fact that an aircraft engine s performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters such as efficiencies and flow capacities related to each major engine module. Through Kalman-filter-based estimation techniques, the level of engine performance degradation can be estimated, given that there are at least as many sensors as health parameters to be estimated. However, in an aircraft engine, the number of sensors available is typically less than the number of health parameters, presenting an under-determined estimation problem. A common approach to address this shortcoming is to estimate a subset of the health parameters, referred to as model tuning parameters. The problem/objective is to optimally select the model tuning parameters to minimize Kalman-filterbased estimation error. A tuner selection technique has been developed that specifically addresses the under-determined estimation problem, where there are more unknown parameters than available sensor measurements. A systematic approach is applied to produce a model tuning parameter vector of appropriate dimension to enable estimation by a Kalman filter, while minimizing the estimation error in the parameters of interest. Tuning parameter selection is performed using a multi-variable iterative search routine that seeks to minimize the theoretical mean-squared estimation error of the Kalman filter. This approach can significantly reduce the error in onboard aircraft engine parameter estimation applications such as model-based diagnostic, controls, and life usage calculations. The advantage of the innovation is the significant reduction in estimation errors that it can provide relative to the conventional approach of selecting a subset of health parameters to serve as the model tuning parameter vector. Because this technique needs only to be performed during the system design process, it places no additional computation burden on the onboard Kalman filter implementation. The technique has been developed for aircraft engine onboard estimation applications, as this application typically presents an under-determined estimation problem. However, this generic technique could be applied to other industries using gas turbine engine technology.

Simon, Donald L.↗

Generating Vegetation Leaf Area Index Earth System Data Record from Multiple Sensors: Theory - Part 1

The generation of multi-decade long Earth System Data Records (ESDRs) of Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation absorbed by vegetation (FPAR) from remote sensing measurements of multiple sensors is key to monitoring long-term changes in vegetation due to natural and anthropogenic influences. Challenges in developing such ESDRs include problems in remote sensing science (modeling of variability in global vegetation, scaling, atmospheric correction) and sensor hardware (differences in spatial resolution, spectral bands, calibration, and information content). In this paper, we develop a physically based approach for deriving LAI and FPAR products from the Advanced Very High Resolution Radiometer (AVHRR) data that are of comparable quality to the Moderate resolution Imaging Spectroradiometer (MODIS) LAI and FPAR products, thus realizing the objective of producing a long (multi-decadal) time series of these products. The approach is based on the radiative transfer theory of canopy spectral invariants which facilitates parameterization of the canopy spectral bidirectional reflectance factor (BRF). The methodology permits decoupling of the structural and radiometric components and obeys the energy conservation law. The approach is applicable to any optical sensor, however, it requires selection of sensor-specific values of configurable parameters, namely, the single scattering albedo and data uncertainty. According to the theory of spectral invariants, the single scattering albedo is a function of the spatial scale, and thus, accounts for the variation in BRF with sensor spatial resolution. Likewise, the single scattering albedo accounts for the variation in spectral BRF with sensor bandwidths. The second adjustable parameter is data uncertainty, which accounts for varying information content of the remote sensing measurements, i.e., Normalized Difference Vegetation Index (NDVI, low information content), vs. spectral BRF (higher information content). Implementation of this approach indicates good consistency in LAI values retrieved from NDVI (AVHRRmode) and spectral BRF (MODIS-mode). Specific details of the implementation and evaluation of the derived products are detailed in the second part of this two-paper series.

Ganguly, Sangram↗

An Improved Cryosat-2 Sea Ice Freeboard Retrieval Algorithm Through the Use of Waveform Fitting

We develop an empirical model capable of simulating the mean echo power cross product of CryoSat-2 SAR and SAR In mode waveforms over sea ice covered regions. The model simulations are used to show the importance of variations in the radar backscatter coefficient with incidence angle and surface roughness for the retrieval of surfaceelevation of both sea ice floes and leads. The numerical model is used to fit CryoSat-2 waveforms to enable retrieval of surface elevation through the use of look-up tables and a bounded trust region Newton least squares fitting approach. The use of a model to fit returns from sea ice regions offers advantages over currently used threshold retrackingmethods which are here shown to be sensitive to the combined effect of bandwidth limited range resolution and surface roughness variations. Laxon et al. (2013) have compared ice thickness results from CryoSat-2 and IceBridge, and found good agreement, however consistent assumptions about the snow depth and density of sea ice werenot used in the comparisons. To address this issue, we directly compare ice freeboard and thickness retrievals from the waveform fitting and threshold tracker methods of CryoSat-2 to Operation IceBridge data using a consistent set of parameterizations. For three IceBridge campaign periods from March 20112013, mean differences (CryoSat-2 IceBridge) of 0.144m and 1.351m are respectively found between the freeboard and thickness retrievals using a 50 sea ice floe threshold retracker, while mean differences of 0.019m and 0.182m are found when using the waveform fitting method. This suggests the waveform fitting technique is capable of better reconciling the seaice thickness data record from laser and radar altimetry data sets through the usage of consistent physical assumptions.

CryoSat-2↗

The Threat of Uncertainty: Why Using Traditional Approaches for Evaluating Spacecraft Reliability are Insufficient for Future Human Mars Missions

Through the Evolvable Mars Campaign (EMC) study, the National Aeronautics and Space Administration (NASA) continues to evaluate potential approaches for sending humans beyond low Earth orbit (LEO). A key aspect of these missions is the strategy that is employed to maintain and repair the spacecraft systems, ensuring that they continue to function and support the crew. Long duration missions beyond LEO present unique and severe maintainability challenges due to a variety of factors, including: limited to no opportunities for resupply, the distance from Earth, mass and volume constraints of spacecraft, high sensitivity of transportation element designs to variation in mass, the lack of abort opportunities to Earth, limited hardware heritage information, and the operation of human-rated systems in a radiation environment with little to no experience. The current approach to maintainability, as implemented on ISS, which includes a large number of spares pre-positioned on ISS, a larger supply sitting on Earth waiting to be flown to ISS, and an on demand delivery of logistics from Earth, is not feasible for future deep space human missions. For missions beyond LEO, significant modifications to the maintainability approach will be required.Through the EMC evaluations, several key findings related to the reliability and safety of the Mars spacecraft have been made. The nature of random and induced failures presents significant issues for deep space missions. Because spare parts cannot be flown as needed for Mars missions, all required spares must be flown with the mission or pre-positioned. These spares must cover all anticipated failure modes and provide a level of overall reliability and safety that is satisfactory for human missions. This will require a large amount of mass and volume be dedicated to storage and transport of spares for the mission. Further, there is, and will continue to be, a significant amount of uncertainty regarding failure rates for spacecraft components. This uncertainty makes it much more difficult to anticipate failures and will potentially require an even larger amount of spares to provide an acceptable level of safety. Ultimately, the approach to maintenance and repair applied to ISS, focusing on the supply of spare parts, may not be tenable for deep space missions. Other approaches, such as commonality of components, simplification of systems, and in-situ manufacturing will be required.

Stromgren, Chel↗

On the accuracy of shape sensitivity

The calculation of sensitivity of the response of a structure modeled by finite elements to shape variation is known to be subject to numerical difficulties. The accuracy of a given method is typically measured against the yard stick of finite-difference sensitivity calculation. The present paper demonstrates with a simple example that this approach may be flawed because of discretization errors associated with the finite element mesh. Seven methods for calculating sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the mesh is refined, displacement sensitivity derivatives converge more slowly than the displacements. Six of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer rather than reduced accuracy. Additionally, it is observed that small derivatives are particularly prone to accuracy problems.

Haftka, R. T.↗