A New Approach for Solving Weight Functions of Electromagnetic Flowmeters Using Resistive Network Modeling

Abstract

The contribution of the flow signal is generally addressed by the weight function in the researches of the electromagnetic flowmeters, and various mathematical technologies were concentrated on the methodologies for solving the value of the weight function. However, it is still difficult to avoid the abstruse mathematical theories and the complex calculation when the solution domain is irregular in shape. This paper treats the problem within the intuitive physical perspective, and the approach, in which the proportion of the current is considered as the substitute for the weight function with the hypothetic current excitation source, is presented. A simple mathematical modeling of the current is built by means of the resistive network without the redundant assumption, and the strict mathematical derivation for the conventional asymmetric flow in the circular flowmeter is made to verify the feasibility and the correctness of the approach. The distributions of the weight function in various situations are obtained with the simulation employed, using the resistive network modeling, and the advantages of the approach are discussed.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Recommended by Oluwole Daniel Makinde

School of Mechatronics Engineering and Automation, Shanghai University, Shanghai 200072, China

Received 22 February 2013; Accepted 17 June 2013

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

The study on the weight function is regarded as the general and essential focus in the field of the electromagnetic flowmeter (EMF). In 1954, the first analytical study for the EMF was performed by Shercliff [1] based on the two-dimensional rectilinear flow and the uniformly distributed magnetic field. The microscopic characteristics of the induction electromotive force were demonstrated with the weight function method he proposed which is now widely used in analyzing the EMF. In 1970, Bevir [2] made further progress on the weight function theory with the virtual current introduced and extended it to the three-dimensional case. In 1971, Smyth [3] investigated the weight functions for the circular and rectangular channel EMFs by means of Green's function and conformal mapping. In 1983, O'Sullivan and Wyatt [4] researched the weight functions for electrodes with different geometries, and various maps of the distribution were obtained. In 1986, Hemp and Versteeg [5] solved the virtual current potential and the magnetic scalar potential by the expansion of Double Fourier series and simplified the calculation of the weight function significantly. In 1989, Zhang [6] proposed the alternating iteration method for solving the Laplace equation and applied it to achieve the semianalytical solution of the weight function for the EMF, on condition that the mixed boundary exists. In 2009, Hu et al. [7] divided the complex geometry by different auxiliary surfaces to obtain the corresponding analytical solutions of the weight function with the method of separation of variables.

Much effort has been made to solve the weight function for different situations by virtue of the mathematical skills in previous studies, based on the theory of the virtual current. The distribution of the weight function is especially the same as that of the virtual current if the magnetic field is uniform. However, the physical meaning of the weight function was less mentioned for the purpose of the calculation, by which the simple and direct approach to solve the weight function can be obtained.

In this paper, the theoretical background is briefly introduced, and the physical meaning of the weight function is proposed with the innovative modeling of the resistive network built, which results in the analytical solution of the weight function for the conventional asymmetric flow in the circular flowmeter. Besides, the simulations of the weight function are also involved to verify the correctness of the methodology addressed, and the conclusion is presented.

2. Theoretical Background

As is known, the basic measurement principle of the EMF is Faraday's Law, and the potential between two electrodes fulfills the integral equation which can be expressed as follows: [figure omitted; refer to PDF] where W [arrow right] is given by [figure omitted; refer to PDF] Here, W [arrow right] is the weight vector which is the extension of the weight function termed by Bevir and is considered as the general index for the distribution. B [arrow right] is the magnetic flux vector, and j [arrow right] is the current vector, which characterizes the current density on condition that the unit current passes from the positive electrode to the negative electrode. However, the current does not exist in the flow channel with the practical application of the EMF, and thus the current is called "virtual current." With the uniform magnetic field assumed, the distribution of the weight function is equal to that of the virtual current, and its value can be calculated by the partial differential equations.

Generally, there are no electric power and magnetic source inside the sensor, and the electrical potential [straight phi] and magnetic potential F within the flow area satisfy the Laplace' equations, [figure omitted; refer to PDF] and the virtual current density j [arrow right] and magnetic flux vector B [arrow right] can be further denoted by [figure omitted; refer to PDF] According to (2), (3), and (4), the weight function is obtained by the following equation: [figure omitted; refer to PDF] where _{X e} [arrow right] , _{Y e} [arrow right] , and _{Z e} [arrow right] are the unit vectors, respectively, in cartesian coordinate system.

From the previous equations, the solution of the weight function is mainly determined by the distribution of the virtual current which can be presented by the potential function. In previous researches, the methodologies for solving the direct analytical problems of the Laplace's equations to achieve the weight function were focused on, especially in the case of the various solution domains.

Different from the theory of the virtual current, this paper is concerned with the physical understanding of the weight function, which simply presents the proportion of the real current distribution. With the current excitation source assumed, the real current flows from the positive electrode to the negative electrode in the case of the static fluid, and the new modeling of the resistive network is established as the analogue of the flow, through which the weight function can be solved with an explicit physical explanation.

3. Modeling and Derivation

The EMF is generally assumed as the cylinder, and the electrodes are symmetrically placed in the circumference. The voltage signal is induced between the electrodes when the conductive fluid flows through the uniform magnetic field. The field theory analysis is made conventionally to build the virtual current model which depends on the mathematical skill of the partial differential equations other than the explicit physical meaning.

Based on the previous situation, a new approach for modeling is presented. The flow is regarded as the resistive network, and the current and potential in the fluid domain are solved with Kirchhoff's Current Law (KCL) employed.

Two different kinds of the grid system are considered in this model, the cartesian grid and the polar grid. The former grid shown in Figure 1 is established with several vertical and horizontal resistances.

Figure 1: The cartesian grid of the resistive network.

[figure omitted; refer to PDF]

According to the Kirchhoff's Current Law (KCL), in which the sum of the nodal current is zero, the mathematical expression is directly obtained as [figure omitted; refer to PDF] where ρ is the flow resistivity, and Δ x and Δ y are the unit step of the grid along the x axis and the y axis, respectively. Here, let Δ x = Δ y = Δ , with the equal interval step introduced: [figure omitted; refer to PDF]

The left part of the equation above is found to be the Laplacian operator, and (7) is equal to the equation ^{∇ 2} U ( x , y ) = 0 , through which the feasibility and the correctness of the modeling using the resistive network can be verified. Nevertheless, it is difficult to approach the boundary with the cartesian grid utilized because of the pipe section in shape, and the latter grid, that is the polar grid, is then used to make up for the deficiency, demonstrated in Figure 2.

Figure 2: The polar grid of the resistive network.

[figure omitted; refer to PDF]

Similarly, the equation is established on the basis of the Kirchhoff's Current Law (KCL): [figure omitted; refer to PDF]

The grid is generated minutely, such that [figure omitted; refer to PDF] The partial differential equation is obtained with the combination of (8) and (9) [figure omitted; refer to PDF]

With the method of separation of variables, the equation has the particular solution which has the form as follows: [figure omitted; refer to PDF]

Substituting (11) into (10), the equation can be arranged as follows: [figure omitted; refer to PDF] where n is the constant based on the structure of the equation.

Equation (12) is considered separately. The first section of it is the Euler equation, while the second section of it is the second-order linear homogeneous equation of constant coefficients. Then the analytical solution of the potential is [figure omitted; refer to PDF] where _{K 1} and _{K 2} are denoted as the coefficient terms related the constant n .

The boundary condition is used to determine the coefficients. According to Kirchhoff's Current Law (KCL), the expression of the pipe boundary is obtained: [figure omitted; refer to PDF]

The current represented by the first two terms of the equation does not exist because of the insulating material expect the electrodes in the senior, and the continuous form of (14) is [figure omitted; refer to PDF] where a presents the pipe radius.

For the electrodes, it is generally regarded as the point charges, and the expression at the electrodes is conventionally set as a constant by means of the mathematical skill without the physical meaning. In this paper, with the assumed current excitation source, the expression at the electrodes can be obtained: [figure omitted; refer to PDF]

Combined with (13), (15), and (16), the equation is arranged as follows: [figure omitted; refer to PDF]

Multiplying sin n θ on the both sides of the equation, with the definite integration applied, (17) is arranged, and _{K 1} can be solved: [figure omitted; refer to PDF]

It is noticed that _{K 1} is not equal to zero, so n is the odd number. Assume that n = 2 m - 1 , _{K 1} is transformed as [figure omitted; refer to PDF]

Similarly, multiplying cos n θ on both sides of the equation, with the definite integration employed, (17) can also be arranged, and _{K 2} is determined: [figure omitted; refer to PDF]

By means of (13), (19), and (22), the analytical solution of the potential is [figure omitted; refer to PDF]

For simplification, ( 1 + 4 ( 2 m - 1 ^{) 2} + 1 ) / 2 is equivalent to 2 m - 1 in the infinite series, and then the potential is finally obtained as follows: [figure omitted; refer to PDF]

Based on the Kirchhoff's Current Law (KCL), the radial current and the circumferential current of the arbitrary point in the flow domain are addressed as [figure omitted; refer to PDF]

The y axis is selected as the unified direction of the current; thus, the actual current is given by the sum of the component along the y axis, including the radial current and the circumferential current, shown as follows: [figure omitted; refer to PDF]

According to the definition, the weight function not only reflects the degree of the contribution of the flow signal generated from the motion of the flow through the magnetic field, but also describes the attenuation coefficient of the induced voltage caused by the geometric position of the points in the effective domain. From the physical meaning, the weight function is actually the distribution proportion. Thus, it is extended in this paper and regarded as the proportion of the current in the total current. Consider the following: [figure omitted; refer to PDF]

Let N = m - 1 ; then (27) is arranged as [figure omitted; refer to PDF]

On basis of the identities, [figure omitted; refer to PDF] the analytical expression of the weight function is solved as [figure omitted; refer to PDF] where K is the constant.

The normalization processing of the weight function is then made by dividing _{W ^} ( 0,0 ) , where [figure omitted; refer to PDF] Thus, the normalized expression of the weight function is finally obtained: [figure omitted; refer to PDF]

With the transformation to the cartesian coordinate system, the solution is [figure omitted; refer to PDF]

4. Results and Discussion

Depending on (33), the solution of the weight function is exactly the same as the classical one, using the approach of the resistive network modeling, which validates the correctness of the scheme in this paper. From the simulation of the current with the modeling employed, it is demonstrated that the proportion of the current in the total current greatly reflects the weight function and it can be the analysis object for the research on the weight function with its modeling presented. In Figures 3(a) to 3(f), the weight functions for both full pipe and partially filled pipe are numerically identical to the proportion of the current in the total current from the simulation results.

(a)-(f) The distributions of the weight function [8, 9] and the current simulated using the proposed approach.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

(e) [figure omitted; refer to PDF]

(f) [figure omitted; refer to PDF]

In the previous studies on the classical situation, it was generally assumed that the electrodes were two point charges with different polarities, which caused the Neumann boundary condition, and the assumption of the conventional modeling previous brought the complexity to establish the Neumann equation in which the value of the partial derivative component determined by some mathematical skills [1, 6] is not explicit. With the characterization of the weight function by using the proportion of the current introduced in this paper, the mathematical constraint of the natural boundary condition is weakened, and the convenience for the calculation is provided.

With the methodology addressed, the weight functions for the EMFs which possess various kinds of electrodes can also be solved and the problems of the boundary condition are transformed to the difficulties of the modeling of the electrodes contacting the fluid, and thus, the reasonable selection of the model for boundary electrodes is essential by using the resistive network in the research on the weight function in the field of the electromagnetic flowmeter (see Figure 4).

(a)-(b) The solutions of the weight function for EMFs with long-arc electrodes using resistive network modeling.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

5. Conclusions

The approach for solving the weight function is developed with the strict derivation of the mathematical expressions. The principle of the methodology is simple and effective, with the intuitive physical meaning of the weight function introduced, and it also possesses practical value to some extent in the field of the analysis on EMF. However, the boundary condition is various for different kinds of EMF, and the uniform modeling of the electrodes with complex geometric structure is to be further studied with this approach proposed.

Acknowledgments

The authors are grateful to the financial support of the Specialized Research Fund for the Doctoral Program of Higher Education of China (no. 20123108110011), titled "Research on the Weight Function and the Measurement of the EMF in Partially Filled Pipes".

References

[1] J. A. Shercliff The Theory of Electromagnetic Flow-Measurement , Cambridge university press, London, UK, 1962.

[2] M. K. Bevir, "The theory of induced voltage electromagnetic flowmeter," Journal of Fluid Mechanics , vol. 43, pp. 577-590, 1970.

[3] C. C. Smyth, "Derivation of weight functions for the circular and rectangular channel magnetic flowmeters, by means of Green's theorem and conformal mapping," Journal of Physics E , vol. 4, no. 1, pp. 29-34, 1971.

[4] V. T. O'Sullivan, D. G. Wyatt, "Computation of electromagnetic flowmeter characteristics from magnetic field data. III. Rectilinear weight functions," Journal of Physics D , vol. 16, no. 8, pp. 1461-1476, 1983.

[5] J. Hemp, H. K. Versteeg, "Prediction of electromagnetic flowmeter characteristics," Journal of Physics D , vol. 19, no. 8, pp. 1459-1476, 1986.

[6] Z. X. Zhang, "A method for solving Laplace's equation with mixed boundary condition in electro magnetic flow meters," Journal of Physics D , vol. 20, pp. 573-576, 1989.

[7] L. Hu, J. Zou, X. Fu, H. Y. Yang, X. D. Ruan, C. Y. Wang, "Divisionally analytical solutions of Laplace's equations for dry calibration of electromagnetic velocity probes," Applied Mathematical Modelling , vol. 33, no. 7, pp. 3130-3150, 2009.

[8] J. A. Shercliff, "Relation between the velocity profile and the sensitivity of electromagnetic flowmeters," Journal of Applied Physics , vol. 25, no. 6, pp. 817-818, 1954.

[9] X.-Z. Zhang, "The virtual current of an electromagnetic flow meter in partially filled pipes," Measurement Science and Technology , vol. 9, no. 11, pp. 1852-1855, 1998.