Search NASA⌕ Search

SEARCH · Search NASA

Results for “Module scale model”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Asymmetrically Varying Guide Field During Magnetic Reconnection: Particle-In-Cell Simulations

Using fully kinetic particle-in-cell modeling, we investigate how magnetic reconnection responds to a varying guide field in one of the inflow regions. We find that the reconnection rate varies significantly when the orientation of the magnetic field changes between being strictly antiparallel and having a guide field. These variations are fairly consistent with the scaling relation for asymmetric reconnection developed by Cassak and Shay (2007). However, the rate is also found to be nonlinearly modulated by changes in the ion inflow velocity. The spatio-temporal change in the inflow velocity arises as the magnetic forces reconfigure to regions of different magnetic field strengths. The variations in the inflow magnetic field configuration allow for different gradients in the magnetic field, leading to asymmetries in the magnetic tension force. By momentum conservation, this facilitates asymmetries in the inflow velocity, which in turn affects the flux transport into the reconnection site. The outflow is found to be less laminar when the inflow varies, and various signatures of the inflow variations are identified in the outflow.

Susanne F Spinnangr↗
DOE OSTI2022

A new large-scale suspended sediment model and its application over the United States

Abstract. Suspended sediment plays a vital role in the regional and global cycling of carbon and nutrients by carrying carbon and nutrients from headwaters into lowland rivers and the oceans. Sediment transport through river systems is often fundamentally modified by human activities such as reservoir management. However, a physically based representation of sediment transport is still missing in most existing Earth system models (ESMs), which are essential tools for modeling and predicting Earth system changes. Here, we introduce a multi-process river sediment module for ESMs, which includes (1) hillslope soil erosion and sediment discharge into streams, (2) sediment transport processes through river networks, (3) reservoir operation based on the inflows from upstream areas and water demand from downstream areas, and (4) sediment trapping by reservoirs. All model parameters are estimated a priori without calibration. We apply this new sediment modeling framework to the contiguous United States and validate it against historical observations of monthly streamflow and sediment discharges at 35 river gauges. The model reasonably well captures the long-term balance and seasonal variations of suspended sediment in large river systems. Furthermore, our model results show that suspended sediment discharge in managed rivers is affected more by reservoirs' direct trapping of sediment particles than by their flow regulation. This new sediment module enables future modeling of the transportation and transformation of carbon and nutrients carried by the fine sediment along the river–ocean continuum to close the global carbon and nutrient cycles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Development and Rigorous Validation of Modules for Predicting Organic Aerosol CCN Activity from Molecular Properties

The indirect interaction of anthropogenic aerosols with clouds (aerosol indirect forcing) represents one of the main uncertainties in climate projections, and reducing this uncertainty requires improved understanding of the aerosol life cycle, including atmospheric aging and the ability to nucleate cloud droplets or crystals. A key parameter that links the chemical properties of aerosol to its nucleating properties, and ultimately to aerosol indirect forcing, is its hygroscopicity (ability to uptake water). One way to compare and systematize the ability of various chemical compounds or particle types to nucleate cloud drops is via the hygroscopicity parameter κ, which can be used directly in global models to simulate aerosol indirect forcing. In this project we further investigated the relationships between the chemical composition of organic aerosol particles, which are known to be important components of atmospheric particulate matter, and their ability to serve as cloud condensation nuclei (CCN). In particular, we isolated the role of viscosity on κ by applying methods similar to those we used previously in our DOE/ASR research. Our motivation for the focus on viscosity was through its link to particle phase state. Particle phase affects water uptake on timescales relevant to the laboratory and also likely relevant to atmospheric processes, with reduced apparent hygroscopicities for semisolid and solid particles, compared with liquid particles. Water uptake itself affects viscosity, an important feedback that we explored in this project. The objectives of the project were as follows: (1) measure the effects of molecular structure, temperature, and relative humidity on the viscosity of organic aerosol formed via controlled laboratory experiments, (2) measure the effects of organic aerosol viscosity on water uptake and CCN activity above and below water saturation, (3) further develop the κ parameterization of organic aerosol hygroscopicity to include viscosity/phase, and (4) experimentally validate the κ parameterization for organic aerosol with chemical composition ranging from relatively simple to complex. It was proposed that particle viscosity exerts a kinetic control on equilibrium water uptake. High viscosity was expected to delay water uptake based on diffusional limitations. This would manifest itself in reduced observed CCN activity or κ. Understanding the interplay between water content and viscosity (composition control on particle viscosity) and between viscosity and equilibrium water content (diffusional control on water uptake) were the overall science objectives of this grant. Targeted experiments were carried out to test the kinetic control hypothesis. These experiments included studies that explicitly measured the equilibration time scale of water uptake, and experiments that tested the role of viscous particle coatings on the observed CCN activity. From the combined work, it is concluded that mass transfer limitations by glassy organic particles are unlikely to affect cloud droplet activation near laboratory temperatures. This implies that viscosity does not need to be included in parameterizations that link particle chemical composition and κ. A parameterization to predict κ through the UNIFAC group contribution method was developed with previous ASR funding. As planned, this computational module was further developed and validated here through targeted laboratory studies. The computational module has been made available as free software under the GNU public license. A database of κ-values to aid model development has been compiled and is available online via an open access model. In conducting this work, significant effort was made to understand how aerosol composition, including the associated aerosol water content, affects viscosity. This included a study detailing the influence of functional group composition on organic aerosol viscosity. We identified a strong link between volatility and viscosity. Functional groups responsible for low vapor pressure also promote high viscosity. Specifically, carboxylic acid and hydroxyl groups were identified as major contributors to high viscosity. Several studies were carried out to characterize the link between composition and the glass-to-liquid phase transition. Compositions tested include models of organic/inorganic particles and models of complex organic aerosol formed through secondary processes. A simple computational model to predict the amorphous phase diagram was developed and tested. Excellent agreement between model and measurements was obtained for organic/inorganic mixtures. The computational model is simple enough to be included in global model simulations. However, constraints of input parameters such as the Gordon-Taylor constant for binary mixtures for a wide range of complex mixtures have not yet been fully evaluated. Combined, the performed work provides new constraints on the amorphous phase state of organic aerosols.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Nonlinear models for coupling the effects of stimulated Raman scattering to inertial confinement fusion codes

Laser plasma instabilities (LPI) reduce driver-target coupling, alter implosion symmetry, and therefore can fundamentally limit fusion performance in inertial confinement fusion (ICF). Developing a predictive modeling capability for LPI effects can critically advance the success of the field. We perform vector particle-in-cell simulations of multi-speckled laser beams undergoing stimulated Raman scattering (SRS) at various densities and intensities relevant to mainly indirectly driven and a subset of parameter space for directly driven ICF systems, focusing on the regimes with intensities above the onset of electron trapping. Based on the wavenumber of the SRS daughter electron plasma wave, we identify several regions with underpinning SRS saturation physics: the electron-trapping dominated region with intermediate kλD values, the strong Landau damping region at larger kλD values, and the region where the Langmuir decay instability arises at lower kλD values. We develop a nonlinear SRS reflectivity model that features the base trapping-dominated scaling of (kλD)−4 and its modifications. Electron trapping deforms the initialized electron distribution functions, and we have developed a new δf-Gaussian-mixture algorithm for an accurate characterization of the trapped hot electron population. With this SRS hot electron description, we construct a nonlinear energy deposition model and a hot electron source model—based on a modified Manley–Rowe relation—suitable for including SRS effects as a sub-grid module in a high-fidelity ICF design code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Cloud Atlas: Unraveling the Vertical Cloud Structure with the Time-series Spectrophotometry of an Unusually Red Brown Dwarf

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352—an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR 8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the HST and the Spitzer time-series observations, as well as the time-averaged spectra of WISE 0047. Additionally, in the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/CH 4 chemistry, our models also reproduce the redder $J-{K}_{{\rm{s}}}$ color of WISE 0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1963

Characteristics of a Lunar Landing Configuration Having Various Multiple-Leg Landing-Gear Arrangements

An experimental investigation has been made of some lunar-landing characteristics of a 1/6-scale dynamic model of a landing module having multiple-leg landing-gear systems. Symmetric four-point and five-point systems and an asymmetric four-point system were investigated. The landing-gear legs were inverted tripod arrangements having a telescoping main strut which incorporated a yielding-metal strap for energy dissipation, hinged V-struts, and circular pads. The landing tests were made by launching a free model onto an impenetrable hard surface (concrete) and onto a powdered-pumice overlay of various depths. Landing motion and acceleration data were obtained for a range of touchdown speeds, touchdown speeds, touch attitudes, and landing-surface conditions. Symmetric four-point and five-point systems and an Maximum normal acceleration experienced at the module center of gravity during landings on hard surface or pumice was 2g (full-scale lunar value in terms of earth's gravity) over a wide range of touchdown conditions. Maximum angular acceleration experienced was 12-1/2 radians/sec(exp 2) and maximum longitudinal acceleration was 1-3/4 g. The module was very stable with all gear configurations during landings on hard surface (coefficient of friction, microns=0.4) at all conditions tested. Some overturn instability occurred during landings on powdered pumice (microns=0.7 to 1.0) depending upon flight path, pitch and yaw attitude, depth of pumice, surface topography, and landing-gear configuration. The effect of stability of roll attitude for the limited amount of roll-attitude landing data obtained was insignificant. Compared with the four-point system, the five-point system with equal maximum gear radius increased landing stability slightly and improved the static stability for subsequent lunar launch. A considerable increase in landing stability in the direction of motion was obtained with an asymmetric four-point gear having two pads offset to increase gear radius by 33 percent in the direction of horizontal flight.

Blanchard, Ulysse J.↗
NASA NTRS2014

Advanced Mirror & Modelling Technology Development

The 2020 Decadal technology survey is starting in 2018. Technology on the shelf at that time will help guide selection to future low risk and low cost missions. The Advanced Mirror Technology Development (AMTD) team has identified development priorities based on science goals and engineering requirements for Ultraviolet Optical near‐Infrared (UVOIR) missions in order to contribute to the selection process. One key development identified was lightweight mirror fabrication and testing. A monolithic, stacked, deep core mirror was fused and replicated twice to achieve the desired radius of curvature. It was subsequently successfully polished and tested. A recently awarded second phase to the AMTD project will develop larger mirrors to demonstrate the lateral scaling of the deep core mirror technology. Another key development was rapid modeling for the mirror. One model focused on generating optical and structural model results in minutes instead of months. Many variables could be accounted for regarding the core, face plate and back structure details. A portion of a spacecraft model was also developed. The spacecraft model incorporated direct integration to transform optical path difference to Point Spread Function (PSF) and between PSF to modulation transfer function. The second phase to the project will take the results of the rapid mirror modeler and integrate them into the rapid spacecraft modeler.

Effinger, Michael↗
NASA NTRS1998

Topography of the South Polar Cap and Layered Deposits of Mars: Viking Stereo Grametry at Regional and Local Scales

Layered deposits and residual polar caps on Mars may record the deposition of ice and sediment modulated by periodic climate change. Topographic information relating to layer thicknesses, erosional processes, and formation of dark spirals within these deposits has been sparce or unreliable until the arrival of MOLA in orbit in September 1997. To assist in evaluating these terrains prior to launch and to assess formation and erosion processes in the polar deposits, we have assembled Viking stereo mosaics of the region and have produced the first reliable DEM models of the south polar deposits using automated stereogrammetry tools. Here we report our preliminary topographic results, pending final image pointing updates. The maximum total thickness of the layered deposits in the south polar region is 2.5 km. The thick layered deposits consist of a series of megaterraces. Each terrace is several tens of kilometers wide and is flat or slopes very gently toward the pole. These terraces step downward from a central plateau near the south pole. Terraces are bounded by relatively steep scarps 100-500 meters high that face toward the equator. These scarps correspond to the pattern of dark spirals observed within the residual cap in southern summer, and are interpreted as ice or frost-free surfaces warmed by solar insolation. Several tongue-shaped troughs, with rounded cirquelike heads, are observed near the margins of the deposit. These troughs are 300-600 meters in deep and may be similar to troughs observed in the northern polar deposit.

Schenk, P.↗
DOE OSTI2024

Advances in statistical methods for cancer surveillance research: an age-period-cohort perspective

Background: Analysis of Lexis diagrams (population-based cancer incidence and mortality rates indexed by age group and calendar period) requires specialized statistical methods. However, existing methods have limitations that can now be overcome using new approaches. Methods: We assembled a “toolbox” of novel methods to identify trends and patterns by age group, calendar period, and birth cohort. We evaluated operating characteristics across 152 cancer incidence Lexis diagrams compiled from United States (US) Surveillance, Epidemiology and End Results Program data for 21 leading cancers in men and women in four race and ethnicity groups (the “cancer incidence panel”). Results: Nonparametric singular values adaptive kernel filtration (SIFT) decreased the estimated root mean squared error by 90% across the cancer incidence panel. A novel method for semi-parametric age-period-cohort analysis (SAGE) provided optimally smoothed estimates of age-period-cohort (APC) estimable functions and stabilized estimates of lack-of-fit (LOF). SAGE identified statistically significant birth cohort effects across the entire cancer panel; LOF had little impact. As illustrated for colon cancer, newly developed methods for comparative age-period-cohort analysis can elucidate cancer heterogeneity that would otherwise be difficult or impossible to discern using standard methods. Conclusions: Cancer surveillance researchers can now identify fine-scale temporal signals with unprecedented accuracy and elucidate cancer heterogeneity with unprecedented specificity. Birth cohort effects are ubiquitous modulators of cancer incidence in the US. The novel methods described here can advance cancer surveillance research.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2020

State of the Profession Considerations: NASA Langley Research Center Capabilities / Technologies for Autonomous In-Space Assembly and Modular Persistent Assets

Successfully implementing OSAM into next generation revolutionary observatories requires integrating expertise and technologies in modular space structures, assembly operations, autonomy, and modeling/simulation. LaRC OSAM technologies/capabilities have been presented to inform the Planetary Science and Astrobiology Decadal Survey community of the robust and mature existing capability to support an OSAM based architecture for their next observatory. LaRC Structures and Assembly capabilities enable; a modular telescope architecture, high-performance structural modules, and robotic assembly techniques. LaRC Autonomy capabilities ensure that the robotic assembly will be accomplished in a safe and robust manner and only require humans in a supervisory role. The LaRC toolbox of Modeling and Simulation capabilities that is calibrated using module-level ground testing, will ensure that the performance of the fully assembled observatory, a very large zero-g system that will never be assembled/tested in a gravity environment, meets all performance requirements when it enters into service. Integrating all three LaRC capabilities and including embedded metrology, will enable servicing, repair, instrument upgrades (and/or replacement) while ensuring a very long lifetime for the observatory and providing a return-on-investment that is substantially greater than the initial cost. Further confidence will be achieved as OSAM technologies are validated in a new LaRC OSAM laboratory that allows large-scale collaborative testing of modular hardware, simulation software and algorithms, and autonomous agents.

Large space structures↗
DOE OSTI2021

A Comprehensive Method to Mitigate Forced Oscillations in Large Interconnected Power Grids

Multiple severe forced oscillation events have recently been observed across the North American interconnections and around the world. These forced oscillations have caused power swings, limited power transfer capability, damaged equipment, and persisted indefinitely until the driving source was located and removed. This paper proposes a comprehensive method to mitigate forced oscillations. Once a forced oscillation is detected, a new source location algorithm can be used to locate the source based on the oscillation mode angle without requiring system topology information. If the source cannot be quickly located and removed, a control strategy can be activated to modulate the active power of utility-scale inverterbased battery energy storage systems (BESSs) to reduce the energy of forced oscillations to a safe level and allow sufficient time for locating the exact source. The proposed source location algorithm is validated using measurements collected during the January 11, 2019 forced oscillation event in North America and other actual grid events, while the proposed control strategy is verified using the Eastern Interconnection dynamic model under the replicated January 11, 2019 forced oscillation event. The simulation results demonstrated that the proposed source location algorithm can accurately identify the forced oscillation source, and the proposed control strategy can significantly reduce forced oscillation energy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

U.S. Solar Photovoltaic System and Energy Storage Cost Benchmark: Q1 2021 [Slides]

Our benchmarking method includes bottom-up accounting for all necessary system and projectdevelopment costs incurred when installing residential, commercial, and utility-scale systems, and it models the Q1 2021 costs for such systems, excluding any previous supply agreements or contracts. In general, we attempt to model the typical installation techniques and business operations from an installed-cost perspective, and our benchmarks are national averages. The residential PV-only benchmark and the commercial rooftop PV-only benchmark average costs by inverter type (string inverters, string inverters with direct current [DC] optimizers, and microinverters), weighted by inverter market share. The residential PV-only benchmark is further averaged across small installer and national integrator business models, weighted by market share. All benchmarks include variations—accounting for the differences in size, equipment, and operational use (particularly for storage)—that are currently available in the marketplace. All benchmarks assume nonunionized construction labor; residential and commercial PV systems predominantly use nonunionized labor, and the type of labor required for utility-scale PV systems depends heavily on the development process. All benchmarks assume the use of monofacial monocrystalline silicon PV modules. Benchmarks using cadmium telluride or bifacial modules could result in significantly different results. The data in this annual benchmark report inform the formulation of and track progress toward the U.S. Department of Energy (DOE) Solar Energy Technologies Office’s Government Performance and Reporting Act cost targets.

14 SOLAR ENERGY↗
NASA NTRS1980

Ocean wave-radar modulation transfer functions from the West Coast experiment

Short gravity-capillary waves, the equilibrium, or the steady state excitations of the ocean surface are modulated by longer ocean waves. These short waves are the predominant microwave scatterers on the ocean surface under many viewing conditions so that the modulation is readily measured with CW Doppler radar used as a two-scale wave probe. Modulation transfer functions (the ratio of the cross spectrum of the line-of-sight orbital speed and backscattered microwave power to the autospectrum of the line-of-sight orbital speed) were measured at 9.375 and 1.5 GHz (Bragg wavelengths of 2.3 and 13 cm) for winds up to 10 m/s and ocean wave periods from 2-18 s. The measurements were compared with the relaxation-time model; the principal result is that a source of modulation other than straining by the horizontal component of orbital speed, possibly the wave-induced airflow, is responsible for most of the modulation by waves of typical ocean wave period (10 s). The modulations are large; for unit coherence, spectra of radar images of deep-water waves should be proportional to the quotient of the slope spectra of the ocean waves by the ocean wave frequency.

Wright, J. W.↗
NASA NTRS1987

Galactic cosmic rays and anomalous components in the heliosphere

The progress made during 1983-1986 by U.S scientists on modulation of Galactic cosmic rays by the sun and on the origin, acceleration, and modulation of the anomalous components in the heliosphere is reviewed. In these years, simultaneous observations of modulation, covering a radial range from 1 to more than 30 AU and a heliospheric latitude range from 0 to greater than 20 deg, made it possible to investigate directly the dynamics and spatial scale of modulation during both the increasing and declining phases of solar activity and to make in situ correlations between changes in the interplanetary medium and changes in the modulation of cosmic rays. Particular attention is given to the use of the observations related to time-dependent modulation effects, curvature drifts, and anomalous components, as well as of the measurements of the cosmic-ray radial gradients and latitudinal intensity variations, in the modeling of solar modulation effects.

Mckibben, R. B.↗
DOE OSTI2024

Preliminary searches for spin-dependent interactions using sidebands of nuclear spin-precession signals

Various theories beyond the Standard Model predict new particles with masses in the sub-eV range with very weak couplings to ordinary matter. A new P-odd and T-odd interaction between polarized and unpolarized nucleons proportional to s⃗⋅r̂ is one such possibility, where r⃗=rr̂ is the spatial vector connecting the nucleons, and s⃗ is the spin of the polarized nucleon. Such an interaction involving a scalar coupling gsN at one vertex and a pseudoscalar coupling gpn at the polarized nucleon vertex can be induced by the exchange of spin-0 pseudoscalar bosons. We describe a new technique to search for interactions of this form and present the first measurements of this type. We show that future improvements to this technique can improve the laboratory upper bound on the product gsNgpn by two orders of magnitude for interaction ranges at the 100 micron scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Phase Behavior of Mixed Polymer Brushes Grown from Ultrathin Coatings

Experimental validation of the predicted melt phase behavior of A/B mixed or binary brush on planar substrate is presented using poly(methyl methacrylate) (A)/ polystyrene (B) (PMMA/PS) with equal number of A/B chains as an example. Specifically, well-defined mixed A/B brushes are synthesized using a single component inimer coating to achieve high grafting density (0.9 chains/ nm 2 ), uniformity of grafting sites, and predictable chain length. The inimer coating is a copolymer ofnitroxide-mediated polymerization (NMP) inimer, atom transfer radical polymerization (ATRP) inimer, styrene, and glycidyl methacrylate (GMA). Cross-linking of the film provides the required stability to probe the melt morphology. Our studies show that even with equal grafting density of the A and B the morphology can be modulated by varying the length of B chains while keeping that of A fixed. We show the very first example of transition of self-assembled structures from disorder to cylinder to ripple phase at sub-30 nm length scale on a planar surface by thermal annealing of mixed brushes. These results are supported by a new phase diagram established through Monte-Carlo simulation using a coarse-grained particle-based model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A multi-scale pipeline linking drug transcriptomics with pharmacokinetics predicts in vivo interactions of tuberculosis drugs

Tuberculosis (TB) is the deadliest infectious disease worldwide. The design of new treatments for TB is hindered by the large number of candidate drugs, drug combinations, dosing choices, and complex pharmaco-kinetics/dynamics (PK/PD). Here we study the interplay of these factors in designing combination therapies by linking a machine-learning model, INDIGO-MTB, which predicts in vitro drug interactions using drug transcriptomics, with a multi-scale model of drug PK/PD and pathogen-immune interactions called GranSim. We calculate an in vivo drug interaction score (iDIS) from dynamics of drug diffusion, spatial distribution, and activity within lesions against various pathogen sub-populations. The iDIS of drug regimens evaluated against non-replicating bacteria significantly correlates with efficacy metrics from clinical trials. Our approach identifies mechanisms that can amplify synergistic or mitigate antagonistic drug interactions in vivo by modulating the relative distribution of drugs. Our mechanistic framework enables efficient evaluation of in vivo drug interactions and optimization of combination therapies.

59 BASIC BIOLOGICAL SCIENCES↗