SEARCH · Engineering Papers
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.
The relationship between length of day (LOD) variations and the Southern Oscillation (SO), and with the associated El Nino phenomenon, is studied. Interannual changes in the LOD are shown to be linearly related to a SO sea level pressure index. A possible physical mechanism for this relationship is the thermal winds caused by changes in the equator to pole heat transports associated with El Nino events and the SO.
To date modeling of crystal growth of optoelectronic materials using Physical Vapor Transport has been limited to the study of the fluid phase. To achieve it, the equations of coupled heat, mass and momentum transfer in the gas have to be solved. The first objective of this study is to examine the effect of heat conduction in the crystal on the fluid flow in the neighborhood of the interface. Heat transfer boundary conditions on both interfaces were modified to take into account the additional heat flux between gas and solid. It is proved that heat conduction does not affect the fluid flow. In the presence of gravity, density gradients in the fluid phase generate convection responsible for the problem of a nonplanar growth of the interface. The second objective is to study systematically under one-g the different possible flows in order to solve this problem. Depending on the parameters, a diffusive mode and three convective modes (thermal, solutal and thermo-solutal) are observed. The competition between thermal and solutal convections leads to a mathematical condition which can be used to achieve a planar growth. It is proven that, under the physical conditions chosen, this mathematical condition cannot be thermodynamically satisfied.
Techniques for extracting boundary layer parameters from measurements of a short pulse CO2 Doppler Lidar are described. The radial velocity measurements have a range resolution of 150 m. With a pulse repetition rate of 20 Hz, it is possible to perform scannings in two perpendicular vertical planes in approx. 72 s. By continuously operating the Lidar for about an hour, one can extract stable statistics of the radial velocities. Assuming that the turbulence is horizontally homogeneous, the mean wind, its standard deviations, and the momentum fluxes were estimated. From the vertically pointing beam, the first, second, and third moments of the vertical velocity were also estimated. Spectral analysis of the radial velocities is also performed from which, by examining the amplitude of the power spectrum at the inertial range, the kinetic energy dissipation was deduced. Finally, using the statistical form of the Navier-Stokes equations, the surface heat flux is derived as the residual balance between the vertical gradient of the third moment of the vertical velocity and the kinetic energy dissipation.
Techniques for extraction boundary layer parameters from measurements of a short-pulse CO2 Doppler lidar are described. The measurements are those collected during the First International Satellites Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE). By continuously operating the lidar for about an hour, stable statistics of the radial velocities can be extracted. Assuming that the turbulence is horizontally homogeneous, the mean wind, its standard deviations, and the momentum fluxes were estimated. Spectral analysis of the radial velocities is also performed from which, by examining the amplitude of the power spectrum at the inertial range, the kinetic energy dissipation was deduced. Finally, using the statistical form of the Navier-Stokes equations, the surface heat flux is derived as the residual balance between the vertical gradient of the third moment of the vertical velocity and the kinetic energy dissipation. Combining many measurements would normally reduce the error provided that, it is unbiased and uncorrelated. The nature of some of the algorithms however, is such that, biased and correlated errors may be generated even though the raw measurements are not. Data processing procedures were developed that eliminate bias and minimize error correlation. Once bias and error correlations are accounted for, the large sample size is shown to reduce the errors substantially. The principal features of the derived turbulence statistics for two case studied are presented.
The method of Bozzolo, Ferrante and Smith (BFS) is applied for the calculation of the heat of segregation of single substitutional impurities in fcc metals. A simple equation for predicting the heat of segregation is derived for the rigid case (no atomic relaxations). The results of including atomic relaxation using a Monte Carlo method are also presented and the results compared with a number of experimental and theoretical results.
A Reynolds stress turbulence model has been implemented in the COMMIX code, together with transport equations describing turbulent heat fluxes, variance of temperature fluctuations, and dissipation of turbulence kinetic energy. The model has been verified partially by simulating homogeneous turbulent shear flow, and stable and unstable stratified shear flows with strong buoyancy-suppressing or enhancing turbulence. This article outlines the model, explains the verifications performed thus far, and discusses potential applications of the COMMIX-RSM code in several domains, including, but not limited to, analysis of thermal striping in engineering systems, simulation of turbulence in combustors, and predictions of bubbly and particulate flows.
Accurate numerical simulation of chemical vapor deposition (CVD) processes requires a general purpose computational fluid dynamics package combined with specialized capabilities for high temperature chemistry. In this report, we describe the implementation of these specialized capabilities in the spectral element code NEKTON. The thermal expansion of the gases involved is shown to be accurately approximated by the low Mach number perturbation expansion of the incompressible Navier-Stokes equations. The radiative heat transfer between multiple interacting radiating surfaces is shown to be tractable using the method of Gebhart. The disparate rates of reaction and diffusion in CVD processes are calculated via a point-implicit time integration scheme. We demonstrate the use above capabilities on prototypical CVD applications.
A model is presented to compute the turbulent kinetic energy dissipation length scale l(sub epsilon) in a stably stratified shear flow. The expression for l(sub epsilon) is derived from solving the spectral balance equation for the turbulent kinetic energy. The buoyancy spectrum entering such equation is constructed using a Lagrangian timescale with modifications due to stratification. The final result for l(sub epsilon) is given in algebraic form as a function of the Froude number Fr and the flux Richardson number R(sub f), l(sub epsilon) = l(sub epsilon)(Fr, R(sub f). The model predicts that for R(sub f) less than R(sub fc), l(sub epsilon) decreases with stratification. An attractive feature of the present model is that it encompasses, as special cases, some seemingly different models for l(sub epsilon) that have been proposed in the past by Deardorff, Hunt et al., Weinstock, and Canuto and Minotti. An alternative form for the dissipation rate epsilon is also discussed that may be useful when one uses a prognostic equation for the heat flux. The present model is applicable to subgrid-scale models, which are needed in large eddy simulations (LES), as well as to ensemble average models. The model is applied to predict the variation of l(sub epsilon) with height z in the planetary boundary layer. The resulting l(sub epsilon) versus z profile reproduces very closely the nonmonotonic profile of l(sub epsilon) exhibited by many LES calculations, beginning with the one by Deardorff in 1974.
Almost one hundred years ago Korteweg published a theory of how stresses could be induced in miscible fluids by concentration gradients, causing phenomena that would appear to be the same as with immiscible fluids. Miscible fluids could manifest a transient or effective interfacial tension (EIT). To this day, there has been no definitive experiment to confirm Korteweg's model but numerous fascinating and suggestive experiments have been reported. The goal of TIPMPS is to answer the question: Can concentration and temperature gradients in miscible materials induce stresses that cause convection? Many polymer processes involving miscible monomer and polymer systems could be affected by fluid flow and so this work could help understand miscible polymer processing, not only in microgravity, but also on earth. Demonstrating the existence of this phenomenon in miscible fluids will open up a new area of study for materials science. The science objectives of TIPMPS are: (1) Determine if convection can be induced by variation of the width of a miscible interface; (2) Determine if convection can be induced by variation of temperature along a miscible interface; (3) Determine if convection can be induced by variation of conversion along a miscible interface An interface between two miscible fluids can best be created via a spatially-selective photopolymerization of dodecyl acrylate with a photoinitiator, which allows the creation of precise and accurate concentration gradients between polymer and monomer. Optical techniques will be used to measure the refractive index variation caused by the resultant temperature and concentration fields. The viscosity of the polymer will be measured from the increase in the fluorescence of pyrene. Because the large concentration and temperature gradients cause buoyancy-driven convection that prevents the observation of the predicted flows, the experiment must be done in microgravity. In this report, we will consider our efforts to estimate the square gradient parameter, k, and our use of the estimates in modeling of the planned TIPMPS experiments. We developed a model consisting of the heat and diffusion equations with convective terms and of the Navier-Stokes equations with an additional volume force written in the form of the Korteweg stresses arising from nonlocal interaction in the fluid. The fluid's viscosity dependence on polymer conversion and temperature was taken from measurements of poly(dodecyl acrylate). Numerical modeling demonstrated that significant flows would arise for conditions corresponding to the planned experiments.
This paper describes recent improvements to the LAURA and HARA codes. LAURA is a CFD code for aerothermodynamics, and HARA evaluates the shock-layer radiation that provides the radiative source term for the flowfield energy equations and radiative heating to a surface. The next release of LAURA and HARA includes a variety of new capabilities. These new capabilities include an automated uncertainty quantification workflow for radiative heat transfer, options for specifying surface roughness and turbulent transition location in the algebraic turbulence models, and improved grid and solution interpolation techniques. Additionally, the computational efficiency of both LAURA and HARA have been improved. Optimization of the MPI communication routines in LAURA are shown to improve the parallel efficiency of the primary flow when running with multiple processes per block, and recent optimization of HARA leverage graphics processing unit (GPU) acceleration in the radiation calculations. Using GPU acceleration of HARA is shown to decrease the cost of the radiation line-of-sight calculation by approximately one order of magnitude for a 10.5 km/s Earth entry simulation.
Mathematical analysis of shock wave structure incorporating Maxwell moment methods and Boltzmann equation for stress and heat flux constitutive relations
This paper describes a method of evaluating heat transfer coefficient from a velocity distribution in the boundary layer. The power law velocity profile and universal velocity distribution in a smooth pipe have been used to evaluate the wall shear stress, and the heat transfer coefficient was evaluated from the wall shear stress using the Chilton-Colburn analogy. The predicted pressure drop was compared with the pressure drop computed by the Colebrook equation, and the predicted heat transfer coefficient was compared with the heat transfer coefficient calculated by the Dittus-Boelter equation. The accuracy of pressure drop prediction was within 1%, and the accuracy of heat transfer coefficient prediction was within 8%. This method has the potential to use a Navier-Stokes based CFD solution to evaluate the heat transfer coefficient in nodal or network flow codes.
This paper describes a method of evaluating heat transfer coefficient from a velocity distribution in the boundary layer. The power law velocity profile and universal velocity distribution in a smooth pipe have been used to evaluate the wall shear stress, and the heat transfer coefficient was evaluated from the wall shear stress using the Chilton-Colburn analogy. The predicted pressure drop was compared with the pressure drop computed by the Colebrook equation, and the predicted heat transfer coefficient was compared with the heat transfer coefficient calculated by the Dittus-Boelter equation. The accuracy of pressure drop prediction was within 1%, and the accuracy of heat transfer coefficient prediction was within 8%. This method has the potential to use a Navier-Stokes based CFD solution to evaluate the heat transfer coefficient in nodal or network flow codes.
The linearized global primitive equation model of Salby and Garcia (1987) was modified to study tropospheric circulations and their interaction with tropical convection. The vertical resolution was increased to approximately 1.5 km, with the effective vertical domain extending the approximately 40 km (higher altitudes in the calculation are used to ensure the radiation upper boundary condition is satisfied). A formal wave conditional instability of the second kind (CISK) parameterization was introduced. That parameterization allows low-frequency disturbances to interact with convection by organizing surface moisture. Boundary layer friction has two important effects in these calculations: it exaggerates surface convergence near the equator by driving the flow out of geostrophic balance; and it introduces a phase shift between the convergence pattern and the circulation, which drives the heating out of quadrature with the temperature field and allows a positive feedback to take place. The observational complement of this theoretical work involves Global Cloud Imagery (GCI) composited from six satellites in the International Satellite Cloud Climatology Project (ISCCP) which were simultaneously observing the Earth's cloud field. Having global coverage, horizontal resolution of half a degree, and produced synoptically every 3 hours, the GCI affords an unprecendented view of global convection. A full year (1983 to 1984) of GCI has now been successfully created. Among other phenomena, that imagery reveals a regular progression of easterly waves in the convective pattern of the Western Pacific. Those disturbances have odd symmetry about the equator and organize convection into a pattern with similar symmetry. In time-mean maps of cloud cover, these propagating disturbances produce a split Inter-Tropical Convergence Zone (ITCZ) in the western Pacific.
While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.
Interactions at the interface between atmospheric pressure plasmas and liquids are being investigated to address applications ranging from nanoparticle synthesis to decontamination and fertilizer production. Furthermore, many of these applications involve activation of droplets wherein the droplet is fully immersed in the plasma and synergistically interacts with the plasma. To better understand these interactions, two-dimensional modeling of radio frequency (RF) glow discharges at atmospheric pressure operated in He with an embedded lossy dielectric droplet (tens of microns in size) was performed. The properties of the sheath that forms around the droplet were investigated over the RF cycle. The electric field in the bulk plasma polarizes the dielectric droplet while the electron drift in the external electric field is shadowed by the droplet. The interaction between the bulk and sheath electric fields produces a maximum in E/N (electric field/gas number density) at the equator on one side of the droplet where the bulk and sheath fields are aligned in the same direction and a minimum along the opposite equator. Due to resistive heating, the electron temperature T e is maximum 45° above and below the equator of the droplet where power deposition per electron is the highest. Although the droplet is, on the average, negatively charged, the charge density on the droplet is positive on the poles and negative on the equator, as the electron motion is primarily due to diffusion at the poles but due to drift at the equator.
Temperature distributions and heat generation in viscoelastic isotropic solid resulting from mechanical deformations determined by energy equation
Closed form expressions for inelastic state-to-state and state-specific dissociative rate coefficients for utilization in vibrational master equation studies of shock heated CO, N2, and O2 highly dilute in Ar are considered. The master equation is linearized by neglecting diatom-diatom collisions and recombination. Master equation results indicate that the most significant contribution to dissociation comes from low and mid lying vibrational levels.