Search NASASearch

SEARCH · Search NASA

Results for “Non-flammable Polymer Separator”

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.

641 records · Page 26

Battery Failure Databank

The Battery Failure Databank contains thermal runaway results gathered from nearly 300 small format fractional thermal runaway calorimetry (S-FTRC) experiments. A majority of these experiments were conducted at synchrotron facilities where high-speed x-ray videography was conducted of the cell while tested inside of the S-FTRC. The databank is a two-component system which consists of a Microsoft ExcelTM spreadsheet which provides S-FTRC results in tabular format and a radiographic video library containing the high-speed x-ray videos. Overall, the databank provides thermal results from S-FTRC experiments conducted on a mixture of commercially available lithium-ion (Li-ion) cells and specialized Li-ion test cells with varying cell format (18650, 21700, and D-cell) and trigger mechanism (heaters, heaters plus internal short circuiting device, and nail penetration). The radiography video component provides insight into the initiation and propagation of TR in the cells, in addition to consequences of TR measured by the FTRC. Fractions of mass ejected for the cell types are separated into regimes based on the different failure modes of the cells, such as purely venting, partial ejection, and total ejection. With knowledge of the rate and characterization of the internal degradation of the cells during TR, and the amounts of mass ejected and unrecovered, the extent of TR, and therefore the mitigation of TR, is revealed relative to cell types and failure modes.

Lithium-ion battery

Strong gradient neoclassical transport in the plateau regime

Strong gradient regions in tokamaks such as the pedestal or internal transport barriers are regions of reduced turbulence where neoclassical transport can play a dominant role. In pedestals, gradient lengths comparable to the ion poloidal gyroradius have been measured. Standard neoclassical theory can miss important strong gradient effects in these regions because it assumes that the gradient length scales of density, temperature and potential are larger than the ion poloidal gyroradius. We extend plateau regime neoclassical theory into regions of gradients of the order of the ion poloidal gyroradius to capture strong gradient effects on transport processes in the pedestal and internal transport barriers. The fundamental idea behind our new framework is to keep a scale separation between the orbit widths and the gradient length scales by performing a large aspect ratio expansion. In the plateau regime, strong gradients cause poloidal variation that is in–out as well as up–down asymmetric. We study two different test cases assuming either radial force balance or the absence of turbulence and show that strong gradient effects can enhance or reduce standard neoclassical theory predictions in the plateau regime in strong gradient regions.

fusion plasma

Secondary-Side Active Rectifier Synchronization and Control in LCC-LC Compensated Inductive Power Transfer Systems

This paper proposes a communication-less synchronization and output power control strategy for LCC-LC compensated inductive power transfer systems with a secondaryside active rectifier. By replacing the passive rectifier with an active bridge, the proposed architecture eliminates the need for a separate secondary-side DC-DC converter and enables direct regulation of battery power. The key challenge in such a system is achieving robust synchronization of the secondary active bridge with the primary inverter without access to primary-side signals. To address this issue, the secondary-side variables are analyzed and the current through the secondary series inductor is identified as a load-independent synchronization variable with a fixed phase relationship to the primary excitation. A fixed-frequency second order generalized (SOGI)-based phase locked loop (PLL) is used to estimate the switching phase and frequency, and a phase-shift-based output power controller is developed using a phasor-transformer model of the active bridges. A modified control formulation is further introduced to remove phase-dependent loop-gain variation and simplify controller design. Simulation results validate the proposed method under both fixed and varying switching-frequency conditions, demonstrating successful frequency tracking from 82 kHz to 88 kHz and accurate output power regulation over a wide operating range. The proposed method offers a compact and effective solution for high-power IPT systems by enabling secondary-side synchronization and control without wireless communication or an additional DC-DC conversion stage.

Fernandes, Arnold Anthony [ORNL] (ORCID:0009000631

Technoeconomic and Emissions Analysis of the Hybrid Redox Process for the Production of Acetic Acid with CO 2 Utilization

The production of oxygenated hydrocarbons, such as acetic acid, using captured CO 2 is a promising pathway to reduce greenhouse gas emissions in the chemical industry. The use of a chemical looping‐based hybrid redox process (HRP) is proposed to convert CO 2 and natural gas into separate CO and syngas streams that can be used to produce various commodity oxygenates, while allowing the beneficial utilization of captured CO 2 . Here, a detailed technoeconomic analysis of HRP applied to the production of acetic acid is presented. Emissions and energy analyses show the ability of HRP to lower the CO 2 emissions for acetic acid synthesis by 74% compared to a conventional steam and autothermal reforming route. HRP also offers a potential 34% reduction in capital costs. Compared to a dry reforming based acetic acid production route, HRP has the potential for significantly lower costs. If integrated with a low carbon energy source, HRP has the potential to achieve a negative emission of greenhouse gas (‐0.50 kg CO 2 per kg acetic acid).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Stationary Testing of a Pulse Tube Cryocooler Designed to Rotate at 6800 RPM with the Superconducting Rotor of an Electric Machine

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Most of these electric machines employ superconductors on their rotor, thereby requiring cryogenic cooling of the rotating system. An attractive solution is to conductively cool the rotor using a cryocooler that rotates with the rotor. NASA has been developing a 50 K pulse tube cryocooler for its High Efficiency Megawatt Motor (HEMM) that can operate while rotating at 6,800 rpm. The design of this cryocooler and testing of its linear motor have been discussed in prior work. This presentation presents the first measurements of the complete HEMM cryocooler. The design and results of a stationary (non-rotating) test of the cryocooler are described.

electrified aircraft propulsion

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely

The occurrence rate of galaxies with polar structures may be significantly underestimated

Polar-ring galaxies are photometrically and kinematically decoupled systems that are highly inclined to the major axis of the host. These galaxies have been explored since the 1970s, but the rarity of these systems has made such systematic studies difficult. However, over 250 good candidates have been identified. In this work, we examine a sample of over 18 000 galaxies from the Sloan Digital Sky Survey (SDSS) Stripe 82 for the presence of galaxies with polar structures. Using deep SDSS Stripe 82, DESI Legacy Imaging Surveys, and Hyper Suprime-Cam Subaru Strategic Program, we selected 53 good candidate galaxies with photometrically decoupled polar rings, 9 galaxies with polar halos, 6 galaxies with polar bulges, and 34 possibly forming polar-ring galaxies, versus 13 polar-ring candidates previously selected in Stripe 82. Our results suggest that the occurrence rate of galaxies with polar structures may be significantly underestimated, as revealed by the deep observations, and may amount to 1–3% of non-dwarf galaxies.

Astronomy & Astrophysics

Broadband and Tunable Microwave Absorption Properties from Large Magnetic Loss in Ni–Zn Ferrite

Highly effective electromagnetic (EM) wave absorber materials with strong reflection loss (RL) and a wide absorption bandwidth (EBW) in gigahertz (GHz) frequencies are crucial for advanced wireless applications and portable electronics. Traditional microwave absorbers lack magnetic loss and struggle with impedance matching, while ferrites are stable, exhibit excellent magnetic and dielectric losses, and offer better impedance matching. However, achieving the desired EBW in ferrites remains a challenge, necessitating further composition design. In this study, impedance matching is successfully enhanced and EBW in Ni–Zn ferrite is broadened by successive doping with Mn and Co , without incorporation of any polymer filler. It is found that Ni 0.4 Co 0.1 Zn 0.5 Fe 1.9 Mn 0.1 O 4 material exhibits exceptional EM wave absorption, with a maximum RL of −48.7 dB. It also featured a significant EBW of 10.8 GHz, maintaining a 90% absorption rate (RL < −10 dB) for a thickness of 4.5 mm. These outstanding properties result from substantial magnetic losses and favorable impedance matching. These findings represent a significant step forward in the development of microwave absorber materials, addressing EM wave pollution concerns within GHz frequencies, including the frequency band used in popular 5G technology.

36 MATERIALS SCIENCE

Interactions Between the Space Station and the Environment: A Preliminary Assessment of EMI

A review of the interactions between proposed Space Station systems/payloads and the environment that contribute to electromagnetic interference was performed. Seven prime sources of interference have been identified. These are: The Space Station power system; active experiments such as beam injection; ASTROMAG; ram and wake density gradients; pick up ions produced by vented or offgassed clouds; waves produced by current loops that include the plasma and structure; arcing from high voltage solar arrays (or possible ESD in polar orbit). This review indicates that: minimizing leakage current from the 20 kHz power system to the structure; keeping the surfaces of the Space Station structure, arrays, and radiators nonconducting; minimizing venting of payloads or systems to non-operational periods; careful placement of payloads sensitive to magnetic field perturbations or wake noise; and designing an operational timeline compatible with experiment requirement are the most effective means of minimizing the effects of this interference. High degrees of uncertainty exist in the estimates of magnitudes of gas emission induced EMI, radiation of 20 kHz and harmonics, ASTROMAG induced interference, and arc threshold/frequency of the solar array. These processes demand further attention so that mitigation efforts are properly calibrated.

G B Murphy

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE

Osmotic and phoretic competition explains chemotaxic assembly and sorting

Microscale objects responding to chemical gradients by migrating toward or away from a preferred species is a simple yet constitutive mechanism by which transport occurs in biological organisms. Synthetic chemotaxis provides key physical descriptions of simplified systems that can be used in biological models, or in the creation of advanced responsive material systems. In this article, we provide a quantitative framework for understanding synthetic chemotaxis of microparticles which involves a competition between phoresis and osmosis. We present separate quantitative measurements of phoresis and osmosis acting on individual taxing particles, finding that phoresis follows the long-predicted v ∼ 1 / r 2 scaling while the osmotic contribution depends on the geometry and details of the system, and must be solved on a case-by-case basis. Through this, we are able to develop a more accurate picture of particle transport at the single particle level. Equipped with this approach, we go on to describe how high concentrations of particles in a symmetric chemical gradient grow close-packed hives that reach a steady-state size tunable through light intensity or particle size. Last, we demonstrate that mixed particles experiencing the same chemical gradient will selectively migrate toward or away depending on the nature of the particle surface, thereby locally sorting out a particular species. We anticipate these results will be important in describing both biological and synthetic chemotaxis in phoretic systems and should bring a wealth of studies that take advantage of competing osmotic flows to illicit unexpected dynamic active behavior.

Science & Technology - Other Topics

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Gorazd Medic

A Probabilistic Approach to Load Modeling for Central HVAC Systems in Large Commercial Buildings for Retrofit Decisions Under Uncertainty

Retrofitting central HVAC systems in large commercial buildings with advanced technologies like heat recovery chillers (HRCs) offers a significant opportunity to enhance energy efficiency. However, analyzing these retrofits is challenging with traditional whole-building simulation tools, which require intensive calibration and struggle to model innovative system configurations and controls. To overcome these limitations, this study proposes a load profilebased retrofit analysis framework that provides better decisions under uncertainty. The main focus of this paper is the development of a probabilistic load profile model that can be used in the framework by using exploratory data analysis (EDA) of measured building data to properly quantify its inherent variability. A non-parametric Gaussian Process (GP) model was employed to capture the time- and weather-dependent characteristics of the heating load while explicitly modeling its uncertainty. The model's effectiveness is demonstrated through strong predictive performance on unseen data and physically interpretable insights into load behavior. This data-driven, probabilistic load profile serves as a robust and flexible input for subsequent system simulations, enabling a more confident and statistically sound analysis of retrofit potential.

Ham, S W

Test Data Analysis of the Thermodynamic Vent System-Augmented Top Spray Injector Liquid Nitrogen Transfer Experiments

Traditionally, a cryogenic tank must be pre-chilled to some “target” temperature before the main vent valve can be closed to attempt a non-vented fill (NVF) of cryogenic liquid propellant. This methodology is particularly attractive for performing in-space transfer of cryogens due to the unknown location of the liquid/vapor interface in microgravity and the high likelihood of venting liquid if the vent valve is opened during transfer. This paper presents in-depth test data analysis of a Thermodynamic Vent System (TVS) augmented injector used for cryogenic tank chilldown and fill experiments of a thin-walled Titanium tank. Eight tests were conducted using liquid nitrogen across a range of inlet conditions and boundary conditions, and three different chilldown/fill methods. For four of the tests, the injector sprays liquid into the tank as normal, but also uses a TVS heat exchanger to cool the metallic injector itself as well as the main incoming liquid stream. Results show that using the TVS augmented injector simplifies transfer operation via enhanced condensation at the injector surface at the cost of sacrificing only a small amount of propellant.

Thermodynamic Vent System

SuperHERO X-ray Telescope Balloon Mission

The SuperHERO hard X-ray telescope will observe extended sources at unprecedented imaging resolution, revealing the origins of non-thermal emission in extreme astrophysical environments. Employing NASA Marshall Space Flight Center’s replicated NiCo full-shell X-ray mirror technology to achieve better than 10-arcsecond half-power diameter (HPD) angular resolution on a balloon platform, SuperHERO consists of seven identical, co-aligned telescopes with a combined effective area of 45 cm² at 30 keV. SuperHERO was selected in the APRA 2023 call as a five-year mission, with an inaugural two-day flight scheduled from Fort Sumner, New Mexico, in the fall of 2028. For its initial target, SuperHERO will observe the Crab pulsar wind nebula, localizing hard X-ray emission. Mission development is currently underway, with updates and status to be provided in this presentation.

Nick Thomas

High-Throughput Uniformity and Defect Monitoring in Low-Temperature Electrolysis Porous Transport Layers Using X-Ray Radiography

Effective quality control (QC) for manufacturing proton exchange membrane water electrolysis (PEMWE) components is critical to enabling widespread adoption of the technology for hydrogen generation. This study investigates X-ray radiography as a novel, high-throughput, potentially in-line QC technique for detecting defects and assessing material property distributions in titanium-based porous transport layers (PTLs) which constitute a crucial component of low temperature PEMWE stacks. We obtain radiographs of a set of fifteen PTLs and model their absorbance of the broadband radiation as a second-order polynomial to account for the non-monoenergetic radiation source used in this study. The resulting model serves as a basis for predicting the areal density and porosity distributions of the PTLs. We find radiography successful in detecting multiple instances of defects, including holes/depressions, cracks, and excess material on the surface or in the pores of the material, demonstrating its potential as a robust in-line QC tool for PTL manufacturing.

08 HYDROGEN

Nuclear Physics Made Very, Very Easy

The fundamental approach to nuclear physics was prepared to introduce basic reactor principles to various groups of non-nuclear technical personnel associated with NERVA Test Operations. NERVA Test Operations functions as the field test group for the Nuclear Rocket Engine Program. Nuclear Engine for Rocket Vehicle Application (NERVA) program is the combined efforts of Aerojet-General Corporation as prime contractor, and Westinghouse Astronuclear Laboratory as the major subcontractor, for the assembly and testing of nuclear rocket engines. Development of the NERVA Program is under the direction of the Space Nuclear Propulsion Office, a joint agency of the U. S. Atomic Energy Commission and the National Aeronautics and Space Administration. This report is being reprinted for use in the U. S. Atomic Energy Commission and National Aeronautics and Space Administration educational and technology utilization programs.

22 GENERAL STUDIES OF NUCLEAR REACTORS

SODAs: sparse optimization for the discovery of differential and algebraic equations

Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.

DAE