Numerical analysis consisted largely of an analysis of the change in magnetic flux density due to the shape of the sensor’s core and an analysis evaluating the effect of collecting ferrous particles in the flow. Through improvement of the performance of the sensor, it is intended to numerically show whether the sensitivity is improved by improving the collection effect of the sensor in the flow. Because this sensor measures the amount of ferrous particles after attaching ferrous particles to the sensor with a permanent magnet, the collection effect of the sensor is important.
2.1. Analysis of Magnetic Flux Density for Changes in Sensor’s Internal Design Paramters
First of all, a numerical model was developed to analyze the performance of the existing ferrous particle sensor.
Figure 2 shows the numerical model of the ferrous particle sensor with permanent magnetic. It is hard to obtain exact design parameters, including the sensor’s materials. Thus, primary design parameters were roughly estimated by dismantling the senor. These approximations were used in the analysis. Additionally, the magnetic core was made of low-carbon steel M-50. The magnetic flux density (B)- magnetic field intensity (H) curve was derived from data provided in the analysis program, as shown in
Figure 3. The number of turns of the coil was 300 and the applied current was 0.1 A.
Figure 4 shows the mesh used in numerical model of the sensor. In the analysis, triangle and quad meshes were mixed. The total number of meshes was 11,605. Meshes were densely applied around the sensor to increase the accuracy of the analysis.
Figure 5 shows four sensor models in which the shape of the core onside the sensor is changed. A-model is similar to the existing sensor model. In B-model, the shape of the core was obliquely inclined to concentrate the magnet flux density in the center. In C-model, the magnetic flux density was reduced to the side of the sensor to focus the magnetic flux density on the upper part compared to the B-model. In D-model, the area of the core below the coil was removed and the shape was changed to further improve the magnetic flux density on the upper part of the sensor. Geometries for the four models of the sensor are shown in
Table 1.
This study utilized a multi-physics analysis method. COMSOL 6.0, a commercial multi-physics software was used for numerical calculations. The sensor’s magnetic field was calculated based on Maxwell’s equations. Formulations used for calculating the magnetic field are shown in Equations (1) - (3):
where
H,
J,
B, and
D are magnetic field intensity [A/m], current density [A/m
2], magnetic flux density [T], and magnetic vector potential [A], respectively. The force on a charge σ [C] moving with velocity
v [m/s] in the presence of an electric and magnetic field
E [V/m],
B, is called the Lorentz force and
Je [A/m
2] is an externally generated current density. In this program, the AC/DC module contains the electromagnetic field interface model, which calculates the magnetic flux of the ferrous particle sensor. For cases where currents and electromagnetic field vary slowly, the induced displacement current can be ignored. This assumption, referred to as the quasistatic approximation, is widely used in modeling low-frequencies electromagnetic where the dimensions of the structure is small compared to the wave length [
41].
Magnetic flux density distributions for the four sensor models are shown in
Figure 6. B-model had a higher maximum magnetic flux density than A-model. However, some magnetic flux density was generated at the side part (red dotted area) of the sensor. Thus, the sensitivity of the sensor was not greatly improved. Compared to B-model, C-model removed the area of the core below the coil. As a result, the magnetic flux density was concentrated in the upper part and less magnetic flux density was generated in the side part. Therefore, D-model improved the sensitivity of the sensor by changing the shape of core. The D-model increased the maximum magnetic flux density by about 210% compared to the conventional sensor type (A-model) and improved the sensitivity by lowering the magnetic flux density of the side part.
2.2. Evaluating Sensitivity of the Sensor in the Flow Field
In
Section 2.1, magnetic flux density of the sensor was improved. Although it is essential to fabricate an actual sensor and conduct an experiment for evaluating the performance, a method for verifying whether the sensitivity of the sensor has been improved is suggested using a numerical method. Evaluating the sensitivity of the sensor numerically is economical in terms of cost and time. However, verification through experiments for test devices and lubrication systems is absolutely necessary. This study focused on how to evaluate sensitivity of the sensor with a numerical method. Among the four models of the sensor, only the A-model and D-model were evaluated for sensitivity numerically. This is because A-model is a case in which the existing sensor is described and D-model is a case in which the magnetic flux density of the sensor is greatly improved.
The numerical analysis employed Navier-Stokes equations, the electromagnetic field interface model, and particle tracing module. The AC/DC module contains the electromagnetic field interface model (equations (1)-(3)), which calculates the magnetic flux of the ferrous particle sensor. Navier-Stokes equations for rotating domains are shown below:
where, Ω,
I,
τ, and
F mean angular velocity [1/s], identity matrix, shear stress [N/m
2], volume force [N/m
3], respectively.
The particle tracing module is utilized to compute individual particles’ paths by solving their equations of motion over time, which allows for evaluation of discrete trajectories:
where
mp,
v1,
FD, and
Fext are particle mass [kg], velocity vector of the particle [m/s], drag force [N], and magnetophoretic force [N], respectively. In addition,
τp,
dp,
rp,
ρp,
μ,
μ0,
μr, and
K mean particle velocity response time [s], particle diameter [m], particle radius [m], particle density [kg/m
3], viscosity of fluid [Pa∙s], vacuum permeability [kg∙m∙s/A
2], relative permeability, and nondimentional parameter, respectively.
The motion of particles in a fluid follows Newton’s second law, which states that the net force on an object is equal to the time derivative of its linear momentum in an inertial reference frame as shown in equation (6).
This analysis was conducted using a planar symmetric model as shown in
Figure 7. Various meshes such as tetrahedral, pyramid, prism, and hexahedral were used. The total number of elements was 872,406 in A-model and 1,158,494 in D-model. Moreover, a dense mesh was applied to around the sensor to ensure analytical accuracy. The size of the flow channel was 140 mm (x) × 80 mm (y) × 80 mm (z). The flow was laminar. Working conditions for numerical calculations are shown in
Table 4. Particles used in the analysis were spherical and the material was iron with a density of 8030 kg/m
3. During the initial five seconds of calculation, 1,500 particles were injected from the particle injection area (blue shade area) at intervals of 0.5 seconds. The total number of injected particles was 15,000. The density and viscosity of the lubricant used in the analysis were 870 kg/m
3 and 0.04 Pa∙s, respectively.
In this analysis, not only trajectories of the particles, but also the number of particles collected by the sensor were evaluated. The number of particles collected by the sensor was evaluated while the fluid velocity changed from 0.002 m/s to 0.1 m/s. That is, the collecting effect of the sensor with variation of the fluid velocity was evaluated. Since this sensor uses a permanent magnet to collect ferrous particles and then measures the amount of ferrous wear particles through a change in magnetic field, the sensor’s collecting effect is the most important factor of the sensor’s sensitivity. Therefore, the sensitivity of the sensor was evaluated by determining how many ferrous particles adhered to the sensor under several flow conditions.
Figure 8 shows particle trajectories with time in A-model and D-model when the velocity of fluid is 0.002 m/s. A total of 15,000 ferrous particles flowed from inlet side to the outlet, and among them, some particles started to attach to the top of the sensor from about 30 seconds when A-model was applied. In the case of D-model, some particles began to adhere to the top of the sensor at a time similar to that of A-model. However, the number of particles attached to the sensor of the D-model was greater than that attached to the sensor of the A-model.
Figure 9 shows collected particles on the top of the sensor in A-model and D-model when the fluid velocity is 0.002 m/s. As shown in
Figure 9, the sensitivity of the sensor was evaluated by the number of particles collected on the top of the sensor. In this
Figure 9, the color of particles indicates the magnitude of the magnetophoretic force. Particles collected in the central part of the top surface of the sensor are subjected to a large magnetophoretic force.
Compared to the existing sensor model (A-model), the improved sensor (D-model) investigated how the particle collecting effect appeared according to the change of fluid velocity as shown in
Figure 10. When the fluid velocity was 0.002 m/s, the number of particles attached to the sensor was 3313 in the A-model and 3470 in the D-model. That is, the case of D-model increased about 4.7% compared to the case of A-model in terms of the number of attached particles. When the fluid velocity increased to 0.02 m/s and 0.04 m/s, results of D-model increased about 9.2% and 44%, respectively, compared to those of A-model in terms of the number of attached particles. When the fluid velocity was 0.1 m/s, the D-model had 22 particles attached to the sensor, but the A-model had no particle attached to the sensor. It is difficult to attach particles to the sensor under conditions where the velocity of the fluid is higher. This is because the inertia force of the fluid increases as the velocity of the fluid increases. This result confirms that the sensor of the improved model shows a distinct improvement in sensitivity in a situation where the fluid velocity increases.
Figure 11 shows magnetic flux density and magnetic force lines around the sensor in a fluid field where the fluid velocity is 0.1 m/s. The maximum magnetic flux density of A-model was 0.436 T. However, the maximum magnetic flux density of D-model was 0.913 T, which was about 209 % higher than that of A-model. In addition, it can be seen that it is more advantage to collect particles in the sensor because it is formed toward the upper end of the core in the D-model than in the A-model through lines of magnetic force around the core of the sensor. Therefore, a method for evaluating the sensitivity of the ferrous particle sensor with permanent magnet using a numerical analysis based on multi-physics was proposed. Through this method, it was shown that the sensitivity of the design-changed model was improved in the flow field.