Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic parameter 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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Enhancements to the RAVEN code in FY16

The RAVEN code has been under development at the Idaho National Laboratory since 2012. Its main goal is to create a multi-purpose platform for the deploying of all the capabilities needed for Probabilistic Risk Assessment, uncertainty quantification, data mining analysis and optimization studies. RAVEN has currently reached a good level of maturity in terms of deployed state-of-art and advanced capabilities. The main subject of this report is to show the activities that have been recently accomplished: • Implementation of ensemble modeling for time-series, and • initial implementation of model validation for surrogate models, and • advanced visualization capability for topology based data analysis The development of ensemble modeling for time-series has been performed in order to begin tackling the needs of those RISMC applications that need to communicate 1-Dimensional information (e.g. power histories, etc.) among different models. In this document the implementation details and an application example is reported. The second subject of this report is about the initial development of methods, within the RAVEN framework, to assess the validity of the predictive capabilities of surrogate models. Indeed, after the construction of a surrogate tight to a certain physical model, it is crucial to assess the goodness of its representation, in order to be confident with its prediction. In this initial activity, a cross-validation technique has been employed. This report will highlight the implementation details and proof its correct implementation by an application example. The final subject of this report is about the implementation of advanced visualization capability in RAVEN, for interactive data analysis. Indeed, RAVEN offers several post-processing capabilities that can structurally decompose data extracted from experimental results offering both data clustering/partitioning and dimensionality reduction techniques. A disadvantage of the workflow available in RAVEN is that it treats these as black box operations and the user is expected to know specific information about their data including the number of partitions to expect in some cases or the “correct” parameter settings for a particular algorithm. In order to overcome this limitation, it has been added an interactive user interface that can be run in RAVEN to explore different parameter settings on the fly specifically for the topological post-processor and allows the user to explore the data interactively. The design of this user interface is generalizable to other algorithms and can be used as a model to generate more dynamic and user-friendly visualization capabilities within RAVEN.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Mapping Thermal Conductivity at the Atomic Scale: A Step toward the Thermal Design of Materials

We describe a spatial decomposition of the thermal conductivity, termed site-projected thermal conductivity, which quantifies the thermal conduction activity at each atomic site—a critical parameter for the thermal design of materials. The method is based on the Green–Kubo formula and the harmonic approximation and requires the force-constant and dynamical matrices, as well as a relaxed structural model. Throughout the study, it uses high quality models previously tested and compared to many experiments. It discusses the method and underlying approximations for amorphous silicon, carries the detailed analysis for amorphous silicon, and then examines an amorphous-crystal silicon interface and representative carbon materials. This study identifies the sites and local structures that reduce heat transport, and quantify these (estimate the spatial range) over which these “thermal defects” are effective. It identifies filamentary structures in the amorphous silicon network which impact heat transport and electronic structure (the Urbach edge) and electronic transport.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An efficient surrogate model of secondary electron formation and evolution

This work extends the adjoint-deep learning framework for runaway electron (RE) evolution, developed by McDevitt et al. [Phys. Plasmas 32, 042503 (2025)], to account for large-angle collisions. By incorporating large-angle collisions, the framework allows the avalanche of REs to be captured, an essential component of RE dynamics. This extension is accomplished by using a Rosenbluth–Putvinski approximation to estimate the distribution of secondary electrons generated by large-angle collisions. By evolving both the primary and multiple generations of secondary electrons, the present formulation can capture both the detailed temporal evolution of a RE population beginning from an arbitrary initial momentum space distribution, along with providing approximations to the saturated growth and decay rates of the RE population. Predictions of the adjoint-deep learning framework are verified against a traditional RE solver, with good agreement present across a broad range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Triggering collapse of the presolar dense cloud core and injecting short-lived radioisotopes with a shock wave. III. Rotating three-dimensional cloud cores

A key test of the supernova triggering and injection hypothesis for the origin of the solar system's short-lived radioisotopes is to reproduce the inferred initial abundances of these isotopes. We present here the most detailed models to date of the shock wave triggering and injection process, where shock waves with varied properties strike fully three-dimensional, rotating, dense cloud cores. The models are calculated with the FLASH adaptive mesh hydrodynamics code. Three different outcomes can result: triggered collapse leading to fragmentation into a multiple protostar system; triggered collapse leading to a single protostar embedded in a protostellar disk; or failure to undergo dynamic collapse. Shock wave material is injected into the collapsing clouds through Rayleigh-Taylor fingers, resulting in initially inhomogeneous distributions in the protostars and protostellar disks. Cloud rotation about an axis aligned with the shock propagation direction does not increase the injection efficiency appreciably, as the shock parameters were chosen to be optimal for injection even in the absence of rotation. For a shock wave from a core-collapse supernova, the dilution factors for supernova material are in the range of ∼10{sup –4} to ∼3 × 10{sup –4}, in agreement with recent laboratory estimates of the required amount of dilution for {sup 60}Fe and {sup 26}Al. We conclude that a type II supernova remains as a promising candidate for synthesizing the solar system's short-lived radioisotopes shortly before their injection into the presolar cloud core by the supernova's remnant shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Toward Development of a Stochastic Wake Model: Validation Using LES and Turbine Loads

Wind turbines within an array do not experience free-stream undisturbed flow fields. Rather, the flow fields on internal turbines are influenced by wakes generated by upwind unit and exhibit different dynamic characteristics relative to the free stream. The International Electrotechnical Commission (IEC) standard 61400-1 for the design of wind turbines only considers a deterministic wake model for the design of a wind plant. This study is focused on the development of a stochastic model for waked wind fields. First, high-fidelity physics-based waked wind velocity fields are generated using Large-Eddy Simulation (LES). Stochastic characteristics of these LES waked wind velocity field, including mean and turbulence components, are analyzed. Wake-related mean and turbulence field-related parameters are then estimated for use with a stochastic model, using Multivariate Multiple Linear Regression (MMLR) with the LES data. To validate the simulated wind fields based on the stochastic model, wind turbine tower and blade loads are generated using aeroelastic simulation for utility-scale wind turbine models and compared with those based directly on the LES inflow. The study’s overall objective is to offer efficient and validated stochastic approaches that are computationally tractable for assessing the performance and loads of turbines operating in wakes.

17 WIND ENERGY↗
DOE OSTI2018

Millisecond Pulsar Origin of the Galactic Center Excess and Extended Gamma-Ray Emission from Andromeda: A Closer Look

A new measurement of a spatially extended gamma-ray signal from the center of the Andromeda galaxy (M31) has recently been published by the Fermi -LAT collaboration, reporting that the emission broadly resembles the so-called Galactic center excess (GCE) of the Milky Way (MW). The weight of the evidence is steadily accumulating on a millisecond pulsar (MSPs) origin for the GCE. These elements prompt us to compare these observations with what is, perhaps, the simplest model for an MSP population, which is solely obtained by rescaling of the MSP luminosity function that is determined in the local MW disk via the respective stellar mass of the systems. Remarkably, we find that without free fitting parameters, this model can account for both the energetics and the morphology of the GCE within uncertainties. For M31, the estimated luminosity due to primordial MSPs is expected to only contribute about a quarter of the detected emission, although a stronger contribution cannot be excluded given the large uncertainties. If correct, the model predicts that the M31 disk emission due to MSPs is not far below the present upper bound. We also discuss additional refinements of this simple model. Using the correlation between globular cluster gamma-ray luminosity and stellar encounter rate, we gauge the dynamical MSP formation in the bulge. This component is expected to contribute to the GCE only at a level of ≲5%, it could affect the signal's morphology. In conclusion, we also comment on the limitations of our model and on future perspectives for improved diagnostics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

AI-Science for Performance Optimization and Diagnosis of Science Instrument Federations

Next generation of science workflows are expected to be executed over complex federations composed of supercomputers, science instruments, storage systems and networks, with new additions of the edge and cloud systems and services. The sheer complexity of these multi-domain federations makes it hard to manage them and optimize their performance, as small impedance mismatches (that can dynamically develop between systems) could drastically degrade the entire federation performance. Recent proliferation of Software Defined Everything (SDX) technologies combined with containerization frameworks provide custom instruments that can monitor and collect critical measurements at various levels to support diagnoses and performance optimization; but their data too enormous for human operators and analysts to process and generate decisions. Machine Learning (ML) methods that extract critical parameters, relationships and trends from the data offer general solutions. Artificial Intelligence (AI) and ML methods must be custom-developed for these problems based on solid, rigorous foundations, since black-box approaches are often ineffective and unsound.We propose to develop comprehensive AI-Science for the performance of science federations to (i) monitor and control storage, networks, experiments, and computing systems across multiple domains via softwarization layers, at speeds and scales orders of magnitude superior to current practice, (ii) optimally realize and orchestrate complex workflows with high performance by using dynamic state and performance estimation methods, and (iii) aggregate measurements across sites and time to develop infrastructure-level profiles, optimizations and diagnoses using AI-Science based on foundational principles from ML, game theory, and information fusion areas.

Rao, Nageswara↗
DOE OSTI2023

A phased workflow to define permit‐ready locations for large volume CO 2 injection and storage

Abstract To‐date, only two UIC Class VI permits have been issued by the US Environmental Protection Agency. We illustrate a four‐phase workflow to first identify regional storage resources and then down‐select sites to yield permit‐ready locations that can accept and store large volumes of CO 2 . Specific permit requirements should guide objectives and define deliverables of respective workflow phases. In the first phase we used available regional data and screened structure and injection zones to locate resources that match CO 2 volumes planned to be captured. Available data were also used to assess presence and depth of usable groundwater, the key resource being protected via permitting. We then used advanced, closed‐form, analytical solutions (EASiTool) to estimate CO 2 injectivity into each hydrologically connected injection compartment. In the second phase we acquired and conditioned additional wireline logs and leased available seismic datasets. We interpreted the depositional systems from wireline well‐log character and mapped sandbody geometry to interpolate injection and confining‐zone distribution. Using available data, we mapped faults and locations of freshwater and overpressure (or other capacity‐limiting geologic parameters) in more detail. In the third phase, we used the augmented geologic data to develop a static model for the selected area, extracted the areas of highest interest, and generated and ran dynamic (flow) models. In a fourth phase, we reduced major uncertainties identified in earlier phases. Our case study indicates that to complete preparation of a permit application requires (1) improved lithologic characterization information (thicknesses and horizontal and vertical connectivity) and (2) better definition of poorly defined local faults. © 2023 The Authors. Greenhouse Gases: Science and Technology published by Society of Chemical Industry and John Wiley & Sons Ltd.

58 GEOSCIENCES↗
DOE OSTI2020

Impact analysis of DERs on bulk power system stability through the parameterization of aggregated DER_a model for real feeders

With an ever increasing percentage of distributed energy resources (DERs) connected behind the meter in the distribution system, it is becoming increasingly important to equip transmission planners with the visibility of DER dynamic performance in distribution system. Not having visibility of the disconnection of DERs with the occurrence of transmission events, could result in an erroneous view of the stability of bulk power system. Here, in this paper, a parameterized aggregated model (DER_a) serves as a representation of the distribution-level dynamics of real residential feeders, which is used for analysis of bulk power system stability. The parameters are obtained by executing dynamic Monte Carlo simulations. Faults are then induced at the substation level causing the DERs to trip which subsequently enables the parameterization of the low and high voltage breakpoints (v l0 , v l1 , v h0 , and v h1 ) of the DER_a model’s partial voltage trip block. These parameters are then utilized to study the effectiveness of the DER_a model to represent the behavior of the aggregated DERs’ response and their impact on the bulk power system. The case study shows the ability of the positive-sequence DER_a model to provide an accurate estimation of DERs that are susceptible to trip due to 3-$\phi$ and 1-$\phi$ faults of transmission.

14 SOLAR ENERGY↗
DOE OSTI2023

Encoding nonlinear and unsteady aerodynamics of limit cycle oscillations using nonlinear sparse Bayesian learning

This article investigates the applicability of a recently proposed, nonlinear sparse Bayesian learning (NSBL) algorithm to identify and estimate the complex aerodynamics of limit cycle oscillations. NSBL provides a semi-analytical framework for determining the data-optimal sparse model nested within a (potentially) over-parameterized model. This is particularly relevant to nonlinear dynamical systems where modelling approaches involve the use of physics-based and data-driven components. In such cases, the data-driven components, where analytical descriptions of the physical processes are not readily available, are often prone to overfitting, meaning that the empirical aspects of these models will often involve the calibration of an unnecessarily large number of parameters. While an overparameterized model may fit the observed data well, such models may be inadequate for making predictions in regimes that are different from those wherein the data were recorded. In view of this, it is desirable to not only calibrate the model parameters, but also identify the optimal compromise between data fit and model complexity. In this article, we exhibit the optimal model discovery for an aeroelastic system wherein the structural dynamics are well-known and described by a differential equation model, coupled with a semi-empirical aerodynamic model for laminar separation flutter, resulting in low-amplitude limit cycle oscillations (LCO). To illustrate the performance of the algorithm, in this article, we use synthetic data and demonstrate the ability of the algorithm to correctly rediscover the optimal model and model parameters, given a known data-generating model. The synthetic data are generated from a forward simulation of a known differential equation model with parameters selected so as to mimic the dynamics observed in wind-tunnel experiments. Subsequently, we demonstrate the performance of the algorithm for model selection using noisy LCO data from wind tunnel experiments. As there is no ground truth available for the experimental data case, we provide a comparison between NSBL and Bayesian model selection to validate the results, and demonstrate the use of NSBL as an efficient alternative to traditional methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. Finally, this versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Using field observations to inform thermal hydrology models of permafrost dynamics with ATS (v0.83)

Climate change is profoundly transforming the carbon-rich Arctic tundra landscape, potentially moving it from a carbon sink to a carbon source by increasing the thickness of soil that thaws on a seasonal basis. Thus, the modeling capability and precise parameterizations of the physical characteristics needed to estimate projected active layer thickness (ALT) are limited in Earth system models (ESMs). In particular, discrepancies in spatial scale between field measurements and Earth system models challenge validation and parameterization of hydrothermal models. A recently developed surface–subsurface model for permafrost thermal hydrology, the Advanced Terrestrial Simulator (ATS), is used in combination with field measurements to achieve the goals of constructing a process-rich model based on plausible parameters and to identify fine-scale controls of ALT in ice-wedge polygon tundra in Barrow, Alaska. An iterative model refinement procedure that cycles between borehole temperature and snow cover measurements and simulations functions to evaluate and parameterize different model processes necessary to simulate freeze–thaw processes and ALT formation. After model refinement and calibration, reasonable matches between simulated and measured soil temperatures are obtained, with the largest errors occurring during early summer above ice wedges (e.g., troughs). The results suggest that properly constructed and calibrated one-dimensional thermal hydrology models have the potential to provide reasonable representation of the subsurface thermal response and can be used to infer model input parameters and process representations. The models for soil thermal conductivity and snow distribution were found to be the most sensitive process representations. However, information on lateral flow and snowpack evolution might be needed to constrain model representations of surface hydrology and snow depth.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Thermodynamics of Restricted Boltzmann Machines and Related Learning Dynamics

We investigate the thermodynamic properties of a restricted Boltzmann machine (RBM), a simple energy-based generative model used in the context of unsupervised learning. Assuming the information content of this model to be mainly reflected by the spectral properties of its weight matrix W, we try to make a realistic analysis by averaging over an appropriate statistical ensemble of RBMs. First, a phase diagram is derived. Otherwise similar to that of the Sherrington–Kirkpatrick (SK) model with ferromagnetic couplings, the RBM’s phase diagram presents a ferromagnetic phase which may or may not be of compositional type depending on the kurtosis of the distribution of the components of the singular vectors of W. Subsequently, the learning dynamics of the RBM is studied in the thermodynamic limit. A “typical” learning trajectory is shown to solve an effective dynamical equation, based on the aforementioned ensemble average and explicitly involving order parameters obtained from the thermodynamic analysis. In particular, this let us show how the evolution of the dominant singular values of W, and thus of the unstable modes, is driven by the input data. At the beginning of the training, in which the RBM is found to operate in the linear regime, the unstable modes reflect the dominant covariance modes of the data. In the non-linear regime, instead, the selected modes interact and eventually impose a matching of the order parameters to their empirical counterparts estimated from the data. Finally, we illustrate our considerations by performing experiments on both artificial and real data, showing in particular how the RBM operates in the ferromagnetic compositional phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Plant-wide modeling and techno-economic analysis of a direct non-oxidative methane dehydroaromatization process via conventional and microwave-assisted catalysis

Direct non-oxidative methane dehydroaromatization (DHA) process via conventional and microwave (MW)-assisted thermo-catalytic catalysis is studied. Rate models for methane DHA reactions, including the effect of catalyst deactivation, are developed by using the in-house experimental data. Model results for gas concentration profile and catalyst deactivation are in good agreement with the experimental data. This rate model is then used for the development of dynamic multi-scale, multi-physics commercial-scale reactor models. Total number of fixed bed reactors desired for a cyclic steady state process is estimated. Plant-wide models are then developed for conventional and MW-assisted processes for producing products of desired specifications. Techno-economic analysis of the methane DHA process is undertaken. Economics of these methane DHA processes are compared with the typical multi-step natural gas to aromatics production process via methanol synthesis. Sensitivity of internal rate of return (IRR) and net present value (NPV) to various economic and process parameters such as plant scale, desired rate of return, reactor cost, feedstock and utility cost, catalyst variable cost, and MW reactor cost is studied. Here, electric equivalent efficiency of the conventional methane DHA process is found to be 69.2 % and 67.3 % at 750 °C and 800 °C, respectively, while the MW-assisted methane DHA process has the electric equivalent efficiency of 48.9 % at 800 °C. IRRs of the conventional methane DHA process at 750 °C and 800 °C, and MW-assisted process are 15.2 %, 17.5 %, and 18.8 %, respectively for a methane feed flowrate of 19,782 kg/h, while the IRR of the multi-step natural gas to aromatics production process is estimated to be 0 % for the same plant scale. Impact of change in the methane price, electricity price, and catalyst cost is found to be considerable on the process economics, while the cost of the MW reactor is found to have negligible impact.

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

Towards modeling phage therapy

Patients infected with life-threatening multi-drug resistant (MDR) bacteria have been treated with cocktails of bacteriophages. This is a complicated form of personalized medicine as the phages given to a patient have to be selected beforehand on the basis of their lytic capacity of the infecting bacteria. Because bacteria rapidly become resistant, the evolution of resistance to a diverse cocktail of phages is a complicated dynamical process, during which competing bacterial strains replace one another by accumulating several resistance mechanisms, each of which may involve a fitness cost. As a consequence, it is typically not known why a particular phage therapy succeeded or failed, and how one can optimize the composition of the cocktails to maximize the rate of success. To improve upon this, we extend an existing in vivo -calibrated mouse model into a novel mathematical model for the human situation, and include multiple phages infecting multiple bacterial strains, differing in their resistance to each of the phages. We adjust several parameter estimates of the bacterial model to the human situation, and use the model to describe a successful case of phage therapy involving several cocktails, each containing several phages. In the model, treatment success crucially depended on pretreatment resistance levels, and on the diversity and the timing of the cocktails. Once an appropriate cocktail is found, it is less important to further optimize the infection rates of the phages. Resistant bacterial strains expand rapidly when sensitive strains decline, and the higher the infectivity of the phages, the faster resistant strains expand. Because resistance evolves rapidly, it is best to provide a diverse set of phages right from the start of therapy, i.e., to hit hard and early, and create a high genetic barrier to bacterial resistance.

59 BASIC BIOLOGICAL SCIENCES↗