SEARCH · Search NASA

Results for “Reproducibility”

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.

SourceResearch connectionFrom yearThrough year

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

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2020

CMB distance priors revisited: effects of dark energy dynamics, spatial curvature, primordial power spectrum, and neutrino parameters

As a physical and sufficient compression of the full CMB data, the CMB distance priors, or shift parameters, have been widely used and provide a convenient way to include CMB data when obtaining cosmological constraints. In this paper, we revisit this data vector and examine its stability under different cosmological models. We find that the CMB distance priors are an accurate substitute for the full CMB data when probing dark energy dynamics. This is true when the primordial power spectrum model is directly generalized from the power spectrum of the model used in the derivation of the distance priors from the CMB data. We discover a difference when a non-flat model with the untilted primordial inflation power spectrum is used to measure the distance priors. This power spectrum is a radical change from the more conventional tilted primordial power spectrum and violates fundamental assumptions for the reliability of the CMB shift parameters. We also investigate the performance of CMB distance priors when the sum of neutrino masses Σm ν and the effective number of relativistic species N eff are allowed to vary. Our findings are consistent with earlier results: the neutrino parameters can change the measurement of the sound horizon from CMB data, and thus the CMB distance priors. We find that when the neutrino model is allowed to vary, the cold dark matter density ω c and N eff need to be included in the set of parameters that summarize CMB data, in order to reproduce the constraints from the full CMB data. As a result, we present an updated and expanded set of CMB distance priors which can reproduce constraints from the full CMB data within 1σ, and are applicable to models with massive neutrinos, as well as non-standard cosmologies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Static response and Love numbers of Schwarzschild black holes

We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild-(A)dS background in arbitrary dimensions. Here, we then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers — in particular that they vanish for all types of perturbation in four spacetime dimensions — but also leads to new results. For instance, we find that neutral Schwarzschild black holes polarize in the presence of an electromagnetic background in any number of spacetime dimensions except four. Moreover, we calculate for the first time black hole Love numbers for vector-type gravitational perturbations in higher dimensions and find that they generically do not vanish. Along the way, we shed some light on an apparent discrepancy between previous results in the literature, and clarify some aspects of the matching between perturbative calculations of static response on a Schwarzschild background and the point-particle effective theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Synthetic spectra for Lyman- α forest analysis in the Dark Energy Spectroscopic Instrument

Synthetic data sets are used in cosmology to test analysis procedures, to verify that systematic errors are well understood and to demonstrate that measurements are unbiased. In this work we describe the methods used to generate synthetic datasets of Lyman-α quasar spectra aimed for studies with the Dark Energy Spectroscopic Instrument (DESI). In particular, we focus on demonstrating that our simulations reproduces important features of real samples, making them suitable to test the analysis methods to be used in DESI and to place limits on systematic effects on measurements of Baryon Acoustic Oscillations (BAO). We present a set of mocks that reproduce the statistical properties of the DESI early data set with good agreement. Additionally, we use a synthetic dataset to forecast the BAO scale constraining power of the completed DESI survey through the Lyman-α forest.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Effect of sawtooth crashes on fast ion distribution in NSTX-U

During the 2016 experimental campaign of NSTX-U (Menard et al 2012 Nucl. Fusion 52 083015), long L-mode and reproducible sawtoothing plasmas have been achieved that were previously not accessible before the upgrade. This provides a good opportunity to study the effect of sawtooth crashes on fast ion confinement and redistribution in spherical tokamaks. The solid-state neutral particle analyzer (SSNPA) and fast ion D-alpha (FIDA) diagnostics on NSTX-U each have a tangentially-viewing instrument and a radially- or vertically-viewing instrument, which are mainly sensitive to passing and trapped fast ions, respectively. It has been observed on both diagnostics that passing particles are strongly expelled from the plasma core to the edge during sawtooth crashes, while trapped fast ions are weakly affected. The tangentially-viewing SSNPA observes large signal spikes at the sawtooth crashes because fast ions move to the edge and charge exchange with edge neutrals. The radially-viewing SSNPA data suggest that there is a small drop of trapped particles in the core. The tangentially-viewing FIDA (t-FIDA) system observes a depletion as large as 25% in the region inside the inversion radius, while an increase at the outer region. There is almost no change in the signals of the vertically-viewing FIDA system. The neutron emission can drop as much as 15% at the sawtooth crashes, accompanied by an increase of edge D α light. Simulations with the Kadomstev and Porcelli sawtooth models have also been performed and compared with the measurements. The full reconnection Kadomtsev model overestimates the neutron rate drop at each sawtooth crash. When tuning the sawtooth input parameters, the partial reconnection Porcelli sawtooth model can qualitatively reproduce the neutron rate drop and t-FIDA signal drop in the core, but it fails to predict the t-FIDA signal increase at the edge.

43 PARTICLE ACCELERATORS↗

DOE OSTI2019

Modification of the Alfvén wave spectrum by pellet injection

Alfvén eigenmodes driven by energetic particles are routinely observed in tokamak plasmas. These modes consist of poloidal harmonics of shear Alfvén waves coupled by inhomogeneity in the magnetic field. In this work, further coupling is introduced by 3D inhomogeneities in the ion density during the assimilation of injected pellets. This additional coupling modifies the Alfvén continuum and discrete eigenmode spectrum. We find that the frequencies of Alfvén eigenmodes drop dramatically when a pellet is injected in JET. From these observations, information about the changes in the ion density caused by a pellet can be inferred. To use Alfvén eigenmodes for MHD spectroscopy of pellet injected plasmas, the 3D MHD codes Stellgap and AE3D were generalised to incorporate 3D density profiles. Additionally, a model for the expansion of the ionised pellet plasmoid along a magnetic field line was derived from the fluid equations. Thereby, the time evolution of the Alfvén eigenfrequency is reproduced. By comparing the numerical frequency drop of a toroidal Alfvén eigenmode (TAE) to experimental observations, the initial ion density of a cigar-shaped ablation region of length 4 cm is estimated to be n * =6.8×10 22 m -3 at the TAE location (r/a≈0.75). The frequency sweeping of an Alfvén eigenmode ends when the ion density homogenises poloidally. Modelling suggests that the time for poloidal homogenisation of the ion density at the TAE position is τ h =18±4 ms for inboard pellet injection, and τ h =26±2 ms for outboard pellet injection. By reproducing the frequency evolution of the elliptical Alfvén eigenmode (EAE), the initial ion density at the EAE location (r/a≈0.9) can be estimated to be n * =4.8×10 22 m -3 . Poloidal homogenisation of the ion density takes 2.7 times longer at the EAE location than at the TAE location for both inboard and outboard pellet injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Modeling the sharp deuterium release from beryllium co-deposits

Many studies have shown a sharp deuterium desorption from beryllium at around 450 K. In an attempt to model this sharp release, a new model was developed, describing the formation and decomposition of beryllium deuteride. Such a model was found to reproduce the sharp desorption peak in the thermal desorption spectrum. Moreover, in contrast to other models, the new model was able to reproduce the constant deuterium desorption from beryllium held at a constant elevated temperature. Further, fitting the modeling results to experimental data allowed us to determine the deuterium solubility limit in beryllium, which indicated an endothermic formation of beryllium deuteride. New values of the deuterium surface recombination coefficient were also obtained, differing from the ones found in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Progress from ASDEX Upgrade experiments in preparing the physics basis of ITER operation and DEMO scenario development

An overview of recent results obtained at the tokamak ASDEX Upgrade (AUG) is given. A work flow for predictive profile modelling of AUG discharges was established which is able to reproduce experimental H-mode plasma profiles based on engineering parameters only. In the plasma center, theoretical predictions on plasma current redistribution by a dynamo effect were confirmed experimentally. For core transport, the stabilizing effect of fast ion distributions on turbulent transport is shown to be important to explain the core isotope effect and improves the description of hollow low-Z impurity profiles. The L–H power threshold of hydrogen plasmas is not affected by small helium admixtures and it increases continuously from the deuterium to the hydrogen level when the hydrogen concentration is raised from 0 to 100%. One focus of recent campaigns was the search for a fusion relevant integrated plasma scenario without large edge localised modes (ELMs). Results from six different ELM-free confinement regimes are compared with respect to reactor relevance: ELM suppression by magnetic perturbation coils could be attributed to toroidally asymmetric turbulent fluctuations in the vicinity of the separatrix. Stable improved confinement mode plasma phases with a detached inner divertor were obtained using a feedback control of the plasma β. The enhanced D α H-mode regime was extended to higher heating power by feedback controlled radiative cooling with argon. The quasi-coherent exhaust regime was developed into an integrated scenario at high heating power and energy confinement, with a detached divertor and without large ELMs. Small ELMs close to the separatrix lead to peeling-ballooning stability and quasi continuous power exhaust. Helium beam density fluctuation measurements confirm that transport close to the separatrix is important to achieve the different ELM-free regimes. Based on separatrix plasma parameters and interchange-drift-Alfvén turbulence, an analytic model was derived that reproduces the experimentally found important operational boundaries of the density limit and between L- and H-mode confinement. Feedback control for the X-point radiator (XPR) position was established as an important element for divertor detachment control. Stable and detached ELM-free phases with H-mode confinement quality were obtained when the XPR was moved 10 cm above the X-point. Investigations of the plasma in the future flexible snow-flake divertor of AUG by means of first SOLPS-ITER simulations with drifts activated predict beneficial detachment properties and the activation of an additional strike point by the drifts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

The role of ETG modes in JET–ILW pedestals with varying levels of power and fuelling

Abstract We present the results of GENE gyrokinetic calculations based on a series of JET–ITER-like-wall (ILW) type I ELMy H-mode discharges operating with similar experimental inputs but at different levels of power and gas fuelling. We show that turbulence due to electron-temperature-gradient (ETGs) modes produces a significant amount of heat flux in four JET–ILW discharges, and, when combined with neoclassical simulations, is able to reproduce the experimental heat flux for the two low gas pulses. The simulations plausibly reproduce the high-gas heat fluxes as well, although power balance analysis is complicated by short ELM cycles. By independently varying the normalised temperature gradients ( ω T e ) and normalised density gradients ( ω n e ) around their experimental values, we demonstrate that it is the ratio of these two quantities η e = ω T e / ω n e that determines the location of the peak in the ETG growth rate and heat flux spectra. The heat flux increases rapidly as η e increases above the experimental point, suggesting that ETGs limit the temperature gradient in these pulses. When quantities are normalised using the minor radius, only increases in ω T e produce appreciable increases in the ETG growth rates, as well as the largest increases in turbulent heat flux which follow scalings similar to that of critical balance theory. However, when the heat flux is normalised to the electron gyro-Bohm heat flux using the temperature gradient scale length L T e , it follows a linear trend in correspondence with previous work by different authors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Parameter dependencies of the separatrix density in low triangularity L-mode and H-mode JET-ILW plasmas

Abstract The midplane electron separatrix density, n e,sep , in JET-ILW L-mode and H-mode low triangularity deuterium fuelled plasmas exhibits a strong explicit dependence on the averaged outer divertor target electron temperature, n e,sep ∼ T e,ot −1/2 . This dependence is reproduced by analytic reversed two point model (rev-2PM), and arises from parallel pressure balance, as well as the ratio of the power and momentum volumetric loss factors, (1 − f cooling )/(1 − f mom-loss ). Quantifying the influence of the (1 − f cooling ) and (1 − f mom-loss ) loss factors on n e,sep has been enabled by measurement estimates of these quantities from L-mode density (fueling) ramps in the outer horizontal, VH(C), and vertical target, VV, divertor configurations. Rev-2PM n e,sep estimates from the extended H-mode and more limited L-mode datasets are recovered to within ±25% of the measurements, with a scaling factor applied to account for use of T e,ot , an averaged quantity, rather than flux tube resolved target values. Both the (1 − f cooling ) and (1 − f mom-loss ) trends and recovery of n e,sep using the rev-2PM formatting are reproduced in EDGE2D-EIRENE L-mode-like and H-mode-like density scan simulations. The general lack of a divertor configuration effect in the JET-ILW n e,sep trends can be attributed to a significant influence of main chamber recycling, which has been shown in the EDGE2D-EIRENE results to moderate n e,sep with respect to changes in divertor neutral leakage imposed by changes in the divertor configuration. The unified n e,sep vs T e,ot trends can, however, be broken if large modifications to the divertor geometry (e.g. complete removal of the outer divertor baffle structure) are introduced in the model. The more pronounced high-field side high density region formation in the VH(C) configuration with reduced clearance to the separatrix does not appear to have a significant influence on the outer midplane separatrix and pedestal parameters when mapped to T e,ot , although conditions at the inner midplane could not be assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas

Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE↗

DOE OSTI2020

Heavy flavor dynamics across system size at the LHC

One of the fundamental signatures of the Quark Gluon Plasma has been the suppression of heavy flavor (specifically D mesons), which has been measured via the nuclear modification factor, R AA and azimuthal anisotropies, v n , in large systems. However, multiple competing models can reproduce the same data for R AA to v n . In this talk we break down the competing effects that conspire together to successfully reproduce R AA and v n in experimental data using Trento+v-USPhydro+DAB-MOD. Then using our best fit model we make predictions for R AA and v n across system size for 208 PbPb, 129 XeXe, 40 ArAr, and 16 OO collisions. We find that 0–10% centrality has a non-trivial interplay between the system size and eccentricities such that system size effects are masked in v 2 whereas in 30–50% centrality the eccentricities are approximately constant across system size and, therefore, is a better centrality class to study D meson dynamics across system size.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Wind Farm Wakes and Farm-to-Farm Interactions: Lidar and Wind Tunnel Tests

Recent experimental and numerical evidence has shown that the cumulative wake generated from the overlapping of multiple wakes within a wind farm could reduce power performance and enhance fatigue loads of wind turbines installed in neighboring downstream wind farms and may also extend up to distances one order of magnitude larger than those typically considered for intra-farm wake interactions. Similar to individual wind turbine wakes, wind farm wakes have a velocity deficit and added turbulence intensity, both affected by the turbine rotor thrust forces and the incoming turbulence intensity. Therefore, the evolution of wind farm wakes will vary for different operational and atmospheric conditions. In this paper, lidar measurements collected during the American WAKE experimeNt (AWAKEN) and wind tunnel tests of wind farms reproduced by porous disks are leveraged to investigate wind farm wakes.

17 WIND ENERGY↗

DOE OSTI2026

Industrial deposition of wavelength-shifting films for liquid argon photon detection systems

The Deep Underground Neutrino Experiment (DUNE) Phase-II Far Detector is considering an approximately 2000 m 2 photon detection system to achieve a target mean light yield of 180 PE/MeV. Meeting this requirement demands scalable, cost-effective, and high-quality wavelength-shifter (WLS) coatings capable of converting 127 nm liquid-argon scintillation light into visible photons with controlled and reproducible optical performance. We report on the successful realization of an industrial physical vapor deposition (PVD) process forp-terphenyl (pTP) coatings, adapted from vacuum deposition techniques developed for OLED display manufacturing, to produce uniform WLS layers on large-area inorganic substrates, a task traditionally challenged by adhesion and uniformity issues at organic-inorganic interfaces. Surface characterization by profilometry and spectroscopic measurements demonstrates edge-region thickness variation below 10% and emission spectra consistent with high-quality pTP reference samples. The industrial process demonstrates reproducibility, scalability, and significantly reduced production time compared to laboratory-based methods, while maintaining optical characteristics consistent with established pTP reference samples. These results establish a viable pathway for mass production of high-performance pTP coatings for DUNE FD3 and future neutrino experiments, from a coating manufacturing and process standpoint. Detector-level performance validation, including quantitative VUV conversion efficiency measurements at 127 nm, is identified as future work.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

The Thomas theorem and the Efimov States within a generalized Bohr model

Using a generalized Bohr model and the hyper-spherical formalism for a three-body system, we derive the universal energy functions. We apply our model to 4He atom where only the Coulomb potential is dominant and find that its binding energy is well reproduced (less than 5.5% difference). Later, we focus on the equal mass three-body systems and derive the Thomas theorem assuming a simple interaction depending on the range of the potential. We discuss the conditions for which an unbound two-body system produces a bound three-body system and apply our model to 4He and triton atoms as well as to the triton nucleus. Using their scattering lengths and effective ranges, we are able to reproduce the two-body or the three-body binding energies (less than 5% difference) with only one parameter fitted. Prediction for excited (Efimov) levels are also given and in particular we demonstrate that for some hyper-angles two equal minima appear which indicate a phase (shape) transition similar to the Landau's theory of phase transition. We suggest that the observed excited levels in two different experiments for the triton nucleus are indeed Efimov levels and there may be more surprises.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Best practices for photovoltaic performance loss rate calculations

Abstract The performance loss rate (PLR) is a vital parameter for the time-dependent assessment of photovoltaic (PV) system performance and health state. Although this metric can be calculated in a relatively straightforward manner, it is challenging to achieve accurate and reproducible results with low uncertainty. Furthermore, the temporal evolution of PV system performance is usually nonlinear, but in many cases a linear evaluation is preferred as it simplifies the assessment and it is easier to evaluate. As such, the search for a robust and reproducible calculation methodology providing reliable linear PLR values across different types of systems and conditions has been the focus of many research activities in recent years. In this paper, the determination of PV system PLR using different pipelines and approaches is critically evaluated and recommendations for best practices are given. As nonlinear PLR assessments are fairly new, there is no consent on how to calculate reliable values. Several promising nonlinear approaches have been developed recently and are presented as tools to evaluate the PV system performance in great detail. Furthermore, challenges are discussed with respect to the PLR calculation but also opportunities for differentiating individual performance losses from a generic PLR value having the potential of enabling actionable insights for maintenance.

14 SOLAR ENERGY↗

DOE OSTI2020

An open-source high-content analysis workflow for CFTR function measurements using the forskolin-induced swelling assay

Abstract Motivation The forskolin-induced swelling (FIS) assay has become the preferential assay to predict the efficacy of approved and investigational CFTR-modulating drugs for individuals with cystic fibrosis (CF). Currently, no standardized quantification method of FIS data exists thereby hampering inter-laboratory reproducibility. Results We developed a complete open-source workflow for standardized high-content analysis of CFTR function measurements in intestinal organoids using raw microscopy images as input. The workflow includes tools for (i) file and metadata handling; (ii) image quantification and (iii) statistical analysis. Our workflow reproduced results generated by published proprietary analysis protocols and enables standardized CFTR function measurements in CF organoids. Availability and implementation All workflow components are open-source and freely available: the htmrenamer R package for file handling https://github.com/hmbotelho/htmrenamer; CellProfiler and ImageJ analysis scripts/pipelines https://github.com/hmbotelho/FIS_image_analysis; the Organoid Analyst application for statistical analysis https://github.com/hmbotelho/organoid_analyst; detailed usage instructions and a demonstration dataset https://github.com/hmbotelho/FIS_analysis. Distributed under GPL v3.0. Supplementary information Supplementary data are available at Bioinformatics online.

Hagemeijer, Marne C.↗