Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematical biology”

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 757 records · Page 42

Structural inference of networked dynamical systems with universal differential equations

Networked dynamical systems are common throughout science in engineering; e.g., biological networks, reaction networks, power systems, and the like. For many such systems, nonlinearity drives populations of identical (or near-identical) units to exhibit a wide range of nontrivial behaviors, such as the emergence of coherent structures (e.g., waves and patterns) or otherwise notable dynamics (e.g., synchrony and chaos). Here, we seek to infer (i) the intrinsic physics of a base unit of a population, (ii) the underlying graphical structure shared between units, and (iii) the coupling physics of a given networked dynamical system given observations of nodal states. These tasks are formulated around the notion of the Universal Differential Equation, whereby unknown dynamical systems can be approximated with neural networks, mathematical terms known a priori (albeit with unknown parameterizations), or combinations of the two. We demonstrate the value of these inference tasks by investigating not only future state predictions but also the inference of system behavior on varied network topologies. The effectiveness and utility of these methods are shown with their application to canonical networked nonlinear coupled oscillators.

97 MATHEMATICS AND COMPUTING↗

Regional Relative Risk, a Physics-Based Metric for Characterizing Airborne Infectious Disease Transmission

Airborne infectious disease transmission events occur over a wide range of spatial scales and can be an important means of disease transmission. Physics- and biology-based models can assist in predicting airborne transmission events, overall disease incidence, and disease control strategy efficacy. We describe a new theory that extends current approaches for the case in which an individual is infected by a single airborne particle, including the scenario in which numerous infectious particles are present in the air but only one causes infection. A single infectious particle can contain more than one pathogenic microorganism and be physically larger than the pathogen itself. This approach allows robust relative risk estimates even when there is wide variation in (i) individual exposures and (ii) the individual response to that exposure (the pathogen dose-response function can take any mathematical form and vary by individual). Based on this theory, we propose the regional relative risk—a new metric, distinct from the traditional relative risk metric, that compares the risk between two regions. In theory, these regions can range from individual rooms to large geographic areas. In this paper, we apply the regional relative risk metric to outdoor disease transmission events over spatial scales ranging from 50 m to 20 km, demonstrating that in many common cases minimal input information is required to use the metric. Also, we demonstrate that the model predictions are consistent with data from prior outbreaks. Future efforts could apply and validate this theory for other spatial scales, such as transmission within indoor environments. This work provides context for (i) the initial stages of an airborne disease outbreak and (ii) larger-scale disease spread, including unexpected low-probability disease “sparks” that potentially affect remote populations, a key practical issue in controlling airborne disease outbreaks.

54 ENVIRONMENTAL SCIENCES↗

Fully energy-dependent HZETRN (a galactic cosmic-ray transport code)

For extended manned space missions, the radiation shielding design requires efficient and accurate cosmic-ray transport codes that can handle the physics processes in detail. The Langley Research Center galactic cosmic-ray transport code (HZETRN) is currently under development for such design use. The cross sections for the production of secondary nucleons in the existing HZETRN code are energy dependent only for nucleon collisions. The approximation of energy-independent, heavy-ion fragmentation cross section is now removed by implementing a mathematically simplified energy-dependent stepping formalism for heavy ions. The cross section at each computational grid is obtained by linear interpolation from a few tabulated data to minimize computing time. Test runs were made for galactic cosmic-ray transport through a liquid hydrogen shield and a water shield at solar minimum. The results show no appreciable change in total fluxes or computing time compared with energy-independent calculations. Differences in high LET (linear energy transfer) spectra are noted, however, because of the large variation in cross sections at the low-energy region. The high LET components are significantly higher in the new code and have important implications on biological risk estimates for heavy-ion exposure.

Shinn, Judy L.↗

NASA/MSFC/NSSTC Science Communication Roundtable

For the last several years the Science Directorate at Marshall Space Flight Center has carried out a diverse program of Internet-based science communication. The Directorate's Science Roundtable includes active researchers, NASA public relations, educators, and administrators. The Science@NASA award-winning family of Web sites features science, mathematics, and space news. The program includes extended stories about NASA science, a curriculum resource for teachers tied to national education standards, on-line activities for students, and webcasts of real-time events. Science stories cover a variety of space-related subjects and are expressed in simple terms everyone can understand. The sites address such questions as: what is space weather, what's in the heart of a hurricane, can humans live on Mars, and what is it like to live aboard the International Space Station? Along with a new look, the new format now offers articles organized by subject matter, such as astronomy, living in space, earth science or biology. The focus of sharing real-time science related events has been to involve and excite students and the public about science. Events have involved meteor showers, solar eclipses, natural very low frequency radio emissions, and amateur balloon flights. In some cases broadcasts accommodate active feedback and questions from Internet participants. Information will be provided about each member of the Science@NASA web sites.

Adams, Mitzi L.↗

Using emergent order to shape a space society

A fast-growing movement in the scientific community is reshaping the way that we view the world around us. The short-hand name for this movement is 'chaos'. Chaos is a science of the global, nonlinear nature of systems. The center of this set of ideas is that simple, deterministic systems can breed complexity. Systems as complex as the human body, ecology, the mind or a human society. While it is true that simple laws can breed complexity, the other side is that complex systems can breed order. It is the latter that I will focus on in this paper. In the past, nonlinear was nearly synonymous with unsolvable because no general analytic solutions exist. Mathematically, an essential difference exists between linear and nonlinear systems. For linear systems, you just break up the complicated system into many simple pieces and patch together the separated solutions for each piece to form a solution to the full problem. In contrast, solutions to a nonlinear system cannot be added to form a new solution. The system must be treated in its full complexity. While it is true that no general analytical approach exists for reducing a complex system such as a society, it can be modeled. The technical involves a mathematical construct called phase space. In this space stable structures can appear which I use as analogies for the stable structures that appear in a complex system such as an ecology, the mind or a society. The common denominator in all of these systems is that they rely on a process called feedback loops. Feedback loops link the microscopic (individual) parts to the macroscopic (global) parts. The key, then, in shaping a space society, is in effectively using feedback loops. This paper will illustrate how one can model a space society by using methods that chaoticists have developed over the last hundred years. And I will show that common threads exist in the modeling of biological, economical, philosophical, and sociological systems.

Graps, Amara L.↗

50 Years of metabolic control analysis

It has been 50 years since the publication of the seminal works by Savageau, Kacser & Burns and Heinrich & Rapoport, on the control of biochemical networks. What is remarkable about these three publications (especially the two papers by Kacser et al. and Rapoport et al.) is that they independently developed almost the same approach and came to the same conclusions. This approach is today called metabolic control analysis (MCA), and is a body of work that uses mathematics to reason about the properties of biochemical pathways in a deductive manner. Although originally focused on metabolic pathways, it was realized early on that it could be equally applied to genetic and signalling networks. In fact anything that one could associate a stoichiometry matrix with. In recent years, others have reinvented MCA, using different notation, highlighting the fact that there has been a good deal of conceptual convergence in this area.

59 BASIC BIOLOGICAL SCIENCES↗

Stereology techniques in radiation biology

Clinicians involved in conventional radiation therapy are very concerned about the dose-response relationships of normal tissues. Before proceeding to new clinical protocols, radiation biologists involved with conformal proton therapy believe it is necessary to quantify the dose response and tolerance of the organs and tissues that will be irradiated. An important focus is on the vasculature. This presentation reviews the methodology and format of using confocal microscopy and stereological methods to quantify tissue parameters, cell number, tissue volume and surface area, and vessel length using the microvasculature as a model tissue. Stereological methods and their concepts are illustrated using an ongoing study of the dose response of the microvessels in proton-irradiated hemibrain. Methods for estimating the volume of the brain and the brain cortex, the total number of endothelial cells in cortical microvessels, the length of cortical microvessels, and the total surface area of cortical microvessel walls are presented step by step in a way understandable for readers with little mathematical background. It is shown that stereological techniques, based on a sound theoretical basis, are powerful and reliable and have been used successfully.

Non-NASA Center↗

Nonlinear model of infection wavy oscillation of COVID-19 in Japan based on diffusion kinetics

The infectious propagation of SARS-CoV-2 is continuing worldwide, and specifically, Japan is facing severe circumstances. Medical resource maintenance and action limitations remain the central measures. An analysis of long-term follow-up reports in Japan shows that the infection number follows a unique wavy oscillation, increasing and decreasing over time. However, only a few studies explain the infection wavy oscillation. This study introduces a novel nonlinear mathematical model of the new infection wavy oscillation by applying the macromolecule diffusion theory. In this model, the diffusion coefficient that depends on population density gives nonlinearity in infection propagation. As a result, our model accurately simulated infection wavy oscillations, and the infection wavy oscillation frequency and amplitude were closely linked with the recovery rate of infected individuals. In conclusion, our model provides a novel nonlinear contact infection analysis framework.

60 APPLIED LIFE SCIENCES↗

Generative diffusion model surrogates for mechanistic agent-based biological models

Mechanistic, multicellular, agent-based models are commonly used to investigate tissue, organ, and organism-scale biology at single-cell resolution. The Cellular-Potts Model (CPM) is a powerful and popular framework for developing and interrogating these models. CPMs become computationally expensive at large space- and time- scales making application and investigation of developed models difficult. Surrogate models may allow for the accelerated evaluation of CPMs of complex biological systems. However, the stochastic nature of these models means each set of parameters may give rise to different model configurations, complicating surrogate model development. In this work, we leverage denoising diffusion probabilistic models (DDPMs) to train a generative AI surrogate of a CPM used to investigate in vitro vasculogenesis. We describe the use of an image classifier to learn the characteristics that define unique areas of a 2-dimensional parameter space. We then apply this classifier to aid in surrogate model selection and verification. Our CPM model surrogate generates model configurations 20,000 timesteps ahead of a reference configuration and demonstrates approximately a 22x reduction in computational time as compared to native code execution. Our work represents a step towards the implementation of DDPMs to develop digital twins of stochastic biological systems.

97 MATHEMATICS AND COMPUTING↗

Position Papers for the 2024 ASCR Workshop on Analog Computing for Science

Analog computing capabilities have existed since the dawn of science, but modern techniques potentially allow for the construction, verification, and characterization of complicated analog computing systems in a wide variety of contexts, from high-performance computational accelerators to nanorobotics and synthetic biology. In short, advances in analog computing can enable the creation of physical systems with complex behaviors that meet sophisticated requirements. Meeting our nation's needs, from needs in computing and modeling, to needs for advanced materials and energy technologies, continues to motivate pursuing novel kinds of complex systems, and thus analog computing techniques, in all of these spaces.

97 MATHEMATICS AND COMPUTING↗

Uncovering heterogeneous intercommunity disease transmission from neutral allele frequency time series

The COVID-19 pandemic has underscored the need for accurate epidemic forecasting to predict pathogen spread, evolution, and evaluate intervention strategies. Forecast reliability hinges on detailed knowledge of disease transmission across population segments, which may be inferred from contact surveys or mobility data. However, these indirect approaches make it difficult to estimate rare transmissions between socially or geographically distant communities. We show that the steep ramp-up of genome sequencing surveillance during the pandemic can be leveraged to directly identify transmission patterns between geographically defined communities. Our approach uses a hidden Markov model to infer the fraction of infections a community imports from others based on how rapidly allele frequencies in the focal community converge to those in the donor communities. Applying this method to SARS-CoV-2 sequencing data from England and the United States, we uncover networks of intercommunity transmission that reflect geographical relationships while exposing significant long-range interactions. The scaling of importation rate with distance is consistent across both countries, yet weaker than expected based on mobility data, highlighting limitations of indirect inference. We show that transmission patterns can change between waves of variants of concern and analyze how the inferred heterogeneity in intercommunity transmission impacts evolutionary forecasts. While applied here to geographically defined communities, our approach could be applied to those defined by other traits (e.g., age, socioeconomic status), provided time-series data can be stratified accordingly. Overall, our study highlights population genomic time series data as a crucial record of epidemiological interactions, which can be deciphered using tree-free inference methods.

Okada, Takashi [Department of Physics; University ↗

Surrogate modeling of Cellular-Potts agent-based models as a segmentation task using the U-Net neural network architecture

The Cellular-Potts model is a powerful and ubiquitous framework for developing computational models for simulating complex multicellular biological systems. Cellular-Potts models (CPMs) are often computationally expensive due to the explicit modeling of interactions among large numbers of individual model agents and diffusive fields described by partial differential equations (PDEs). In this work, we develop a convolutional neural network (CNN) surrogate model using a U-Net architecture that accounts for periodic boundary conditions. We use this model to accelerate the evaluation of a mechanistic CPM previously used to investigate in vitro vasculogenesis. The surrogate model was trained to predict 100 computational steps ahead (Monte-Carlo steps, MCS), accelerating simulation evaluations by a factor of 562 times compared to single-core CPM code execution on CPU. Over short timescales of up to 3 recursive evaluations, or 300 MCS, our model captures the emergent behaviors demonstrated by the original Cellular-Potts model such as vessel sprouting, extension and anastomosis, and contraction of vascular lacunae. This approach demonstrates the potential for deep learning to serve as a step toward efficient surrogate models for CPM simulations, enabling faster evaluation of computationally expensive CPM simulations of biological processes.

97 MATHEMATICS AND COMPUTING↗

A Modular System for Increasing Predictiveness for Extreme Climate Predictions

We know that climate change is poised to reshape our world, but we lack clear enough predictions about precisely how. The preponderance of these changes is associated with human activity, specifically the emission of CO 2 and other greenhouse gases. Problematically, projections of climate change continue to be marred by unacceptably large uncertainties which hamper informed decision-making and cost society a chance to adapt proactively and effectively. These uncertainties stem from deficiencies in predictions of future greenhouse gas emissions, but also from inaccuracies in the representation of the physical models used to predict the climate response to such emissions. The uncertainties in projections associated with the inaccurate representation of climate physics, chemistry and biology are similar to those that plagued the first global climate models developed fifty years ago, despite more than a factor 10 8 increase in computer performance. Our transformational question is then, how can the accuracy of climate projections be dramatically improved by applying recent advances in the computational and data sciences to train the models with the wealth of data being constantly collected about the ongoing changes in the climate system?

54 ENVIRONMENTAL SCIENCES↗

Land Surface Modeling 2.0 for agricultural climate change impact assessments

Focus Area: This white paper addresses DOE AI4ESP focus area #2 by providing a sketch blueprint for a next generation, hybrid AI/process-based global Land Surface Modeling (LSM) framework to improve projections of climate change impacts on the land surface system including agriculture and the hydrological cycle. Science Challenge: Climate change impacts on agriculture are highly uncertain: -50% to +150% global production changes for major grains under high-end climate change, for example. The core challenge in land surface system predictability is the large gap between the scale at which the relevant biological processes act and the scale at which the risks need to be assessed. Conventional empirical and process-based modeling approaches are insufficient. A new multi-scale modeling paradigm which employs AI/ML to learn from new streams of remote sensing data and targeted ’gene-to-global’ simulations is needed.

54 ENVIRONMENTAL SCIENCES↗

Considering Cell Proliferation to Optimize Detection of Radiation-induced 53BP1+ Foci in 15 Mouse Strains ex vivo

Due to high metabolic activity, proliferating cells continuously generate free radicals, which induce DNA double strand breaks (DSB). Fluorescently tagged nuclear foci of DNA repair protein 53 binding protein-1(53BP1) are used as a standard metric for measuring DSB formation at baseline and in response to environmental insults such as radiation. Here we demonstrate that the background level of spontaneous 53BP1+ foci formation can be modeled mathematically as a function of cell confluence, which is a metric of proliferation rate. This model was validated using spontaneous 53BP1+ foci data from 68 cultures of primary skin fibroblasts derived from 15 different strains of mice, showing a ~10 fold decrease from low to full confluence that is independent of mouse strain. We developed an online open access tool to correct for the impact of cell confluence on the detection of radiation-induced 53BP1+ foci (RIF). This tool provides guidelines for the number of cells required to reach statistical significance for the detection of excess foci induced by low doses of ionizing radiation as a function of confluence and time post-irradiation. We hope that this quantification tool will help the radiation biology community in the design of future experiments that utilize 53BP1+ foci-based quantification of radiation responses in vitro. Our “tool for enhanced results of RIF in cells” (terRIFic) can be found at: https://radbiolab.shinyapps.io/terrific/

Radiation, DNA damage, confluence↗

Development of Hydropower Biological Evaluation Toolset (HBET): V2.1.9 Release Notes for HBET

The following release notes reflect changes made to HBET for proposed changes to be released in July 2024. Notes are broken up into three sections: 1) Key Improvements, 2) Bug Fixes, and 3) Data Changes • Key Improvements: primary features added and changes to existing features that affect the user experience. • Bug Fixes: Issues discovered or reported that were fixed in the proposed work to be released. • Data Changes: Any work done on the databases directly or the process to calculate data for the system.

13 HYDRO ENERGY↗

Statistical relationships across epigenomes using large-scale hierarchical clustering

Recent advances in genomics and sequencing platforms have revolutionized our ability to create immense data sets, particularly for studying epigenetic regulation of gene expression. However, the avalanche of epigenomic data is difficult to parse for biological interpretation given nonlinear complex patterns and relationships. This attractive challenge in epigenomic data lends itself to machine learning for discerning infectivity and susceptibility. In this study, we explore over 3000 epigenomes of uninfected individuals and provide a framework to characterize the relationships among epigenetic modifiers, their modifiers, genetic loci, and specific immune cell types across all chromosomes using hierarchical clustering. Hierarchical clustering of epigenomic data revealed consistent epigenetic patterns across chromosomes, demonstrating that variation due to epigenetic modifiers is greater than variation between cell types. Gene Ontology and KEGG pathway analyses indicated significant enrichment of genes involved in chromatin remodeling, mRNA splicing, immune responses, and the regulation of microRNAs and snoRNAs. Epigenetic modifiers frequently formed biologically relevant clusters, including the cohesin complex, RNA Polymerase II transcription factors, and PRC2 complex members. These clustering behaviors remained consistent across all chromosomes, supported by entropy analysis and high Adjusted Rand Index scores, indicating robust cross-chromosomal similarity. Co-occurrence analysis further revealed specific sets of modifiers that consistently appeared together within clusters, reflecting shared biological functions and interactions. Validation using another dataset confirmed the reproducibility of these clustering patterns and modifier co-occurrence relationships, underscoring the reliability and generalizability of the methodology.

97 MATHEMATICS AND COMPUTING↗

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Furthermore, our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

97 MATHEMATICS AND COMPUTING↗