A method comprising: obtaining, by a computing system, a decomposed fluid pipeline network comprising a fluid pipeline network decomposed into a plurality of subnetworks (310, 312, 314), the plurality of subnetworks comprising a plurality of upstream subnetworks (312, 314) and a downstream subnetwork, and a plurality of boundary nodes where an upstream subnetwork (312, 314) overlaps the downstream subnetwork, the boundary node comprising a sink (302, 308) of the upstream subnetwork (312,314) and a source (304,306) of the downstream subnetwork individually executing, by the computing system, and using a first set of control device settings and a first set of boundary conditions, a network solver on the plurality of upstream subnetworks (312, 314) in parallel, wherein the first set of control device settings are initial settings for selected control devices in the upstream subnetwork that control the flow of fluid through the fluid pipeline network and where the first set of boundary conditions comprise fluid property values at a source of the each upstream subnetwork (312, 314), wherein individually executing the network solver comprises repetitively performing a simulation using a pipeline network model to obtain a simulation result comprising a revised set of boundary conditions and control device settings until consecutive simulation results converge to within a threshold and outputting the final simulation result comprising a second set of boundary conditions and a second set of control device settings; executing, by the computing system and using the second set of boundary conditions and second set of control settings, the network solver on the downstream subnetwork by repetitively performing a simulation using the pipeline network model to obtain a simulation result comprising a revised set of boundary conditions and control device settings until the consecutive simulation results converge to within a threshold and outputting the final simulation result comprising a third set of control device settings; testing, by the computing system, whether a threshold has been reached, by simulating the pipeline network using the second and third set of control device settings and evaluating the simulation result with an objective function, wherein the threshold is based on a test for convergence with previous evaluation results, wherein the objective function is selected so as to optimize so as to optimize exploration and production operations of the wellsite system; and presenting, by the computing system, the third set of control device settings based on the result satisfying the threshold. obtaining, by a computing system, a decomposed fluid pipeline network comprising a fluid pipeline network decomposed into a plurality of subnetworks (310, 312, 314), the plurality of subnetworks comprising a plurality of upstream subnetworks (312, 314) and a downstream subnetwork, and a plurality of boundary nodes where an upstream subnetwork (312, 314) overlaps the downstream subnetwork, the boundary node comprising a sink (302, 308) of the upstream subnetwork (312,314) and a source (304,306) of the downstream subnetwork individually executing, by the computing system, and using a first set of control device settings and a first set of boundary conditions, a network solver on the plurality of upstream subnetworks (312, 314) in parallel, wherein the first set of control device settings are initial settings for selected control devices in the upstream subnetwork that control the flow of fluid through the fluid pipeline network and where the first set of boundary conditions comprise fluid property values at a source of the each upstream subnetwork (312, 314), wherein individually executing the network solver comprises repetitively performing a simulation using a pipeline network model to obtain a simulation result comprising a revised set of boundary conditions and control device settings until consecutive simulation results converge to within a threshold and outputting the final simulation result comprising a second set of boundary conditions and a second set of control device settings; executing, by the computing system and using the second set of boundary conditions and second set of control settings, the network solver on the downstream subnetwork by repetitively performing a simulation using the pipeline network model to obtain a simulation result comprising a revised set of boundary conditions and control device settings until the consecutive simulation results converge to within a threshold and outputting the final simulation result comprising a third set of control device settings; testing, by the computing system, whether a threshold has been reached, by simulating the pipeline network using the second and third set of control device settings and evaluating the simulation result with an objective function, wherein the threshold is based on a test for convergence with previous evaluation results, wherein the objective function is selected so as to optimize so as to optimize exploration and production operations of the wellsite system; and presenting, by the computing system, the third set of control device settings based on the result satisfying the threshold.