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,729 records · Page 96

Modeling dust as component minerals in the Community Atmosphere Model: development of framework and impact on radiative forcing

The mineralogy of desert dust is important due to its effect on radiation, clouds and biogeochemical cycling of trace nutrients. This study presents the simulation of dust radiative forcing as a function of both mineral composition and size at the global scale, using mineral soil maps for estimating emissions. Externally mixed mineral aerosols in the bulk aerosol module in the Community Atmosphere Model version 4 (CAM4) and internally mixed mineral aerosols in the modal aerosol module in the Community Atmosphere Model version 5.1 (CAM5) embedded in the Community Earth System Model version 1.0.5 (CESM) are speciated into common mineral components in place of total dust. The simulations with mineralogy are compared to available observations of mineral atmospheric distribution and deposition along with observations of clear-sky radiative forcing efficiency. Based on these simulations, we estimate the all-sky direct radiative forcing at the top of the atmosphere as + 0.05 Wm −2 for both CAM4 and CAM5 simulations with mineralogy. We compare this to the radiative forcing from simulations of dust in release versions of CAM4 and CAM5 (+0.08 and +0.17 Wm −2 ) and of dust with optimized optical properties, wet scavenging and particle size distribution in CAM4 and CAM5, −0.05 and −0.17 Wm −2 , respectively. The ability to correctly include the mineralogy of dust in climate models is hindered by its spatial and temporal variability as well as insufficient global in situ observations, incomplete and uncertain source mineralogies and the uncertainties associated with data retrieved from remote sensing methods.

54 ENVIRONMENTAL SCIENCES↗

Automated Resonance Fitting for Nuclear Data Evaluation

Global and national efforts to deliver high-quality nuclear data to users have a wide-ranging impact, affecting applications in national security, reactor operations, basic science, medicine, and more. Cross section evaluation is a major part of this effort, combining theory and experimentation to produce recommended values and uncertainties for reaction probabilities. Resonance region evaluation is a specialized type of nuclear data evaluation that can require significant manual effort and months of time from expert scientists. In this article, non-convex non-linear optimization methods are combined with concepts of inferential statistics to infer a resonance model from experimental data in an automated manner that is not dependent on prior evaluation(s). This methodology aims to enhance the workflow of a resonance evaluator by minimizing time, effort, and the potential for bias from prior assumptions, while enhancing reproducibility and documentation, thereby addressing well-known challenges in the field.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assimilating All-Sky Microwave Radiance Data to Improve NASA GEOS Forecasts and Analysis

The NASA Global Modeling and Assimilation Office (GMAO) has been pursuing efforts to utilize all-sky (clear+cloudy+precipitating) MW radiance data and has developed a system to assimilate all-sky GPM Microwave Imager (GMI) radiance data in the Goddard Earth Observing System (GEOS) during the last PMM funding period. The system provides additional constraints on the analysis process near the storm regions and adjusts the geophysical parameters such as precipitation, cloud, moisture, surface pressure, and wind by combining information from GMI radiance measurements and model forecasts in an optimal manner. The system proved that assimilating the GMI all-sky radiance data improve the GEOS atmospheric analyses and forecasts. This all-sky data framework has been included in the GEOS Forward Processing (FP) system since July 11, 2018 and assimilates all-sky GMI data in real-time for GEOS global analysis and forecast production at the GMAO. We are currently extending this all-sky GMI radiance data assimilation system to assimilate more all-sky MW radiance data from other sensors such as the Microwave Humidity Sounder (MHS), the Advanced Technology Microwave Sounder (ATMS), the Special Sensor Microwave Imager/Sounder (SSMIS), Advanced Microwave Scanning Radiometer 2 (AMSR2), and the Sounder for Atmospheric Profiling of Humidity in the Intertropics by Radiometery (SAPHIR) onboard the GPM constellation spacecrafts. Preliminary results from this extended all-sky system show increased benefit from cloud- and precipitation-affected MW radiances with much larger spatial and temporal coverages compared to the all-sky system assimilating GMI alone and improved GEOS forecast skills especially for lower tropospheric humidity fields.

Kim, Min-Jeong↗

A nonsmooth nonconvex optimization algorithm for two-stage optimization problems

An optimization algorithm for a group of nonsmooth nonconvex problems inspired by two-stage stochastic programming problems is proposed. The main challenges for these problems include (1) the problems lack the popular lower-type properties such as prox-regularity assumed in many nonsmooth nonconvex optimization algorithms, (2) the objective can not be analytically expressed and (3) the evaluation of function values and subgradients are computationally expensive. To address these challenges, this report first examines the properties that exist in many two-stage problems, specifically upper-C 2 objectives. Then, we show that quadratic penalty method for securityconstrained alternating current optimal power flow (SCACOPF) contingency problems can make the contingency solution functions upper-C 2 . Based on these observations, a simplified bundle algorithm that bears similarity to sequential quadratic programming (SQP) method is proposed. It is more efficient in implementation and computation compared to conventional bundle methods. Global convergence analysis of the algorithm is presented under novel and reasonable assumptions. The proposed algorithm therefore fills the gap of theoretical convergence for smoothed SCACOPF problems. The inconsistency that might arise in our treatment of the constraints are addressed through a penalty algorithm whose convergence analysis is also provided. Finally, theoretical capabilities and numerical performance of the algorithm are demonstrated through numerical examples.

97 MATHEMATICS AND COMPUTING↗

An integrated approach for the knowledge discovery in computer simulation models with a multi-dimensional parameter space

In production industries, parameter identification, sensitivity analysis and multi-dimensional visualization are vital steps in the planning process for achieving optimal designs and gaining valuable information. Sensitivity analysis and visualization can help in identifying the most-influential parameters and quantify their contribution to the model output, reduce the model complexity, and enhance the understanding of the model behavior. Typically, this requires a large number of simulations, which can be both very expensive and time consuming when the simulation models are numerically complex and the number of parameter inputs increases. There are three main constituent parts in this work. The first part is to substitute the numerical, physical model by an accurate surrogate model, the so-called metamodel. The second part includes a multi-dimensional visualization approach for the visual exploration of metamodels. In the third part, the metamodel is used to provide the two global sensitivity measures: i) the Elementary Effect for screening the parameters, and ii) the variance decomposition method for calculating the Sobol indices that quantify both the main and interaction effects. The application of the proposed approach is illustrated with an industrial application with the goal of optimizing a drilling process using a Gaussian laser beam.

97 MATHEMATICS AND COMPUTING↗

Comparative Studies of Optoelectrical Properties of Prominent PV Materials: Halide Perovskite, CdTe, and GaAs

We compare three representative high performance PV materials: halide perovskite MAPbI3, CdTe, and GaAs, in terms of photoluminescence (PL) efficiency, PL lineshape, carrier diffusion, and surface recombination and passivation, over multiple orders of photo-excitation density or carrier density appropriate for different applications. An analytic model is used to describe the excitation density dependence of PL intensity and extract the internal PL efficiency and multiple pertinent recombination parameters. A PL imaging technique is used to obtain carrier diffusion length without using a PL quencher, thus, free of unintended influence beyond pure diffusion. Our results show that perovskite samples tend to exhibit lower Shockley-Read-Hall (SRH) recombination rate in both bulk and surface, thus higher PL efficiency than the inorganic counterparts, particularly under low excitation density, even with no or preliminary surface passivation. PL lineshape and diffusion analysis indicate that there is considerable structural disordering in the perovskite materials, and thus photo-generated carriers are not in global thermal equilibrium, which in turn suppresses the nonradiative recombination. This study suggests that relatively low point-defect density, less detrimental surface recombination, and moderate structural disordering contribute to the high PV efficiency in the perovskite. This comparative photovoltaics study provides more insights into the fundamental material science and the search for optimal device designs by learning from different technologies.

14 SOLAR ENERGY↗

An Evaluation of the Effect of Network Cost Optimization for Leadership Class Supercomputers

Dragonfly-based networks are an extensively deployed network topology in large-scale high-performance computing due to their cost-effectiveness and efficiency. The US will soon have three Exascale supercomputers for leadership class workloads deployed using dragonfly networks. Compared to indirect networks of similar scale, the dragonfly network has considerably reduced cable lengths, cable counts, and switch counts, resulting in significant network cost savings for a given system size, however, these cost reductions result in reduced global minimal paths and more challenging routing. Additionally, large scale dragonfly networks often require a taper at the global link level, resulting in less bisection bandwidth than is achievable in other traditional non-blocking topologies of equivalent scale. While dragonfly networks have been extensively studied, they have yet to be fully evaluated in an extreme scale (i.e., exascale) system that targets capability workloads. In this paper, we present the results of the first large scale evaluation of a dragonfly network on an exascale system (Frontier) and compare its behavior to a similar scale fat-tree network on a previous generation TOP500 system (Summit). This evaluation aims to determine the effect of network cost optimizations by measuring a tapered topology’s impact on capability workloads. Our evaluation is based on a collection of synthetic microbenchmarks, mini-apps, and full scale applications. It compares the scaling efficiencies of each benchmark between the dragonfly-based Frontier and the fat-tree-based Summit systems. Our results show that a dragonfly network is $\sim \mathbf{3 0 \%}$ more cost efficient than a fat-tree topology, which amortizes to $\sim 3 \%$ of an exascale system cost. Furthermore, while tapered dragonfly networks impose significant tradeoffs, the impacts are not as broad as initially thought and are mostly seen in applications with global communication patterns, particularly all-to-all (e.g., FFT-based algorithms), but also local communication patterns (e.g., nearest-neighbor algorithms) that are sensitive to network performance variability.

Khan, Awais↗

Applying Constrained Bayesian Optimization to the Design of Critical Experiments

Often when planning a criticality experiment, many design configurations are iteratively investigated with a Monte Carlo transport code. The goal is that the experiment will be optimal with respect to some variable, like the fraction of fissions occurring at a certain energy range, while simultaneously being critical. Unfortunately, the Monte Carlo transport simulations are expensive, which can ultimately limit the number of configurations that can be explored. In this work, we present how Gaussian processes (GPs) can be used as a reduced-order model in a constrained Bayesian optimization (CBO) algorithm to design a criticality experiment. The GPs replace the Monte Carlo transport simulations that explore the design space. The CBO algorithm efficiently identifies new points in the design space to run the Monte Carlo transport code while respecting the criticality constraint. It does so in a manner that both improves the accuracy of the GP and finds the approximate global optimum. We demonstrate the performance of CBO with the design of a Thermal Epithermal eXperiment (TEX) for the criticality safety validation of nuclear waste models of the Hanford Tank Farm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Molecular structure calculations: A unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation

We elaborate on the theory for the variational solution of the Schroedinger equation of small atomic and molecular systems without relying on the Born-Oppenheimer paradigm. The all-particle Schroedinger equation is solved in a numerical procedure using the variational principle, Cartesian coordinates, parameterized explicitly correlated Gaussian functions with polynomial prefactors, and the global vector representation. As a result, non-relativistic energy levels and wave functions of few-particle systems can be obtained for various angular momentum, parity, and spin quantum numbers. A stochastic variational optimization of the basis function parameters facilitates the calculation of accurate energies and wave functions for the ground and some excited rotational-(vibrational-)electronic states of H{sub 2}{sup +} and H{sub 2}, three bound states of the positronium molecule, Ps{sub 2}, and the ground and two excited states of the {sup 7}Li atom.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-E-J-23: An Accurate Algorithm to Match Imperfectly Matched Images for Lung Tumor Detection Without Markers

Purpose: In order to locate lung tumors on projection images without internal markers, digitally reconstructed radiograph (DRR) is created and compared with projection images. Since lung tumors always move and their locations change on projection images while they are static on DRRs, a special DRR (background DRR) is generated based on modified anatomy from which lung tumors are removed. In addition, global discrepancies exist between DRRs and projections due to their different image originations, scattering, and noises. This adversely affects comparison accuracy. A simple but efficient comparison algorithm is reported. Methods: This method divides global images into a matrix of small tiles and similarities will be evaluated by calculating normalized cross correlation (NCC) between corresponding tiles on projections and DRRs. The tile configuration (tile locations) will be automatically optimized to keep the tumor within a single tile which has bad matching with the corresponding DRR tile. A pixel based linear transformation will be determined by linear interpolations of tile transformation results obtained during tile matching. The DRR will be transformed to the projection image level and subtracted from it. The resulting subtracted image now contains only the tumor. A DRR of the tumor is registered to the subtracted image to locate the tumor. Results: This method has been successfully applied to kV fluoro images (about 1000 images) acquired on a Vero (Brainlab) for dynamic tumor tracking on phantom studies. Radiation opaque markers are implanted and used as ground truth for tumor positions. Although, other organs and bony structures introduce strong signals superimposed on tumors at some angles, this method accurately locates tumors on every projection over 12 gantry angles. The maximum error is less than 2.6 mm while the total average error is 1.0 mm. Conclusion: This algorithm is capable of detecting tumor without markers despite strong background signals.

60 APPLIED LIFE SCIENCES↗

Integrative SP3 Workflow for Multi-PTM Proteomics Profiling (TZ-DP0)

The goal of the experiment was to demonstrate that the optimized multiplexed multi-PTM profiling workflow can comprehensively and quantitatively capture dynamic changes in protein abundance, cysteine oxidation, phosphorylation, and acetylation in cytokine-induced inflammatory stress in mouse pancreatic ß-cells. Global proteomic, redox proteomic, phosphoproteomic, and acetylomic were data collected from mouse Beta-TC-6 pancreatic Beta-cells, untreated (mock) and cytokine-treated Beta-cells at 4, 8, and 24 hours with 4 biological replicates. Samples were digested with trypsin and Lys-C, then analyzed by LC-MS/MS. Data were searched with MS-GF+, MASIC, and MaxQuant using PNNL's DMS processing pipeline.

59 BASIC BIOLOGICAL SCIENCES↗

Thermoanaerosceptrum fracticalcis gen. nov. sp. nov., a Novel Fumarate-Fermenting Microorganism From a Deep Fractured Carbonate Aquifer of the US Great Basin

Deep fractured rock ecosystems across most of North America have not been studied extensively. However, the US Great Basin, in particular the Nevada National Security Site (NNSS, formerly the Nevada Test Site), has hosted a number of influential subsurface investigations over the years. This investigation focuses on resident microbiota recovered from a hydrogeologically confined aquifer in fractured Paleozoic carbonate rocks at 863 – 923 meters below land surface. Analysis of the microorganisms living in this oligotrophic environment provides a perspective into microbial metabolic strategies required to endure prolonged hydrogeological isolation deep underground. Here we present a microbiological and physicochemical characterization of a deep continental carbonate ecosystem and describe a bacterial genus isolated from the ecosystem. Strain DRI-13 T is a strictly anaerobic, moderately thermophilic, fumarate-respiring member of the phylum Firmicutes . This bacterium grows optimally at 55°C and pH 8.0, can tolerate a concentration of 100 mM NaCl, and appears to obligately metabolize fumarate to acetate and succinate. Culture-independent 16S rRNA gene sequencing indicates a global subsurface distribution, while the closest cultured relatives of DRI-13 T are Pelotomaculum thermopropionicum (90.0% similarity) and Desulfotomaculum gibsoniae (88.0% similarity). The predominant fatty acid profile is iso-C 15:0 , C 15:0 , C 16:0 and C 14:0 . The percentage of the straight-chain fatty acid C 15:0 is a defining characteristic not present in the other closely related species. The genome is estimated to be 3,649,665 bp, composed of 87.3% coding regions with an overall average of 45.1% G + C content. Strain DRI-13 T represents a novel genus of subsurface bacterium isolated from a previously uncharacterized rock-hosted geothermal habitat. The characterization of the bacterium combined with the sequenced genome provides insights into metabolism strategies of the deep subsurface biosphere. Based on our characterization analysis we propose the name Thermoanaerosceptrum fracticalcis (DRI-13 T = DSM 100382 T = ATCC TSD-12 T ).

16S rRNA↗

A generalized method for optimization of active noise controllers in three-dimensional spaces

In this investigation, the formulation of an optimal active noise controller for harmonic, enclosed sound fields is derived using an indirect boundary element method. The optimal active noise controller is defined as the volume velocity secondary source strengths which minimize the sound pressure level at a number of discrete interior locations. In addition, the formulation of an optimal active noise controller which minimizes the free field power radiated from a generalized, distributed noise source is presented. Representative results are provided for each formulation. For the enclosure problem, results are given for global control, local control, and passive wall treatments using a rectangular cavity. For the free field radiation problem, the effects of secondary source location and quantity on active noise controller performance are investigated. The relative merits of each formulation are discussed.

Mollo, C. G.↗

Recent advances in approximation concepts for optimum structural design

The basic approximation concepts used in structural optimization are reviewed. Some of the most recent developments in that area since the introduction of the concept in the mid-seventies are discussed. The paper distinguishes between local, medium-range, and global approximations; it covers functions approximations and problem approximations. It shows that, although the lack of comparative data established on reference test cases prevents an accurate assessment, there have been significant improvements. The largest number of developments have been in the areas of local function approximations and use of intermediate variable and response quantities. It also appears that some new methodologies are emerging which could greatly benefit from the introduction of new computer architecture.

Barthelemy, Jean-Francois M.↗

Independent Constituents in Observatory Time Series in Characterizing the Secular Variation

By an ansatz of a finite set of N superimposed oscillations, observatory time series are fitted. This start model is optimized for the independence of the used 2N-dimensional functional system as a criterion. Here, the contributions of each of the dimensions to the volume of the corresponding 2N-dimensional parallelpiped of the functional space is to be evaluated. Independence and an equal standard for the constituents are achieved by the optimizing procedure varying the start values of frequencies in the fitting ansatz when quotients of successively constructed Gram-determinants are beyond a definite threshold and of nearly the same numerical value. Results for globally distributed observatories characterize the secular variation.

Webers, W.↗

Recent advances in approximation concepts for optimum structural design

The basic approximation concepts used in structural optimization are reviewed. Some of the most recent developments in that area since the introduction of the concept in the mid-seventies are discussed. The paper distinguishes between local, medium-range, and global approximations; it covers functions approximations and problem approximations. It shows that, although the lack of comparative data established on reference test cases prevents an accurate assessment, there have been significant improvements. The largest number of developments have been in the areas of local function approximations and use of intermediate variable and response quantities. It also appears that some new methodologies are emerging which could greatly benefit from the introduction of new computer architecture.

Barthelemy, J.-F. M.↗

Transition-Independent Decentralized Markov Decision Processes

There has been substantial progress with formal models for sequential decision making by individual agents using the Markov decision process (MDP). However, similar treatment of multi-agent systems is lacking. A recent complexity result, showing that solving decentralized MDPs is NEXP-hard, provides a partial explanation. To overcome this complexity barrier, we identify a general class of transition-independent decentralized MDPs that is widely applicable. The class consists of independent collaborating agents that are tied up by a global reward function that depends on both of their histories. We present a novel algorithm for solving this class of problems and examine its properties. The result is the first effective technique to solve optimally a class of decentralized MDPs. This lays the foundation for further work in this area on both exact and approximate solutions.

Becker, Raphen↗

Receiver Would Control Phasing of a Phased-Array Antenna

In a proposed digital signal-processing technique, a radio receiver would control the phasing of a phased-array antenna to aim the peaks of the antenna radiation pattern toward desired signal sources while aiming the nulls of the pattern toward interfering signal sources. The technique was conceived for use in a Global Positioning System (GPS) receiver, for which the desired signal sources would be GPS satellites and typical interference sources would be terrestrial objects that cause multipath propagation. The technique could also be used to optimize reception in spread-spectrum cellular-telephone and military communication systems. During reception of radio signals in a conventional phased-array antenna system, received signals at their original carrier frequencies are phase-shifted, then combined by analog circuitry. The combination signal is then subjected to down-conversion and demodulation. In a system according to the proposed technique (see figure), the signal received by each antenna would be subjected to down-conversion, spread-spectrum demodulation, and correlation; this processing would be performed separately from, and simultaneously with, similar processing of signals received by the other antenna elements. Following analog down-conversion to baseband, the signals would be digitized, and all subsequent processing would be digital. In the digital process, residual carriers would be removed and each signal would be correlated with a locally generated model pseudorandum-noise code, all following normal GPS procedure. As part of this procedure, accumulated values would be added in software and the resulting signals would be phase-shifted in software by the amounts necessary to synthesize the desired antenna directional gain pattern of peaks and nulls. The principal advantage of this technique over the conventional radio-frequency-combining technique is that the parallel digital baseband processing of the signals from the various antenna elements would be a relatively inexpensive and flexible means for exploiting the inherent multiple-peak/multiple-null aiming capability of a phased-array antenna. In the original intended GPS application, the peaks and nulls could be directed independently for each GPS signal being tracked by the GPS receiver. The technique could also be applied to other code-division multiple-access communication systems.

Dunn, Charles E.↗