Search NASA⌕ Search

SEARCH · Search NASA

Results for “energy equations”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗
NASA NTRS1993

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 3: Programmer's reference

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗
NASA NTRS1993

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗
NASA NTRS1996

Measurements and Modeling of Soot Formation and Radiation in Microgravity Jet Diffusion Flames

This is a computational and experimental study for soot formation and radiative heat transfer in jet diffusion flames under normal gravity (1-g) and microgravity (0-g) conditions. Instantaneous soot volume fraction maps are measured using a full-field imaging absorption technique developed by the authors. A compact, self-contained drop rig is used for microgravity experiments in the 2.2-second drop tower facility at NASA Lewis Research Center. On modeling, we have coupled flame structure and soot formation models with detailed radiation transfer calculations. Favre-averaged boundary layer equations with a k-e-g turbulence model are used to predict the flow field, and a conserved scalar approach with an assumed Beta-pdf are used to predict gaseous species mole fraction. Scalar transport equations are used to describe soot volume fraction and number density distributions, with formation and oxidation terms modeled by one-step rate equations and thermophoretic effects included. An energy equation is included to couple flame structure and radiation analyses through iterations, neglecting turbulence-radiation interactions. The YIX solution for a finite cylindrical enclosure is used for radiative heat transfer calculations. The spectral absorption coefficient for soot aggregates is calculated from the Rayleigh solution using complex refractive index data from a Drude- Lorentz model. The exponential-wide-band model is used to calculate the spectral absorption coefficient for H20 and C02. It is shown that when compared to results from true spectral integration, the Rosseland mean absorption coefficient can provide reasonably accurate predictions for the type of flames studied. The soot formation model proposed by Moss, Syed, and Stewart seems to produce better fits to experimental data and more physically sound than the simpler model by Khan et al. Predicted soot volume fraction and temperature results agree well with published data for a normal gravity co-flow laminar flames and turbulent jet flames. Predicted soot volume fraction results also agree with our data for 1-g and 0-g laminar jet names as well as 1-g turbulent jet flames.

Ku, Jerry C.↗
NASA NTRS1997

A Thermal and Electrical Analysis of Power Semiconductor Devices

The state-of-art power semiconductor devices require a thorough understanding of the thermal behavior for these devices. Traditional thermal analysis have (1) failed to account for the thermo-electrical interaction which is significant for power semiconductor devices operating at high temperature, and (2) failed to account for the thermal interactions among all the levels involved in, from the entire device to the gate micro-structure. Furthermore there is a lack of quantitative studies of the thermal breakdown phenomenon which is one of the major failure mechanisms for power electronics. This research work is directed towards addressing. Using a coupled thermal and electrical simulation, in which the drift-diffusion equations for the semiconductor and the energy equation for temperature are solved simultaneously, the thermo-electrical interactions at the micron scale of various junction structures are thoroughly investigated. The optimization of gate structure designs and doping designs is then addressed. An iterative numerical procedure which incorporates the thermal analysis at the device, chip and junction levels of the power device is proposed for the first time and utilized in a BJT power semiconductor device. In this procedure, interactions of different levels are fully considered. The thermal stability issue is studied both analytically and numerically in this research work in order to understand the mechanism for thermal breakdown.

Vafai, Kambiz↗
NASA NTRS2001

Plasma Reactor Modeling and Validation Experiments

Plasma processing is a key processing stop in integrated circuit manufacturing. Low pressure, high density plum reactors are widely used for etching and deposition. Inductively coupled plasma (ICP) source has become popular recently in many processing applications. In order to accelerate equipment and process design, an understanding of the physics and chemistry, particularly, plasma power coupling, plasma and processing uniformity and mechanism is important. This understanding is facilitated by comprehensive modeling and simulation as well as plasma diagnostics to provide the necessary data for model validation which are addressed in this presentation. We have developed a complete code for simulating an ICP reactor and the model consists of transport of electrons, ions, and neutrals, Poisson's equation, and Maxwell's equation along with gas flow and energy equations. Results will be presented for chlorine and fluorocarbon plasmas and compared with data from Langmuir probe, mass spectrometry and FTIR.

Meyyappan, M.↗
NASA NTRS2005

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Numbers

A new turbulence model suited for calculating the turbulent Prandtl number as part of the solution is presented. The model is based on a set of two equations: one governing the variance of the enthalpy and the other governing its dissipation rate. These equations were derived from the exact energy equation and thus take into consideration compressibility and dissipation terms. The model is used to study two cases involving shock wave/boundary layer interaction at Mach 9.22 and Mach 5.0. In general, heat transfer prediction showed great improvement over traditional turbulence models where the turbulent Prandtl number is assumed constant. It is concluded that using a model that calculates the turbulent Prandtl number as part of the solution is the key to bridging the gap between theory and experiment for flows dominated by shock wave/boundary layer interactions.

Gaffney, R. L., Jr.↗
NASA NTRS2007

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Numbers

A new turbulence model suited for calculating the turbulent Prandtl number as part of the solution is presented. The model is based on a set of two equations: one governing the variance of the enthalpy and the other governing its dissipation rate. These equations were derived from the exact energy equation and thus take into consideration compressibility and dissipation terms. The model is used to study two cases involving shock wave/boundary layer interaction at Mach 9.22 and Mach 5.0. In general, heat transfer prediction showed great improvement over traditional turbulence models where the turbulent Prandtl number is assumed constant. It is concluded that using a model that calculates the turbulent Prandtl number as part of the solution is the key to bridging the gap between theory and experiment for flows dominated by shock wave/boundary layer interactions.

Xiao, X.↗
NASA NTRS1983

MHD stability of incompressible coronal loops with radiative energy loss

Previous studies of the magnetohydrodynamic (MHD) stability of solar coronal loops have not taken into account the effects of radiative energy loss in the energy equation. However, since coronal loops continuously lose energy by radiation and heat conduction, it is important to understand how these energy loss mechanisms affect MHD stability. We investigate the problem assuming that a magnetic loop has cylindrical geometry. As a first step, stability is studied for a localized mode, and the result is applied to a specific equilibrium. We find that the radiative energy loss effect not only changes the growth rate of ideally unstable modes, but also alters the stability boundary predicted by ideal MHD theory.

An, C.-H.↗
NASA NTRS1988

Fluid-thermal-structural study of aerodynamically heated leading edges

A finite element approach for integrated fluid-thermal-structural analysis of aerodynamically heated leading edges is presented. The Navier-Stokes equations for high speed compressible flow, the energy equation, and the quasi-static equilibrium equations for the leading edge are solved using a single finite element approach in one integrated, vectorized computer program called LIFTS. The fluid-thermal-structural coupling is studied for Mach 6.47 flow over a 3-in diam cylinder for which the flow behavior and the aerothermal loads are calibrated by experimental data. Issues of the thermal-structural response are studied for hydrogen-cooled, super thermal conducting leading edges subjected to intense aerodynamic heating.

Deuchamphai, Pramote↗
NASA NTRS1988

Flow-thermal-structural study of aerodynamically heated leading edges

A finite element approach for integrated fluid-thermal-structural analysis of aerodynamically heated leading edges is presented. The Navier-Stokes equations for high speed compressible flow, the energy equation, and the quasi-static equilibrium equations for the leading edge are solved using a single finite element approach in one integrated, vectorized computer program called LIFTS. The fluid-thermal-structural coupling is studied for Mach 6.47 flow over a 3-inch diameter cylinder for which the flow behavior and the aerothermal loads are calibrated by experimental data. Issues of the thermal-structural response are studied for hydrogen cooled, super thermal conducting leading edges subjected to intense aerodynamic heating.

Dechaumphai, Pramote↗
NASA NTRS1986

Efficient numerical techniques for complex fluid flows

The central feature in any flow prediction method is the treatment of the coupling between the momentum and continuity equations. In natural-convection flows, the energy equation also becomes strongly coupled with the momentum equations. Because of the nonlinear nature of the coupling, these equations are solved iteratively. Iterative methods are often prone to slow convergence, divergence, and extreme sensitivity to underrelaxation factors. The aim of the present research is to develop more efficient and reliable solution schemes for the coupled flow equations. Such schemes will significantly reduce the expense of computing complex flows encountered in combustion chambers, gas turbines, heat exchangers, and other practical equipment. In the work completed so far, a technique employing norm reduction in conjunction with the successive-substitution and Newton-Raphson techniques was developed. Also, a block-correction procedure for the flow equations is currently being formulated and tested.

Patankar, Suhas V.↗
NASA NTRS1989

Fluid-thermal-structural interaction of aerodynamically heated leading edges

A two-dimensional finite element approach is presented for the integrated fluid-thermal-structural analysis of aerodynamically heated leading edges. The approach is combined with an adaptive unstructured remeshing technique to solve the Navier-Stokes equations for high speed compressible flow, the energy equation for the structure thermal response, and the quasi-static equilibrium equations for the structural response. Coupling and interaction between the three disciplines are demonstrated using two applications for high speed flow over a cylinder and a simulated engine leading edge verification test.

Dechaumphai, Pramote↗
NASA NTRS1989

Coupled fluid-thermal-structural analysis for aerodynamically heated structures

A finite element approach for integrated fluid-thermal-structural analysis of aerodynamically heated structures is described. The approach is used to solve: (1) the Navier-Stokes equations for the aerodynamic flow, (2) the energy equation of the structure for the temperature response, and (3) the quasi-static equilibrium equations of the structure for the structural response. Coupling between the three disciplines is demonstrated using two applications: Mach 7 flow over a panel and Mach 8 shock-shock interference on a cylinder.

Dechaumphai, Pramote↗
DOE OSTI2026

Cosmological constraints from the DESI DR1 joint power spectrum and bispectrum analysis

We derive cosmological parameter constraints from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) galaxy clustering data, based on a joint full-shape analysis of the power spectrum multipoles and the bispectrum monopole using the ShapeFit framework. This is the follow-up of our previous work, in which we obtained for the first time constraints on the ShapeFit parameters using the bispectrum of DESI DR1. Here we present the first ShapeFit cosmological inference results using the bispectrum of DESI DR1. We recover values for the matter density parameter and Hubble constant of respectively $Ω_m=0.310\pm0.012$ and $H_0=[68.92\pm0.97]\,\mathrm{km\, s^{-1} Mpc^{-1}}$, consistent with previous results from the full DESI DR1 dataset that did not use the bispectrum signal. The inclusion of the bispectrum significantly tightens the constraints on the amplitude of fluctuations, reducing the error-bars in $\ln(A_s\times10^{10})$ by approximately 20%, compared to using the power spectrum alone. We also explore extended cosmological models by performing fits for the evolving dark energy equation of state $w_0w_a$, and the sum of neutrino masses $\sum m_ν$. In these cases, we obtain constraints slightly larger than the ones from previous works from the DESI collaboration, due to not combining the full-shape results with other probes in all tracers. We find no strong evidence of deviations from standard $Λ$CDM, with the dark energy equation-of-state remaining within 2$σ$ from a cosmological constant $Λ$, and the neutrino mass being consistent with the normal hierarchy, $\sum m_ν<0.1\,[eV]$ at 95% confidence limit. These constraints are broadly consistent with other DESI DR1 analyses, thus validating the robustness of the ShapeFit compression approach and the inclusion of the bispectrum for cosmological inference.

Novell-Masot, S. [ICC, Barcelona U.; Geneva U., De↗
DOE OSTI2026

Cosmological constraints from the DESI DR1 joint power spectrum and bispectrum analysis

We derive cosmological parameter constraints from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) galaxy clustering data, based on a joint full-shape analysis of the power spectrum multipoles and the bispectrum monopole using the ShapeFit framework. This is the follow-up of our previous work, in which we obtained for the first time constraints on the ShapeFit parameters using the bispectrum of DESI DR1. Here we present the first ShapeFit cosmological inference results using the bispectrum of DESI DR1. We recover values for the matter density parameter and Hubble constant of respectively $Ω_m=0.310\pm0.012$ and $H_0=[68.92\pm0.97]\,\mathrm{km\, s^{-1} Mpc^{-1}}$, consistent with previous results from the full DESI DR1 dataset that did not use the bispectrum signal. The inclusion of the bispectrum significantly tightens the constraints on the amplitude of fluctuations, reducing the error-bars in $\ln(A_s\times10^{10})$ by approximately 20%, compared to using the power spectrum alone. We also explore extended cosmological models by performing fits for the evolving dark energy equation of state $w_0w_a$, and the sum of neutrino masses $\sum m_ν$. In these cases, we obtain constraints slightly larger than the ones from previous works from the DESI collaboration, due to not combining the full-shape results with other probes in all tracers. We find no strong evidence of deviations from standard $Λ$CDM, with the dark energy equation-of-state remaining within 2$σ$ from a cosmological constant $Λ$, and the neutrino mass being consistent with the normal hierarchy, $\sum m_ν<0.1\,[eV]$ at 95% confidence limit. These constraints are broadly consistent with other DESI DR1 analyses, thus validating the robustness of the ShapeFit compression approach and the inclusion of the bispectrum for cosmological inference.

Novell-Masot, S. [ICC, Barcelona U.; Geneva U., De↗