Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimal start”

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 631 records · Page 35

Arcjet thruster research and technology. Phase 1: Executive summary

The principal objective of this two phase program is to conduct the development research required to make the low power arcjet a flight ready technology. Many important results were obtained during Phase 1 to move closer to this objective. Fundamental analyses were performed of the arcjet nozzle, the gas kinetic reaction effects, the thermal environment, and the arc stabilizing vortex. These aided the conceptual understanding of the arcjet and guided design work. A hydrazine (N2H4) arcjet was designed that combined a flight qualified catalyst bed with a modular arcjet. Extensive testing was performed which demonstrated the feasibility of using this propellant in an arcjet for the first time. Startup techniques were developed, stability maintained, material compatibility tests conducted, and performance mapping tests performed. Specific impulse values from 400 to 730 seconds were produced with a non-optimized design. These levels are higher than were originally thought possible and proved that extremely high enthalpy values can be obtained with constricted arc technology. Erosion rate data are promising for lifetime extensions to meet flight application requirements. Power control unit (PCU) development was started with the design and fabrication of a laboratory high switching frequency supply. Valuable data were obtained on PCU operation and on the interaction with the dynamic arc. Phase 2 efforts presently underway are resolving key issues for multi-hundred hour lifetimes, are continuing to investigate arcjet/PCU interactions, and will demonstrate duty cycle N2H4 arcjet/PCU operation in a simulated flight mode for lifetimes consistent with initial applications.

Source record↗

Numerical Investigations of High Pressure Acoustic Waves in Resonators

This presentation presents work on numerical investigations of nonlinear acoustic phenomena in resonators that can generate high-pressure waves using acoustic forcing of the flow. Time-accurate simulations of the flow in a closed cone resonator were performed at different oscillation frequencies and amplitudes, and the numerical results for the resonance frequency and fluid pressure increase match the GRC experimental data well. Work on cone resonator assembly simulations has started and will involve calculations of the flow through the resonator assembly with and without acoustic excitation. A new technique for direct calculation of resonance frequency of complex shaped resonators is also being investigated. Script-driven command procedures will also be developed for optimization of the resonator shape for maximum pressure increase.

Athavale, Mahesh↗

DSN Beowulf Cluster-Based VLBI Correlator

The NASA Deep Space Network (DSN) requires a broadband VLBI (very long baseline interferometry) correlator to process data routinely taken as part of the VLBI source Catalogue Maintenance and Enhancement task (CAT M&E) and the Time and Earth Motion Precision Observations task (TEMPO). The data provided by these measurements are a crucial ingredient in the formation of precision deep-space navigation models. In addition, a VLBI correlator is needed to provide support for other VLBI related activities for both internal and external customers. The JPL VLBI Correlator (JVC) was designed, developed, and delivered to the DSN as a successor to the legacy Block II Correlator. The JVC is a full-capability VLBI correlator that uses software processes running on multiple computers to cross-correlate two-antenna broadband noise data. Components of this new system (see Figure 1) consist of Linux PCs integrated into a Beowulf Cluster, an existing Mark5 data storage system, a RAID array, an existing software correlator package (SoftC) originally developed for Delta DOR Navigation processing, and various custom- developed software processes and scripts. Parallel processing on the JVC is achieved by assigning slave nodes of the Beowulf cluster to process separate scans in parallel until all scans have been processed. Due to the single stream sequential playback of the Mark5 data, some ramp-up time is required before all nodes can have access to required scan data. Core functions of each processing step are accomplished using optimized C programs. The coordination and execution of these programs across the cluster is accomplished using Pearl scripts, PostgreSQL commands, and a handful of miscellaneous system utilities. Mark5 data modules are loaded on Mark5 Data systems playback units, one per station. Data processing is started when the operator scans the Mark5 systems and runs a script that reads various configuration files and then creates an experiment-dependent status database used to delegate parallel tasks between nodes and storage areas (see Figure 2). This script forks into three processes: extract, translate, and correlate. Each of these processes iterates on available scan data and updates the status database as the work for each scan is completed. The extract process coordinates and monitors the transfer of data from each of the Mark5s to the Beowulf RAID storage systems. The translate process monitors and executes the data conversion processes on available scan files, and writes the translated files to the slave nodes. The correlate process monitors the execution of SoftC correlation processes on the slave nodes for scans that have completed translation. A comparison of the JVC and the legacy Block II correlator outputs reveals they are well within a formal error, and that the data are comparable with respect to their use in flight navigation. The processing speed of the JVC is improved over the Block II correlator by a factor of 4, largely due to the elimination of the reel-to-reel tape drives used in the Block II correlator.

Rogstad, Stephen P.↗

The Utilization Profiles of the CCSDS Unified Space Link Protocol (USLP)

The purpose of this paper is to identify the utilization profiles for interfacing the Data Protocol Sublayer using the Unified Space Link Protocols (USLP) (reference 1) with the space link coding procedures as specified in the CCSDS Coding & Synchronization Blue Books (references 2 through 5), used in both telecommand and telemetry applications. This paper describes how the USLP Protocol utilizes the coding and synchronization sublayer to support: a. Direct to Earth (DTE) telemetry links for engineering and science data b. Direct to Earth (DTE) telemetry links for very high rate science data c. Direct from Earth (DFE) command, sequencing and flight software loads d. Space to Space Links (Proximity) utilized by orbiters for data exchange to/from surface bound assets. The CCSDS has divided the functions of the Data Link Layer into two sublayers: the Data Link Protocol Sublayer (DLP-SL) and the Coding and Synchronization Sublayer (CS-SL). The Data Link Protocol Sublayer (DLP-SL) interfaces to the users, accepting the data that is to be transported, on the sending side of the link, and delivering that data on the receiving end. The Transfer Frame is the data unit that is transferred across the Data Link Protocol Sublayer and the Coding and Synchronization Sublayer boundary. The Coding and Synchronization Sublayer (CS-SL) provides the encoding, randomization, and frame synchronization functions that prepares the USLP Transfer Frame for transport across the space link. The CS-SL is divided into 2 processes: 1) The Frame Interface Processes (FIP) performs the interface functions required to prepare the data for delivery to the Coding/Decoding Process (CDP). This process includes prepending a Frame Start Marker to the provided frame, when management has designated that the frame is not to be aligned to the codeblock or when there is no block code used. 2) The Coding/Decoding Process (CDP) performs the forward error correction processes that are used to optimize the performance of the link and minimize the error rate. The CDP creates the symbol stream that is delivered to the Physical Layer. The transfer of the USLP transfer frames across different types of space links is the focus of this paper. The Protocol Data Unit (PDU) that is passed in both directions between the Data Link Protocol Sublayer (DLP-SL) and Coding and Synchronization Sublayer (CS-SL) is the transfer frame. The USLP frame structure provides flexibility that can be constrained by the functions utilized within the CS-SL that prepare the transfer frame for transit. For example, the USLP transfer frame contains a length field that enables the frame to be of variable length but CS-SL under certain conditions may constrain the frame to be fixed in length. This paper describes 5 operational modes available for use by the Data Link Layer to provide data exchange across the USLP space link. These modes are different because different operational requirements apply to vastly different types of space links and thus the communications implementation requirements differ. The environmental issues include the power or energy available, the distance between the end points of the link, the complexity of the equipment available at those end points, the atmospheric conditions and radiometric frequency selection. The CS-SL utilizes different forward error correcting codes supported by specific operational modes to configure the data for transit. This paper describes all of the operational modes in a series of data models which decompose the functionality between the Data Link Protocol Sublayer and the Coding and Synchronization sublayer. The operational modes described are: 1. Uncoded Mode: has been used for short links that contain significant available power to provide an acceptable frame error rate. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 2. Convolutional Only Mode: is currently the prime forward error correction coding used for the proximity links. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 3. Variable Length Frame Aligned to Variable Length Codeblock (TC): is used for Direct from Earth links were power levels are high and the simple, least complex code i.e., the BCH code is used. This mode has been in use since the early 1970s. The BCH code is a short code and the decoder is easy to implement. 4. Fixed Length Frame Aligned to Fixed Length Codeblock (AOS/TM): was introduced when the concatenated Convolutional and Reed-Solomon Code was formulated to provide significant reduction in link data error rate and the ability to determine if there was an error in the decoded codeblock. The frame is aligned to the codeblock so that there is a one to one relationship of frame errors to codeblock errors without additional error detection coding being added. This mode requires the protocol frames to be the exact size of the message portion of the codeblock. 5. Frames Unaligned to Fixed Length Codeblocks (Currently used for very high rates and space to space links): This mode is currently used for missions that have a very high data rate that can be controlled adaptively as the environment changes and as the next generation operating mode for the proximity link. This mode from a coded data stream point of view is exactly like that described in 4. above, except that the frame need not be aligned to the codeblock. There is no requirement on frame length when using this mode. Thus when using USLP it can be used to support links that require short or long frames. There is also no mandatory requirement that frames cannot be separated by idle data reducing the tight data rate connection requirements between the data link protocol sublayer and the coding & synchronization sublayer. In conclusion, how these operational modes can be put to use in mission operational scenarios is described for Direct from Earth links (DFE), Direct to Earth links (DTE), and Proximity links.

Greenberg, E.↗

Jump-Start U.S. Manufacturing of High-Performance Permanent Magnets

Gas atomized feedstock powders for Nd-Fe-B based magnets are more robust and cost effective than the strip-cast sheet currently used as feedstock for milling by the magnet industry. However, gas atomized powder’s crystallinity and grain boundary phases are less ideal than the strip cast flake. We proposed to optimize the atomization process as well as the powder chemistry to address this issue and to eliminate heavy rare earth elements from these magnet alloys. Our goal is to help GKN-Hoeganaes, a US metal powder company, who used to make Nd-Fe-B magnets, to bring back US manufacturing of these high energy density permanent magnets.

36 MATERIALS SCIENCE↗

Optimization of the design parameters for a wide-band radiometric system

The optimun design parameters for a swept frequency wide-band radiometric antenna system for spacecraft applications are studied. Wide band antenna systems are needed to observe layered surfaces which are frequency sensitive and require multiple measurements for interpretation. The lowest frequency band of interest is between 1.4 to 2.8 Ghz. Starting with a given size reflector fed in the offset mode by a corrugated horn located at the focus of the parabola, the primary performance indexes; e.g., half power beamwidth, cross polarization level, and overall beam efficiency were calculated over a wide frequency range (two to one) for different physical horn dimensions and for different values of f/D ratio. These data are used to find the best design under given restriction of reflector size and blockage.

Agrawal, P. K.↗

Towards On-Chip Learning for Low Latency Reasoning with End-to-End Synthesis

The Software Defined Architectures (SODA) Synthesizer is an open-source compiler-based tool able to automatically generate domain-specialized systems targeting Application-Specific Integrated Circuits (ASICs) or Field Programmable Gate Arrays (FPGAs) starting from high-level programming. SODA is composed of a frontend, SODA-OPT, which leverages the multilevel intermediate representation (MLIR) framework to interface with productive programming tools (e.g., machine learning frame-works), identify kernels suitable for acceleration, and perform high-level optimizations, and of a state-of-the-art high-level synthesis backend, Bambu from the PandA framework, to generate custom accelerators. One specific application of the SODA Synthesizer is the generation of accelerators to enable ultra-low latency inference and control on autonomous systems for scientific discovery (e.g., electron microscopes, sensors in particle accelerators, etc.). This paper provides an overview of the flow in the context of the generation of accelerators for edge processing to be integrated in transmission electron microscopy (TEM) devices, focusing on use cases from precision material synthesis. We show the tool in action with an example of design space exploration for inference on reconfigurable devices with a conventional deep neural network model (LeNet). Finally, we discuss the research directions and opportunities enabled by SODA in the area of autonomous control for scientific experimental workflows.

Castellana, Vito G.↗

A Self-Sustained CPS Design for Reliable Wildfire Monitoring

Continuous monitoring of areas nearby the electric grid is critical for preventing and early detection of devastating wildfires. Existing wildfire monitoring systems are intermittent and oblivious to local ambient risk factors, resulting in poor wildfire awareness. Ambient sensor suites deployed near the gridlines can increase the monitoring granularity and detection accuracy. However, these sensors must address two challenging and competing objectives at the same time. First, they must remain powered for years without manual maintenance due to their remote locations. Second, they must provide and transmit reliable information if and when a wildfire starts. The first objective requires aggressive energy savings and ambient energy harvesting, while the second requires continuous operation of a range of sensors. To the best of our knowledge, this paper presents the first self-sustained cyber-physical system that dynamically co-optimizes the wildfire detection accuracy and active time of sensors. The proposed approach employs reinforcement learning to train a policy that controls the sensor operations as a function of the environment (i.e., current sensor readings), harvested energy, and battery level. Here, the proposed cyber-physical system is evaluated extensively using real-life temperature, wind, and solar energy harvesting datasets and an open-source wildfire simulator. In long-term (5 years) evaluations, the proposed framework achieves 89% uptime, which is 46% higher than a carefully tuned heuristic approach. At the same time, it averages a 2-minute initial response time, which is at least 2.5× faster than the same heuristic approach. Furthermore, the policy network consumes 0.6 mJ per day on the TI CC2652R microcontroller using TensorFlow Lite for Micro, which is negligible compared to the daily sensor suite energy consumption.

54 ENVIRONMENTAL SCIENCES↗

Genetic Algorithm for Optimization of Neural Networks for Bayesian Inference of Model Uncertainty

The objective of this work was to develop a genetic optimization algorithm that can design a neural network capable of producing uncertainty estimates along with predictions. This algorithm is necessary because the inclusion of uncertainty modeling in a neural network greatly complicates the network’s design space, making the development of a converging model extremely difficult and time consuming. The genetic algorithm presented in this work uses a number of value ranges for various configurable neural network parameters to create a randomly generated population of network architectures. The initially generated population is then evolved over the course of several generations, with the best performing models breeding to produce novel network configurations. Mutations are randomly applied to the network designs to facilitate the development of adaptations beneficial to the task being performed. An experiment was conducted to validate the proposed algorithm, in which the genetic optimizer was tasked with producing a neural network capable of predicting the sound pressure level (SPL) resulting from jet-surface interaction (JSI) noise. The data used for this task was generated at the NASA Glenn Research Center in the Aero-Acoustic Propulsion Laboratory. Starting with an initial population size of 35 randomly generated networks, and evolved over the course of 10 generations, the genetic algorithm produced a design able to predict SPL as a result of JSI noise within 0.272 dB, on average.

Genetic algorithm↗

Future of the Search for Life: Workshop Report

The 2-week, virtual Future of the Search for Life science and engineering workshop brought together more than 100 scientists, engineers, and technologists in March and April 2022 to provide their expert opinion on the interconnections between life-detection science and technology. Participants identified the advances in measurement and sampling technologies they believed to be necessary to perform in situ searches for life elsewhere in our Solar System, 20 years or more in the future. Among suggested measurements for these searches, those pertaining to three potential indicators of life termed “dynamic disequilibrium,” “catalysis,” and “informational polymers” were identified as particularly promising avenues for further exploration. For these three indicators, small breakout groups of participants identified measurement needs and knowledge gaps, along with corresponding constraints on sample handling (acquisition and processing) approaches for a variety of environments on Enceladus, Europa, Mars, and Titan. Despite the diversity of these environments, sample processing approaches all tend to be more complex than those that have been implemented on missions or envisioned for mission concepts to date. The approaches considered by workshop breakout groups progress from nondestructive to destructive measurement techniques, and most involve the need for fluid (especially liquid) sample processing. Sample processing needs were identified as technology gaps. These gaps include technology and associated sampling strategies that allow the preservation of the thermal, mechanical, and chemical integrity of the samples upon acquisition; and to optimize the sample information obtained by operating suites of instruments on common samples. Crucially, the interplay between science-driven life-detection strategies and their technological implementation highlights the need for an unprecedented level of payload integration and extensive collaboration between scientists and engineers, starting from concept formulation through mission deployment of life-detection instruments and sample processing systems.

Marc Neveu↗

Strength optimization of alpha-SiC by improved processing

Silicon carbide is of great interest for structural use in aircraft and automobile engines. This ceramic combines high thermal conductivity and low coefficient of thermal expansion, and consequently has good thermal shock resistance. However, like other ceramics, silicon carbide shows strength variability due to processing flaws such as large voids, shrinkage cracks, and inclusions. Agglomerates in the starting powder seem to be the predominant cause of such defects. Improved processing techniques such as slurry pressing and hot isostatic pressing (HIP) were employed to minimize these defects and to improve strength and reliability in the fabricated material. For this purpose 2-inch diameter disks were fabricated by various consolidation techniques. These include: dry pressing and sintering, slurry pressing and sintering, and slurry pressing and HIPing. The results are discussed.

Dutta, Sunil↗

Project Title: Demonstration High Temperature Superconducting NonPlanar Stellarator Magnet with Advanced Manufactured Assemblies

This is the final report for the project “Demonstration High Temperature Superconducting Non- Planar Stellarator Magnet with Advanced Manufactured Assemblies”, funded by DOE, and performed by Type One Energy from September, 2020, to March 2024 involving the Fusion Technology Institute at the University of Wisconsin–Madison, the Plasma Science and Fusion Center (PSFC) at the Massachusetts Institute of Technology (MIT) and Commonwealth Fusion Systems (CFS) to design and fabricate the first non-planar HTS (REBCO) coil for a high-field stellarator based on the SPARC tokamak’s VIPER cable concept. Stellarators at high fields make high-temperature superconducting magnets necessary for a compact fusion device. But the asymmetric and non-planar nature of its components, especially the magnets make it difficult for scalable producibility. To address these challenges, two promising technologies have emerged: advanced manufacturing (AM) for the supporting plates for forming the magnets, and high-temperature superconducting (HTS) cables inside the plates. AM has advanced enough to produce stellarator components with the necessary geometric complexity, size, and the precision, leading to potentially significant reductions in production time, cost, and waste. The cost of HTS tape has decreased dramatically, and progress in HTS planar magnet development has reached a point where it can be proposed for application to complex 3D non-planar magnets. The main objective of this project is to develop, demonstrate and pre-commercialize a novel, non-planar HTS coil shape that remains superconducting to achieve production scalable reductions in time and cost and performance. The proposed technology is based on the novel concept of a precision sub-scale HTS nonplanar coil assembly. This project focuses on the design, fabrication, material optimization of cable design, and validation and demonstration of the high current carrying capability of superconducting magnets and their support in a complex 3D shape needed for application to stellarator magnetic plasma confinement. The specific objectives of this research program include: (1) The successful application of metal AM to build a precision sub-scale HTS nonplanar coil, (2) An HTS cable and cross-section design that can conform to the required nonplanar coil shape (bend radii as tight as 10-cm) and remains superconducting at an engineering current density of 1.35 kA/cm 2 at 77 K and 1 tesla at the conductor (5 kA in the cable). To achieve the above challenging goals, we have formed a multidisciplinary research team consisting of members from Type One Energy and UW-Madison, MIT PSFC and CFS with complementary skills and strong facilities. The team worked collaboratively on fundamental and applied research on the following three major technical areas: (1) Design, fabrication, and optimization of non-planar HTS Cable The ultimate goal of the project is to determine if commercial REBCO tapes and additive manufacturing can be used to fabricate high field (≥ 10T) non-planar coils with tight bending radii (≃ 100mm) and with a degradation of the critical current (Ic) smaller than 20% with respect to the expected performance. We started with shorter length cable to evaluate the scalability of the production process and eventually reached multiple turns for higher magnetic fields. Our findings suggest that a stellarator coil system of a relevant size, characterized by its asymmetric and non-planar components, can be fabricated using a formed cable in plate method. This system can be simulated using a large-scale modeling approach. The use of hybrid modeling 3 techniques will be pivotal in reducing the complexity of the model and in assessing expected performance in designs. (2) Modeling and simulation of the non-planar HTS Cable Multiphysics simulations are performed using the commercial software and are carried out in self-field conditions, involving 2D and 3D models and twisted around one slot of twist-pitched VIPER cable. Multiphysics simulations are mainly focused on the coil for the critical current evaluation, the magnetic field map, self-Lorentz forces and mechanical, and magnetothermal behavior and the quench dynamics. The detailed model and prediction of the superconducting performance of a stellarator-relevant demonstration cable from numerical simulations supports the results from the actual testing backing the results. A detailed description and results are provided in the later sections. (3) Design, fabrication, and optimization of support for the non-planar HTS Cable The team developed an additive manufactured (AM) coil positioning plate that formed into the required non-planar geometry (with bend radii as tight as 10-cm) and to acceptable tolerances required for a stellarator magnet: (+0.25-mm from ideal on dimensions of coil positioning plates and up to +1-mm from ideal for position of wound coil). The plate materials is also included in this selection process from fabrication and 3D printing perspective and commensurate with eventual application to a fusion reactor. From the cost effectiveness point of view, the HTS coil and plate has the potential to cost less than that made in conventional methods with less waste (<75% waste) reducing time (<50%) and cost (<50%), especially as the AM field matures. The application of advanced manufacturing in the construction of the support plates will also lead to cost reduction, as the cables can be easily replaced, thereby making the assembly modular. With such high primary cost and time savings, high current densities and magnetic field, the funded R&D work has validated the designs, proven the feasibility, and characterized the performance of the HTS coil and plate assembly, paving the way for a relevant-size stellarator coil system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tutorial review on structure – dendrite growth relations in metal battery anode supports

This tutorial review explains the emerging understanding of the surface and bulk chemistry – electrochemical performance relations in anode supports (aka secondary current collectors, substrates, templates, hosts) for lithium, sodium and potassium metal batteries (LMBs, SMBs or NMBs, and KMBs or PMBs). In relation to each section, the possible future research directions that may yield both new insight and improved cycling behavior are explored. Representative case studies from Li, Na and K metal anode literature are discussed. The tutorial starts with an overview of the solid electrolyte interphase (SEI), covering both the “classic” understanding of the SEI structure and the “modern” insights obtained by site-specific cryogenic stage TEM analysis. Next, the multiple roles of supports in promoting cycling stability are detailed. Without an optimized support architecture, the metal–electrolyte interface becomes geometrically unstable at a lower current density and cycle number. Taking into consideration the available literature on LMBs, SMBs and KMBs, it is concluded that effective architectures are geometrically complex and electrochemically lithiophilic, sodiophilic or potassiophilic, so as to promote conformal electrochemical wetting of the metal during plating/stripping. One way that philicity is achieved is through support oxygen surface chemistry, which yields a reversibly reactive metal–support interface. Examples of this include the well-known oxygen–carbon moieties in reduced graphene oxide (rGO), as well as classic ion battery reversible conversion reaction oxides such as SnO 2 . Unreactive surfaces lead to dewetted island growth of the metal, which is a precursor to dendrites, and possibly to non-uniform dissolution. Surveying the literature on various Li, Na and K metal supports, it is concluded that the key bulk thermodynamic property that will predict electrochemical wetting behavior is the enthalpy of infinite solution (Δ sol H ∞ ) of the metal (solute) into the support (solvent). Large and negative Δ sol H ∞ promotes uniform metal wetting on the support surface, corresponding to relatively low plating overpotential. Positive Δ sol H ∞ promotes dewetted islands and a relatively high overpotential. This simple rule explains a broad range of studies on Li, Na and K metal – support interactions, including the previously reported correlation between mutual solubility and wetting.

25 ENERGY STORAGE↗

Flex-PCB Quench Antenna Developments at FNAL

We report flexible printed-circuit-board-based quench antenna arrays have been developed at FNAL since 2017. They proved to be low-noise, low-cost, sensitive, and non-invasive sensors for quench characterization in superconducting magnets, having the potential to extend even to quench detection. Starting with applications of small area sensors of this type, we have explored feasibility for operations at high pressure, and have invested in multi-channel quench antenna arrays covering the complete area of the innermost magnet layer of the conductor. We've also tested and optimized sensor geometry and relative orientations to increase sensitivity and spatial resolution. Beyond developments of warm and cold bore antennas, we've progressed to concepts which attach quench antennas to the conductor surfaces of accelerator magnets. Multiple programs and projects at FNAL and collaborating institutions benefited from this work already and our plans target further extensions of capabilities and applicability of these devices. In this paper we present the logic behind the development steps taken so far, status of our work and direction of near-term research on quench antennas at Fermilab.

43 PARTICLE ACCELERATORS↗

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↗

A method for designing robust multivariable feedback systems

A new methodology is developed for the synthesis of linear, time-invariant (LTI) controllers for multivariable LTI systems. The aim is to achieve stability and performance robustness of the feedback system in the presence of multiple unstructured uncertainty blocks; i.e., to satisfy a frequency-domain inequality in terms of the structured singular value. The design technique is referred to as the Causality Recovery Methodology (CRM). Starting with an initial (nominally) stabilizing compensator, the CRM produces a closed-loop system whose performance-robustness is at least as good as, and hopefully superior to, that of the original design. The robustness improvement is obtained by solving an infinite-dimensional, convex optimization program. A finite-dimensional implementation of the CRM was developed, and it was applied to a multivariate design example.

Milich, David Albert↗

Minimum time attitude slewing maneuvers of a rigid spacecraft

The minimum time attitude slewing motion of a rigid spacecraft with its controls provided by torques and forces, which have their upper and lower limits prescribed, is considered. The two-point boundary-value problem is derived by applying the Pontriagin's Maximum Principle to the system and solved by using a quasi-linearization algorithm. The nominal solutions to the problem as well as the starting values of the total slewing time and the unknown initial costates for this algorithm are generated by using Euler's eigenaxis rotation theorem. It is pointed out that one of the four initial costates associated with the quaternions can be arbitrarily selected without affecting the optimal controls and, thus, simplifying the computation. The minimum slewing time is determined by shortening the total slewing time until at least one of the controls becomes a bang-bang type. Several numerical tests for the rigidized SCOLE model are presented to show the applications of the methods.

Li, Feiyue↗

A method for designing robust multivariable feedback systems

A new methodology is developed for the synthesis of linear, time-invariant (LTI) controllers for multivariable LTI systems. The aim is to achieve stability and performance robustness of the feedback system in the presence of multiple unstructured uncertainty blocks; i.e., to satisfy a frequency-domain inequality in terms of the structured singular value. The design technique is referred to as the causality recovery methodology (CRM). Starting with an initial (nominally) stabilizing compensator, the CRM produces a closed-loop system whose performance-robustness is at least as good as, and hopefully superior to, that of the original design. The robustness improvement is obtained by solving an infinite-dimensional, convex optimization program. A finite-dimensional implementation of the CRM was developed, and it was applied to a multivariate design example.

Milich, David A.↗