Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 1,603 records · Page 89

Experimental investigation of carbonate formation on Mars

If Mars had an early dense CO2 atmosphere (greater than 1 bar), then the transition to the present low pressure (7 mbar) may require carbonate formation without liquid water (e.g., although the early greenhouse has been questioned on the basis that modeled CO2 condensation prevents buildup of such a pressure). Our experiments have focused on the kinetics of such 'dry' reactions. Booth and Kieffer suggested rates consistent with removal of greater than 1 bar of CO2 over geologic time, but their results were questionable since their reactions did not proceed beyond one monolayer of carbonate formed on powdered silicate under martian conditions. Theoretical considerations implied that if nonporous rinds formed, then the rate of carbonate formation would be limited by diffusion, not reaction. Therefore, our experiments optimize reaction conditions to test for continued carbonate production, as well as to verify previous conclusions for initial growth of carbonate.

Stephens, Stuart K.↗

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗

Heterogeneous Combustion of Porous Graphite Particles in Microgravity

Recent theoretical investigations on graphite particle combustion have employed several levels of heterogeneous reaction models, ranging from global to elementary models, to describe the oxidation of carbon to gaseous products. Unlike the counterpart homogeneous reaction models, these heterogeneous reaction models are not well developed because of the difficulties associated with decoupling the physical characteristics of the solid (e.g. surface area taking part in combustion) from the chemical kinetic data. This is certainly true for porous graphite particle combustion, where heterogeneous and homogeneous reactions occur within the pores and play an important role in the overall oxidation process. As a result, there are considerable uncertainties of physical phenomena predicted using different heterogeneous kinetic models available in the literature. A good example, discussed later in this paper, is the predicted critical particle size below which the mass burning rate becomes exponentially small. The main goal of this study is to understand the basic mechanism controlling such rapid changes in burning rates, by developing a model where physical contributions are decoupled from chemical rate constants in a consistent manner. Another important goal of the proposed study is to develop a truly intrinsic, detailed heterogeneous reaction model for porous graphite combustion at high-temperatures, and to derive a systematically reduced heterogeneous reaction model in terms of the elementary reaction rate constants of the detailed model. The validation of chemical kinetic models describing the heterogeneous and homogeneous combustion in and around a spherically symmetric porous graphite particle can be considerably simplified by experimental measurements obtained under microgravity conditions. A vital component of this study is to conduct such supporting experiments on particle burning rate and surface temperature using NASA microgravity facilities, in close coordination with the theoretical effort. The basic understanding obtained and models developed as part of this project will be useful for optimal design of coal combustion devices. These models can also be extended to investigate the role of heterogeneous chemistry on pollutant formation pathways in combustion devices. The theoretical approach developed here, with pore diffusion effects decoupled from the chemical effects, can also be extended to understand the heterogeneous combustion of other porous fuels, for example, combustion of magnesium in a CO2 environment for propulsion in the Martian atmosphere.

Chelliah, Harsha K.↗

Structure of Laminar Permanently Blue, Opposed-Jet Ethylene-Fueled Diffusion Flames

The structure and state relationships of laminar soot-free (permanently blue) ethylene-fueled diffusion flames at various strain rates were studied both experimentally and computationally using an opposed-jet configuration. Measurements of gas velocities, temperatures, and compositions were carried out along the stagnation stream line. Corresponding predictions of flame structure were obtained, based on numerical simulations using several contemporary reaction mechanisms for methane oxidation. Flame conditions studied included ethylene-fueled opposed-jet diffusion flames having stoichiometric mixture fractions of 0.7 with measurements involving strain rates of 60-240/s and predictions involving strain rates of 0-1140/s at normal temperature and pressure. It was found that measured major gas species concentrations and temperature distributions were in reasonably good agreement with predictions using mechanisms due to GRI-Mech and Peters and that effects of preferential diffusion significantly influence flame structure even when reactant mass diffusivities are similar. Oxygen leakage to fuel-rich conditions and carbon monoxide leakage to fuel-lean conditions both increased as strain rates increased. Furthermore, increased strain rates caused increased fuel concentrations near the flame sheet, decreased peak gas temperatures, and decreased concentrations of carbon dioxide and water vapor throughout the flames. State relationships for major gas species and gas temperatures were found to exist over a broad range of strain rates, providing potential for significant computational simplifications for modeling purposes in some instances.

Lin, K.-C.↗

Structure of Laminar Permanently Blue, Opposed-Jet Ethylene-Fueled Diffusion Flames

The structure and state relationships of laminar soot-free (permanently blue) ethylene-fueled diffusion flames at various strain rates were studied both experimentally and computationally using an opposed-jet configuration. Measurements of gas velocities, temperatures, and compositions were carried out along the stagnation stream line. Corresponding predictions of flame structure were obtained, based on numerical simulations using several contemporary reaction mechanisms for methane oxidation. Flame conditions studied included ethylene-fueled opposed-jet diffusion flames having stoichiometric mixture fractions of 0.7 with measurements involving strain rates of 60-240/s and predictions involving strain rates of 0-1140/s at normal temperature and pressure. It was found that measured major gas species concentrations and temperature distributions were in reasonably good agreement with predictions using mechanisms due to GRI-Mech and Peters and that effects of preferential diffusion significantly influence flame structure even when reactant mass diffusivities are similar. Oxygen leakage to fuel-rich conditions and carbon monoxide leakage to fuel-lean conditions both increased as strain rates increased. Furthermore, increased strain rates caused increased fuel concentrations near the flame sheet, decreased peak gas temperatures, and decreased concentrations of carbon dioxide and water vapor throughout the flames. State relationships for major gas species and gas temperatures were found to exist over a broad range of strain rates, providing potential for significant computational simplifications for modeling purposes in some instances.

Lin, K.-C.↗

Cu Migration in Polycrystalline CdTe Solar Cells

An impurity reaction-diffusion model is applied to Cu defects and related intrinsic defects in polycrystalline CdTe for a better understanding of Cu’s role in the cell level reliability of CdTe PV devices. The simulation yields transient Cu distributions in polycrystalline CdTe during solar cell processing and stressing. Preliminary results for Cu migration using available diffusivity and solubility data show that Cu accumulates near the back contact, a phenomena that is commonly observed in devices after back-contact processing or stress conditions.

14 SOLAR ENERGY↗

Soot formation and precursor chemistry in Counterflow flames of aviation fuel surrogates

To meet market demands, the aviation sector is interested in utilizing drop-in Synthetic Aviation Turbine Fuels (SATF), either as neat fuels or in blends with conventional Jet A. SATF currently approved in standard specifications may have lower aromatic content with significant fractions of normal, branched, and cyclo-alkanes. Fundamental studies on soot formation from aviation fuels (Jet A, SATF) and their surrogate components are essential to understand how fuel composition influences soot and soot precursor formation. Here, this study reports new measurements of polycyclic aromatic hydrocarbons (PAH) and soot in counterflow diffusion flames (CDFs) of aviation fuel surrogates. Both intrusive and non-intrusive diagnostics are employed to determine the profiles of temperature, gas phase species, PAHs (up to C16), and soot volume fraction (SVF) in CDFs of iso-octane and surrogate mixtures. These measurements shed light on the transition of soot precursors to primary soot particles. In addition to serving as a common surrogate component in Jet A surrogate mixtures, iso-octane is a template species for larger, less volatile branched alkanes found in SATF mixtures. The newly developed Lawrence Livermore National Laboratory (LLNL) PAH and soot model successfully captures temperature, precursor species, and SVF profiles for the mixtures and conditions discussed in this work. Finally, a high-fidelity surrogate for Jet A is proposed that matches targeted physical and chemical properties well, while leveraging the wide range of candidate fuel molecules available in the LLNL detailed chemical model. The proposed surrogate formulation is validated against newly acquired measurements of the surrogate and literature measurements of Jet A. These new experiments and simulations provide critical insights into the PAH and soot formation from aviation fuels. Reaction pathways which require further investigation are highlighted, such that future work may bridge the remaining quantitative gaps in predicting soot formation from aviation fuel surrogates and surrogate components.

Aviation fuels↗

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗

Mercury contaminated sediment sites—An evaluation of remedial options

Mercury (Hg) is a naturally-occurring element that is ubiquitous in the aquatic environment. Though efforts have been made in recent years to decrease Hg emissions, historically-emitted Hg can be retained in the sediments of aquatic bodies where they may be slowly converted to methylmercury (MeHg). Consequently, Hg in historically-contaminated sediments can result in high levels of significant exposure for aquatic species, wildlife and human populations consuming fish. Even if source control of contaminated wastewater is achievable, it may take a very long time, perhaps decades, for Hg-contaminated aquatic systems to reach relatively safe Hg levels in both water and surface sediment naturally. It may take even longer if Hg is present at higher concentration levels in deep sediment. Hg contaminated sediment results from previous releases or ongoing contributions from sources that are difficult to identify. Due to human activities or physical, chemical, or biological processes (e.g. hydrodynamic flows, bioturbation, molecular diffusion, and chemical transformation), the buried Hg can be remobilized into the overlying water. Hg speciation in the water column and sediments critically affect the reactivity (i.e. conversion of inorganic Hg(II) to MeHg), transport, and its exposure to living organisms. Also, geochemical conditions affect the activity of methylating bacteria and its availability for methylation. This review paper discusses remedial considerations (e.g. key chemical factors in fate and transport of Hg, source characterization and control, environmental management procedures, remediation options, modeling tools) and includes practical case studies for cleaning up Hg-contaminated sediment sites. -- Highlights: ► Managing mercury-contaminated sediment sites are challenging to remediate. ► Remediation technologies are making a difference in managing these sites. ► Partitioning plays a dominant role in the distribution of mercury species. ► Mathematical models can be used to help us understand the chemistry and processes.

60 APPLIED LIFE SCIENCES↗

No double detonations but core carbon ignitions in high-resolution, grid-based simulations of binary white dwarf mergers

We study the violent phase of the merger of massive binary white dwarf systems. Our aim is to characterize the conditions for explosive burning to occur, and identify a possible explosion mechanism of Type Ia supernovae. The primary components of our model systems are carbon–oxygen (C/O) white dwarfs, while the secondaries are made either of C/O or of pure helium. We account for tidal effects in the initial conditions in a self-consistent way, and consider initially well-separated systems with slow inspiral rates. We study the merger evolution using an adaptive mesh refinement, reactive, Eulerian code in three dimensions, assuming symmetry across the orbital plane. We use a corotating reference frame to minimize the effects of numerical diffusion, and solve for self-gravity using a multigrid approach. We find a novel detonation mechanism in C/O mergers with massive primaries. Here, the detonation occurs in the primary's core and relies on the combined action of tidal heating, accretion heating, and self-heating due to nuclear burning. The exploding structure is compositionally stratified, with a reverse shock formed at the surface of the dense ejecta. The existence of such a shock has not been reported elsewhere. The explosion energy (1.6 × 10 51 erg) and 56 Ni mass (0.86 M ⊙ ) are consistent with an SN Ia at the bright end of the luminosity distribution, with an approximated decline rate of Δm 15 (B) ≈ 0.99. Our study does not support double-detonation scenarios in the case of a system with a 0.6 M ⊙ helium secondary and a 0.9 M ⊙ primary. Although the accreted helium detonates, it fails to ignite carbon at the base of the boundary layer or in the primary's core.

79 ASTRONOMY AND ASTROPHYSICS↗

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/$\mathrm{cm^3}$ to 16 g/$\mathrm{cm^3}$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. Here, we find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crystallization dynamics in glass-forming systems

Crystallization under far-from-equilibrium conditions is investigated for two different scenarios: crystallization of the metallic glass alloy Cu 50 Zr 50 and solidification of a transparent organic compound, o-terphenyl. For Cu 50 Zr 50 , crystallization kinetics are quanti ed through a new procedure that directly fits thermal analysis data to the commonly utilized JMAK model. The phase evolution during crystallization is quantified through in-situ measurements (HEXRD, DSC) and ex-situ microstructural analysis (TEM, HRTEM). The influence of chemical partitioning, diffusion, and crystallographic orientation on this sequence are examined. For o-terphenyl, the relationship between crystal growth velocity and interface undercooling is systematically studied via directional solidification.

36 MATERIALS SCIENCE↗

Study of diamond film growth and properties

The objective was to study diamond film growth and its properties in order to enhance the laser damage threshold of substrate materials. Calculations were performed to evaluate laser induced thermal stress parameter, R(sub T) of diamond. It is found that diamond has several orders of magnitude higher in value for R(sub T) compared to other materials. Thus, the laser induced damage threshold (LIDT) of diamond is much higher. Diamond films were grown using a microwave plasma enhanced chemical vapor deposition (MPECVD) system at various conditions of gas composition, pressure, temperature, and substrate materials. A 0.5 percent CH4 in H2 at 20 torr were ideal conditions for growing of high quality diamond films on substrates maintained at 900 C. The diamond films were polycrystalline which were characterized by scanning electron microscopy (SEM) and Raman scattering spectroscopy. The top surface of the growing film is always rough due to the facets of polycrystalline film while the back surface of the film replicates the substrate surface. An analytical model based on two dimensional periodic heat flow was developed to calculate the effective in-plane (face parallel) diffusivity of a two layer system. The effective diffusivity of diamond/silicon samples was measured using a laser pulse technique. The thermal conductivity of the films was measured to be 13.5 W/cm K, which is better than that of a type Ia natural diamond. Laser induced damage experiments were performed on bare Si substrates, diamond film coated Si, and diamond film windows. Significant improvements in the LIDT were obtained for diamond film coated Si compared to the bare Si.

Albin, Sacharial↗

Effect of Marangoni Convection Generated by Voids on Segregation During Low-G and 1-G Solidification

Solidification experiments, especially microgravity solidification experiments, are often compromised by the evolution of unwanted voids or bubbles in the melt. Although these voids and/or bubbles are highly undesirable, there is currently no effective means of preventing their formation or of eliminating their adverse effects, particularly during microgravity experiments. Marangoni convection caused by these voids can drastically change the transport processes in the melt. Recent microgravity experiments by Matthiesen (1) Andrews (2) and Fripp (3) are perfect examples of how voids and bubbles can affect the outcome of costly space experiments and significantly increase the level of difficulty in interpreting their results. Formation of bubbles have caused problems in microgravity experiments for a long time. Even in the early Skylab mission an unexpectedly large number of bubbles were detected in the four materials processing experiments reported by Papazian and Wilcox (4). They demonstrated that while during ground-based tests bubbles were seen to detach from the interface easily and float to the top of the melt, in low-gravity tests no detachment from the interface occurred and large voids were grown in the crystal. More recently, the lead-tin-telluride crystal growth experiment of Fripp et al.(3) flown aboard the USMP-3 mission has provided very interesting results. The purpose of the study was to investigate the effect of natural convection on the solidification process by growing the samples at different orientations with respect to the gravitational field. Large pores and voids were found in the three solid crystal samples processed in space. Post-growth characterization of the compositional profiles of the cells indicated considerable levels of mixing even in the sample grown in the hot-on-top stable configuration. The mixing was attributed to thermocapillary convection caused by the voids and bubbles which evolved during growth. Since the thermocapillary convection is orientation-independent, diffusion-controlled growth was not possible in any of the samples, even the top-heated one. These results are consistent with recent studies of thermocapillary convection generated by a bubble on a heated surface undertaken by Kassemi and Rashidnia (5-7) where it is numerically and experimentally shown that the thermocapillary flow generated by a bubble in a model fluid (silicone oil) can drastically modify the temperature field through vigorous mixing of the fluid around it, especially under microgravity conditions.

Kassemi, M.↗

Parfume (particle Fuel Model)

PARFUME is an advanced gas-cooled reactor fuel performance modeling and analysis code. It has been developed as an integrated mechanistic code that evaluates the thermal, mechanical, and physio-chemical behavior of fuel particles during irradiation and postulated accident conditions to determine the failure probability of a population of fuel particles given the particle-to-particle statistical variations in physical dimensions and material properties that arise from the fuel fabrication process, accounting for all viable mechanisms that can lead to particle failure. The code also determines the diffusion of fission products from the fuel through the particle coating layers, and through the fuel matrix to the coolant boundary. The subsequent release of fission products is calculated at the fuel element level (release of fission products from the compact or pebble). PARFUME calculates the release fraction as the ratio of the number of atoms released from the fuel element to the amount.

Gougar, Hans D.↗

Probing soft X-ray induced photoreduction of a model Mn-complex at cryogenic conditions

Soft X-ray absorption spectroscopy of first row transition elements at their respective L -edges provides important information about the oxidation and spin states of the metal centers. However, the associated sample damage in radiation-sensitive samples substantially alters the electronic and chemical structures of redox-active metal centers. Here, we measure the soft X-ray spectrum of the model Mn III (acac) 3 complex containing a redox-active Mn III metal center in an octahedral environment with a superconducting transition-edge sensor detector. To reduce the secondary damage resulting primarily from the diffusion of radicals and electrons, the spectra are collected at 30 K and 80 K on solid samples. Starting from the first scan, we detect the contribution of X-ray induced sample damage leading to a change in the Mn II intensity. However, at low temperatures, particularly at 30 K, we do not observe a gradual increase in the radiation damage with successive scans with the X-ray beam at the same spot. At our estimated dose of 90 kGy, we find 62% of Mn III (acac) 3 is still intact at 30 K. However, at room temperature, we see a gradual increase in radiation damage with increasing numbers of scans at the same spot, which is consistent with the possibility of increased diffusion rates of secondary radicals and electrons as noted in other studies.

L-edge spectroscopy↗

NiAl Oxidation Reaction Processes Studied In Situ Using MEMS-Based Closed-Cell Gas Reaction Transmission Electron Microscopy

The nanoscale oxidation mechanisms and kinetics of a model β-NiAl system were investigated using in situ closed-cell gas reaction scanning transmission electron microscopy (STEM). Here, we directly visualize the dynamic structural and chemical changes that occur during high-temperature oxidation at a high spatial resolution of 50.3Ni–49.7Al (at.%) nanoparticles under static air conditions at 730 Torr with heating up to 750 °C at 5 °C/s. A MEMS-based gas cell system, with microfabricated heater devices and a gas delivery system, was used to reveal site-specific oxidation initiation sites. Through time-resolved annular dark-field STEM imaging, we tracked the nanoscale oxidation kinetics of Al 2 O 3 . After oxidation at 750 °C, nucleation of voids at the Ni/Al 2 O 3 interface was observed along a NiAl grain boundary, followed by the formation of faceted NiO crystals. Small faceted cubic crystals of NiO were formed at the initial stage of oxidation at high PO 2 due to the outward self-diffusion of Ni 2+ ions, followed by the formation of a mixture of metastable and stable α-Al 2 O 3 at the oxide/metal interface that is attributed to a PO 2 decrease with oxidation time, which agreed with thermodynamic modeling calculations. Furthermore, the results from these in situ oxidation experiments in the β-NiAl system are in agreement with the established oxidation mechanisms; however, with in situ closed-cell gas microscopy it is now feasible to investigate nanoscale oxidation mechanisms and kinetics in real time and at high spatial resolution and can be broadly applied to understand the basic high-temperature oxidation mechanisms for a wide range of alloy compositions.

36 MATERIALS SCIENCE↗

Hepatocytes Determine the Hypoxic Microenvironment and Radiosensitivity of Colorectal Cancer Cells Through Production of Nitric Oxide That Targets Mitochondrial Respiration

Purpose: To determine whether host hepatocytes may reverse hypoxic radioresistance through nitric oxide (NO)-induced oxygen sparing, in a model relevant to colorectal cancer (CRC) liver metastases. Methods and Materials: Hepatocytes and a panel of CRC cells were incubated in a tissue-mimetic coculture system with diffusion-limited oxygenation, and oxygen levels were monitored by an oxygen-sensing fluorescence probe. To activate endogenous NO production, cocultures were exposed to a cytokine mixture, and the expression of inducible nitric oxide synthase was analyzed by reverse transcription–polymerase chain reaction, Western blotting, and NO/nitrite production. The mitochondrial targets of NO were examined by enzymatic activity. To assess hypoxic radioresponse, cocultures were irradiated and reseeded for colonies. Results: Resting hepatocytes consumed 10-40 times more oxygen than mouse CT26 and human DLD-1, HT29, HCT116, and SW480 CRC cells, and thus seemed to be the major effectors of hypoxic conditioning. As a result, hepatocytes caused uniform radioprotection of tumor cells at a 1:1 ratio. Conversely, NO-producing hepatocytes radiosensitized all CRC cell lines more than 1.5-fold, similar to the effect of selective mitochondrial inhibitors. The radiosensitizing effect was associated with a respiratory self-arrest of hepatocytes at the level of aconitase and complex II, which resulted in profound reoxygenation of tumor cells through oxygen sparing. Nitric oxide–producing hepatocytes were at least 10 times more active than NO-producing macrophages to reverse hypoxia-induced radioresistance. Conclusions: Hepatocytes were the major determinants of the hypoxic microenvironment and radioresponse of CRC cells in our model of metabolic hypoxia. We provide evidence that reoxygenation and radiosensitization of hypoxic CRC cells can be achieved through oxygen sparing induced by endogenous NO production in host hepatocytes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗