Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 829 records · Page 46

Enabling Grid-Aware Market Participation of Aggregate Flexible Resources

Increasing integration of distributed energy resources (DERs) within distribution feeders provides unprecedented flexibility at the distribution-transmission interconnection. With the new FERC 2222 order, DER aggregations are allowed to participate in energy market. To enable market participation, these virtual power plants need to provide their generation cost curves. This paper proposes efficient optimization formulations and solution approaches for the characterization of hourly as well as multi-time-step generation cost curves for a distribution system with high penetration of DERs. Network and DER constraints are taken into account when deriving these cost curves, and they enable active distribution systems to bid into the electricity market. The problems of deriving linear and quadratic cost curves are formulated as robust optimization problems and tractable reformulation/solution algorithm are developed to facilitate efficient calculations. The proposed formulations and solution algorithm are validated on a realistic test feeder with high penetration of flexible resources.

aggregated distributed energy resources↗

Precision X-Band Linac Technologies for Nuclear Photonics Gamma-Ray Sources

Nuclear photonics is an emerging field of research requiring new tools, including high spectral brightness, tunable gamma-ray sources; high photon energy, ultrahigh-resolution crystal spectrometers; and novel detectors. This presentation focuses on the precision linac technology required for Compton scattering gamma-ray light sources, and on the optimization of the laser and electron beam pulse format to achieve unprecedented spectral brightness. Within this context, high-gradient X-band technology will be shown to offer optimal performance in a compact package, when used in conjunction with the appropriate pulse format, and photocathode illumination and interaction laser technologies. The nascent field of nuclear photonics is enabled by the recent maturation of new technologies, including high-gradient X-band electron acceleration, robust fiber laser systems, and hyper-dispersion CPA. Recent work has been performed at LLNL to demonstrate isotope-specific detection of shielded materials via NRF using a tunable, quasi-monochromatic Compton scattering gamma-ray source operating between 0.2 MeV and 0.9 MeV photon energy. This technique is called Fluorescence Imaging in the Nuclear Domain with Energetic Radiation (or FINDER). This work has, among other things, demonstrated the detection of {sup 7}Li shielded by Pb, utilizing gamma rays generated by a linac-driven, laser-based Compton scattering gamma-ray source developed at LLNL. Within this context, a new facility is currently under construction at LLNL, with the goal of generating tunable {gamma}-rays in the 0.5-2.5 MeV photon energy range, at a repetition rate of 120 Hz, and with a peak brightness in the 10{sup 20} photons/(s x mm{sup 2} x mrad{sup 2} x 0.1% bw).

42 ENGINEERING↗

A matheuristic for design and dispatch of a utility-connected distributed energy system

Modeling distributed power generation systems often requires complicated mathematical expressions that present challenges for commercial optimization solvers. Here, this paper presents a matheuristic to solve a mixed-integer optimization model that informs decisions regarding the design and dispatch of a utility-connected microgrid. We deploy a genetic algorithm to search the system design space and a linear program to solve the economic dispatch problem. The model is a component of a web tool that requires solutions within a few minutes. Our method yields objective function values within 5% of an exogenously produced optimal in fewer than 30 seconds for 90% of our test cases compared to only 10% of our test cases by a traditional optimization solver in the same amount of time.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A unified RANS–LES model: Computational development, accuracy and cost

Large eddy simulation (LES) is computationally extremely expensive for the investigation of wall-bounded turbulent flows at high Reynolds numbers. A way to reduce the computational cost of LES by orders of magnitude is to combine LES equations with Reynolds-averaged Navier–Stokes (RANS) equations used in the near-wall region. A large variety of such hybrid RANS–LES methods are currently in use such that there is the question of which hybrid RANS-LES method represents the optimal approach. The properties of an optimal hybrid RANS–LES model are formulated here by taking reference to fundamental properties of fluid flow equations. It is shown that unified RANS–LES models derived from an underlying stochastic turbulence model have the properties of optimal hybrid RANS–LES models. The rest of the paper is organized in two parts. First, a priori and a posteriori analyses of channel flow data are used to find the optimal computational formulation of the theoretically derived unified RANS–LES model and to show that this computational model, which is referred to as linear unified model (LUM), does also have all the properties of an optimal hybrid RANS–LES model. Second, a posteriori analyses of channel flow data are used to study the accuracy and cost features of the LUM. The following conclusions are obtained. (i) Compared to RANS, which require evidence for their predictions, the LUM has the significant advantage that the quality of predictions is relatively independent of the RANS model applied. (ii) Compared to LES, the significant advantage of the LUM is a cost reduction of high-Reynolds number simulations by a factor of 0.07Re{sup 0.46}. For coarse grids, the LUM has a significant accuracy advantage over corresponding LES. (iii) Compared to other usually applied hybrid RANS–LES models, it is shown that the LUM provides significantly improved predictions.

97 MATHEMATICS AND COMPUTING↗

Hydropower flexibility valuation tool for flow requirement evaluation

Timing of generation is becoming more and more valuable. This creates greater potential tension between environmental and power system objectives since both systems require their own flow patterns. Identifying win–win outcomes in this context requires being able to discuss the value of flexibility across stakeholder groups. This research proposes a two-stage optimization method to understand hydropower flexibility to meet both environmental and power system requirements. The tool simulates the two-settlement market process in the U.S. by maximizing revenues from both the day-ahead and real-time markets, subject to plant operational limits, regulatory flow and ramping requirements, and uncertainties associated with water availability and market prices. The model is formulated as linear programming problems and solved using IBM ILOG CPLEX optimizer. By examining a range of flow requirements, ramping constraints, and storage capacities, the proposed tool shows how to make more informed decisions to weigh the cost of specific flow requirements in the context of the overall license requirements. Results from the case study show that revenue is more sensitive to the ramping constraints than the minimum flow constraints. We also demonstrate that removing flow constraints in a dry month increases monthly revenue by up to 118%, as opposed to only 1% in a wet month. In addition, our results suggest that using a learning-based water flow forecast results in an increase of monthly revenue up to 6.4% compared with persistence forecast.

13 HYDRO ENERGY↗

Atacama Cosmology Telescope: Component-separated maps of CMB temperature and the thermal Sunyaev-Zel’dovich effect

Optimal analyses of many signals in the cosmic microwave background (CMB) require map-level extraction of individual components in the microwave sky, rather than measurements at the power spectrum level alone. To date, nearly all map-level component separation in CMB analyses has been performed exclusively using satellite data. In this paper, we implement a component separation method based on the internal linear combination (ILC) approach which we have designed to optimally account for the anisotropic noise (in the 2D Fourier domain) often found in ground-based CMB experiments. Using this method, we combine multi-frequency data from the Planck satellite and the Atacama Cosmology Telescope Polarimeter (ACTPol) to construct the first wide-area, arcminute-resolution component-separated maps (≈ 2100 sq. deg.) of the CMB temperature anisotropy and the thermal Sunyaev-Zel'dovich (tSZ) effect sourced by the inverse-Compton scattering of CMB photons off hot, ionized gas. Our ILC pipeline allows for explicit deprojection of various contaminating signals, including a modified blackbody approximation of the cosmic infrared background (CIB) spectral energy distribution. The cleaned CMB maps will be a useful resource for CMB lensing reconstruction, kinematic SZ cross-correlations, and primordial non-Gaussianity studies. Here, the tSZ maps will be used to study the pressure profiles of galaxies, groups, and clusters through cross-correlations with halo catalogs, with dust contamination controlled via CIB deprojection. The data products described in this paper are available on LAMBDA.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-linear Model Predictive Control of Wave Energy Converters with Realistic Power Take-off Configurations and Loss Model

Model Predictive Control (MPC) has been recognized as a well-developed control strategy for optimizing the performance of wave energy converters (WECs). The standard problem involves maximizing the absorbed power from oncoming waves while respecting motion constraints of the device and force constraints of the power take-off (PTO). The generated electrical power is then calculated by assuming an ideal PTO with no power conversion losses or limitations to power flow. In this article, we remove the assumption of an ideal PTO and present four different control options that reflect unique combinations of PTO control and power flow constraints. We also recast the MPC problem to maximize the generated power instead of the absorbed power by accounting for the losses in the power train. Using a WEC with a hydraulic PTO, we develop a loss model suitable for controls optimization. Lastly, we illustrate the process of setting up the non-linear MPC problem and compare the performance of competing control strategies using a common benchmark based on annual energy captured at a reference site.

16 TIDAL AND WAVE POWER↗

Nature of size effects in compact models of field effect transistors

Investigations have shown that in the local approximation (for sizes L < 100 μm), AlGaN/GaN high electron mobility transistor (HEMT) structures satisfy to all properties of chaotic systems and can be described in the language of fractal geometry of fractional dimensions. For such objects, values of their electrophysical characteristics depend on the linear sizes of the examined regions, which explain the presence of the so-called size effects—dependences of the electrophysical and instrumental characteristics on the linear sizes of the active elements of semiconductor devices. In the present work, a relationship has been established for the linear model parameters of the equivalent circuit elements of internal transistors with fractal geometry of the heteroepitaxial structure manifested through a dependence of its relative electrophysical characteristics on the linear sizes of the examined surface areas. For the HEMTs, this implies dependences of their relative static (A/mm, mA/V/mm, Ω/mm, etc.) and microwave characteristics (W/mm) on the width d of the sink-source channel and on the number of sections n that leads to a nonlinear dependence of the retrieved parameter values of equivalent circuit elements of linear internal transistor models on n and d. Thus, it has been demonstrated that the size effects in semiconductors determined by the fractal geometry must be taken into account when investigating the properties of semiconductor objects on the levels less than the local approximation limit and designing and manufacturing field effect transistors. In general, the suggested approach allows a complex of problems to be solved on designing, optimizing, and retrieving the parameters of equivalent circuits of linear and nonlinear models of not only field effect transistors but also any arbitrary semiconductor devices with nonlinear instrumental characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗

Integrate FARM with PID controllers: IES Simulation Ecosystem Control System Development

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM was designed to support the HERON software module in the solution of the optimal dispatch problem for IES units. As the result of HERON-FARM dispatch simulation, the set-point trajectories are optimized to meet constraints on both the production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and the process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.) at a coarse time resolution (every 10 or 100 seconds) over long time horizons (several days or weeks). In case the operational constraints need to be met at finer time resolution, the computational burden of HERON-FARM would linearly increase with the sampling rate, and sub-optimal solutions might be obtained. System responses characterized by overshoots and damped oscillations temporarily violating the imposed constraints might occur during abrupt power transients. In this report, a hierarchical control system architecture for the operation of the Dynamic Energy Technology and Integration Laboratory (DETAIL) facility constructed at INL was proposed. First, the preliminary studies on the proposed control strategy for operating the facility and the designed PI controllers were reviewed. In particular, the current approach for generating the set-point trajectories was studied, and its limits were identified. To this aim, the inclusion of a Supervisory Control layer embedding a modified version of the FARM algorithm for preserving the system safe operation over both long and real-time horizons was proposed. In this way, FARM would be applied twice, i.e., the original version (“FARM Validator”) aiding the solution of the power dispatch problem, and the modified version (“FARM Supervisory” coordinating the PI controllers to address the real-time control tasks. Despite the kernel of the two modules is the same algorithm, their roles, tasks, and capabilities are quite different. A detailed description of the role of FARM at addressing low-level control tasks is provided, along with tentative operational procedures for training the models embedded into the algorithm by using the collected experimental data.

42 ENGINEERING↗

Diffusion bonding of commercially pure Ni using Cu interlayer

The concentration dependence of diffusivity in a multi-component diffusion system makes it complicated to predict the concentration profiles of diffusing species. This so called chemical diffusivity can be expressed as a function of thermodynamic and kinetic data. DICTRA software can calculate the concentration profiles using appropriate mobility and thermodynamic data. It can also optimize the diffusivity data using experimental diffusivity data. Then the optimized diffusivity data is stored as mobility data which is a linear function of temperature. In this work, diffusion bonding of commercially pure Ni using Cu interlayers is reported. The mobility parameters of Ni-Cu alloy binary systems were optimized using DICTRA/Thermocalc software from the available self-, tracer and chemical diffusion coefficients. The optimized mobility parameters were used to simulate concentration profiles of Ni-Cu diffusion joints using DICTRA/Thermocalc software. The calculated and experimental concentration profiles agreed well at 1100 Degree-Sign C. Agreement between the simulated and experimental profiles was less good at 1050 Degree-Sign C due to the grain boundary contribution to the overall diffusion. - Highlights: Black-Right-Pointing-Pointer The concentration profiles of Cu in Ni-Cu diffusion joints are modeled. Black-Right-Pointing-Pointer Interdiffusion coefficients in Ni-Cu system are optimized. Black-Right-Pointing-Pointer Optimized interdiffusion coefficients are expressed as mobility parameters. Black-Right-Pointing-Pointer Simulated profiles are comparable with experimental profiles.

36 MATERIALS SCIENCE↗

Revealing Decision Conservativeness Through Inverse Distributionally Robust Optimization

This paper introduces Inverse Distributionally Robust Optimization (I-DRO) as a method to infer the conservativeness level of a decision-maker, represented by the size of a Wasserstein metric-based ambiguity set, from the optimal decisions made using Forward Distributionally Robust Optimization (F-DRO). By leveraging the Karush-Kuhn-Tucker (KKT) conditions of the convex F-DRO model, we formulate I-DRO as a bi-linear program, which can be solved using off-the-shelf optimization solvers. Additionally, this formulation exhibits several advantageous properties. We demonstrate that I-DRO not only guarantees the existence and uniqueness of an optimal solution but also establishes the necessary and sufficient conditions for this optimal solution to accurately match the actual conservativeness level in F-DRO. Furthermore, we identify three extreme scenarios that may impact I-DRO effectiveness. Our case study applies F-DRO for power system scheduling under uncertainty and employs I-DRO to recover the conservativeness level of system operators. Numerical experiments based on an IEEE 5-bus system and a realistic NYISO 11-zone system demonstrate I-DRO performance in both normal and extreme scenarios. An extended version of this paper with additional analyses is available at li2024revealing.

distributionally robust optimization↗

Quantum dot spontaneous emission control in a ridge waveguide

We investigate the spontaneous emission (SE) of self-assembled InAs quantum dots (QDs) embedded in GaAs ridge waveguides that lay on a low index substrate. In thin enough waveguides, the coupling to the fundamental guided mode is vanishingly small. A pronounced anisotropy in the coupling to non-guided modes is then directly evidenced by normal-incidence photoluminescence polarization measurements. In this regime, a measurement of the QD decay rate reveals a SE inhibition by a factor up to 4. In larger wires, which ensure an optimal transverse confinement of the fundamental guided mode, the decay rate approaches the bulk value. Building on the good agreement with theoretical predictions, we infer from calculations the fraction β of SE coupled to the fundamental guided mode for some important QD excitonic complexes. For a charged exciton (isotropic in plane optical dipole), β reaches 0.61 at maximum for an on-axis QD. In the case of a purely transverse linear optical dipole, β increases up to 0.91. This optimal configuration is achievable through the selective excitation of one of the bright neutral excitons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SimCCS: An open-source tool for optimizing CO 2 capture, transport, and storage infrastructure

Commercial-scale carbon capture and storage (CCS) technology will involve deploying infrastructure on a massive and costly scale. This effort will require careful and comprehensive planning to ensure that capture locations, storage sites, and the dedicated CO 2 distribution pipelines are selected in a robust and cost-effective manner. Presented in 2009, SimCCS is an optimization model for integrated system design that enables researchers, stakeholders, and policy makers to design CCS infrastructure networks. SimCCS 2.0 is a complete, ground-up redesign that is now a portable software package, useable and shareable by the CCS research, industrial, policy, and public communities. SimCCS 2.0 integrates multiple new capabilities including a refined optimization model, novel candidate network generation techniques, and optional integration with high-performance computing platforms. Accessing user-provided CO 2 source, sink, and transportation data, SimCCS 2.0 creates candidate transportation routes and formalizes an optimization problem that determines the most cost-effective CCS system design. This optimization problem is further solved either through a high-performance computing interface, or through third-party software on a local desktop computing platform. Finally, SimCCS 2.0 employs an open-access geographic information system framework to enable analysis and visualization capabilities. SimCCS 2.0 is written in Java and is publicly available via GitHub to encourage collaboration, modification, and community development.

42 ENGINEERING↗

Performance evaluation of ZnO–CuO hetero junction solid state room temperature ethanol sensor

Graphical abstract: Sensor response (resistance) curves of time were changed from 150 ppm to 250 ppm alcohol concentration of ZnO–CuO 1:1. The response and recovery times were measured to be 62 and 83 s, respectively. The sensing material ZnO–CuO is a high potential alcohol sensor which provides a simple, rapid and highly sensitive alcohol gas sensor operating at room temperature. Highlights: ► The main advantages of the ethanol sensor are as followings. ► Novel materials ZnO–CuO ethanol sensor. ► The optimized ZnO–CuO hetero contact system. ► A good sensor response and room working temperature (save energy). -- Abstract: A semiconductor ethanol sensor was developed using ZnO–CuO and its performance was evaluated at room temperature. Hetero-junction sensor was made of ZnO–CuO nanoparticles for sensing alcohol at room temperature. Nanoparticles were prepared by hydrothermal method and optimized with different weight ratios. Sensor characteristics were linear for the concentration range of 150–250 ppm. Composite materials of ZnO–CuO were characterized using X-ray diffraction (XRD), temperature-programmed reduction (TPR) and high-resolution transmission electron microscopy (HR-TEM). ZnO–CuO (1:1) material showed maximum sensor response (S = R{sub air}/R{sub alcohol}) of 3.32 ± 0.1 toward 200 ppm of alcohol vapor at room temperature. The response and recovery times were measured to be 62 and 83 s, respectively. The linearity R{sup 2} of the sensor response was 0.9026. The sensing materials ZnO–CuO (1:1) provide a simple, rapid and highly sensitive alcohol gas sensor operating at room temperature.

36 MATERIALS SCIENCE↗

The optimal use of silicon pixel charge information for particle identification

Particle identification using the energy loss in silicon detectors is a powerful technique for probing the Standard Model (SM) as well as searching for new particles beyond the SM. Traditionally, such techniques use the truncated mean of the energy loss on multiple layers, in order to mitigate heavy tails in the charge fluctuation distribution. We show that the optimal scheme using the charge in multiple layers significantly outperforms the truncated mean. Furthermore, truncation itself does not significantly degrade performance and the optimal classifier is well-approximated by a linear combination of the truncated mean and truncated variance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Network Modeling of the U.S. Petrochemical Industry under Raw Material and Hurricane Harvey Disruptions

A geographically resolved network model of the U.S. chemical industry in 2017 is designed, and optimal chemical flows between units are calculated using linear programming. A baseline solution and three disruption scenarios (primary raw material disruptions, reported Hurricane Harvey ethylene cracker disruptions, and assumed capacity disruptions based on the Hurricane Harvey geographic storm track) are reported on to determine how the structure of the industry is modified to adapt to widespread and geographically specific disruptions. The calculated impacts of the assumed Hurricane Harvey disruption include 170 chemical units in 26 states that change production level as a result of supply chain disruptions during the storm. The systemic impact for the assumed Hurricane Harvey disruption is 19.3 million tonnes of gross chemical production. The day with the largest impact on gross chemical production shows a reduction from baseline operations of 1.3 million tonnes (42% of baseline). This model can be utilized for analysis of future disruption scenarios and to test resilience strategies, including impacts of new manufacturing configurations or technologies.

02 PETROLEUM↗