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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Finite-amplitude thermal convection in a spherical shell

The properties of finite-amplitude thermal convection for a Boussinesq fluid contained in a spherical shell are investigated. All nonlinear terms are retained in the equations, and both axisymmetric and nonaxisymmetric solutions are studied. The velocity is expanded in terms of poloidal and toroidal vectors. Spherical surface harmonics resolve the horizontal structure of the flow, but finite differences are used in the vertical. With a few modifications, the transform method developed by Orszag (1970) is used to calculate the nonlinear terms, while Green's function techniques are applied to the poloidal equation and diffusion terms.

Young, R. E.↗

DOE OSTI2021

A vortex sheet based analytical model of the curled wake behind yawed wind turbines

Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1981

Zero-gravity aerosol behavior

The feasibility and scientific benefits of a zero gravity aerosol study in an orbiting laboratory were examined. A macroscopic model was devised to deal with the simultaneous effects of diffusion and coagulation of particles in the confined aerosol. An analytical solution was found by treating the particle coagulation and diffusion constants as ensemble parameters and employing a transformation of variables. The solution was used to carry out simulated zero gravity aerosol decay experiments in a compact cylindrical chamber. The results demonstrate that the limitations of physical space and time imposed by the orbital situation are not prohibitive in terms of observing the history of an aerosol confined under zero gravity conditions. While the absence of convective effects would be a definite benefit for the experiment, the mathematical complexity of the problem is not greatly reduced when the gravitational term drops out of the equation. Since the model does not deal directly with the evolution of the particle size distribution, it may be desirable to develop more detailed models before undertaking an orbital experiment.

Edwards, H. W.↗

NASA NTRS1974

Solutions of the Navier-Stokes equations for vortex breakdown

Steady solutions of the Navier-Stokes equations, in terms of velocity and pressure, for breakdown in an unconfined viscous vortex are obtained numerically using the artificial compressibility technique of Chorin combined with an ADI finite-difference scheme. Axisymmetry is assumed and boundary conditions are carefully applied at the boundaries of a large finite region in an axial plane while resolution near the axis is maintained by a coordinate transformation. The solutions, which are obtained for Reynolds numbers up to 200 based on the free-stream axial velocity and a characteristic core radius, show that breakdown results from the diffusion and convection of vorticity away from the vortex core which, because of the strong coupling between the circumferential and axial velocity fields in strongly swirling flows, can lead to stagnation and reversal of the axial flow near the axis.

Grabowski, W. J.↗

NASA NTRS1978

Turbulent convection within rapidly rotating superadiabatic fluids with horizontal temperature gradients

A model is developed for small scale thermally driven turbulent convection in rapidly rotating systems. The Boussinesq equations of adiabatic and inviscid flow are perturbed about a basic state which is superadiabatic and which contains a horizontal temperature gradient balanced by a geostrophic wind. The growth rates and length scales of the most unstable modes are determined. It is found that these modes are either zonally or latitudinally symmetric, depending on the latitude and the horizontal temperature gradient. The transports of heat and momentum by the modes are computed. The linear amplitudes are assumed to be limited by shear instability. The model predicts a latitude dependent eddy diffusivity. The applicability of the model to astrophysical and planetary bodies is briefly discussed.

Flasar, F. M.↗

NASA NTRS1994

Steepened structures in equatorial spread F. 2: Theory

This paper investigates the properties of a one-dimensional fluid model of plasma convection in the equatorial F region ionosphere. The model equations are similar in form to Burgers equation except for additional higher-order spatial derivatives. Like Burgers equation, solution to the model have the form of propagating, shocklike structures. Numerical simulations of the model closely resemble the steepened structures observed by sounding rocket plasma density probes within equatorial spread F. Simulated denstiy power spectra, like the spectra computed from in situ data, seem to possess power law forms with a break at wavelengths of about 100 m. The precise wavenumber of the spectral break is determined by the ambipolar diffusion coefficient. The model predicts that electric field fluctuations perpendicular to the direction of plasma steepening should be proportional to the plasma density fluctuations. Electric field fluctuations parallel to the steepening will be due primarily to the ambipolar field and have a Boltzmann relationship with density (square of the absolute value of delta E) approximately equal to (K(exp 2))(square of the absolute value of (delta n/n)). At wavelengths less than about 300 m, the ambipolar field should be the dominant component of the total field intensity.

Hysell, D. L↗

NASA NTRS1987

Energetic particle acceleration in spherically symmetric accretion flows and shocks

Steady state, spherically symmetric solutions of the cosmic-ray transport equation describing the acceleration of energetic particles in galactic accretion flows onto neutron stars, black holes, white dwarfs, and protostars are studied. The results indicate that astrophysical accretion flows can be partitioned into distinct classes depending upon whether the accretion rate lies above or below a critical value of a few times 10 to the -7th stellar masses/yr. When the particle transport is convection-dominated, both classes of accretion flows exhibit a spectral index appropriate for first-order Fermi acceleration at a plane shock in the absence of losses. As the particle transport becomes diffusion-dominated, both classes show a break and precipitous falloff in the particle spectrum due to the escape of these particles from the accretion flow. The precise nature of the spectrum depends on the relationship between the particle momentum and the spatial diffusion coefficient.

Webb, G. M.↗

NASA NTRS1989

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Yung, Chain-Nan↗

NASA NTRS1994

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗

NASA NTRS2003

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

DOE OSTI2021

Magnetic-field generation and its effect on ablative Rayleigh–Taylor instability in diffusive ablation fronts

The effects of self-generated magnetic fields on the ablative Rayleigh–Taylor (RT) instability are investigated in the linear regime. The main governing parameters are the Froude number (Fr), which stands for the ratio between ablative convection and acceleration of the target, and the Mach number at the ablation front (Ma), assumed to be small (isobaricity). During the development of the RT instability, magnetic fields are generated due to misalignment between pressure and density gradients (Biermann-battery effect). They accumulate at the section of the ablation front where the Nernst and the plasma velocities cancel each other. The magnetic field modifies the dynamics of the instability through the Righi–Leduc term, which acts as a heat source in the energy equation. It is found that the B fields affect perturbations with short wavelengths up to the most unstable wave in the spectrum. Here, the B field plays a destabilizing role for moderate Froude numbers and becomes stabilizing for large Froude numbers. For plastic ablators, the Fr threshold is found to be Fr = 5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1987

A comparison of numerical flux formulas for the Euler and Navier-Stokes equations

Numerical flux formulas for the convection terms in the Euler or Navier-Stokes equations are analyzed with regard to their accuracy in representing steady nonlinear and linear waves (shocks and entropy/shear waves, respectively). Numerical results are obtained for a one-dimensional conical Navier-Stokes flow including both a shock and a boundary layer. Analysis and experiments indicate that for an accurate representation of both layers the flux formula must include information about all different waves by which neighboring cells interact, as in Roe's flux-difference splitting. In comparison, Van Leer's flux-vector splitting, which ignores the linear waves, badly diffuses the boundary layer. The results of MacCormack's scheme, if properly tuned, are significantly better. The use of a sufficiently detailed flux formula appears to reduce the number of cells required to resolve a boundary layer by a factor 1/2 to 1/4 and thus pays off.

Van Leer, Bram↗

NASA NTRS1976

Similarity theory of the buoyantly interactive planetary boundary layer with entrainment

A similarity model is developed for the vertical profiles of turbulent flow variables in an entraining turbulent boundary layer of arbitrary buoyant stability. In the general formulation the vertical profiles, internal rotation of the velocity vector, discontinuities or jumps at a capping inversion and bulk aerodynamic coefficients of the boundary layer are given by solutions to a system of ordinary differential equations in the similarity variable. To close the system, a formulation for buoyantly interactive eddy diffusivity in the boundary layer is introduced which recovers Monin-Obukhov similarity near the surface and incorporates a hypothesis accounting for the observed variation of mixing length throughout the boundary layer. The model is tested in simplified versions which depend only on roughness, surface buoyancy, and Coriolis effects by comparison with planetary-boundary-layer wind- and temperature-profile observations, measurements of flat-plate boundary layers in a thermally stratified wind tunnel and observations of profiles of terms in the turbulent kinetic-energy budget of convective planetary boundary layers. On balance, the simplified model reproduced the trend of these various observations and experiments reasonably well, suggesting that the full similarity formulation be pursued further.

Hoffert, M. I.↗

NASA NTRS2001

Our Sun IV: The Standard Model and Helioseismology: Consequences of Uncertainties in Input Physics and in Observed Solar Parameters

Helioseismic frequency observations provide an extremely accurate window into the solar interior; frequencies from the Michaelson Doppler Imager (MDI) on the Solar and Heliospheric Observatory (SOHO) spacecraft, enable the adiabatic sound speed and adiabatic index to be inferred with an accuracy of a few parts in 10(exp 4) and the density with an accuracy of a few parts in 10(exp 3). This has become a Serious challenge to theoretical models of the Sun. Therefore, we have undertaken a self-consistent, systematic study of the sources of uncertainties in the standard solar models. We found that the largest effect on the interior structure arises from the observational uncertainties in the photospheric abundances of the elements, which affect the sound speed profile at the level of 3 parts in 10(exp 3). The estimated 4% uncertainty in the OPAL opacities could lead to effects of 1 part in 10(exp 3); the approximately 5%, uncertainty in the basic pp nuclear reaction rate would have a similar effect, as would uncertainties of approximately 15% in the diffusion constants for the gravitational settling of helium. The approximately 50% uncertainties in diffusion constants for the heavier elements would have nearly as large an effect. Different observational methods for determining the solar radius yield results differing by as much as 7 parts in 10(exp 4); we found that this leads to uncertainties of a few parts in 10(exp 3) in the sound speed int the solar convective envelope, but has negligible effect on the interior. Our reference standard solar model yielded a convective envelope position of 0.7135 solar radius, in excellent agreement with the observed value of 0.713 +/- 0.001 solar radius and was significantly affected only by Z/X, the pp rate, and the uncertainties in helium diffusion constants. Our reference model also yielded envelope helium abundance of 0.2424, in good agreement with the approximate range of 0.24 to 0.25 inferred from helioseismic observations; only extreme Z/X values yielded envelope helium abundance outside this range. We found that other current uncertainties, namely, in the solar age and luminosity, in nuclear rates other than the pp reaction, in the low-temperature molecular opacities, and in the low-density equation of state, have no significant effect on the quantities that can be inferred from helioseismic observations. The predicted pre-main-sequence lithium depletion is uncertain by a factor of 2. The predicted neutrino capture rate is uncertain by approximately 30% for the Cl-27 experiment and by approximately 3% for Ga-71 experiments, while the B-8 neutrino flux is uncertain by approximately 30%.

Boothroyd, Arnold I.↗

NASA NTRS1972

An experimental and theoretical evaluation of increased thermal diffusivity phase change devices

This study was to experimentally evaluate and mathematically model the performance of phase change thermal control devices containing high thermal conductivity metal matrices. Three aluminum honeycomb filters were evaluated at five different heat flux levels using n-oct-adecane as the test material. The system was mathematically modeled by approximating the partial differential equations with a three-dimensional implicit alternating direction technique. The mathematical model predicts the system quite well. All of the phase change times are predicted. The heating of solid phase is predicted exactly while there is some variation between theoretical and experimental results in the liquid phase. This variation in the liquid phase could be accounted for by the fact that there are some heat losses in the cell and there could be some convection in the experimental system.

White, S. P.↗

DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

NASA NTRS1997

Instability Mechanisms of Thermally-Driven Interfacial Flows in Liquid-Encapsulated Crystal Growth

During the past year, a great deal of effort was focused on the enhancement and refinement of the computational tools developed as part of our previous NASA grant. In particular, the interface mollification algorithm developed earlier was extended to incorporate the effects of surface-rheological properties in order to allow the study of thermocapillary flows in the presence of surface contamination. These tools will be used in the computational component of the proposed research in the remaining years of this grant. A detailed description of the progress made in this area is provided elsewhere. Briefly, the method developed allows for the convection and diffusion of bulk-insoluble surfactants on a moving and deforming interface. The novelty of the method is its grid independence: there is no need for front tracking, surface reconstruction, body-fitted grid generation, or metric evaluations; these are all very expensive computational tasks in three dimensions. For small local radii of curvature there is a need for local grid adaption so that the smearing thickness remains a small fraction of the radius of curvature. A special Neumann boundary condition was devised and applied so that the calculated surfactant concentration has no variations normal to the interface, and it is hence truly a surface-defined quantity. The discretized governing equations are solved subsequently using a time-split integration scheme which updates the concentration and the shape successively. Results demonstrate excellent agreement between the computed and exact solutions.

Haj-Hariri, Hossein↗

NASA NTRS1990

Chemical heterogeneity and imperfect mixing in the solar nebula

The 'standard' accretion-disk model of the solar nebula is examined together with its basic equations, which are solved for the steady-state case. It is shown that, for this model of the solar nebula, the efficiency of mixing the products of thermochemical processing from small to large disk radii depends not on the magnitude of the eddy diffusivity but only on the ratio k between the eddy diffusivity and the eddy viscosity. In most cases, this ratio k is probably of order of unity or less, implying that most of the solar nebula was not contaminated by the consequences of thermochemical equilibria that were established at 'small' disk radii. It is noted that, while this condition is almost certainly satisfied if the physical process responsible for viscosity is thermal convection or waves, it has uncertain validity during the early phases of disk evolution where accretion-induced shear instabilities may dominate.

Stevenson, David J.↗