Engineering topics
Publications and source records attributed to Radhakrishnan, Krishnan.
Our hypothesis predicts that retinal blood vessels increase in density during early-stage progression to moderate nonproliferative diabetic retinopathy (NPDR). The renin-angiotensin system (RAS) is implicated in the pathogenesis of DR and in the function of circulating angiogenic cells (CACs), a critical bone marrow-derived population that is instrumental in vascular repair.
Our hypothesis predicts that retinal blood vessels increase in density during early-stage progression to moderate nonproliferative diabetic retinopathy (NPDR). The prevailing paradigm of NPDR progression is that vessels drop out prior to abnormal, vision-impairing regrowth at late-stage proliferative diabetic retinopathy (DR). However, surprising results for our previous preliminary study 1 with NASA's VESsel GENeration Analysis (VESGEN) software showed that vessels proliferated considerably during moderate NPDR compared to drop out at both mild and severe NPDR. Validation of our hypothesis will support development of successful early-stage regenerative therapies such as vascular repair by circulating angiogenic cells (CACs). The renin-angiotensin system (RAS)is implicated in the pathogenesis of DR and in the function of CACs, a critical bone marrow-derived population that is instrumental in vascular repair.
Long-term exposure to space microgravity poses significant risks for visual impairment. Evidence suggests such vision changes are linked to cephalad fluid shifts, prompting a need to directly quantify microgravity-induced retinal vascular changes. The quality of retinal images used for such vascular remodeling analysis, however, is dependent on imaging methodology. For our exploratory study, we hypothesized that retinal images captured using fluorescein imaging methodologies would be of higher quality in comparison to images captured without fluorescein. A semi-automated image quality assessment was developed using Vessel Generation Analysis (VESGEN) software and MATLAB® image analysis toolboxes. An analysis of ten images found that the fluorescein imaging modality provided a 36% increase in overall image quality (two-tailed p=0.089) in comparison to nonfluorescein imaging techniques.
The effects of microvascular branching on blood flow were investigated in vivo by microscopic particle imaging velocimetry (micro-PIV). We use micro-PIV to measure blood flow by tracking red blood cells (RBC) as the moving particles. Velocity flow fields, including flow pulsatility, were analyzed for the first four branching orders of capillaries, postcapillary venules and small veins of the microvascular network within the developing avian yolksac at embryonic day 5 (E5). Increasing volumetric flowrates were obtained from parabolic laminar flow profiles as a function of increasing vessel diameter and branching order. Maximum flow velocities increased approximately twenty-fold as the function of increasing vessel diameter and branching order compared to flow velocities of 100 - 150 micron/sec in the capillaries. Results from our study will be useful for the increased understanding of blood flow within anastomotic, heterogeneous microvascular networks.
An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.
Vascular endothelial growth factor-165 (VEGF(sub 165)) stimulated angiogenesis in the quail chorioallantoic membrane (CAM) by vessel expansion from the capillary network. However, lymphangiogenesis was stimulated by the filopodial guidance of tip cells located on blind-ended lymphatic sprouts. As quantified by fractal/generational branching analysis using the computer code VESGEN, vascular density increased maximally at low VEGF concentrations, and vascular diameter increased most at high VEGF concentrations. Increased vascular density and diameter were statistically independent events (r(sub s), -0.06). By fluorescence immunohistochemistry of VEGF receptors VEGFR-1 and VEGFR-2, alpha smooth muscle actin ((alpha) SMA) and a vascular/lymphatic marker, VEGF(sub 165) increased the density and diameter of sprouting lymphatic vessels guided by tip cells (accompanied by the dissociation of lymphatics from blood vessels). Isolated migratory cells expressing (alpha)SMA were recruited to blood vessels, whereas isolated cells expressing VEGFR-2 were recruited primarily to lymphatics. In conclusion, VEGF(sub 165) increased lymphatic vessel density by lymphatic sprouting, but increased blood vessel density by vascular expansion from the capillary network.
An acute reduction in oxygen delivery to a tissue is associated with metabolic changes aimed at maintaining ATP homeostasis. However, given the complexity of the human bioenergetic system, it is difficult to determine quantitatively how cellular metabolic processes interact to maintain ATP homeostasis during stress (e.g., hypoxia, ischemia, and exercise). In particular, we are interested in determining mechanisms relating cellular oxygen concentration to observed metabolic responses at the cellular, tissue, organ, and whole body levels and in quantifying how changes in tissue oxygen availability affect the pathways of ATP synthesis and the metabolites that control these pathways. In this study, we extend a previously developed mathematical model of human bioenergetics, to provide a physicochemical framework that permits quantitative understanding of oxygen as a metabolic regulator. Specifically, the enhancement--sensitivity analysis--permits studying the effects of variations in tissue oxygenation and parameters controlling cellular respiration on glycolysis, lactate production, and pyruvate oxidation. The analysis can distinguish between parameters that must be determined accurately and those that require less precision, based on their effects on model predictions. This capability may prove to be important in optimizing experimental design, thus reducing use of animals.
Quantitative analysis of vascular generational branching demonstrated that transforming growth factor-beta1 (TGF-beta1), a multifunctional cytokine and angiogenic regulator, strongly inhibited angiogenesis in the arterial tree of the developing quail chorioallantoic membrane (CAM) by inhibition of the normal increase in the number of new, small vessels. The cytokine was applied uniformly in solution at embryonic day 7 (E7) to the CAMs of quail embryos cultured in petri dishes. After 24 h the rate of arterial growth was inhibited by as much as 105% as a function of increasing TGF-beta1 concentration. Inhibition of the rate of angiogenesis in the arterial tree by TGF-beta1 relative to controls was measured in digital images by three well-correlated, computerized methods. The first computerized method, direct measurement by the computer code VESGEN of vascular morphological parameters according to branching generations G(sub 1) through G(sub greater than or equal to 5), revealed that TGF-beta1 selectively inhibited the increase in the number density of small vessels, N(sub v greater than or equal to 5), (382 plus or minus 85 per square centimeter) for specimens treated with 1 microgram TGF-beta1/CAM for 24 h, compared to 583 plus or minus 99 per square centimeter for controls), but did not significantly affect other parameters such as average vessel length or vessel diameter. The second and third methods, the fractal dimension (D(sub f)) and grid intersection (rho (sub v)), are statistical descriptors of spatial pattern and density. According to D(sub f) and rho(sub v), arterial density increased in control specimens from 1.382 plus or minus 0.007 and 662 plus or minus 52 per square centimeters at E7 (0 h) to 1.439 plus or minus 0.013 and 884 plus or minus 55 per square centimeters at E8 (24 h), compared to 1.379 plus or minus 0.039 and 650 plus or minus 111 per square centimeter for specimens treated with 1 microgram TGF-beta1/CAM for 24 h. TGF-beta1 therefore regulates vascular pattern and the rate of angiogenesis in a unique fingerprint manner, as do other major angiogenic regulators that include VEGF, FGF-2 (bFGF), and angiostatin. TGF-beta1 did not stimulate angiogenesis significantly at low cytokine concentrations, which suggests that this quail CAM model of angiogenesis is not associated with an inflammatory response.
An acute reduction in oxygen delivery to skeletal muscle is generally associated with profound derangements in substrate metabolism. Given the complexity of the human bioenergetic system and its components, it is difficult to quantify the interaction of cellular metabolic processes to maintain ATP homeostasis during stress (e.g., hypoxia, ischemia, and exercise). Of special interest is the determination of mechanisms relating tissue oxygenation to observed metabolic responses at the tissue, organ, and whole body levels and the quantification of how changes in oxygen availability affect the pathways of ATP synthesis and their regulation. In this study, we apply a previously developed mathematical model of human bioenergetics to study effects of ischemia during periods of increased ATP turnover (e.g., exercise). By using systematic sensitivity analysis the oxidative phosphorylation rate was found to be the most important rate parameter affecting lactate production during ischemia under resting conditions. Here we examine whether mild exercise under ischemic conditions alters the relative importance of pathways and parameters previously obtained.
An acute reduction in oxygen delivery to a tissue is associated with metabolic changes aimed at maintaining ATP homeostasis. However, given the complexity of the human bio-energetic system, it is difficult to determine quantitatively how cellular metabolic processes interact to maintain ATP homeostasis during stress (e.g., hypoxia, ischemia, and exercise). In particular, we are interested in determining mechanisms relating cellular oxygen concentration to observed metabolic responses at the cellular, tissue, organ, and whole body levels and in quantifying how changes in tissue oxygen availability affect the pathways of ATP synthesis and the metabolites that control these pathways. In this study; we extend a previously developed mathematical model of human bioenergetics, to provide a physicochemical framework that permits quantitative understanding of oxygen as a metabolic regulator. Specifically, the enhancement - sensitivity analysis - permits studying the effects of variations in tissue oxygenation and parameters controlling cellular respiration on glycolysis, lactate production, and pyruvate oxidation. The analysis can distinguish between parameters that must be determined accurately and those that require less precision, based on their effects on model predictions. This capability may prove to be important in optimizing experimental design, thus reducing use of animals.
Although all tissues in the body can adapt to varying physiological/pathological conditions, muscle is the most adaptable. To understand the significance of cellular events and their role in controlling metabolic adaptations in complex physiological systems, it is necessary to link cellular and system levels by means of mechanistic computational models. The main objective of this work is to improve understanding of the regulation of energy metabolism during skeletal/cardiac muscle ischemia by combining in vivo experiments and quantitative models of metabolism. Our main focus is to investigate factors affecting lactate metabolism (e.g., NADH/NAD) and the inter-regulation between carbohydrate and fatty acid metabolism during a reduction in regional blood flow. A mechanistic mathematical model of energy metabolism has been developed to link cellular metabolic processes and their control mechanisms to tissue (skeletal muscle) and organ (heart) physiological responses. We applied this model to simulate the relationship between tissue oxygenation, redox state, and lactate metabolism in skeletal muscle. The model was validated using human data from published occlusion studies. Currently, we are investigating the difference in the responses to sudden vs. gradual onset ischemia in swine by combining in vivo experimental studies with computational models of myocardial energy metabolism during normal and ischemic conditions.
In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.
The goals of the trace chemistry group were to identify the processes relevant to aerosol and aerosol precursor formation occurring within aircraft gas turbine engines; that is, within the combustor, turbine, and nozzle. The topics of discussion focused on whether the chemistry of aerosol formation is homogeneous or heterogeneous; what species are important for aerosol and aerosol precursor formation; what modeling/theoretical activities to pursue; what experiments to carry out that both support modeling activities and elucidate fundamental processes; and the role of particulates in aerosol and aerosol precursor formation. The consensus of the group was that attention should be focused on SO2, SO3, and aerosols. Of immediate concern is the measurement of the concentration of the species SO3, SO2, H2SO4 OH, HO2, H2O2, O, NO, NO2, HONO, HNO3, CO, and CO2 and particulates in various engines, both those currently in use and those in development. The recommendation was that concentration measurements should be made at both the combustor exit and the engine exit. At each location the above species were classified into one of four categories of decreasing importance, Priority I through IV, as follows: Combustor exit: Priority I species - SO3:SO2 ratio, SO3, SO2, and particulates; Priority II species: OH and O; Priority III species - NO and NO2; and Priority IV species - CO and CO2. For the Engine exit: Priority I species - SO3:SO2 ratio, SO3, SO2,H2SO4, and particulates; Priority II species: OH,HO2, H2O2, and O; Priority III species - NO, NO2, HONO, and HNO3; and Priority IV species - CO and CO2. Table I summarizes the anticipated concentration range of each of these species. For particulate matter, the quantities of interest are the number density, size distribution, and composition. In order to provide data for validating multidimensional reacting flow models, it would be desirable to make 2-D, time-resolved measurements of the concentrations of the above species and, in addition, of the pressure, temperature, and velocity. A near term goal of the experimental program should be to confirm the nonlinear effects of sulfur speciation, and if present, to provide an explanation for them. It is also desirable to examine if the particulate matter retains any sulfur. The recommendation is to examine the effects on SOx production of variations in fuel-bound sulfur and aromatic content (which may affect the amount of particulates formed). These experiments should help us to understand if there is a coupling between particulate formation and SO, concentration. Similarly, any coupling with NOx can be examined either by introducing NOx into the combustion air or by using fuel-bound nitrogen. Also of immediate urgency is the need to establish and validate a detailed mechanism for sulfur oxidation/aerosol formation, whose chemistry is concluded to be homogeneous, because there is not enough surface area for heterogeneous effects. It is envisaged that this work will involve both experimental and theoretical programs. The experimental work will require, in addition to the measurements described above, fundamental studies in devices such as flow reactors and shock tubes. Complementing this effort should be modeling and theoretical activities. One impediment to the successful modeling of sulfur oxidation is the lack of reliable data for thermodynamic and transport properties for several species, such as aqueous nitric acid, sulfur oxides, and sulfuric acid. Quantum mechanical calculations are recommended as a convenient means of deriving values for these properties. Such calculations would also help establish rate constants for several important reactions for which experimental measurements are inherently fraught with uncertainty. Efforts to implement sufficiently detailed chemistry into computational fluid dynamic codes should be continued. Zero- and one-dimensional flow models are also useful vehicles for elucidating the minimal set of species and reactions that must be included in two- and three-dimensional modeling studies.
An acute reduction in oxygen (O2) delivery to a tissue is generally associated with a decrease in phosphocreatine, increases in ADP, NADH/NAD, and inorganic phosphate, increased rates of glycolysis and lactate production, and reduced rates of pyruvate and fatty acid oxidation. However, given the complexity of the human bioenergetic system and its components, it is difficult to determine quantitatively how cellular metabolic processes interact to maintain ATP homeostasis during stress (e.g., hypoxia, ischemia, and exercise). Of special interest is the determination of mechanisms relating tissue oxygenation to observed metabolic responses at the tissue, organ, and whole body levels and the quantification of how changes in tissue O2 availability affect the pathways of ATP synthesis and the metabolites that control these pathways. In this study, we extend a previously developed mathematical model of human bioenergetics to provide a physicochemical framework that permits quantitative understanding of O2 as a metabolic regulator. Specifically, the enhancement permits studying the effects of variations in tissue oxygenation and in parameters controlling the rate of cellular respiration on glycolysis, lactate production, and pyruvate oxidation. The whole body is described as a bioenergetic system consisting of metabolically distinct tissue/organ subsystems that exchange materials with the blood. In order to study the dynamic response of each subsystem to stimuli, we solve the ordinary differential equations describing the temporal evolution of metabolite levels, given the initial concentrations. The solver used in the present study is the packaged code LSODE, as implemented in the NASA Lewis kinetics and sensitivity analysis code, LSENS. A major advantage of LSENS is the efficient procedures supporting systematic sensitivity analysis, which provides the basic methods for studying parameter sensitivities (i.e., changes in model behavior due to parameter variation). Sensitivity analysis establishes relationships between model predictions and problem parameters (i.e., initial concentrations, rate coefficients, etc). It helps determine the effects of uncertainties or changes in these input parameters on the predictions, which ultimately are compared with experimental observations in order to validate the model. Sensitivity analysis can identify parameters that must be determined accurately because of their large effect on the model predictions and parameters that need not be known with great precision because they have little or no effect on the solution. This capability may prove to be important in optimizing the design of experiments, thereby reducing the use of animals. This approach can be applied to study the metabolic effects of reduced oxygen delivery to cardiac muscle due to local myocardial ischemia and the effects of acute hypoxia on brain metabolism. Other important applications of sensitivity analysis include identification of quantitatively relevant pathways and biochemical species within an overall mechanism, when examining the effects of a genetic anomaly or pathological state on energetic system components and whole system behavior.
Noise generated in gas turbine combustors can exist in several forms-broadband noise, sharp resonant peaks, and regular or intermittent nonlinear pulsing. In the present study, dynamic pressure measurements were made in several JP-5-fueled combustor configurations, at various mean pressures and temperatures. The fluctuating pressure was measured at mean pressures from 6 to 14 atm and inlet temperatures from 550 K to 850 K. The goal of the present work was to study the effect of changes in mean flow conditions on combustor noise: both broadband noise and sharp tones were considered. In general, the shape of the broadband noise spectrum was consistent from one configuration to another. The shape of the spectrum was influenced by the acoustic filtering of the combustion zone. This filtering ensured the basic consistency of the spectra. In general, the trends in broadband noise observed at low mean pressures were also seen at high mean pressures; that is, the total sound level decreased with both increasing equivalence ratio and increasing inlet temperature. The combustor configurations without a central pilot experienced higher broadband noise levels and were more susceptible to narrow peak resonances than configurations with a central pilot. The sharp peaks were more sensitive to the mean flow than was the broadband noise, and the effects were not always the same. In some situations, increasing the equivalence ratio made the sharp peaks grow, while at other conditions, increasing the equivalence ratio made the sharp peaks shrink. Thus, it was difficult to predict when resonances would occur; however, they were reproducible. Acoustic coupling between the upstream and downstream regions of the combustor may play a role in the sharp-peaked oscillations. Noise was also observed near lean blow out. As with other types of noise, lean blow out noise was affected by the combustion chamber acoustics, which apparently maintains the fluctuations at a uniform frequency. However, the actual conditions when this type of noise was experienced appeared to simply follow the lean blow out limit as it varied with mean temperature and pressure.
Combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system. The problem of combustion-acoustic instability is a concern in many devices for various reasons, as each device may have a unique mechanism causing unsteady heat release rates and many have unique boundary conditions. To accurately predict and quantify combustion-acoustic stabilities, the unsteady heat release rate and boundary conditions need to be accurately determined. The present review brings together work performed on a variety of practical combustion devices. Many theoretical and experimental investigations of the unsteady heat release rate have been performed, some based on perturbations in the fuel delivery system particularly for rocket instabilities, while others are based on hydrodynamic processes as in ramjet dump combustors. The boundary conditions for rocket engines have been analyzed and measured extensively. However, less work has been done to measure acoustic boundary conditions in many other combustion systems.
Lean, prevaporized, premixed combustors are susceptible to combustion-acoustic instabilities. A model was developed to predict eigenvalues of axial modes for combustion-acoustic interactions in a premixed combustor. This work extends previous work by including variable area and detailed chemical kinetics mechanisms, using the code LSENS. Thus the acoustic equations could be integrated through the flame zone. Linear perturbations were made of the continuity, momentum, energy, chemical species, and state equations. The qualitative accuracy of our approach was checked by examining its predictions for various unsteady heat release rate models. Perturbations in fuel flow rate are currently being added to the model.
Premixed combustors, which are being considered for low NOx engines, are susceptible to instabilities due to feedback between pressure perturbations and combustion. This feedback can cause damaging mechanical vibrations of the system as well as degrade the emissions characteristics and combustion efficiency. In a lean combustor instabilities can also lead to blowout. A model was developed to perform linear combustion-acoustic stability analysis using detailed chemical kinetic mechanisms. The Lewis Kinetics and Sensitivity Analysis Code, LSENS, was used to calculate the sensitivities of the heat release rate to perturbations in density and temperature. In the present work, an assumption was made that the mean flow velocity was small relative to the speed of sound. Results of this model showed the regions of growth of perturbations to be most sensitive to the reflectivity of the boundary when reflectivities were close to unity.