Search results

The mathematical complexity of the B_J(x) Brillouin function makes it unsuitable for most calculations and its application difficult for computer programming in magnetism. Here, its approximation with the tanh function is proposed to ease the mathematical operations for most cases. The approximation works with good accuracy, acceptable in practical calculations. This approximation has already formed the foundation of the “hyperbolic model” in magnetism for the study of hysteretic phenomena. The reversal of the Brillouin function is an important but difficult mathematical problem for practical purposes. Here, a proposal has been put forward for an easy approximation using an analytical expression. This provides a good workable solution for the B_J(x)⁻¹ function dependent on J, the angular momentum quantum number of the material used. The proposed approximation is applicable within the working range of practical applications. The paper aims to discuss these issues.

The multi-variant Brillouin function is closely approximated by the tanh function to ease calculations. Its mathematically unsolved reversed function is approximated by a simple analytical expression with a good working accuracy.

The Brillouin function and its reversal can be approximated for practical users mostly for professionals working in Magnetism.

Most if not all practical problems in Magnetism can be solved within the limitations of the two approximations.

Both proposed functions can ease the mathematical problems faced by researchers and other users in Magnetism.

Ease the frustration of most users working in the field of Magnetism.

The application of the tanh function for replacing the Brillouin function led to the creation of the hyperbolic model of hysteresis. To the author's knowledge, the reverse function was mathematically only solved in 2015 with a vastly complicated mathematics, and is hardly suitable for practical calculations in Magnetism. The proposed simple expression can be very useful for theorists and experimental scientists.

This paper starts with the description of a purely mathematical model of the saturation curve and the hysteresis loop based on the fundamental similarities between the Langevin function the specified T(x) function and the sigmoid shape. The T(x) function which is composed of tangent hyperbolic and linear functions with its free parameters can describe the regular anhysteretic magnetisation curve. Developed from this function the model describes not only the regular hysteresis loop but also the biased and other minor loops like the ones produced by the interrupted and reversed magnetisation process and the open “loops” created by a piecewise monotonic magnetising field input of diminishing amplitude. The remanent magnetism as the function of the interrupted field co‐ordinates is predicted by the model in this mathematical form for the first time. The model presented here is based on the principle that all processes follow the shape of the T(x) function describing the shape of the major hysteresis loop of the ferromagnetic specimen under investigation. The model is also applicable to hysteretic processes in other fields.

Improperly fitted parameters for the Jiles–Atherton (JA) hysteresis model can lead to non-physical hysteresis loops when ferromagnetic materials are simulated. This can be remedied by including a proper physical constraint in the parameter-fitting optimization algorithm. This paper aims to implement the constraint in the meta-heuristic simulated annealing (SA) optimization and Nelder–Mead simplex (NMS) algorithms to find JA model parameters that yield a physical hysteresis loop. The quasi-static B(H)-characteristics of a non-oriented (NO) silicon steel sheet are simulated, using existing measurements from a single sheet tester. Hysteresis loops received from the JA model under modified logistic function and piecewise cubic spline fitted to the average M(H) curve are compared against the measured minor and major hysteresis loops.

A physical constraint takes into account the anhysteretic susceptibility at the origin. This helps in the optimization decision-making, whether to accept or reject randomly generated parameters at a given iteration step. A combination of global and local heuristic optimization methods is used to determine the parameters of the JA hysteresis model. First, the SA method is applied and after that the NMS method is used in the process.

The implementation of a physical constraint improves the robustness of the parameter fitting and leads to more physical hysteresis loops. Modeling the anhysteretic magnetization by a spline fitted to the average of a measured major hysteresis loop provides a significantly better fit with the data than using analytical functions for the purpose. The results show that a modified logistic function can be considered a suitable anhysteretic (analytical) function for the NO silicon steel used in this paper. At high magnitude excitations, the average M(H) curve yields the proper fitting with the measured hysteresis loop. However, the parameters valid for the major hysteresis loop do not produce proper fitting for minor hysteresis loops.

The physical constraint is added in the SA and NMS optimization algorithms. The optimization algorithms are taken from the GNU Scientific Library, which is available from the GNU project. The methods described in this paper can be applied to estimate the physical parameters of the JA hysteresis model, particularly for the unidirectional alternating B(H) characteristics of NO silicon steel.

The purpose of this paper is to investigate an alternative to established hysteresis models.

Different mathematical representations of the magnetic hysteresis are compared and some differences are briefly discussed. After this, the application of the T(x) function is presented and an inductor model is developed. Implementation details of the used transient circuit simulator code are further discussed. From real measurement results, parameters for the model are extracted. The results of the final simulation are finally discussed and compared to measurements.

The T(x) function possesses a fast mathematical formulation with very good accuracy. It is shown that this formulation is very well suited for an implementation in transient circuit simulator codes. Simulation results using the developed model are in very good agreement with measurements.

For the purpose of this paper, only soft magnetic materials were considered. However, literature suggests, that the T(x) function can be extended to hard magnetic materials. Investigations on this topic are considered as future work.

While the mathematical background of the T(x) function is very well presented in the referenced papers, the application in a model of a real device is not very well discussed yet. The presented paper is directly applicable to typical problems in the field of power electronics.

This paper aims to present an analytical way of formulating the vital parameters of an equivalent hysteresis loop of a composite, multi-component magnetic substance. By using the hyperbolic model, the only model, which separates the constituent parts of the composite magnetic materials, an equivalent loop can be composed analytically. So far, it was only possible to superimpose the tanh functions by numerical method. With this transformation, all multi-component composite substances can be treated mathematically as a single-phase material, as in the T(x) model, and include it in mathematical operations. The transformation works with good accuracy for major and minor loops and provides an easy analytical way to arrive to the vital parameters. This also shows an analytical way to the easy solution of some of the difficult problems in magnetism for multi-component ferrous materials, such as Fourier and Laplace transforms, accommodation and energy loss, already solved for the T(x) model.

The mathematical single loop formulation of hysteresis loop of a multi-phase substance shows the way in good approximation of the sum of constituent loops, described by tanh functions. That was so far only possible by numerical methods. By doing so, it becomes equivalent to the T(x) model for mathematical operations.

The described method gives an analytical formulation [identical to the T(x) model] of multi-component hysteresis loops described by hyperbolic model, leading to simple solution of difficult problems in magnetism such as loop reversal.

Although the method is an approximation, its accuracy is good enough for use in magnetic research and practical applications in industries engaged in application of magnetic materials.

The hyperbolic model is the only one which separates the magnetic substance, used in practice, to constituent components by describing its multi-component state. Superimposing the components was only possible so far by numerical means. The transformation shown is an analytical approximation applicable in mathematical calculations. The transformation described here enables the user to apply all rules applicable to the T(x) model.

This study equally helps researchers and practical users of the hyperbolic model.

This novel analytical approach to the problem provides an acceptable mathematical solution for practical problems in research and manufacturing. It shows a way to solutions of many difficult problems in magnetism.

The paper sets out to formulate the intermolecular forces leading to Barkhausen instability. In the approach the known concept of effective field is used within the framework of the T(x) model. The aim is to provide a mathematical tool to theoreticians and applied scientists in magnetism that is easier to use than those of other models. At the same time to demonstrate the easy applicability of the T(x) model to hysteretic phenomena.

With the combination of the effective and the external field the model is applied to hysteresis loops as well as to the anhysteretic state showing in both cases the local development of unstable conditions at beyond a critical point, leading to local hysteresis loops.

The paper formulates the critical conditions for the hysteretic and the anhysteretic process and calculates the susceptibility as the functions of magnetisation and the applied field.

Experimental verification will be required to prove the applicability to the various magnetic materials and to the accuracy of the model.

The paper provides an easy mathematical and visual method to show the conditions before and after the Barkhausen instability sets in during the magnetisation process.

The paper provides an easy mathematical tool for theoreticians and experimental scientists with a visual presentation of processes leading to Barkhausen instability and magnetic behaviour beyond that by using the T(x) model.

The purpose of this paper is the modelling and estimation of inrush currents while energising power devices under no load conditions. An analytical representation of the nonlinear B-H curve serves for considering the hysteresis behaviour in the numerical model.

The model is implemented into a standard finite element formulation to compute transient problems.

Inrush currents behave like faults in power distribution facilities. Its prior estimation helps to distinguish between operating conditions and faults.

The magnetic cores may become extremely magnetised. At such high material saturations, the material characteristics are not measurable accurately. Hence, the results depend on the extrapolation of the B-H curves.

The use of first-order reversal curves within the major hysteresis loops helps in a convenient way to estimate peak and shape of the inrush currents.

The aim of the paper is to present a simple approach to modelling minor hysteresis loops in grain‐oriented steel sheets under quasi‐static and dynamic conditions. The hysteresis phenomenon is described with a recently developed hybrid model, which combines ideas inherent in the product Preisach model and the Jiles‐Atherton description. The dynamic effects due to eddy currents are taken into account in the description using a lagged response with respect to the input.

It is assumed that some model parameters might be dependent on the level of relative magnetization within the material. Their dependencies could be given as power laws. The values of scaling coefficients in power laws are determined.

A satisfactory agreement of experimental and modelled quasi‐static and dynamic hysteresis loops is obtained.

The present study provides a starting point for further verification of the approach for other classes of soft magnetic materials, which could be described with the developed model. At present, the approach to model minor loops by the update of model parameters is verified for the B‐sine excitation case.

The “branch‐and‐bound” optimization algorithm is a useful tool for recovery of the values of both model parameters and scaling coefficients as well.

The recently developed hybrid description of hysteresis phenomenon can be successfully extended to take into account symmetric minor loops. The developed approach could be a framework to develop a comprehensive description of magnetization phenomena in the future.

The purpose of this paper is to demonstrate an effective and faster numerical solution for nonlinear‐coupled differential equations describing fiber amplifiers which have no explicit solution. MATLAB boundary value problem (BVP) solver of bvp6c function is addressed for the solution.

Coding method with the bvp6c is introduced, signal evolution, threshold calculation method is introduced, gain and noise figure are plotted and superiority of the bvp6c solver is compared with the Newton‐Raphson method.

bvp6c function appears to be an effective tool for the solution fiber amplifier equations and can be used for different pump configurations of BFAs and RFAs. The excellent agreement between the proposed and reported results shows the reliability of the proposed threshold power calculation method.

The paper eases the work of the fiber optic research community, who suffer from two point BVPs. Moreover, the stiffness of the signal evolution which is faced with high pump powers and/or long fiber lengths can be solved with continuation. This superiority of the solver can be used to overcome any stiff changes of the signals for the future studies.

The main outcome of this paper is the numerically calculation of the threshold values of fiber amplifiers without the necessity of the experiment. The robustness improvement of the solution is that the solver is able to solve the equations even with the poor guess values and the solution can be obtained without the necessity of analytical Jacobian matrix.

MATLAB bvp6c solver has proven to be effective for the numerical solution of nonlinear‐coupled intensity differential equations describing fiber amplifiers with two‐point boundary values. Beside the signal evolution, thresholds of Brillouin and Raman fiber amplifiers can also be calculated by using the proposed solver. This is a notable and promising improvement of the paper, at least from a fiber optic amplifier designer point of view.

The purpose of this paper is to present the method of lines (MOL) solution of the stimulated Brillouin scattering (SBS) equations (a system of three first-order hyperbolic partial differential equations (PDEs)), describing the three-wave interaction resulting from a coupling between light and acoustic waves. The system has complex numbers and boundary values.

System of three first-order hyperbolic PDEs are first transformed and then spatially discretized. Superbee flux limiter is proposed to offset numerical damping and dispersion, brought on by the low order approximation of spatial derivatives in the PDEs. In order to increase computational efficiency, the structured structure of the PDE Jacobian matrix is identified and a sparse integration algorithm option of the ordinary differential equation (ODE) solvers is used. The flux limiter based on higher order approximations eliminates numerical oscillation. Examples are presented, and the performance of the Matlab ODE solvers is evaluated by comparison.

This type of solution provides a rapid means of investigating SBS as a tool in fiber optic sensing.

To the best of the authors' knowledge, MOL solution is proposed for the first time for the modeling of three-wave interaction in a SBS-based fiber optic sensor.