Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scenario 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 253 records · Page 14

Tightest Mixed-Integer Programming Formulations for Quadratic SCUC Optimization

In this project, we developed new, tighter Mixed-Integer Programming (MIP) formulations for the combined Alternating Current (AC) Security-Constrained Unit Commitment (SCUC) and Security-Constrained Optimal Power Flow (SCOPF). The work addresses a critical challenge in power system operations: efficiently determining which generation units to commit and how to optimally dispatch them while maintaining network reliability constraints for both normal and contingency scenarios. Our efforts: 1. Advance the Understanding of SCUC/SCOPF Modeling: By introducing tighter MIP formulations and leveraging cutting-edge optimization tools (Julia/JuMP, PowerModels.jl), this project has pushed forward the state of the art in efficient power systems scheduling. 2. Enhance Technical and Economic Feasibility: The methods developed provide more accurate and potentially faster solutions to large-scale, realistic scheduling and dispatch problems in electric power systems, which can translate into improved reliability and potentially lower costs for grid operations. 3. Benefit to the Public: Greater efficiency in power system operations leads to cost savings for utilities and end-users. Improved reliability and integration of advanced modeling approaches can facilitate the adoption of clean energy resources and better accommodate uncertainties in renewable generation. Because this technology could impact bulk power markets and reliability, these innovations have far-reaching public benefits in terms of cost savings, reliability, and sustainability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗

A Sustainable Wind–Biogas Hybrid System for Remote Areas in Jordan: A Case Study of Mobile Hospital for a Zaatari Syrian Refugee Camp

Access to reliable and sustainable energy in remote areas remains a pressing global challenge, significantly affecting economic development and the quality of life. This study focuses on the implementation of fully off-grid wind–biogas hybrid power systems to address this issue, with a focus on remote healthcare camp facilities. This paper investigates the performance of a hybrid renewable energy system within the context of one of Jordan’s northern remote areas, the Zaatari Syrian Refugee Camp, assessing its efficiency and environmental impact by taking the Zaatari hospital as the case study. Simulations were conducted to evaluate system components, including wind turbines, biogas generators, and diesel generators. A comprehensive evaluation was conducted, encompassing both the operational efficiency of the system and its impact on the environment. This study also considered various scenarios (SC#), including grid availability and autonomy levels, to optimize system configurations. The techno-economic assessment employed the levelized cost of energy (LCOE) as a key performance indicator, and sensitivity analyses explored the impact of diesel costs and wind power fluctuations on the system. Additionally, environmental assessment was conducted to evaluate the environmental effects of hybrid systems, with a specific focus on reducing greenhouse gas emissions. This investigation involved an examination of emissions in three different scenarios. The results indicate that the lowest LCOE that could be achieved was 0.0734 USD/kWh in SC#1 with 72.42% autonomy, whereas achieving 100% autonomy increased the LCOE to 0.1756 USD/kWh. Additionally, the results reveal that in scenarios SC#2 and SC#3, which have a higher proportion of diesel generator usage, there were elevated levels of NO x and CO 2 emissions. Conversely, in SC#1, which lacks diesel generators, emissions were notably lower. The proposed hybrid system demonstrates its potential to provide a reliable energy supply to healthcare facilities in remote regions, emphasizing both economic feasibility and environmental benefits. These findings contribute to informed decision making for sustainable energy solutions in similar contexts, promoting healthcare accessibility and environmental sustainability.

09 BIOMASS FUELS↗

Sensitivity and Uncertainty Quantification of Transition Scenario Simulations

This report documents the first collective attempt at developing and applying capabilities to quantify uncertainties, assess parametric sensitivities, and optimize multiple parameters and metrics in fuel cycle simulations generated by the SA&I Campaign. To do this, external codes that were designed to perform sensitivity analysis and uncertainty quantification (SA&UQ) needed to be coupled to the SA&I Campaign’s nuclear fuel cycle simulators (NFCS). In FY20, two approaches were pursued: 1) coupling Cyclus to an ORNL-internal code called MOT (Metaheuristic Optimization Tool) and 2) coupling DYMOND to the opensource SA&UQ tool kit Dakota. The primary objective of having these NFCS/SA&UQ coupled capabilities is to better inform DOE-NE and other stakeholders on the results generated from the NFCS. For a given set of fuel cycle strategies, policies, and technology assumptions that make up a fuel cycle scenario, these NFCS have traditionally been used by the SA&I Campaign to provide quantitative answers in terms of year-by-year mass flows, infrastructure requirements, costs, etc. With these newly developed coupled capabilities, the SA&I Campaign can now efficiently simulate hundreds or thousands of these scenarios, sample large ranges of parameters and assumptions, and use the unique features of the SA&UQ tools to process the data. This enables providing answers with known and propagated uncertainties, determining the sensitivity of important metrics to different parameters and assumptions, quantifying how much fuel cycle and technology parameters impact each other, and producing optimized fuel cycle strategies for single and multiple variables. To demonstrate these new capabilities, the Cyclus/MOT was used to model several scenarios ranging from simple fleet retirements to transitions to advanced reactors. Specifically, for a transition scenario from LWRs to SFRs and advanced LWRs, uncertainty quantification, sensitivity analysis, and optimization studies were applied to cases involving single and multiple parameter (input) and single and multiple metric (output) variations. In addition, a similar transition scenario was modeled to demonstrate how to optimize the reprocessing capacity parameter to minimize two performance metrics while taking into account uncertainties from two other parameters. Lastly, a depletion module based on SCALE/ORIGEN was added in Cyclus to simulate the third scenario that was designed to quantify the impact of the modeling assumption that all LWR used nuclear fuel have the same burnup. The newly developed DYMOND/Dakota capability was also applied to a transition scenario from the existing fleet to small modular reactors and fast reactors. This particular scenario involves not only explicit isotopic depletion via ORIGEN-2, but also includes multirecycling and utilizing the criticality search feature to determine the fresh fuel composition of recycled fuel, a feature unique to the DYMOND NFCS. A large database of simulations were run with 4 main parameters that were sampled: start date of reprocessing, reprocessing capacity, energy demand growth rate, and advanced reactor share of the fleet. The 4 main metrics were uranium consumption, enrichment requirements, waste generation, and levelized cost of electricity using data from the Cost Basis Report. The demonstrated SA&UQ results include those that inform on how to choose parameters to avoid “failed” scenarios, Sobol’ indices that inform on the importance of various parameters individually and synergistically, and Analysis of Variance (ANOVA) studies that decompose parameter ranges into groups and informs on whether variations are statistically significant.

Feng, B.↗

Compatibility of divertor detachment and ELM suppression in DIII-D high- β p plasmas with ITER-similar shape

Abstract Integration of transient and steady-state divertor heat fluxes control with a high-performance core is necessary for future fusion reactors. In recent DIII-D high- β p experiments, divertor detachment and simultaneous edge localized mode (ELM) suppression are demonstrated while the plasma confinement quality is maintained high in ITER-similar shape. By optimizing the neon injection in high- β p scenario with ITER-similar shape, deep detachment and ELM suppression are achieved with a high-performance core ( β N ∼ 2.8, β p ∼ 2.3) at q 95 ∼ 7.5. Partial divertor detachment and suppression of large ELMs are achieved at q 95 ∼ 6. The stability analyses suggest that with low neon injection, the density pedestal becomes higher and steeper and the T i profile also increases, therefore the increased edge pressure and higher current density destabilize the Peeling-Ballooning mode (PBM), which would lead to a large ELM collapse. With strong neon gas puffing, the significantly reduced pedestal pressure and current density, due to the degraded T e pedestal, lead to the stabilization of PBM and ELMs are suppressed. For both cases, the coupling between the large radius internal transport barrier (ITB) and edge pedestal is the key reason for maintaining high global performance. The formation of large radius ITB compensates for pedestal degradation. Such results could provide an attractive scenario to well control the transient and steady-state heat flux onto the divertor plates while maintaining good plasma performance, which is an important step toward the steady-state operation of future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessing the Implications of Automated Merging Control in a Mixed and Heterogeneous Traffic Environment

Previous efforts to explore the implications of partial market penetration of connected and automated vehicles (CAVs) show a consensus on the benefits of higher market penetration rates (MPR) of vehicles enabled with connectivity and/or automation. There is, however, a level of uncertainty regarding the effects of lower market penetration rates and the consideration of heterogeneous vehicle fleets. Using VISSIM to perform microscopic traffic simulation and, vehicle simulation models, we assess the impacts of different CAVs market penetration rates on fuel consumption considering a heterogeneous traffic environment. The results show that the fuel efficiency benefits of optimal coordination control are maximized in moderate congested scenarios when the CAVs MPR exceeds 40%.

Rios Torres, Jackeline↗

An Optimization-Based HVAC Load and PV Disaggregation Methodology With Contextual Supervision

Here, this letter presents an optimization-based Heating, Ventilation, and Air Conditioning (HVAC) and PV disaggregation approach. This letter builds on the previous works of authors, which discuss HVAC disaggregation strategy for aggregated levels without requiring sub-metered information. This letter expands on previous works to discuss optimization-based HVAC disaggregation by varying PV penetration scenarios with joint disaggregation of HVAC load and PV.

Analytics↗

Variability in Wind Energy Generation across the Contiguous United States

ERA5 provides high-resolution, high-quality hourly wind speeds at 100 m and is a unique resource for quantifying temporal variability in likely wind-derived power production across the United States. Gross capacity factors (CF) in seven independent system operators (ISOs) are estimated using the location and rated power of each wind turbine, a simplified power curve, and ERA5 output from 1979 to 2018. Excluding the California ISO, the marginal probability of a calm (zero power production) is less than 0.1 in any ERA5 grid cell. When a calm occurs, the mean co-occurrence across wind-turbine-containing grid cells ranges from 0.38 to 0.39 for ISOs in the Midwest and central plains [Midcontinent (or Midwest) ISO (MISO), Southwest Power Pool (SPP), and the Electric Reliability Council of Texas (ERCOT) region], increasing to 0.54–0.58 for ISOs in the eastern United States [Pennsylvania–New Jersey–Maryland interconnection (PJM), New York ISO (NYISO), and New England ISO (NEISO)]. Periods with low gross CF have a median duration of ≤6 h, except in California, and are most likely during summer. Additionally, gross CF exhibit highest variance at periods of 1 day in ERCOT and SPP; on synoptic scales in MISO, NEISO, and NYISO; and on interannual time scales in PJM. This implies differences in optimal strategies for ensuring resilience of supply. Theoretical scenarios show adding wind energy capacity near existing wind farms is advantageous even in areas with high existing installed capacity (IC), while expanding into areas with lower IC is more beneficial to reducing ramps and the probability of gross CF falling below 20%. These results emphasize the benefits of large balancing areas and aggregation in reducing wind power variability and the likelihood of wind droughts.

17 WIND ENERGY↗

Optimal resource allocation for flexible-grid entanglement distribution networks

We use a genetic algorithm (GA) as a design aid for determining the optimal provisioning of entangled photon spectrum in flex-grid quantum networks with arbitrary numbers of channels and users. After introducing a general model for entanglement distribution based on frequency-polarization hyperentangled biphotons, we derive upper bounds on fidelity and entangled bit rate for networks comprising one-to-one user connections. Simple conditions based on user detector quality and link efficiencies are found that determine whether entanglement is possible. We successfully apply a GA to find optimal resource allocations in four different representative network scenarios and validate features of our model experimentally in a quantum local area network in deployed fiber. Our results show promise for the rapid design of large-scale entanglement distribution networks.

97 MATHEMATICS AND COMPUTING↗

Cyber Secure Sensor Network for Fossil Fuel Power Generation Assets Monitoring

The energy sector is undergoing digital transformation which is to say more and more power generation assets have become automated and connected to the internet. The connected sensors can tap into plants to monitor the health of assets and manage fleets remotely. These are just some of the benefits digitalization is bringing to the power industry. Within a power plant, control systems are no longer concerned with one system or one piece of equipment, but rather whole fleet of assets inter-connected with smart sensors which have the function of continuously monitoring and transmitting real-time operational data to operators. These connected systems will form a part of the industrial internet of things (IIoT). The big data scenarios provide benefits of performing system prognostics and optimization which is a key selling point of power plant digitalization.

20 FOSSIL-FUELED POWER PLANTS↗

Arbitrage and Capacity Firming in Coordination with Day-Ahead Bidding of a Hybrid PV Plant

A hybrid PV plant (HPP) combines a photovoltaic (PV) plant with a battery energy storage system (BESS), which is considered a promising step towards the future of renewable power plants by the U.S. Department of Energy. When the renewable penetration reaches a significant level, a hybrid PV plant can bid in as a controllable thermal plant in the future electricity market. In this study, a bidding and BESS scheduling model is proposed for the HPP. The robust optimization (RO) technique has been utilized to identify the worst-case scenario of uncertainties during the bidding process. To address the overly conservative issue of the single-stage RO, we have decoupled the BESS schedule for arbitrage and PV capacity firming by a two-stage RO formulation. By comparing the output of single-stage RO and two-stage RO, the two-stage RO bids and schedules in a more aggressive manner, which increases the income of HPP. Also, the penalty of under-generation is considered in our model so that the day-ahead bidding decision and arbitrage schedules can be adjusted based on the potential UNDER-GENERATION penalty. Because the proposed model is non-convex and contains multi-stages, the Column-and-Constraint Generation (C&CG) algorithm is applied to the model as the solution. The proposed model has shown better economic performance compared to a state-of-art single-stage bidding method in case studies.

BESS scheduling↗

Incorporation of market signals for the optimal design of post combustion carbon capture systems

Recent studies have shown that fossil generators equipped with post-combustion carbon capture (PCC) systems are needed to reduce the cost of deep decarbonization. Such generators need to be flexible and responsive to grid conditions, particularly in a high variable renewable energy (VRE) environment. In this work, we evaluate the net present value (NPV) of retrofitting an existing natural gas combined cycle (NGCC) unit with a flexible PCC system while incorporating market signals from a high VRE grid. We use our industrial partner’s NGCC configuration as representative of existing NGCC units and Svante’s rapid-temperature swing adsorption (TSA) for PCC. Because of its ability to rapidly startup/shutdown and ramp-up/ramp-down, the chosen capture technology is very attractive for load-following operations. For a given set of market signals, we formulate a two-stage stochastic multi-period optimization problem, under the price-taker assumption, to simultaneously optimize the design of the capture system and operation of the entire plant. Rigorous models for the NGCC unit, PCC system, and compression system are developed using commercial process simulators and validated with either plant or vendor data. For computational tractability, we develop surrogate/reduced-order models for use in the optimization problem. The surrogate model for the NGCC plant is constructed by linearizing the rigorous dynamic model at 75% load, while data-driven nonlinear surrogate models for the capture and compression systems are constructed using simulation data from the rigorous models. The optimization problem, formulated as a mixed integer bilinear program, is implemented in the IDAES® integrated platform and solved to global optimality using Gurobi 9.5. Using this formulation, we determine the profitability of retrofitting an existing NGCC unit with the chosen capture system for multiple regions in the U.S. under two scenarios with different carbon prices. Importantly, the results show that the optimal decision strongly depends on the region and on the carbon price, thereby demonstrating the importance of the inclusion of market signals in the design process.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Scalable Truck Charging Demand Simulation for Cost-Optimized Infrastructure Planning

This project developed a scalable, high-resolution model to simulate medium- and heavy-duty (MHD) electric truck charging demand and assess its impact on grid infrastructure. Using generative modeling, simulation, and cost optimization, the project delivered an end-to-end software pipeline and a library of 96 real-world scenarios for the Dallas–Houston megaregion. We demonstrated a modular architecture for transportation and grid modeling, implemented cost-optimized infrastructure planning methods, and quantified grid capital, operational, and environmental costs across a wide range of truck electrification scenarios. The results have been adopted by major utility stakeholders and contributed to regional planning efforts.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DEEP Solar: Data DrivEn Modeling and Analytics for Enhanced System Layer ImPlementation

Realizing the SETO 2030 mission of reducing solar energy costs to 3-5 c/kWh will require innovative enabling research on effective, cost-efficient integration of local PV within distribution systems. However, the intermittent and variable nature of PVs compels operators to impose conservative hosting capacity constraints. Given the extremely high variability of (intermittent and unpredictable) solar energy generation, relaxing the capacity constraints (which are currently around 15%) and achieving 100% or greater integration of renewables will require a fundamental transformation of the power grid via the utilization of exponentially larger amounts of AMI enabled fine-grained data. To address the challenges in increasing the penetration of renewable energy based DERs, this project envisions an Enhanced System Layer (ESL) at the distribution network level that is reliable, cost-effective and scalable to millions of Distributed Energy Resources (DERs)/devices. This includes developing: 1) Transformative and highly scalable machine learning based predictive analytics tools that plug into distribution system planning and provide real-time situational awareness at the distribution level for short and long-term operational planning. The tools will be built using novel data-driven energy models of millions of active nodes with AMI, 2) Adaptive stochastic analysis and optimization algorithms for real-time grid operations, 3) Dynamic Scenario Analysis using parallel Cloudenabled implementations with < 1 minute computational cycle times.

14 SOLAR ENERGY↗

Study of Storage Requirements and Costs for Shaping Renewables and Nuclear Energy (FY23 Summary Report)

To decarbonize electricity generation primarily by using wind and solar resources, it is likely that additional low (or zero) carbon dioxide (CO 2 ) alternatives will be required for ensuring a reliable electricity supply. This study extends the work and modeling framework that we consolidated in FY 2022, by assessing the zero-CO 2 pathways of a Texas power system in which the availability of variable renewable energy (VRE), nuclear energy, and energy storage types (i.e., battery and thermal energy storage [TES]) are considered the sole resources available for expansion. Using the Risk Analysis Virtual Environment (RAVEN) and Holistic Energy Resource Optimization Network (HERON) frameworks, this study investigates the market of two nuclear energy technologies (i.e., large light-water reactors [LWRs] and LWR-type small modular reactors [SMRs]) under two plausible storage coupling scenarios (i.e., electric coupling and direct thermal coupling). We also explore optimal environments for nuclear energy deployment, and identify key performance characteristics that make nuclear energy economically viable in areas with significant intermittent energy sources. Our analysis encompasses 24 cases. We model the least-cost grid systems, highlight the potential for a coupled LWR-TES approach, and utilize SMRs to achieve deep decarbonization in an affordable, technically feasible manner. We also examine seasonal variations in balancing electricity supply and demand, and account for varying performance, costs, and grid constraints. Overall, the findings of our modeling afford valuable insights for supporting future technology investment decisions in the energy sector.

25 ENERGY STORAGE↗

Determining spatially varying profit‐maximizing management practices for miscanthus and switchgrass production in the rainfed United States

Abstract Determining optimal management practices for the profitable production of perennial energy crops is critical for scaling up production beyond experimental levels. Although many experimental field studies have examined the effects of management practices on the performance of miscanthus and switchgrass, there are no recommendations for economically optimal nitrogen (N) application rates and how they should vary spatially and with the age of the energy crop as well as on optimal rotation age of the energy crop to maximize profits. We develop a modeling framework to determine economically optimal crop management decisions and simulate the variability under various scenarios for miscanthus and switchgrass production across 2287 counties in the rainfed United States. We find that profit‐maximizing N recommendations for these crops vary across maturity stages and regions and can increase the landowner's profits compared with a uniform N rate across ages and regions. We also find that the optimal rotation for these crops is shorter than the productive physical lifespan (15–20 and 10 years for miscanthus and switchgrass, respectively). Specifically, the N rate that maximizes the economic returns is negligible for miscanthus and 111 kg ha −1 for switchgrass production at age 2. The mean profit‐maximizing N rate increases with age for miscanthus, peaking at 151 kg ha −1 at age 11 before declining to 114 kg ha −1 at the optimal rotation age of 13 years while that for switchgrass is 150 kg ha −1 for middle‐aged stands and declines to 114 kg ha −1 at the optimal rotation of 8–9 years. We find that miscanthus is the most profitable energy crop in the northern region of the rainfed United States while switchgrass is most profitable in the south of the rainfed United States. Our findings are useful for improving assessments of the profitability of energy crops and guiding future management decisions by landowners.

54 ENVIRONMENTAL SCIENCES↗

Parametric raytracing modeling for NSTX-U scenario development with high harmonic fast waves and neutral beam injection

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker–Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k φ phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Algebraic Multigrid with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization

Large-scale contact mechanics simulations are crucial in many engineering fields such as structural design and manufacturing. In the frictionless case, contact can be modeled by minimizing an energy functional; however, these problems are often nonlinear, nonconvex, and increasingly difficult to solve as mesh resolution increases. In this work, we employ a Newton-based interior-point (IP) filter line-search method, an effective approach for large-scale constrained optimization. While this method converges rapidly, each iteration requires solving a large saddle-point linear system that becomes ill-conditioned as the optimization process converges, largely due to IP treatment of the contact constraints. Such ill-conditioning can hinder solver scalability and increase iteration counts with mesh refinement. Here, to address this, we introduce a novel preconditioner, algebraic multigrid with filtering (AMGF), tailored to the Schur complement of the saddle-point system. Building on the classical AMG solver, commonly used for elasticity, we augment it with a specialized subspace correction that filters near null space components introduced by contact interface constraints. Through theoretical analysis and numerical experiments on a range of linear and nonlinear contact problems, we demonstrate that the proposed solver achieves mesh independent convergence and maintains robustness against the ill-conditioning that notoriously plagues IP methods. These results indicate that AMGF makes contact mechanics simulations more tractable and broadens the applicability of Newton-based IP methods in challenging engineering scenarios. More broadly, AMGF is well suited for problems, optimization or otherwise, where solver performance is limited by a low-dimensional subspace, such as those arising from localized constraints, interface conditions, or model heterogeneities. This makes the method widely applicable beyond contact mechanics and constrained optimization.

Mathematics and Computing↗