Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global Optimization”

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,297 records · Page 72

Execution time support for scientific programs on distributed memory machines

Optimizations are considered that are required for efficient execution of code segments that consists of loops over distributed data structures. The PARTI (Parallel Automated Runtime Toolkit at ICASE) execution time primitives are designed to carry out these optimizations and can be used to implement a wide range of scientific algorithms on distributed memory machines. These primitives allow the user to control array mappings in a way that gives an appearance of shared memory. Computations can be based on a global index set. Primitives are used to carry out gather and scatter operations on distributed arrays. Communications patterns are derived at runtime, and the appropriate send and receive messages are automatically generated.

Berryman, Harry↗

A Self-Tuning Kalman Filter for Autonomous Spacecraft Navigation

Most navigation systems currently operated by NASA are ground-based, and require extensive support to produce accurate results. Recently developed systems that use Kalman Filter and Global Positioning System (GPS) data for orbit determination greatly reduce dependency on ground support, and have potential to provide significant economies for NASA spacecraft navigation. Current techniques of Kalman filtering, however, still rely on manual tuning from analysts, and cannot help in optimizing autonomy without compromising accuracy and performance. This paper presents an approach to produce a high accuracy autonomous navigation system fully integrated with the flight system. The resulting system performs real-time state estimation by using an Extended Kalman Filter (EKF) implemented with high-fidelity state dynamics model, as does the GPS Enhanced Orbit Determination Experiment (GEODE) system developed by the NASA Goddard Space Flight Center. Augmented to the EKF is a sophisticated neural-fuzzy system, which combines the explicit knowledge representation of fuzzy logic with the learning power of neural networks. The fuzzy-neural system performs most of the self-tuning capability and helps the navigation system recover from estimation errors. The core requirement is a method of state estimation that handles uncertainties robustly, capable of identifying estimation problems, flexible enough to make decisions and adjustments to recover from these problems, and compact enough to run on flight hardware. The resulting system can be extended to support geosynchronous spacecraft and high-eccentricity orbits. Mathematical methodology, systems and operations concepts, and implementation of a system prototype are presented in this paper. Results from the use of the prototype to evaluate optimal control algorithms implemented are discussed. Test data and major control issues (e.g., how to define specific roles for fuzzy logic to support the self-learning capability) are also discussed. In addition, architecture of a complete end-to-end candidate flight system that provides navigation with highly autonomous control using data from GPS is presented.

Truong, Son H.↗

Algorithms for Maneuvering Spacecraft Around Small Bodies

A document describes mathematical derivations and applications of autonomous guidance algorithms for maneuvering spacecraft in the vicinities of small astronomical bodies like comets or asteroids. These algorithms compute fuel- or energy-optimal trajectories for typical maneuvers by solving the associated optimal-control problems with relevant control and state constraints. In the derivations, these problems are converted from their original continuous (infinite-dimensional) forms to finite-dimensional forms through (1) discretization of the time axis and (2) spectral discretization of control inputs via a finite number of Chebyshev basis functions. In these doubly discretized problems, the Chebyshev coefficients are the variables. These problems are, variously, either convex programming problems or programming problems that can be convexified. The resulting discrete problems are convex parameter-optimization problems; this is desirable because one can take advantage of very efficient and robust algorithms that have been developed previously and are well established for solving such problems. These algorithms are fast, do not require initial guesses, and always converge to global optima. Following the derivations, the algorithms are demonstrated by applying them to numerical examples of flyby, descent-to-hover, and ascent-from-hover maneuvers.

Acikmese, A. Bechet↗

Dryness controls temperature-optimized gross primary productivity across vegetation types

Temperature response of gross primary productivity (GPP) is a well-known property of ecosystem, but GPP at the optimum temperature (GPP_T opt ) has not been fully discussed. Our understanding of how GPP_T opt responds to warming and water availability is highly limited. Here, in this study, we analyzed data at 326 globally distributed eddy covariance sites (79°N-37°S), to identify controlling factors of GPP_T opt . Although GPP_T opt was significantly influenced by soil moisture, global solar radiation, mean annual temperature, and vapor pressure deficit in a non-linear pattern (R 2 = 0.47), the direction and magnitude of these climate variables’ effects on GPP_T opt depend on the dryness index (DI), a ratio of potential evapotranspiration to precipitation. The spatial pattern showed that soil moisture did not affect GPP_T opt across energy-limited sites with DI < 1 while dominated GPP_T opt across water-limited sites with DI >1. The temporal pattern showed that GPP_T opt was lowered by warming or low precipitation in water-limited sites while energy-limited sites tended to maintain a stable GPP_T opt regardless of changes in air temperature. Vegetation types in humid climates tended to have higher GPP_T opt and were more likely to benefit from a warmer climate since it was not restricted by water conditions. This study highlights that the response of GPP_T opt to global warming depends on the dryness conditions, which explains the nonlinear control of water and temperature over GPP_T opt . Our finding is essential to realistic prediction of terrestrial carbon uptake under future climate and vegetation conditions.

54 ENVIRONMENTAL SCIENCES↗

Creating Near Real-time Hourly VIIRS Radiance Mosaic over Polar Regions from SatCORPS Global Cloud Composite (GCC)

While polar orbiting satellites are essential in providing the needed observations over polar regions for data assimilation, weather diagnoses/forecasting, and flight guidance and situation awareness in airborne field campaigns, there is limited temporal coverage. Often, polar mosaics are created by overlaying multiple satellite imagers on top of each other which can lead to transition lines and discontinuities in the data products. To help ease some of these shortcomings, we utilize the NASA Langley Satellite ClOud and Radiation Property retrieval System (SatCORPS) Mosaic algorithm for optimally fusing and merging selected imager radiances from the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard the Suomi National polar-orbiting Partnership (SNPP), Joint Polar Satellite System (JPSS-1), and JPSS-2 satellites, and the MODerate resolution Imaging Spectroradiometer (MODIS) instruments onboard AQUA and TERRA satellites. Since multiple satellites overpasses and instruments are used in creating hourly 3-km composites over the polar regions, selection of the best satellite data for each 3-km composite pixel is based on aggregated ranking of satellite resolution, time, viewing zenith angle, and distance from terminator. This ranking system provides a smoother transition and improves data discontinuity in the merged output resulting in a more seamless polar region dataset. We will discuss the use of this near real-time dataset to support polar airborne field campaigns and SatCORPS satellite-derived Global Cloud Composite (GCC) generation.

Satellite composites↗

Transcriptome and metabolome integration in sugarcane through culm development

Abstract Sugarcane (Saccharum sp.) is a tropical and subtropical C4 plant with a high photosynthetic and carbon assimilation efficiency that stores sucrose. Culm biomass is also composed of bagasse fibre, a by‐product of the sugarcane industry. This high‐yielding grass, high in sucrose and lignocellulosic biomass, is considered an optimal feedstock as an alternative to fossil fuels and to produce a broad range of high‐value biomaterials. The ideal sugarcane production system would optimise the relative production of sugar and these new products. Multi‐omics correlation analysis was used to generate a global view of the essential metabolic pathways identifying critical genes involved in carbon partitioning during different stages of development. This research employed an unprecedented metabolic and transcriptomic dataset of 360 samples from a selection of 1440 culms of 24 genotypes at five different development stages. Chemical composition and metabolome analysis showed an increase through the culm development of lignin, sucrose, carbon, and amino acids such as aspartic acid, serine, alanine, methionine, threonine 3‐cyano‐L‐alanine, and citric acid. Transcriptome analysis revealed functionalities such as transcription, nucleotide transport and metabolism, and the biosynthesis of amino acids that are highly activated during the immature stage and highly down‐regulated during the most mature age.

Perlo, Virginie↗

A Type II Hamiltonian Variational Principle and Adjoint Systems for Lie Groups

We present a novel Type II variational principle on the cotangent bundle of a Lie group which enforces Type II boundary conditions, i.e., fixed initial position and final momentum. In general, such Type II variational principles are only globally defined on vector spaces or locally defined on general manifolds; however, by left translation, we are able to define this variational principle globally on cotangent bundles of Lie groups. Type II boundary conditions are particularly important for adjoint sensitivity analysis, which is our motivating application. As such, we additionally discuss adjoint systems on Lie groups, their properties, and how they can be used to solve optimization problems subject to dynamics on Lie groups.

97 MATHEMATICS AND COMPUTING↗

In-situ resonant band engineering of solution-processed semiconductors generates high performance n-type thermoelectric nano-inks

Thermoelectric devices possess enormous potential to reshape the global energy landscape by converting waste heat into electricity, yet their commercial implementation has been limited by their high cost to output power ratio. No single “champion” thermoelectric material exists due to a broad range of material-dependent thermal and electrical property optimization challenges. While the advent of nanostructuring provided a general design paradigm for reducing material thermal conductivities, there exists no analogous strategy for homogeneous, precise doping of materials. Here, we demonstrate a nanoscale interface-engineering approach that harnesses the large chemically accessible surface areas of nanomaterials to yield massive, finely-controlled, and stable changes in the Seebeck coefficient, switching a poor nonconventional p- type thermoelectric material, tellurium, into a robust n -type material exhibiting stable properties over months of testing. These remodeled, n -type nanowires display extremely high power factors (~500 µW m -1 K -2 ) that are orders of magnitude higher than their bulk p -type counterparts.

42 ENGINEERING↗

An experimental characterization of core turbulence regimes in Wendelstein 7-X

First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, power balance analysis has shown that turbulent transport poses a serious limitation to the global performance of the machine. Furthermore, several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of T i ≃ 1.7 keV, except in the few scenarios in which turbulence can be suppressed. In order to understand the precise turbulent mechanisms at play and thus design improved performance scenarios, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. As a first step in this direction, in this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ( k ⊥ ρ i ≃ 1) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients and . Furthermore, we discuss the importance of other potentially relevant parameters such as T e / T i , E r or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio η i = , as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which T i,core > 1.7 keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mars gravity field based on a short-arc technique

The magnitudes of 92 surface mass points at designated locations were estimated from the radio tracking data of the Mariner Mars 1971 (M9) orbiter. This result is the first mass point model of a global field. The derived surface mass distribution correlates positively with the visible topography. The Hellas basin contains a mass deficiency, in contrast to some of the lunar basins which contain mass excesses. The Mars gravity field represented by the four parameters of an optimally located mass point (superimposed on an oblate spheroid) has third- and fourth-degree harmonics comparable to those of the complete model.

Sjogren, W. L.↗

The potential application of the blackboard model of problem solving to multidisciplinary design

Problems associated with the sequential approach to multidisciplinary design are discussed. A blackboard model is suggested as a potential tool for implementing the multilevel decomposition approach to overcome these problems. The blackboard model serves as a global database for the solution with each discipline acting as a knowledge source for updating the solution. With this approach, it is possible for engineers to improve the coordination, communication, and cooperation in the conceptual design process, allowing them to achieve a more optimal design from an interdisciplinary standpoint.

Rogers, J. L.↗

Computer aiding for low-altitude helicopter flight

A computer-aiding concept for low-altitude helicopter flight was developed and evaluated in a real-time piloted simulation. The concept included an optimal control trajectory-generated algorithm based on dynamic programming, and a head-up display (HUD) presentation of a pathway-in-the-sky, a phantom aircraft, and flight-path vector/predictor symbol. The trajectory-generation algorithm uses knowledge of the global mission requirements, a digital terrain map, aircraft performance capabilities, and advanced navigation information to determine a trajectory between mission waypoints that minimizes threat exposure by seeking valleys. The pilot evaluation was conducted at NASA Ames Research Center's Sim Lab facility in both the fixed-base Interchangeable Cab (ICAB) simulator and the moving-base Vertical Motion Simulator (VMS) by pilots representing NASA, the U.S. Army, and the U.S. Air Force. The pilots manually tracked the trajectory generated by the algorithm utilizing the HUD symbology. They were able to satisfactorily perform the tracking tasks while maintaining a high degree of awareness of the outside world.

Swenson, Harry N.↗

Computer aiding for low-altitude helicopter flight

A computer-aiding concept for low-altitude helicopter flight was developed and evaluated in a real-time piloted simulation. The concept included an optimal control trajectory-generated algorithm based on dynamic programming, and a head-up display (HUD) presentation of a pathway-in-the-sky, a phantom aircraft, and flight-path vector/predictor symbol. The trajectory-generation algorithm uses knowledge of the global mission requirements, a digital terrain map, aircraft performance capabilities, and advanced navigation information to determine a trajectory between mission waypoints that minimizes threat exposure by seeking valleys. The pilot evaluation was conducted at NASA Ames Research Center's Sim Lab facility in both the fixed-base Interchangeable Cab (ICAB) simulator and the moving-base Vertical Motion Simulator (VMS) by pilots representing NASA, the U.S. Army, and the U.S. Air Force. The pilots manually tracked the trajectory generated by the algorithm utilizing the HUD symbology. They were able to satisfactorily perform the tracking tasks while maintaining a high degree of awareness of the outside world.

Swenson, Harry N.↗

Assessment of the 1997-1998 Asian Monsoon Anomalies

Using State-of-the-art satellite-gauge monthly rainfall estimate and optimally interpolated sea surface temperature (SST) data, we have assessed the 1997-98 Asian monsoon anomalies in terms of three basic causal factors: basin-scale SST, regional coupling, and internal variability. Singular Value Decomposition analysis of rainfall and SST are carried out globally over the entire tropics and regionally over the Asian monsoon domain. Contributions to monsoon rainfall predictability by various factors are evaluated from cumulative anomaly correlation with dominant regional SVD modes. Results reveal a dominant, large-scale monsoon-El Nino coupled mode with well-defined centers of action in the near-equatorial monsoon regions. it is noted that some subcontinental regions such as all-India, or arbitrarily chosen land regions over East Asia, while important socio-economically, are not near the centers of influence from El Nino, hence are not necessarily representative of the response of the entire monsoon region to El Nino. The observed 1997-98 Asian monsoon anomalies are found to be very complex with approximately 34% of the anomalies attributable to basin- scale SST influence associated with El Nino. Regional coupled processes contribute an additional 19%, leaving about 47% due to internal dynamics. Also noted is that the highest monsoon predictability is not necessary associated with major El Nino events (e.g. 1997, 1982) but rather in non-El Nino years (e.g. 1980, 1988) when contributions from the regional coupled modes far exceed those from the basin-scale SST. The results suggest that in order to improve monsoon seasonal-to-interannual predictability, there is a need to exploit not only monsoon-El Nino relationship, but also monsoon regional coupled processes and their modulation by long-term climate change.

Lau, William K.-M.↗

Decadal Prediction Skill in the GEOS-5 Forecast System

A suite of decadal predictions has been conducted with the NASA Global Modeling and Assimilation Office?s GEOS-5 Atmosphere-Ocean General Circulation Model (AOGCM). The hindcasts are initialized every December from 1959 to 2010 following the CMIP5 experimental protocol for decadal predictions. The initial conditions are from a multi-variate ensemble optimal interpolation ocean and sea-ice reanalysis, and from the atmospheric reanalysis (MERRA, the Modern-Era Retrospective Analysis for Research and Applications) generated using the GEOS-5 atmospheric model. The forecast skill of a three-member-ensemble mean is compared to that of an experiment without initialization but forced with observed CO2. The results show that initialization acts to increase the forecast skill of Northern Atlantic SST compared to the uninitialized runs, with the increase in skill maintained for almost a decade over the subtropical and mid-latitude Atlantic. The annual-mean Atlantic Meridional Overturning Circulation (AMOC) index is predictable up to a 5-year lead time, consistent with the predictable signal in upper ocean heat content over the Northern Atlantic. While the skill measured by Mean Squared Skill Score (MSSS) shows 50% improvement up to 10-year lead forecast over the subtropical and mid-latitude Atlantic, however, prediction skill is relatively low in the subpolar gyre, due in part to the fact that the spatial pattern of the dominant simulated decadal mode in upper ocean heat content over this region appears to be unrealistic. An analysis of the large-scale temperature budget shows that this is the result of a model bias, implying that realistic simulation of the climatological fields is crucial for skillful decadal forecasts.

Ham, Yoo-Geun↗

Linear Subpixel Learning Algorithm for Land Cover Classification from WELD using High Performance Computing

In this work, we use a Fully Constrained Least Squares Subpixel Learning Algorithm to unmix global WELD (Web Enabled Landsat Data) to obtain fractions or abundances of substrate (S), vegetation (V) and dark objects (D) classes. Because of the sheer nature of data and compute needs, we leveraged the NASA Earth Exchange (NEX) high performance computing architecture to optimize and scale our algorithm for large-scale processing. Subsequently, the S-V-D abundance maps were characterized into 4 classes namely, forest, farmland, water and urban areas (with NPP-VIIRS-national polar orbiting partnership visible infrared imaging radiometer suite nighttime lights data) over California, USA using Random Forest classifier. Validation of these land cover maps with NLCD (National Land Cover Database) 2011 products and NAFD (North American Forest Dynamics) static forest cover maps showed that an overall classification accuracy of over 91 percent was achieved, which is a 6 percent improvement in unmixing based classification relative to per-pixel-based classification. As such, abundance maps continue to offer an useful alternative to high-spatial resolution data derived classification maps for forest inventory analysis, multi-class mapping for eco-climatic models and applications, fast multi-temporal trend analysis and for societal and policy-relevant applications needed at the watershed scale.

Subpixel↗

Quantum Dynamics at Scale: Ultrafast Control of Emergent Functional Materials

Confluence of extreme-scale quantum dynamics simulations (i.e. quantum@scale) and cutting-edge x-ray free-electron laser experiments are revolutionizing materials science. An archetypal example is the exciting concept of using picosecond light pulses to control emergent material properties on demand in atomically-thin layered materials. This paper describes efforts to scale our quantum molecular dynamics engine toward the United States' first exaflop/s computer, under an Aurora Early Science Program project named "Metascalable layered material genome". Key algorithmic and computing techniques incorporated are: (1) globally-scalable and locally-fast solvers within a linear-scaling divide-conquer-recombine algorithmic framework; (2) algebraic `BLASification' of computational kernels; and (3) data alignment and loop restructuring, along with register and cache blocking, for enhanced vectorization and efficient memory access. The resulting weak-scaling parallel efficiency was 0.93 on 131,072 Intel Xeon Phi cores for a 56.6 million atom (or 169 million valence-electron) system, whereas the various code transformations achieved 5-fold speedup. The optimized simulation engine allowed us for the first time to establish a significant effect of substrate on the dynamics of layered material upon electronic excitation.

Algebraic BLASification↗

Performance of an Optimized Eta Model Code on the Cray T3E and a Network of PCs

In the year 2001, NASA will launch the satellite TRIANA that will be the first Earth observing mission to provide a continuous, full disk view of the sunlit Earth. As a part of the HPCC Program at NASA GSFC, we have started a project whose objectives are to develop and implement a 3D cloud data assimilation system, by combining TRIANA measurements with model simulation, and to produce accurate statistics of global cloud coverage as an important element of the Earth's climate. For simulation of the atmosphere within this project we are using the NCEP/NOAA operational Eta model. In order to compare TRIANA and the Eta model data on approximately the same grid without significant downscaling, the Eta model will be integrated at a resolution of about 15 km. The integration domain (from -70 to +70 deg in latitude and 150 deg in longitude) will cover most of the sunlit Earth disc and will continuously rotate around the globe following TRIANA. The cloud data assimilation is supposed to run and produce 3D clouds on a near real-time basis. Such a numerical setup and integration design is very ambitious and computationally demanding. Thus, though the Eta model code has been very carefully developed and its computational efficiency has been systematically polished during the years of operational implementation at NCEP, the current MPI version may still have problems with memory and efficiency for the TRIANA simulations. Within this work, we optimize a parallel version of the Eta model code on a Cray T3E and a network of PCs (theHIVE) in order to improve its overall efficiency. Our optimization procedure consists of introducing dynamically allocated arrays to reduce the size of static memory, and optimizing on a single processor by splitting loops to limit the number of streams. All the presented results are derived using an integration domain centered at the equator, with a size of 60 x 60 deg, and with horizontal resolutions of 1/2 and 1/3 deg, respectively. In accompanying charts we report the elapsed time, the speedup and the Mflops as a function of the number of processors for the non-optimized version of the code on the T3E and theHIVE. The large amount of communication required for model integration explains its poor performance on theHIVE. Our initial implementation of the dynamic memory allocation has contributed to about 12% reduction of memory but has introduced a 3% overhead in computing time. This overhead was removed by performing loop splitting in some of the high demanding subroutines. When the Eta code is fully optimized in order to meet the memory requirement for TRIANA simulations, a non-negligeable overhead may appear that may seriously affect the efficiency of the code. To alleviate this problem, we are considering implementation of a new algorithm for the horizontal advection that is computationally less expensive, and also a new approach for marching in time.

Kouatchou, Jules↗