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

Correspondence between discrete and continuous models of excitable media: trigger waves

We present a theoretical framework for relating continuous partial differential equation (PDE) models of excitable media to discrete cellular automata (CA) models on a randomized lattice. These relations establish a quantitative link between the CA model and the specific physical system under study. We derive expressions for the CA model's plane wave speed, critical curvature, and effective diffusion constant in terms of the model's internal parameters (the interaction radius, excitation threshold, and time step). We then equate these expressions to the corresponding quantities obtained from solution of the PDEs (for a fixed excitability). This yields a set of coupled equations with a unique solution for the required CA parameter values. Here we restrict our analysis to "trigger" wave solutions obtained in the limiting case of a two-dimensional excitable medium with no recovery processes. We tested the correspondence between our CA model and two PDE models (the FitzHugh-Nagumo medium and a medium with a "sawtooth" nonlinear reaction source) and found good agreement with the numerical solutions of the PDEs. Our results suggest that the behavior of trigger waves is actually controlled by a small number of parameters.

Non-NASA Center↗
NASA NTRS2004

Control of gravitropic orientation. I. Non-vertical orientation by primary roots of maize results from decay of competence for orthogravitropic induction

Plant organs may respond to gravity by vertical (orthogravitropic), oblique (plagiogravitropic) or horizontal (diagravitropic) growth. Primary roots of maize (Zea mays L.) provide a good system for studying such behaviours because they are reportedly capable of displaying all three responses. In current work using maize seedlings of the Silver Queen cultivar, stabilisation of growth at an oblique orientation was commonplace. Hypothetically, plagiogravitropism may be accomplished either by a process we call graded orthogravitropism or by hunting about a sensed non-vertical setpoint. In graded orthotropism primary bending is unidirectional and depends on facilitative stimuli that determine its extent. The hallmark of the setpoint mechanism is restorative curvature of either sign following a displacement; both diagravitropism and orthogravitropism are based on setpoints. Roots settled in a plagiogravitropic orientation were tested with various illumination and displacement protocols designed to distinguish between these two hypotheses. The tests refuted the setpoint hypothesis and supported that of graded orthotropism. No evidence of diagravitropism could be found, thus, earlier claims were likely based on inadequately controlled observations of graded orthotropism. We propose that orthotropism is graded by the sequential action of dual gravity receptors: induction of a vectorial gravitropic response requires gravitational induction of a separate facilitative response, whose decay in the absence of fresh stimuli can brake gravitropism at plagiotropic angles.

NASA Discipline Plant Biology↗
DOE OSTI2023

Catalyst: Fast and flexible modeling of reaction networks

We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst’s broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Control of the optical properties of gum and dentin tissue of a human tooth at laser spectral lines in the range of 200 – 800 nm

A change in the optical properties of the tissue of the gums and dentin of a human tooth is experimentally investigated in vitro with the application of a 99.5 % glycerol solution. The kinetics of the total transmittance in the spectral range of 200 – 800 nm on characteristic laser lines is studied. It is shown that the effectiveness of optical clearing, expressed as a relative increase in the total transmittance, for gum tissue in the laser radiation range of 200 – 800 nm can vary significantly (by an order of magnitude). In this case, the absolute values of optical transmittance in the UV region are relatively small due to the strong absorption by the endogenous chromophores of the studied tissues, and in the visible and NIR regions, on the contrary, they are significant, which is due to the absence of strongly absorbing chromophores. The effectiveness of dentin clearing in the entire visible and NIR spectral regions is on average 25 % for 300 minutes of immersion of the sample in a 99.5 % glycerol solution. The presented results are necessary for increasing the efficiency of laser therapy and methods for optical monitoring of morphological changes in biological tissues for the purpose of early diagnosis of diseases. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Isotopically nonstationary metabolic flux analysis of plants: recent progress and future opportunities

Metabolic flux analysis (MFA) is a valuable tool for quantifying cellular phenotypes and to guide plant metabolic engineering. By introducing stable isotopic tracers and employing mathematical models, MFA can quantify the rates of metabolic reactions through biochemical pathways. Recent applications of isotopically nonstationary MFA (INST–MFA) to plants have elucidated nonintuitive metabolism in leaves under optimal and stress conditions, described coupled fluxes for fast–growing algae, and produced a synergistic multi–organ flux map that is a first in MFA for any biological system. Finally, these insights could not be elucidated through other approaches and show the potential of INST–MFA to correct an oversimplified understanding of plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗
NASA NTRS1996

Noninvasive Determination of Bone Mechanical Properties using Vibration Response: A Refined Model and Validation in vivo

Accurate non-invasive mechanical measurement of long bones is made difficult by the masking effect of surrounding soft tissues. Mechanical Response Tissue Analysis (MRTA) offers a method for separating the effects of the soft tissue and bone; however, a direct validation has been lacking. A theoretical analysis of wave propagation through the compressed tissue revealed a strong mass effect dependent on the relative accelerations of the probe and bone. The previous mathematical model of the bone and overlying tissue system was reconfigured to incorporate the theoretical finding. This newer model (six-parameter) was used to interpret results using MRTA to determine bone cross-sectional bending stiffness, EI(sub MRTA). The relationship between EI(MRTA) and theoretical EI values for padded aluminum rods was R(exp 2) = 0.999. A biological validation followed using monkey tibias. Each bone was tested in vivo with the MRTA instrument. Postmortem, the same tibias were excised and tested to failure in three-point bending to determine EI(sub 3-PT) and maximum load. Diaphyseal Bone Mineral Density (BMD) measurements were also made. The relationship between E(sub 3-PT) and in vivo EI(sub MRTA) using the six-parameter model is strong (R(exp 2) = 0.947) and better than that using the older model (R(exp 2) = 0.645). EI(MRTA) and BMD are also highly correlated (R(exp 2) = 0.853). MRTA measurements in vivo and BMD ex vivo are both good predictors of scaled maximum strength (R(exp 2) = 0.915 and R(exp 2) = 0.894, respectively). This is the first biological validation of a non-invasive mechanical measurement of bone by comparison to actual values. The MRTA technique has potential clinical value for assessing long-bone mechanical properties.

Roberts, S. G.↗
NASA NTRS1996

Noninvasive Determination of Bone Mechanical Properties Using Vibration Response: A Refined Model and Validation in vivo

Accurate non-invasive mechanical measurement of long bones is made difficult by the masking effect of surrounding soft tissues. Mechanical response tissue analysis (MRTA) offers a method for separating the effects of the soft tissue and bone; however, a direct validation has been lacking. A theoretical analysis of wave propagation through the compressed tissue revealed a strong mass effect dependent on the relative accelerations of the probe and bone. The previous mathematical model of the bone and overlying tissue system was reconfigured to incorporate the theoretical finding. This newer model (six-parameter) was used to interpret results using MRTA to determine bone cross-sectional bending stiffness, EI(sub MRTA). The relationship between EI(sub MRTA) and theoretical EI values for padded aluminum rods was R(sup 2) = 0.999. A biological validation followed using monkey tibias. Each bone was tested in vivo with the MRTA instrument. Postmortem, the same tibias were excised and tested to failure in three-point bending to determine EI(sub 3-PT) and maximum load. Diaphyseal bone mineral density (BMD) measurements were also made. The relationship between EI(sub 3-PT) and in vivo EI(sub MRTA) using the six-parameter model is strong (R(sup 2) = 0.947) and better than that using the older model (R(sup 2) = 0.645). EI(sub MRTA) and BMD are also highly correlated (R(sup 2) = 0.853). MRTA measurements in vivo and BMD ex vivo are both good predictors of scaled maximum strength (R(sup 2) = 0.915 and R(sup 2) = 0.894, respectively). This is the first biological validation of a non- invasive mechanical measurement of bone by comparison to actual values. The MRTA technique has potential clinical value for assessing long-bone mechanical properties.

Roberts, S. G.↗
NASA NTRS1994

The excitability of plant cells: with a special emphasis on characean internodal cells

This review describes the basic principles of electrophysiology using the generation of an action potential in characean internodal cells as a pedagogical tool. Electrophysiology has proven to be a powerful tool in understanding animal physiology and development, yet it has been virtually neglected in the study of plant physiology and development. This review is, in essence, a written account of my personal journey over the past five years to understand the basic principles of electrophysiology so that I can apply them to the study of plant physiology and development. My formal background is in classical botany and cell biology. I have learned electrophysiology by reading many books on physics written for the lay person and by talking informally with many patient biophysicists. I have written this review for the botanist who is unfamiliar with the basics of membrane biology but would like to know that she or he can become familiar with the latest information without much effort. I also wrote it for the neurophysiologist who is proficient in membrane biology but knows little about plant biology (but may want to teach one lecture on "plant action potentials"). And lastly, I wrote this for people interested in the history of science and how the studies of electrical and chemical communication in physiology and development progressed in the botanical and zoological disciplines.

NASA Discipline Plant Biology↗
NASA NTRS2009

System and method for deriving a process-based specification

A system and method for deriving a process-based specification for a system is disclosed. The process-based specification is mathematically inferred from a trace-based specification. The trace-based specification is derived from a non-empty set of traces or natural language scenarios. The process-based specification is mathematically equivalent to the trace-based specification. Code is generated, if applicable, from the process-based specification. A process, or phases of a process, using the features disclosed can be reversed and repeated to allow for an interactive development and modification of legacy systems. The process is applicable to any class of system, including, but not limited to, biological and physical systems, electrical and electro-mechanical systems in addition to software, hardware and hybrid hardware-software systems.

Hinchey, Michael Gerard↗
DOE OSTI2021

A Mathematical Model of the Dynamics of Cytokine Expression and Human Immune Cell Activation in Response to the Pathogen Staphylococcus aureus

Cell-based mathematical models have previously been developed to simulate the immune system in response to pathogens. Mathematical modeling papers which study the human immune response to pathogens have predicted concentrations of a variety of cells, including activated and resting macrophages, plasma cells, and antibodies. This study aims to create a comprehensive mathematical model that can predict cytokine levels in response to a gram-positive bacterium, S. aureus by coupling previous models. To accomplish this, the cytokines Tumor Necrosis Factor Alpha (TNF-α), Interleukin 6 (IL-6), Interleukin 8 (IL-8), and Interleukin 10 (IL-10) are included to quantify the relationship between cytokine release from macrophages and the concentration of the pathogen, S. aureus, ex vivo. Partial differential equations (PDEs) are used to model cellular response and ordinary differential equations (ODEs) are used to model cytokine response, and interactions between both components produce a more robust and more complete systems-level understanding of immune activation. In the coupled cellular and cytokine model outlined in this paper, a low concentration of S. aureus is used to stimulate the measured cellular response and cytokine expression. Results show that our cellular activation and cytokine expression model characterizing septic conditions can predict ex vivo mechanisms in response to gram-negative and gram-positive bacteria. Our simulations provide new insights into how the human immune system responds to infections from different pathogens. Novel applications of these insights help in the development of more powerful tools and protocols in infection biology.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1974

An interdisciplinary study of the estuarine and coastal oceanography of Block Island Sound and adjacent New York coastal waters

The author has identified the following significant results. Photo-optical additive color quantitative measurements were made of ERTS-1 reprocessed positives of New York Bight and Block Island Sound. Regression of these data on almost simultaneous ship sample data of water's physical, chemical, biological, and optical properties showed that ERTS bands 5 and 6 can be used to predict the absolute value of the total number of particles and bands 4 and 5 to predict the relative extinction coefficient in New York Bight. Water masses and mixing patterns in Block Island Sound heretofore considered transient were found to be persistent phenomena requiring revision of existing mathematical and hydraulic models.

Yost, E. F.↗
DOE OSTI2022

Random Knotting in Fractal Ring Polymers

Many ring polymer systems of physical and biological interest exhibit both pronounced topological effects and nontrivial self-similarity, but the relationship between these two phenomena has not yet been clearly established. Here, we use theory and simulation to formulate such a connection by studying a fundamental topological property—the random knotting probability—for ring polymers with varying fractal dimension, d f . Using straightforward scaling arguments, we generalize a classic mathematical result, showing that the probability of a trivial knot decays exponentially with chain size, N, for all fractal dimensions: P 0 (N) ∝ exp(–N/N 0 ). However, no such simple considerations can account for the dependence of the knotting length, N 0 , on d f , necessitating a more involved analytical calculation. This analysis reveals a complicated double-exponential dependence, which is well supported by numerical data. By contrast, functional forms typical of simple scaling theories fail to adequately describe the observations. These findings are equally valid for two-dimensional ring polymer systems, where “knotting” is defined as the intersection of any two segments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗