Interactive Simulation Model Research Articles

Efficient two-way fluid–structure interaction simulation for performance prediction of pressure-compensating emitter

Drip irrigation using a high-performance pressure-compensating (PC) emitter is one of the essential components for precision agriculture, and it is necessary to accurately predict its performance prior to design. In this study, an efficient two-way fluid–structure interaction (FSI) simulation model was developed and verified through an enlarged model experiment. The computational fluid dynamics (CFD) and computational solid mechanics (CSM) models of the FSI simulation were systematically verified, and a calibration method for the overestimated flow rate in the re-rising range was applied. The CFD model was determined to be the shear stress transport turbulence model, and the CSM model was determined to be the Ogden hyperelastic model for the PC emitter. The minimum prediction error for the flow rate was 7.93%, which was within 10% for all cases. The simulation model demonstrated its efficiency by analysing the performance of a single PC emitter with an average total analysis time of 18.6 h. In addition, by comparing various cases according to the design parameters, it is considered that the hardness of the diaphragm has a significant impact on the design of low-pressure PC emitters. The simulation model of this study can accurately predict the performance of PC emitter under specific conditions, yet improvement of simulation model is required to be applied in design optimisation. Future studies may benefit from combining an improved FSI simulation with a surrogate model to further enhance optimisation efforts.

Publisher's Note: “A modified weakly compressible smoothed particle hydrodynamics mixture model for accurate simulation of wave and porous structure interaction” [Phys. Fluids 36, 047128 (2024)

Publisher's Note: “A modified weakly compressible smoothed particle hydrodynamics mixture model for accurate simulation of wave and porous structure interaction” [Phys. Fluids <b>36</b>, 047128 (2024)

Use of virtual reality in physical rehabilitation: A narrative review

Virtual reality (VR) has emerged as an innovative technology in various fields. It transforms how we experience and interact with the world, and it has also developed into physical rehabilitation. Rehabilitation has also evolved from mere exercises to using technology while giving exercises to patients. VR is a technology that can enhance treatment and improve outcomes. VR uses computer-generated three-dimensional interactive simulation models, allowing users to engage in environments resembling real-world objects and events. An extensive literature search was conducted on SCOPUS, PubMed, the Web of Science and relevant rehabilitation databases. The keywords ‘virtual reality’, ‘rehabilitation’ and ‘physical therapy’ were employed to identify pertinent studies. The inclusion criteria encompassed research investigating the use of VR in rehabilitating musculoskeletal, neurological and cardiovascular conditions. VR has been used in various rehabilitation domains. It is effectively used in balance and gait training, cognition and memory training, pain management, muscle strengthening, etc. It has also been shown to improve adherence to exercises. VR training in physical therapy represents a transformative advancement in rehabilitation. Integrating VR technology offers many benefits, including enhanced engagement, customised therapy regimens and a safe environment for patients to practice and improve their physical abilities. The evidence from various studies indicates its potential to improve outcomes for individuals with a wide range of physical impairments.

A modified weakly compressible smoothed particle hydrodynamics mixture model for accurate simulation of wave and porous structure interaction

In this study, a modified weakly compressible smoothed particle hydrodynamics (WCSPH) mixture model was developed to more accurately simulate the interaction between waves and porous structures. In this model, we enhanced the governing equations of the traditional WCSPH mixture model by introducing Darcy velocity, apparent density, and an adjustable smoothing length. This refinement ensures that the modified model effectively maintains the conservation of fluid volume in seepage simulations. Additionally, this paper proposes a permeable interface treatment technique that replaces traditional smoothed particle hydrodynamics interpolation with finite element shape function interpolation, significantly enhancing computational efficiency. At the same time, we also introduced and revised a particle shifting technique, which further increases the computational precision of the model. The modified WCSPH mixture model was then applied to simulate several physical experiments, including the dam-break wave propagation in a permeable dam, the attenuation of solitary waves on a permeable riverbed, the propagation of the solitary wave on a submerged porous structure, and the breaking process of waves passing through permeable breakwaters. Through comparison with the experimental data and other numerical results, the current model was comprehensively verified from various aspects, such as fluid volume conservation, wave evolution in and around the porous structure, and pressure distribution characteristics. The results confirm the excellent performance of the current model in simulating the interaction between waves and porous structures.

An assessment of the impact of traditional rice cooking practice and eating habits on arsenic and iron transfer into the food chain of smallholders of Indo-Gangetic plain of South-Asia: Using AMMI and Monte-Carlo simulation model

The current study was designed to investigate the consequences of rice cooking and soaking of cooked rice (CR) with or without arsenic (As) contaminated water on As and Fe (iron) transfer to the human body along with associated health risk assessment using additive main-effects and multiplicative interaction (AMMI) and Monte Carlo Simulation model. In comparison to raw rice, As content in cooked rice (CR) and soaked cooked rice (SCR) enhanced significantly (at p < 0.05 level), regardless of rice cultivars and locations (at p < 0.05 level) due to the use of As-rich water for cooking and soaking purposes. Whereas As content in CR and SCR was reduced significantly due to the use of As-free water for cooking and soaking purposes. The use of As-free water (AFW) also enhanced the Fe content in CR. The overnight soaking of rice invariably enhanced the Fe content despite the use of As-contaminated water in SCR however, comparatively in lesser amount than As-free rice. In the studied area, due to consumption of As-rich CR and SCR children are more vulnerable to health hazards than adults. Consumption of SCR (prepared with AFW) could be an effective method to minimize As transmission and Fe enrichment among consumers.

Multi-objective optimisation of heat transfer and structural strength of aero-piston air-cooled engine cylinder based on orthogonal test

Owing to their compact and lightweight structure, air-cooled aviation piston engines are used in general aviation, agriculture, and military fields. To address the kind of engine thermal load problem, this study used a horizontal four-cylinder four-stroke aviation piston air-cooled gasoline engine as the research object. Tested the engine cylinder temperature and pressure under calibrated power conditions, built a one-dimensional simulation model and a fluid–structure interaction simulation model. Four main structural parameters—fin thickness, fin spacing, fin length, and cylinder wall thickness—were considered as factors. Each factors have five levels to construct orthogonal tests to analyse their influence on the cylinder heat transfer and structural strength performance. A range analysis of the orthogonal test showed that the fin thickness had the most significant influence on the thermal stress, while the fin length most significantly influenced the temperature and cylinder hole deformation. A comprehensive frequency analysis method combining intuitive and range analyses is used to determine the optimal combination scheme, and the heat transfer and strength performance of the engine cylinder before and after optimization is compared. The results indicated that the fin thickness and length had the most influence on cylinder performance. After optimisation, the temperature decreased by 3.3 % from 198.3 ℃ to 191.8 ℃, while the thermal stress decreased by 21.9 % from 100.5 MPa to 78.5 MPa. In addition, the cylinder hole deformation decreased by 11 %, from 0.273 mm to 0.243 mm. These research results can provide a reference for the design and optimisation of heat fins of air-cooled aviation piston engines.

Agent-based simulation model to evaluate government policies for farmers’ adoption and synergy in improving irrigation systems: A case study of Lake Urmia basin

The complexities of water issues, in which human factors directly affect the qualitative and quantitative indicators of water resources in a variable, extensive, and undefined manner, make social interactions necessary for the water resources sustainable development. In this situation, local government admins (local divisions of the Water Institutions and Ministry of Agriculture) (LGA) can play a fundamental role in attracting the participation of independent farmers in sustainable agriculture and preserving water resources by modifying support policies. In this study, an agent-based simulation model of social interactions and economic interests of farmers is presented along with different scenarios of support policies as well as random supervision and training by LGA agents to evaluate its impact on the independent decision making of farmers in the form of a complex adaptive system. Questionnaires and expert interviews were used to determine the data and scenarios. Farmers’ adaptation for development technology in irrigation management (DTIM), simulated output of WP, field application efficiency (FAE), and their adaptation speed are influenced by input factors such as assumed cost of sustainable water supply (SWS), mean government expenditure (subsidy count), and the level of supervision and training by LGA. The results showed that the assumed cost of SWS and government supervision and training are the most influential indicators of the farmer agent’s propensity to DTIM with F-value=4375.59 and F-value=1055.10, respectively. The results confirm that appropriate policies in strengthening LGA supervision and training and increasing farmers’ awareness of the importance of long-term water resource sustainability can help achieve more adaptability.

Investigation of the Dynamic Performance of High-Speed Railway Bridges in Cold Regions Based on Vehicle–Bridge Coupling Vibrations

Abstract This work uses numerical simulations to systematically study the dynamic performance of high-speed railroad bridges in cold areas based on coupled vibration analyses of axles and vehicles during freeze–thaw cycles. The freezing and thawing cycles of perennial frozen soil were established with comsol software, and the changes of the soil components during freezing and thawing cycles were analyzed with an indirect solution method and reasonable hydrothermal parameters and freezing and thawing boundary conditions. With the establishment of a pile–soil interaction simulation model, the change rule for pile foundation displacement and the force characteristics of pile foundations during freeze–thaw cycles were studied. abaqus and um software were used to simulate the coupled vibrations of axles, select the indexes for evaluation of train safety and the smoothness of the high-speed railroad trains, and investigate the influence of foundation freezing and pulling displacements on the dynamic performance of high-speed trains.

Establishment of the Interaction Simulation Model between Plug Seedlings and Soil

Currently, the simulation parameters for the model of the interaction between the transplanter, the plug seedlings, the soil, and the pot damage mechanism still need to be clarified. The optimization design of the planters and the improvement of planting quality are still urgent issues that need to be solved. In this paper, the simulation parameters of the pot and the soil were calibrated based on the pressure distribution measurement technology. With the actual collision impact force and matrix loss rate as the targets, a four-factor, three-level orthogonal test was designed to obtain the optimal parameters. Through the optimization analysis of the experimental results, it could be concluded that the pot–soil restitution coefficient, the pot–soil static friction coefficient, the pot–soil rolling friction coefficient, and the surface energy were 0.31, 0.88, 0.35, and 1.07 J/m2, respectively. The experimental verification of the optimal parameter combination showed that the relative error of the collision impact force was 1.65% and that the relative error of the matrix loss rate was 2.32%, which verified the model’s reliability. Based on the optimal parameters, the movement law of the hole tray seedlings was studied at different positions during the transplanting process. The plug seedlings collided not only with the planter but also with the soil, which led to the breakage and looseness of the pot structure. The relative error between the matrix loss rate of the transplanter inserting soil, the matrix loss rate of the transplanter that did not enter the soil, and the simulated matrix loss rate was less than 10%, which further proved the accuracy of the simulation model.

A numerical model for simulation of two-phase flows interaction with flexible slender bodies

This paper introduces a fluid–structure interaction (FSI) model for simulation of the coupled dynamics between two-phase flows and elastic slender structures. This model is extended from the early FSI model for single-phase flows by Wang et al. [“A coupled flow and beam model for fluid–slender body interaction,” J. Fluids Struct. 115, 103781 (2022)] to two-phase flows with interfaces of two liquids or free surfaces between gas and liquid. To better capture the interface movement of the two-phase flows and the interaction with the deformable structure, a consistent mass and momentum flux scheme is developed to reduce the spurious oscillation of fluid velocities near the interface, especially in the lower density (e.g., air) region. The proposed model is validated by a series of two-dimensional laboratory experiments of flow impact on a deformable plate, demonstrating that the model has good capabilities of conserving mass and momentum during the process of plate deformation by impulsive flow forces. The model is also applied to the investigation of three-dimensional dam break flow impact on a column of elastic plates. The complex interaction between the plates and the flow is discussed based on the simulation results.

Research on system of ultra-flat carrying robot based on improved PSO algorithm

Ultra-flat carrying robots (UCR) are used to carry soft targets for functional safety road tests of intelligent driving vehicles and should have superior control performance. For the sake of analyzing and upgrading the motion control performance of the ultra-flat carrying robot, this paper develops the mathematical model of its motion control system on the basis of the test data and the system identification method. Aiming at ameliorating the defects of the standard particle swarm optimization (PSO) algorithm, namely, low accuracy, being susceptible to being caught in a local optimum, and slow convergence when dealing with the parameter identification problems of complex systems, this paper proposes a refined PSO algorithm with inertia weight cosine adjustment and introduction of natural selection principle (IWCNS-PSO), and verifies the superiority of the algorithm by test functions. Based on the IWCNS-PSO algorithm, the identification of transfer functions in the motion control system of the ultra-flat carrying robot was completed. In comparison with the identification results of the standard PSO and linear decreasing inertia weight (LDIW)-PSO algorithms, it indicated that the IWCNS-PSO has the optimal performance, with the number of iterations it takes to reach convergence being only 95 and the fitness value being only 0.117. The interactive simulation model was constructed in MATLAB/Simulink, and the critical proportioning method and the IWCNS-PSO algorithm were employed respectively to complete the tuning and optimization of the Proportional-Integral (PI) controller parameters. The results of simulation indicated that the PI parameters optimized by the IWCNS-PSO algorithm reduce the adjustment time to 7.99 s and the overshoot to 13.41% of the system, and the system is significantly improved with regard to the control performance, which basically meets the performance requirements of speed, stability, and accuracy for the control system. In conclusion, the IWCNS-PSO algorithm presented in this paper represents an efficient system identification method, as well as a system optimization method.

Reliability analysis of 15MW horizontal axis wind turbine rotor blades using fluid-structure interaction simulation and adaptive kriging model

Over the course of the last four decades, the rotor diameter of Horizontal Axis Wind Turbines (HAWTs) has undergone a substantial increase, expanding from 15 m (30 kW) to an impressive 240 m (15MW), primarily aimed at enhancing their power generation capacity. This growth in blade swept area, however, gives rise to heightened loads, stresses and deflections, imposing more rigorous demands on the structural robustness of these components. To prevent sudden failure and to plan effective inspection, maintenance, and repair activities, it is vital to estimate the reliability of the rotor blades by considering all the forces (aerodynamic and structural dynamics) acting on them over the turbine's lifespan. This research proposes a comprehensive methodology that seamlessly combines fluid-structure interaction (FSI) simulation, Kriging model/algorithm and Adaptive Kriging Monte Carlo Simulation (AKMCS) to assess the reliability of the HAWT rotor blades. Firstly, high-fidelity FSI simulations are performed to investigate the dynamic response of the rotor blade under varying wind conditions. Recognizing the computationally intensive nature and time-consuming aspects of FSI simulations, a judicious approach involves harnessing an economical Kriging model as a surrogate. This surrogate model adeptly predicts blade deflection along its length, utilizing training and testing data derived from FSI simulations. Impressively, the Kriging model predicts blade deflection 400 times faster than the FSI simulations, showcasing its enhanced efficiency. The optimized surrogate model is then used to estimate the flap wise blade tip deflection for one million wind speed samples generated using Weibull distribution. Thereafter, to evaluate the reliability of the blades, statistical modeling using methods such as Monte Carlo Simulation (MCS), AKMCS is performed. The results demonstrate the faster convergence of AKMCS requiring only 21 samples, as opposed to 1 million samples for MCS with minimal reduction in the precision of the estimated probability of failure (Pf) and reliability index (β). Demonstrated on the backdrop of an IEA-15MW offshore reference WT rotor blade, the proposed methodology underscores its potential to be seamlessly incorporated into the creation of WT digital twins, due to its near real-time predictive capabilities for Pf and β assessments.

Simulation model of contact interaction during surface strengthening of steel parts

In the processes of surface strengthening of steel parts, the stress-strain state is decisive for explaining the physical processes of strengthening, forming the dimensions of the contact area. Analytical dependences of contact parameters are quite approximate. In this work, based on the Ansys software complex, a simulated model of the contact of a truncated torus with a cylinder is proposed, which demonstrates the kinetics of the process of pressing a hard alloy tool into a steel workpiece - a cylinder. The experiment was conducted for 4 seconds in order to determine the maximum level of stresses, the distribution of stresses and the amount of residual stresses after removing the load. The clamping force was applied mainly in the zone of elastic deformations. The results showed an uneven stress distribution with a maximum in the center of the contact spot of 1082 MPa. After changing the load direction, small residual deformations at the level of 0.00311 μm were observed in the center of the contact patch. This indicates a violation of the elastic region on a small contact area, which does not affect the general nature of the stress distribution and can be removed during the finishing process. The results of simulation of the stressed state are used for the correlation with the observed structural changes of the material during the action of thermal and power stresses. The stress peak was formed at a distance of 200 μm, which contributes to the formation of maximum values of microhardness at this depth.

A Study on the Adsorption Mechanism and Compactness of the TFS Coating Interfacial Layer

Chrome-plated plates, also known as tin-free plates (TFS), are the latest substrates for coating plates. The coating plate cannot be separated from the TFS during the stamping and extension process, and the interface layer of the TFS coating plate cannot produce pores to ensure good corrosion resistance and the appearance of the metal packaging cans. This requires the TFS coating plate interfacial layer to have good adsorption and compactness. In this paper, the molecular simulation model of the interfacial layer interaction of the TFS coating plate was established by using molecular mechanics simulation, Monte Carlo simulation, and molecular dynamics simulation, and the influential rules of chromium oxide crystalline structure, coating functional group type, and coating pressure on the adsorption and compactness of interfacial layer were analyzed and verified by experiments. The results show that the adsorption is stronger when the surface of the TFS is a chromium oxide (110) crystalline surface and contains hydroxide ions. The adsorption of polyester polyurethane coating and polyether polyurethane coating for and the adsorption of polyester polyurethane coating functional groups is stronger than polyether functional groups, and the adsorption of other functional groups is ranked by the same method. The interfacial layer compactness increases with an increase in coating pressure. For this experimental sample, the value of the film pressure sensor is 18,940 g when meeting the requirements of adsorption and compactness of the interfacial layer of the TFS coating plate, which can be extended for other coating plates.

VOICE INFORMATION IN JSC “RUSSIAN RAILWAYS”

Background. The work is devoted to the study of the processes of interaction between the company and consumers within the framework of the information services provided in order to improve or innovate the methods of voice communication for further operation and automation; as well as to substantiate the need for local research through the growth of automation of interaction and the speed of obtaining and perception of information.&#x0D; Materials and methods. Materials from open sources on the Internet, simulation modeling, data analysis, systems theory, system approach.&#x0D; Results. Information systems at Russian Railways using voice information methods have been studied, the area within which automation of user interaction with the system is possible has been analyzed, a simulation model of customer interaction with the system using voice information has been built to ensure the issuance of information on request.&#x0D; Conclusion. Based on the results of the study, a solution was proposed for automation in the field of event-based informing of customers using the transport services of Russian Railways, company employees in order to improve information processing and automation of work with it. The proposed model will demonstrate and justify the issuance, processing and transmission of information for the transition to the introduction of voice information, through an assessment of the perception and reaction of the company’s employees to the implementation of audio accompaniment and interaction with the system through voice. As a result, it is proposed to develop a model reflecting the possibility of working with an information system using voice guidance and information.

INTERDISCIPLINARY ANALYSIS OF THE AGE CHARACTERISTICS OF TEENAGERS IN THE CONTEXT OF INTERCULTURAL COMMUNICATION

The article substantiates the relevance of the problem of forming teenagers' readiness for intercultural communication in the age context from the point of view of interdisciplinary analysis. The reasons for the relevance of studying such issues are indicated: communication with representatives of other cultures, implementation of academic, professional or personal mobility. The purpose of the presented research is to do the interdisciplinary analysis of the age characteristics of teenagers` in the context of intercultural communication from pedagogical, social, psychological, physiological and axiological aspects. Analysis and systematization of scientific literature, classification, generalization were chosen as the leading methods. Positions are presented that are fundamental in the context of the research problem, namely: 1) the process of their socialization has a strong influence on teenagers, primarily in their own environment; 2) a person's life motivation begins to form at the teenager`s age; 3) taking into account the capabilities of the teenager`s body and giving feasible tasks. It was found that in order to engage all students at the lesson, it is desirable to use interactive collective learning methods that allow the teacher to establish emotional and intellectual relationships with teenagers aimed at achieving joint success. The variability of modern methods and technologies, which involve the use of technological, software and technical means, is revealed, namely: training methods, interactive collective learning methods, web quests, project methods, role-playing games, dialogue technologies, simulation and interactive modeling methods. It is noted that the indicated methods and technologies of forming teenagers' readiness for intercultural communication ensure awareness of their own cultural identity, interest and respect for universal and national cultural values.

Port-of-Entry Simulation Model for Potential Wait Time Reduction and Air Quality Improvement: A Case Study at the Gateway International Bridge in Brownsville, Texas, USA

The mathematical study known as queueing theory has recently become a major point of interest for many government agencies and private companies for increasing efficiency. One such application is vehicle queueing at an international port-of-entry (POE). When queueing, fumes from idling vehicles negatively affect the overall health and well-being of the community, especially the U.S. Customs and Border Protection (CBP) agents that work at the POEs. As such, there is a need to analyze and optimize the border crossing queuing operations to minimize wait times and number of vehicles in the queue and, thus, reduce the vehicle emissions. For this research, the U.S.–Mexico POE located at The Gateway International Bridge in Brownsville, Texas, is used as a case study. Due to data privacy concerns, the hourly wait times for vehicles arriving at the border had to be extracted manually each day using a live wait time tracker online. The data extraction was performed for the month of March 2022. Using these wait times, a queueing simulation software, SIMIO, was used to develop an interactive simulation model and calibrate the service rates. The output from the SIMIO model was then used to develop an artificial neural network (ANN) to predict hourly particulate matter content with an R2 of 0.402. From the ANN, a predictive equation has been developed, which may be used by CBP to make operational decisions and improve the overall efficiency of this POE. Thus, lowering the average wait times and the emissions from idling vehicles in the queue.

Mitigation of the blast load on RC column by hanging granite slabs

Reinforced concrete (RC) column, which is the critical axial load carrying structural component, is vulnerable to blast loads, especially to close-in blast loads. The blast damage of the RC column might lead to the progressive collapse of the whole building structure. Therefore, it is necessary to ensure that the RC column could survive under design blast loads once the building structure is at the risk of terrorist or accidental explosions. One of the solutions is that the RC column is designed to be strong enough to resist the blast loads, which is very expensive. While the other method is to protect the RC column against blast loads with energy absorb materials. The granite slab, which is often used as the decoration hanging outside the structural components, is utilized to mitigate the blast loads applied on the RC column. In this study, firstly the concept of hanging granite slab for mitigating blast load on RC columns is proposed. Furthermore, the finite element (FE) model for simulation of the blast wave interaction with the dry-hanging granite slab is established to reveal the protection mechanism of RC column hanging with granite slab based on the platform LS-DYNA. The accuracy of FE model is validated by comparing the numerical results and the test results from blast field test. The influences of slab thickness T and hanging distance D on blast mitigation are investigated. Finally, the blast load mitigation factors are proposed and their empirical model are established as a function of scaled distance Z, hanging slab thickness T and hanging standoff distance D.

Design of emotional branding communication model based on system dynamics in social media environment and its influence on new product sales.

In the current social media environment, emotional branding communication has become a common marketing tool for brand owners, and therefore it has become particularly important and urgent to study it. Based on the perspective of brand equity theory, combined with the new characteristics of marketing communication in the social media environment, this paper constructed an emotional branding communication model in the social media environment. The system dynamics (SD) method was used to simulate and analyze the new product marketing system to assess whether it could stir the emotional needs of the consumers and resonate within their hearts. This paper discusses the asymmetric communication of different brands regarding the same commodity to determine the impact of this exchange mechanism, that is, only the weak brands in the market initially adopt marketing methods, while the strong brands do not participate in social marketing activities. It was found that the influence of marketing frequency and marketing intensity on symmetric and asymmetric communications was different. In the face of different types of competitors, the marketing strategy of weak brands needs emphasis. Through unit consistency test, structure verification test, effectiveness, and rationality test, it was proven that the emotional branding communication model and new product sales interaction simulation model established in this paper were reasonable and effective.

Developing intelligence in food security: An agent-based modelling approach of Qatar's food system interactions under socio-economic and environmental considerations

Maintaining the consistent provision of food is becoming a persistent challenge in the light of increasing demands for resources and the multiple risks governing food systems. These include population growth, market volatilities, changing climate and the lack of coordination amongst involved stakeholders. Dynamic and interactive simulation models have demonstrated their utility for developing possible scenarios related to the food sector's subsystems, allowing stakeholders to understand their outcomes, should they occur, and make efficient and sustainable decisions that potentially alleviate the risks. This study proposes an Agent-based model (ABM) that supports informed decision-making through allowing the simulation of multiple scenarios that can potentially influence the food system's strategies and actions. The framework is illustrated considering the economic actors within Qatar's food sector, with a special focus on delivering products in high demand, such as the tomato crop. The model was implemented using the Python-based Mesa library, and it simulates interactions between domestic farmers, local importers, and international exporters to assess the performance of the food sector considering economic, environmental and social scenarios. Scenario (1) describes the baseline case wherein local farmers and importers are acting independently with the aim to maximise their individual profits. Scenario's (2) and (3) represent enhancement plans that foster collaboration, water security and improved food quality. Findings of this study assert that from an economic standpoint, the food system can supply 75% of the market with imported tomatoes, as for the remaining 25%, it can be provided through national farms, with a 76% contribution of greenhouses and 24% open fields. This option offers the minimum economic cost and the lowest Greenhouse Gas (GHG) emissions. However, this perspective implies massive food losses and relatively high-water consumption. Scenario 2, which advocates a full reliance on imports offers a costly monthly strategy in terms of economic expenses and environmental emissions, however, it results in considerable water savings that relieve the strain on the local water supply system. When implementing a plan driven by quality consideration (scenario 3), the reliance on imports reaches up to 82%, and the remaining local produce is mostly supplied by greenhouses, with approximately 78%, leading to quality improvement and a substantial decrease in food losses yet at the expense of both the economic and environmental performance. These results confirm that aspiring water security and social welfare while provisioning food products requires tremendous investment and might lead to significant GHG emissions. The framework developed in this study provides a generic and interactive decision-making scheme that can assist in undertaking future decisions considering economic and environmental constraints and market-driven specifications.