A Fluid-Solid-Magnetic Coupling Algorithm of Internal Crack Growth in the Weld of Oil and Gas Pipelines

In order to characterize the dynamic process of the crack growth in the weld of oil and gas pipelines, a mathematical model of fluid-solid-magnetic multifield coupling was constructed in this paper. Based on this model, the bidirectional fluid-solid coupling and unidirectional magnetic structure coupling caused by the weld deformation were achieved by dynamic application of the fluid permeation pressure, calculating the internal crack growth in the pipe weld, reconstructing the computational grid near the internal crack, and discussing the characteristics of the magnetic leakage field in the process of the internal crack growth in pipe weld.Thus, a fluid-solid-magnetic coupling algorithm for the internal crack growth in pipe welds considering fluid permeation pressure is established. According to the characteristics of the internal crack opening distance, internal crack growth length, crack tip energy release rate, peak values of magnetic induction intensity level, and vertical component, the process of the internal crack growth is measured.The results show that the fluidosmotic pressure accelerates the process of the internal crack growth and this algorithmcan solve the problem of the characterization and evaluation of crack growth in pipe welds under fluid-solid-magnetic coupling action.


Introduction
Oil and gas pipelines are important ties to connect the production, transportation, and marketing as the transportation tool of petroleum and natural gas.Petroleum and natural gas both belong to state key supervision hazardous chemicals, which are inflammable, explosive, and toxic substances.Hence, it is easy for them to cause the accident of fire and explosion upon leakage and ignition source.In addition, the oil and gas pipeline used for transporting oil and gas resources are a major hazard because of the long pipeline, high pressure, and changeable environment.Hence, the safety management, detection, and protection of oil and gas pipelines are key technologies that are involved in the prevention and control of national production accidents.In the leakage and explosion accidents of oil and gas pipeline, one of the most critical risk factors affecting its safety is the cracking of the pipe weld.Many catastrophic accidents originate from the rapid expansion of the crack under the internal pressure of the pipeline fluid, and the cracking of the pipeline is a dynamic process.Therefore, the research on detecting and evaluating the in-service inspection of weld crack propagation is very important in the field of engineering practice.Besides, the magnetic flux leakage (MFL) technology which was originated from 1960s integrates nondestructive testing, magnetic physics, and metallography and has been widely used to detection in the storage tank and pipelines because of its high accuracy, reliability, and easy automation [1,2].
Because the fluid penetration pressure has a direct impact on the crack in the inner wall of the pipe, it is more dangerous than that caused by the outer wall crack.When there is defect damage on the inner wall of the weld of pipelines, the propagation of the crack defect is aggravated under the erosion of the fluid due to the space of the pipeline weld joint where the leakage defect of the conveying fluid is located, which is the fluid-solid coupling effect between the defect-space contact surface of the inner wall of the pipeline weld seam and the conveying fluid.During MFL detection, such defect expansion is caused by fluid-structure interaction between the defect-space contact surface in inner wall of the weld and the transporting fluid, and this is also accompanied with multifield coupling problem of magnetic force fields.Many researches about the fluid-solid coupling in pipelines have been reported.Professor Zheng Jinyang et al. proposed a fluid-solid coupling algorithm with computational stability to investigate the fluid-solid coupling problem of deformation and fracture in pipeline explosion consequences simulation [3].Professor Feng Weimin et al. analyzed the fluid-solid coupling phenomenon of pressure pipelines and constructed the finite-element mathematical model of fluid-solid coupling, simulated the transition process caused by valve switch by using the finite-element software ADINA, and performed the numerical simulation calculation to fluid-solid coupling phenomenon of straight pressure pipelines under different constraints [4].However, there are few research papers reported on the fluid-solid-magnetic multifield coupling studies by using the application of magnetic force field.E. Aursand et al. developed a coupled fluid-structure interaction model, which consisted of a homogeneous equilibrium flow model, coupled with the Span-Wagner equation of state and including solid-phase formation, and a finite-element model of the pipe taking into account large deformations and fracture propagation through a local fracture criterion [5].H. O. Nordhagen et al. modeled the pipe material and fracture propagation using the finite-element method with a ductile fracture criterion.The finite-volume method has been employed to simulate the fluid flow inside the pipe, and the resulting pressure profile was applied as a load in the finiteelement model.Choked-flow theory was used to calculate the flow through the pipe crack [6].
The preliminary experimental research on different degrees and positions of typical cracking defect in the weld of oil and gas pipelines was carried out, and the magnetic induction intensity distribution in the leakage magnetic field caused by the pipe weld cracks was obtained, and then the surface remanufacturing was studied [7,8].However, this study only involved the static MFL testing of the weld cracks of pipelines.The research on the dynamic MFL in the crack growth process was not discussed in detail.Therefore, a simulation algorithm of the internal crack growth in the weld of oil and gas pipelines based on the fluid-solid-magnetic multifield coupling was established in this paper.And this algorithm could achieve the detection of internal cracks in the weld of serving pipelines by using MFL technique and characterize and evaluate the internal crack growth in the weld of pipelines under the fluid-solid-magnetic multifield coupling effect.The realization of the algorithm can provide a theoretical basis for the detection and evaluation of weld cracks in-service pipeline by using MFL method and the remanufacture, repair, or replacement of oil and gas pipelines.

Fluid-Solid-Magnetic Multiphysical Field Model
According to the principle of MFL detection of pipeline weld [7], the fluid-solid-magnetic multiphysical field model of internal crack growth in the weld of pipelines was constructed.Firstly, the structural model of weld was set up.According to the material and parameters of the pipe weld (such as pipe diameter D, pipe wall thickness , weld reinforcement height   , and weld width   ), the stress analysis model of the weld structure of pipelines was established with the plain strain model and the two-dimensional PLANE182 structural unit in finite-element software.Secondly, the initial internal crack was defined in advance.The initial length of internal crack was L b and the position of crack tip was T. The expansion route was defined in advance from T to the weld center on the outer wall of pipelines along the radial direction.The preexpansion route length was L p .As shown in Figure 1, the internal crack distribution along the radial inner wall was calculated by taking the position of 1/2 arc length from the fusion line to the center of the welding path as an example. p was discretized into INTER202 unit.And at INTER202 unit, TARGE169 was chosen as the target unit and CONTA171 was chosen as the contact unit to create the contact relationship.
After that, the excitation structural model outside the pipelines was constructed.According to the structural size, the weld and excitation structural model of pipelines is shown in Figure 2. In the excitation structure, Nd-Fe-B was chosen as the permanent magnet and ingot iron was used as armature and pole shoes.The two-dimensional PLANE13 magnetic field unit was applied to construct the magnetic field analysis model for the excitation structure according to the excitation structural size.Fluids at internal cracks in pipelines and weld were discretized into the two-dimensional PLANE13 magnetic field unit, and the outer spaces of the weld and surrounding air domain of excitation structure were discretized into the two-dimensional PLANE13 magnetic field unit.Subsequently, the fluid-solid-magnetic multifield coupling finite-element gridding map of internal crack growth in the weld of pipelines was established (Figure 3).

Mathematical Model of
Fluid-Solid-Magnetic Multifield Coupling

Mathematical Equations.
In the multiphysical field model, unit discretions of structure and magnetic field were independent mutually.The fluid-solid-magnetic multifield coupling was constructed to realize subsequent coupling calculation of internal crack growth in the weld.The construction process was as follows.
(a) The Fluid-Solid-Magnetic Multifield Coupling Finite-Element Equation Was Established where  is the stiffness matrix of structure,  is the displacement vector of nodes,  () is the load vector of nodes,  m() is the magnetic stiffness matrix,  is the magnetic potential vector of nodes, and  is the magnetic flux density matrix. () in ( 1) is expressed as follows: where  () is the fluid pressure load on the inner wall and at the internal crack of pipelines and  mg() is the magnetic load.
(b) The Bidirectional Fluid-Solid Coupling Form Was Proposed and the Fluid Pressure Load Was Applied Dynamically.
With the internal crack growth on inner wall of pipelines, fluid penetrated into cracks and fluid osmotic pressure was generated.The application of  () in (2) accelerated the internal crack growth, which caused deformation of the weld structure.The displacement vector of nodes () was changed accordingly.As a result, the internal cracks propagated synchronously after structural deformation and further permeation of fluid, which influenced distribution of  () and further reaccelerated growth in the internal crack.The pressure distribution of the fluid was a dynamic process with the growth of the internal crack,  () was the function of pipeline displacement (), and  () was also applied dynamically.This was the bidirectional fluid-solid coupling process between internal fluid and weld structure of pipelines, which was also known as the strong coupling relationship.
(c) The Unidirectional Magnetic Structural Coupling Form Was Proposed.During internal crack growth in the weld of pipelines, fluid domain in internal cracks in the pipelines and the weld as well as outer spaces of the weld and surrounding air domain of excitation structure changed.Changes of  affected magnetic field distribution in the weld of pipelines during the internal crack growth. m() in (1) was the function of .Magnetic force of the magnetic field influenced the weld structure of pipelines slightly, and  mg() in (2) could be neglected.This was the unidirectional magnetic coupling process from the weld structure to the magnetic field, and it was known as the weak coupling relationship.

VCCT Technique.
There are three cracking modes or types of structural material cracks under the external loads, which are opening mode, sliding mode, and shearing mode.
In the opening mode (usually called Mode-I crack), the external load is perpendicular to the crack plane, and the crack tip opens along the direction which is vertical to the external load.Longitudinal crack of a cylindrical shell under hoop tensile stress is one of Mode-I cracks.In this study, the longitudinal direction weld crack of the oil and gas pipeline which is developed from the circumferential tension stress caused by the fluid pressure is a Mode-I crack.
In fracture mechanics, the fracture criterion includes stress intensity factor criterion   <   ( = I, II, III), energy release rate criterion   <   ( = I, II, III), j-integral criterion   <   , and CTOD (crack tip opening displacement) criterion CTOD < .The energy release rate criterion was applied to judge the crack propagation in this study [9].Nowadays, many scholars all over the world have proposed several FEM-based numerical methods for computing the strain energy release rate, and virtual crack closure technique (VCCT) is the most typical one based on fracture mechanics and it is an important tool for studying crack growth.In the FEM model of the cracked pipe weld in this study, VCCT was applied to compute the strain energy release rate , which was a metric parameter of fracture toughness denoting the energy rate of crack growth [10][11][12][13][14].

Mathematical Problems in Engineering
In this paper, the pipe weld is simplified into a plane strain model.For Mode-I crack, the parameter   under critical plane strain is expressed as The basic idea of VCCT is that the work of the crack developing from a to  + Δ is equivalent to the work of the crack closing from  + Δ to a. Figure 4 shows a crack with the length "a" and the virtual increment "Δ".The closure integral of Mode-I crack is where  (1)  yy is the stress component at crack tip when the length of the crack is d, Δ] (2) is the relative opening displacement of the point in the virtual crack plane, Δ is the virtual crack increment, and   is the component of the strain release rate of Mode-I crack.
In consideration of the FEM meshes near the crack tip, the work done by the stress in the virtual crack line equals the product of the node force and its displacement.The basic hypothesis of VCCT is that the opening displacement at the virtual crack tip is close to that at the initial crack tip.Therefore, for the FEM meshing model in Figure 5, the strain energy release rate is where F y1 is the force at node 1 and Δ] 3,4 is the relative vertical displacement between node 3 and node 4.
By comparing the energy release rate   obtained by FEM computation with the critical energy release rate of the material   , the crack propagation in pipe welds could be predicted: if   ≥   , then the crack will grow.

Dynamic Application of the Fluid Permeation Pressure.
(1) The fluid pressure in the final loading step which was applied on the inner wall of pipeline was  end ( end = 30MPa).
(2) During the solution process, the automatic load step is turned on, and the fluid pressure applied at each load step was  i .In the first load step, apply the fluid pressure   =  end /1000.If it converges easily,   =  end /(100 ∼ 1000), while if it converges slowly,   =  end /(1000 ∼ 10 7 ).
(3)  end was applied on the TARGE169 unit and CONTA171 unit.As shown in Figure 6, the initial points (A and B) and final points (C and D) of the fluid osmotic pressure were defined.TARGE169 unit and CONTA171 unit were at the same position.
(4) In the solving process, the internal crack growth length ( s ) was verified and the fluid osmotic pressure ( s ) was applied.
If  s = 0, no internal growth was observed.The fluid osmotic pressure at the initial internal crack length ( b ) was as follows: where  is the radial distance variable from the inner wall at the internal crack to the weld center at the outer wall of pipelines.
The fluid osmotic pressure at  p was If  s > 0, the internal crack began to propagate.The fluid osmotic pressure at the internal crack  b +  s was The fluid osmotic pressure at the nongrowth position ( p −  s ) of preexpansion route was In ( 6)∼( 9), the fluid osmotic pressure ( s ) and fluid pressure on the inner pipeline wall (  ) formed the fluid pressure load ( () ) in the process of internal crack growth.The magnetic load ( mg() ) was overlooked, and the load vector of the node ( () ) in ( 1) was formed again in each loading step.The dynamic application of fluid pressure load was accomplished, which would further influence the weld structural deformation ().In this way, the bidirectional fluid-solid coupling between the fluid in the pipeline and weld structure was realized.

Application of Displacement Boundary Conditions.
The displacement constraint was applied.Radial and circumferential degrees of freedom (DOFs) were restrained on 0 ∘ and 180 ∘ circumferential directions of pipelines, while the annual DOF was restrained on the 270 ∘ circumferential direction of pipelines.

Application of Magnetic Field Loads and Magnetic Field Boundary Conditions. (1)
The coercivity of permanent magnet was 919000 A/m and it was applied along the height of permanent magnet.
(2) The far-field unit was used, and only computational domain with finite size was required to be established.INFIN110 far-field unit was chosen and far-field boundaries were applied surrounding the computational domain.The condition of convergence was L2 norm ‖{R}‖ 2 = (∑ R 2 i ) 1/2 , as defaulted in FEM software.

Fluid-Solid-Magnetic Multifield Coupling Algorithm
The flowchart of the proposed fluid-solid-magnetic coupling algorithm is shown in Figure 7.The major steps are as follows.
(1) Construct the fluid-solid-magnetic multiphysical field model of internal crack growth in the weld of pipelines.
(2) Construct the fluid-solid-magnetic multifield mathematical model of internal crack growth in the weld of pipelines.
(3) If   ≥   , the crack grows,  s > 0; then calculate internal crack growth in the weld of pipelines.
(4) Reconstruct computational domain grid surrounding internal cracks in the weld of pipelines.
(5) Analyze the MFL field in the internal crack growth process.
(6) Compare the internal crack expansion length ( s ) and preexpansion route length ( p ), and judge whether the weld in pipelines is broken.The sketch map of  s and  p is shown in Figure 8.
(b) If  s ≥  p , the weld is broken.End the calculation.
(7) Postprocess the fluid-solid-magnetic multifield coupling results of internal crack growth in the weld of pipelines.

Calculation Parameters.
The excitation structure contains a closed magnetic circuit consisting of the permanent magnet, armature, pole shoe, air gap, and pipeline weld.The dimensions and materials of the magnetization structure are presented in Table 1.The parameters of the pipe weld and the initial parameters of the internal crack are shown in Table 2.In this paper, the internal crack distribution along the radial inner wall was calculated by taking the position of 1/2 arc length from the fusion line to the center of the welding path as a numerical example (see Figure 1).Q235 steel was used as the material of the pipe, and the constitutive model of the pipe weld material was a linear strengthening elastic-plastic model, and the mechanical properties of this material are shown in Table 3.

Internal Crack
Growth.The internal crack growth was computed by the VCCT internal crack growth algorithm and the energy release rate criterion in fracture mechanics.At the current loading step , the energy release rate (  ) at crack tip was calculated.And then, compared with the critical energy release rate ( IC ) of the selected pipeline material, if   ≥  IC , the internal cracks propagated.
Figure 9(a) presents the distribution of the fluid osmotic pressure in the calculation process of the internal crack propagation.The load step is  = 68, and the fluid osmotic pressure is   =  s = 13.7511MPa.However, when the fluid osmotic pressure is ignored ( s = 0), the distribution of the fluid pressure is displayed in Figure 9(b).And the load step is  = 58 and the fluid pressure is   = 14.6475MPa.
The result of the load step (progressive application of the fluid osmotic pressure) in the crack propagation was extracted, and the node coordinates were updated according to the deformation amount of the inner wall structure in pipelines during the internal crack propagation.The internal crack growth process is exhibited in Figure 10.
With the internal crack growth on inner wall of pipelines, the fluid penetrated into cracks and fluid osmotic pressure was generated, resulting in the acceleration of the internal crack propagation and the deformation of the pipeline weld structure.Once the structure of the pipe weld was deformed, the internal cracks were also extended synchronously and the fluid penetrated further, which affected the distribution of the fluid permeation pressure and further accelerated the expansion of the internal cracks.Thus, the bidirectional fluidsolid coupling process was achieved between the internal fluid and weld structure of pipelines.

Gridding Reconstruction for the Internal Crack.
The key points generated at the crack profile nodes according to the opening distance of internal cracks  o and extended length  s were connected into a surface, which was then meshed and given the properties of air units.Thus, the air gridding located at the crack was reconstructed.The gridding reconstruction is a process of dynamic internal crack growth and dynamic gridding reconstruction, and it is the key technology of the fluid-solid-magnetic multifield coupling   algorithm for the internal crack growth in the weld of oil and gas pipelines.The gridding of the fluid domain inside the pipelines and in the crack of the weld was reconstructed according to the deformation volume of the weld and inner wall structure during the internal crack growth under the action of fluid-structure interaction.Besides, the gridding of the outer space of the weld and surrounding air domain of excitation structure was also reconstructed according to the same deformation under the action of fluid-structure interaction.The reconstructed gridding maps are shown in Figure 11.After the mesh reconstruction for fluid domain in pipelines and air domain outside the pipelines, the magnetic stiffness matrix ( m() ) in (1) was regenerated in each loading step, implying that  m() is related to .Thus, the coupling of unidirectional magnetic structural was implemented.

Analysis on the Influence of Structural Deformation.
Based on the internal crack growth and reconstruction of each loading step, the magnetic pole direction in the excitation structure was set.Thus, the magnetic return path between the weld and excitation structure was formed.The MFL field during internal crack growth in the weld was computed by the nonlinear vector potential method.The opening distance of internal cracks  o and extended length  s at each loading step () from the growth beginning to the buckling failure was extracted, and Figures 12 and 13 indicated   −  o and   −  s curves obtained at the internal cracks could be used to compare the influence of the fluid osmotic pressure  s on the internal crack growth under condition of considering and not considering the fluid osmotic pressure.It was observed from Figures 12 and 13 that the opening distance of internal cracks  o , extended length  s obtained at the last loading step before the buckling failure was   S  = 3.6506mm and   S s = 3.1892mm, respectively, when the fluid osmotic pressure  s was applied, whereas the opening distance of internal cracks  o and extended length  s obtained without the fluid osmotic pressure  s was   i  = 3.3909mm and   i s = 3.0009mm, respectively.This phenomenon illustrated that the opening distance of internal cracks and extended length was large, when the fluid osmotic pressure  s was considered.
Besides, the energy release rate   at the crack tip was calculated, namely, the command flow programmed by using the finite-element software under different pressures and growth lengths, and the results were displayed in Tables 4 and  5 and Figure 14.It was found that   <   (63.02N/mm), and the internal crack could not expand, when the internal crack entered the stage of unexpanded growth without considering the fluid osmotic pressure.However, if   ≥   , the internal crack began to expand and entered the stage of steady-state growth.When the   value continued to increase, the crack entered the stage of accelerating growth, and it would soon enter the unstable fracture stage.In this stage, only a small load was applied to keep the pipe damaged until it broke.Consequently, the curvilinear trends in Figures 12-14 can characterize the process of crack growth and the risk level of pipelines.
In order to further compare the growth degree of internal cracks under condition of considering and not considering the fluid osmotic pressure, the  s values of steady-state and accelerating growth stages were fitted from Figures 12 and 13.The fitting curve and equation are obtained in Figure 15 and fitting equation ( 7)∼ (13).
Fitting equation of fitting curve4 without considering  s in the stage of accelerating growth is  = −1.18E−71  (−/0.0992)+ 3.2031E −72  (/0.0888)+ 0.5551 (14.2563 ≤ x ≤ 14.6808) (13) In order to compare the  s growth ratios in the same pressure range, ten further ranges were divided from 12.8034MPa to 13.7703 MPa, and the pressure increment delta was set at  = 0.0967MPa.Thus, the data of  s growth ratios were obtained at the same pressure and the results were listed in Table 6, and then the comparative histogram was drawn in Figure 16.It was found from Table 6 that  s values  were approximately 1.1∼10.7,when the pressure increased gradually.
The  s growth ratio was defined as follows: where     is the  s growth of the internal crack with considering the fluid osmotic pressure and     is the  s   growth of the internal crack without considering the fluid osmotic pressure.welds during the process of internal crack growth was affected by the deformation of pipe weld structure.Consequently, the horizontal component of the magnetic induction intensity   and vertical component of the magnetic induction intensity   showed an increasing trend of whether the fluid osmotic pressure  s was considered or not.

Analysis on the Influence of Magnetic Field
In order to discuss the influence of the fluid osmotic pressure  s at internal cracks on magnetic induction intensity of MFL field, the peak horizontal component of the magnetic induction intensity   and the peak vertical component of the magnetic induction intensity   were extracted at each loading step () from Figure 17.Then the   −   peak and   −   peak curves were obtained as shown in Figure 18.
It was observed from Figure 18 that the peak horizontal component of the magnetic induction intensity   and peak vertical component of the magnetic induction intensity   obtained at the fluid osmotic pressure  s was   S  = 0.2712T and   S  = 0.1093T, respectively.However,   and   obtained without considering the fluid osmotic pressure  s were   i  = 0.2395T and   i  = 0.0776T, respectively.Consequently, the peak horizontal component of the magnetic induction intensity   and peak vertical component of the magnetic induction intensity   detected by the excitation structure was relatively large, when the fluid osmotic pressure  s was considered.Figures 12,13, and 18 that the opening distance of internal cracks  o , extended length  s , peak horizontal component of the magnetic induction intensity   , and peak vertical component of the magnetic induction intensity   increased slowly during the steady-state growth stage, while the four eigen values ( o ,  s ,   , and   ) increased greatly during the accelerating growth stage.Only a small load was applied to keep the pipe damaged until it broke in the last loading step before the buckling failure.In addition, when the fluid pressure (  ) was the same,  o ,  s ,   , and   with considering the fluid osmotic pressure were higher than those without considering the fluid osmotic pressure.The four eigen values ( o ,  s ,   , and   ) could be used to measure the internal crack growth.And the fluid osmotic pressure accelerated the internal crack growth in pipelines.

Further Discussion. It was revealed from
It was seen from Tables 4 and 5, the fluid pressure at the beginning of the internal crack growth was 12.8034 MPa without considering the fluid osmotic pressure, whereas the fluid pressure at the beginning of the internal crack growth was 12.5127 MPa with considering the fluid osmotic pressure.This result indicated that the application of the fluid osmotic pressure was more likely to cause the growth of internal cracks.In addition, when the fluid osmotic pressure was considered, the fluid pressure difference (from the beginning of the extended loading step to the loading step before the pipeline was buckling failure) was  It was also found from Tables 4 and 5 that when the fluid osmotic pressure was considered, the  s values difference (from the beginning of the extended loading step to the loading step before the pipeline was buckling failure) was   S s d = 3.1892 − 0.2391 = 2.9501mm.However, when the fluid osmotic pressure was elided, the  s values difference was   i s d = 3.0009 − 0.2379 = 2.763mm.
In order to measure the speed of crack growth in detail, we defined the internal crack propagation velocity  c , which was expressed as follows: where    is the  s values difference (from the beginning of the extended loading step to the loading step before the pipeline was buckling failure) and   is the fluid pressure difference (from the beginning of the extended loading step to the loading step before the pipeline was buckling failure).It was clear that the internal crack propagation velocity   could be calculated by using (15).When the fluid osmotic pressure was considered, the internal crack propagation velocity was     =   s s d / When the fluid osmotic pressure was considered, the  o and  s values were large, thus resulting in the higher peaks of   and   .This phenomenon explained the curve variations in Figure 18.When the fluid osmotic pressure was considered, peaks of   and   detected by the excitation structure were relatively large due to the large  o and  s values with considering the fluid osmotic pressure (see Figure 12).

Conclusions
(1) When there are some internal cracks in the pipe weld, the fluid carried by the pipe invades the pipe weld where the internal cracks are located.In order to solve this problem, a fluid-solid-magnetic multiphysical field model was established by using the finite-element software.The constitutive model of the pipe weld material was a linear strengthening elastic-plastic model.Based on the mathematical model of fluid-solid-magnetic multifield coupling, the bidirectional fluid-solid coupling and unidirectional magnetic structure coupling caused by the weld deformation were achieved by using the grid reconstruction and dynamic application of the fluid permeation pressure in this paper.And a new fluid-solid-magnetic coupling algorithm for calculating the crack growth in weld of oil and gas pipelines was built.
(2) The fluid-solid-magnetic coupling algorithm considered that the influence of the fluid osmotic pressure was more suitable to predict the engineering practices.According to the coupling results of this example, through the calculation of the energy release rate at the crack tip, the internal crack propagation could be divided into three stages: the stage of unexpanded growth, the stage of steadystate growth, and the stage of accelerating growth.Based on the calculation results, the  s values of  s growth ratios obtained at the same pressure were approximately 1.1∼10.7,and the internal crack propagation velocity with considering the fluid osmotic pressure was 1.6 times larger than those without considering the fluid osmotic pressure.Thus, it can be seen that the velocity of the internal crack growth would increase when the fluid osmotic pressure was considered.
(3) The fluid-solid-magnetic coupling algorithm could be used to characterize the extent of the internal crack growth, and it could also be applied to other conditions, such as internal cracks at different positions of the pipeline weld, different lengths of internal cracks, and multiple cracks in the engineering practice.The algorithm provides a theoretical basis for the detection and evaluation of internal cracks in weld of in-service pipeline by the MFL method and for the reproduction and repair of the pipeline.

Figure 1 :
Figure 1: Diagram of the initial internal crack of pipe weld.

Figure 4 :
Figure 4: Diagram of the virtual crack closure method.

Figure 5 :
Figure 5: Diagram of the release rate of strain energy calculated by the virtual crack closure method.

Figure 6 :
Figure 6: Diagram of the starting and ending points for the application of osmotic pressure.
Without considering the fluid osmotic pressure

Figure 9 :
Figure 9: With considering and without considering the fluid osmotic pressure.

Figure 10 :Figure 11 :
Figure 10: Process of the internal crack growth.

Figure 12 :
Figure 12:   −  o and   −  s curves with considering the fluid osmotic pressure.

Figure 13 :
Figure 13:   −  o and   −  s curves without considering the fluid osmotic pressure.
crack tip (N/mm)  s -Stage of unexpanded growth   -without considering  s -Stage of unexpanded growth   -with considering  s -Stage of steady-state growth   -with considering  s -Stage of accelerating growth   -with considering  s -Stage of steady-state growth   -without considering  s -Stage of accelerating growth  

Figure 14 :Figure 15 :
Figure 14:   −   curves obtained under condition of considering and not considering the fluid osmotic pressure.
= 1.2576MPa.However, when the fluid osmotic pressure was elided, the fluid pressure difference was   i d = 14.6808−12.8034= 1.8774MPa.Therefore, the fluid pressure difference was small when the fluid pressure penetration was considered due to the result of  s d <   i d .Besides, this result could lead to a large crack propagation velocity and a few extension steps.

𝑃 s d = 2 .
9501/1.2576 = 2.35mm/MPa.However, when the fluid osmotic pressure was elided, the internal crack propagation velocity was    =    s d /   d = 2.763/1.8774= 1.47mm/MPa.Therefore, the internal crack propagation velocity with considering the fluid osmotic pressure was 1.6 times larger than those without considering the fluid osmotic pressure.

Table 1 :
Magnetized structural dimensions and materials.

Table 3 :
Mechanical properties of Q235 steel.

Table 4 :
Energy release rate   with considering the fluid osmotic pressure.

Table 5 :
Energy release rate   without considering the fluid osmotic pressure.

Table 6 :
s values of  s growth ratios at the same pressure.