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

Data Assimilative Optimization of WSA Source Surface and Interface Radii using Particle Filtering

TheWang-Sheeley-Arge (WSA) model estimates solar wind speed and interplanetary magnetic field polarity in the inner heliosphere using global photospheric magnetic field maps. WSA employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to determine the Sun's global coronal magnetic field. The PFSS and SCS models are connected through two radial parameters, the source surface and interface radii, which specify the overlap region between the inner SCS and outer PFSS models. Though both radii values are adjustable, they have typically been fixed to 2.5 solar radii. Our work highlights how solar wind predictions improve when the radii are allowed to vary over time. Data assimilation using particle filtering (sequential Monte Carlo) is used to infer optimal values over a fixed time window. Solar wind model predictions and satellite observations are compared with a newly developed quality-of-agreement prediction metric. The agreement metric between the model and observations is assumed to correspond to the probability of the two key WSA model parameters, the source surface and interface radii, where the highest metric value implies the optimal radii.We find that the optimal particle filter values of solar radii can perform twice as well as standard values for an exploratory period during Carrington Rotation 1901, with these values also reducing nonphysical kinking effects seen in solar magnetic field lines. Data assimilation choices of input realization and time frame have implications for variation in the solar wind over time. We present this work's theoretical context and practical applications for prediction accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

Monte Carlo simulations for water adsorption in porous materials: Best practices and new insights

Technologies based on water adsorption such as water harvesting from air have tremendous potential in mitigating important global crises such as water scarcity. An important challenge to the deployment of such technologies is finding optimal adsorbent materials. Given the large materials space of available adsorbents, large-scale computational screening can be extremely helpful for this task. This work explores the methods and details associated with such screening procedures and recommends best practices. We also shed light on the limitations of traditionally used and inexpensive to compute prescreening approaches involving geometric and energetic features to predict water adsorption behavior of porous materials. Such approaches can provide general trends to predict adsorption behavior but may lead to the overlook of potentially important structures due to the complex nature of water adsorption. Finally, this study offers insights for future water adsorption simulations to facilitate the development of optimal water adsorbents.

42 ENGINEERING↗

Toward a least-effort principle for evaluating prices of elements as indicators of sustainability

In this article, we explore simplifying mechanisms in materials economics for conceptionally organizing the flows of chemical elements into society. We use the concept of rank, which plays a crucial role in information theory where power laws are involved, and is a useful measure for seeking extremal principles. Borrowing from an analogy to Zipf’s law leads to a power-law relationship between abundance and rank that—pressing the Zipf’s law analogy further—may suggest that the price paid for elemental materials optimizes efficient utilization. Factors such as weight, volume, and moles in the earth's crust or in global mining production of elements are used in our analysis. For this research, 42 elements with price and production data for the year 2010 (a relatively typical year) were considered. We observe that the role of the volume from global mining production dominates in materials economics and indicate that elements such as Te and He should be sold at prices at leas ~ 100–1000 times more than their present price for efficient utilization.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Correcting the PFNS for more consistent fission modeling

For FY20, we had a deliverable to write a report detailing efforts to simultaneously evaluate both the prompt neutron multiplicity, $\overline{ν}$, and the prompt neutron fission neutron spectrum, PFNS, using CGMF. CGMF is the LANL-developed fission fragment decay code that consistently evaporates prompt neutrons and γ rays using the Hauser-Feshbach statistical theory of compound reactions. The decay begins by constructing the initial conditions of the fission fragments, then decaying each one from the excited state by neutrons and γ rays, conserving energy, momentum, spin, and parity in each step of the emission. The initial conditions of the fragments, along with the multiplicity, energy, and direction of each emitted neutron or γ ray, are recorded, allowing for the full reconstruction of the fission event. These event histories allow us to reconstruct average quantities, as well as correlations between observables, that can be compared with experimental or evaluated data. In that initial report, although there was already a favorable comparison between $\overline{ν}$ from CGMF, experiment, and the current ENDF/B-VIII.0 evaluation, we showed that there was still significant work to be done to improve the PFNS from CGMF. Historically, the PFNS is calculated too soft by Hauser-Feshbach fission models, and CGMF is no exception. The incorrect shape presents a significant challenge in fission modeling, including for our understanding of the fission process and for our ability to consistently calculate and predict a variety of prompt fission observables (such as fission fragment initial conditions, neutron and γ-ray multiplicities and energies, and the correlations between all observables). In our companion report, we detail our success in using CGMF to evaluate $\overline{ν}$. Although not included in the optimization explicitly, we also keep the initial conditions of the fission fragments physical, along with reproducing reasonably well the neutron multiplicity distribution. As we would expect from the sensitivities calculations from, the average neutron energies change very little from the $\overline{ν}$ optimization along with the PFNS (as will be shown in Sec. 2.6). The conclusion was that the global and statistical models would have to be investigated instead of just the fission fragment initial conditions (as is sufficient for $\overline{ν}$). This report details those efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗

Perspectives on biorefineries in microbial production of fuels and chemicals

Microbes drive our complex biosphere by regulating the global ecosystem through cycling elements and energy. Humankind has barely begun leveraging this biotransformation capacity to impact global economies and ecologies. Advances in genetic engineering, molecular analysis, metabolic flux modeling, microbial consortia/biome mapping and engineering, cell-free bioproduction, artificial intelligence/machine learning and the ever expanding -omics frontiers have set the stage for paradigm changes to how humankind produces, uses, transforms, and recycles carbon and energy through microbes. Harnessing this enormous potential could drive a global bioeconomy and manage carbon at a planetary level but requires understanding and application at a grand scale across a broad range of science and engineering disciplines. The penultimate manifestation of these advances is the “bio-refinery”, which is often referenced, but is a long way from being fully developed as a global carbon management platform. Broadening the feed stocks, processing operations, and product portfolio to a sequential cascade optimizing the conversion as a whole instead of limited outputs could greatly advance deployment and stability of a bioeconomy.

09 BIOMASS FUELS↗

The modification and application of Kuo's parameterization of cumulus convection in middle latitudes

An evaluation to reveal physical insight into the involved mechanisms through a combination of the large scale heat and moisture budgets and formulation of the latent heat released was performed. To improve upon the treatment of the interaction of deep convection with the environment, modifications consisting of considerations of the large scale moisture supply and of the vertical transport of moisture and of dry static energy inside the cloud are made. Also, a two layer cloud ensemble model is combined with the modified scheme. An examination of the modified scheme is performed by combining the large scale heat and moisture budgets and the modified formulation of the latent heat released. The modified parameterization procedures are compared with the original scheme, and the results of tests of both schemes are discussed. A variational optimization scheme is developed to analyze the wind field such that the integrated continuity equation and the global boundary condition are satisfied.

Liang, W. J.↗

Global analysis of a buck regulator

Sufficient conditions for global stability of a buck regulator using a discrete control law are found. The method of paired systems and Liapunov functions are used to establish global stability and to study the convergence of the regulator. A heuristic argument is given that the optimal switching curves associated with the paired continuous systems approximate the optimal switching curves of the discrete systems.

Edwards, D. B.↗

Space shuttle propulsion parameter estimation using optimal estimation techniques

The first twelve system state variables are presented with the necessary mathematical developments for incorporating them into the filter/smoother algorithm. Other state variables, i.e., aerodynamic coefficients can be easily incorporated into the estimation algorithm, representing uncertain parameters, but for initial checkout purposes are treated as known quantities. An approach for incorporating the NASA propulsion predictive model results into the optimal estimation algorithm was identified. This approach utilizes numerical derivatives and nominal predictions within the algorithm with global iterations of the algorithm. The iterative process is terminated when the quality of the estimates provided no longer significantly improves.

Source record↗

Technoeconomic assessment of distillation heat transfer intensification using aeroelastically fluttering reeds

Dry cooling, where forced air is the heat transfer medium, is a preferable cooling method in arid locations lacking readily available process water. However, such locations often experience high ambient temperatures that limit the effectiveness of air cooling. The objective of this study is to quantify the economic and energetic benefits of heat transfer intensification via the implementation of aeroelastically fluttering reeds to the air-cooled condenser of a methanol distillation column. Condenser size and performance, regarding recovered methanol and required fan power, is evaluated across condenser operating temperatures (T cond ) from 60 to 62°C and heat transfer coefficients (U) U base –2U base for a range of inlet air temperatures based on ambient temperature data from Yuma, Arizona. Under typical design sizing, condenser capital cost was reduced by 6%–35% (1.3U base –2U base ) and nominal methanol recovery was increased from 0.26% to 0.38% (T cond = 62–60°C). At optimized condenser size, all enhanced U and T cond pairs increase methanol recovery and reduce fan power costs compared to the optimal Ubase reference. Overall, using enhanced heat transfer to maintain condenser temperature under a wider range of inlet conditions, rather than to reduce operation temperature, produces more favorable performance. Methanol price is not a determining factor in which pairs are profitable. Analysis was repeated for a global warming scenario, revealing more valuable improvements under elevated temperatures. Energy savings from condenser improvement to a methanol production system are not significant with respect to an optimized conventional system. Furthermore, unit economic and energetic incentives suggest implementation of fluttering reeds may be justified in other applications.

42 ENGINEERING↗

Preliminary Design Optimization For A Supersonic Turbine For Rocket Propulsion

In this study, we present a method for optimizing, at the preliminary design level, a supersonic turbine for rocket propulsion system application. Single-, two- and three-stage turbines are considered with the number of design variables increasing from 6 to 11 then to 15, in accordance with the number of stages. Due to its global nature and flexibility in handling different types of information, the response surface methodology (RSM) is applied in the present study. A major goal of the present Optimization effort is to balance the desire of maximizing aerodynamic performance and minimizing weight. To ascertain required predictive capability of the RSM, a two-level domain refinement approach has been adopted. The accuracy of the predicted optimal design points based on this strategy is shown to he satisfactory. Our investigation indicates that the efficiency rises quickly from single stage to 2 stages but that the increase is much less pronounced with 3 stages. A 1-stage turbine performs poorly under the engine balance boundary condition. A portion of fluid kinetic energy is lost at the turbine discharge of the 1-stage design due to high stage pressure ratio and high-energy content, mostly hydrogen, of the working fluid. Regarding the optimization technique, issues related to the design of experiments (DOE) has also been investigated. It is demonstrated that the criteria for selecting the data base exhibit significant impact on the efficiency and effectiveness of the construction of the response surface.

Papila, Nilay↗

Multi-objective Bayesian active learning for MeV-ultrafast electron diffraction

Ultrafast electron diffraction using MeV energy beams(MeV-UED) has enabled unprecedented scientific opportunities in the study of ultrafast structural dynamics in a variety of gas, liquid and solid state systems. Broad scientific applications usually pose different requirements for electron probe properties. Due to the complex, nonlinear and correlated nature of accelerator systems, electron beam property optimization is a time-taking process and often relies on extensive hand-tuning by experienced human operators. Algorithm based efficient online tuning strategies are highly desired. Here, we demonstrate multi-objective Bayesian active learning for speeding up online beam tuning at the SLAC MeV-UED facility. The multi-objective Bayesian optimization algorithm was used for efficiently searching the parameter space and mapping out the Pareto Fronts which give the trade-offs between key beam properties. Such scheme enables an unprecedented overview of the global behavior of the experimental system and takes a significantly smaller number of measurements compared with traditional methods such as a grid scan. This methodology can be applied in other experimental scenarios that require simultaneously optimizing multiple objectives by explorations in high dimensional, nonlinear and correlated systems.

43 PARTICLE ACCELERATORS↗

Improving Prediction of Peroxide Value of Edible Oils Using Regularized Regression Models

We present four unique prediction techniques, combined with multiple data pre-processing methods, utilizing a wide range of both oil types and oil peroxide values (PV) as well as incorporating natural aging for peroxide creation. Samples were PV assayed using a standard starch titration method, AOCS Method Cd 8-53, and used as a verified reference method for PV determination. Near-infrared (NIR) spectra were collected from each sample in two unique optical pathlengths (OPLs), 2 and 24 mm, then fused into a third distinct set. All three sets were used in partial least squares (PLS) regression, ridge regression, LASSO regression, and elastic net regression model calculation. While no individual regression model was established as the best, global models for each regression type and pre-processing method show good agreement between all regression types when performed in their optimal scenarios. Furthermore, small spectral window size boxcar averaging shows prediction accuracy improvements for edible oil PVs. Best-performing models for each regression type are: PLS regression, 25 point boxcar window fused OPL spectral information RMSEP = 2.50; ridge regression, 5 point boxcar window, 24 mm OPL, RMSEP = 2.20; LASSO raw spectral information, 24 mm OPL, RMSEP = 1.80; and elastic net, 10 point boxcar window, 24 mm OPL, RMSEP = 1.91. The results show promising advancements in the development of a full global model for PV determination of edible oils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global trajectory targeting via computer graphics

A technique is described in which the two-point boundary value problem (TPBVP) may be solved with the aid of interactive computer graphics. The particular TPBVP considered is the optimal electric propulsion space trajectory problem. An appropriate two-dimensional projection of the TPBVP mapping, or trajectory, is displayed on the computer's television screen, and a man-in-the-loop varies selected trajectory starting conditions in the fashion of a nonlinear walk until the viewed trajectory endpoint lies near a displayed target. Once global targeting is accomplished in this manner, program internal logic can easily handle local targeting to strongly solve the TPBVP.

Mann, F. I.↗

Design optimization of integrated cooling inserts in modular Fischer-Tropsch reactors

Sustainable production of liquid fuels and feedstocks from atmospheric CO 2 through carbon recycling technologies is necessary to broaden decarbonization efforts and further reduce global emissions. Fischer–Tropsch (FT) reactors can address this challenge by providing storable, high-value liquid hydrocarbon fuels and feedstocks from syngas generated through reductive CO 2 utilization. FT technology, however, is most cost effective at large-scale, fixed-site plants and does not effectively address the smaller, globally distributed CO 2 point sources. Alternatively, smaller, modular reactors can be composed together to the specific scale of the emission source, yielding a flexible solution that enables more widespread deployment. In these modular reactors, thermal management using a finned cooling insert embedded within the catalyst matrix is critical to maintaining the performance and viability. Thus, in this work, topology optimization is used to determine the cooling insert geometry that maximizes reactor productivity while preventing auto-thermal runaway. Optimal designs are generated for varying number of constituent fins and over a range of maximum operating temperatures. Constraints including minimum feature length-scales and prescribed cooling insert material are imposed on the designs to enhance manufacturability. The impact of design features such as insert tapering and increased length scale hierarchy is automatically revealed by the systematic design framework employed. Here, this approach thus generates novel optimal geometries while automating, accelerating, and enhancing the design process compared to traditional heuristic approaches for cooling insert design in modular FT reactors.

30 DIRECT ENERGY CONVERSION↗

Repository of HydroSMADE: Hydropower Site-level Monthly Availability Data Ensemble for 1950-2100 at Existing and Potential Global Sites

This repository presents HydroSMADE—Hydropower Site-level Monthly Availability Data Ensemble, a new open dataset that provides monthly hydropower availability for 1,593 existing and 124,333 potential sites worldwide over the period 1950–2100. The dataset is generated by using a global hydrologic model (Xanthos) with explicit representation of hydropower operation. Specifically, HydroSMADE distinguishes between storage and diversion sites, applies optimized operating rules, and incorporates site-specific characteristics such as generation capacity, maximum turbine flow, and reservoir storage. Driven by bias-corrected meteorological inputs, the data is provided for 30 alternative future scenarios. The scenarios consist of the full factorial combination of three standard CMIP6 atmospheric forcing pathways (SSP1-2.6, SSP3-7.0, and SSP5-8.5) and ten CMIP6 General Circulation Models (GCMs): GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, EC-Earth3, CanESM5, MIROC6, CNRM-ESM2-1, UKESM1-0-LL, and CNRM-CM6-1. The repository contains a total of 122 files: a text file (readme.txt) containing a brief description of the included data, a CSV file containing site attributes, and the remaining 120 files (in CSV) containing site-level monthly hydropower availability. Example Jupyter Notebooks to explore the HydroSMADE dataset are available on GitHub at https://github.com/kamal0013/HydroSMADE More details on the methods and technical validation of HydroSMADE are available in the following paper by the same authors: Chowdhury, A. K., Abeshu, G. W., Zhao, M., Wild, T. B., Hassan, N., Ying, Z., Kim, G. J., Matthew, B., Jonathan, L., & Li, H.-Y. (Submitted). Hydropower Site-level Monthly Availability Data Ensemble for 1950-2100 at Existing and Potential Global Sites.

Existing and Potential Sites↗

Optimal message log reclamation for independent checkpointing

Independent (uncoordinated) check pointing for parallel and distributed systems allows maximum process autonomy but suffers from possible domino effects and the associated storage space overhead for maintaining multiple checkpoints and message logs. In most research on check pointing and recovery, it was assumed that only the checkpoints and message logs older than the global recovery line can be discarded. It is shown how recovery line transformation and decomposition can be applied to the problem of efficiently identifying all discardable message logs, thereby achieving optimal garbage collection. Communication trace-driven simulation for several parallel programs is used to show the benefits of the proposed algorithm for message log reclamation.

Wang, Yi-Min↗