Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convection diffusion equation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Analogies Between Colloidal Sedimentation and Turbulent Convection at High Prandtl Numbers

A new set of coarse-grained equations of motion is proposed to describe concentration and velocity fluctuations in a dilute sedimenting suspension of non-Brownian particles. With these equations, colloidal sedimentation is found to be analogous to turbulent convection at high Prandtl numbers. Using Kraichnan's mixing-length theory, we obtain scaling relations for the diffusive dissipation length delta(sub theta), the velocity variance delta u, and the concentration variance delta phi. The obtained scaling laws over varying particle radius alpha and volume fraction phi(sub ) are in excellent agreement with the recent experiment by Segre, Herbolzheimer, and Chaikin. The analogy between colloidal sedimentation and turbulent convection gives a simple interpretation for the existence of a velocity cut-off length, which prevents hydrodynamic dispersion coefficients from being divergent. It also provides a coherent framework for the study of sedimentation dynamics in different colloidal systems.

Tong, P.↗
NASA NTRS1981

State of the art in subgrid scale modeling

The present state of the art in subgrid scale modeling for large eddy simulation is reviewed. It is found that the eddy viscosity models are sufficient for simulating homogeneous flows but are probably insufficient for inhomogeneous shear flows. One- and two-equation subgrid scale models probably will not be of great value in subgrid scale modeling, and full Reynolds stress models may be required. In the latter, the convection and diffusion terms are expected to be of greater importance than they are in time-average modeling, but the pressure redistribution terms may be of less importance. A scale similarity model may provide an alternative.

Ferziger, J. H.↗
NASA NTRS1985

The equilibrium structure of thin magnetic flux tubes. I

A model atmosphere within a thin magnetic flux tube that is embedded in an arbitrarily stratified medium is presently constructed by solving the radiative transfer equation in the two-stream approximation for gray opacity, under the assumption that the magnetic field is sufficiently strong to warrant the neglect of both thermal conduction and convective diffusion; energy inside the flux tube therefore being transported solely by radiation. The structure of the internal atmosphere is determined on the basis of the hydrostatic and radiative equilibrium conditions of the tube embedded in an external atmosphere. The gas temperature along the axis of the tube is computed, and the geometry of the flux tube is determined on the basis of magnetic flux conservation and total pressure equilibrium.

Ferrari, A.↗
NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗
NASA NTRS1976

Solutions of the Navier-Stokes equations for vortex breakdown

A numerical investigation of vortex breakdown has been undertaken in an attempt to understand its properties, and the mechanisms responsible for it. Solutions of the full steady axisymmetric Navier-Stokes equations for breakdown in an unconfined viscous vortex have been obtained for core Reynolds numbers up to 200, for a two-parameter family of assumed upstream velocity distributions. Diffusion and convection of vorticity away from the vortex core, and the strong coupling between the circumferential and axial velocity fields in highly-swirling flows, are shown to lead to stagnation and reversal of the axial flow near the axis. The various theories of vortex breakdown are considered in light of the present numerical solutions.

Grabowski, W. J.↗
NASA NTRS1978

Numerical method for boundary layers with blowing - The exponential box scheme

The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.

El-Mistikawy, T. M.↗
NASA NTRS1993

Flux-difference split parabolized Navier-Stokes algorithm for non-equilibrium chemically reacting flows

A flux-difference split explicit finite-difference algorithm is presented for solving the parabolized form of the equations governing three-dimensional nonequilibrium chemically reacting flows. The algorithm is based on an explicit noniterative, upwind space-marching scheme developed by Korte, but differs in that the unsteady Riemann problem, rather than the steady Riemann problem, is solved. The algorithm allows either a second or an approximately third-order accurate upwind treatment of the convection terms by employing the unsteady approximate Riemann solver of Roe. The source terms of the species transport equations are treated in either an explicit or implicit manner, and the species diffusion terms are modeled with either a Fickian or a multicomponent model. A validation of the algorithm is performed by comparing computational results with the 2-D Mach 14, 15 degree compression-corner data of Holden. The three-dimensional capability of the algorithm is demonstrated by computing Mach 2.7 flow over a swept wedge scramjet fuel injector, and three-dimensional reacting flow capability is demonstrated by a computing a shock-jet interaction concept for mixing and combustion enhancement.

White, J. A.↗
NASA NTRS2005

Removal of organic contaminants by RO and NF membranes

Rejection characteristics of organic and inorganic compounds were examined for six reverse osmosis (RO) membranes and two nanofiltration (NF) membranes that are commercially available. A batch stirred-cell was employed to determine the membrane flux and the solute rejection for solutions at various concentrations and different pH conditions. The results show that for ionic solutes the degree of separation is influenced mainly by electrostatic exclusion, while for organic solutes the removal depends mainly upon the solute radius and molecular structure. In order to provide a better understanding of rejection mechanisms for the RO and NF membranes, the ratio of solute radius (r(i,s)) to effective membrane pore radius (r(p)) was employed to compare rejections. An empirical relation for the dependence of the rejection of organic compounds on the ratio r(i,s)/r(p) is presented. The rejection for organic compounds is over 75% when r(i,s)/r(p) is greater than 0.8. In addition, the rejection of organic compounds is examined using the extended Nernst-Planck equation coupled with a steric hindrance model. The transport of organic solutes is controlled mainly by diffusion for the compounds that have a high r(i,s)/r(p) ratio, while convection is dominant for compounds that have a small r(i,s)/r(p) ratio. c2005 Elsevier B.V. All rights reserved.

NASA Discipline Life Sciences Technologies↗
NASA NTRS1993

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗
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.↗
DOE OSTI2023

Learning subgrid-scale models with neural ordinary differential equations

We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.

97 MATHEMATICS AND COMPUTING↗
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

Vapor Transport Within the Thermal Diffusion Cloud Chamber

A review of the equations used to determine the 1-D vapor transport in the thermal diffusion cloud chamber (TDCC) is presented. These equations closely follow those of the classical Stefan tube problem in which there is transport of a volatile species through a noncondensible, carrier gas. In both cases, the very plausible assumption is made that the background gas is stagnant. Unfortunately, this assumption results in a convective flux which is inconsistent with the momentum and continuity equations for both systems. The approximation permits derivation of an analytical solution for the concentration profile in the Stefan tube, but there is no computational advantage in the case of the TDCC. Furthermore, the degree of supersaturation is a sensitive function of the concentration profile in the TD CC and the stagnant background gas approximation can make a dramatic difference in the calculated supersaturation. In this work, the equations typically used with a TDCC are compared with very general transport equations describing the 1-D diffusion of the volatile species. Whereas no pressure dependence is predicted with the typical equations, a strong pressure dependence is present with the more general equations given in this work. The predicted behavior is consistent with observations in diffusion cloud experiments. It appears that the new equations may account for much of the pressure dependence noted in TDCC experiments, but a comparison between the new equations and previously obtained experimental data are needed for verification.

Ferguson, Frank T.↗
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.↗
NASA NTRS1986

Diffusion flame extinction in slow convection 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 Damkchler 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, Chiun-Hsun↗
NASA NTRS1992

Intensity variations in the interplanetary magnetic field measured by Voyager 2 and the 11-year solar cycle modulation of galactic cosmic rays

New evidence is presented to support the hypothesis that the 11-yr solar cycle modulation of galactic cosmic rays is caused by strong diffusion inside long-lived merged interaction regions. To test this hypothesis, the 1D force-field approximation of the cosmic ray modulation equation is solved. It is assumed that a constant solar wind speed convects magnetic field compressions and rarefactions unchanged through a model heliosphere. The result is a reasonable simulation of the integrated high-energy cosmic ray intensity profile from about 1982 to mid-1989. This period encompasses both the full recovery portion of the last profile from about 1982 to mid-1989. This model responds to the Voyager 2 magnetic field data by correctly timing the beginning of the new modulation cycle in late 1987. It is concluded that the present hypothesis is consistent with the results of this simulation.

Perko, J. S.↗
NASA NTRS1975

Scatter-free collimated convection and cosmic-ray transport at 1 AU

Observations of solar cosmic rays between 0.1 MeV and 1 GeV are more consistent with scatter-free than diffusive propagation near 1 AU. Since these particles also are constrained to interplanetary magnetic field lines, the parallel and transverse flow are decoupled by collimated convection. Local propagation is therefore more appropriately described by Liouville's equation, and transport properties can be examined from the exact zeroeth and first moments of the equation averaged over all directions in momentum space. Observational results which may be explained directly without invoking local scattering include the decay of solar flare events as well as the radial gradient, diurnal and semidiurnal variation of galactic cosmic rays.

Roelof, E. C.↗
NASA NTRS1986

Vapor transport processes of Ge-chalcogenides under normal and reduced gravity conditions

Experimental and theoretical studies of the GeSe-GeI4 system revealed the multicomponent, multireaction nature of this system and the existence of a diffusion boundary layer; they also yielded equations for the prediction of diffusion limited mass transport rates. Vapor transport and crystal growth experiments of the GeSe-GeI4 and GeSe-Xe systems on earth and in a microgravity environment demonstrated the effects of convection and of microgravity on mass flux and crystal morphology. The surface and bulk morphology of GeSe crystals obtained for both systems in microgravity is considerably improved relative to ground-based specimens. Unexpected crystal growth phenomena were observed for the GeSe-Xenon system on recent Shuttle flights.

Wiedemeier, H.↗