Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical homogenization”

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 271 records · Page 15

A modified restricted Euler equation for turbulent flows with mean velocity gradients

The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This modele is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence.

Girimaji, Sharath S.↗

Plasma instabilities of a finite-radius electron beam in a uniform plasma

To study electromagnetic radiation induced by electronn beam injection from the Space Shuttle, the electromagnetic dispersion equation or a finite-radius cold electron beam in a neutralizing background was solved numerically. The numerical solutions indicate that a keV electron beam can drive the beam and whistler modes unstable, regardless of whether the beam is homogeneous or has a finite radius. The results obtained were applied to explain the whistler waves radiated from the keV electron beam injected from Spacelab 2.

Wong, H. K.↗

Physics Based Model for Online Fault Detection in Autonomous Cryogenic Loading System

We report the progress in the development of the chilldown model for rapid cryogenic loading system developed at KSC. The nontrivial characteristic feature of the analyzed chilldown regime is its active control by dump valves. The two-phase flow model of the chilldown is approximated as one-dimensional homogeneous fluid flow with no slip condition for the interphase velocity. The model is built using commercial SINDAFLUINT software. The results of numerical predictions are in good agreement with the experimental time traces. The obtained results pave the way to the application of the SINDAFLUINT model as a verification tool for the design and algorithm development required for autonomous loading operation.

Cryogenics↗

Double-Diffusive Convection During Growth of Halides and Selenides

Heavy metal halides and selenides have unique properties which make them excellent materials for chemical, biological and radiological sensors. Recently it has been shown that selenohalides are even better materials than halides or selenides for gamma-ray detection. These materials also meet the strong needs of a wide band imaging technology to cover ultra-violet (UV), midwave infrared wavelength (MWIR) to very long wavelength infrared (VLWIR) region for hyperspectral imager components such as etalon filters and acousto-optic tunable filters (AO). In fact AOTF based imagers based on these materials have some superiority than imagers based on liquid crystals, FTIR, Fabry-Perot, grating, etalon, electro-optic modulation, piezoelectric and several other concepts. For example, broadband spectral and imagers have problems of processing large amount of information during real-time observation. Acousto-Optic Tunable Filter (AOTF) imagers are being developed to fill the need of reducing processing time of data, low cost operation and key to achieving the goal of covering long-wave infrared (LWIR). At the present time spectral imaging systems are based on the use of diffraction gratings are typically used in a pushbroom or whiskbroom mode. They are mostly used in systems and acquire large amounts of hyperspectral data that is processed off-line later. In contrast, acousto-optic tunable filter spectral imagers require very little image processing, providing new strategies for object recognition and tracking. They are ideally suited for tactical situations requiring immediate real-time image processing. But the performance of these imagers depends on the quality and homogeneity of acousto-optic materials. In addition for many systems requirements are so demanding that crystals up to sizes of 10 cm length are desired. We have studied several selenides and halide crystals for laser and AO imagers for MWIR and LWIR wavelength regions. We have grown and fabricated crystals of several materials such as mercurous chloride, mercurous bromide, mercurous iodide, lead chloride lead bromide, lead iodide, thallium arsenic selenide, gallium selenide, zince sulfide zinc selenide and several crystals into devices. We have used both Bridgman and physical vapor transport (PVT) crystal growth methods. In the past have examined PVT growth numerically for conditions where the boundary of the enclosure is subjected to a nonlinear thermal profile. Since past few months we have been working on binary and ternary materials such as selenoiodides, doped zinc sulfides and mercurous chloro bromide and mercurous bromoiodides. In the doped and ternary materials thermal and solutal convection play extremely important role during the growth. Very commonly striations and banding is observed. Our experiments have indicated that even in highly purified source materials, homogeneity in 1-g environment is very difficult. Some of our previous numerical studies have indicated that gravity level less than 10-4 (-g) helps in controlling the thermosolutal convection. We will discuss the ground based growth results of HgClxBr(1-x) and ZnSe growth results for the mm thick to large cm size crystals. These results will be compared with our microgravity experiments performed with this class of materials. For both HgCl-HgBr and ZnS-ZnSe the lattice parameters of the mixtures obey Vagard's law in the studied composition range. The study demonstrates that properties are very anisotropic with crystal orientation, and performance achievement requires extremely careful fabrication to utilize highest figure of merit. In addition, some parameters such as crystal growth fabrication, processing time, resolution, field of view and efficiency will be described based on novel solid solution materials. It was predicted that very similar to the pure compounds solid solutions also have very large anisotropy, and very precise oriented and homogeneous bulk and thin film crystals is required to achieve maximum performance of laser or imagers. Some of the parameters controlling the homogeneity such as thermos-solutal convection driven forces can be controlled in microgravity environments to utilize the benefits of these unique materials.

Singh, N. B.↗

Spacecraft detumbling using movable telescoping appendages

The dynamics of detumbling a randomly spinning spacecraft using externally mounted, movable telescoping appendages are studied both analytically and numerically. Two types of telescoping appendages are considered: (1) where an end mass is mounted at the end of an (assumed) massless boom; and (2) where the appendage is assumed to consist of a uniformly distributed homogeneous mass throughout its length. From an application of Liapunov's second method boom extension maneuvers can be determined to approach either of two desired final states: close to a zero inertial angular velocity state, and a final spin rate about only one of the principal axes. Recovery dynamics are evaluated analytically for the case of symmetrical deployment. Numerical examination of other asymmetrical cases verifies the practicality of using movable appendages to recover a randomly tumbling spacecraft.

Bainum, P. M.↗

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗

The parallel diffusion of cosmic rays in a random magnetic field.

Within the quasi-linear approximation, the existence of the parallel diffusion coefficient for cosmic rays in a random magnetic field (homogeneous, isotropic), despite the slow decay of the interaction between particles and random field, is demonstrated. As an example, the results of a numerical calculation of the parallel diffusion coefficient for a Gaussian random-field correlation function are presented. The numerical results are corroborated by asymptotic analysis and are compared to those of other theories.

Klimas, A.↗

Higher level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and various approaches to this simulation are compared. The subgrid model required by large eddy simulation is discussed as well as the use of this type of simulation in the development of models for the Reynolds-averaged equations and the application of direct simulation to the testing of both subgrid scale and Reynolds-averaged models. Numerical methods used in large eddy and direct simulation are described with emphasis on special purpose methods. Topics covered include the simulation of homogeneous flow, free shear flows, the mixing layer, wakes, and wall-bounded flows including channel flow and the boundary layer. Applications of large eddy simulation in the laboratory as in meteorological and other environmental flows are examined. Directions in which the work is proceeding and what can be expected from higher levels simulated on are examined.

Ferziger, J. H.↗

Optimal airframe synthesis for gust loads

An optimization capability for sizing airframe structures that are subjected to a combination of deterministic and random flight loads was established. The random vibration environment introduces the need for selecting a statistical process that best describes the random loads and permits computation of the dynamic response parameters of interest. Furthermore, it requires a formulation of design constraints that would minimize the conservativeness in the design and retain computational viability. The random loads are treated as a stationary, homogeneous process with a Gaussian probability distribution. The formulation of the analysis problem, the structure of the optimization programming system and a representative numerical example are discussed.

Hajela, P.↗

Hierarchic models for laminated plates

The research conducted in the formulation of hierarchic models for laminated plates is described. The work is an extension of the work done for laminated strips. The use of a single parameter, beta, is investigated that represents the degree to which the equilibrium equations of three dimensional elasticity are satisfied. The powers of beta identify members in the hierarchic sequence. Numerical examples that were analyzed with the proposed sequence of models are included. The results obtained for square plates with uniform loading and homogeneous boundary conditions are very encouraging. Several cross-ply and angle-ply laminates were evaluated and the results compared with those of the fully three dimensional model, computed using MSC/PROBE, and with previously reported work on laminated strips.

Szabo, Barna A.↗

Passive Microwave Measurements Over Conifer Forests at L-Band and C-Band

Measurements have been made at L-band and C-band over conifer forests in Virginia to study the response of passive microwave instruments to biomass and soil moisture. A series of aircraft measurements were made in July, August and November, 1999 over relatively homogenous conifer forests of varying biomass. Three radiometers participated in these measurements. These were: 1) the L-band radiometer ESTAR, a horizontally polarized synthetic aperture radiometer which has been used extensively in past measurements of soil moisture; 2) the L-band radiometer SLFMR, a vertically polarized cross-track scanner which has been used successfully in the past for mapping sea surface salinity; and 3) The ACMR, a new C-band radiometer which operates at V- and H-polarization and in the configuration for these experiments did not scan. All three radiometers were flown on the NASA P-3 aircraft based at the Goddard Space Flight Center's Wallops Flight Facility. The ESTAR and SLFMR were mounted in the bomb bay of the P-3 and imaged across track whereas the ACMR was mounted to look aft at 54 degrees up from nadir. Data was collected at altitudes of 915 meters and 457 meters. The forests consisted of relatively homogeneous "managed" stands of conifer located near Waverly, Virginia. This is a relatively flat area about 30 miles southeast of Richmond, VA with numerous stands of trees being grown for the forestry industry. The stands selected for study consisted of areas of regrowth and mature stands of pine. In addition, a small stand of very large trees was observed. Soil moisture sampling was done in each stand during the aircraft over flights. Data was collected on July 7, August 27, November 15 and November 30, 1999. Measurements were made with ESTAR on all days. The ACMR flew on the summer missions and the SLFMR was present only on the August 27 flight. Soil moisture varied from quite dry on July 7 to quite moist on November 30 (which was shortly after a period of rain). The microwave images clearly distinguish between the different forest stands. Research is continuing to seek a quantitative correlation with biomass and surface soil moisture.

LeVine, D. M.↗

Modeling and Validation of Damped Plexiglas Windows for Noise Control

Windows are a significant path for structure-borne and air-borne noise transmission in general aviation aircraft. In this paper, numerical and experimental results are used to evaluate damped plexiglas windows for the reduction of structure-borne and air-borne noise transmitted into the interior of an aircraft. In contrast to conventional homogeneous windows, the damped plexiglas windows were fabricated using two or three layers of plexiglas with transparent viscoelastic damping material sandwiched between the layers. Transmission loss and radiated sound power measurements were used to compare different layups of the damped plexiglas windows with uniform windows of the same nominal thickness. This vibro-acoustic test data was also used for the verification and validation of finite element and boundary element models of the damped plexiglas windows. Numerical models are presented for the prediction of radiated sound power for a point force excitation and transmission loss for diffuse acoustic excitation. Radiated sound power and transmission loss predictions are in good agreement with experimental data. Once validated, the numerical models were used to perform a parametric study to determine the optimum configuration of the damped plexiglas windows for reducing the radiated sound power for a point force excitation.

Buehrle, Ralph D.↗

Coherent Eigenmodes in Homogeneous MHD Turbulence

The statistical mechanics of Fourier models of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence is discussed, along with their relevance for dissipative magnetofluids. Although statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation, i.e., we have coherent structure. We use eigenanalysis of the modal covariance matrices in the probability density function to explain this phenomena in terms of `broken ergodicity', which is defined to occur when dynamical behavior does not match ensemble predictions on very long time-scales. We provide examples from 2-D and 3-D magnetohydrodynamic simulations of homogeneous turbulence, and show new results from long-time simulations of MHD turbulence with and without a mean magnetic field

Shebalin, John V.↗

Electromagnetic electron beam instabilities - Hot, isotropic beams

This paper considers the linear theory of electromagnetic instabilities driven by an electron beam in a homogeneous, nonrelativistic, Vlasov plasma. The beam is relatively hot, isotropic in its own frame, and streams parallel or antiparallel to a magnetic field B. Numerical solutions of the full dispersion equation for propagation parallel or antiparallel to B are presented, and the linear properties of the whistler heat flux and electron beam firehose instabilities are exhibited and compared. Under a broad range of parmameters the former mode has the lower beam speed threshold, and the larger maximum growth rate. In addition, it is demonstrated that, for a sufficiently large relative beam density, relative beam temperature, and plasma beta the whistler heat flux instability has a much lower beam speed threshold than the electrostatic electron beam instability. The application of these instabilities to first-order Fermi acceleration of electrons at space plasma shocks is discussed.

Gary, S. P.↗

Direct numerical study of leading-edge contamination

Instability, turbulence, and relaminarization in the attachment-line region of swept and unswept cylindrical bodies are studied by numerical solution of the full Navier-Stokes equations. The flow is simulated over a strip containing the attachment-line and treated as homogeneous in the spanwise direction; the disturbances decay exponentially upstream. Transpiration through the wall may be prescribed. The new method, which admits completely general disturbance, agrees with published linear-stability results, which were limited to an apparently restrictive form of disturbance. Fully developed turbulent solutions with sweep are generated and compare well with the experiment. The turbulence is subcritical (except for blowing), resulting in large hysteresis loops. By lowering the sweep Reynolds number, or increasing the suction, the turbulent flow is made to relaminarize. The relaminarization Reynolds number is much less sensitive to suction than the linear-stability Reynolds number. Extensive attempts to detect the postulated nonlinear instability of the unswept flow failed, suggesting that this flow is linearly and nonlinearly stable.

Spalart, P. R.↗

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗

Adhesion of Mineral and Soot Aerosols can Strongly Affect their Scattering and Absorption Properties

We use the numerically exact superposition T-matrix method to compute the optical cross sections and the Stokes scattering matrix for polydisperse mineral aerosols (modeled as homogeneous spheres) covered with a large number of much smaller soot particles. These results are compared with the Lorenz-Mie results for a uniform external mixture of mineral and soot aerosols. We show that the effect of soot particles adhering to large mineral particles can be to change the extinction and scattering cross sections and the asymmetry parameter quite substantially. The effect on the phase function and degree of linear polarization can be equally significant.

Mishchenko, Michael I.↗

General chemical kinetics computer program for static and flow reactions, with application to combustion and shock-tube kinetics

A general chemical kinetics program is described for complex, homogeneous ideal-gas reactions in any chemical system. Its main features are flexibility and convenience in treating many different reaction conditions. The program solves numerically the differential equations describing complex reaction in either a static system or one-dimensional inviscid flow. Applications include ignition and combustion, shock wave reactions, and general reactions in a flowing or static system. An implicit numerical solution method is used which works efficiently for the extreme conditions of a very slow or a very fast reaction. The theory is described, and the computer program and users' manual are included.

Bittker, D. A.↗