C $********************************************************************************$

C

C CO$2$ Injection Well Simulator Final Version $(4/4/2024)$

C

C $********************************************************************************$

DIMENSION D$(100),$TF$(100),$TIME$(100),$DEN$(100,100),$VIS$(100,100)$

DIMENSION V$(100,100),$P$(100,100),$T$(100,100),$SH$(100,100)$

OPEN$(2,$FILE$=$'Resulpm.txt$\text{'},$STATUS$=$'UNKNOWN'$)$

OPEN$(3,$FILE$=$'Resultm.txt$\text{'},$STATUS$=$'UNKNOWN'$)$

C

C Parameters definition: Injection pressure $($MPa$)$ and temperature $($C$),$formation surface temperature $($C$),$

C geothermal gradient $($C$/100$m$),$well depth $($m$),$ injection rate $($m$3/$min$),$ tubing internal diameter $($m$),$

C heat transfer coefficient $($W$/($m$2.$C$)),$tubing section area $($m$2),$ velocity $($m$/$s$),$ element length $($m$),$

C number of elements, time step size $($s$),$ specific heat $($J$/($kg$.$C$))$

C

C In the matrices, first dimension is depth and the second time

PINJ$=80.$

TINJ$=-20.$

TS$=24.$

GG$=2\mathrm{.}95/100$

WD$=2800$

QM$=4\mathrm{.}0$

DI$=0\mathrm{.}076$

U$=100.$

A$=3\mathrm{.}1416*$DI$**2/4$

VI$=$QM$/($A$*60)$

DX$=50.$

N$=$WD$/$DX

DT$=3$

TIME$(1)=0.$

C Setting grid depths and temperatures

Do $10$ I$=1,$N$+1$

D$($I$)=$DX$*($I$-1)$

TF$($I$)=$TS$+$GG$*$D$($I$)$

$10$ Continue

C Defining the initial conditions

P$(1,1)=$PINJ

T$(1,1)=$TS

V$(1,1)=$VI

Call PROPERTIES$($P$(1,1),$T$(1,1),$DEN$(1,1),$VIS$(1,1),$SH$(1,1))$

Do $40$ I$=1,$N$+1$

PA$=$P$(1,1)$

TA$=$T$(1,1)$

$30$ Call PROPERTIES$($PA$,$TA$,$DEN$($I$+1,1),$VIS$($I$+1,1),$SH$($I$+1,1))$

DENA$=($DEN$($I$,1)+$DEN$($I$+1,1))/2$

VISA$=($VIS$($I$,1)+$VIS$($I$+1,1))/2$

C Mass balance equation

V$($I$+1,1)=$V$($I$,1)*$DEN$($I$,1)/$DEN$($I$+1,1)$

VA$=($V$($I$,1)+$V$($I$+1,1))/2$

C Momentun balance equation

Call FRICG$($VA$,$VISA,DI$,0.,$DENA,GFP$)$

GPH$=9\mathrm{.}81*$DENA$/1000000.$

DP$=($GPH$-$GFP$)*$DX

P$($I$+1,1)=$P$($I$,1)+$DP

####PLEASE WRITE PYTHON CODE FOR BELOW

###I have a Fortran code that I need to convert to python. Please do not Use CHATGPT or AI tools, I need only code with output...

C Heat Transfer equation

A$1=$A$*$DEN$($I$+1,1)*$V$($I$+1,1)*$SH$($I$+1,1)/$DX

B$1=3\mathrm{.}1416*$U$*$DI

T$($I$+1,1)=($A$1*$T$($I$,1)+$B$1*$TF$($I$+1))/($A$1+$B$1)$

If$($Abs$($P$($I$+1,1)-$PA$).$LT$\mathrm{.}0\mathrm{.}1.$and.Abs$($T$($I$+1,1)-$TA$).$LT$\mathrm{.}0\mathrm{.}1)$goto $50$

PA$=$P$($I$+1,1)$

TA$=$T$($I$+1,1)$

Goto $30$

$50$ Write$(*,*)$D$($I$),$P$($I$,1),$T$($I$,1),$TF$($I$)$

$40$ Continue

C Time loop

Do $60$ J$=2,90$

TIME$($J$)=$TIME$($J$-1)+$DT$/60.$

P$(1,$J$)=$PINJ

T$(1,$J$)=$TINJ

V$(1,$J$)=$VI

Call PROPERTIES$($P$(1,$J$),$T$(1,$J$),$DEN$(1,$J$),$VIS$(1,$J$),$SH$(1,$J$))$

DO $45$ I$=1,$N$+1$

PA$=$P$(1,$J$)$

TA$=$T$(1,$J$)$

$35$ Call PROPERTIES$($PA$,$TA$,$DEN$($I$+1,$J$),$VIS$($I$+1,$J$),$SH$($I$+1,$J$))$

C Mass balance equation

V$($I$+1,$J$)=$V$($I$,$J$)*$DEN$($I$,$J$)/$DEN$($I$+1,$J$)$

C Call Velc$($DEN$($I$,$J$-1),$DEN$($I$+1,$J$-1),$DEN$($I$,$J$),$DEN$($I$+1,$J$),$V$($I$,$J$-1),$V$($I$+1,$J$-1),$V$($I$,$J$),$DX$,$DT$,$V$($I$+1,$J$))$

C Momentun balance equation

DENA$=($DEN$($I$,$J$)+$DEN$($I$+1,$J$))/2$

VISA$=($VIS$($I$,$J$)+$VIS$(1,$J$))/2$

VA$=($V$($I$,$J$)+$V$($I$+1,$J$))/2$

Call FRICG$($VA$,$VISA,DI$,0.,$DENA,GFP$)$

GPH$=9\mathrm{.}81*$DENA$/1000000.$

DP$=($GPH$-$GFP$)*$DX

P$($I$+1,$J$)=$P$($I$,$J$)+$DP

C Heat transfer equation

Call Temp$($DEN$($I$,$J$-1),$DEN$($I$+1,$J$-1),$DEN$($I$,$J$),$DEN$($I$+1,$J$),$V$($I$,$J$-1),$V$($I$+1,$J$-1),$V$($I$,$J$),$V$($I$+1,$J$),$SH$($I$,$J$-1),$

SH$($I$+1,$J$-1),$SH$($I$,$J$),$SH$($I$+1,$J$),$T$($I$+1,$J$-1),$T$($I$,$J$),$TF$($I$+1),$U$,$A$,$DI$,$DX$,$DT$,$T$($I$+1,$J$))$

If$($Abs$($P$($I$+1,$J$)-$PA$).$LT$\mathrm{.}0\mathrm{.}1.$and.Abs$($T$($I$+1,$J$)-$TA$).$LT$\mathrm{.}0\mathrm{.}1)$goto $55$

PA

Step by Step Solution

There are 3 Steps involved in it

1 Expert Approved Answer

Step: 1 Unlock

Question Has Been Solved by an Expert!

Get step-by-step solutions from verified subject matter experts

Step: 2 Unlock

Step: 3 Unlock

Students Have Also Explored These Related Programming Questions!

Q:

Instructor: Jones Student ID Major 1 ECON 2 FIN 3 MIS 4 MRKT 5 MNGT 6 ACCT 7 ACCT 8 ACCT 9 MIS 10 MIS 11 MIS 12 MIS 13 MIS 14 FIN 15 FIN 16 FIN 17 ECON 18 MRKT 19 MRKT 20 ECON Correct Answer...

Q:

RFLPs can include sequence changes created by single-base change, single-base insertion, single-base deletion, or tetranucleotide repeat, as shown in Figure 21.21. For each polymorphism, show a...

Q:

Why do you suppose that forensic DNA analysis relies principally on short tandem repeats (repeat polymorphisms), rather than single- nucleotide polymorphisms, such as that described on page 669 and...

Q:

In the figure are shown three capacitors with capacitances c 1 = 7.46?10 ?6 F , c 2 = 2.48?10 ?6 F , and c 3 = 3.8?10 ?6 F . The capacitor network is connected to an applied potential V a b . After...

Q:

answer the question clearly You are building a flight-control system for which a convincing safety case must be made. Would you assign the tasks of safety requirements engineering, test case...

Q:

Question 1: Describe a principal-agent relationship. a. a principal makes a decision that affects an agent b. an agent makes a decision that affects a principal c. an employee acting on behalf of its...

Q:

A researcher at a marketing firm examines whether the age of a consumer matters when buying athletic clothing. Her initial feeling is that Brand A attracts a younger customer, whereas the more...

Q:

MATHEMATICS FOR MACHINE LEARNING Marc Peter Deisenroth A. Aldo Faisal Cheng Soon Ong Contents Foreword 1 Part I Mathematical Foundations 9 1 Introduction and Motivation 11 1.1 Finding Words for...

Q:

1 (1 point) Electronic options for withdrawals or payments include Question 1 options: A) automatic payments. B) debit and ATM cards C) direct deposits. D) a. and b. E) a., b., and c. SaveQuestion 2...

Q:

Debra wants to know if her sister can be her qualifying child.Which of the following is not a test to be a Qualifying Child? Select one: a. The child must be under age 19 at the end of the year or...

Q:

Estimates for the number of computers still running Windows Server 2003 range from hundreds of thousands to several million even though Microsoft has stopped supporting the product. What are the...

Q:

Scherer Corporation is preparing a bid for a special order that would require 720 liters of U48N material. The company already has 560 liters of this raw material in stock, which originally cost...

Q:

Which of the following $ 1 . 0 0 0 face - value securities has the highest yieds to matulity? A 5 percent coupon bond with a price of $ 1 , 2 0 0 A 5 percent coupon bond with a price of $ 1 . 0 0 0 A...

Q:

Exercise 2-1 Compute a Predetermined Overhead Rate [LO2-1] Harris Fabrics computes its plantwide predetermined overhead rate annually on the basis of direct labor-hours. At the beginning of the year,...

Previous Question Next Question