## The Boundary Element Method Simplifies Corrosion Simulation

##### Bertil Nistad February 17, 2016

In version 5.2 of COMSOL Multiphysics, we offer a new feature for simulating corrosion in slender structures. This significantly speeds up the total time spent when working with structures such as oil platforms. By using the boundary element method (BEM) and specialized beam elements in the *Current Distribution on Edges, BEM* interface, there is no longer a need for a finite element mesh to resolve the whole 3D structure, saving time for large corrosion problems consisting of slender components.

### Comparing the Finite Element and Boundary Element Methods

When solving physical equations on a real, 3D geometry, numerical methods require you to discretize the geometry. This involves breaking up the geometry into elements. At certain points in each element, the equations are solved and the relationship between the elements is described so that a finite set of equations is created.

If you choose to use a finite element method (FEM), you have to discretize the volume between the boundaries of your geometry, since the finite elements describe the relationship between the neighboring elements only. If you choose to use a boundary element method (BEM), the equation system is set up so that you only need to discretize your boundaries in the geometry. There are then many fewer elements; however, the physical equations connect every boundary element with every other boundary element, unlike finite elements on a volume, which are connected by equations only to their neighbors. Both methods have advantages and disadvantages. This blog post discusses an example where the BEM demonstrates its efficiency.

### Studying Corrosion in a Complex Structure

We wrote about the importance of corrosion management in oil platforms in a previous blog post. Such structures usually stand on a seabed with a jacket structure or float above the oil field as a floating production unit. Such jacket structures are generally quite large and complex and include a lot of pipes, steel legs, and hopefully a number of sacrificial anodes, which oxidize to protect the structure from corrosion.

When working with large structures, such as a jacket structure, there may be several hundred or even thousands of anodes that are placed quite close to the structural members. If you want to use a volumetric (3D) mesh on such structures, you will end up with a lot of surfaces and small mesh elements due to the curved surfaces and the thin gaps between the anodes and the structural parts. In order to quickly assess the cathodic protection scheme of a large jacket structure, it can be useful to utilize the boundary element method with the edge elements provided in the Corrosion Module.

*Electric potential of steel versus the Ag/AgCl electrochemical potential on a jacket structure submerged in seawater, calculated with the* Current Distribution on Edges, BEM *interface.*

The figure above demonstrates results from a model containing more than 3000 different edges. Each edge is given a radius and an appropriate current density is specified depending on whether it is made from steel (protected cathode) or represents a sacrificial anode. It is easy to combine several different dimensions and current densities using the selection capabilities in COMSOL Multiphysics.

By using the *Current Distribution on Edges, BEM* interface, you only need to create a discretization for the edges. The rest of the calculation is based on the boundary element method, which implies that no discretization of the surrounding volume is needed. Creating an edge mesh for this structure will take only a few seconds, but for a full-sized solid structure with a volumetric mesh, much more work and computational resources would be required.

### Utilizing the Capabilities of the Boundary Element Method

When building a model with the *Current Distribution on Edges, BEM* physics interface, you need to select it from the *Physics Interfaces* selection list under the *Electrochemistry* branch.

*A screenshot showing the* Select Physics *window, with access to the Corrosion Module.*

Once you have built your geometry consisting of edges, you can add the *Edge Radius* of the different edges as an attribute in the *BEM Geometrical Properties* node.

*Define the radius of your edge or tube as an* Edge Radius *attribute in the* BEM Geometrical Properties *feature.*

By specifying the edge radius, you may also specify that it should compensate for the tube volume in the electrolyte. In turn, the edges are implemented as edge sources in an infinite electrolyte. However, if there are many edges or tubes and they are close together, they will reduce the effective conductivity of the electrolyte. Thus, if you check the box for compensation for the tube volume, each tube is treated as a solid tube at the cost of two more degrees of freedom at each element.

*Compensation for tube volume, which displaces the electrolyte, may be enabled by the check box in the* Settings *window.*

Now you are ready to specify the potential or current conditions on each edge. By right-clicking on the *Current Distribution on Edges, BEM* node, you’ll find the selection list of your available conditions.

*Define boundary conditions on the edges as usual.*

The default plot that appears after solving includes a *Line* plot along the edges, where the thickness of the tube is taken from the specified tube radius in the physics interface.

*The potential on the edges with a tube plot, showing the virtual dimensions of the tubes.*

In this blog post, we have shown a new way of handling large slender structures with the *Current Distribution on Edges, BEM* interface. By utilizing the boundary element method and specialized elements, we are able to reduce the time spent on meshing and model setup.

### Enhance Your Corrosion Simulation Process with These Resources

- Take a closer look at the Corrosion Module to learn more about the capabilities of COMSOL Multiphysics for cathodic protection and corrosion simulation
- Check out the COMSOL Blog for in-depth stories and demonstrations on both corrosion and electrochemistry topics

## Comments

EH OoiFebruary 19, 2016 1:05 amMay I know if the BEM feature is available only for the Electrochemistry Module?

Thank you.

Bertil NistadFebruary 19, 2016 9:22 amThe Current distribution on edges, BEM interface is only available in the electrochemistry products. (Batteries and Fuel Cells-, Corrosion-, and Electrochemistry Module)

Let me know if you have any additional questions.

Best regards

Bertil

Djedjiga BOUKHLEFMarch 30, 2017 9:03 amI need to distance training for simulation of the cathodic protection for the buried pipes by the methode of protection by galvanic anode.

Cordially.

Caty FaircloughApril 13, 2017 5:09 pmHi,

Thanks for your comment!

For questions related to your modeling, please contact our support team.

Online support center: https://www.comsol.com/support

Email: support@comsol.com

Best,

Caty