Total Current Research Articles

Identification of hypoxene metabolites in urine samples using gas chromatography–tandem mass spectrometry for anti-doping control

Objectives. Hypoxen is a drug which possesses antioxidant and antihypoxic effects. It achieves this by increasing the utilization of oxygen by mitochondria, intensifying oxidative phosphorylation, and as a result, improving tissue respiration. Athletes take it during prolonged exercise, in order to increase efficiency and reduce physical overwork. Since 2023, the World Anti-Doping Agency has included the drug in the monitoring program, in the belief that it can be used to gain a competitive advantage. It is thus a candidate for inclusion in the Prohibited List as a potential regulator of the human metabolism. Currently, there are no studies or scientific publications focusing on the identification of hypoxene in biofluids for the purpose of anti-doping control. The aim of this study is to determine the possible metabolites of the drug and their chromato-mass spectrometric characteristics in urine samples using gas chromatography–tandem mass spectrometry (GC–MS/MS) for doping control screening purposes.Methods. Sample preparation of urine samples was carried out using enzymatic hydrolysis, liquid–liquid extraction and derivatization. The GC–MS/MS method was used for analysis. Screening of hypoxene metabolites was carried out in the mode of total ion current after fragmentation of selected parent ions.Results. Three specific metabolites of hypoxene (m/z 342, 300, and 346, including trimethylsilyl derivatives) were identified in urine samples of volunteers (n = 3). They can act as markers for taking the target antihypoxant, and their possible structural formulas are given. The excretion curves of two metabolites with an m/z of 300 and 346 respectively in urine were studied. The maximum concentration is reached after 8–14 and 1.5–6 h, respectively. It was established, that these metabolites are reliably identified in urine 90 h or more after a single dose of the drug.Conclusions. Possible structures of hypoxene metabolites in urine samples from volunteers were determined for the first time and their chromato-mass spectrometric characteristics were established. The approach developed in this study can be used for screening analysis of hypoxene for the purpose of anti-doping control.

Chemical Bath Deposition of NiFe Alloy Anode for Efficient Alkaline Water Electrolyzer Integration.

Green hydrogen production from water splitting is a feasible way for intermittent renewable energy storage and utilization, where the exploration and scale-up preparation of high-performance anodic oxygen evolution electrocatalysts are critical prerequisites for its industrial-level applications. Herein, a chemical bath deposition of FeNi3 intermetallic alloys onto Ni mesh support is performed, which delivers a current density of 0.62 Acm-2 at 1.72 V versus reversible hydrogen electrode for alkaline water oxidation in 1 m KOH and an excellent electrolysis stability at 0.2 A cm-2 for over 300 h. Moreover, via 3D computational fluid dynamics simulation and flow field optimization, a homogeneous deposition of ≈5400 cm2 NiFe anode is demonstrated within 4 min using the developed flow bath reactor. Once integrating the as-prepared NiFe anodes into alkaline electrolyzer stack, the voltage variation between each unit cell is below 40 mV at a total operation current of 71 A, or ca. current density of 0.2 A cm-2, confirming the uniformity of this batch synthesis protocol and its great potential for industrial alkaline water electrolysis.

Numerical simulation of the main stage of a lightning

We present a numerical model of the main stage of a lightning discharge. Within the framework of the developed model, evolution of parameters of the current channel upon the return stroke (the lightning main stage) is described by the system of equations governing conservation of mass, momentum, total energy, along with the transmission-line equations for determining the electric potential and the total current in each channel cross section. The main characteristics of lightning at the stage of the return stroke detectable experimentally, such as gas heating in the channel to temperatures in the range of 10–40 kK, the fundamental possibility of propagation of the potential-gradient wave at a speed varying from several hundredth to several tenths of the speed of light, and the possibility of the return-stroke wave propagating a relatively long distance without substantial attenuation, are demonstrated numerically. The conclusion that the developed physical and numerical model of the lightning discharge describes physical processes that occur under real conditions qualitatively correctly can be drawn based on the results on simulation of lightning discharges of various intensity.

The electromagnetic field of a circular parallel plate capacitor excited by an AC current and AC voltage supply

Abstract This study investigates the electromagnetic field (EMF) distribution of an ideal circular parallel plate capacitor excited by a time-harmonic power source. Considering the lead wire and capacitor as a charged whole, we formulate the boundary value problems of the Helmholtz equation for the EMF in the lead wire space and capacitor space, respectively. First, we solve for the EMF generated by the total current in the lead wire space of an AC current. Following this, we solve for the EMF boundary value problem of a capacitor filled with linear, uniform, isotropic, and non-magnetic lossy dielectric under an AC current excitation, using the continuity of the total current as a basis. Second, the EMF distributions in the capacitor space and the lead wire space under an AC voltage excitation are provided, following the principles of the generalized Helmholtz theorem. Third, the EMF distribution is discussed when the electromagnetic ``standing-wave phenomenon" occurs in the ideal dielectric capacitor, and identify the ``resonance phenomenon" and the ``current-stopping phenomenon" of the capacitor's EMF when excited by an AC current and AC voltage, respectively. We also present the corresponding ``resonance frequency" and ``current-stopping frequency". Finally, we explore the quasi-stable solution of the capacitor's EMF under low-frequency condition and the static solutions of the electric and magnetic fields under DC excitation and static voltage excitation. Our findings suggest that the existing formula for the capacitor's EMF approximates our analytical solution under quasi-stable condition.

Control-Oriented Free-Boundary Equilibrium Solver for Tokamaks

A free-boundary equilibrium solver for an axisymmetric tokamak geometry was developed based on the finite difference method and Picard iteration in a rectangular computational area. The solver can run either in forward mode, where external coil currents are prescribed until the converged magnetic flux function ψ(R,Z) map is achieved, or in inverse mode, where the desired plasma boundary, with or without an X-point, is prescribed to determine the required coil currents. The equilibrium solutions are made consistent with prescribed plasma parameters, such as the total plasma current, poloidal beta, or safety factor at a specified flux surface. To verify the mathematical correctness and accuracy of the solver, the solution obtained using this numerical solver was compared with that from an analytic fixed-boundary equilibrium solver based on the EAST geometry. Additionally, the proposed solver was benchmarked against another numerical solver based on the finite-element and Newton-iteration methods in a triangular-based mesh. Finally, the proposed solver was compared with equilibrium reconstruction results in DIII-D experiments.

Conservation of the Charge in Signal from Drift Tube Ion Mobility Spectrometers.

Quantitative information obtained from drift tube detectors used in ion mobility spectrometry is contained in the area of peaks forming the drift time spectrum. The area of all peaks corresponds to the total charge of ions entering the drift section of the spectrometer. It was found that this charge is not conserved when the ion composition changes. This work is devoted to studying the causes of this phenomenon. Experimental research consisted of recording drift time spectra for 2-pentanone and n-heptanone, at various analyte concentrations and different opening times of the shutter grid. Measurements of the total ion current were also performed in static mode with an open grid. The research results indicated that the reasons for the lack of ion charge conservation in the drift time spectrum are ion recombination, mutual repulsion, and mobility-dependent transmission of ions through the shutter grid. The explanation of the relationships obtained experimentally was based on a simple theoretical model, which considered the phenomenon of ion transport along the reaction section and the penetration of ions through the shutter. The developed model provides a good description of the measurement results and allows the estimation of ion currents and ion concentrations in the reaction section upstream of the grid. This information is important for proper quantitative analysis as well as when the ion mobility spectrometer is used in quantitative studies of chemical ionization processes.

Scalable Electro‐Biosynthesis of Ectoine from Greenhouse Gases

Converting greenhouse gases into valuable products has become a promising approach for achieving a carbon‐neutral economy and sustainable development. However, the conversion efficiency depends on the energy yield of the substrate. In this study, we developed an electro‐biocatalytic system by integrating electrochemical and microbial processes to upcycle CO2 into a valuable product (ectoine) using renewable energy. This system initiates the electrocatalytic reduction of CO2 to methane, an energy‐dense molecule, which then serves as an electrofuel to energize the growth of an engineered methanotrophic cell factory for ectoine biosynthesis. The scalability of this system was demonstrated using an array of ten 25 cm2 electrochemical cells equipped with a high‐performance carbon‐supported isolated copper catalyst. The system consistently generated methane at the cathode under a total partial current of approximately −37 A (~175 mmolCH4 h–1) and O2 at the anode under a total partial current of approximately 62 A (~583 mmolO2 h–1). This output met the requirements of a 3‐L bioreactor, even at maximum CH4 and O2 consumption, resulting in the high‐yield conversion of CO2 to ectoine (1146.9 mg L–1). This work underscores the potential of electrifying the biosynthesis of valuable products from CO2, providing a sustainable avenue for biomanufacturing and energy storage.

Scalable Electro-Biosynthesis of Ectoine from Greenhouse Gases.

Converting greenhouse gases into valuable products has become a promising approach for achieving a carbon-neutral economy and sustainable development. However, the conversion efficiency depends on the energy yield of the substrate. In this study, we developed an electro-biocatalyticsystem by integrating electrochemical and microbial processes to upcycle CO2 into a valuable product (ectoine) using renewable energy. This system initiates the electrocatalytic reduction of CO2 to methane, an energy-dense molecule, which then serves as an electrofuel to energize the growth of an engineered methanotrophic cell factory for ectoine biosynthesis. The scalability of this system was demonstrated using an array of ten 25 cm2 electrochemical cells equipped with a high-performance carbon-supported isolated copper catalyst. The system consistently generated methane at the cathode under a total partial current of approximately -37 A (~175 mmolCH4 h-1) and O2 at the anode under a total partial current of approximately 62 A (~583 mmolO2 h-1). This output met the requirements of a 3-L bioreactor, even at maximum CH4 and O2 consumption, resulting in the high-yield conversion of CO2 to ectoine (1146.9 mg L-1). This work underscores the potential of electrifying the biosynthesis of valuable products from CO2, providing a sustainable avenue for biomanufacturing and energy storage.

Numerical analysis of Gaussian potential patches model depending on the substrate doping in inhomogeneous Schottky barrier diodes over a wide temperature range

The current-voltage (I-V-T) characteristics of an inhomogeneous n-type GaAs Schottky barrier diode have been investigated by numerical analysis using the modified thermionic emission (TE) current equation by Tung in the 40–320 K range at 40 K intervals. This total current (TC) equation consists of TE current and the patch current components. The patch current dominates through the low Schottky barrier height patches at low temperatures. From the I-V-T characteristics given for three different standard deviations (σ) at each substrate doping value Nd, we have determined the temperatures at which the patch current begins to dominate. The starting temperature of the patch current has decreased as the σ and Nd values decrease. It has been seen that the temperature at which the patch current component begins to dominate is about 120, 80, and 60 K for σ4, σ3, and σ2 at Nd=1.0×1014cm−3 or Nd=1.0×1015cm−3, respectively; 160, 120, and 80 K at Nd=5.0×1015cm−3; and 200, 160, and 80 K at Nd=1.0×1016cm−3, respectively. Moreover, for the substrate with high doping, it has been observed that the I-V curve of the patch current component or the TC shifts toward higher voltages than the expected position at low temperatures. Thus, from the I-V-T characteristics, it has appeared that Tung’s pinch-off model tends to be more applicable to lightly doped semiconductors. Moreover, the TC equation should be used at high temperatures because the I-V curves at high temperatures belong to the TE component, and the patch current expression without the TE component should be especially used for fit to the experimental curves at low temperatures.

The extremely strong non-neutralized electric currents of the unique solar active region NOAA 13664

Context. In May 2024, the extremely complex active region National Oceanic and Atmospheric Administration (NOAA) 13664 produced the strongest geomagnetic storm since 2003. Aims.The aim of this study is to explore the development of the extreme magnetic complexity of NOAA 13664 in terms of its photospheric electric current. Methods. The non-neutralized electric current was derived from photospheric vector magnetograms, provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The calculation method is based on image processing, thresholding, and error analysis. The spatial and temporal evolution of the non-neutralized electric current of the region as well as its constituent subregions was examined. For context, a comparison with other complex, flare-prolific active regions is provided. Results. Active region NOAA 13664 was formed by the emergence and interaction of three subregions, two of which were of notable individual complexity. It consisted of numerous persistent, current-carrying magnetic partitions that exhibited periods of conspicuous motions and strongly increasing electric current at many locations within the region. These periods were followed by intense and repeated flaring. The total unsigned non-neutralized electric currents and average injection rates reached 5.95 ⋅ 1013 A and 1.5 ⋅ 1013 A/day, and are the strongest observed so far, significantly surpassing other super-active regions of Solar Cycle 24 and 25. Conclusions. Active region NOAA 13664 presents a unique case of complexity. Further scrutiny of the spatial and temporal variation of the net electric currents during the emergence and development of super-active regions is paramount to understand the origin of complex regions and adverse space weather.

A variational method for the sheath potential of hypersonic leading edges with space-charge limitations

Electron transpiration cooling for the leading edges (LE) of hypersonic aircraft utilizes thermionic emission; however, space-charge effects limit the electron emission rate, potentially diminishing the efficiency of this cooling mechanism. We develop a variational weak form of the Poisson equation that describes the sheath potential and then numerically solve it using the finite element method. This formulation has two main benefits: (1) the space-charge limit condition can be incorporated as a constraint and (2) it allows for the analysis of three-dimensional geometries with complex boundary conditions. We demonstrate that the current emitted from the surface of an LE is generally a small fraction of the Child–Langmuir limit due to space charge. We then propose several methods to enhance the emitted current from the surface and to boost the cooling effect of thermionic emission. These include increasing the plasma density, applying a negative surface potential, and using fringe fields under suitable geometric conditions. For a LaB6 emitting LE, the total emitted current is shown to be minimal and independent of the temperature of a surface with floating potential. However, when a negative potential is applied and the surface is heated, the emitted current follows the Richardson–Dushman relationship up to a critical temperature, beyond which it remains constant. At an applied surface potential of −5 V, the critical temperature is around 1700 K.

A data-driven model for the field emission from broad-area electrodes

Electron emission from cathodes in high field gradients is a quantum tunneling effect. The 1928 Fowler–Nordheim field emission (FE) equation and the 1956 Murphy–Good FE equation have traditionally been key in describing cold field emissions, offering estimates for emitters for almost a century. Nevertheless, applying FE theory in practice is often constrained by the lack of data on the distribution and geometry of the emission sites. Predictions become more challenging with an uneven electric field distribution at the cathode surface. Consequently, FE formulations are frequently calibrated using current–voltage data after test, limiting their efficacy as true predictive models.This study develops an alternative model for field emission using a data-driven predictive approach based on (1) vast experimental data, (2) electrostatic simulations of the cathode surface, and (3) detailed material and geometry properties, which together overcome these limitations. The objective of this work is to develop and harness this comprehensive dataset to train a machine learning model capable of providing precise predictions of the cathode current in order to further the understanding and application of field emission phenomena. More than 259 h of experimental data have been processed to train and benchmark some of the well-known machine learning models. After two stages of optimization, a coefficient of determination >98% is achieved in the prediction total field emission current using ensemble models.

Моделирование суточной вариации атмосферного электрического поля в турбулентном приземном слое

An electrodynamic model of the electric field dynamics behavior in the atmospheric surface layer due to the action of local factors is presented. The model was obtained by reducing the differential equations system of the electrode effect to the so-called total current equation, which is a second-order equation of parabolic type, considered in a two-dimensional space-time region. The total current equation allows us to connect the set of main factors influencing on the atmospheric surface layer electric field state: conduction current, turbulent current and current arising as a result of convective processes in the atmosphere with the so-called total current in the nearground layer, reflecting changes in the ionospheric potential. The described method provides significant advantages in research, since within the framework of one model it allows the formulation of various problems of electrodynamics of the surface layer and a comparative analysis of the influence on the behavior of the electric field in the surface layer of both individual factors and their combinations. Based on meteorological data at the Peak Cheget mountain station in the Elbrus region, a daily variation in the values of the turbulent diffusion coefficient was constructed. The resulting analytical dependence was used to solve the equation for the total current in the surface layer, provided that it is constant at the upper boundary of the turbulent electrode layer. Solutions obtained by the Fourier method describe the daily variation of the electric field strength depending on the degree of turbulent mixing in the atmosphere. The appearance of a time shift in daily extrema, a change in their amplitude, and the appearance of additional extrema depending on the elec-tric field values have been established. All of these effects are comparable to the global unitary variation and increase with increasing electric field.

Evolution of Characteristics of Vertical Electric Current and Magnetic Field in Active Regions of the Sun and their Relation to Powerful Flares

The study of evolution of magnetic field and electric currents in active regions of the Sun over a long-time interval is of interest for understanding the processes of accumulation and release of energy in them, leading to various phenomena that affect space weather. In this work, based on the photospheric vector magnetograms of the Helioseismic and Magnetic Imager instrument aboard the Solar Dynamics Observatory, an analysis was made of the evolution of a number of characteristics of the magnetic field and vertical electric current in three active regions 11158, 11675, and 12673 that produced class M and X flares, during the time from their origin in the Eastern hemisphere, during the passage through the solar disk, and until the disappearance near the Western limb with a step of 2 hours. The characteristics under consideration included: the power-law exponent of the probability density function of the absolute value of the vertical electric current density, the maximum of the absolute value of the vertical current density, the signed and unsigned total vertical currents and the unsigned total vertical and horizontal magnetic fluxes, the energy of the nonlinear force-free and potential magnetic fields, the free magnetic energy, and the number of islands with strong vertical current. Some regularities in the behaviour of the characteristics under consideration are found, in particular regarding the occurrence of solar flares. The correlation coefficients between pairs of these characteristics are calculated. Additionally, M. Aschwanden’s approach is shown to be promising for predicting the maximum X-ray class of a flare based on the calculation of the energy of the potential magnetic field in active regions. The results obtained can be used to predict powerful solar flares.

CO2-mediated bicarbonate conversion to concentrated formate in a CEM-based electrolyzer

Renewable energy-driven bicarbonate conversion to valuable chemicals presents an attractive strategy for mitigating CO2 emissions, as bicarbonate can be efficiently generated from the capture of atmospheric CO2 using alkaline solutions with reactive absorption. In this work, we present a CO2-mediated bicarbonate conversion to pure formate using a cation exchange membrane-based electrolyzer with a 25 cm2 electrode area. Our electrolysis achieved selectivities exceeding 75% for formate at a total current of 2.5 A, achieving formate concentrations up to 1.2 M and yields as high as 95% over extended periods. The techno-economic assessment confirmed the economic viability of the process, highlighting the potential for bicarbonate electrolysis as a sustainable method for producing valuable chemicals.

Comprehensive analysis of charge carriers dynamics through the honeycomb structure of graphite thin films and polymer graphite with applications in cold field emission and scanning tunneling microscopy

Polymer graphite electron sources have performed satisfactorily as field emission emitters and scanning tunneling microscopy probes in the past few years. However, the emission process was characterized by limited total emission currents. This paper introduces the elemental, vibrational, electronic structure, and optical analysis of polymer graphite and glass-graphite composite field emission cathodes to study these limitations. Moreover, the field emission characteristics are studied including the changes in the potential energy barrier of the used materials and structures. Among the studied structures, the cathodes prepared from graphite thin films deposited on a micropointed glass substrate (film-GMF) showed superior performance as random field emission arrays. This includes obtaining much higher emission current values ≈ 20 times) and lower threshold voltages ≈ 1/2) compared to the results obtained from polymer graphite samples. The enhancement factor in such emitters is believed to be the three-dimensional honeycomb structure of graphite. Moreover, the study includes applying graphite coatings to tungsten nano-field emission cathodes and scanning tunneling microscopy probes, which improves the performance of such cathodes/probes in both microscopic techniques.

A new airborne system for simultaneous high-resolution ocean vector current and wind mapping: first demonstration of the SeaSTAR mission concept in the macrotidal Iroise Sea

Abstract. Coastal seas, shelf seas and marginal ice zones are dominated by small-scale ocean surface dynamic processes that play a vital role in the transport and exchange of climate-relevant properties such as carbon, heat, water and nutrients between land, ocean, ice and atmosphere. Mounting evidence indicates that ocean scales below 10 km have far-ranging impacts on air–sea interactions, lateral ocean dispersion, vertical stratification, ocean carbon cycling and marine productivity – governing exchanges across key interfaces of the Earth system, the global ocean, and atmosphere circulation and climate. Yet, these processes remain poorly observed at the fine spatial and temporal scales necessary to resolve them. The Ocean Surface Current Airborne Radar (OSCAR) is a new airborne instrument with the capacity to inform these questions by mapping vectorial fields of total ocean surface currents and winds at high resolution over a wide swath. Developed for the European Space Agency (ESA), OSCAR is the airborne demonstrator of the satellite mission concept SeaSTAR, which aims to map total surface current and ocean wind vectors with unprecedented accuracy, spatial resolution and temporal revisit across all coastal seas, shelf seas and marginal ice zones. Like SeaSTAR, OSCAR is an active microwave synthetic aperture radar along-track interferometer (SAR-ATI) with optimal three-azimuth sensing enabled by unique highly squinted beams. In May 2022, OSCAR was flown over the Iroise Sea, France, in its first scientific campaign as part of the ESA-funded SEASTARex project. The campaign successfully demonstrated the capabilities of OSCAR to produce high-resolution 2D images of total surface current vectors and near-surface ocean vector winds, simultaneously, in a highly dynamic, macrotidal coastal environment. OSCAR current and wind vectors show excellent agreement with ground-based X-band-radar-derived surface currents, numerical model outputs and NovaSAR-1 satellite SAR imagery, with root mean square differences from the X-band radar better than 0.2 m s−1 for currents at 200 m resolution. These results are the first demonstration of simultaneous retrieval of total current and wind vectors from a high-squint three-look SAR-ATI instrument and the first geophysical validation of the OSCAR and SeaSTAR observing principle. OSCAR presents a remarkable new ocean observing capability to support the study of small-scale ocean dynamics and air–sea interactions across the Earth's coastal, shelf and polar seas.

An Improved Collaborative Control Scheme to Resist Grid Voltage Unbalance for BDFG-Based Wind Turbine

This article presents an improved collaborative control to resist grid voltage unbalance for brushless doubly fed generator (BDFG)-based wind turbine (BDFGWT). The mathematical model of grid-connected BDFG including machine side converter (MSC) and grid side converter (GSC) in the αβ reference frame during unbalanced grid voltage condition is established. On this base, the improved collaborative control between MSC and GSC is presented. Under the control, the control objective of the whole BDFGWT system, including canceling the pulsations of electromagnetic torque and the unbalance of BDFGWT’s total currents, pulsations of BDFGWT’s total powers are capable of being realized; therefore, the control capability of BDFGWT to resist unbalanced grid voltage is greatly improved. Moreover, improved single-loop current controllers adopting PR regulators are proposed for both MSC and GSC where the sequence extractions for both MSC and GSC currents are not needed any more, and hence the proposed control is much simpler. In addition, the transient characteristics are also improved. Moreover, in order to achieve the decoupling control of current and average power, current controller also adopts the feedforward control approach. Case studies for a two MW BDFGWT system are implemented, and the results verify that the presented control is capable of effectively improving the control capability for BDFGWT to resist grid voltage unbalance and exhibit good stable and dynamic control performances.

Mapping Kinetic Energy Hotspots in the Persian Gulf and Oman Sea Using Surface Current Derived by Geodetic Observations and Data Assimilation

Harnessing ocean kinetic energy has emerged as a promising renewable energy solution in recent years. However, identifying optimal locations for extracting this energy remains a significant challenge. This study presents a novel scheme to estimate the total surface current (TSC) as permanent surface current by integrating geodetic data and in-situ measurements. The TSC is typically a combination of the geostrophic current, derived from dynamic topography, and the Ekman current. We utilize NOAA’s Ekman current data to complement the geostrophic current and obtain the TSC. To further enhance the accuracy of the TSC estimates, we employ a 3DVAR data assimilation method, incorporating local current meter observations. The results are verified against two control current meter stations. The data-assimilation process resulted in an improvement of 4 to 15 cm/s in the precision of calculated TSC. Using the assimilated TSC data, we then assess the kinetic energy potential and identify six regions with the most significant promise for marine kinetic energy extraction. This innovative approach can assist researchers and policymakers in targeting the most suitable locations for harnessing renewable ocean energy.

A numerical method for calculating the driven current of neutral beam injection in tokamaks

The current driven by neutral beam injection in tokamak is calculated, and the slowing-down distribution function of the fast ion is obtained by the backward Euler iteration method, including the pitch angle scattering collision. This study reveals that when the pitch-angle cosine is small, the trapped fast-ion current significantly contributes to the total driven current, particularly when the neutral beam is injected perpendicularly. In such cases, the current densities of passing and trapped ions are of the same order of magnitude, with the trapped fast-ion current contributing over 10% to the total neutral beam-driven current. This results in a parabolic profile of the total current in the radial direction, promoting the formation of a negative shear equilibrium structure in the core of the tokamak plasma. The numerical approach was validated against the NUBEAM code while considering electron shielding effects and applied to calculate the neutral beam-driven current in multiple tokamaks. The influence of pitch-angle cosine and neutral beam injection power on the driven current was studied at different radial positions.