Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function Estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Measuring the Alfvénic nature of the interstellar medium: Velocity anisotropy revisited

The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel and Lazarian method to estimate the Alfvén Mach number using the structure function anisotropy in velocity centroid data from Position-Position-Velocity maps. We utilize three-dimensional magnetohydrodynamic simulations of fully developed turbulence, with a large range of sonic and Alfvénic Mach numbers, to produce synthetic observations of velocity centroids with observational characteristics such as thermal broadening, cloud boundaries, noise, and radiative transfer effects of carbon monoxide. In addition, we investigate how the resulting anisotropy-Alfvén Mach number dependency found in Esquivel and Lazarian might change when taking the second moment of the Position-Position-Velocity cube or when using different expressions to calculate the velocity centroids. We find that the degree of anisotropy is related primarily to the magnetic field strength (i.e., Alfvén Mach number) and the line-of-sight orientation, with a secondary effect on sonic Mach number. If the line of sight is parallel to up to ≈45 deg off of the mean field direction, the velocity centroid anisotropy is not prominent enough to distinguish different Alfvénic regimes. The observed anisotropy is not strongly affected by including radiative transfer, although future studies should include additional tests for opacity effects. These results open up the possibility of studying the magnetic nature of the ISM using statistical methods in addition to existing observational techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

CONSTRAINTS ON SPATIAL VARIATIONS IN THE FINE-STRUCTURE CONSTANT FROM PLANCK

We use the cosmic microwave background (CMB) anisotropy data from Planck to constrain the spatial fluctuations of the fine-structure constant α at a redshift of 1100. We use a quadratic estimator to measure the four-point correlation function of the CMB temperature anisotropies and extract the angular power spectrum fine-structure constant spatial variations projected along the line of sight at the last scattering surface. At tens of degree angular scales and above, we constrain the fractional rms fluctuations of the fine-structure constant to be (δα/α){sub rms} < 3.4 × 10{sup –3} at the 68% confidence level. We find no evidence for a spatially varying α at a redshift of 10{sup 3}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

JET LUMINOSITY OF GAMMA-RAY BURSTS: THE BLANDFORD–ZNAJEK MECHANISM VERSUS THE NEUTRINO ANNIHILATION PROCESS

A neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engine of gamma-ray bursts (GRBs). Two mechanisms, i.e., the Blandford–Znajek (BZ) mechanism and the neutrino annihilation process, are generally considered to power GRBs. Using the analytic solutions from Xue et al. and ignoring the effects of the magnetic field configuration, we estimate the BZ and neutrino annihilation luminosities as functions of the disk masses and BH spin parameters to contrast the observational jet luminosities of GRBs. Our results show that although the neutrino annihilation processes could account for most GRBs, the BZ mechanism is more effective, especially for long-duration GRBs. Actually, if the energy of the afterglows and flares of GRBs is included, then the distinction between these two mechanisms is more significant. Furthermore, massive disk mass and high BH spin are beneficial for powering the high luminosities of GRBs. Finally, we discuss possible physical mechanisms that could enhance the disk mass or neutrino emission rate of NDAFs and the relevant difference between these two mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Determining the anisotropy and exchange parameters of polycrystalline spin-1 magnets

Abstract Although low-dimensional S = 1 antiferromagnets remain of great interest, difficulty in obtaining high-quality single crystals of the newest materials hinders experimental research in this area. Polycrystalline samples are more readily produced, but there are inherent problems in extracting the magnetic properties of anisotropic systems from powder data. Following a discussion of the effect of powder-averaging on various measurement techniques, we present a methodology to overcome this issue using thermodynamic measurements. In particular we focus on whether it is possible to characterise the magnetic properties of polycrystalline, anisotropic samples using readily available laboratory equipment. We test the efficacy of our method using the magnets [Ni(H 2 O) 2 (3,5-lutidine) 4 ](BF 4 ) 2 and Ni(H 2 O) 2 (acetate) 2 (4-picoline) 2 , which have negligible exchange interactions, as well as the antiferromagnet [Ni(H 2 O) 2 (pyrazine) 2 ](BF 4 ) 2 , and show that we are able to extract the anisotropy parameters in each case. The results obtained from the thermodynamic measurements are checked against electron-spin resonance and neutron diffraction. We also present a density functional method, which incorporates spin–orbit coupling to estimate the size of the anisotropy in [Ni(H 2 O) 2 (pyrazine) 2 ](BF 4 ) 2 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Fine-scale vegetation composition and structure shape spatiotemporal variation in surface albedo across a low Arctic tundra landscape

The unprecedented rate of warming in the Arctic is driving changes in the structure and composition of tundra vegetation. Increases in deciduous tall shrub cover, height, and density are of particular concern, as these changes alter local surface albedo in ways that could amplify effects on the regional surface energy budget (SEB). Despite this importance, significant uncertainties remain in understanding the interplay between fine-scale vegetation patterns and emergent albedo dynamics across space and time. Here, we address these uncertainties by (1) quantifying spatiotemporal variation in surface shortwave albedo and (2) determining the relative influence of fine-scale vegetation composition, structure, and environmental conditions on albedo across a representative low-Arctic tundra landscape on Alaska’s Seward Peninsula. To do this, we synthesized multi-scale, multi-platform remote sensing observations, including a novel Landsat-derived albedo time series, a fine-scale map of Arctic plant functional type (PFT) fractional cover, and airborne LiDAR estimates of canopy height and topography. We show that there are substantial reductions in winter albedo for pixels dominated by tall, woody PFTs (28.13%) relative to pixels dominated by non-woody vegetation, but almost no change in summer albedo (3% increase). Further, we identified a unimodal trend in the relationship between canopy height and the timing of the springtime transition from high (snowy) to low (leafy) albedo (peak at 5.5 m), possibly because of competing ‘snow-fence’ and ‘protrusion’ snow-shrub interactions. To explore the primary drivers of albedo, we constructed a random forest model and found that canopy height and the fractional cover of woody PFTs were as- or more important predictors of winter albedo than topographic features. These findings provide strong evidence for the impacts of local vegetation characteristics on regional surface albedo, highlighting the need for better quantification of snow-shrub interactions to accurately predict the Arctic’s SEB under future environmental change.

Arctic↗
DOE OSTI2021

Ionizing radiation from AGNs at z > 3.3 with the Subaru Hyper Suprime-Cam Survey and the CFHT Large Area U-band Deep Survey (CLAUDS)

ABSTRACT We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum, LyC) escape fraction from active galactic nuclei (AGNs) at z ∼ 3–4. For 94 AGNs with spectroscopic redshifts at 3.3 < z < 4.0, we use their U/i-band flux ratios to estimate LyC transmission of individual AGNs. The distribution of their LyC transmission shows values lower than the range of LyC transmission values for IGM of the same redshift range, which suggests that LyC escape fraction of AGNs at z > 3.3 is considerably lower than unity in most cases. We do not find any trend in LyC transmission values depending on their UV luminosities. Based on the photometry of stacked images we find the average flux ratio of LyC and non-ionizing UV photons escaping from the objects (fLyC/fUV)out = 0.182 ± 0.043 for AGNs at 3.3 < z < 3.6, which corresponds to LyC escape fraction fesc = 0.303 ± 0.072 if we assume a fiducial intrinsic SED of AGN. Based on the estimated LyC escape fraction and the UV luminosity function (UVLF) of AGNs, we argue that UV-selected AGNs’ contribution to the LyC emissivity at the epoch is minor, although the size of their contribution largely depends on the shape of the UVLF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Non-local contribution from small scales in galaxy–galaxy lensing: comparison of mitigation schemes

ABSTRACT Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3 × 2pt, had to discard a lot of signal to noise from small scales due to our inability to accurately model non-linearities and baryonic effects. Galaxy–galaxy lensing, or the position–shear correlation between lens and source galaxies, is one of the three two-point correlation functions that are included in such analyses, usually estimated with the mean tangential shear. However, tangential shear measurements at a given angular scale θ or physical scale R carry information from all scales below that, forcing the scale cuts applied in real data to be significantly larger than the scale at which theoretical uncertainties become problematic. Recently, there have been a few independent efforts that aim to mitigate the non-locality of the galaxy–galaxy lensing signal. Here, we perform a comparison of the different methods, including the Y-transformation, the point-mass marginalization methodology, and the annular differential surface density statistic. We do the comparison at the cosmological constraints level in a combined galaxy clustering and galaxy–galaxy lensing analysis. We find that all the estimators yield equivalent cosmological results assuming a simulated Rubin Observatory Legacy Survey of Space and Time (LSST) Year 1 like set-up and also when applied to DES Y3 data. With the LSST Y1 set-up, we find that the mitigation schemes yield ∼1.3 times more constraining S8 results than applying larger scale cuts without using any mitigation scheme.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Constraining primordial non-Gaussianity from the large scale structure two-point and three-point correlation functions

Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by f NL ⁠, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spectrum. We propose a framework to measure f NL using the configuration space two-point correlation function (2pcf) monopole and three-point correlation function (3pcf) monopole of survey tracers. Our model estimates the effect of the scale-dependent bias induced by the presence of PNG on the 2pcf and 3pcf from the clustering of simulated dark matter haloes. We describe how this effect may be scaled to an arbitrary tracer of the cosmological matter density. The 2pcf and 3pcf of this tracer are measured to constrain the value of f NL ⁠. In LSS surveys, the effect of imaging systematics on two-point statistics is often degenerate with the PNG signal. Our proposed model employs three-point statistics primarily to break this degeneracy. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our method. Our forecast indicates the ability to constrain f NL to a precision of σf NL ≈ 22 with one year of DESI survey data, as well as the ability to constrain the imaging systematic weights in situ.

early Universe↗
DOE OSTI2019

Likelihood non-Gaussianity in large-scale structure analyses

The standard present-day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference – that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. power spectrum multipole ($P_ℓ$) analysis and the Sinha et al. group multiplicity function ($ζ$) analysis. Using non-parametric divergence estimators on mock catalogues originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the $P_ℓ$ and $ζ$ likelihoods. We then use Gaussian mixture density estimation and independent component analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian pseudo-likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. and Sinha et al. parameters. Likelihood non-Gaussianity shifts the fσ8 constraint by -0.44σ, but otherwise does not significantly impact the overall parameter constraints of Beutler et al. For the $ζ$ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of the Sinha et al. halo occupation parameters. For log $M_1$ and α, the posteriors are shifted by +0.43σ and -0.51σ and broadened by 42 per cent and 66 per cent, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Mass functions, luminosity functions, and completeness measurements from clustering redshifts

Abstract This paper presents stellar mass functions and i-band luminosity functions for Sloan Digital Sky Survey (SDSS) galaxies with i < 21 using clustering redshifts. From these measurements, we also compute targeting completeness measurements for the Baryon Oscillation Spectroscopic Survey (BOSS). Clustering redshifts is a method of obtaining the redshift distribution of a sample of galaxies with only photometric information by measuring the angular cross-correlation with a spectroscopic sample in different redshift bins. We construct a spectroscopic sample containing data from the BOSS + eBOSS surveys, allowing us to recover redshift distributions from photometric data out to z ≃ 2.5. We produce k-corrected i-band luminosity functions and stellar mass functions by applying clustering redshifts to SDSS DR8 galaxies in small bins of colour and magnitude. There is little evolution in the mass function between 0.2 < z < 0.8, implying that the most massive galaxies form most of their mass before z = 0.8. These mass functions are used to produce stellar mass completeness estimates for the BOSS, giving a stellar mass completeness of $80{{\ \rm per\ cent}}$ above M⋆ > 1011.4 between 0.2 < z < 0.7, with completeness falling significantly at redshifts higher than 0.7, and at lower masses. Large photometric data sets will be available in the near future (DECaLS, DES, Euclid), so this and similar techniques will become increasingly useful in order to fully utilize these data.

Bates, Dominic J.↗
DOE OSTI2023

Magnetic ground state and perturbations of the distorted kagome Ising metal TmAgGe

Here, we present the magnetic orders and excitations of the distorted kagome intermetallic magnet TmAgGe. Using neutron single crystal diffraction we identify the propagation vectors k = ($\frac{1}{2}$ 0 0) and k = (0 0 0) and determine the magnetic structures of the zero-field and magnetic field-induced phases for H along the $\mathcal{a}$ and [–110] crystal directions. We determine the experimental magnetic field- temperature (H, T)-phase diagram and reproduce it by Monte Carlo simulations of an effective spin exchange Hamiltonian for one distorted kagome layer. Our model includes a strong axial single-ion anisotropy and significantly smaller exchange couplings, which span up to the third-nearest neighbors within the layer. Single crystal inelastic neutron scattering (INS) measurements reveal an almost flat, only weakly dispersive mode around 7 meV that we use alongside bulk magnetization data to deduce the crystal-electric field (CEF) scheme for the Tm 3+ ions. Random phase approximation (RPA) calculations based on the determined CEF wave functions of the two lowest quasidoublets enable an estimation of the interlayer coupling that is compatible with the experimental INS spectra. No evidence for low-energy spin waves associated to the magnetic order was found, which is consistent with the strongly Ising nature of the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Lateral cracks during sliding indentation on various optical materials

A series of static and sliding indentation (ie, scratching) was performed and characterized on a wide range of optical workpiece materials [single crystals of Al 2 O 3 (sapphire), SiC, Y 3 Al 5 O 12 (YAG), CaF 2 , and LiB 3 O 5 (LBO); a SiO 2 –Al 2 O 3 –P 2 O 5 –Li 2 O glass ceramic (Zerodur); and glasses of SiO 2 :TiO 2 (ULE), SiO 2 (fused silica), and P 2 O 5 –Al 2 O 3 –K 2 O–BaO (Phosphate)] at various applied loads using various indenters (Vickers, 10 µm conical, and 200 µm conical). Despite having different load dependencies, the lateral crack depth formed during sliding indentation quantitatively scales with that formed during static indentation, explaining why static indentation has been historically effective in describing various grinding parameters. Depending on the indenter geometry, the amount of residual trench damage (plastic deformation and local fracturing) during sliding indentation was often enhanced by more than an order of magnitude compared with static indentation. A simple ploughing scratch model, which considers both tangential and normal stresses (where the tangential stress is amplified by relatively small tangential contact area), explains this enhancement and other observed trends. Accounting for the high correlation between residual trench depth and volumetric fracturing, the model is extended to estimate the amount of fracture damage as a function of the material properties of the workpiece, indenter geometry, and applied load. Such a model has utility in the design of optimized grinding processes, particularly the abrasive geometry. Finally, at higher loads (>1 N), lateral cracks were often observed to preferentially propagate in the forward scratching direction, as opposed to perpendicular to the scratch as typically observed. High-speed imaging of the scratch process confirms that these cracks propagate ahead of the sliding indenter during the scratching event. Finite element stress analysis suggests the ploughing frictional forces increase the mode I tensile stresses at the leading edge of the sliding indenter explaining the direction of crack propagation of such cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2014

SU-F-19A-04: Dosimetric Evaluation of a Novel CT/MR Compatible Fletcher Applicator for Intracavitary Brachytherapy of the Cervix Uteri

Purpose: To validate a Monte Carlo model and evaluate the dosimetric capabilities of a novel commercial CT/MR compatible Fletcher applicator for cervical cancer brachytherapy. Methods: MCNPX 2.7.0 was used to model the Fletcher CT/MR shielded applicator (FA) and 192Ir HDR source. Energy deposition was calculated with a track length estimator modified by an energy-dependent heating function. A high density polystyrene phantom was constructed with three film pockets for validation of the MCNPX model. Three planes of data were calculated with the MCNPX model corresponding to the three film planes in phantom. The planes were located 1 cm from the most anterior, posterior, and medial extents of the FA right ovoid. Unshielded distributions were calculated by modeling the shielded cells as air instead of the tungsten alloy. A third order polynomial fit to the OD to dose curve was used to convert OD of the three film planes to dose. Each film and MCNPX plane dose distribution was normalized to a point 2 cm from the center of the film plane and in a region of low dose gradient. MCNPX and film were overlaid and compared with a distance-to-agreement criterion of (±2%/±2mm). Shielded and unshielded distributions were overlaid and a percent shielded plot was created. Results: 85.2%, 97.1%, and 96.6% of the MCNPX points passed the (±2%/±2mm) criterion respectively for the anterior, lateral, and posterior film comparison planes. A majority of the points in the anterior plane that exceeded the DTA criterion were either along edges of where the film was cut or near the terminal edges of the film. The percent shielded matrices indicated that the maximum % shielding was 50%. Conclusion: These data confirm the validity of the FA Monte Carlo model. The FA ovoid can shield up to 50% of the dose in the anteroposterior direction.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-G-204-01: BEST IN PHYSICS (IMAGING): Effect of Image Processing Parameters On Nodule Detectability in Chest Radiography

Purpose: To investigate the effect of varying system image processing parameters on lung nodule detectability in digital radiography. Methods: An anthropomorphic chest phantom was imaged in the posterior-anterior position using a GE Discovery XR656 digital radiography system. To simulate lung nodules, a polystyrene board with 6.35mm diameter PMMA spheres was placed adjacent to the phantom (into the x-ray path). Due to magnification, the projected simulated nodules had a diameter in the radiographs of approximately 7.5 mm. The images were processed using one of GE’s default chest settings (Factory3) and reprocessed by varying the “Edge” and “Tissue Contrast” processing parameters, which were the two user-configurable parameters for a single edge and contrast enhancement algorithm. For each parameter setting, the nodule signals were calculated by subtracting the chest-only image from the image with simulated nodules. Twenty nodule signals were averaged, Gaussian filtered, and radially averaged in order to generate an approximately noiseless signal. For each processing parameter setting, this noise-free signal and 180 background samples from across the lung were used to estimate ideal observer performance in a signal-known-exactly detection task. Performance was estimated using a channelized Hotelling observer with 10 Laguerre-Gauss channel functions. Results: The “Edge” and “Tissue Contrast” parameters each had an effect on the detectability as calculated by the model observer. The CHO-estimated signal detectability ranged from 2.36 to 2.93 and was highest for “Edge” = 4 and “Tissue Contrast” = −0.15. In general, detectability tended to decrease as “Edge” was increased and as “Tissue Contrast” was increased. A human observer study should be performed to validate the relation to human detection performance. Conclusion: Image processing parameters can affect lung nodule detection performance in radiography. While validation with a human observer study is needed, model observer detectability for common tasks could provide a means for optimizing image processing parameters.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Damage Identification of Wind Turbine Blades Using Piezoelectric Transducers

This paper presents the experimental results of active-sensing structural health monitoring (SHM) techniques, which utilize piezoelectric transducers as sensors and actuators, for determining the structural integrity of wind turbine blades. Specifically, Lamb wave propagations and frequency response functions at high frequency ranges are used to estimate the condition of wind turbine blades. For experiments, a 1 m section of a CX-100 blade is used. The goal of this study is to assess and compare the performance of each method in identifying incipient damage with a consideration given to field deployability. Overall, these methods yielded a sufficient damage detection capability to warrant further investigation. This paper also summarizes the SHM results of a full-scale fatigue test of a 9 m CX-100 blade using piezoelectric active sensors. This paper outlines considerations needed to design such SHM systems, experimental procedures and results, and additional issues that can be used as guidelines for future investigations.

Choi, Seong-Won↗
DOE OSTI2017

Quinoa, Version 0.2

Quinoa is a set of computational tools that enables research and numerical analysis in fluid dynamics. At this time it remains a test-bed to experiment with various algorithms using fully asynchronous runtime systems. Currently, Quinoa consists of the following tools: (1) Walker, a numerical integrator for systems of stochastic differential equations in time. It is a mathematical tool to analyze and design the behavior of stochastic differential equations. It allows the estimation of arbitrary coupled statistics and probability density functions and is currently used for the design of statistical moment approximations for multiple mixing materials in variable-density turbulence. (2) Inciter, an overdecomposition-aware finite element field solver for partial differential equations using 3D unstructured grids. Inciter is used to research asynchronous mesh-based algorithms and to experiment with coupling asynchronous to bulk-synchronous parallel code. Two planned new features of Inciter, compared to the previous release (LA-CC-16-015), to be implemented in 2017, are (a) a simple Navier-Stokes solver for ideal single-material compressible gases, and (b) solution-adaptive mesh refinement (AMR), which enables dynamically concentrating compute resources to regions with interesting physics. Using the NS-AMR problem we plan to explore how to scale such high-load-imbalance simulations, representative of large production multiphysics codes, to very large problems on very large computers using an asynchronous runtime system. (3) RNGTest, a test harness to subject random number generators to stringent statistical tests enabling quantitative ranking with respect to their quality and computational cost. (4) UnitTest, a unit test harness, running hundreds of tests per second, capable of testing serial, synchronous, and asynchronous functions. (5) MeshConv, a mesh file converter that can be used to convert 3D tetrahedron meshes from and to either of the following formats: Gmsh, http://www.geuz.org/gmsh, Netgen, http://sourceforge.net/apps/mediawiki/netgen-mesher, ExodusII, http://sourceforge.net/projects/exodusii, HyperMesh, http://www.altairhyperworks.com/product/HyperMesh.

Bakosi, Jozsef↗
DOE OSTI2026

Quantum Stochastic Programming [SWR-26-040]

The Quantum Stochastic Programming tool contains quantum computing algorithms for two-stage stochastic optimization, with a focus on the Unit Commitment (UC) problem in power systems. The algorithms combine Discrete Quantum Annealing (DQA) with Quantum Amplitude Estimation (QAE) to compute expected-value objective functions over a probability distribution of wind-power scenarios. Based on: arXiv 2402.15029 - "Quantum algorithms for the two-stage stochastic unit commitment problem"

Maack, Jonathan [National Laboratory of the Rockie↗
DOE OSTI2023

Incidence and predictors of tuberculosis among children receiving antiretroviral therapy in the Wolaita Zone: A retrospective cohort study

Tuberculosis is the leading cause of morbidity and mortality among children living with the human immunodeficiency virus (HIV), mainly in sub-Saharan Africa, including Ethiopia. Tuberculosis remains a significant health concern for HIV-positive children in Ethiopia. There is a paucity of data on the incidence and predictors of tuberculosis among children living with HIV on antiretroviral therapy in the Wolaita zone. Hence, this study aimed to assess the incidence and predictors of tuberculosis among children living with HIV on antiretroviral therapy in the Wolaita zone between January 2010 to December 2020. A retrospective cohort study was conducted among 389 children receiving antiretroviral therapy in Wolaita zone health facilities between January 2010 to December 2020. The checklist was adapted from the standardized antiretroviral treatment (ART) follow-up form currently used by the institutions’ ART clinics. The Kaplan-Meier survival function and Log-rank were used to estimate the survival for each categorical variable to compare the survival between different exposure groups. Both bivariable and multivariable parametric survival Gompertz models were fitted to identify predictors of tuberculosis among HIV-positive children. The association was summarized using an adjusted hazard ratio (AHR), and statistical significance was declared at 95% CI and p-value < 0.05. The goodness of the model fit was assessed using a Cox-Snell residual plot. The incidence rate of tuberculosis among children living with HIV was 3.5 (95% CI 2.7–4.5) per 100 child years. World Health Organization clinical stage III or IV (AHR = 2.31, 95% CI [1.26, 4.22]), hemoglobin level <10 g/dL (AHR = 2.87, 95% CI [1.51, 5.45]), fair or poor ART adherence (AHR = 4.4, 95% CI[2.18, 9.05]), underweight (AHR = 2.55, 95% CI [1.45, 4.51]), age >10 years (AHR = 3.62; 95% CI [1.29, 10.0]), and cotrimoxazole preventive therapy (AHR = 0.23; 95% CI [0.08, 0.60]) were among the independent predictors of TB occurrence. The incidence of tuberculosis among children on ART was high. HIV-positive children presenting with advanced disease staging (III and IV), anemia, “fair” and “poor” ART adherence, underweight, age above ten years, and not receiving cotrimoxazole preventive therapy were at higher risk of TB. Therefore, counseling on ART adherence, early diagnosis, and prompt treatment of anemia and malnutrition are recommended to avert tuberculosis.

60 APPLIED LIFE SCIENCES↗