Search NASA⌕ Search

SEARCH · Search NASA

Results for “solute transport modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Iterative solution of the multistream electron transport equation. I - Comparison with laboratory beam injection experiments

The Neumann iteration method presently used for solving the electron transport equation in which energy, attitude, and pitch angle are independent variables is fast, and can compute numerical point-response-function solutions of the electron transport equation. Because both the inelastic cross sections and angular elastic cross sections of the model are empirically based, the solutions obtained represent a test of compatibility between various sets of cross sections and energy deposition measurements. The use of a numerical quadrature based on analytic phase function forms yields accurate phase function integrals at low computational cost.

Porter, H. S.↗
NASA NTRS1989

Diffusional aspects of the high-temperature oxidation of protective coatings

The role of diffusional transport associated with the high-temperature oxidation of coatings is examined, with special attention given to the low-pressure plasma spraying MCrAl-type overlay coatings and similar Ni-base alloys which form protective AlO3 scales. The use of diffusional analysis to predict the minimum solute concentration necessary to form and grow a solute oxide scale is illustrated. Modeling procedures designed to simulate the diffusional transport in coatings and substrates are presented to show their use in understanding coating degradation, predicting the protective life of a coating, and evaluating various coating parameters to guide coating development.

Nesbitt, J. A.↗
NASA NTRS2023

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

Solar energetic particles↗
NASA NTRS2019

Effects of the Serber First Step in 3DHZETRN-v2.1

3DHZETRN-v2 includes a detailed three dimensional (3D) treatment of neutron/light-ion transport based on a quasi-elastic/multiple-production assumption allowing improved agreement of the neutron/light-ion fluence compared with results of three Monte Carlo (MC) codes in the sense that the variance with respect to the individual MC results is less than the variance among the MC code results. The current numerical methods are no longer the main limitation to HZETRN code development and further changes in the nuclear model are required. In a prior study, an improved quasi-elastic spectrum based on a solution of the transport approximation to nuclear media effects showed promise, but the remaining multiple-production spectrum was based on a database derived from the Ranft model that used Bertini multiplicities. In the present paper, we will implement a more complete Serber first step into the 3DHZETRN-v2 code, but we retain the Bertini-Ranft branching ratios and evaporation multiplicities. It is shown that the new Serber model in the 3HZETRN-v2 code reduces the variance with individual MC codes, which are largely due to nuclear cross section model differences.

Wilson, John W.↗
NASA NTRS1999

Interaction of Porosity with an Advancing Solid/Liquid Interface: a Real-Time Investigation

Problems associated with formation of porosity during solidification continue to have a daily impact on the metal forming industry. Several past investigations have dealt with the nucleation and growth aspects of porosity. However, investigations related to the interaction of porosity with that of a solidification front has been limited mostly to organic analogues. In this paper we report on real time experimental observations of such interactions in metal alloys. Using a state of the art X-Ray Transmission Microscope (XTM) we have been able to observe and record the dynamics of the interaction. This includes distortion of the solid/liquid interface near a poro.sity, solute segr,egation patterns surrounding a porosity and the change in shape of the porosity during interaction with an advancing solid/liquid interface. Results will be presented for different Al alloys and growth conditions. The experimental data will be compared to theory using a recently developed 2D numerical model. The model employs a finite difference approach where the solid/liquid interface is defined through the points at which the interface intersects the grid lines. The transport variables are calculated at these points and the motion of the solidification front is determined by the magnitude of the transport variables. The model accounts for the interplay of the thermal and solutal field and the influence of capilarity to predict the shape of the solid/liquid interface with time in the vicinity of porosity. One can further calculate the perturbation of the solutal field by the presence of porosity in the melt.

Sen, S.↗
NASA NTRS1971

Steady-state temperature distribution in living tissue modeled as cylindrical shells.

Closed form, analytical solutions to the problem of steady-state heat transfer in living tissue modeled as cylindrical shells are presented and discussed. These solutions are particularly useful for the study of temperature distributions in the extremities. Metabolic heat generation, conduction, and heat transported by the blood perfusing the tissue are considered in the model. The results demonstrate the important role that the blood stream plays in the transfer of heat inside living tissue. Solutions are also presented for the limiting cases of diminishing blood flow that would occur during vasoconstriction or occlusion of blood by external means.

Shitzer, A.↗
NASA NTRS1981

An evaluation of computer codes for simulating the Galileo Probe aerothermal entry environment

The approaches of three computer flow field codes (HYVIS, COLTS, and RASLE), used to determine the Galileo Probe aerothermal environment and its effect on the design of the thermal protection system, are analyzed in order to resolve differences in their predicted results. All three codes account for the hypersonic, massively blown, radiation shock layers, characteristic of Jupiter entry. Significant differences, however, are evident in their solution procedures: the governing conservation equations, the numerical differencing methods, the governing physics (chemical, radiation, diffusion, and turbulence models), and the basic physical data (thermodynamic, transport, chemical, and spectral properties for atomic and molecular species). Solutions are compared for two near peak heating entry conditions for a Galileo Probe baseline configuration, having an initial mass of 242 kg and simulating entry into the Orton nominal atmosphere. The modern numerical methodology of COLTS and RASLE appear to provide an improved capability for coupled flow-field solutions.

Menees, G. P.↗
NASA NTRS2021

Assessment of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon transition correlations that were derived from a limited number of low-speed experiments and these models often fail to account for all of the relevant transition mechanisms and/or the variation in those mechanisms with respect to changes in the significant flow parameters. Available data from a recent experiment on the Common Research Model with Natural Laminar Flow (CRM-NLF) in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.3b solver for a swept wing configuration with nonzero taper and transonic cruise conditions. Specifically, the OVERFLOW solutions are used to evaluate the accuracy and robustness of the transport-equation-based transition models. Results highlight that the Spalart-Allmaras-based amplification factor transport (AFT-2017b) equation model and Menter’s shear-stress transport equation (SST2003)-based Langtry-Menter transition models (either with or without the modeling of crossflow transition) significantly underpredict the reported extent of laminar flow region over the entire span of the wing, irrespective of which instability mechanism(s) is expected to dominate the onset of the transition process. We show that the transition correlations underlying these models fail to account for the stabilizing effect of compressibility on the Tollmien-Schlichting transition, which is likely to be a major contributor to the underprediction of the laminar flow region on the CRM-NLF. The SST-2003-based Langtry-Menter model also appears to inaccurately predict the chordwise pressure variation along the majority of the wing span at all the flow conditions studied herein, due to how the turbulence intensity levels were enforced in the computations and how that was interfering with the functioning of the underlying turbulence model within the boundary layer. The AFT and Langtry-Menter models appear to be sensitive to the level of the freestream turbulence intensity, but the degree of sensitivity varies across the models.

CFD modeling↗
NASA NTRS1973

Study of helium emissions from active solar regions

A theoretical study is made of the visible and UV line radiation of He I atoms and He II ions from a plane-parallel model flare layer. Codes were developed for the solution of the statistically steady state equation for a 30 level He I - II - III model, and the line and continuum transport equations. These codes are described and documented in the report along with sample solutions. Optical depths and some line intensities are presented for a 1000 km thick layer. Solutions of the steady state equations are presented for electron temperatures 10,000 to 50,000 K and electron densities 10 to the 10th power to 10 to the 14th power/cu cm.

Kulander, J. L.↗
NASA NTRS2004

Minnowbrook IV: 2003 Workshop on Transition and Unsteady Aspects of Turbomachinery Flows

This Minnowbrook IV 2003 workshop on Transition and Unsteady Aspects of Turbomachinery Flows includes the following topics: 1) Current Issues in Unsteady Turbomachinery Flows; 2) Global Instability and Control of Low-Pressure Turbine Flows; 3) Influence of End Wall Leakage on Secondary Flow Development in Axial Turbines; 4) Active and Passive Flow Control on Low Pressure Turbine Airfoils; 5) Experimental and Numerical Investigation of Transitional Flows as Affected by Passing Wakes; 6) Effects of Freestream Turbulence on Turbine Blade Heat Transfer; 7) Bypass Transition Via Continuous Modes and Unsteady Effects on Film Cooling; 8) High Frequency Surface Heat Flux Imaging of Bypass Transition; 9) Skin Friction and Heat Flux Oscillations in Upstream Moving Wave Packets; 10) Transition Mechanisms and Use of Surface Roughness to Enhance the Benefits of Wake Passing in LP Turbines; 11) Transient Growth Approach to Roughness-Induced Transition; 12) Roughness- and Freestream-Turbulence-Induced Transient Growth as a Bypass Transition Mechanism; 13) Receptivity Calculations as a Means to Predicting Transition; 14) On Streamwise Vortices in a Curved Wall Jet and Their Effect on the Mean Flow; 15) Plasma Actuators for Separation Control of Low Pressure Turbine Blades; 16) Boundary-Layer Separation Control Under Low-Pressure-Turbine Conditions Using Glow-Discharge Plasma Actuators; 17) Control of Separation for Low Pressure Turbine Blades: Numerical Simulation; 18) Effects of Elevated Free-Stream Turbulence on Active Control of a Separation Bubble; 19) Wakes, Calming and Transition Under Strong Adverse Pressure Gradients; 20) Transitional Bubble in Periodic Flow Phase Shift; 21) Modelling Spots: The Calmed Region, Pressure Gradient Effects and Background; 22) Modeling of Unsteady Transitional Flow on Axial Compressor Blades; 23) Challenges in Predicting Component Efficiencies in Turbomachines With Low Reynolds Number Blading; 24) Observations on the Causal Relationship Between Blade Count and Developing Rotating Stall in a Four Stage Axial Compressor; 25) Experimental and Numerical Study of Non-Linear Interactions in Transonic Nozzle Flow; 26) Clocking Effects on a Modern Stage and One-Half Transonic Turbine; 27) DNS and LES of Transition on Turbine Blades; 28) The Use of Cellular Automata in Modeling the Transition; 29) Predicting Unsteady Buffet Onset Using RANS Solutions; 30) Transition Modelling With the SST Turbulence Model and an Intermittency Transport; and 31) Equation Workshop Summary Transcript

LaGraff, John E.↗
NASA NTRS2002

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗
NASA NTRS1973

Preliminary numerical analysis of improved gas chromatograph model

A mathematical model for the gas chromatograph was developed which incorporates the heretofore neglected transport mechanisms of intraparticle diffusion and rates of adsorption. Because a closed-form analytical solution to the model does not appear realizable, techniques for the numerical solution of the model equations are being investigated. Criteria were developed for using a finite terminal boundary condition in place of an infinite boundary condition used in analytical solution techniques. The class of weighted residual methods known as orthogonal collocation is presently being investigated and appears promising.

Woodrow, P. T.↗
NASA NTRS1976

Modeling of three-dimensional mixing and reacting ducted flows

A computer code, based upon a finite element solution algorithm, was developed to solve the governing equations for three-dimensional, reacting boundary region, and constant area ducted flow fields. Effective diffusion coefficients are employed to allow analyses of turbulent, transitional or laminar flows. The code was used to investigate mixing and reacting hydrogen jets injected from multiple orifices, transverse and parallel to a supersonic air stream. Computational results provide a three-dimensional description of velocity, temperature, and species-concentration fields downstream of injection. Experimental data for eight cases covering different injection conditions and geometries were modeled using mixing length theory (MLT). These results were used as a baseline for examining the relative merits of other mixing models. Calculations were made using a two-equation turbulence model (k+d) and comparisons were made between experiment and mixing length theory predictions. The k+d model shows only a slight improvement in predictive capability over MLT. Results of an examination of the effect of tensorial transport coefficients on mass and momentum field distribution are also presented. Solutions demonstrating the ability of the code to model ducted flows and parallel strut injection are presented and discussed.

Zelazny, S. W.↗
NASA NTRS1978

Effects of particle drift on cosmic ray transport. II - Analytical solution to the modulation problem with no latitudinal diffusion

An analytical solution to a model of the modulation of galactic cosmic rays in the presence of particle drifts is presented and discussed. The solution assumes an energy-independent radial diffusion coefficient proportional to distance and no latitudinal diffusion, and includes energy-independent particle drift velocities similar to those expected in a Parker spiral magnetic field with an equatorial current sheet. The solutions clearly demonstrate the large effects of drifts on the modulated cosmic-ray intensity. For values of the radial diffusion coefficient and particle drift velocity which are plausible for 1-GV-rigidity protons, the logarithmic radial gradient in the inner solar system is reduced by more than a factor of 5 over the value calculated in the absence of drifts. It is found that even for much smaller values of particle drift velocity and radial diffusion coefficient, such as might be expected for protons with energies of the order of 10 MeV, the effects of the drifts can be substantial.

Isenberg, P. A.↗
NASA NTRS2016

Hall2De Simulations with a First-Principles Electron Transport Model Based on the Electron Cyclotron Drift Instability

Several years of work have combined plasma measurements in a laboratory Hall thruster and r-z numerical simulations with Hall2De to isolate the spatial variation of the anomalous collision frequency needed in Ohm’s law to produce the observed thruster behavior. This numerical solution is used here to test the validity of a first-principles model of the anomalous transport in these devices before such model is implemented self-consistently in r-z fluid codes like Hall2De. The first-principles model employs quasi-linear theory and is based on the hypothesis that the Electron Cyclotron Drift Instability (ECDI) excites ion acoustic turbulence that, in turn, enhances the effective collision frequency in these devices. We find that a model of the ECDI that assumes Maxwellian velocity distributions for electrons and singly-charged, main-beam, cold ions (Ti=0.07 eV) is insufficient to explain the expected variation of the anomalous collision frequency both in the interior and exterior of the acceleration channel. When warm ions (~0.5-3 eV) are accounted for, the ECDI model in the channel interior appears more promising but fails by orders of magnitude in the near plume region due to the much higher Landau damping of the ion acoustic waves there. This implies that either (a) some process allows the ECDI instability to remain uninhibited by classical Landau damping or, (b) that a different instability (or instabilities) altogether, also insusceptible to Landau damping, is/are active in this region. A previous hypothesis, that convection of wave energy generated by the ECDI in the channel plays a significant role in the near plume, is not supported by the results of the simulations.

Mikellides, Ioannis G.↗
NASA NTRS1994

Process modelling for materials preparation experiments

The main goals of the research under this grant consist of the development of mathematical tools and measurement techniques for transport properties necessary for high fidelity modelling of crystal growth from the melt and solution. Of the tasks described in detail in the original proposal, two remain to be worked on: development of a spectral code for moving boundary problems, and development of an expedient diffusivity measurement technique for concentrated and supersaturated solutions. We have focused on developing a code to solve for interface shape, heat and species transport during directional solidification. The work involved the computation of heat, mass and momentum transfer during Bridgman-Stockbarger solidification of compound semiconductors. Domain decomposition techniques and preconditioning methods were used in conjunction with Chebyshev spectral methods to accelerate convergence while retaining the high-order spectral accuracy. During the report period we have further improved our experimental setup. These improvements include: temperature control of the measurement cell to 0.1 C between 10 and 60 C; enclosure of the optical measurement path outside the ZYGO interferometer in a metal housing that is temperature controlled to the same temperature setting as the measurement cell; simultaneous dispensing and partial removal of the lower concentration (lighter) solution above the higher concentration (heavier) solution through independently motor-driven syringes; three-fold increase in data resolution by orientation of the interferometer with respect to diffusion direction; and increase of the optical path length in the solution cell to 12 mm.

Rosenberger, Franz↗
NASA NTRS1982

Turbulent transport modelling of separating and reattaching shear flows

The improvement of capabilities for computer simulation of turbulent recirculating flows was investigated. Attention has been limited to two dimensional flows and principally to statistically stationary motion. Improvement of turbulence modeling explored the treatment of the near wall sublayer and of the exterior fully turbulent region, working within the framework of turbulence closures requiring the solution of transport equations for the turbulence energy and its dissipation rate. The work on the numerical procedure, based on the Gosman-Pun program TEACH, addressed the problems of incorporating the turbulence model as well as the extension to time dependent flows, the incorporation of a third order approximation of convective transport, and the treatment of non-orthogonal boundaries.

Launder, B. E.↗
NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗