Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maximization”

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,351 records · Page 75

Design optimization of an ethanol heavy-duty engine using design of experiments and bayesian optimization

Diesel-fueled engines still hold a large market share in the medium and heavy-duty transportation sector. However, the increase in fossil fuel prices and the strict emission regulations are leading engine manufacturers to seek cleaner alternatives without a compromise in performance. Alcohol-based fuels, such as ethanol, offer a promising alternative to diesel fuel in meeting regulatory demands. Ethanol provides cleaner combustion and lower levels of soot due to its chemical properties, in particular its lower level of carbon content. In addition, the stoichiometric operating conditions of alcohol fueled engines enable the mitigation of NOx emissions in aftertreatment stage. With the promise of retrofitting diesel engines to run on ethanol to reduce emissions, the thermal efficiency of these engines remains the primary optimization target. In order to find the optimal ethanol-fueled engine design that maximizes the thermal efficiency, a large design space needs to be investigated using engineering tools. In this study, previous research by the authors on optimizing the design of a single-cylinder ethanol-fueled engine was extended to explore the design space for a heavy-duty multi-cylinder engine configuration. A heavy-duty engine setup with multiple operating conditions at different engine speeds and loads were considered. A design optimization analysis was performed to identify the potential designs that maximize the indicated thermal efficiency in an ethanol-fueled compression ignition engine. First, a computational fluid dynamics (CFD) model of the engine was validated using experimental data for four drive cycle points. Using a design of experiments (DoE) approach and a parameterized piston bowl geometry, the model was then exercised to explore the relationship among geometric features of the piston bowl and spray targeting angle and indicated thermal efficiency across all tested operating conditions. After evaluating 165~candidate designs, a piston bowl geometry was identified that yielded an increase between 1.3 to 2.2 percentage points in indicated thermal efficiency for all tested conditions, while satisfying the operational design constraints for peak pressure and maximum pressure rise rate. The increased performance was attributed to enhanced mixing that led to the formation of a more homogeneous distribution of in-cylinder temperature and equivalence ratio, higher combustion temperatures, and shorter combustion duration. Finally, a Bayesian optimization (BOpt) analysis was employed to find the optimal piston bowl geometry with a fixed spray injector angle for one of the operating conditions. Using BOpt, a piston candidate was identified that resulted in a 1.9~percentage point increase in thermal efficiency from the baseline design, yet only required 65\% of the design samples investigated using the DoE approach.

Tekgul, Bulut↗

Accelerated CO2 Storage Optimization Using Multi-Resolution Fourier Neural Operator at the Illinois Basin Decatur Project (IBDP)

This paper presents a deep learning-based approach for optimizing CO2 injection in carbon capture and storage (CCS) operations. We developed a multi-resolution machine learning model to significantly reduce data generation costs. Utilizing this proxy model, we implemented a multi-objective genetic algorithm to optimize well control during the CO2 injection process. The proposed approach was applied to the Illinois Basin Decatur Project (IBDP), successfully optimizing the CO2 injection schedule based on three key objectives: maximizing the amount of CO2 stored, maximizing sweep efficiency, and minimizing pressure increase. The use of the proxy model accelerated the optimization workflow by two orders of magnitude, while the cost of data generation for the proxy model was reduced by 90% by utilizing a coarse-scale model.

accelerated CO2 storage optimization↗

DEVELOPMENT AND APPLICATION OF RISK ANALYSIS TOOLKIT FOR PLANT RESOURCE OPTIMIZATION

This paper presents the development of methods and tools that are being designed to optimize plant operations (e.g., maintenance/replacement schedules and optimal maintenance postures for plant components) in a manner that is more cost effective than current approaches and makes better use of available component health and cost data. These methods include both data- and model-based optimization methods. Model-based optimization methods directly include reliability and cost models to determine an optimal plant operational strategy. We consider gradient-based and evolutionary (based on genetic algorithms) optimization methods. The second class of methods target more specific use cases (e.g., project schedule optimization) and are not based on reliability models directly, but they require specific component reliability and cost data. This class of methods is based on variants of the knapsack problem with an aim to determine an optimal project schedule that maximizes the overall NPV. This paper also presents multi-objective methods designed to identify an optimal maintenance posture based on a Pareto frontier analysis. Rather than dictating the “right” tradeoff (i.e., identify the absolute best posture), we show how it is possible to perform a trade space exploration approach (i.e., identify value and costs of several postures and let the analysis account for desired value and cost metrics). This is performed by identifying maintenance postures that maximize value (e.g., system availability) and minimize operational costs, i.e., the Pareto frontier in a value-cost trade space. For all these methods we present detailed applicative examples that show their validity from a decision-making perspective.

97 - MATHEMATICS AND COMPUTING↗

Biconjugate Impedance Matching Control of the HERO WEC: Preprint

Effective control strategies are essential for maximizing the power output of wave energy converters (WECs). In this study, we explore the methodology of port-based modelling which can be used for bi-conjugate impedance matching, a type of reactive control that involves maximizing power transfer by setting the load (PTO) impedance to be the complex conjugate of the source impedance (WEC). Despite its foundational role in electrical engineering, this method has seen limited application in wave energy and has not yet been demonstrated in wave-powered desalination systems. We represent our complex multidisciplinary system as an equivalent circuit to explore frequency-based relations. The HERO WEC, a small-scale, modular point absorber (PA) designed for reverse osmosis desalination is used as a case study to demonstrate port-based modelling for control applications. Improving the performance of wave-powered water desalination systems could play a critical role in addressing freshwater scarcity for coastal communities, especially those in remote or disaster-prone regions.

16 TIDAL AND WAVE POWER↗

Optimization for Bioenergy Systems

The Sustainable Aviation Fuel (SAF) Grand Challenge (Langholtz, 2024 ) seeks to generate 35 billion gallons of SAF each year by 2050, with corn stover, an agricultural byproduct, playing a key role as a feedstock. This study develops an optimization framework to enhance the quality and quantity of corn stover while ensuring economic and environmental viability. Using the Decision Support System for Agrotechnology Transfer (DSSAT) crop model, we simulate the effects of cover crops on rotation yield, soil moisture balance, and nitrogen cycling across diverse climates and soils. The model outputs, including yield data and soil quality changes, inform a Mixed-Integer Linear Programming (MILP) optimization model. This model aims to maximize economic and environmental returns by incorporating production costs, direct and indirect income, and environmental incentives. The optimization model evaluates 280 agriculture management plans composed of various crop management strategies, including corn stover removal rates, cover crop adoption, and fertilization practices. It seeks to identify the optimal combination of crop and tillage decisions for each subfield, maximizing profits while enhancing soil carbon sequestration and reducing greenhouse gas emissions. Outputs include detailed subfield locations, optimal management plans, and profits per hectare and per acre, allowing for comparison with literature values on farm profits. This study provides a robust optimization framework supporting the SAF Grand Challenge by proposing economically viable and environmentally sustainable strategies for corn stover utilization. The findings highlight corn stover's potential as a sustainable feedstock for SAF production, offering practical solutions to enhance its quality and quantity while maintaining soil health. Idaho is used as a case study to demonstrate the framework's applicability and effectiveness in real-world scenarios. Langholtz, M. H., Davis, M., Hellwinckel, C., De La Torre Ugarte, D., Efroymson, R., Jacobson, R., Milbrandt, A., Coleman, A., Davis, R., Kline, K. L., Badgett, A., Curran, S., Schmidt, E., Theiss, T., Fried, J., English, B., Lambert, L., Cook, H., Field, J., ... Walker, L. (2024). 2023 Billion-Ton Report: An Assessment of U.S. Renewable Carbon Resources. https://doi.org/10.2172/2441098 DSSAT Foundation. (2025). Decision Support System for Agrotechnology Transfer (DSSAT). Retrieved from https://dssat.net/

09 - BIOMASS FUELS↗

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗

3D LiMn 2 O 4 nThin Film Deposited by ALD: A Road toward High-Capacity Electrode for 3D Li-Ion Microbatteries

Miniaturized electronics suffer from a lack of energy autonomy. In that context, the fabrication of lithium-ion solid-state microbatteries with high performance is mandatory for powering the next generation of portable electronic devices. Here, the fabrication of a thin film positive electrode for 3D Li-ion microbatteries made by the atomic layer deposition (ALD) method and in situ lithiation step is demonstrated. The 3D electrodes based on spinel LiMn 2 O 4 films operate at high working potential (4.1 V vs Li/Li + ) and are capable of delivering a remarkable surface capacity (approximate to 180 mu Ah cm -2 ) at low C-rate while maintaining more than 40 mu Ah cm -2 at C/2 (time constant = 2 h). Both the thickness of the electrode material and the 3D gain of the template are carefully tuned to maximize the electrode performance. Further, advanced characterization techniques such as transmission electron and X-ray transmission microscopies are proposed as perfect tools to study the conformality of the deposited films and the interfaces between each layer: no interdiffusion or segregation are observed. This work represents a major issue towards the fabrication of 3D-lithiated electrode by ALD-without any prelithiation step by electrochemical technique-making it an attractive solution for the fabrication of 3D Li-ion solid-state microbatteries with semiconductor processing methods.

3D Li-ion micro-batteries↗

An alternative approach for MLE calculation in nonlinear continuous dynamic systems

As an important metric to tell if a nonlinear dynamic system has a singular attractor or divergent trajectory, the maximal Lyapunov exponent (MLE) can be calculated from either system models or time series of state variable measurement. However, in the real world, due to inaccurate models, measurement noise, and the fact that sometimes state variables cannot be measured directly, it is very difficult to get an accurate MLE, which limits its application in, for example, in prediction of a nonlinear physical system (e. g. power systems, etc.) behavior. To overcome these factors, this paper proposed a trajectory-estimation-based MLE calculation approach. The proposed approach addressed how to calculate the MLE when state variables cannot be accessed directly, and uncertainties in system models, as well as noise in measurements. The simulation results shows that the proposed approach is able to handle well the nonlinear measurement functions between state variables and measurements, and get better results than pure model-based approaches or measurement-based approaches in front of measurement noise and model uncertainties.

maximal Lyapunov exponents, nonlinear differential↗

Re-evaluating probable maximum precipitation estimates: sensitivity to transposition domains and storm rotation using modern datasets

This study examines the sensitivity of Probable Maximum Precipitation (PMP) estimates to key methodological decisions embedded in the legacy approach adopted in the U.S. National Weather Service Hydrometeorological Reports No. 51 and No. 52. Although widely used for infrastructure design and risk regulation, fundamental aspects of PMP estimation—such as storm sample size, transposition domain, maximization procedures, and storm rotation—remain poorly constrained and lack formal guidance. Using the Red Rock watershed in Iowa as a case study, and leveraging the 2002–2023 NOAA Analysis of Record for Calibration (AORC) precipitation dataset, we systematically evaluate how each methodological choice, individually and in combination, influences PMP estimates. Our findings demonstrate that PMP is not a fixed physical upper bound but rather a modeling construct shaped heavily by user-defined assumptions. Notably, PMP values derived from modern gridded rainfall datasets can be substantially higher than the legacy estimate used in the original spillway design for Red Rock Dam. Decisions regarding storm sample size, domain extent, climatological window, and particularly storm rotation all contributed to higher PMP estimates. Storm rotation alone—a loosely constrained element in the current PMP practice—can amplify PMP by more than 25%. These results reveal the lack of standardized bounds in current PMP workflows and the need for systematic sensitivity and uncertainty analysis. As PMP estimation shifts toward probabilistic approaches, incorporating physically meaningful storm attributes will be key to developing more transparent, defensible methods for dam safety and climate-resilient infrastructure.

Probable maximum precipitation↗

Energy based fracture initiation criterion for strain-crystallizing rubber-like materials with pre-existing cracks

Fracture prediction is indispensable for polymers, like rubbers, which have a broad range of applications mainly due to their high extensibility. The phenomenon known as strain-induced crystallization further contributes to the fracture toughness of certain rubbers. Here in this study, a criterion based on internal bond energy, incorporating the effects of crystallization, is proposed to predict fracture initiation in rubber-like materials with pre-existing cracks. First, a multi-scale mechanical model is developed for characterizing the behavior of rubber when subjected to both uniaxial and biaxial deformation states. At the microscale, both the amorphous and crystalline chain segments are modeled as elastic in order to consider the energy contribution by the molecular bond distortions. This internal energy is considered along with the entropic and crystalline free energy for each chain. In the chain model, the effects of loading condition and the relative orientation of a chain on its crystallinity are taken into account. At the macroscopic scale, an existing crystallinity distribution function is adapted and a mixed finite element formulation with an augmented Lagrangian multiplier is utilized to impose the incompressibility constraint. A non-affine maximal advance path constraint based homogenization model is utilized for bridging the two scales. Its potential to account for anisotropy in the stretched network compels the model to be preferable due to its physical significance, for the purpose of fracture modeling. The rigidity of the crystallites is accounted for by proposing a crystallite distortion energy in addition to the critical bond dissociation energy, for fracture initiation to occur. The model is validated by comparison with existing experimental results for both crystallizing and non-crystallizing rubbers. In addition to its potential to predict the material behavior when subjected to uniaxial and biaxial loading, the capability of the model to quantitatively estimate the effect of crystallization on fracture initiation is also verified.

36 MATERIALS SCIENCE↗

A biradical-tagged phospholipid as a polarizing agent for solid-state MAS Dynamic Nuclear Polarization NMR of membrane proteins

We report a novel Dynamic Nuclear Polarization (DNP) NMR polarizing agent ToSMTSL-PTE representing a phospholipid with a biradical TOTAPOL tethered to the polar head group has been synthesized, characterized, and employed to enhance solid-state Nuclear Magnetic Resonance (SSNMR) signal of a lipid-reconstituted integral membrane protein proteorhodopsin (PR). A matrix-free PR formulation for DNP improved the absolute sensitivity of NMR signal by a factor of ca. 4 compared to a conventional preparation with TOTAPOL dispersed in a glassy glycerol/water matrix. DNP enhancements measured at 400 MHz/263 GHz and 600 MHz/395 GHz showed a strong field dependence but remained moderate at both fields, and comparable to those obtained for PR covalently modified with ToSMTSL. Additional continuous wave (CW) X-band electron paramagnetic resonance (EPR) experiments with ToSMTSL-PTE in solutions and in lipid bilayers revealed that an unfavorable conformational change of the linker connecting mononitroxides could be one of the reasons for moderate DNP enhancements. Further, differential scanning calorimetry (DSC) and CW EPR experiments indicated an inhomogeneous distribution and/or a possibility of a partial aggregation of ToSMTSL-PTE in DMPC:DMPA bilayers when the concentration of the polarizing agent was increased to 20 mol% to maximize the DNP enhancement. Thus, conformational changes and an inhomogeneous distribution of the lipid-based biradicals in lipid bilayers emerged as important factors to consider for further development of this matrix-free approach for DNP of membrane proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the Relationship of Microstructure and Conductivity in Metal Halide Perovskites via Active Learning-Driven Automated Scanning Probe Microscopy

Electronic transport and hysteresis in metal halide perovskites (MHPs) are key to the applications in photovoltaics, light emitting devices, and light and chemical sensors. These phenomena are strongly affected by the materials microstructure including grain boundaries, ferroic domain walls, and secondary phase inclusions. Here, we demonstrate an active machine learning framework for “driving” an automated scanning probe microscope (SPM) to discover the microstructures responsible for specific aspects of transport behavior in MHPs. In our setup, the microscope can discover the microstructural elements that maximize the onset of conduction, hysteresis, or any other characteristic that can be derived from a set of current–voltage spectra. Furthermore, this approach opens new opportunities for exploring the origins of materials functionality in complex materials by SPM and can be integrated with other characterization techniques either before (prior knowledge) or after (identification of locations of interest for detail studies) functional probing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Orientation and Stability of Asymmetric Magnetic Reconnection X Line

Here, the orientation and stability of the reconnection x line in asymmetric geometry is studied using three-dimensional (3-D) particle-in-cell simulations. We initiate reconnection at the center of a large simulation domain to minimize the boundary effect. The resulting x line has sufficient freedom to develop along an optimal orientation, and it remains laminar. Companion 2-D simulations indicate that this x line orientation maximizes the reconnection rate. The divergence of the nongyrotropic pressure tensor breaks the frozen-in condition, consistent with its 2-D counterpart. We then design 3-D simulations with one dimension being short to fix the x line orientation but long enough to allow the growth of the fastest growing oblique tearing modes. This numerical experiment suggests that reconnection tends to radiate secondary oblique tearing modes if it is externally (globally) forced to proceed along an orientation not favored by the local physics. Finally, the development of oblique structure easily leads to turbulence inside small periodic systems.

3-D magnetic reconnection↗

Exploration of lattice Hamiltonians for functional and structural discovery via Gaussian process-based exploration–exploitation

Statistical physics models ranging from simple lattice to complex quantum Hamiltonians are one of the mainstays of modern physics that have allowed both decades of scientific discovery and provided a universal framework to understand a broad range of phenomena from alloying to frustrated and phase separated materials to quantum systems. Traditionally, exploration of the phase diagrams corresponding to multidimensional parameter spaces of Hamiltonians was performed using a combination of basic physical principles, analytical approximations, and extensive numerical modeling. However, exploration of complex multidimensional parameter spaces is subject to the classic dimensionality problem, and the behaviors of interest concentrated on low dimensional manifolds remain undiscovered. Here, we demonstrate that a combination of exploration and exploration–exploitation with Gaussian process modeling and Bayesian optimization allows effective exploration of the parameter space for lattice Hamiltonians and effectively maps the regions at which specific macroscopic functionalities or local structures are maximized. We argue that this approach is general and can be further extended well beyond the lattice Hamiltonians to effectively explore the parameter space of more complex off-lattice and dynamic models.

42 ENGINEERING↗

Improved microwave SQUID multiplexer readout using a kinetic-inductance traveling-wave parametric amplifier

We report on the use of a kinetic-inductance traveling-wave parametric amplifier (KITWPA) as the first amplifier in the readout chain of a microwave superconducting quantum interference device (SQUID) multiplexer (μmux). This μmux is designed to multiplex signals from arrays of low temperature detectors, such as superconducting transition-edge sensor microcalorimeters. When modulated with a periodic flux-ramp to linearize the SQUID response, the flux noise improves, on average, from 1.6 μ $\Phi$ 0 /$\sqrt{Hz}$ with the KITWPA off to 0.77 μ $\Phi$ 0 /$\sqrt{Hz}$ with the KITWPA on. When statically biasing the μmux to the maximally flux-sensitive point, the flux noise drops from 0.45 μ $\Phi$ 0 /$\sqrt{Hz}$ to 0.2 μ $\Phi$ 0 /$\sqrt{Hz}$⁠. We validate this new readout scheme by coupling a transition-edge sensor microcalorimeter to the μmux and detecting background radiation. In conclusion, the combination of μmux and KITWPA provides a variety of new capabilities including improved detector sensitivity and more efficient bandwidth utilization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cohort organized learning: clustering through agreement

In this article we describe cohort organized learning (CoOL), a method for clustering data without explicit distance or similarity computations. Herein, we will describe CoOL, derive the gradients determined by expectation maximization to train the networks, show how to monitor convergence during training and evaluate the clusters after training, and discuss a series of examples and use cases. We also discuss CoOL’s limitations and future prospects on related tasks. Because CoOL uses neural networks to estimate the clusters, it can be used to cluster any data that can be made compatible and we illustrate this on vector data and images.

clustering↗

Multi-Timescale Three-Phase Unbalanced Distribution System Operation With Variable Renewable Generations

This paper proposes a multi-timescale operational approach for a three-phase unbalanced distribution system, which contains hourly scheduling at substation level and minutes power flow operation at feeder level. In the substation scheduling model, the objective is to minimize the system cost with variable renewable generations. The given error distribution model of renewable generation is formulated as a chance constraint, and derived into a deterministic form by Gaussian mixture model (GMM) with genetic algorithm-based expectation-maximization (GAEM). In the feeder scheduling model, the system cost is further reduced with the optimal power flow (OPF) at a timescale of minutes. Considering the nonconvexity of the three-phase unbalanced OPF problem in distribution systems, the semidefinite programming (SDP) is used to relax the problem into a convex problem, and a distributed computation approach is built based on alternating direction method of multiplier (ADMM). The IEEE 123-bus distribution system, DU (University of Denver) campus distribution system, and the IEEE 8500-bus distribution system are used as the test bench for the proposed approach. And the numerical results demonstrate the effectiveness and validity of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗