Mixed Integer Linear Programming Optimization Research Articles

Comparative Analysis of Conventional (Financial, Network) and Optimized Mixed Integer Linear Programming Scenarios for Telecom Network Investment Decision Making

The telecommunications industry in Indonesia faces numerous challenges in allocating resources for network development, including network complexity, technological changes, market competition, and customer expectations. To achieve optimal financial outcomes and customer satisfaction, investments must be well-managed, well-placed, and well-executed. This study analyses network investment portfolio selection outcomes at Telkomsel; Indonesia's largest cellular operator, comparing conventional scenarios (financial and network) with an optimized scenario using Mixed Integer Linear Programming (MILP). The research evaluates these scenarios based on total portfolio score, incremental revenue, and Internal Rate of Return (IRR). Financial indicators such as net present value (NPV), IRR, EBIT margin, incremental revenue, and network indicators such as capacity and customer satisfaction (download/upload throughput, latency, packet loss, jitter) assess investment feasibility. The MILP optimization scenario strongly correlates with incremental revenue, NPV, and IRR, indicating higher financial performance. Sites selected in the MILP optimization outperscenario formed others in total score (22% better than the financial scenario, 63% better than the network scenario), incremental revenue (3.5% better than the financial scenario, 90.2% better than the network scenario), and portfolio IRR (4% better than the financial scenario, 70% better than the network scenario).

Cross regional online food delivery: Service quality optimization and real-time order assignment

Online food delivery (OFD) represents a rapidly evolving e-business application that leverages cloud computing data centers, playing a crucial role in meeting the demands of urban lifestyles. With diverse order fulfillment features and increasing expectations for service quality, the task of effectively assigning riders for timely long-distance, cross-regional deliveries presents a significant engineering challenge. Previous studies often relied on traditional rider allocation methods that fail to account for varying capacities, or they utilized non-intelligent systems that did not adequately address fluctuating order demands and service delays. In this study, we introduce a robust Mixed Integer Linear Programming (MILP) optimization framework designed to minimize the total service time and delivery cost for cross-regional orders. This framework divides a large OFD area into multiple regions and utilizes both transfer vehicles and riders to optimize deliveries. To enhance the predictive accuracy of our model, we incorporate advanced machine learning techniques. Specifically, we employ the Long Short-Term Memory (LSTM) model to forecast regional order demands accurately, reflecting the dynamic nature of the marketplace. Additionally, Extreme Gradient Boosting (XGBoost) is tailored to dynamically predict travel times from restaurants to customer locations, facilitating more precise scheduling and resource allocation within the MILP framework. These machine learning techniques significantly bolster the MILP framework by providing detailed, accurate predictions that improve decision-making processes and adaptability to real-time conditions. Acknowledging the complexity of this optimization problem, we further enhance our approach by integrating a meta-heuristic algorithm, Adaptive Large Neighbor Search (ALNS), which efficiently assigns orders to the appropriate transfer vehicles and riders within polynomial time. Our Cross Regional Online Food Delivery (XROFD) system is meticulously designed to optimize both customer satisfaction and rider incentives. Simulation experiments confirm that the XROFD system not only reduces service times and delivery costs but also markedly enhances customer satisfaction and provides superior incentives for riders, outperforming existing state-of-the-art methods.

Optimizing supply chain networks using mixed integer linear programming (MILP)

Mixed Integer Linear Programming (MILP) has emerged as a powerful tool for optimizing complex supply chain networks. This paper explores the theoretical foundations of MILP, including the integration of integer variables and advanced solution techniques such as branch-and-bound and branch-and-cut algorithms. Through detailed modeling of production planning, network design, and transportation logistics, MILP enables companies to achieve significant cost reductions and operational efficiencies. We present case studies from retail, manufacturing, and pharmaceutical sectors to illustrate the practical applications of MILP. These examples demonstrate how MILP optimization can lead to reductions in production and inventory costs, improved customer satisfaction, and enhanced service levels. The findings underscore the value of MILP in addressing the multifaceted challenges of modern supply chain management.

Sustainable development of the water-energy-Co2 nexus in the refining sector: A stochastic multi-objective optimization under emissions trading systems

Refineries consume substantial energy and water while emitting significant CO2. Solvent-based Carbon Capture (CC) offers a viable CO2 mitigation solution but increases water and energy consumption and requires considerable investments. This research aims to overcome these barriers and improve the Iranian refining sector towards Sustainable Development Goals (SDGs) 6, 7, and 13 by mitigating CO2 emissions, producing low-carbon hydrogen, and utilizing unconventional water resources for CC plants. A comprehensive nexus approach was incorporated to account for the intricate interconnections between SDGs and refineries' water, energy, and CO2 networks. We evaluate the implementation of CO2 market regulations corresponding to seven international Emissions Trading Systems (ETSs) of the EU, UK, New Zealand (NZ), Regional Greenhouse Gas Initiative (RGGI), South Korea, California, and China in the Iranian refining sector to incentivize CC integration. The problem was formulated using Stochastic Multi-Objective Mixed-Integer Linear Programming optimization, with K-means Clustering applied to account for ETS CO2 price fluctuations and uncertainties. Results showed intense competition between economic and environmental objectives without ETS regulations, but appropriate ETS implementation could transform this conflict into a cooperative compromise. Under EU and UK scenarios, 60.3% and 67% of CO2 emissions could be avoided with fully covered capture costs, potentially leading to net economic profit through surplus allowance sales. Implementation of the other five ETSs could reduce capture costs by 18–67%, with a 16.4% CO2 avoidance.

Design and investment strategy optimization of district cooling system during the ramp-up phase

District Cooling (DC) systems, known for their efficiency and integration of renewable energy sources, encounter substantial design and investment complexities during the initial ramp-up phase. This critical phase is characterized by increasing cooling demand, necessitating a scientific approach to equipment selection and sizing that spans multiple temporal scales. A key yet frequently underestimated element in DC systems is the inclusion of backup chiller. While traditionally reserved to emergency use, its proactive interactions with regular system operations can significantly diminish operational expenses, despite the inherent complexities. Our study presents an innovative mixed-integer linear programming (MILP) optimization framework. This framework is designed to manage and optimize the intricate DC system design and investment strategy during the ramp-up phase, with a particular focus on the strategic inclusion of backup chillers. To demonstrate this framework, our research undertakes a comprehensive analysis of a real-world DC project in Kuala Lumpur. The findings show that the MILP optimization delivers substantial reductions in Life Cycle Costs (LCC) up to 3.33% across various system configurations. This is of significant value for energy efficiency and sustainable urban energy systems.

Integration of MILP and discrete-event simulation for flow shop scheduling using Benders cuts

Optimization-based scheduling in the chemical industry is highly beneficial but also highly difficult due to its combinatorial complexity. Different modeling and optimization techniques exist, each with individual strengths. We propose Benders decomposition to integrate mixed-integer linear programming (MILP) and discrete-event simulation (DES) to solve flow shop scheduling problems with makespan minimization objective. The basic idea is to generate valid Benders cuts based on sensitivity information of the DES sub problem, which can be found in the critical paths of DES solutions. For scaled literature flow shops, our approach requires at least an order of magnitude fewer iterations than a genetic algorithm and provides optimality gap information. For a real-world case study, our approach finds good solutions very quickly, making it a powerful alternative to established methods. We conclude that the Benders-DES algorithm is a promising approach to combine rigorous MILP optimization capabilities with high-fidelity DES modeling capabilities.

Towards sustainable smart cities: Integration of home energy management system for efficient energy utilization

Smart Homes (SHs) play a crucial role in the broader context of smart cities, contributing to the overall efficiency, sustainability, and quality of life for residents. This work presents a computationally efficient mathematical model for Home Energy Management Systems (HEMS) to activate the Demand Side Management (DSM) while minimizing the energy consumption costs of SHs. The proposed model is represented as a standard Mixed-Integer Linear Programming (MILP) optimization problem. The proposed HEMS provides the optimal scheduling of home appliances and plugging time for charging Electric Vehicles (EV) at home. The main objective of this optimization problem is to minimize the electricity bills while improving the load profile of the distribution system in the context of smart cities. In this study, load management for a typical household has been studied considering several electricity tariffs in Portugal. The simulation results confirm that in the presence of the proposed HEMS, the electricity bill will be effectively reduced while consumer behavior is changed. The cost savings that can be attained through only shifting loads is 5.7 %, while it can reach 67.65 % by installing a hybrid energy system at the studied smart home. The weekly cost reduction percentage due to charging the EV at home is around 16.24 % compared to charging the EV using public charging stations regardless of the traveling time and waiting at the public charging stations.

A peer-to-peer energy trading market framework for single buyer model: An integration of decentralised finance instruments and cooperative game model

Single Buyer (SB) model has gained traction in many Southeast Asia, Middle Eastern, and North African countries for its blend of centralisation and controlled competition. The semi-liberalised regulatory environment has hindered the integration of cogeneration or combined heat and power (CHP) as a prosumer within the electricity value chain, despite its pivotal role in enhancing energy efficiency. There is a disconnection in the existing literature as electric supply configurations and enablement of peer-to-peer (P2P) electricity trading are individually studied. There also lacks the studies addressing electricity as a special commodity and its compatibility with the market. This paper bridges the research gaps by proposing a P2P energy trading market framework that integrates decentralised finance (DeFi) instruments and a cooperative game model. Malaysia, which implements the SB model in its electric supply structure, is selected for a case study. The comprehensive framework develops (1) a stableswap hybrid function automated market maker (HFMM) decentralised exchange (DEX) for P2P trading, (2) a mixed-integer linear programming (MILP) optimisation model to determine the optimal operational strategy with minimised total annualised cost (TAC) of the system, and (3) an allocation scheme based on a cooperative game model to gauge the individual and coalitional cost savings allocations. Overall, this study presents a novel approach to addressing the challenges with innovating the SB model, while also fostering competitions among prosumers. The SB structure of Malaysia is briefly inspected, and the proposed framework is expected to contribute to the growth of macro-level CHP integration in the nation's energy planning.

Pricing and Replenishment Prediction Model for Vegetable Commodities Based on Mixed Integer Linear Programming Optimization XGBoost

In fresh food superstores, the replenishment volume and pricing strategy of vegetables are affected by various factors such as seasonal seasons, supply demand, purchasing costs, and profits of vegetable categories. In this paper, we use the data related to the sales flow details and wholesale prices of each commodity from July 1, 2020 to June 30, 2023 as a dataset, and build a model based on the gradient boosting tree algorithm, XGBoost algorithm [1], to predict the total daily replenishment volume and pricing strategy of each vegetable category in the coming week (July 1-7, 2023) [2]. It is also optimized using a mixed-variable 0-1 linear programming (MILP) [3] model with decision variables including a binary variable of whether to order the vegetable category or not and a continuous variable of the number of vegetable categories ordered to compute the optimal replenishment volume and pricing strategy for July 1, 2023 for each item. Through the model in this paper, under the premise of meeting the market demand for each category of vegetable goods, while keeping the control within a specific range, the superstore can get the optimal replenishment plan and pricing strategy, reasonably allocate the resources of the warehouse space, improve the market competitiveness of the superstore, and maximize the revenue.

Stochastic optimization framework for hybridization of existing offshore wind farms with wave energy and floating photovoltaic systems

Due to the complementarity of renewable energy sources, there has been a focus on technology hybridization in recent years. In the area of hybrid offshore power plants, the current research projects mostly focus on the combinational implementation of wind, solar, and wave energy technologies. Accordingly, considering the already existing offshore wind farms, there is the potential for the implementation of hybrid power plants by adding wave energy converters and floating photovoltaics. In this work, a stochastic sizing model is developed for the hybridization of existing offshore wind farms using wave energy converters and floating photovoltaics considering the export cable capacity limitation. The problem is modeled from an investor perspective to maximize the economic profits of the hybridization, while the costs and revenues regarding the existing units and the export cable are excluded. Furthermore, to tackle the uncertainties of renewable energy generation, as well as the energy price, a scenario generation method based on copula theory is proposed to consider the dependency structure between the different random variables. Altogether, the hybridization study is modeled in a mixed integer linear programming optimization framework considering the net present value of the project as the objective function. The results showed that hybrid-sources-based energy generation provided the highest economic profit in the studied cases in the different geographical locations. Furthermore, the technical specifications of the farms have also been considerably improved providing more stable energy generation, guaranteeing a minimum level of power in a high share of the time, and with a better utilization of the capacity of the cable while the curtailment of energy is maintained within the acceptable range.

Impact of tariff structures on energy community and grid operational parameters

This paper presents a mixed-integer linear programming optimization model of a renewable energy community comprised of members with local generators, battery energy storage systems, electric vehicles, and heat pumps and thermal energy storage, thus representing a local multi-energy system. The goal of the paper is to analyze the impact of different tariff structures for incentivizing energy sharing within energy communities on the distribution grid operational parameters. To achieve this, a comprehensive analysis of a hypothetical renewable energy community is conducted using the presented optimization model. Various tariff structures are evaluated, including flat and dynamic tariffs, reflecting the day-ahead market prices for energy and supply, as well as various network tariff designs, such as regulated charge reduction, time-of-use network tariffs and capacity payments. The results indicate that sharing energy within communities reduces energy losses in the grid, peak loads and reverse power flows and results in a smoother voltage profile. It is also found that tariffs that include dynamic electricity and supply costs can further increase these benefits, while including capacity payments as an incentive mechanism yields most favourable results for all analyzed criteria.

Potentials of Mixed-Integer Linear Programming (MILP)-Based Optimization for Low-Carbon Hydrogen Production and Development Pathways in China

Hydrogen (H2) is considered one of the main pillars for transforming the conventional “dark” energy system to a net-zero carbon or “green” energy system. This work reviewed the potential resources for producing low-carbon hydrogen in China, as well as the possible hydrogen production methods based on the available resources. The analysis and comparison of the levelized cost of hydrogen (LCOH) for different hydrogen production pathways, and the optimal technology mixes to produce H2 in China from 2020 to 2050 were obtained using the mixed-integer linear programming (MILP) optimization model. The results were concluded as three major ones: (a) By 2050, the LCOH of solar- and onshore-wind-powered hydrogen will reach around 70–80 $/MWh, which is lower than the current H2 price and the future low-carbon H2 price. (b) Fuel costs (>40%) and capital investments (~20%) of different hydrogen technologies are the major cost components, and also are the major direction to further reduce the hydrogen price. (c) For the optimal hydrogen technology mix under the higher renewable ratio (70%) in 2050, the installed capacities of the renewable-powered electrolysers are all more than 200 GW, and the overall LCOH is 68.46 $/MWh. This value is higher than the LCOH (62.95 $/MWh) of the scenario with higher coal gasification with carbon capture and the storage (CG-CCS) ratio (>50%). Overall, this work is the first time that hydrogen production methods in China has been discussed comprehensively, as well as the acquisition of the optimal H2 production technology mix by the MILP optimization model, which can provide guidance on future hydrogen development pathways and technology development potential in China.

Multi-objective optimization of straw-based bio-natural gas supply chains considering cost, CO2 emission, and safety

Bio-natural gas (BNG) has received recent interest as a renewable alternative to natural gas (NG). BNG supply chains (BSC), including production, transportation, and utilization, are compatible and can overlap with existing NG infrastructure. The goal is to develop a set of models that achieve the various trade-offs between economic, environmental and safety risks by combining life cycle assessment (LCA), mixed-integer linear programming (MILP) and fuzzy optimization, expanding the utilization of BSC while quantifying the risk of the whole process, so as to improve the safety awareness of the project operation. BNG technologies were considered: anaerobic digestion (AD), gasification, and syngas methanation. This model can analyze the changes in the different objectives of each module and the distribution of biomass mass flow by comparing different straw-based harvesting rates. Results indicate that AD technology possesses the lowest economic cost, and BNG as transportation fuel will generate the most benefits. The unit cost of BNG for different technologies ranges from 0.36 to 0.67 $/m³，the emission reduction is from 5.18 kg/m³ to 10.28 kg/m³，and the maximum safety risk value is 1.96. The proposed model was applied to a case study in China to illustrate the optimization of a regional straw-based BSC. The results illustrate how BNG supply chains can be systematically planned as part of an integrated decarbonization policy while mitigating safety risks that may occur from handling this fuel.

Quantum Circuit Reconstruction from Power Side-Channel Attacks on Quantum Computer Controllers

The interest in quantum computing has grown rapidly in recent years, and with it grows the importance of securing quantum circuits. A novel type of threat to quantum circuits that dedicated attackers could launch are power trace attacks. To address this threat, this paper presents first formalization and demonstration of using power traces to unlock and steal quantum circuit secrets. With access to power traces, attackers can recover information about the control pulses sent to quantum computers. From the control pulses, the gate level description of the circuits, and eventually the secret algorithms can be reverse engineered. This work demonstrates how and what information could be recovered. This work uses algebraic reconstruction from power traces to realize two new types of single trace attacks: per-channel and total power attacks. The former attack relies on per-channel measurements to perform a brute-force attack to reconstruct the quantum circuits. The latter attack performs a single-trace attack using Mixed-Integer Linear Programming optimization. Through the use of algebraic reconstruction, this work demonstrates that quantum circuit secrets can be stolen with high accuracy. Evaluation on 32 real benchmark quantum circuits shows that our technique is highly effective at reconstructing quantum circuits. The findings not only show the veracity of the potential attacks, but also the need to develop new means to protect quantum circuits from power trace attacks. Throughout this work real control pulse information from real quantum computers is used to demonstrate potential attacks based on simulation of collection of power traces.

Closed loop model predictive control of a hybrid battery-hydrogen energy storage system using mixed-integer linear programming

The derivation of an efficient operational strategy for storing intermittent renewable energies using a hybrid battery-hydrogen energy storage system is a difficult task. One approach for deriving an efficient operational strategy is using mathematical optimization in the context of model predictive control. However, mathematical optimization derives an operational strategy based on a non-exact mathematical system representation for a specified prediction horizon to optimize a specified target. Thus, the resulting operational strategies can vary depending on the optimization settings.This work focuses on evaluating potential improvements in the operational strategy for a hybrid battery-hydrogen energy storage system using mathematical optimization. To investigate the operation, a simulation model of a hybrid energy storage system and a tailor-made mixed integer linear programming optimization model of this specific system are utilized in the context of a model predictive control framework. The resulting operational strategies for different settings of the model predictive control framework are compared to a rule-based controller to show the potential benefits of model predictive control compared to a conventional approach. Furthermore, an in-depth analysis of different factors that impact the effectiveness of the model predictive controller is done. Therefore, a sensitivity analysis of the effect of different electricity demands and resource sizes on the performance relative to a rule-based controller is conducted. The model predictive controller reduced the energy consumption by at least 3.9 % and up to 17.9% compared to a rule-based controller. Finally, Pareto fronts for multi-objective optimizations with different prediction and control horizons are derived and compared to the results of a rule-based controller. A cost reduction of up to 47 % is achieved by a model predictive controller with a prediction horizon of 7 days and perfect foresight.

Comparative modeling of cost-optimal energy system flexibility for Swedish and Austrian regions

This study develops a reproducible method for estimating the cost-efficient flexibility potential of a local or regional energy system. Future scenarios that achieve ambitious climate targets and estimate the cost-efficient flexibility potential of demonstration sites were defined. Flexible potentials for energy system assessment are upscaled from the demonstration sites in Eskilstuna (Sweden) and Lower Austria (Austria). As heat pumps (HPs) and district heating (DH) are critical for future heat demand, these sites are representative types of DH networks in terms of size and integration with the electricity grid. In both regions a TIMES model is used for energy system optimization, while for upscaling, Eskilstuna uses the building-stock model ECCABS, whereas Lower Austria uses a mixed integer linear programming optimization model, and the BALMOREL power system model. According to the modeling, HPs will dominate Eskilstuna’s heating sector by 2040. In Lower Austria, DH becomes more prevalent, in combination with wood biomass and HPs. These findings are explained by the postulated technological-economic parameters, energy prices, and CO2 prices. We conclude that future electricity prices will determine future heating systems: either a high share of centralized HPs (if electricity prices are low) or a high share of combined heat-and-power (if electricity prices are high). Large-scale energy storage and biomass can be essential solutions as may deliver increased cost-effectiveness, if available and under certain conditions.

A three-stage reliability-centered framework for critical feeder identification, failure modes prioritization, and optimal maintenance strategy assignment in power distribution system

In this paper, a three-stage reliability-centered maintenance (RCM) framework for critical feeder identification, failure modes prioritization, and optimal maintenance strategy assignment is presented for the power distribution system. This three-stage methodology addresses a notable gap in the literature by incorporating these three stages together. The first stage ranks critical feeders using the BWM and TOPSIS methods. In the second stage, the failure modes are prioritized using feeder rank, severity, occurrence, and detection factors and are assigned a score representing their priority. Finally, in the third stage, a maintenance strategy assignment, formulated as a mixed-integer linear programming (MILP) optimization problem is proposed. The optimization problem considers three strategies, cold-line, hot-line, and run-to-failure (RTF), assigned to each failure mode. Besides, the possibility of implementing cold-line repairs overlapping with the upstream grid's scheduled outages and the availability of several maintenance crews are taken into account. The objective function is to minimize the total maintenance costs, which include operation, equipment, energy not supplied (ENS), expected energy not supplied (EENS), and future corrective maintenance costs. The proposed approach is implemented on a real distribution system. The results indicate that the proposed RCM framework leads to a 26.32 % reduction in the EENS compared to selecting the RTF strategy for all components. Besides, the proposed optimization-based method leads to a 29 % lower objective function than the business-as-usual maintenance planning, and the proposed three-stage RCM framework outperforms a previously published state-of-the-art method by decreasing the objective function and the EENS.

The Optimal Size of a Heterogeneous Air Taxi Fleet in Advanced Air Mobility: A Traffic Demand and Flight Scheduling Approach

Introducing Advanced Air Mobility (AAM) as a novel transportation mode poses unique challenges due to limited practical and empirical data. One of these challenges involves accurately estimating future passenger demand and the required number of air taxis, given uncertainties in modal shift dynamics, induced traffic patterns, and long-term price elasticity. In our study, we use mobility data obtained from a Dresden traffic survey and modal shift rates to estimate the demand for AAM air taxi operations for this regional use case. We organize these operations into an air taxi rotation schedule using a Mixed Integer Linear Programming (MILP) optimization model and set a tolerance for slight deviations from the requested arrival times for higher productivity. The resulting schedule aids in determining the AAM fleet size while accounting for flight performance, energy consumption, and battery charging requirements tailored to three distinct types of air taxi fleets. According to our case study, the methodology produces feasible and high-quality air taxi flight rotations within an efficient computational time of 1.5 h. The approach provides extensive insights into air taxi utilization, charging durations at various locations, and assists in fleet planning that adapts to varying, potentially uncertain, traffic demands. Our findings reveal an average productivity of 12 trips per day per air taxi, covering distances from 13 to 99 km. These outcomes contribute to a sustainable, business-focused implementation of AAM while highlighting the interaction between operational parameters and overall system performance and contributing to vertiport capacity considerations.

A Mixed‐Integer Linear Programming Framework for Optimization of Water Network Operations Problems

AbstractWater distribution systems (WDSs) are critical infrastructure used to convey water from sources to consumers. The mathematical framework governing the distribution of flows and heads in extended period simulations of WDSs lends itself to application in a wide range of optimization problems. Applying the classical mixed integer linear programming (MILP) approach to model WDSs hydraulics within an optimization framework can contribute to higher solution accuracy with lower computational effort. However, adapting WDSs models to conform to a MILP formulation has proven challenging because of the intrinsic non‐linearity of system hydraulics and the complexity associated with modeling hydraulic devices that influence the state of the WDS. This paper introduces MILPNet, an adjustable framework for WDSs that can be used to build and solve an extensive array of MILP optimization problems. MILPNet includes constraints that represent the mass balance and energy conservation equations, hydraulic devices, control rules, and status checks. To conform to MILP structure, MILPNet employs piece‐wise linear approximation and integer programming. MILPNet was implemented and tested using Gurobi Python API. Modeling accuracy was shown to be comparable to EPANET, a public domain software for hydraulic modeling, and sensitivity analyses were conducted to examine the impacts of the modeling assumptions on the performance of MILPNet. Additionally, application of the framework was demonstrated using pump scheduling optimization examples in single and rolling horizon scenarios. Our results show that MILPNet can facilitate the construction and solution of optimization problems for a range of applications in WDSs operations.

Fertilizer Logistics in Brazil: Application of a Mixed-Integer Programming Mathematical Model for Optimal Mixer Locations

Background: Brazil is one of the largest consumers of fertilizers and is highly dependent on the international market to meet its demand for agricultural production inputs. The complexity of the fertilizer supply chain motivated us to carry out this study on redesigning the fertilizer logistics chain and evaluate strategies for reducing logistics costs by redesigning the fertilizer mixing network in Brazil, a country that is heavily dependent on imported fertilizers for agriculture. Methods: We introduce a multi-product mixed-integer linear programming optimization model encompassing the logistics network, from import ports to mixing factories and agricultural fertilizer supply centers. This model includes logistics infrastructure and taxes, accounting for greenhouse gas emissions (specifically carbon dioxide) in fertilizer logistics. Results: The results indicate that expanding the port capacity for fertilizer importation can significantly reduce logistics costs and greenhouse gas emissions by up to 22.5%, decreasing by 23.9% compared to the baseline. We also observed that removing taxes on fertilizer importation can reduce logistics costs by approximately 11%, but it increases greenhouse gas emissions by 2.25% due to increased reliance on road transport. We identified 15 highly resilient regions for establishing mixing factories, evaluated various scenarios and determined the importance of these locations in optimizing the fertilizer supply network in the country. Moreover, the results suggest a significant potential to enhance the role of Brazil’s Northern Arc region in fertilizer import flows. Conclusions: Public policies and private initiatives could be directed toward encouraging the establishment of mixing factories in the identified regions and increasing transport capacity in the Northern Arc region. Improving the logistical conditions of the fertilizer network would contribute to food security by reducing the costs of essential inputs in food production and promoting sustainability by reducing greenhouse gas emissions.