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 127 records · Page 7

The small ice cap instability in diffusive climate models

Simple climate models which invoke diffusive heat transport and ice cap albedo feedback have equilibrium solutions with no stable ice cap smaller than a radius of about 20 deg on a great circle. Attention is presently given to a solution of this phenomenon which is physically appealing. The ice-free solution has a thermal minimum, and if the minimum temperature is just above the critical value for ice formation, then the artificial addition of a patch of ice leads to a widespread depression of the temperature below the critical freezing temperature. A second stable solution will then exist whose spatial extent is determined by the range of the influence function of a point sink of heat, due to the albedo shift in the patch.

North, G. R.↗

The vagaries of variable viscosity convection

The wavelength dependence of convection with a temperature-dependent viscosity is studied using two-dimensional finite difference models. A steady-state solution with a horizontal wavelength of eight times the depth is found to transport heat more efficiently than a solution with exactly the same parameters but with a wavelength of twice the depth. Only at the longer wavelengths is it possible for the top high viscosity boundary layer to participate more fully in the flow. Using models with a constant heat flux resulted in a decrease in the interior temperature of approximately 30% for the longer wavelength. It is also demonstrated that a parameterization using only a constant viscosity thermal structure and an interior viscosity is inadequate to parameterize the variable viscosity calculations described. This suggests that parameterized convection schemes may not accurately resolve the thermal evolution of planetary interior unless they can correctly include the effects of vertical viscosity structure.

Daly, S. F.↗

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

A stability theorem for energy-balance climate models

The paper treats the stability of steady-state solutions of some simple, latitude-dependent, energy-balance climate models. For north-south symmetric solutions of models with an ice-cap-type albedo feedback, and for the sum of horizontal transport and infrared radiation given by a linear operator, it is possible to prove a 'slope stability' theorem, i.e., if the local slope of the steady-state iceline latitude versus solar constant curve is positive (negative) the steady-state solution is stable (unstable). Certain rather weak restrictions on the albedo function and on the heat transport are required for the proof, and their physical basis is discussed.

Cahalan, R. F.↗

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗

Three-dimensional numerical modelling of transport in horizontal MOCVD reactors

A systematic three-dimensional study for the MOCVD of gallium arsenide is presented. Numerical solutions, obtained for the system of equations governing the steady state heat, momentum, and mass transfer in horizontal reactors with large and small aspect ratio, are compared with experimental growth rate distributions. The effects on layer thickness uniformity of lateral thermal boundary conditions and susceptor tilt are explored. The validity of a parabolic flow approximation of the equations is also investigated by comparison with the complete three-dimensional treatment.

Ouazzani, Jalil↗

On the explicit versus family solution of the fully diurnal photochemical equations of the stratosphere

The family approach used to solve chemical equations is examined by integrating a heirarchy of versions of a 2D model of stratospheric chemistry and transport on an isentropic surface. Errors of up to 20 percent in radical concentrations can result from the use of a small number of families; the use of more families increases the accuracy of the solutions. The minimum number of families which can be used to provide an accurate solution is discussed, and the parameterization of long-time-scale species is recommended for making the model more efficient. The most pervasive errors are found in calculations of the concentrations of radical species, a problem which can lead to errors in predicting O3 over the long run. Methods for improving the schemes are presented, emphasizing the possibility of producing an improved 3D model.

Austin, John↗

Martian thermal boundary layers: Subhourly variations induced by radiative-conductive heat transfer within the dust-laden atmosphere-ground system

Radiative-conductive heat transfer has been investigated for the ground-atmosphere system of the planet Mars. The basic goal was the quantitative determination of time dependent vertical distributions of temperature and static stability for Southern-Hemispheric summer season and middle and polar latitudes, for both dust-free and dust-laden atmospheric conditions. The numerical algorithm which models at high spatial and temporal resolution the thermal energy transports in the dual ground-atmosphere system, is based on solution of the applicable heating rate equation, including radiative and molecular-conductive heat transport terms. The two subsystems are coupled by an internal thermal boundary condition applied at the ground-atmosphere interface level. Initial data and input parameters are based on Mariner 4, 6, 7, and 9 measurements and the JPL Mars Scientific Model. Numerical experiments were run for dust-free and dust-laden conditions in the midlatitudes, as well as ice-free and ice-covered polar regions. Representative results and their interpretation are presented. Finally, the theoretical framework of the generalized problem with nonconservative Mie scattering and explicit thermal-convective heat transfer is formulated, and applicable solution algorithms are outlined.

Pallmann, A. J.↗

Transport models of the turbulent velocity-temperature products for computations of recirculating flows

Performance tests of the third-order turbulence closure for predictions of separating and recirculating flows in backward-facing steps were studied. Computations of the momentum and temperature fields in the flow domain being considered entail the solution of time-averaged transport equations containing the second-order turbulent fluctuating products. The triple products, which are responsible for the diffusive transport of the second-order products, attain greater significance in separating and reattaching flows. The computations are compared with several algebraic models and with the experimental data. The prediction was improved considerably, particularly in the separated shear layer. Computations are further made for the temperature-velocity double products and triple products. Finally, several advantages were observed in the usage of the transport equations for the evaluation of the turbulence triple products; one of the most important features is that the transport model can always take the effects of convection and diffusion into account in strong convective shear flows such as reattaching separated layers while conventional algebraic models cannot account for these effects in the evaluation of turbulence variables.

Amano, R. S.↗

Designing a Decontamination Solution for the Low-Earth-Orbit, Cryogenic SPHEREx Mission

To address the scientific goals of NASA’s astrophysics division, the JPL and Caltech’s SPHEREx mission conducts the first near-infrared all-sky spectral survey in low-earth-orbit using a passively-cooled, cryogenic telescope. Several unique water contamination issues arise due to the combination of cryogenic operation temperatures and the limited temperature control authority of passive cooling. Because water molecules have adsorption bands within the wavelengths of interest for SPHEREx’s survey, it is imperative that water contamination be controlled and minimized. In this work, a model is developed for the SPHEREx mission to predict the transport and accumulation of outgassed water onto sensitive payload components. The model utilizes the time dependent thermal profiles of the individual payload components, initial water content, and spacecraft geometry to calculate the time dependent water diffusion, transport, and adsorption. This model is subsequently used to predict water contamination risks and design decontamination solutions for each risk. Two major water contamination risks were predicted, water accumulation during cooldown and on cryogenic surfaces during the mission. The first risk occurs during the cooldown of the payload to cryogenic temperatures, predicting water accumulation on optical surfaces in excess of the allowable levels. To mitigate this accumulation, the temperature of the optical surfaces is controlled during cooldown through a combination of heaters and spacecraft pointing. The second risk is that over the course of the mission, is unavoidable water accumulation onto the payload thermal system’s cryogenic surfaces, possibly jeopardizing thermal performance and temperature stability, both of which are required for science success. To decontaminate any water that had collected onto cryogenic thermal system surfaces, a decontamination maneuver was designed. In this decontamination maneuver, cryogenic surfaces of the thermal system are warmed to a temperature where water will desorb by spacecraft pointing while still meeting all avoidance constraints, and thereby avoiding undue risk to the observatory hardware. Through application of the developed analytical model and inclusion of the decontamination maneuver in the mission design, SPHEREx can confidently demonstrate that it is able to decontaminate, at the start of and during the mission as needed to meet its end of life science performance requirements.

Soares, Carolos E.↗

Numerical modeling of pollutant transport using a Lagrangian marker particle technique

A derivation and code were developed for the three-dimensional mass transport equation, using a particle-in-cell solution technique, to solve coastal zone waste discharge problems where particles are a major component of the waste. Improvements in the particle movement techniques are suggested and typical examples illustrated. Preliminary model comparisons with analytic solutions for an instantaneous point release in a uniform flow show good results in resolving the waste motion. The findings to date indicate that this computational model will provide a useful technique to study the motion of sediment, dredged spoils, and other particulate waste commonly deposited in coastal waters.

Spaulding, M.↗

Process modelling for materials preparation experiments

The development is examined of mathematical tools and measurement of transport properties necessary for high fidelity modeling of crystal growth from the melt and solution, in particular for the Bridgman-Stockbarger growth of mercury cadmium telluride (MCT) and the solution growth of triglycine sulphate (TGS). The tasks include development of a spectral code for moving boundary problems, kinematic viscosity measurements on liquid MCT at temperatures close to the melting point, and diffusivity measurements on concentrated and supersaturated TGS solutions. A detailed description is given of the work performed for these tasks, together with a summary of the resulting publications and presentations.

Rosenberger, Franz↗

Process modelling for materials preparation experiments

The main goals of the research under this grant consist of the development of mathematical tools and measurement of transport properties necessary for high fidelity modeling of crystal growth from the melt and solution, in particular, for the Bridgman-Stockbarger growth of mercury cadmium telluride (MCT) and the solution growth of triglycine sulphate (TGS). Of the tasks described in detail in the original proposal, two remain to be worked on: (1) development of a spectral code for moving boundary problems; and (2) diffusivity measurements on concentrated and supersaturated TGS solutions. Progress made during this seventh half-year period is reported.

Rosenberger, Franz↗

Process modelling for materials preparation experiments

The main goals of the research consist of the development of mathematical tools and measurement of transport properties necessary for high fidelity modeling of crystal growth from the melt and solution, in particular for the Bridgman-Stockbarger growth of mercury cadmium telluride (MCT) and the solution growth of triglycine sulphate (TGS). 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 diffusivity measurements on concentrated and supersaturated TGS solutions. During this eighth half-year period, good progress was made on these tasks.

Rosenberger, Franz↗

Analysis of Residual Acceleration Effects on Transport and Segregation During Directional Solidification of Tin-Bismuth in the MEPHISTO Furnace Facility

The objective of this work is to approach the problem of determining the transport conditions (and effects of residual acceleration) during the plane-front directional solidification of a tin-bismuth alloy under low gravity conditions. The work involves using a combination of 2- and 3-D numerical models, scaling analyses, ID models and the results of ground-based and low-gravity experiments. The latter are to be conducted during the MEPHISTO experiment scheduled for USMP-3 in early 1996. The models will be used to predict the response of the transport conditions and consequent solute segregation in directionally solidifying tin-bismuth melt. Real-time Seebeck voltage variations across a Sn-Bi melt during directional solidification in MEPHISTO on USMP-1 show a distinct variation which can be correlated with thruster firings. The Seebeck voltage measurement is related to the response of the instantaneous average melt composition at the melt-solid interface. This allows a direct comparison of numerical simulations with the Seebeck signals obtained on USMP-1. The effects of such accelerations on composition for a directionally solidifying Sn-Bi alloy have been simulated numerically. USMP-1 acceleration data was used to assist in our choice of acceleration magnitude and orientation. The results show good agreement with experimental observations. The USMP-3 experiments took place earlier this year (February 22 through March 6). There were several differences between the USMP-3 experiments as compared to USMP-1. Firstly a more concentrated alloy was solidified and, secondly, Primary Reaction Control System thruster burns were requested at particular times during four separate growth runs. This allowed us to monitor the response Seebeck response under well-characterized growth conditions. In addition, we carried out simulations during the experiment in order to interpret the Seebeck signal. Preliminary results are described here.

Alexander J. Iwan D.↗

Performance of a Bounce-Averaged Global Model of Super-Thermal Electron Transport in the Earth's Magnetic Field

In this paper, we report the results of our recent research on the application of a multiprocessor Cray T916 supercomputer in modeling super-thermal electron transport in the earth's magnetic field. In general, this mathematical model requires numerical solution of a system of partial differential equations. The code we use for this model is moderately vectorized. By using Amdahl's Law for vector processors, it can be verified that the code is about 60% vectorized on a Cray computer. Speedup factors on the order of 2.5 were obtained compared to the unvectorized code. In the following sections, we discuss the methodology of improving the code. In addition to our goal of optimizing the code for solution on the Cray computer, we had the goal of scalability in mind. Scalability combines the concepts of portabilty with near-linear speedup. Specifically, a scalable program is one whose performance is portable across many different architectures with differing numbers of processors for many different problem sizes. Though we have access to a Cray at this time, the goal was to also have code which would run well on a variety of architectures.

McGuire, Tim↗

Survey of Turbulence Models for the Computation of Turbulent Jet Flow and Noise

The report presents an overview of jet noise computation utilizing the computational fluid dynamic solution of the turbulent jet flow field. The jet flow solution obtained with an appropriate turbulence model provides the turbulence characteristics needed for the computation of jet mixing noise. A brief account of turbulence models that are relevant for the jet noise computation is presented. The jet flow solutions that have been directly used to calculate jet noise are first reviewed. Then, the turbulent jet flow studies that compute the turbulence characteristics that may be used for noise calculations are summarized. In particular, flow solutions obtained with the k-e model, algebraic Reynolds stress model, and Reynolds stress transport equation model are reviewed. Since, the small scale jet mixing noise predictions can be improved by utilizing anisotropic turbulence characteristics, turbulence models that can provide the Reynolds stress components must now be considered for jet flow computations. In this regard, algebraic stress models and Reynolds stress transport models are good candidates. Reynolds stress transport models involve more modeling and computational effort and time compared to algebraic stress models. Hence, it is recommended that an algebraic Reynolds stress model (ASM) be implemented in flow solvers to compute the Reynolds stress components.

Nallasamy, N.↗

Start-On-The-Part Transient Model for In-Situ Automated Tape Placement of Thermoplastic Composites

Fabrication of a complex part by automated tape placement (ATP) can require starting up a new tape-end in the part interior, termed start-on-the-part. Careful thermal management of the starting transient is needed to achieve uniform crystallinity and inter-laminar weld strength - which is the objective of this modeling effort. The transient is modeled by a Fourier-Laplace transform solution of the time-dependent thermal transport equation in two spatial dimensions. The solution is subject to a quasi-steady approximation for the speed and length of the consolidation head. Sample calculations are done for the Langley ATP robot applying PEEK/carbon fiber composite and for two upgrades in robot performance. The head starts out almost at rest which meets an engineering requirement for accurate placement of the new tape-end. The head then rapidly accelerates until it reaches its steady state speed. This rapid acceleration, however, violates the quasi-steady approximation, so uniform weld strength and crystallinity during the starting transient are not actually achieved. The solution does give the elapsed time and distance from start-up to validity of the quasi-steady approximation - which quantifies the length of the non-uniform region. The elapsed time was always less than 0.1 s and the elapsed distance less than 1 cm. This quantification would allow the non-uniform region to be either trimmed away or compensated for in the design of a part. Such compensation would require experiments to measure the degree of non-uniformity, because the solution does not provide this information. The rapid acceleration suggests that the consolidation roller or belt be actively synchronized to avoid abrading the tape.

Costen, Robert c.↗