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 901 records · Page 50

Simultaneous Retrieval of Aerosol and Cloud Properties During the MILAGRO Field Campaign

Estimation of Direct Climate Forcing (DCF) due to aerosols in cloudy areas has historically been a difficult task, mainly because of a lack of appropriate measurements. Recently, passive remote sensing instruments have been developed that have the potential to retrieve both cloud and aerosol properties using polarimetric, multiple view angle, and multi spectral observations, and therefore determine DCF from aerosols above clouds. One such instrument is the Research Scanning Polarimeter (RSP), an airborne prototype of a sensor on the NASA Glory satellite, which unfortunately failed to reach orbit during its launch in March of 2011. In the spring of 2006, the RSP was deployed on an aircraft based in Veracruz, Mexico, as part of the Megacity Initiative: Local and Global Research Observations (MILAGRO) field campaign. On 13 March, the RSP over flew an aerosol layer lofted above a low altitude marine stratocumulus cloud close to shore in the Gulf of Mexico. We investigate the feasibility of retrieving aerosol properties over clouds using these data. Our approach is to first determine cloud droplet size distribution using the angular location of the cloud bow and other features in the polarized reflectance. The selected cloud was then used in a multiple scattering radiative transfer model optimization to determine the aerosol optical properties and fine tune the cloud size distribution. In this scene, we were able to retrieve aerosol optical depth, the fine mode aerosol size distribution parameters and the cloud droplet size distribution parameters to a degree of accuracy required for climate modeling. This required assumptions about the aerosol vertical distribution and the optical properties of the coarse aerosol size mode. A sensitivity study was also performed to place this study in the context of future systematic scanning polarimeter observations, which found that the aerosol complex refractive index can also be observed accurately if the aerosol optical depth is larger than roughly 0.8 at a wavelength of (0.555 m).

Knobelspiesse, K.↗

Cubesat Constellation Design for Air Traffic Monitoring

Suitably equipped global and local air traffic can be tracked. The tracking information may then be used for control from ground-based stations by receiving the Automatic Dependent Surveillance-Broadcast (ADS-B) signal. The ADS-B signal, emitted from the aircraft's Mode-S transponder, is currently tracked by terrestrial based receivers but not over remote oceans or sparsely populated regions such as Alaska or the Pacific Ocean. Lack of real-time aircraft time/location information in remote areas significantly hinders optimal planning and control because bigger "safety bubbles" (lateral and vertical separation) are required around the aircraft until they reach radar-controlled airspace. Moreover, it presents a search-and-rescue bottleneck. Aircraft in distress, e.g. Air France AF449 that crashed in 2009, take days to be located or cannot be located at all, e.g. Malaysia Airlines MH370 in 2014. In this paper, we describe a tool for designing a constellation of small satellites which demonstrates, through high-fidelity modeling based on simulated air traffic data, the value of space-based ADS-B monitoring and provides recommendations for cost-efficient deployment of a constellation of small satellites to increase safety and situational awareness in the currently poorly-served surveillance area of Alaska. Air traffic data has been obtained from the Future ATM Concepts Evaluation Tool (FACET), developed at NASA Ames Research Center, simulated over the Alaskan airspace over a period of one day. The simulation is driven by MATLAB with satellites propagated and coverage calculated using AGI's Satellite ToolKit(STK10).

cubesats↗

Applications of Dynamic Land Surface Information for Passive Microwave Precipitation Retrieval

Accurate, physically based precipitation retrieval over global land surfaces is an important goal of the NASA/JAXA Global Precipitation Measurement Mission (GPM). This is a difficult problem for the passive microwave constellation, as the signal over radiometrically warm land surfaces in the microwave frequencies means that the measurements used are indirect and typically require inferring some type of relationship between an observed scattering signal and precipitation at the surface. GPM, with collocated radiometer and dual-frequency radar, is an excellent tool for tackling this problem and improving global retrievals. In the years following the launch of the GPM Core Observatory satellite, physically based passive microwave retrieval of precipitation over land continues to be challenging. Validation efforts suggest that the operational GPM passive microwave algorithm, the Goddard profiling algorithm (GPROF), tends to overestimate precipitation at the low (<5 mm h−1) end of the distribution over land. In this work, retrieval sensitivities to dynamic surface conditions are explored through enhancement of the algorithm with dynamic, retrieved information from a GPM-derived optimal estimation scheme. The retrieved parameters describing surface and background characteristics replace current static or ancillary GPROF information including emissivity, water vapor, and snow cover. Results show that adding this information decreases probability of false detection by 50% and, most importantly, the enhancements with retrieved parameters move the retrieval away from dependence on ancillary datasets and lead to improved physical consistency.

Sarah Ringerud↗

Wavelet Approximation in Data Assimilation

Estimation of the state of the atmosphere with the Kalman filter remains a distant goal because of high computational cost of evolving the error covariance for both linear and nonlinear systems. Wavelet approximation is presented here as a possible solution that efficiently compresses both global and local covariance information. We demonstrate the compression characteristics on the the error correlation field from a global two-dimensional chemical constituent assimilation, and implement an adaptive wavelet approximation scheme on the assimilation of the one-dimensional Burger's equation. In the former problem, we show that 99%, of the error correlation can be represented by just 3% of the wavelet coefficients, with good representation of localized features. In the Burger's equation assimilation, the discrete linearized equations (tangent linear model) and analysis covariance are projected onto a wavelet basis and truncated to just 6%, of the coefficients. A nearly optimal forecast is achieved and we show that errors due to truncation of the dynamics are no greater than the errors due to covariance truncation.

Tangborn, Andrew↗

Robotically Assembled Aerospace Structures: Digital Material Assembly using a Gantry-Type Assembler

This paper evaluates the development of automated assembly techniques for discrete lattice structures using a multi-axis gantry type CNC machine. These lattices are made of discrete components called digital materials. We present the development of a specialized end effector that works in conjunction with the CNC machine to assemble these lattices. With this configuration we are able to place voxels at a rate of 1.5 per minute. The scalability of digital material structures due to the incremental modular assembly is one of its key traits and an important metric of interest. We investigate the build times of a 5x5 beam structure on the scale of 1 meter (325 parts), 10 meters (3,250 parts), and 30 meters (9,750 parts). Utilizing the current configuration with a single end effector, performing serial assembly with a globally fixed feed station at the edge of the build volume, the build time increases according to a scaling law of n4, where n is the build scale. Build times can be reduced significantly by integrating feed systems into the gantry itself, resulting in a scaling law of n3. A completely serial assembly process will encounter time limitations as build scale increases. Automated assembly for digital materials can assemble high performance structures from discrete parts, and techniques such as built in feed systems, parallelization, and optimization of the fastening process will yield much higher throughput.

Mechanical Structure↗

Robotically Assembled Aerospace Structures: Digital Material Assembly using a Gantry-Type Assembler

This paper evaluates the development of automated assembly techniques for discrete lattice structures using a multi-axis gantry type CNC machine. These lattices are made of discrete components called "digital materials." We present the development of a specialized end effector that works in conjunction with the CNC machine to assemble these lattices. With this configuration we are able to place voxels at a rate of 1.5 per minute. The scalability of digital material structures due to the incremental modular assembly is one of its key traits and an important metric of interest. We investigate the build times of a 5x5 beam structure on the scale of 1 meter (325 parts), 10 meters (3,250 parts), and 30 meters (9,750 parts). Utilizing the current configuration with a single end effector, performing serial assembly with a globally fixed feed station at the edge of the build volume, the build time increases according to a scaling law of n4, where n is the build scale. Build times can be reduced significantly by integrating feed systems into the gantry itself, resulting in a scaling law of n3. A completely serial assembly process will encounter time limitations as build scale increases. Automated assembly for digital materials can assemble high performance structures from discrete parts, and techniques such as built in feed systems, parallelization, and optimization of the fastening process will yield much higher throughput.

Manufacturing↗

High quality Gaussian basis sets for fourth-row atoms

Energy optimized Gaussian basis sets of triple-zeta quality for the atoms Rb-Xe have been derived. Two series of basis sets are developed: (24s 16p 10d) and (26s 16p 10d) sets which were expanded to 13d and 19p functions as the 4d and 5p shells become occupied. For the atoms lighter than Cd, the (24s 16p 10d) sets with triple-zeta valence distributions are higher in energy than the corresponding double-zeta distribution. To ensure a triple-zeta distribution and a global energy minimum, the (26s 16p 10d) sets were derived. Total atomic energies from the largest basis sets are between 198 and 284 (mu)E(sub H) above the numerical Hartree-Fock energies.

Partridge, Harry↗

Inspection planning development: An evolutionary approach using reliability engineering as a tool

This paper proposes an evolutionary approach for inspection planning which introduces various reliability engineering tools into the process and assess system trade-offs among reliability, engineering requirement, manufacturing capability and inspection cost to establish an optimal inspection plan. The examples presented in the paper illustrate some advantages and benefits of the new approach. Through the analysis, reliability and engineering impacts due to manufacturing process capability and inspection uncertainty are clearly understood; the most cost effective and efficient inspection plan can be established and associated risks are well controlled; some inspection reductions and relaxations are well justified; and design feedbacks and changes may be initiated from the analysis conclusion to further enhance reliability and reduce cost. The approach is particularly promising as global competitions and customer quality improvement expectations are rapidly increasing.

Graf, David A.↗

Use of MODIS Cloud Top Pressure to Improve Assimilation Yields of AIRS Radiances in GSI

Radiances from hyperspectral sounders such as the Atmospheric Infrared Sounder (AIRS) are routinely assimilated both globally and regionally in operational numerical weather prediction (NWP) systems using the Gridpoint Statistical Interpolation (GSI) data assimilation system. However, only thinned, cloud-free radiances from a 281-channel subset are used, so the overall percentage of these observations that are assimilated is somewhere on the order of 5%. Cloud checks are performed within GSI to determine which channels peak above cloud top; inaccuracies may lead to less assimilated radiances or introduction of biases from cloud-contaminated radiances.Relatively large footprint from AIRS may not optimally represent small-scale cloud features that might be better resolved by higher-resolution imagers like the Moderate Resolution Imaging Spectroradiometer (MODIS). Objective of this project is to "swap" the MODIS-derived cloud top pressure (CTP) for that designated by the AIRS-only quality control within GSI to test the hypothesis that better representation of cloud features will result in higher assimilated radiance yields and improved forecasts.

Zavodsky, Bradley↗

Exploring Network-Related Optimization Problems Using Quantum Heuristics

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics.One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Wang, Zhihui↗

Study network-related optimization problems using quantum alternating optimization ansatz

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics. One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Zhihui Wang↗

Simultaneous Characterization of Wildfire Smoke and Surface Properties With Imaging Spectroscopy During the FIREX-AQ Field Campaign

We introduce and evaluate an approach for the simultaneous retrieval of aerosol and surface properties from Airborne Visible/Infrared Imaging Spectrometer Classic (AVIRIS-C) data collected during wildfires. The joint National Aeronautics and Space Administration (NASA) National Oceanic and Atmospheric Administration Fire Influence on Regional to Global Environments and Air Quality field campaign took place in August 2019, and involved two aircraft and coordinated ground-based observations. The AVIRIS-C instrument acquired data from onboard NASA's high altitude ER-2 research aircraft, coincident in space and time with aerosol observations obtained from the Aerosol Robotic Network (AERONET) DRAGON mobile platform in the smoke plume downwind of the Williams Flats Fire in northern Washington in August 2019. Observations in this smoke plume were used to assess the capacity of optimal-estimation based retrievals to simultaneously estimate aerosol optical depth (AOD) and surface reflectance from Visible Shortwave Infrared (VSWIR) imaging spectroscopy. Radiative transfer modeling of the sensitivities in spectral information collected over smoke reveal the potential capacity of high spectral resolution retrievals to distinguish between sulfate and smoke aerosol models, as well as sensitivity to the aerosol size distribution. Comparison with ground-based AERONET observations demonstrates that AVIRIS-C retrievals of AOD compare favorably with direct sun AOD measurements. Our analyses suggest that spectral information collected from the full VSWIR spectral interval, not just the shortest wavelengths, enables accurate retrievals. We use this approach to continuously map both aerosols and surface reflectance at high spatial resolution across heterogeneous terrain, even under relatively high AOD conditions associated with wildfire smoke.

FIREX-AQ↗

Climate Hazard Assessment for Stakeholder Adaptation Planning in New York City

This paper describes a time-sensitive approach to climate change projections, developed as part of New York City's climate change adaptation process, that has provided decision support to stakeholders from 40 agencies, regional planning associations, and private companies. The approach optimizes production of projections given constraints faced by decision makers as they incorporate climate change into long-term planning and policy. New York City stakeholders, who are well-versed in risk management, helped pre-select the climate variables most likely to impact urban infrastructure, and requested a projection range rather than a single 'most likely' outcome. The climate projections approach is transferable to other regions and consistent with broader efforts to provide climate services, including impact, vulnerability, and adaptation information. The approach uses 16 Global Climate Models (GCMs) and three emissions scenarios to calculate monthly change factors based on 30-year average future time slices relative to a 30- year model baseline. Projecting these model mean changes onto observed station data for New York City yields dramatic changes in the frequency of extreme events such as coastal flooding and dangerous heat events. Based on these methods, the current 1-in-10 year coastal flood is projected to occur more than once every 3 years by the end of the century, and heat events are projected to approximately triple in frequency. These frequency changes are of sufficient magnitude to merit consideration in long-term adaptation planning, even though the precise changes in extreme event frequency are highly uncertain

Horton, Radley M.↗

Stacking sequence optimization of simply supported laminates with stability and strain constraints

An integer programming formulation for the design of symmetric and balanced rectangular composite laminates with simply supported boundary conditions subject to buckling and strain constraints is presented. The design variables that define the stacking sequence of the laminate are ply-identity zero-one integers. The buckling constraint is linear in terms of the ply-identity design variables, but strains are nonlinear functions of these variables. A linear approximation is developed for the strain constraints so that the problem can be solved by sequential linearization using the branch and bound algorithm. Examples of graphite-epoxy plates under biaxial compression are presented. Optimum stacking sequences obtained using the linear approximation are compared with global optimum designs obtained using a genetic search procedure.

Nagendra, S.↗

Estimating Evapotranspiration Using an Observation Based Terrestrial Water Budget

Evapotranspiration (ET) is difficult to measure at the scales of climate models and climate variability. While satellite retrieval algorithms do exist, their accuracy is limited by the sparseness of in situ observations available for calibration and validation, which themselves may be unrepresentative of 500m and larger scale satellite footprints and grid pixels. Here, we use a combination of satellite and ground-based observations to close the water budgets of seven continental scale river basins (Mackenzie, Fraser, Nelson, Mississippi, Tocantins, Danube, and Ubangi), estimating mean ET as a residual. For any river basin, ET must equal total precipitation minus net runoff minus the change in total terrestrial water storage (TWS), in order for mass to be conserved. We make use of precipitation from two global observation-based products, archived runoff data, and TWS changes from the Gravity Recovery and Climate Experiment satellite mission. We demonstrate that while uncertainty in the water budget-based estimates of monthly ET is often too large for those estimates to be useful, the uncertainty in the mean annual cycle is small enough that it is practical for evaluating other ET products. Here, we evaluate five land surface model simulations, two operational atmospheric analyses, and a recent global reanalysis product based on our results. An important outcome is that the water budget-based ET time series in two tropical river basins, one in Brazil and the other in central Africa, exhibit a weak annual cycle, which may help to resolve debate about the strength of the annual cycle of ET in such regions and how ET is constrained throughout the year. The methods described will be useful for water and energy budget studies, weather and climate model assessments, and satellite-based ET retrieval optimization.

Rodell, Matthew↗

Performance of a Bounce-Averaged Global Model of Super-Thermal Electron Transport in the Earth's Magnetic Field

In this paper, we report the results of our recent research on the application of a multiprocessor Cray T916 supercomputer in modeling super-thermal electron transport in the earth's magnetic field. In general, this mathematical model requires numerical solution of a system of partial differential equations. The code we use for this model is moderately vectorized. By using Amdahl's Law for vector processors, it can be verified that the code is about 60% vectorized on a Cray computer. Speedup factors on the order of 2.5 were obtained compared to the unvectorized code. In the following sections, we discuss the methodology of improving the code. In addition to our goal of optimizing the code for solution on the Cray computer, we had the goal of scalability in mind. Scalability combines the concepts of portabilty with near-linear speedup. Specifically, a scalable program is one whose performance is portable across many different architectures with differing numbers of processors for many different problem sizes. Though we have access to a Cray at this time, the goal was to also have code which would run well on a variety of architectures.

McGuire, Tim↗

Advanced Modeling Strategies for the Analysis of Tile-Reinforced Composite Armor

A detailed investigation of the deformation mechanisms in tile-reinforced armored components was conducted to develop the most efficient modeling strategies for the structural analysis of large components of the Composite Armored Vehicle. The limitations of conventional finite elements with respect to the analysis of tile-reinforced structures were examined, and two complementary optimal modeling strategies were developed. These strategies are element layering and the use of a tile-adhesive superelement. Element layering is a technique that uses stacks of shear deformable shell elements to obtain the proper transverse shear distributions through the thickness of the laminate. The tile-adhesive superelement consists of a statically condensed substructure model designed to take advantage of periodicity in tile placement patterns to eliminate numerical redundancies in the analysis. Both approaches can be used simultaneously to create unusually efficient models that accurately predict the global response by incorporating the correct local deformation mechanisms.

Davila, Carlos G.↗

Convair 990 transonic flow-field simulation about the forward fuselage

A three-dimensional, viscous flow code was used to calculate the transonic flow about the forebody of the Convair CV-990 (Galileo II) research aircraft stationed at NASA Ames Research Center. The computations were used to determine the location for a differential pressure system. In addition, attitude sensor placements were verified. These instruments comprise a meteorological measurement system, which will be used for global determination of three-dimensional wind data. The code solves the thin layer form of the Reynolds-averaged Navier-Stokes equations using an implicit numerical procedure. The governing equations are written in a generalized, nonorthogonal coordinate system, and are cast in a strong conservation law form. Laminar boundary layer results are presented for free stream Mach number of 0.8 and angles of attack of zero and 2 deg. Use of this computational tool reduced the development time for the location of the sensors and aided in the optimal placement on the aircraft of these devices.

Chaussee, D. S.↗