IoT Task Offloading in Edge Computing Using Non-Cooperative Game Theory for Healthcare Systems

2024-03-02 01:31DineshMavaluruChettupallyAnilCarieAhmedAlutaibiSatishAnamalamudiBayapaReddyNarapureddyMuraliKrishnaEnduriandMdEzazAhmed

Dinesh Mavaluru,Chettupally Anil Carie,Ahmed I.Alutaibi,Satish Anamalamudi,Bayapa Reddy Narapureddy,Murali Krishna Enduri and Md Ezaz Ahmed

1School of Computing and Informatics,Saudi Electronic University,Riyadh,Saudi Arabia

2Department of Computer Engineering,SRM University AP,Guntur,India

3College of Computer and Information Sciences,Majmaah University,Al Majma’ah,Saudi Arabia

4Department of Public Health,College of Applied Medical Sciences,King Khalid University,Abha,Saudi Arabia

ABSTRACT In this paper, we present a comprehensive system model for Industrial Internet of Things (IIoT) networks empowered by Non-Orthogonal Multiple Access (NOMA) and Mobile Edge Computing (MEC) technologies.The network comprises essential components such as base stations, edge servers, and numerous IIoT devices characterized by limited energy and computing capacities.The central challenge addressed is the optimization of resource allocation and task distribution while adhering to stringent queueing delay constraints and minimizing overall energy consumption.The system operates in discrete time slots and employs a quasi-static approach,with a specific focus on the complexities of task partitioning and the management of constrained resources within the IIoT context.This study makes valuable contributions to the field by enhancing the understanding of resourceefficient management and task allocation,particularly relevant in real-time industrial applications.Experimental results indicate that our proposed algorithm significantly outperforms existing approaches,reducing queue backlog by 45.32%and 17.25%compared to SMRA and ACRA while achieving a 27.31%and 74.12%improvement in QnO.Moreover,the algorithm effectively balances complexity and network performance,as demonstrated when reducing the number of devices in each group(Ng)from 200 to 50,resulting in a 97.21%reduction in complexity with only a 7.35%increase in energy consumption.This research offers a practical solution for optimizing IIoT networks in real-time industrial settings.

KEYWORDS Internet of Things;edge computing;offloading;NOMA

1 Introduction

In recent years, the convergence of Internet of Things (IoT) and edge computing has brought transformative changes across industries, notably in the healthcare sector.The widespread use of connected medical devices and the growing demand for real-time data analysis have underscored the need for innovative solutions capable of efficiently managing and processing the vast volumes of data generated by these devices [1,2].A key strategy in addressing this challenge lies in leveraging edge computing,where data is processed in proximity to its source,reducing latency and enhancing overall system performance.

With advancements in wireless communication, the contemporary distributed network offers expanded coverage,higher capacity,and improved connectivity,facilitating the realization of IoT.IoT devices are increasingly tasked with computationally intensive and time-sensitive operations,especially with the development of sophisticated load control and real-time computing.This has given rise to the need for efficient resource management solutions, given the constraints of limited computation and storage resources on the device side[3].While IoT devices can offload tasks to the nearest edge servers to mitigate delays and conserve energy,the burgeoning number of terminal devices in IoT has made computational tasks more challenging.

As the demand for communication and spectrum resources continues to surge,the development of resource allocation and management algorithms has become urgent.Non-Orthogonal Multiple Access(NOMA) technology presents a practical solution for efficiently using scarce spectrum resources[4,5].However, several challenging problems persist.Jointly optimizing resource unit allocation,computing resource offloading,and task division algorithms is computationally complex due to the coupling between these entities.Moreover, striking a balance between long-term optimization goals and real-time decision-making amplifies the computational complexity of the model.The intricacies of multidimensional task allocation without future data,such as task arrival times and channel state information, add to the challenge.Additionally, the dynamic interference among devices when the same resource units are allocated to different devices creates externalities, leading to exponentially complex resource unit allocation problems.

In response to these challenges,the research explores a multi-time scale multidimensional resource allocation approach for NOMA-based Mobile Edge Computing(NOMA-MEC)in IoT.The proposed algorithm aims to minimize long-term device consumption by jointly optimizing resource unit allocation and task splitting.This approach dissects the long-term stochastic joint optimization problem into three short-term deterministic problems using the Lyapunov technique.It focuses on wireless spectrum allocation,task splitting,and task computation.To enhance efficiency and reduce complexity, IoT devices and resource units are grouped based on clustering schemes [6,7].Wireless spectrum allocation is treated as a one-to-many matching problem,resolved on a larger time scale and executed on the base station(BS)side.The suggestion is to use group-switching matching to allocate resource units.Devices and resource units are initially categorized into several groups,with switching matching performed within each group.Task splitting and computing resource allocation are then handled on the device side in a distributed manner on a smaller time scale[8].

The research presents a novel algorithm for resource unit allocation, which optimizes energy consumption while taking into account the dynamic nature of IoT device preferences and the scarcity of spectrum resources[9].

-Task Splitting Efficiency:The study introduces a task splitting approach that improves the overall system performance by jointly optimizing work partitioning and computing resource allocation,reducing queue backlogs.

- Multi-Time Scale Solution: By addressing long-term queuing delay constraints, the research offers a multi-time scale, multidimensional resource allocation strategy for NOMA-MEC in IoT,striking a balance between complexity and network performance.

- Group-Based Approach: The proposed approach utilizes clustering schemes to group IoT terminals and resource units,reducing computational complexity and enhancing efficiency in wireless spectrum allocation.

2 Related Works

With the advent of Industry 4.0,Industrial IoT(IIoT)has become a pivotal technology,addressing various challenges in smart factories, including the management of industrial big data, equipment monitoring, and maintenance [10,11].IIoT devices generate a substantial amount of data, and some industries require real-time data processing.Conventionally, industrial data is transferred to centralized cloud servers for processing due to their powerful computing capabilities [12].However,the sheer number of IIoT devices can lead to network congestion,causing delays in data transmission and impacting real-time tasks.Therefore,a network infrastructure capable of real-time and efficient data processing is essential for IIoT.

Edge computing is emerging as a valuable complement to cloud computing,utilizing a distributed computing approach to alleviate network congestion and reduce transmission delays associated with cloud computing’s centralized processing[13].By distributing computing tasks across multiple servers throughout the network,edge computing ensures a quicker response to user requests and tasks[14].To minimize data transmission delays in long-distance communication,edge servers are strategically placed closer to the devices.

Several ongoing studies explore the application of edge computing in IIoT environments.These studies optimize time and energy costs to achieve the efficient distribution of green and energy-saving computing resources[15,16].They addressed offloading challenges for multi-hop computing tasks in hybrid edge cloud computing environments,utilizing game theory to meet service quality requirements[17].In [18], authors introduced an innovative approach that integrates Wireless Power Transfer(WPT) with Mobile Edge Computing (MEC) for efficient computation offloading, achieving realtime performance in large-scale networks and addressing energy and latency constraints in wireless environments.In[19], authors optimized task dependencies in IoT edge computing using a directed cyclic graph model and priority-aware scheduling,outperforming other offloading methods in terms of throughput and task satisfaction rate.In[20],authors addressed the growing need for data analysis in the context of IoT systems, where the volume of generated data poses performance challenges.It introduces the IoT-SCOM model, focusing on minimizing transmission latency in edge-cloudhybrid systems.Experimental results show that IoT-SCOM offers improved accuracy and efficiency compared to existing methods, enhancing data-intensive service element deployment in the edgecloud environment.Moreover,a hybrid computing architecture with intelligent resource planning is proposed to meet real-time requirements[21].

In summary, edge computing offers solutions to meet quality of service demands and reduce system overhead by deploying at the network edge,closer to field devices,and providing suitable computing resources.This research investigates computing resource allocation and task offloading in IIoT environments using edge computing and NOMA.It aims to resolve issues related to communication delay, energy consumption, and system scalability by efficiently allocating computer resources and offloading responsibilities from edge servers.The contributions of this paper are twofold.

1.Decomposition based on multi-dimensional optimization: Utilizing Lyapunov optimization,the complex,long-term stochastic multi-dimensional optimization problem is decomposed into three deterministic problems that can be solved effectively.

2.Resource Allocation Problem based on group Switching Matching:The proposed resource unit allocation algorithm,based on group switched matching,provides a more efficient,less complex,and flexible approach to handle interdependencies among IIoT devices and various resources.

The remainder of this paper is structured as follows.Our system model is introduced in the second section.The introduction of the problem model,transformation,and queuing delay limitations follows in Section 3.The proposed algorithm is further explained in Section 4.The simulation results and data analysis are presented in the fifth section.The paper is finally summarized in the sixth section.

3 System Model

This paper examines a common NOMA MEC-based IIoT example,which includes a base station,an edge server, and a significant number of IIoT devices with low energy and computing capacities.The BS is place d where the edge server is,as stated in Fig.1,and it offers radio access and computing services forNIIoT devices.The node symbolizes the attack state,and its set is N,with the set being defined asN= {1,2,...,n,...,N}.The IIoT device either partially processes its tasks or offloads them to the BS.Comparable to[22],we investigate a discrete timeslot form that splits the optimization process intoTtime slots, each with a durationτ, as displayed in Fig.2.The timeslot collection is indicated by the expressionT= {1,...,t,...,T}.Consider a quasi-static example where the CSI varies between time slots but stays the same in one time slots.Each subsequentT′timeslot is connected to a time epoch denoted by the symbolss∈S,S={1,2,...,S}.The definition of the s-th time epoch isT (s)={(s-1)T′+1,...,sT′,sT′+1,...,(s+1)T′}.

The resource is divided intoMtime-frequency resources units with bandwidthBand a time period ofT′τ,the union of which isM={1,...,m,...,M}.The multi-dimensional resource allocation and splitting problems are studied.At the start of each time interval,the infinite resource allocation scheme is optimized over a major timescale.The resource allocation scheme is represented through the binary indicatory(s)={,n∈N,m∈M},where= 1 dicates that resource unitmis allocated to the devicenin the s-th time interval,otherwise=0.Next,computing resource allocation and splitting are jointly optimized on a small-time scale based on resource unit allocation strategy in each timeslot.

Figure 1:System model

Figure 2:Discrete timeslot model

A.Traffic Model on the Device

The notations used in the mathematical analysis is represented in Table 1 whereas description of each equation is mentioned in Table 2.The task partitioning system is used in this paper[23],and it allows for the division of each task into separate sub-tasks of sizeX0(bit).Assume thatsubtasks reach the devicenat thet-th timeslot.The arriving job can be split intoan(t)sub-tasks for local operation and(t)sub-tasks for task offloading,which are two different and parallel models.Given the foregoing,the task splitting at devicenin timeslot t is defined as

Table 1: Notations and equations in the system model

Table 2: Equations with meanings/descriptions

whichis the task length of local processing of devicenat timet,andrepresents the task length when the device n performs computing offloading at timet.

whichandrespectively indicate the amounts of data departureand.

B.Locally Task Data Processing Scheme

The quantity of locally executed data for thet-th timeslot is defined as

wherefn(t)denotes the number of CPU cycles that were allocated to deviceNduring thet-th timeslot,this is based on the 3GPP channel model.The required number of CPU cycles per bit,or computation intensity,isln.At thet-th timeslots,the local processing computing delay and the corresponding energy consumption generated by devicenare defined as

where the chip structure regulates a constant power coefficient calledκn.

C.Based on Task Offloading Scheme

Typically,BS uses successive interference cancellation(SIC)in NOMA-based systems to lessen the acquired signal’s interference.All other signals are viewed as interference,but the BS decodes signals from equipment with higher channel gains in the right order.The channel gain of the uplink from device n to resource unit m in thet-th timeslot is written ashmn (t), and the signal noise ratio (SNR)received at BS is given as follows:

The transmission power isp.Intracellular interference caused by other low channel gain devices makes up the first part of the denominator,while Additive White Gaussian Noise Power makes up the second.The volume of data from a task that can be offloaded in thet-th timeslot and the transmission speed of device n using resource unit m are therefore determined below:

The corresponding energy consumption produced by device n in thet-th timeslot is obtained as follows:

4 Problem Definition and Formulation

In this section,the queueing delay constraint is introduced first.Next,the optimization problem of multi-dimensional resource allocation and task splitting is presented.Later on, the proposed Algorithm is addressed.

A.Queued Delay Constraint

The queueing delay constraint is controlled to assure the effectiveness and timeliness of the offloading task.Little’s Law [24] states that the queuing delay ofandis calculated as

B.Problem Description

By collectively optimizing resource unit allocation,splitting,and computing task scheduling,in the case of queue delay restrictions, the goal is to decrease the overall cumulative long-term energy consumption of all devices.It is possible to express the multidimensional resource allocation and split optimization strategy as

wherey= {y(s)},s∈Srepresents the resource unit allocation vector,a= {a(t)},t∈Trepresents the radio resource splitting vector, anda(t)= {an(t),n∈N},f= {f (t)=fn(t),n∈N,t∈T}represents the computation resource allocation vector.C1is rhe resource splitting constraint.The computational resource allocation restriction on the device is indicated byC2.C3-C5denote that each device can use a maximum of one resource unit.The resource unit m can be allocated to devices up to Nm.According toC6, the resource unit allocated to devicenmust make sure that the SNR obtained at BS is more than the required minimumSNRn.

C.Transformation of the Problem

P1 is an non-deterministic polynomial (NP)-hard issue that is challenging to directly solve.By using Lyapunov optimization [25], the original long term stochastic optimization issue is split up into several deterministic subproblems in the short term.Formulas(11)and(12)can be transformed into the queue stability constraint problem for the virtual queue concept.It is possible to convert the matching virtual queuesandinto

Whenandare average rate stable, the formulas(11) and (12) will be maintained automatically.So P1 can be converted to

It can be concluded that P2 will be decomposed into three optimization sub-tasks,namely SP1:resource unit allocation sub-problem,SP2:task splitting sub-problem,and SP3:computing resource allocation sub-problem.

5 Resource Allocation and Task Splitting

In this section,we first introduce the three decomposed sub-schemes and their related responses and then conclude the raised scheme.

A.Resource Unit Allocation Optimization

At the start of each interval in SP1,devicenand BS choose the resource unit allocation strategy.As,andvary dynamically in timeslot,their empirical averages,namely,andare considered.So SP1 is going to be

A 1-to-n matching method between devices and resource unit is the default setting for SP1.The 1-to-n matching function eta that we are modeling here performs as follows:

Performances 1 and 2 dealt with the constraintsC4andC5,while performance 3 suggests that if devicemand resource unitmmatch,then resource unit m is assigned to devicen,and the opposite is t rue.In particular,there is a strong link between the resource unit allocation indexandλ.

The following are the utility functions for devicenand resource unitm:

It is not practical to match all devices and resource unit due to the extensive IIoT networks’high matching complexity.So, we start by dividing devices and resource unit into sets.Especially,based on the clustering strategy [26], devices and resource unit are partitioned intoKsets, i.e.,N= {NG1,...,NGk,...,NGK}andM= {MG1,...,MGk,...,MGK},Each union hasNg=N/kdevices andMg=M/kresource unit,respectively.The devices in every union usually have regulated resource unit,i.e.,MGKallocated to theNGKfor offloading.Then,switching matching is carried out in each set in a semi-distributed mode,allowing resource unit and devices to build their priorities in descending order based on their utilities.Therefore,two-end group switching matching is used to settle the resource unit assignment problem.

Definition 1.Define matchesλand two device resource unit pairs(n,m),(v,l)∈λ, namely,λ(m)=n,andλ(v)=l,∀n/=vandn,v∈NGk, ∀m/=l,m,l∈MGk,if they meet

={λ(n,m),(v,l)}∪{(n,l),(v,m)}is depicted as the variablesλand≻λswapping process.

Definition 2.when there is no exchange match, the matchedλis bilateral exchange stable.The proposed algorithm summarizes the details of resource unit allocation based on group switching matching.During initialization,devices and resource unit are divided intoKunions,and each device union is assigned resource unit set,namelyNGk←MGk.Subsequently,if only all the restrictions in formula(15)are met,the devices and resource unit in the allocation group will randomly match each other.Each device and resource unit produces a preference based on the formula(20).

In the switching matching period,each deviceninNGKcurrently matching the resource unitminMGKtakes a program to its most desired resource unitlinMGK.For each devicevinNGKthat existing matchlinMGK,If≻λand satisfies all the requirements in formula(15),the old matchλis replaced with the new match.Otherwise,λstays the same.The stage ends until no matches are exchanged.

At the end stage,the formula(19)transforms the finalλinto the resource unit allocation indexy∗(s).

B.Task Splitting and Resource Allocation Optimization

Task splitting issue with the following formula, in the t-th timeslot, SP2 distributes the task splitting decision between local operation and offloading.

Computing resource allocation sub-scheme SP3 controls the number of CPU cycle frequencies each device allocates for local operation in thet-th timeslot,which is gained via the following formula:

Both SP2 and SP3 are convex optimization problems,and Lagrange duality decomposition is easy to resolve.Reference can be made to similar specific derivation process[27].

C.Our Proposed Algorithm

The proposed algorithm consists of three steps,namely

Step 1:Initialize all queue backlog and resource unit allocation strategy metrics.

Step 2:The best resource unit allocationy∗(s)is achieved by each device in a partially distributed fashion according to Algorithm 1 and sends data using the allocated resource unit.

Step 3:Each device learns the ideal resource allocation and task distribution approach.The efficiency of data transmission, power usage, queue overhang, queue delay, and updatesandaccording to the formulas(2), (3), (13),and(14),are then taken into account by each device.Up untilt>T,the iteration between the second and third steps continues.

6 Experiments

This part uses simulation to assess the suggested algorithm.In this study, a unitary cell with a 2000 m radius is considered.Each 100 group contains 20 devices,dividing all equipment and resource unit in the same way.Table 3 lists the exact simulation parameters.

Table 3: Simulation parameters

ACRA(Access Control and Resource Allocation Algorithm):ACRA is an advanced algorithm designed for resource allocation and access control in the context of edge computing and IoT(Internet of Things).It addresses the challenge of efficiently managing resources and allocating them to IoT devices.The key features of ACRA include:

Resource Allocation:ACRA optimizes the allocation of computing and communication resources to IoT devices, ensuring efficient utilization of available resources.Access Control: It implements an access control mechanism, allowing IoT devices to submit tasks to the nearest edge server for processing.This minimizes task delays and conserves energy.Real-time Decision-Making: ACRA balances the need for long-term optimization goals with the requirement for real-time decisionmaking.This is crucial for ensuring that tasks are processed efficiently and in a timely manner.Multi-Time Scale Optimization:ACRA decomposes the complex,long-term stochastic optimization problem into short-term deterministic sub-problems.This approach simplifies resource allocation and improves system performance.Clustering Schemes:The algorithm employs clustering schemes to group IoT terminals and resource units,reducing computational complexity and enhancing efficiency in wireless spectrum allocation.

SMRA(Switch Matching Resource Allocation Algorithm):SMRA is another resource allocation algorithm that focuses on optimizing energy consumption within IoT networks.It is designed to allocate resources efficiently, particularly in terms of energy usage.Key characteristics of SMRA include: Energy Optimization: SMRA primarily aims to minimize energy consumption in IoT networks.It is especially well-suited for scenarios where energy efficiency is a critical consideration.Task Processing:It ensures that devices submit tasks to edge servers when they cannot be completed within the required timeframe.This reduces task delays and saves energy.Subtask Offloading:SMRA emphasizes offloading subtasks from devices to edge servers to reduce queue backlogs and enhance overall system performance.

Comparison and Evaluation: To provide a clear understanding of their respective roles and effectiveness,ACRA and SMRA should be compared based on various evaluation parameters,such as:Energy Efficiency:A comparison should be made to determine which algorithm is more effective in optimizing energy consumption.Task Processing Speed:Evaluation should assess which algorithm minimizes task delays and enhances real-time performance.Queue Backlogs: The impact of ACRA and SMRA on queue backlogs and system performance should be compared.Resource Utilization:An analysis of how efficiently these algorithms allocate and utilize computing and communication resources.Complexity: Assess the computational complexity of each algorithm, which can affect practical implementation.

We compared the two most advanced algorithms.The first is the resource unit allocation algorithm based on switch matching (SMRA) proposed in [28].It should be noted that in the simulation,minimization of energy consumption replaces SMRA maximum of energy efficiency.The second is due to the Lyapunov optimization and pricing matching based Access control and resource allocation algorithm(ACRA)developed.The task split phase is determined randomly,and the local computing resource are set to maximum values in the SMRA and ACRA.

Figs.3–5 show the average power consumption and the average backlog ofandwith timeslots.SMRA outperforms other algorithms in terms of energy consumption because only energy optimization is considered,but at the expense of much inferior backlog performance.Our suggested technique reduces QLM by 45.32%and 17.25%percentage compared to SMRA and ACRA because it jointly optimizes work partitioning and computing resource allocation.In addition, due to the consideration of externalities, our proposed algorithm betters than SMRA and ACRA by 27.31%and 74.12%on

Figure 3:Average energy consumption

Figure 4:Average backlog of

As a result,queue backlogs can be decreased by offloading more subtasks from the device to the edge server.

Fig.6 illustrates how the number of devices in each group ofNgaffects the complexity and average power consumption of our suggested approach.Define the whole number of devices to N=200.Complexity is set as the number of iterations of the exchange demanded to fulfill the match between all resource unit and devices.WhenNgis reduced from 200 to 50, the complexity of our proposed algorithm is decreased by 97.21%, whereas the energy consumption is just increased by 7.35%.We can draw the conclusion that our suggested algorithm can effectively balance complexity and network performance.

Figure 5:Average backlog of

Figure 6:Impact of Ng

7 Discussion

The proposed algorithm is superior to the existing algorithms, SMRA and ACRA, based on several critical metrics:

Energy Consumption: The primary metric indicating the superiority of the proposed algorithm is energy consumption.The results show that the proposed algorithm effectively optimizes energy consumption compared to SMRA.This is a crucial metric for IoT systems where energy efficiency is paramount.The exact percentage reduction in energy consumption should be provided to quantify this improvement.

Queue Backlog: Another important metric is the queue backlog (QLM), specifically for QnLand QnO.The proposed algorithm significantly reduces QLM,as indicated by a 45.32%and 17.25%reduction compared to SMRA and ACRA.A lower queue backlog is essential for real-time and delaysensitive applications,ensuring that tasks are processed promptly.

Complexity: The complexity of the algorithm is also a vital consideration.As demonstrated in the experiments, when the number of devices (Ng) is reduced from 200 to 50, the complexity of the proposed algorithm decreases by 97.21%.Lower complexity is desirable for practical implementation,as it reduces the computational burden and improves system efficiency.

Externalities: The proposed algorithm outperforms both SMRA and ACRA by 27.31% and 74.12% in QnOdue to the consideration of externalities.Externalities are crucial in a real-world IoT environment where devices may interfere with each other.Minimizing externalities is essential for efficient resource allocation and network performance.

8 Conclusion

In conclusion, this study has introduced a comprehensive framework for optimizing resource allocation and task offloading in Industrial IoT (IIoT) environments, leveraging edge computing and Non-Orthogonal Multiple Access(NOMA).Through the use of directed cyclic graphs,priorityaware scheduling,and Markov decision processes,the proposed approach addresses the challenges of network congestion,delay sensitivity,and efficient resource management in IIoT.The results indicate that our strategy outperforms various offloading methods in terms of both throughput and task satisfaction rate.By considering deadlines, user mobility, and dynamic task dependencies, we have shown the potential for real-time,efficient,and responsive IIoT systems.

For future directions,research can explore:

1.Dynamic Network Topologies:Investigating how the proposed framework performs in scenarios with changing network topologies,which is common in IIoT environments.

2.Security and Privacy:Expanding the framework to incorporate enhanced security and privacy measures to protect sensitive IIoT data.

3.Energy-Efficient Edge Servers:Exploring methods to optimize the energy consumption of edge servers,further enhancing the sustainability of IIoT systems.

4.Scalability:Extending the framework’s scalability to accommodate a growing number of IIoT devices and applications.

5.Hybrid Approaches: Combining edge computing with other emerging technologies, such as blockchain or quantum computing,to explore synergies for IIoT applications.

Acknowledgement:Authors would like to thank the anonymous reviewers for their valuable time and feedback in reviewing the manuscript.

Funding Statement:The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through large group research project under Grant Number RGP2/474/44.

Author Contributions:The authors confirm contribution to the paper as follows:Problem statement and proposed idea:Dinesh Mavaluru,Anil Carie Chettupally,Satish Anamalamudi;Data collection:Ahmed I.Alutaibi,Bayapa Reddy Narapureddy;Mathematical analysis and interpretation of results:Murali Krishna Enduri,Satish Anamalamudi,Ahmed I.Alutaibi;Draft manuscript preparation:Md Ezaz Ahmed,Anil Carie Chettupally,Dinesh Mavaluru.All authors reviewed the results and approved the final version of the manuscript.

Availability of Data and Materials:The simulation code is not yet available to public due to extended project proposal.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.