Results for “Diffusion equation” — Page 68

NASA NTRS1973

Unsteady, combined radiation and conduction in an absorbing, scattering, and emitting medium.

The transient cooldown of a gray, absorbing, isotropic scattering, emitting, and conducting medium bounded by gray, diffusely emitting and reflecting parallel plates is considered. Numerical solutions are obtained for the initial boundary-value problem with a discontinuous decrease in temperature at one boundary. The quasi-steady equation of radiative transfer is solved using Gaussian quadrature and a matrix eigenvector technique together with explicit numerical solution of the unsteady energy equation. Temperature and energy flux distributions are presented for variations of optical thickness, boundary emissivity, albedo, and conduction-radiation parameter.

Weston, K. C.↗

NASA NTRS1982

Small particle transport across turbulent nonisothermal boundary layers

The interaction between turbulent diffusion, Brownian diffusion, and particle thermophoresis in the limit of vanishing particle inertial effects is quantitatively modeled for applications in gas turbines. The model is initiated with consideration of the particle phase mass conservation equation for a two-dimensional boundary layer, including the thermophoretic flux term directed toward the cold wall. A formalism of a turbulent flow near a flat plate in a heat transfer problem is adopted, and variable property effects are neglected. Attention is given to the limit of very large Schmidt numbers and the particle concentration depletion outside of the Brownian sublayer. It is concluded that, in the parameter range of interest, thermophoresis augments the high Schmidt number mass-transfer coefficient by a factor equal to the product of the outer sink and the thermophoretic suction.

Rosner, D. E.↗

NASA NTRS2000

Design and Performance of McRas in SCMs and GEOS I/II GCMs

The design of a prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) for general circulation models (GCMs) will be discussed. McRAS distinguishes three types of clouds: (1) convective, (2) stratiform, and (3) boundary-layer types. The convective clouds transform and merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts into precipitation following the auto-conversion equations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, as well as diffuse both horizontally and vertically with a fully interactive cloud-microphysics throughout the life-cycle of the cloud, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III and 5-ARN CART datasets has shown that together with the rest of the model physics, McRAS can simulate the observed temperature, humidity, and precipitation without many systematic errors. The time history and time mean incloud water and ice distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the convective anvil and towers, and the convective updraft and downdraft velocities and mass fluxes all show a realistic behavior. Performance of McRAS in GEOS 11 GCM shows several satisfactory features but some of the remaining deficiencies suggest need for additional research involving convective triggers and inhibitors, provision for continuously detraining updraft, a realistic scheme for cumulus gravity wave drag, and refinements to physical conditions for ascertaining cloud detrainment level.

Sud, Yogesh C.↗

NASA NTRS2001

Model study of ATP and ADP buffering, transport of Ca(2+) and Mg(2+), and regulation of ion pumps in ventricular myocyte

We extended the model of the ventricular myocyte by Winslow et al. (Circ. Res 1999, 84:571-586) by incorporating equations for Ca(2+) and Mg(2+) buffering and transport by ATP and ADP and equations for MgATP regulation of ion transporters (Na(+)-K(+) pump, sarcolemmal and sarcoplasmic Ca(2+) pumps). The results indicate that, under normal conditions, Ca(2+) binding by low-affinity ATP and diffusion of CaATP may affect the amplitude and time course of intracellular Ca(2+) signals. The model also suggests that a fall in ATP/ADP ratio significantly reduces sarcoplasmic Ca(2+) content, increases diastolic Ca(2+), lowers systolic Ca(2+), increases Ca(2+) influx through L-type channels, and decreases the efficiency of the Na(+)/Ca(2+) exchanger in extruding Ca(2+) during periodic voltage-clamp stimulation. The analysis suggests that the most important reason for these changes during metabolic inhibition is the down-regulation of the sarcoplasmic Ca(2+)-ATPase pump by reduced diastolic MgATP levels. High Ca(2+) concentrations developed near the membrane might have a greater influence on Mg(2+), ATP, and ADP concentrations than that of the lower Ca(2+) concentrations in the bulk myoplasm. The model predictions are in general agreement with experimental observations measured under normal and pathological conditions.

NASA Discipline Cardiopulmonary↗

NASA NTRS2024

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗

NASA NTRS1982

Environmental solid particle effects on compressor cascade performance

The effect of suspended solid particles on the performance of the compressor cascade was investigated experimentally in a specially built cascade tunnel, using quartz sand particles. The cascades were made of NACA 65(10)10 airfoils. Three cascades were tested, one accelerating cascade and two diffusing cascades. The theoretical analysis assumes inviscid and incompressible two dimensional flow. The momentum exchange between the fluid and the particle is accounted for by the interphase force terms in the fluid momentum equation. The modified fluid phase momentum equations and the continuity equation are reduced to the conventional stream function vorticity formulation. The method treats the fluid phase in the Eulerian system and the particle phase in Lagrangian system. The experimental results indicate a small increase in the blade surface static pressures, while the theoretical results indicate a small decrease. The theoretical analysis, also predicts the loss in total pressure associated with the particulate flow through the cascade.

Tabakoff, W.↗

NASA NTRS1982

Flux-vector splitting for the Euler equations

Flux-vector splitting is shown to be the simplest way of implementing upwind differencing for the full or isenthalpic Euler equations combined with the ideal-gas law. The scheme presented produces steady shock profiles with two interior zones. A disadvantage in using flux-vector splitting is that it leads to numerical diffusion of a contact discontinuity at rest. This diffusion, however, can be removed, and present research is aimed at achieving this with minimal computational effort.

Van Leer, B.↗

NASA NTRS1985

Opening of the magnetic field lines in a fast rotating magnetosphere, with an application to Jupiter

A simple model of a magnetosphere rotating around a rotating Jupiter-like object with a spin-aligned dipolar moment mu is considered. A low-energy plasma is released by inner sources located beyond the corotation radius and is diffused outward through closed lines. The lines form a closed thin equatorial disk, maintaining a quasi-static balance between the centrifugal force and the magnetic tension. Solutions are obtained for the analytical equations characterizing the plasma in the closed region; the structure of the magnetic field outside the disk; and the self-consistent value of the critical radius r sub zero. The analytical solutions are compared to the observed properties of the Jovian magnetosphere, and the results are discussed in detail.

Aly, J. J.↗

NASA NTRS1990

Symmetrization of conservation laws with entropy for high-temperature hypersonic computations

Results of Hughes, France, and Mallet are generalized to conservation law systems taking into account high-temperature effects. Symmetric forms of different equation sets are derived in terms of entropy variables. First, the case of a general divariant gas is studied; it can be specialized to the usual Navier-Stokes equations, as well as to situations where the gas is vibrationally excited, and undergoes equilibrium chemical reactions. The case of gas in thermochemical nonequilibrium is considered next. Transport phenomena, and in particular mass diffusion, are examined in the framework of symmetric advective-diffusive systems.

Chalot, F.↗

NASA NTRS1992

Theoretical study of the design of a catalyst for para to ortho hydrogen conversion

The theory of Petzinger and Scalapino (1973) was thoroughly reviewed, and all of the basic equations for paramagnetic para to ortho hydrogen catalysis re-derived. There are only a few minor phase errors and errors of omission in the description of the theory. Three models (described by Petzinger and Scalapino) for the rate of para to ortho H2 catalysis were worked out, and uniform agreement obtained to within a constant factor of 2 pi. The analytical methods developed in the course of this study were then extended to two new models, which more adequately describe the process of surface catalysis including transfer of hydrogen molecules onto and off of the surface. All five equations for the para to ortho catalytic rate of conversion are described. The two new equations describe the catalytic rate for these models: H2 on the surface is a 2-D gas with lifetime tau; and H2 on the surface is a 2-D liquid undergoing Brownian motion (diffusion) with surface lifetime tau.

Coffman, Robert E.↗

NASA NTRS1997

The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems

In this paper, we study the Local Discontinuous Galerkin methods for nonlinear, time-dependent convection-diffusion systems. These methods are an extension of the Runge-Kutta Discontinuous Galerkin methods for purely hyperbolic systems to convection-diffusion systems and share with those methods their high parallelizability, their high-order formal accuracy, and their easy handling of complicated geometries, for convection dominated problems. It is proven that for scalar equations, the Local Discontinuous Galerkin methods are L(sup 2)-stable in the nonlinear case. Moreover, in the linear case, it is shown that if polynomials of degree k are used, the methods are k-th order accurate for general triangulations; although this order of convergence is suboptimal, it is sharp for the LDG methods. Preliminary numerical examples displaying the performance of the method are shown.

Cockburn, Bernardo↗

NASA NTRS2000

Determination of Interfacial Rheological Properties through Microgravity Oscillations of Bubbles and Drops

This report summarizes our derivations of analytical expressions for the frequencies and damping constants for small-amplitude axisymmetric shape oscillations of a liquid drop suspended in an immiscible fluid host in microgravity. In particular, this work addresses large Reynolds number shape oscillations and focuses on the surface rheological effects that arise from the presence of insoluble surfactants at the interface. Parameters characterizing viscous effects from the bulk phases, surface viscous effects, Marangoni effects from the surface advection and diffusion of surfactants, and the Gibbs elasticity are all considered and analyzed to determine the relative importance of each contribution. Supplementing the analytical treatment for small-amplitude oscillations, a numerical boundary integral equation formulation is developed for the study of large-amplittide axisymmetric oscillations of a drop in vacuum. The boundary integral formulation is an extension of classical potential flow theory and approximately accounts for viscous effects in the bulk fluid as well as the surface viscous and Marangoni effects resulting from an insoluble surfactant contaminating the interface. Theoretical and numerical results are presented for four distinct cases. These, range from the case when the effects of the surfactants are 'negligible' to 'large' when compared to the viscous effects in the bulk phases. The feasibility of the non-contact measurement of the surface parameters, using experimental observations for the oscillation frequencies and damping constants of drops and bubbles, is discussed.

Nadim, Ali↗

NASA NTRS2014

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Gorazd Medic↗

NASA NTRS2014

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Gorazd Medic↗

NASA NTRS2017

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Centrifugal compressor↗

NASA NTRS1991

Finite element computation of compressible flows with the SUPG formulation

Finite element computation of compressible Euler equations is presented in the context of the streamline-upwind/Petrov-Galerkin (SUPG) formulation. The SUPG formulation, which is based on adding stabilizing terms to the Galerkin formulation, is further supplemented with a shock capturing operator which addresses the difficulty in maintaining a satisfactory solution near discontinuities in the solution field. The shock capturing operator, which has been derived from work done in entropy variables for a similar operator, is shown to lead to an appropriate level of additional stabilization near shocks, without resulting in excessive numerical diffusion. An implicit treatment of the impermeable wall boundary condition is also presented. This treatment of the no-penetration condition offers increased stability for large Courant numbers, and accelerated convergence of the computations for both implicit and explicit applications. Several examples are presented to demonstrate the ability of this method to solve the equations governing compressible fluid flow.

Le Beau, G. J.↗

NASA NTRS1974

The spin-temperature theory of dynamic nuclear polarization and nuclear spin-lattice relaxation

A detailed derivation of the equations governing dynamic nuclear polarization (DNP) and nuclear spin lattice relaxation by use of the spin temperature theory has been carried to second order in a perturbation expansion of the density matrix. Nuclear spin diffusion in the rapid diffusion limit and the effects of the coupling of the electron dipole-dipole reservoir (EDDR) with the nuclear spins are incorporated. The complete expression for the dynamic nuclear polarization has been derived and then examined in detail for the limit of well resolved solid effect transitions. Exactly at the solid effect transition peaks, the conventional solid-effect DNP results are obtained, but with EDDR effects on the nuclear relaxation and DNP leakage factor included. Explicit EDDR contributions to DNP are discussed, and a new DNP effect is predicted.

Byvik, C. E.↗

NASA NTRS1986

Diffusion flame extinction in slow convenctive flow under microgravity environment

A theoretical analysis is presented to study the extinction characteristics of a diffusion flame near the leading edge of a thin fuel plate in slow, forced convective flows in a microgravity environment. The mathematical model includes two-dimensional Navier-Stokes momentum, energy and species equations with one-step overall chemical reaction using second-order finite rate Arrhenius kinetics. Radiant heat loss on the fuel plate is applied in the model as it is the dominant mechanism for flame extinguishment in the small convective flow regime. A parametric study based on the variation of convective flow velocity, which varies the Damkohler number (Da), and the surface radiant heat loss parameter (S) simultaneously, is given. An extinction limit is found in the regime of slow convective flow when the rate of radiant heat loss from fuel surface outweighs the rate of heat generation due to combustion. The transition from existent envelope flame to extinguishment consists of gradual flame contraction in the opposed flow direction together with flame temperature reduction as the convective flow velocity decreases continuously until the extinction limit is reached. A case of flame structure subjected to surface radiant heat loss is also presented and discussed.

Chen, C. H.↗