What industries does newtsim serve?

newtsim serves oil and gas (flow assurance, drilling fluids, subsea pipeline), pharmaceutical and drug discovery, energy storage (battery materials, electrolytes), semiconductor process simulation, green hydrogen (electrolyser and catalyst design), nuclear decommissioning, aerospace and composites, corrosion and materials, carbon capture, agrochemicals, mining and geotechnical, and fermentation and bioprocessing.

How is newtsim different from standard CFD or simulation software?

Standard simulation tools use empirical models at a single scale. newtsim runs the full physics chain from quantum mechanics and molecular dynamics through to plant-scale CFD and structural mechanics, GPU-native, with direct data flow between every scale and no third-party solver dependencies. Constitutive behaviour is predicted from first principles rather than empirically fitted.

What does a fixed-price simulation study include?

A newtsim fixed-price study has a defined scope, clear deliverables, fixed cost agreed upfront, and is delivered in weeks rather than quarters. Studies typically include a scoping call, physics-grounded simulation using the appropriate newtsim solvers, results analysis, and a written report with actionable recommendations.

What is physics-grounded simulation?

Physics-grounded simulation derives material properties and constitutive behaviour from lower-scale physics (molecular dynamics and quantum mechanics) rather than from empirical fitting. This gives accurate predictions for novel materials and extreme conditions outside the range of experimental data, and enables root-cause analysis at the mechanistic level.

How does newtsim approach oil and gas flow assurance simulation?

newtsim runs a continuous physics chain from molecular wax crystallisation to full pipeline transient multiphase flow, with no empirical handoffs between scales. Wax appearance temperatures, gel strength, and deposition kinetics are predicted from molecular dynamics (newtsim Bond), passed directly to mesoscale morphology (newtsim Pack), and then propagated into multiphase CFD (newtsim Stream), giving flow assurance predictions that hold outside the range of laboratory data.

Can newtsim simulate battery electrolyte behaviour and degradation?

Yes. The Bond solver runs molecular dynamics for Li-ion transport and solid-electrolyte interphase (SEI) formation; Pack handles electrode morphology at the mesoscale; these results couple upward to cell-scale electrochemical models. This gives physics-grounded predictions of electrolyte degradation, capacity fade, and thermal behaviour without relying solely on empirical ageing fits.

How does physics-grounded drug discovery simulation differ from docking software?

Molecular docking software treats the protein and ligand as rigid bodies and scores binding poses by geometry. Physics-grounded simulation uses molecular dynamics to capture conformational flexibility and explicit solvation effects that rigid-body docking misses, predicting binding free energy, ADMET properties, and metabolic activation from first principles. This is particularly valuable for targets where induced-fit binding or solvent-mediated interactions drive selectivity.

What makes GPU-native simulation faster than CPU-based CFD?

GPU parallelism maps naturally to the computational kernels underlying molecular dynamics and CFD: large numbers of independent force or flux calculations that can execute simultaneously. newtsim runs all simulation scales GPU-native, avoiding the CPU-to-GPU data transfers that bottleneck traditional solvers which run physics on CPU and only offload selected operations to the GPU. This removes the transfer overhead and allows the full multiscale chain to run without the latency penalties of heterogeneous pipelines.

Does newtsim replace existing simulation tools like ANSYS or COMSOL?

newtsim fills the multiscale gap, connecting molecular-scale physics to the engineering-scale models that tools like ANSYS operate at. It is not a replacement for FEA-only or single-scale workflows where those tools are the right fit. Where newtsim adds value is in problems where constitutive behaviour, material properties, or interfacial physics need to be derived from lower-scale simulation rather than taken from a handbook or empirical fit.

How long does a fixed-price simulation study take?

Most studies complete in 4 to 8 weeks depending on scope. The standard timeline is: scoping call and problem definition in week 1; simulation execution across the relevant physics scales in weeks 2 to 4; results analysis, validation, and report delivery in weeks 5 to 6. More complex multiscale studies with broader operating envelopes may extend to 8 weeks. Scope and timeline are agreed and fixed before work begins.

What deliverables does a simulation study produce?

A newtsim simulation study delivers: a validated multiscale model calibrated to your system; an operating envelope with physics-based safety margins; a written technical report with findings and actionable recommendations; and optionally CFD visualisations or an ML surrogate model (newtsim Neural) for fast design-space exploration by your engineering team.

Does newtsim work for green hydrogen electrolyser and catalyst design?

Yes. The Stream solver handles proton exchange membrane and alkaline electrolyser CFD: multiphase flow, bubble dynamics, and membrane transport. Bond handles catalyst surface reaction molecular dynamics, giving adsorption energies and reaction barriers from first principles rather than from empirical correlations. Neural then produces fast-to-run surrogate models for design-space exploration and stack optimisation, reducing the iteration cost for new catalyst compositions and electrode architectures.

newtsim is a physics-grounded multiscale simulation consulting firm and software platform. We deliver fixed-price simulation studies and platform licences for industrial clients, covering the full physics chain from quantum mechanics and molecular dynamics through to plant-scale CFD and structural mechanics. All solvers run GPU-native.

Does newtsim use artificial intelligence or machine learning?

Yes. newtsim Neural produces ML-trained surrogate models derived from high-fidelity physics simulations. These surrogates run in milliseconds rather than hours, enabling rapid design-space exploration and real-time control without sacrificing physical accuracy. The ML models are trained on physics simulation data, not empirical measurements, so they inherit the accuracy and extrapolation properties of the underlying physics.

How do I get started with newtsim?

Book a free technical scoping call at newtsim.com/demo. Bring a description of your engineering challenge, the physics you think are important, and any data you have. The call is 45 minutes and produces a problem definition, proposed physics chain, and outline scope that forms the basis of a fixed-price proposal.

Can newtsim simulate non-Newtonian fluids?

Yes. newtsim specialises in non-Newtonian flow: yield-stress fluids, viscoelastic materials, shear-thinning and shear-thickening systems, and complex multiphase flows. The Stream solver handles macroscale CFD for these systems; Bond provides the molecular-scale rheology predictions that replace empirical viscosity fits; and the full chain gives accurate behaviour outside the calibrated range of laboratory data.

Does newtsim work for small companies or startups?

Yes. Fixed-price studies are available to clients of any size. The fixed-price, fixed-scope model means you know the cost before committing, making it accessible for smaller organisations without large simulation teams. Many newtsim clients are startups and scale-ups who need physics-grounded analysis early in product development before investing in laboratory or pilot-scale work.

What is the difference between molecular dynamics and CFD?

Molecular dynamics (MD) simulates the motion of individual atoms and molecules to compute material properties such as viscosity, diffusivity, reaction rates, and phase behaviour, from first principles. Computational fluid dynamics (CFD) models bulk fluid flow at the engineering scale using those material properties as inputs. newtsim connects the two: MD provides accurate constitutive inputs that CFD uses for engineering-scale predictions, rather than relying on empirical correlations.

What is a surrogate model in simulation?

A surrogate model (also called a reduced-order model or emulator) is a fast mathematical approximation of a complex simulation. newtsim Neural builds ML surrogates trained on high-fidelity physics simulation data. Once trained, a surrogate evaluates in milliseconds instead of hours, enabling real-time optimisation and design-space exploration. Because surrogates are trained on physics data, they are more accurate than empirical correlations, especially outside the training range.

What deliverables does a newtsim simulation study produce?

A newtsim simulation study delivers: (1) a validated multiscale model calibrated to your system; (2) an operating envelope with physics-based safety margins; (3) a written technical report with findings and actionable recommendations; and optionally (4) CFD visualisations and an ML surrogate model for fast ongoing design-space exploration by your engineering team.

Oil & Gas

Oil & Gas · Case Study

Pipeline Integrity and Remaining Life Assessment: Prudhoe Bay Transit Line

Executive Summary

On March 2, 2006, at approximately 16:30 local time, a BP field technician discovered crude oil pooling across approximately 1.9 acres of North Slope tundra in the Western Operating Area of Prudhoe Bay, Alaska. The source was a 6 mm × 4 mm through-wall perforation at the 6 o'clock position of a 20-inch crude transit line — a textbook bottom-of-pipe pitting signature. By the time Alaska regulators confirmed the final spill volume in March 2008, it stood at 212,252 US gallons (5,054 barrels): the largest crude oil release ever recorded on the North Slope.

The cause was not a sudden failure. PHMSA's formal investigation and the subsequent Congressional hearings established that the corroded segment had been flagged by inline inspection in 1998, with anomalies of 20--50% wall loss at six locations on the precise segment that later failed. No follow-up ILI was scheduled for another eight years. Corrosion inhibitor injection at the gathering centre had been discontinued or severely curtailed prior to the failure; Congressional testimony estimated that BP saved approximately USD 9.6 million by deferring cleaning and inspection. A 150 mm deep layer of consolidated black sludge — iron sulphide scale, paraffin wax, produced sand, and bacterial biomass — had accumulated at the bottom of pipe, creating a completely anaerobic microenvironment that isolated the steel surface from any inhibitor contact. Sulphate-reducing bacteria colonised this environment and drove pitting at 0.5--1.0 mm/year over fourteen years. BP pleaded guilty to a criminal misdemeanour under the Clean Water Act in 2007. Total cost: USD 20 million in fines and penalties, over USD 25 million in remediation, and approximately USD 100 million in deferred production revenue during the eight-day Western Operating Area shutdown.

Had a coupled pipeline integrity simulation been commissioned in the months before the failure — drawing on the 1998 ILI data already in BP's possession — it would have identified precisely this risk. A multiscale corrosion model, applied to the line's known fluid chemistry (61% water-cut, pH₂S sufficient to sustain SRB activity, 90-day pigging interval), would have computed BOP pitting rates of 0.8--1.1 mm/year at the two lowest-velocity segments, identified that the 30 ppm imidazoline dose provided only 41% BOP film coverage under stratified flow conditions, and placed the deepest existing anomaly at 14--22 months to through-wall failure at P50 confidence. The warning would have been unambiguous. The intervention — inhibitor uplift, biocide schedule revision, targeted pigging, and a composite sleeve — would have cost a fraction of the subsequent cleanup.

This assessment carries out exactly that analysis, applied to a 34-kilometre North Slope transit line with an identical failure mechanism profile: elevated water-cut (61%), established SRB biofilm, insufficient inhibitor coverage, and an ILI programme with multi-year gaps. The petrosim multiscale simulation chain — from molecular-scale iron dissolution kinetics through to field-scale corrosion rate mapping — predicts BOP pitting rates of 0.81--1.14 mm/year at two hotspot locations, with the deepest confirmed anomaly (68% wall loss at km 21.7) assessed at 14--22 months to through-wall perforation under the current treatment programme. The recommended intervention package extends median remaining life from 18 months to 51 months.

The 2006 Prudhoe Bay failure resulted in USD 20 million in regulatory fines, over USD 25 million in remediation and cleanup, and approximately USD 100 million in deferred production revenue — a total direct consequence of approximately USD 145 million — on top of further pipeline replacement commitments of USD 60--80 million. A newtsim simulation, applied to the 1998 ILI data already in BP's possession, would have identified the at-risk segment and prevented the failure. The two identified hotspot locations at km 14.3 and km 21.7 represent precisely the type of spatially defined risk that can be instrumented with continuous monitoring sensors and fed into newtsim livesim for real-time corrosion rate tracking — reducing the probability of an undetected through-wall failure to near-zero while providing fleet-level corrosion intelligence across the operator's broader gathering network.

Scenario Background (illustrative reference case)

The fictional operator is a major integrated upstream company with approximately 1.2 million barrels per day (MMbbl/d) of gross production on Alaska's North Slope, ranking it among the top five onshore US producers by daily throughput. The operator holds working interests in 14 producing pads within the Prudhoe Bay Unit (PBU), the largest conventional oilfield ever discovered in North America, with an estimated ultimate recovery exceeding 13 billion barrels since first production in 1977. The asset base includes approximately 480 km of operated transit and gathering pipelines, ranging from 6-inch field gathering lines to 36-inch trunk lines feeding the Trans-Alaska Pipeline System (TAPS) interconnect at Pump Station 1.

The asset in question — hereafter designated Line TL-7 — is a legacy crude transit line installed in 1978 during the original field build-out phase. The line carries produced fluids (three-phase oil, water, and gas) from a cluster of four satellite pads on the western operating area (WOA) to the WOA Central Processing Facility (CPF) for primary separation, gas compression, and crude stabilisation. The pipeline is a continuous, below-grade, thermally insulated corridor buried in permafrost-overlying gravel pad material, with above-grade caribou crossings at three locations. The nominal pipe material is API 5L Grade B (SMYS 241 MPa) for the older 1970s sections and API 5L X52 (SMYS 358 MPa) for sections replaced or extended during the 1994 WOA expansion programme.

Line TL-7 has a documented history of internal corrosion concerns. The 1991 ILI MFL run identified localised metal loss at two low-point locations, which were subsequently derated and managed under a fitness-for-service monitoring programme. A 1998 MFL run confirmed slow further progression at these sites and identified six additional anomalies in the 14--22 km segment that were classified as requiring enhanced monitoring. No further ILI run was performed for eight years following the 1998 survey. By the time of the assessment, the operator's asset integrity team had flagged a 6.8-kilometre segment (km 12.9 to km 19.7) as requiring closer scrutiny based on anomalous differential pressure trends and corrosion coupon mass-loss data at a monitoring station at km 15.4.

Current production on Line TL-7 is 23,400 STB/d of liquid (oil + water), with a gross GOR of 620 Sm³/Sm³. Water cut has risen from 14% in 2018 to 61% in 2023 — a trajectory entirely consistent with the late-life production profile of North Slope reservoirs. The five-year water-cut increase is the single most important driver of the corrosion regime change identified in this study.

Challenge

Line TL-7 is an API 5L X52 carbon steel pipe, 457 mm (18-inch) nominal diameter, 9.5 mm nominal wall thickness (design wall 10.3 mm, with 0.8 mm corrosion allowance consumed), operating at 3.1--4.4 MPa line pressure and 28--42°C fluid temperature across its length. Fluid temperature is highest at the producing pad headers (42°C) and lowest in the 14-km buried corridor approaching the CPF (28°C at the CPF inlet flange during winter operation). Burial depth averages 1.4 m in gravel pad areas and 2.1 m in the insulated tundra corridor. Insulation is 50 mm closed-cell polyurethane foam jacketed in HDPE, providing an overall heat transfer coefficient of approximately 2.8 W/m²K.

The produced water chemistry at TL-7 inlet is representative of WOA formation water: total dissolved solids (TDS) 38,000 mg/L, chloride 22,000 mg/L, sulphate 180 mg/L, bicarbonate 410 mg/L, iron (total) 12 mg/L, and pH in the range 6.1--6.8 at operating conditions. The partial pressure of CO₂ dissolved in the produced gas phase is 0.09 MPa (90 kPa), the primary driver of sweet corrosion under the de Waard--Milliams framework. H₂S is present at trace levels of 8--15 ppm in the gas phase, contributing a pH₂S of approximately 0.003 MPa — below the NACE MR0175 SSC Region 0 threshold but sufficient to influence SRB metabolic activity and electrochemical corrosion kinetics.

Multiphase flow characterisation using the Taitel-Dukler flow regime map confirms stratified wavy flow as the dominant regime in the 6.8 km flagged segment at current operating conditions. The superficial liquid velocity is 0.4--0.7 m/s (below the critical 1.0 m/s threshold for preventing water dropout and BOP accumulation at 61% water-cut) and the superficial gas velocity is 2.1--3.4 m/s. This flow regime results in sustained water holdup at the BOP position, with the free water layer estimated at 80--120 mm depth over a 200° arc at each pipe cross-section — creating an essentially stagnant aqueous corrosion environment at the BOP even when the bulk fluid is moving at 0.5 m/s.

The ILI programme history is the central data integrity concern for this assessment. Following the 1998 MFL run, no further ILI was conducted on TL-7 for eight years, partially mirroring the 1992--2006 inspection gap that preceded the real-world Prudhoe Bay failure. The most recent MFL run was completed in Q4 2022 using a 457 mm Rosen MFL tool with reported sensitivity of ±5% wall thickness for metal loss features exceeding 25 mm². The run identified 47 metal loss features exceeding 20% wall loss. The critical features are:

Feature ID	Location (km post)	Axial length (m)	Max wall loss (%)	Wall remaining (mm)	Classification
F-001	14.32	3.2	42	5.5	Cluster A — Priority 1
F-002	14.41	0.8	31	6.5	Cluster A — Priority 2
F-003	17.85	1.1	28	6.8	Cluster B — Priority 3
F-004	21.68	0.9	68	3.0	Cluster C — CRITICAL
F-005	21.74	0.4	55	4.3	Cluster C — Priority 1
F-006	24.12	2.1	38	5.9	Cluster D — Priority 2
F-007–F-047	Various	Various	20–27	6.9–7.6	Priority 3–4

The deepest feature, F-004 at km 21.68, with 68% wall loss and only 3.0 mm remaining wall, is the critical-path item for remaining life assessment. At the nominal pipe SMYS of 358 MPa (API 5L X52) and current MAOP of 4.4 MPa, through-wall failure at this location would be analogous to the 2006 Prudhoe Bay spill.

The corrosion inhibitor programme at the time of assessment consisted of continuous injection of film-forming imidazoline compound at 30 ppm at the pad headers, combined with monthly batch biocide treatment at 250 ppm glutaraldehyde. Produced water samples from the TL-7 km 15.4 monitoring station returned SRB plate counts of 10³--10⁴ CFU/mL over the past 12 months — confirming an established, resilient biofilm not being controlled by the current treatment. 16S rRNA sequencing identified Desulfovibrio and Desulfobacter genera capable of generating localised pH drops to 4.5--5.0 and local H₂S concentrations up to 200 ppm in the BOP biofilm.

Real-World Basis

This study draws directly and extensively from the BP Prudhoe Bay transit line failure of March 2, 2006, which constitutes one of the best-documented internal corrosion pipeline failures in the public literature and provides the primary model validation dataset for this assessment.

The incident sequence. At approximately 16:30 local time on March 2, 2006, a BP field technician discovered a crude oil release from a 20-inch pipeline in the WOA of Prudhoe Bay. The release point was a 6 mm x 4 mm through-wall perforation at the 6 o'clock position of the pipe — the classic BOP pitting signature of SRB-driven microbiologically influenced corrosion under accumulated sediment and sludge. Initial estimates placed the spill at 267,000 US gallons (6,357 barrels); Alaska regulators confirmed the volume in March 2008 as 212,252 US gallons (5,054 barrels), representing approximately 1.9 acres of tundra contamination — the largest crude oil release ever recorded on the North Slope.

Root cause investigation findings. PHMSA conducted a formal investigation and issued a Corrective Action Order (CAO) in 2006, with a final determination published in July 2011. The investigation established four key findings. First, the failed line segment (GC-2 transit line, WOA Pad 9) had last been inspected by MFL ILI in 1998 — an eight-year gap. The 1998 smart pig run had detected corrosion anomalies in the range of 20--50% wall loss at six locations on the precise segment that subsequently failed, yet no follow-up ILI was scheduled for eight additional years. Second, corrosion inhibitor injection at the GC-2 gathering centre had been discontinued or significantly reduced for an extended period prior to the failure. Congressional testimony (House Subcommittee on Oversight and Investigations, September 7, 2006) revealed that no inhibitor was being injected at GC-2 prior to the failure, and prosecutors estimated BP saved approximately USD 9.6 million by deferring pipeline cleaning and inspection. Third, autopsy of the failed pipe section revealed a 150 mm deep layer of consolidated black sludge at the BOP position over a 40-metre section surrounding the failure point. This sludge layer — composed of iron sulphide scale, paraffin wax deposits, produced sand, and bacterial biomass — created a completely anaerobic microenvironment at the pipe wall, isolating the steel surface from any inhibitor chemical contact. Fourth, microbiological analysis confirmed SRB populations as the primary corrosion agent. The failure morphology — a single discrete perforation with minimal surrounding general wall loss — is pathognomonic of SRB-driven pitting rather than CO₂ sweet corrosion, which produces broader distributed metal loss. The BOP pitting rate at the failure location was retrospectively estimated at 0.5--1.0 mm/year based on the period between 1992 (no anomaly detected) and 2006 failure.

Regulatory and financial consequences. BP pleaded guilty to a criminal misdemeanour under the Clean Water Act in 2007 and paid USD 20 million in fines and civil penalties. Total cleanup and remediation costs exceeded USD 25 million. Production loss during the WOA shutdown — affecting approximately 400,000 barrels per day of throughput — cost an estimated USD 100 million in deferred revenue over the 8-day shutdown period. BP subsequently committed to replacing approximately 26 miles of North Slope transit pipeline and implementing a revised ILI programme with a maximum 5-year inspection interval.

Published validation data. The BOP pitting rate of 0.5--1.0 mm/year and the MIC pitting morphology established by the investigation serve as primary benchmarks for the petrosim mesoscale model.

Simulation Approach

The petrosim platform applies a three-scale simulation chain for corrosion and integrity assessment. The following diagram illustrates the full simulation workflow:

Molecular scale. The molecular-scale stage models iron dissolution kinetics at the steel--electrolyte interface under TL-7's specific produced water chemistry: chloride 22,000 mg/L, pH 6.1--6.8, pCO₂ 0.09 MPa, and trace H₂S. The purpose is to compute the anodic dissolution rate from first principles rather than relying on empirical correlations that were not calibrated for this combination of conditions.

Inhibitor film behaviour is modelled at the same scale. The imidazoline compound is evaluated across a dose range of 10--110 ppm to determine how much of the steel surface the film actually covers under field conditions. At 30 ppm (the current dose), the model predicts 41% surface coverage (59% efficiency). At the recommended 95 ppm dose, coverage rises to 78% (84% efficiency). This non-linear dose-response is the basis for the inhibitor uplift recommendation.

The MIC submodel captures SRB biofilm growth and metabolic H₂S generation calibrated against North Slope produced water microbiome data. SRB metabolic rate depends on sulphate concentration (180 mg/L), temperature (28--42°C, optimal at 37°C), and biocide inhibition. The model predicts a local pH depression of 1.2--1.8 pH units at the biofilm--steel interface, substantially increasing dissolution rate above the bulk-fluid electrochemical prediction. This is the mechanism that converts a manageable general corrosion problem into a through-wall pitting threat.

Mesoscale. A CO₂ corrosion model extended for MIC is applied at 1-metre resolution over the full 34 km line. At each element, the corrosion rate is computed as the product of the uninhibited electrochemical rate, an iron carbonate scale correction, a flow-velocity mass-transfer correction, the molecular-scale inhibitor efficiency, and a spatially varying MIC acceleration factor (ranging from 1.0 in clean sections to 4.2 in established biofilm zones). The rationale for this multiplicative structure is that each factor modifies the base rate independently and the molecular-scale inputs replace the empirical calibration constants that limit commercial tools.

Water accumulation geometry in low-points and reduced-bore sections is resolved at 0.1-metre wall segment resolution. At km 14.3 and km 21.7, the model predicts sustained bottom-of-pipe water holdup of 85--120 mm depth — sufficient to prevent inhibitor film renewal between pigging passes and to sustain the anaerobic microenvironment required for SRB activity.

Field scale. Field-scale outputs aggregate segment-level corrosion rates into a spatial corrosion map overlaid on the ILI anomaly catalogue. A B31.8S fitness-for-service assessment computes remaining strength factor and failure pressure for each ILI feature. Monte Carlo analysis (10,000 realisations) propagates uncertainty in wall thickness measurement (±0.5 mm MFL tool tolerance) and corrosion rate prediction (±15%) into a probabilistic remaining life distribution at each feature location.

Simulation Caveats

Classification: STRETCH. This study employs a multiscale chain that spans reactive molecular dynamics (sub-nanosecond atomistic timescales) through to field-scale remaining life assessment (multi-year), requiring bridging assumptions at each scale transition. Results should be presented as quantitative guidance informing engineering judgement rather than as deterministic failure predictions.

The molecular-scale force field parameterisation for the Fe/H₂O/Cl/H₂S/CO₂ system is validated against bulk iron dissolution experiments but has limited experimental benchmarking for the simultaneous four-species environment at TL-7 conditions. The predicted imidazoline surface coverage fraction carries an estimated ±15% absolute uncertainty at low concentrations (<40 ppm), which propagates directly into inhibitor efficiency predictions.

The MIC submodel applies a spatially averaged SRB activity factor alpha_MIC(x) derived from bulk plate count and 16S rRNA data at a single monitoring station (km 15.4). SRB colonisation in real BOP biofilms is highly localised; actual pitting rates at km 14.3 and km 21.7 may deviate from predictions by ±25--35% if local biofilm structure differs materially from the km 15.4 reference.

The ±0.5 mm MFL measurement uncertainty (1 sigma) on a 3.0 mm remaining wall at F-004 represents a ±17% fractional uncertainty in the remaining life denominator. At P90, the Monte Carlo remaining life for F-004 compresses to as little as 8 months. This uncertainty dominates the output at the critical feature and should be reduced by UT confirmation measurements before scheduling the composite sleeve installation window.

The mesoscale corrosion rate equation CR(x) = CR_dW x F_scale x F_vel x (1 - eta) x alpha_MIC multiplies five independently computed factors. Error accumulation across the chain means the combined ±15% stated prediction band is likely optimistic; independent ER/coupon validation at km 14.3 and km 21.7 would narrow this uncertainty before the next ILI run.

Recommended framing: Present the remaining life predictions as probabilistic risk assessments (P10/P50/P90 distributions) rather than point estimates, and emphasise that the immediate fitness-for-service exceedance at F-004 (P_f = 2.91 MPa < MAOP 4.4 MPa) is an ILI-derived finding that stands independently of the simulation model. The simulation's primary value is in quantifying the corrosion rate trajectory and treatment optimisation, not in the specific remaining life numbers.

Key Predictions / Results

Corrosion pit depth vs time — F-004 baseline vs post-intervention trajectory (km 21.7, Line TL-7)

The following table summarises the key quantitative predictions versus the baseline for the two critical anomaly clusters:

Parameter	Baseline (current regime)	With recommended intervention	Prudhoe Bay 2006 benchmark
General internal CR, mid-line	0.23–0.31 mm/yr	0.08–0.12 mm/yr	0.15–0.25 mm/yr (inhibited)
BOP pitting rate, km 14.3	0.72–0.96 mm/yr	0.18–0.29 mm/yr	0.5–1.0 mm/yr (failure zone)
BOP pitting rate, km 21.7	0.81–1.14 mm/yr	0.21–0.34 mm/yr	0.5–1.0 mm/yr (failure zone)
Median remaining life F-004 (68% loss)	14–22 months (P50)	48–55 months (P50)	Failed within inspection gap
Failure pressure F-001 (42% loss, km 14.3)	4.08 MPa	4.08 MPa (static)	—
Failure pressure F-004 (68% loss, km 21.7)	2.91 MPa	2.91 MPa (static)	—
MAOP	4.4 MPa	4.4 MPa	—
SRB count, km 15.4 monitoring station	10³–10⁴ CFU/mL	<10¹ CFU/mL (target)	Active at failure location
Inhibitor BOP film coverage at 30 ppm	41% surface coverage	78% (at 95 ppm dose)	Inhibitor lapsed at failure
Biocide MIC for established Desulfovibrio biofilm	550–900 ppm (EPS shielding)	<100 ppm (fresh biofilm)	—

ILI measured wall loss vs petrosim simulation — all critical features, Line TL-7

The failure pressure of F-004 at 2.91 MPa falls below the current MAOP of 4.4 MPa — meaning this feature has already exceeded its RSTRENG fitness-for-service limit. This constitutes an immediate integrity concern requiring either pressure reduction to below 2.5 MPa (with a 20% safety factor) or emergency composite sleeve installation. The failure pressure of F-001 at 4.08 MPa exceeds MAOP by only 7%, well within measurement uncertainty, and also warrants urgent action.

The following risk matrix illustrates the consequence versus likelihood ranking:

Inhibitor dose-response. The MD-derived dose-response curve for imidazoline film coverage shows a non-linear efficiency function. Below 40 ppm, BOP film coverage falls below 50% — effectively inadequate under the current stratified flow regime. The recommended minimum dose for adequate BOP protection is 85--110 ppm continuous injection. Batch treatment at 500 ppm monthly is insufficient to maintain film integrity between pigging runs, as the film degrades within 4--7 days of pig pass in the sludge-accumulation zones.

Biocide programme. The current 250 ppm glutaraldehyde monthly batch treatment has failed to suppress SRB populations. The petrosim MIC submodel predicts that glutaraldehyde minimum inhibitory concentration for the established Desulfovibrio biofilm at km 15.4 is 550--900 ppm due to exopolysaccharide shielding. Upgrading to 400 ppm on a 14-day cycle is projected to reduce SRB viable counts to <10¹ CFU/mL within three treatment cycles, reducing the MIC acceleration factor alpha_MIC from 3.1--4.2 to 1.1--1.4 at the two hotspot locations.

Remaining life intervention timeline:

Comparison Methodology

Model outputs are benchmarked against three independent data sources.

The first validation layer is the ILI delta wall thickness comparison. Corrosion rate predictions are compared against the change in wall thickness measured between the 1998 and 2022 MFL runs at the 47 feature locations. For features detectable in both runs, the 24-year delta provides an average corrosion rate of 0.11--0.28 mm/year for general corrosion in non-sludge zones — consistent with the prediction of 0.23--0.31 mm/year at baseline. The BOP clusters (F-001, F-004) were not detectable in the 1998 run, consistent with the interpretation that they initiated after the ILI gap began and the inhibitor programme deteriorated.

The second layer is corrosion coupon and ER probe data. The km 15.4 monitoring station hosts an electrical resistance (ER) probe and quarterly retrieval coupon station. ER probe data over the 18 months prior to the assessment shows a measured general corrosion rate of 0.26--0.34 mm/year — directly within the predicted range of 0.23--0.31 mm/year. Coupon mass-loss data from the same station confirms 0.29 ± 0.04 mm/year average, validating the model at a location with known flow and chemistry conditions.

The third layer is the BP Prudhoe Bay 2006 published benchmark, which provides secondary experimental confirmation. Investigation documents record BOP pitting rates at the 2006 failure location of 0.5--1.0 mm/year in the 14 years preceding the failure. The prediction of 0.81--1.14 mm/year for the uninhibited BOP pitting rate at TL-7 km 21.7 falls within and slightly above this range, which is physically reasonable given TL-7's higher water cut (61% vs. approximately 35--40% at the 2006 failure location) and comparable produced water chemistry. Agreement is within the expected ±25% uncertainty band for BOP pitting rate prediction in MIC-active environments.

The primary validation is the higher-fidelity molecular and mesoscale models benchmarked against the field-scale electrochemical corrosion predictions. The contemporary ER/coupon data and historical MFL delta rates from TL-7 itself validate the field-scale outputs, while the published Prudhoe Bay failure morphology data provides independent secondary confirmation of the modelled pitting mechanism and rate.

Deliverables

Full corrosion rate spatial map (34 km, 1-metre resolution) in GIS-compatible shapefile format overlaid on the TL-7 pipeline alignment sheet, with colour-coded corrosion severity zones and inhibitor coverage efficiency contours
Feature-by-feature remaining life table for all 47 ILI anomalies: P10/P50/P90 time to through-wall (months), current failure pressure (MPa), RSTRENG fitness-for-service classification (acceptable/monitor/repair), and recommended action
Fitness-for-service assessment report per ASME B31.8S and DNVGL-RP-F101 for the 12 critical features exceeding 40% wall loss, including burst pressure calculations, MAOP derate recommendations, and permanent repair options (composite sleeve vs. full spool replacement)
Emergency action notification for F-004 (km 21.7): formal FFS classification as failed (P_f = 2.91 MPa < MAOP 4.4 MPa), recommended MAOP reduction to 2.5 MPa and urgent composite sleeve installation within 90 days
Inhibitor and biocide treatment optimisation report: MD-derived dose-response curves, recommended injection rates as function of water cut, pigging interval, and flow rate, with cost estimates for chemical uplift
Revised ILI scheduling recommendation: next MFL run recommended within 18 months (no later than Q2 2025), with simultaneous UT grid scan at km 14.3 and km 21.7 to confirm post-repair wall thickness
Monte Carlo uncertainty report: sensitivity ranking of input parameters (most sensitive: wall thickness measurement uncertainty, inhibitor film coverage, SRB activity factor; least sensitive: CO₂ fugacity at current pCO₂)
Executive risk matrix (consequence x probability) for all 47 anomaly clusters, formatted for operator's integrity management plan (IMP) submission to PHMSA
Petrosim simulation data archive (ReaxFF MD trajectories, mesoscale output tables, field-scale GIS files) in standard formats for operator's data management system

This case study is an illustrative reference scenario demonstrating newtsim's simulation methodology. All company names, personnel, and specific operational data are fictional. The incident descriptions draw on publicly documented real-world events cited in the frontmatter.