An in situ (internal) electric field is used as a dosimetric quantity for human protection from low-frequency electromagnetic fields (lower than 5 MHz) under international safety ...standard/;guidelines. The IEEE standard uses a homogenous elliptical cross section to derive external field strength corresponding to an in situ field strength, while the International Committee on Non-Ionizing Radiation Protection (ICNIRP) guidelines use anatomical models to relate them. In the latter, "the 99th percentile value of the in situ electric field averaged over the cube of its side length of 2 mm" is used to represent the maximum in situ electric field. This metric was introduced to suppress computational artifacts that are inherent when using voxelized anatomical models, in which curved boundaries are discretized with a stair-casing approximation. To suppress the error, a few schemes have been proposed for treating the computational artifacts. In this study, the various schemes to suppress the artifacts are reviewed. Subsequently, a postprocessing method for determining the appropriate maximum in situ field strength is proposed. The performance of the proposed scheme is first verified by comparison with an analytical solution in a multilayered sphere. The method is then applied for different exposure scenarios in anatomically realistic human models where the volume under computation is also considered.
The presence of time-varying electromagnetic fields across a neuron cell may cause changes in its electrical characteristics, most notably, in the action potential dynamics. This phenomenon is ...examined by simulating electrophysiology of a single cortex neuron. Magnetic flux is captured by using a polynomial approximation of a memristor embedded into Hodgkin–Huxley model, equivalent electrical circuit of a neuron cell. Bifurcation analysis is carried out for two different electrical modes associated with the nature of the external neuronal stimulus. Aiming to determine the true influence of the variability of ion channels conductivity, the stochastic sensitivity analysis is undertaken post hoc. Additionally, numerical simulations are enriched with uncertainty quantification, observing values of ion channel conductivity as random variables. The aim of the study is to computationally determine the sensitivity of the action potential dynamics with respect to changes in conductivity of each ion channel so that the future experimental procedures, most often medical treatments, may be adapted to different individuals in various environmental conditions.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Goal: The aim of this paper is to provide a rigorous model and, hence, a more accurate description of the transcranial magnetic stimulation (TMS) induced fields and currents, respectively, by taking ...into account the inductive and capacitive effects, as well as the propagation effects, often being neglected when using quasi-static approximation. Methods: The formulation is based on the surface integral equation (SIE) approach. The model of a lossy homogeneous brain has been derived from the equivalence theorem and using the appropriate boundary conditions for the electric field. The numerical solution of the SIE has been carried out using the method of moments. Results: Numerical results for the induced electric field, electric current density, and the magnetic flux density distribution inside the human brain, presented for three typical TMS coils, are in a good agreement with some previous analysis as well as to the results obtained by analytical approach. Conclusion: The future work should be related to the development of a more detailed geometrical model of the human brain that will take into account complex cortical columnar structures, as well as some additional brain tissues. Significance: To the best of authors knowledge, similar approach in modeling TMS has not been previously reported, albeit integral equation methods are seeing a revival in computational electromagnetics community.
The paper deals with an efficient approach to determine complex power generated by a thin wire antenna above a lossy ground. Once the current distribution along a thin wire is obtained by numerically ...solving the Pocklington integro-differential equation the complex power of the antenna can be obtained by solving the integral over the inner product of tangential component of the electric field and current distribution along the wire stemming from Poynting theorem. Numerical calculation procedure uses the vectors and matrices already constructed within the calculation procedure for the current distribution. Some illustrative numerical results for the active power, reactive power and apparent power of the dipole antenna radiating over a lossy half-space are presented in the paper for different values of input parameters.
This paper is on the use of a hybrid boundary element method/finite element method (BEM/FEM) to determine the induced electric field in spherical human head models exposed to high-frequency plane ...electromagnetic (EM) wave. The geometrically simplified models include the homogeneous one and the non-homogeneous one, featuring compartments such as skin, skull, CSF, and brain. Both models are illuminated by plane EM wave at frequencies including 900 and 1800 MHz, 3500 MHz pertaining to 5G communication systems and also 6000 MHz representing the transition frequency related to the EMF safety standards. The numerical results for the electric field induced in both human head models are presented, while the emphasis is on the electric field along the propagation axis. The novelty of this work is related to the subsequent post-processing of the sampled induced field along the model axis by using two different numerical filtering techniques. It is shown that using the peak detection algorithm, the spline interpolation could be used to estimate the signal envelope. The exponentially decaying nature of the envelope allows to assess the penetration depth of the EM radiation within the biological tissue. Moreover, it is shown that the analytically calculated penetration depth, derived for the unbounded medium, is well reproduced by the numerical computation of EM field.
Recent relevant safety guidelines IEEE-Std C95.1- 2019 and ICNIRP-RF Guidelines 2020 have converged towards 6 GHz as a transition frequency from specific absorption rate (SAR), as basic restriction ...quantity, to absorbed power density (APD). Namely, the penetration of electromagnetic waves into the human tissue rapidly decreases as frequency increases, therefore, tissue heating can be considered as superficial above 6 GHz. However, besides the APD, an alternative internal dosimetric quantity transmitted power density or TPD is sometimes computed since its relation to SAR is more obvious and is easier to obtain. This paper deals with an analytical/numerical approach to determine TPD in planar multi-layered model of the human tissue exposed to the dipole antenna radiation. Analytical approach deals with assumed sinusoidal current distribution, while numerical approach pertains to the determination of current by solving the corresponding Pocklington integro-differential equation via Galerkin-Bubnov Indirect Boundary Element Method. The novelty presented in this paper with respect to previous work is a multilayer geometry whose effects are considered via the corresponding Fresnel plane wave reflection/transmission approximation. Some illustrative results for current distribution, transmitted field, volume power density (VPD) and TPD at various frequencies and distances of the antenna from the interface are given.
The paper deals with an efficient procedure to study human exposure to vertical dipole antenna above a flat lossy half space. The closed form expressions for the corresponding irradiated electrical ...field are obtained assuming the sinusoidal and triangular current distribution along the antenna, respectively. The corresponding integral field expressions are evaluated by means of numerical integration and analytical procedures. The computations have been undertaken in the far field zone for various antenna parameters and compared to rigorous numerical model. Provided the field radiated from the vertical dipole is determined, whole body average Specific Absorption Rate (SARWB) is computed in a simple parallelepiped model of the human body and cylindrical model of the human body, respectively.
The paper deals with the uncertainty quantification of the transient axial current induced along the human body exposed to electromagnetic pulse radiation. The body is modeled as a straight wire ...antenna whose length and radius exhibit random nature. The uncertainty is propagated to the output transient current by means of the stochastic collocation method. The stochastic approach is entirely nonintrusive and serves as a wrapper around the deterministic code. The numerical deterministic model is based on the time domain Hallen integral equation solved by means of the Galerkin-Bubnov indirect boundary element method (GB-IBEM). The stochastic moments, i.e., the mean and the variance of the transient current, are calculated. Confidence margins are obtained for the whole duration of the transient response as well as for the maximal current value. The presented approach enables the estimation of the probability for the induced current to exceed the basic restrictions prescribed by regulatory bodies. The sensitivity analysis of the input parameters indicates to which extent the variation of the input parameter set influences the output values which is particularly interesting for the design of the human equivalent antenna.
Full text
Available for:
FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
International organizations have collaborated to revise standards and guidelines for human protection from exposure to electromagnetic fields. In the frequency range of 6-300 GHz, the permissible ...spatially averaged epithelial/absorbed power density, which is primarily derived from thermal modeling, is considered as the basic restriction. However, for the averaging methods of the epithelial/absorbed power density inside human tissues, only a few groups have presented calculated results obtained using different exposure conditions and numerical methods. Because experimental validation is extremely difficult in this frequency range, this paper presents the first intercomparison study of the calculated epithelial/absorbed power density inside a human body model exposed to different frequency sources ranging from 10-90 GHz. This intercomparison aims to clarify the difference in the calculated results caused by different numerical electromagnetic methods in dosimetry analysis from 11 research groups using planar skin models. To reduce the comparison variances caused by various key parameters, computational conditions (e.g., the antenna type, dimensions, and dielectric properties of the skin models) were unified. The results indicate that the maximum relative standard deviation (RSD) of the peak spatially averaged epithelial/absorbed power densities for one- and three-layer skin models are less than 17.49% and 17.39%, respectively, when using a dipole antenna as the exposure source. For the dipole array antenna, the corresponding maximum RSD increases to 32.49% and 42.55%, respectively. Under the considered exposure scenarios, the RSD in the spatially averaged epithelial/absorbed power densities decrease from 42.55% to 16.7% when the frequency is increased from 10-90 GHz. Furthermore, the deviation from the two equations recommended by the exposure guidelines
Numerical artifacts affect the reliability of computational dosimetry of human exposure to low-frequency electromagnetic fields. In the guidelines of the International Commission of Non-Ionizing ...Radiation Protection, a reduction factor of 3 was considered to take into account numerical uncertainties when determining the limit values for human exposure. However, the rationale for this value is unsure. The IEEE International Committee on Electromagnetic Safety has published a research agenda to resolve numerical uncertainties in low-frequency dosimetry. For this purpose, intercomparison of results computed using different methods by different research groups is important. In previous intercomparison studies for low-frequency exposures, only a few computational methods were used, and the computational scenario was limited to a uniform magnetic field exposure. This study presents an application of various numerical techniques used: different finite-element method (FEM) schemes, method of moments, and boundary-element method (BEM) variants, and, finally, by using a hybrid FEM/BEM approach. As a computational example, the induced electric field in the brain by the coil used in transcranial magnetic stimulation is investigated. Intercomparison of the computational results is presented qualitatively. Some remarks are given for the effectiveness and limitations of application of the various computational methods.