Nonlinear combustion instability - Computer experiments using codes COMB and TRDL
Mathematical models for computer calculations of two dimensional hydrodynamic nonlinear combustion instability flow
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.
Mathematical models for computer calculations of two dimensional hydrodynamic nonlinear combustion instability flow
An extensive parametric investigation was performed of the extraction of energy in output gaps of klystron amplifiers, using 3-D computer programs. Due to complexity of the program which used a hydrodynamic, axially and radially deformable disk ring model and the resulting long computing time, the investigation was limited to the output gap, by far the most important and difficult part of the klystron interaction. Results show that, for a confined flow focused beam throughout the penultimate cavity, radial velocities remain very small and the beam is highly laminar. It was, therefore, concluded that possible errors resulting from treating only the output cavity in 3-D would remain small.
A model is proposed in which the strong hydrodynamic shock forming around the Space Shuttle, and occasionally detaching, generates strong hydrodynamic turbulence downstream. The level of this turbulence is found to be sufficiently great to account for the heating of electrons that has been observed. A shock formed in front of the Orbiter is shown to satisfy the Rankine-Hugoniot momentum and pressure balance equations. The hydrodynamic turbulence is coupled through ion-neutral collisions to the O2(+) and N2(+) ions of outgas.
The process of gravitational fragmentation of molecular clouds is investigated through a comparison of numerical and observational results. The expected millimeter wave molecular line emission from a model fragmenting cloud generated by a numerical hydrodynamic simulation is calculated and compared with observations of HCO(+) from the star-forming region W49A. This investigation suggests that the rotating ring of H II regions, the necklace in W49A, may have formed in the gravitational fragmentation of a flattened rotating molecular cloud.
This annual report summarizes the work accomplished on rewetting of monogroove heat pipe in space station. Specifically, theoretical and experimental investigations of the rewetting characteristics of thin liquid films over unheated and heated capillary grooved plates were performed. To investigate the effect of gravity on rewetting, the grooved surface was placed in upward and downward facing positions. Profound gravitational effects were observed as the rewetting velocity was found to be higher in the upward than in the downward facing orientation. The difference was even greater with higher initial plate temperatures. With either orientation, it was found that the rewetting velocity increased with the initial plate temperature. But when the temperature was raised above a rewetting temperature, the rewetting velocity decreased with the initial plate temperature. Hydrodynamically controlled and heat conduction controlled rewetting models were then presented to explain and to predict the rewetting characteristics in these two distinct regions. The predicted rewetting velocities were found to be in good agreement with experimental data with elevated plate temperatures.
The Greenland, Iceland, and Norwegian Seas provide the only link between the North Atlantic and the Arctic Oceans. This is a very dynamic region, rich in mesoscale variability. A two-layer, hydrodynamic version of the Navy layered ocean model is used to simulate the mesoscale frontal features and associated current systems. The model is wind-driven using monthly mean wind stresses and inflow/outflow mass flux from the Arctic and North Atlantic Oceans. The current transports from the model results compare favorably with transport measurements from a number of observational experiments.
A transient moving contact line is investigated experimentally. The dynamic interface shape between 20 and 800 microns from the contact line is compared with theory. A novel experiment is devised, in which the contact line is set into motion by electrically altering the solid-liquid surface tension gamma(sub SL). The contact line motion simulates that of spontaneous wetting along a vertical plate with a maximum capillary number Ca approx. = 4 x 10(exp -2). The images of the dynamic meniscus are analyzed as a funtion of Ca. For comparison, the steady-state hydrodynamic equation based on the creeping flow model in a wedge geometry and the three-region uniform perturbation expansion of Cox (1986) is adopted. The interface shape is well depicted by the uniform solutions for Ca <= 10(exp -3). However, for Ca > 10(exp -3), the uniform solution over-predicts the viscous bending. This over-prediction can be accounted for by modifying the slip coefficient within the intermediate solution. With this correction, the measured interface shape is seen to match the theoretical prediction for all capillary numbers. The amount of slip needed to fit the measurements does not scale with the capillary number.
The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.
The Interface Region Imaging Spectrograph has routinely observed the flaring Mg II near-ultraviolet (NUV) spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward-modeled both the general properties of these lines and specific flare observations. Generally these have forward-modeled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward-modeling Mg II in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II that (1) partial frequency distribution (PRD) is still required during flare simulations despite the increased densities; (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot; (3) including Mg I in non-LTE (NLTE) results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum; (4) only hydrogen and Mg II need to be included in NLTE; (5) ideally the nonequilibrium hydrogen populations, with nonthermal collisional rates, should be used rather than the statistical equilibrium populations; (6) an atom consisting of only the ground state, h and k upper levels, and continuum level is insufficient to model the resonance lines; and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of the transition region and coronal irradiation in flares.
An investigation of a 1/14-scale dynamically similar model of a panto-base version of the Chase C-123 airplane was conducted to evaluate the hydrodynamic characteristics of the airplane. The resistance, longitudinal stability, and spray patterns during take-off and general behavior in calm- and rough-water landings were determined. Brief calm-water tests were made to compare the initial vertical impact accelerations of the model with and without hydro-skis. Take-off stability was satisfactory for calm-water operation. A ratio of gross load to maximum resistance of 3,6 was obtained. Heavy spray reached the propellers only during ski emergence. The landing behavior in calm water and in waves 3 feet by 150 feet (full scale) was satisfactory for a normal range of trim angles. Initial impacts in calmwater landings resulted in vertical accelerations of about 2 1/2 with the hydro-skis installed and about 4g with the hydro-skis removed,
I will present models of the solar atmospheric response to flare heating. The models solve the equations of non-LTE radiation hydrodynamics with an electron beam added as a flare energy source term. Radiative transfer is solved in detail for many important optically thick hydrogen and helium transitions and numerous optically thin EUV lines making the models ideally suited to study the emission that is produced during flares. I will pay special attention to understanding key EUV lines as well the mechanism for white light production. I will also present preliminary results of how the model solar atmosphere responds to Fletcher & Hudson type flare heating. I will compare this with the results from flare simulations using the standard thick target model.
An efficient, flexible, three-dimensional, hydrodynamic, computer code has been developed for a spherical cap geometry. The code will be used to simulate NASA's Atmospheric General Circulation Experiment (AGCE). The AGCE is a spherical, baroclinic experiment which will model the large-scale dynamics of our atmosphere; it has been proposed to NASA for future Spacelab flights. In the AGCE a radial dielectric body force will simulate gravity, with hot fluid tending to move outwards. In order that this force be dominant, the AGCE must be operated in a low gravity environment such as Spacelab. The full potential of the AGCE will only be realized by working in conjunction with an accurate computer model. Proposed experimental parameter settings will be checked first using model runs. Then actual experimental results will be compared with the model predictions. This interaction between experiment and theory will be very valuable in determining the nature of the AGCE flows and hence their relationship to analytical theories and actual atmospheric dynamics.
A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.
Following the breakup of the Space Shuttle Columbia during reentry a NASA/Contractor investigation team was formed to examine the probable damage inflicted on Orbiter Thermal Protection System elements by impact of External Tank insulating foam projectiles. The authors formed a working subgroup within the larger team to apply the Smooth Particle Hydrodynamics code SPHC to the damage estimation problem. Numerical models of the Orbiter's tiles and of the Tank's foam were constructed and used as inputs into the code. Material properties needed to properly model the tiles and foam were obtained from other working subgroups who performed tests on these items for this purpose. Two- and three-dimensional models of the tiles were constructed, including the glass outer layer, the main body of LI-900 insulation, the densified lower layer of LI-900, the Nomex felt mounting layer, and the Aluminum 2024 vehicle skin. A model for the BX-250 foam including porous compression, elastic rebound, and surface erosion was developed. Code results for the tile damage and foam behavior were extensively validated through comparison with Southwest Research Institute foam-on-tile impact experiments carried out in 1999. These tests involved small projectiles striking individual tiles and small tile arrays. Following code and model validation we simulated impacts of larger foam projectiles on the examples of tile systems used on the Orbiter. Results for impacts on the main landing gear door are presented in this paper, including effects of impacts at several angles, and of rapidly rotating projectiles. General results suggest that foam impacts on tiles at about 500 mph could cause appreciable damage if the impact angle is greater than about 20 degrees. Some variations of the foam properties, such as increased brittleness or increased density could increase damage in some cases. Rotation up to 17 rps failed to increase the damage for the two cases considered. This does not rule out other cases in which the rotational energy might lead to an increase in tile damage, but suggests that in most cases rotation will not be an important factor.
The reported investigation has the objective to provide a basis for understanding the hydrodynamic characteristics of the continuous flow electrophoresis chamber and their effects on the separation process. Particular emphasis is placed on the role buoyancy plays in establishing the basic flow and affecting its stability. The flow and temperature field is discussed, taking into account mathematical models for the structure of the temperature field, models for the structure of the axial velocity field, a two-dimensional constant thermal conductivity model, the axial velocity field, a two-dimensional constant transport properties model, and one-dimensional velocity fields. Questions of hydrodynamic stability are considered, giving attention to the stability of a fully developed flow with an axial temperature gradient.
A model of convective turbulence which takes radiative dissipation, rotation, and convective motion anisotropy into account, on the basis of a closure for the nonlinear interactions that employs the growth rates of hydrodynamic instabilities, is used to obtain a theoretical framework for modeling the primordial solar nebula. It is assumed that convection is the sole source of turbulence causing the solar nebula to evolve. Vertical structure equations in the thin disk approximation are developed and a detailed comparison with the previous solar nebula convective models of such workers as Lin et al. (1981, 1982) is undertaken. The present values for the turbulent efficiency are much lower and more sensitive to opacity and surface density, resulting in low turbulent speeds, a more massive disk, a lower accretion rate 'best value', and a longer characteristic dispersal time for the disk. It is concluded that convection may not be the dominant source of turbulence needed to evolve young solar/stellar nebulae.
A series of one and two dimensional hydrodynamic simulations of a ring wave in interstellar gas disks was completed. These calculations included nonlinear source terms to model the effects of interstellar interactions and star formation, as well as the spatial-temporal gas flow. Toomre's kinematical model was merged with the Arnold, Shandarin, and Zeldovich 'pancake' theory of caustics in galaxy formation. The resulting theory can describe almost all the structure in restricted three-body simulations of single-pass collisions, even with multi-component potentials. Off-center galactic collisions were studied to understand the dynamics involved. Multi-color optical and near-infrared observations of faint tidal features were performed in about two dozen interacting galaxies selected from the Arp atlas. This sample provided evidence for ongoing star formation in tidal structures, and even enhancements of star formation in some cases. The task of assembling the data for gas-rich, late-type galaxies, was undertaken to see if a more coherent picture of the gas distribution would emerge from the more complete data. Analytic solutions of the equations with subsonic flows to balance gas consumption for expulsion form a galactic fountain were also derived.
Tests of a PB2Y-3 flying boat were made at the U.S> Naval Air Station, Patuxent River, Md., to determine its hydrodynamic trim limits of stability. Corresponding tests were also made of a 1/8-size powered dynamic model of the same flying boat in Langley tank no. 1. During the tank tests, the full-size testing procedure was reproduced as closely as possible in order to obtain data for a direct correlation of the results. As a nominal gross load of 66,000 pounds, the lower trim limits of the full-size and model were in good agreement above a speed of 80 feet per second. As the speed decreased below 80 feet per second, the difference between the model trim limits and full-scale trim limits gradually became larger. The upper trim limit of the model with flaps deflected 0 deg was higher than that of the full-size, but the difference was small over the speed range compared. At flap deflections greater than 0 deg, it was not possible to trim either the model of the airplane to the upper limit with the center of gravity at 28 percent of the mean aerodynamic chord. The decrease in the lower trim limits with increase in flap deflection showed good agreement for the airplane and model. The lower trim limits obtained at different gross loads for the full-size airplane were reduced to approximately a single curve by plotting trim against the square root of C(sub delta (sub o)) divided by C(sub V).