Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional rate functions”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Variation in hydroclimate sustains tropical forest biomass and promotes functional diversity

The fate of tropical forests under climate change is unclear as a result, in part, of the uncertainty in projected changes in precipitation and in the ability of vegetation models to capture the effects of drought-induced mortality on aboveground biomass (AGB). In this work, we evaluated the ability of a terrestrial biosphere model with demography and hydrodynamics (Ecosystem Demography, ED2-hydro) to simulate AGB and mortality of four tropical tree plant functional types (PFTs) that operate along light- and water-use axes. Model predictions were compared with observations of canopy trees at Barro Colorado Island (BCI), Panama. We then assessed the implications of eight hypothetical precipitation scenarios, including increased annual precipitation, reduced inter-annual variation, El Niño-related droughts and drier wet or dry seasons, on AGB and functional diversity of the model forest. When forced with observed meteorology, ED2-hydro predictions capture multiple BCI benchmarks. ED2-hydro predicts that AGB will be sustained under lower rainfall via shifts in the functional composition of the forest, except under the drier dry-season scenario. Furthermore these results support the hypothesis that inter-annual variation in mean and seasonal precipitation promotes the coexistence of functionally diverse PFTs because of the relative differences in mortality rates. If the hydroclimate becomes chronically drier or wetter, functional evenness related to drought tolerance may decline.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Global land model development: time to shift from a plant functional type to a plant functional trait approach (Final Report)

This project advanced global land models by shifting from the current plant functional type approach to one that better utilizes what is known about the importance and variability of plant traits, within a framework of simultaneously improving fundamental physiological relations that are at the core of model carbon cycling algorithms. A primary goal for Earth system modeling is to make accurate predictions of the future trajectory of the climate system, based on a mechanistic understanding of processes regulating fluxes of mass and energy among system components. Land plays an important role in modifying the Earth’s mass and energy balance, as a critical link in the global cycling of carbon, among others. Land surface models have developed to include mechanistic representations of vegetation physiology, carbon and nutrient dynamics in plants and soils, how they might respond to changing climate and chemistry, and how those changes might feedback to influence changes in atmospheric greenhouse gases themselves. Existing models represent the global distribution of vegetation types using the Plant Functional Type concept. Plant Functional Types are classes of plant species with similar evolutionary and life history with presumably similar responses to environmental conditions like CO 2 , water and nutrient availability. Fixed properties for each Plant Functional Type are specified through a collection of physiological parameters, or traits. These traits, mostly physiological in nature (e.g., leaf nitrogen and longevity) are used in model algorithms to estimate ecosystem properties and/or drive calculated process rates. In most models, 5 to 15 functional types represent terrestrial vegetation; in essence, they assume there are a total of only 5 to 15 different kinds of plants on the entire globe. This assumption of constant plant traits captured within the functional type concept has serious limitations, as a single set of traits does not reflect trait variation observed within and between species and communities. While this simplification was necessary decades past, substantial improvement is now possible. Rather than assigning a small number of constant parameter values to all grid cells in a model, we developed procedures that predict a frequency distribution of values for any given grid cell. Thus, by incorporating the mean and variance model fidelity has been improved. The trait-based approach has improved land modeling by (1) incorporating patterns and heterogeneity of traits into model parameterization, thus evolving away from a framework that considers large areas of vegetation to have near identical trait values and (2) allowing for improved treatment of physiological responses to environmental drivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

SCALE 6.2.4 Validation: Radiation Shielding

For safe and reliable use of computer codes by the community, accuracy must be clearly evaluated. In particular, the nuclear reactor engineering and licensing field needs accurate tools for radiation shielding modeling. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is one such tool, with built-in variance reduction methods distributed within the SCALE code, and its validity is demonstrated in this report for the released version 6.2.4. Representative benchmarks corresponding to shielding analysis are selected for the validation study. Typical experimental results analyzed from those benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma doses, foil neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between experiment and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, MAVRIC agrees well with the experiment results. MAVRIC is also compared to Monte Carlo N-Particle (MCNP) calculations when available, and both codes generally produce good agreements within estimated uncertainties. The selected benchmarks are obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook (ICSBEP Handbook), the Shielding Integral Benchmark Archive & Database (SINBAD), and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to incrementally expand the shielding validation suite incrementally.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1973

Investigations of respiratory control systems simulation

The Grodins' respiratory control model was investigated and it was determined that the following modifications were necessary before the model would be adaptable for current research efforts: (1) the controller equation must be modified to allow for integration of the respiratory system model with other physiological systems; (2) the system must be more closely correlated to the salient physiological functionings; (3) the respiratory frequency and the heart rate should be expanded to illustrate other physiological relationships and dependencies; and (4) the model should be adapted to particular individuals through a better defined set of initial parameter values in addition to relating these parameter values to the desired environmental conditions. Several of Milhorn's respiratory control models were also investigated in hopes of using some of their features as modifications for Grodins' model.

Gallagher, R. R.↗
NASA NTRS1993

Dynamics of the baryonic component in hierarchical clustering universes

I present self-consistent 3-D simulations of the formation of virialized systems containing both gas and dark matter in a flat universe. A fully Lagrangian code based on the Smoothed Particle Hydrodynamics technique and a tree data structure has been used to evolve regions of comoving radius 2-3 Mpc. Tidal effects are included by coarse-sampling the density of the outer regions up to a radius approx. 20 Mpc. Initial conditions are set at high redshift (z greater than 7) using a standard Cold Dark Matter perturbation spectrum and a baryon mass fraction of 10 percent (omega(sub b) = 0.1). Simulations in which the gas evolves either adiabatically or radiates energy at a rate determined locally by its cooling function were performed. This allows us to investigate with the same set of simulations the importance of radiative losses in the formation of galaxies and the equilibrium structure of virialized systems where cooling is very inefficient. In the absence of radiative losses, the simulations can be rescaled to the density and radius typical of galaxy clusters. A summary of the main results is presented.

Navarro, Julio↗
NASA NTRS1997

Kinetics and Thermochemistry of ClCO Formation from the Cl + CO Association Reaction

Laser flash photolysis of Cl2/CO/M mixtures (M = N2, CO, Ar, CO2) has been employed in conjunction with Cl((sup 2)P(sub J)) detection by time-resolved resonance fluorescence spectroscopy to investigate equilibration kinetics in the reactions Cl((sup 2)P(sub J)) + CO ClCO as a function of temperature (185-260 K) and pressure (14-200 Torr). The association and dissociation reactions are found to be in the low-pressure limit over the range of experimental conditions investigated. In N2 and/or CO buffer gases, the temperature dependences of the ClCO formation and dissociation reaction rate constants are described by the Arrhenius expressions k(sub 1) = (1.05 +/- 0.36) x 10(exp -34) exp[(810 +/- 70)/T] cm(exp 6)/molecules(exp 2).s and k(sub -1) = (4.1 +/- 3.1) x 10(exp -10) exp[(-2960 +/- 60)/T]cu cm/(molecule.s) (errors are 2 sigma). Second- and third-law analyses of the temperature dependence of the equilbrium constant (k/k-1) lead to the following thermodynamic parameters for the association reaction: Delta-H(sub 298) = -7.7 +/- 0.6 kcal/mol, Delta-H(sub 0) = -6.9 +/- 0.7 kcal/mol, Delta-S(sub 298) = -23.8 +/- 2.0 cal/mole.K, Delta-H(sub f,298)(ClCO) = 5.2 +/- 0.6 kcal/mol (errors are 2 sigma). The results repported in this study significantly reduce the uncertainties in all reported kinetic and thermodynamic parameters.

Nicovich, J. M.↗
NASA NTRS1998

Attenuated baroreflex control of sympathetic nerve activity after cardiovascular deconditioning in rats

The effect of cardiovascular deconditioning on baroreflex control of the sympathetic nervous system was evaluated after 14 days of hindlimb unloading (HU) or the control condition. Rats were chronically instrumented with catheters and sympathetic nerve recording electrodes for measurement of mean arterial pressure (MAP) and heart rate (HR) and recording of lumbar (LSNA) or renal (RSNA) sympathetic nerve activity. Experiments were conducted 24 h after surgery, with the animals in a normal posture. Baroreflex function was assessed using a logistic function that related HR and LSNA or RSNA to MAP during infusion of phenylephrine and nitroprusside. Baroreflex influence on HR was not affected by HU. Maximum baroreflex-elicited LSNA was significantly reduced in HU rats (204 +/- 11.9 vs. 342 +/- 30.6% baseline LSNA), as was maximum reflex gain (-4.0 +/- 0.6 vs. -7.8 +/- 1.3 %LSNA/mmHg). Maximum baroreflex-elicited RSNA (259 +/- 10.8 vs. 453 +/- 28.0% baseline RSNA), minimum baroreflex-elicited RSNA (-2 +/- 2.8 vs. 13 +/- 4.5% baseline RSNA), and maximum gain (-5.8 +/- 0.5 vs. -13.6 +/- 3.1 %RSNA/mmHg) were significantly decreased in HU rats. Results demonstrate that baroreflex modulation of sympathetic nervous system activity is attenuated after cardiovascular deconditioning in rodents. Data suggest that alterations in the arterial baroreflex may contribute to orthostatic intolerance after a period of bedrest or spaceflight in humans.

NASA Discipline Cardiopulmonary↗
NASA NTRS2014

A Linear Programming Approach to Routing Control in Networks of Constrained Nonlinear Positive Systems with Concave Flow Rates

We consider control design for positive compartmental systems in which each compartment's outflow rate is described by a concave function of the amount of material in the compartment.We address the problem of determining the routing of material between compartments to satisfy time-varying state constraints while ensuring that material reaches its intended destination over a finite time horizon. We give sufficient conditions for the existence of a time-varying state-dependent routing strategy which ensures that the closed-loop system satisfies basic network properties of positivity, conservation and interconnection while ensuring that capacity constraints are satisfied, when possible, or adjusted if a solution cannot be found. These conditions are formulated as a linear programming problem. Instances of this linear programming problem can be solved iteratively to generate a solution to the finite horizon routing problem. Results are given for the application of this control design method to an example problem. Key words: linear programming; control of networks; positive systems; controller constraints and structure.

Positive Systems↗
NASA NTRS1999

In-flight adaptive performance optimization (APO) control using redundant control effectors of an aircraft

Practical application of real-time (or near real-time) Adaptive Performance Optimization (APO) is provided for a transport aircraft in steady climb, cruise, turn descent or other flight conditions based on measurements and calculations of incremental drag from a forced response maneuver of one or more redundant control effectors defined as those in excess of the minimum set of control effectors required to maintain the steady flight condition in progress. The method comprises the steps of applying excitation in a raised-cosine form over an interval of from 100 to 500 sec. at the rate of 1 to 10 sets/sec of excitation, and data for analysis is gathered in sets of measurements made during the excitation to calculate lift and drag coefficients C.sub.L and C.sub.D from two equations, one for each coefficient. A third equation is an expansion of C.sub.D as a function of parasitic drag, induced drag, Mach and altitude drag effects, and control effector drag, and assumes a quadratic variation of drag with positions .delta..sub.i of redundant control effectors i=1 to n. The third equation is then solved for .delta..sub.iopt the optimal position of redundant control effector i, which is then used to set the control effector i for optimum performance during the remainder of said steady flight or until monitored flight conditions change by some predetermined amount as determined automatically or a predetermined minimum flight time has elapsed.

Gilyard, Glenn B.↗
DOE OSTI2025

Nannochloropsis oceanica IMET1 and its bacterial symbionts for carbon capture, utilization, and storage: biomass and calcium carbonate production under high pH and high alkalinity

ABSTRACT To combat the increasing levels of carbon dioxide (CO 2 ) released from the combustion of fossil fuels, microalgae have emerged as a promising strategy for biological carbon capture, utilization, and storage. This study used a marine microalgal strain, Nannochloropsis oceanica IMET1, which thrives in high CO 2 concentrations. A high-pH, high-alkalinity culture was designed for CO 2 capture through algal biomass production as well as permanent sequestration through calcium carbonate (CaCO 3 ) precipitation. This was accomplished by timed pH elevation and the addition of sodium bicarbonate to cultures of N. oceanica grown at lab scale (1 L) and pilot scale (500 L) with 10% and 5% CO 2 , respectively. Our data showed that 0.02 M NaHCO 3 promoted algal growth and that sparging cultures with ambient air after 12 days raised pH and created favorable CaCO 3 formation conditions. At the 1 L scale, we reached 1.52 g L −1 biomass after 12 days and an extra 9.3% CO 2 was captured in the form of CaCO 3 precipitates. At the 500 L pilot scale, an extra 60% CO 2 was captured (Day 40) with a maximum CO 2 capture rate of 63.2 g m −2 day −1 (Day 35). Bacterial communities associated with the microalgae were dominated by two novel Patescibacteria. Functional analysis revealed that genes for several plant growth-promotion traits (PGPTs) were enriched within this group. The microalgal-bacterial coculture system offers advantages for enhanced carbon mitigation through biomass production and simultaneous precipitation of recalcitrant CaCO 3 for long-term CO 2 storage. IMPORTANCE Capturing carbon dioxide (CO 2 ) released from fossil fuel combustion is of the utmost importance as the impacts of climate change continue to worsen. Microalgae can remove CO 2 through their natural photosynthetic pathways and are additionally able to convert CO 2 into a stable, recalcitrant form as calcium carbonate (CaCO 3 ). We demonstrate that microalgae-based carbon capture systems can be greatly improved with high pH and high alkalinity by providing optimal conditions for carbonate precipitation. Our results with the microalga, Nannochloropsis oceanica strain IMET1, show an extra 9.3% CO 2 captured as CaCO 3 at the 1 L scale and an extra 60% CO 2 captured at the 500 L (pilot) scale. Our optimized system provides a novel approach to capture CO 2 through two mechanisms: (i) as organic carbon within microalgal biomass and (ii) as inorganic carbon stored permanently in the form of CaCO 3.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Ectomycorrhizal fungi enhance pine growth by stimulating iron‐dependent mechanisms with trade‐offs in symbiotic performance

Iron (Fe) is crucial for metabolic functions of living organisms. Plants access occluded Fe through interactions with rhizosphere microorganisms and symbionts. Yet, the interplay between Fe addition and plant–mycorrhizal interactions, especially the molecular mechanisms underlying mycorrhiza-assisted Fe processing in plants, remains largely unexplored. Here, we conducted mesocosms in Pinus plants inoculated with different ectomycorrhizal fungi (EMF) Suillus species under conditions with and without Fe coatings. Meta-transcriptomic, biogeochemical, and X-ray fluorescence imaging analyses were applied to investigate early-stage mycorrhizal roots. While Fe addition promoted Pinus growth, it concurrently reduced mycorrhiza formation rate, symbiosis-related metabolites in plant roots, and aboveground plant carbon and macronutrient content. This suggested potential trade-offs between Fe-enhanced plant growth and symbiotic performance. However, the extent of this trade-off may depend on interactions between host plants and EMF species. Interestingly, dual EMF species were more effective at facilitating plant Fe uptake by inducing diverse Fe-related functions than single-EMF species. This subsequently triggered various Fe-dependent physiological and biochemical processes in Pinus roots, significantly contributing to Pinus growth. However, this resulted in a greater carbon allocation to roots, relatively reducing the aboveground plant carbon content. Our study offers critical insights into how EMF communities rebalance benefits of Fe-induced effects on symbiotic partners.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Mercury's Dynamic Magnetic Tail

The Mariner 10 and MESSENGER flybys of Mercury have revealed a magnetosphere that is likely the most responsive to upstream interplanetary conditions of any in the solar system. The source of the great dynamic variability observed during these brief passages is due to Mercury's proximity to the Sun and the inverse proportionality between reconnection rate and solar wind Alfven Mach number. However, this planet's lack of an ionosphere and its small physical dimensions also contribute to Mercury's very brief Dungey cycle, approx. 2 min, which governs the time scale for internal plasma circulation. Current observations and understanding of the structure and dynamics of Mercury's magnetotail are summarized and discussed. Special emphasis will be placed upon such questions as: 1) How much access does the solar wind have to this small magnetosphere as a function of upstream conditions? 2) What roles do heavy planetary ions play? 3) Do Earth-like substorms take place at Mercury? 4) How does Mercury's tail respond to extreme solar wind events such coronal mass ejections? Prospects for progress due to advances in the global magnetohydrodynamic and hybrid simulation modeling and the measurements to be taken by MESSENGER after it enters Mercury orbit on March 18, 2011 will be discussed.

Slavin, James A.↗
DOE OSTI2016

Ab Initio Reaction Kinetics of CH 3 O$\dot{C}$(=O) and $\dot{C}$H 2 OC(=O)H Radicals

The dissociation and isomerization kinetics of the methyl ester combustion intermediates methoxycarbonyl radical (CH3O$\dot{C}$(=O)) and (formyloxy)methyl radical ($\dot{C}$H2OC(=O)H) are investigated theoretically using high-level ab initio methods and Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory. Geometries obtained at the hybrid density functional theory (DFT) and coupled cluster singles and doubles with perturbative triples correction (CCSD(T)) levels of theory are found to be similar. We employ high-level ab initio wave function methods to refine the potential energy surface: CCSD(T), multireference singles and doubles configuration interaction (MRSDCI) with the Davidson–Silver (DS) correction, and multireference averaged coupled-pair functional (MRACPF2) theory. MRSDCI+DS and MRACPF2 capture the multiconfigurational character of transition states (TSs) and predict lower barrier heights than CCSD(T). The temperature- and pressure-dependent rate coefficients are computed using RRKM/ME theory in the temperature range 300–2500 K and a pressure range of 0.01 atm to the high-pressure limit, which are then fitted to modified Arrhenius expressions. Dissociation of CH3O$\dot{C}$(=O) to $\dot{C}$H3 and CO2 is predicted to be much faster than dissociating to CH3$\dot{O}$ and CO, consistent with its greater exothermicity. Isomerization between CH3O$\dot{C}$(=O) and $\dot{C}$H2OC(=O)H is predicted to be the slowest among the studied reactions and rarely happens even at high temperature and high pressure, suggesting the decomposition pathways of the two radicals are not strongly coupled. The predicted rate coefficients and branching fractions at finite pressures differ significantly from the corresponding high-pressure-limit results, especially at relatively high temperatures. Finally, because it is one of the most important CH3$\dot{O}$ removal mechanisms under atmospheric conditions, the reaction kinetics of CH3$\dot{O}$ + CO was also studied along the PES of CH3O$\dot{C}$(=O); the resulting kinetics predictions are in remarkable agreement with experiments.

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

Orbital-enriched flat-top partition of unity method for the Schrödinger eigenproblem

Quantum mechanical calculations require the repeated solution of a Schrödinger equation for the wavefunctions of the system, from which materials properties follow. Recent work has shown the effectiveness of enriched finite element type Galerkin methods at significantly reducing the degrees of freedom required to obtain accurate solutions. However, time to solution has been adversely affected by the need to solve a generalized rather than standard eigenvalue problem and the ill-conditioning of associated system matrices. In this work, we address both issues by proposing a stable and efficient orbital-enriched partition of unity method to solve the Schrödinger boundary-value problem in a parallelepiped unit cell subject to Bloch-periodic boundary conditions. In the proposed partition of unity method, the three-dimensional domain is covered by overlapping patches, with a compactly-supported weight function associated with each patch. A key ingredient in our approach is the use of non-negative weight functions that possess the flat-top property, i.e., each weight function is identically equal to unity over some finite subset of its support. This flat-top property provides a pathway to devise a stable approximation over the whole domain. On each patch, we use $p$th degree orthogonal (Legendre) polynomials that ensure $p$th order completeness, and in addition include eigenfunctions of the radial Schrödinger equation. Furthermore, we adopt a variational lumping approach to construct a (block-)diagonal overlap matrix that yields a standard eigenvalue problem for which there exist efficient eigensolvers. The accuracy, stability, and efficiency of the proposed method is demonstrated for the Schrödinger equation with a harmonic potential as well as a localized Gaussian potential. In conclusion, we show that the proposed approach delivers optimal rates of convergence in the energy, and the use of orbital enrichment significantly reduces the number of degrees of freedom for a given desired accuracy in the energy eigenvalues while the stability of the enriched approach is fully maintained.

42 ENGINEERING↗
NASA NTRS2001

Models of Small-Scale Patchiness

Patchiness is perhaps the most salient characteristic of plankton populations in the ocean. The scale of this heterogeneity spans many orders of magnitude in its spatial extent, ranging from planetary down to microscale. It has been argued that patchiness plays a fundamental role in the functioning of marine ecosystems, insofar as the mean conditions may not reflect the environment to which organisms are adapted. Understanding the nature of this patchiness is thus one of the major challenges of oceanographic ecology. The patchiness problem is fundamentally one of physical-biological-chemical interactions. This interconnection arises from three basic sources: (1) ocean currents continually redistribute dissolved and suspended constituents by advection; (2) space-time fluctuations in the flows themselves impact biological and chemical processes, and (3) organisms are capable of directed motion through the water. This tripartite linkage poses a difficult challenge to understanding oceanic ecosystems: differentiation between the three sources of variability requires accurate assessment of property distributions in space and time, in addition to detailed knowledge of organismal repertoires and the processes by which ambient conditions control the rates of biological and chemical reactions. Various methods of observing the ocean tend to lie parallel to the axes of the space/time domain in which these physical-biological-chemical interactions take place. Given that a purely observational approach to the patchiness problem is not tractable with finite resources, the coupling of models with observations offers an alternative which provides a context for synthesis of sparse data with articulations of fundamental principles assumed to govern functionality of the system. In a sense, models can be used to fill the gaps in the space/time domain, yielding a framework for exploring the controls on spatially and temporally intermittent processes. The following discussion highlights only a few of the multitude of models which have yielded insight into the dynamics of plankton patchiness. In addition, this particular collection of examples is intended to furnish some exposure to the diversity of modeling approaches which can be brought to bear on the problem. These approaches range from abstract theoretical models intended to elucidate specific processes, to complex numerical formulations which can be used to actually simulate observed distributions in detail.

McGillicuddy, D. J.↗
DOE OSTI2017

Formation of solid thorium monoxide at near-ambient conditions as observed by neutron reflectometry and interpreted by screened hybrid functional calculations

Oxidation of a ~1000 Å sputter-deposited thorium thin film at 150 °C in 100 ppm of flowing oxygen in argon produces the long-sought solid form of thorium monoxide. Changes in the scattering length density (SLD) distribution in the film over the 700-min experiment measured by in-situ, dynamic neutron reflectometry (NR) shows the densities, compositions and thickness of the various thorium oxides layers formed. Screened, hybrid density-functional theory calculations of potential thorium oxides aid interpretation, providing atomic-level picture and energetics for understanding oxygen migration. NR provided evidence of the formation of substoichiometric thorium oxide, ThO y (y < 1) at the interface between the unreacted thorium metal and its dioxide overcoat which grows inward, consuming the thorium at a rate of 2.1 Å/min while y increases until reaching 1:1 oxygen-to-thorium. Its presence indicates that kinetically-favored solid-phase ThO can be preferentially generated as a majority phase under the thermodynamically-favored ThO 2 top layer at conditions close to ambient.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Quantifying Shear-Induced Deformation and Detachment of Individual Adherent Sickle Red Blood Cells

According to the article, vaso-occlusive crisis, a common painful complication of sickle cell disease, is a complex process triggered by intercellular adhesive interactions among blood cells and the endothelium in all human organs (e.g., the oxygen-rich lung as well as hypoxic systems such as liver and kidneys). We present a combined experimental-computational study to quantify the adhesive characteristics of sickle mature erythrocytes (SMEs) and irreversibly sickled cells (ISCs) under flow conditions mimicking those in postcapillary venules. We employed an in vitro microfluidic cell adherence assay, which is coated uniformly with fibronectin. We investigated the adhesion dynamics of SMEs and ISCs in pulsatile flow under well-controlled hypoxic conditions, inferring the cell adhesion strength by increasing the flow rate (or wall shear stress (WSS)) until the onset of cell detachment. In parallel, we performed simulations of individual SMEs and ISCs under shear. We introduced two metrics to quantify the adhesion process, the cell aspect ratio (AR) as a function of WSS and its rate of change (the dynamic deformability index). We found that the AR of SMEs decreases significantly with the increase of WSS, consistent between the experiments and simulations. In contrast, the AR of ISCs remains constant in time and independent of the flow rate. The critical WSS value for detaching a single SME in oxygenated state is in the range of 3.9–5.5 Pa depending on the number of adhesion sites; the critical WSS value for ISCs is lower than that of SMEs. Our simulations show that the critical WSS value for SMEs in deoxygenated state is above 6.2 Pa (multiple adhesion sites), which is greater than their oxygenated counterparts. We investigated the effect of cell shear modulus on the detachment process; we found that for the same cell adhesion spring constant, the higher shear modulus leads to an earlier cell detachment from the functionalized surface. These findings may aid in the understanding of individual roles of sickle cell types in sickle cell disease vaso-occlusion.

Biophysics↗