Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Breakup”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Entry Debris Field Estimation Methods and Application to Compton Gamma Ray Observatory Disposal

For public safety reasons, the Compton Gamma Ray Observatory (CGRO) was intentionally deorbited on June 4, 2000. This deorbit was NASA's first intentional controlled deorbit of a satellite, and more will come including the eventual deorbit of the International Space Station. To maximize public safety, satellite deorbit planning requires conservative estimates of the debris footprint size and location. These estimates are needed to properly design a deorbit sequence that places the debris footprint over unpopulated areas, including protection for deorbit contingencies. This paper details a method for estimating the length (range), width (crossrange), and location of entry and breakup debris footprints. This method utilizes a three degree-of-freedom Monte Carlo simulation incorporating uncertainties in all aspects of the problem, including vehicle and environment uncertainties. The method incorporates a range of debris characteristics based on historical data in addition to any vehicle-specific debris catalog information. This paper describes the method in detail, and presents results of its application as used in planning the deorbit of the CGRO.

Mrozinski, Richard B.↗
NASA NTRS2001

Correlation Studies of Sea Ice Concentration with Surface Temperature and Meltponding

The spatial and temporal variability of sea ice concentrations derived from passive microwave data is studied in conjunction with co-registered high resolution infrared and visible satellite data. Cloud free infrared and visible data provide surface temperature and large scale surface characteristics, respectively, that can be used to better understand regional and seasonal fluctuations in ice concentrations. Results from correlation analysis of ice concentration versus surface temperature data show the intuitively expected negative relationship but the strength in the relationship is unexpectedly very strong. In the Antarctic, the correlation is consistently very high spatially when yearly anomalies are used, and not so high in some areas when seasonal anomalies are used, especially during spring and summer. In the monthly anomalies, the correlation is also good, especially in dynamically active regions. The expanse in the anomalies in surface temperature are shown to go way beyond the sea ice regions into the open ocean and continental areas, suggesting strong atmospheric forcing. Weak correlations are normally found in highly consolidated areas, where large changes in temperature do not cause large changes in ice concentration on a short term, and in open ocean polynya areas, where the change in ice concentration may be cause by melt from the underside of the ice. In the Arctic, strong correlations between surface temperature and ice concentration are evident for all seasons except during the summer. In the summer, factors such as meltponding, surface wetness, and ice breakup, as detected by high resolution visible data, contributes to larger uncertainties in the determination of ice concentration and the lack of good correlation of the variables.

Comiso, J. C.↗
NASA NTRS2001

Nonlinear Dynamics of a Diffusing Interface

Excitation of two miscible-viscous liquids inside a bounded enclosure in a microgravity environment has shown the evolution of quasi-stationary waves of various modes for a range of parameters. We examine computationally the nonlinear dynamics of the system as the interface breakup and bifurcates to resonance structures typified by the Rayleigh-Taylor instability mechanism. Results show that when the mean steady field is much smaller than the amplitude of the sinusoidal excitation, the system behaves linearly, and growth of quasi-stationary waves occurs through the Kelvin-Helmholtz instability mechanism. However, as the amplitude of excitation increases, nonlinearity occurs through subharmonic bifurcation prior to broadband chaos.

Duval, Walter M. B.↗
NASA NTRS2001

Grid Convergence of High Order Methods for Multiscale Complex Unsteady Viscous Compressible Flows

Grid convergence of several high order methods for the computation of rapidly developing complex unsteady viscous compressible flows with a wide range of physical scales is studied. The recently developed adaptive numerical dissipation control high order methods referred to as the ACM and wavelet filter schemes are compared with a fifth-order weighted ENO (WENO) scheme. The two 2-D compressible full Navier-Stokes models considered do not possess known analytical and experimental data. Fine grid solutions from a standard second-order TVD scheme and a MUSCL scheme with limiters are used as reference solutions. The first model is a 2-D viscous analogue of a shock tube problem which involves complex shock/shear/boundary-layer interactions. The second model is a supersonic reactive flow concerning fuel breakup. The fuel mixing involves circular hydrogen bubbles in air interacting with a planar moving shock wave. Both models contain fine scale structures and are stiff in the sense that even though the unsteadiness of the flows are rapidly developing, extreme grid refinement and time step restrictions are needed to resolve all the flow scales as well as the chemical reaction scales.

Sjoegreen, B.↗
NASA NTRS2000

Comet Shoemaker-Levy 9 Dust Size and Velocity Distributions

Pre-impact observations of Comet Shoemaker-Levy 9 (S-L9) obtained with the Hubble Space Telescope are examined, and a model of an active, dust-producing comet is fitted to images of fragments G, H, K, and L. The model assumes steady isotropic dust emission from each fragment's sunlit hemisphere. Best-fit results indicate that the dominant light-scatterers in these fragments' comae were relatively large dust grains of radii 10 micrometers < R < 3 mm. The fragments' dust size distributions were rather flat in comparison to other comets, dN(R) proportional to R(sup -2.3 +/- 0.1), and the dust ejection speeds were approximately 0.5-1.5 m/s. The S-L9 fragments themselves were not detected directly, and upper limits on their radii are 1.0-1.5 km assuming an albedo a = 0.04. However, these fragments' vigorous production of dust, which ranges from 6 to 22 kg/s, places a lower limit of approximately 100 m on their radii at the moment of tidal breakup. Any fragments smaller than this limit, yet experiencing similar mass loss rates, would have dissipated prior to impact. Such bodies would fail to leave an impact scar at Jupiter's atmosphere, as was realized by fragments F, J, P(sub 1), P(sub 2), T, and U.

Hahn, Joseph M.↗
NASA NTRS2001

Passive and Active Stabilization of Liquid Bridges in Low Gravity

Tests are planned in the low gravity environment of the International Space Station (ISS) of new methods for the suppression of the capillary instability of liquid bridges. Our suppression methods are unusual in that they are not limited to liquid bridges having very special properties and may impact a variety of low-gravity and earth-based technologies. There are two main approaches to be investigated: (1) Passive Acoustic Stabilization (PAS); and (2) Active Electrostatic Stabilization (AES). In PAS, the suppression of the mode growth is accomplished by placing the bridge in an acoustic field having the appropriate properties such that the acoustic radiation pressure automatically pulls outward on the thinnest portion of the bridge. In AES, the bridge deformation is sensed optically and counteracted by actively adjusting the electrostatic Maxwell stresses via two ring electrodes concentric with the slightly conducting bridge to offset the growth of the unstable mode. While the present work emphasizes cylindrical bridges, the methods need not be restricted to that case. The methods to be explored are relevant to the suppression of capillary instabilities in floating zone crystal growth, breakup of liquid jets and columns, bubbles, and annular films as well as the management of coolants or propellants in low-gravity.

Marston, Philip L.↗
NASA NTRS2001

The Effects of Aircraft Wake Dynamics on Contrail Development

Results of large-eddy simulations of the development of young persistent ice contrails are presented, concentrating on the interactions between the aircraft wake dynamics and the ice cloud evolution over ages front a few seconds to approx. 30 min. The 3D unsteady evolution of the dispersing engine exhausts, trailing vortex pair interaction and breakup, and subsequent Brunt-Vaisala oscillations of the older wake plume are modeled in detail in high-resolution simulations, coupled with it bulk microphysics model for the contrail ice development. The simulations confirm that the early wake dynamics can have a strong influence on the properties of persistent contrails even at late times. The vortex dynamics are the primary determinant of the vertical extent of the contrail (until precipitate ton becomes significant): and this together with the local wind shear largely determines the horizontal extent. The ice density, ice crystal number density, and a conserved exhaust tracer all develop and disperse in different fashions from each other. The total ice crystal number can be significantly reduced due to adiabatic compression resulting from the downward motion of the vortex system, even for ambient conditions that are substantially supersaturated with respect to ice. The fraction of the initial ice crystals surviving, their spatial distribution and the ice mass distribution are all sensitive to the aircraft type, ambient humidity, assumed initial ice crystal number, and ambient turbulence conditions. There is a significant range of conditions for which a smaller transport such as a B737 produces as significant a persistent contrail as a larger transport such as a B747, even though the latter consumes almost five times as much fuel. The difficulties involved in trying to minimize persistent contrail production are discussed.

Lewellen, D. C.↗
NASA NTRS2002

Frequency Analyzed Laser Light Scattering (Falls)

A simple and robust light scattering technique for spray diagnostics was developed. The technique provides a spray's intact length and breakup frequency. Scattered light is collected from a laser sheet that illuminates a portion of the spray. The scattered light is imaged on a fast photodiode whose output is sampled at 20 kHz. Fast Fourier transform (FFT) analysis is performed on the sampled signal. The shape of the FFT produced can be correlated to one of three distinct spray regimes: intact sheet, ligaments, or droplets. Strobe and laser illuminated photographs are also shown for reference.

Rhys, Noah O.↗
NASA NTRS2002

Spin-Up Instability of a Levitated Molten Drop in MHD-Flow Transition to Turbulence

When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.

Abedian, B.↗
NASA NTRS2002

Ice Shelves and Landfast Ice on the Antarctic Perimeter: Revised Scope of Work

Ice shelves respond quickly and profoundly to a warming climate. Within a decade after mean summertime temperature reaches approx. O C and persistent melt pending is observed, a rapid retreat and disintegration occurs. This link was documented for ice shelves in the Antarctic Peninsula region (the Larsen 'A', 'B' and Wilkins Ice shelves) by the results of a previous grant under ADRO-1. Modeling of ice flow and the effects of meltwater indicated that melt pending accelerates shelf breakup by increasing fracture penetration. SAR data supplemented an AVHRR- and SSM/I-based image analysis of extent and surface characteristic changes. This funded grant is a revised, scaled-down version of an earlier proposal under the ADRO-2 NRA. The overall objective remains the same: we propose to build on the previous study by examining other ice shelves of the Antarctic and incorporate an examination of the climate-related characteristics of landfast ice. The study now considers just a few shelf and fast ice areas for study, and is funded for two years. The study regions are the northeastern Ross Ice Shelf, the Larsen 'B' and 'C' shelves, fast ice and floating shelf ice in the Pine Island Glacier area, and fast ice along the Wilkes Land coast. Further, rather than investigating a host of shelf and fast ice processes, we will home in on developing a series of characteristics associated with climate change over shelf and fast ice areas. Melt pending and break-up are the end stages of a response to a warming climate that may begin with increased melt event frequency (which changes both albedo and emissivity temporarily), changing firn backscatter (due to percolation features), and possibly increased rifting of the shelf surface. Fast ice may show some of these same processes on a seasonal timescale, providing insight into shelf evolution.

Scambos, Ted↗
NASA NTRS2002

Experimental Investigation of Water Droplet Impingement on Airfoils, Finite Wings, and an S-duct Engine Inlet

Validation of trajectory computer codes, for icing analysis, requires experimental water droplet impingement data for a wide range of aircraft geometries as well as flow and icing conditions. This report presents improved experimental and data reduction methods for obtaining water droplet impingement data and provides a comprehensive water droplet impingement database for a range of test geometries including an MS(1)-0317 airfoil, a GLC-305 airfoil, an NACA 65(sub 2)-415 airfoil, a commercial transport tail section, a 36-inch chord natural laminar flow NLF(1)-0414 airfoil, a 48-inch NLF(1)-0414 section with a 25 percent chord simple flap, a state-of-the-art three-element high lift system, a NACA 64A008 finite span swept business jet tail, a full-scale business jet horizontal tail section, a 25 percent-scale business jet empennage, and an S-duct turboprop engine inlet. The experimental results were obtained at the NASA Glenn Icing Research Tunnel (IRT) for spray clouds with median volumetric diameter (MVD) of 11, 11.5, 21, 92, and 94 microns and for a range of angles of attack. The majority of the impingement experiments were conducted at an air speed of 175 mph corresponding to a Reynolds number of approximately 1.6 million per foot. The maximum difference of repeated tests from the average ranged from 0.24 to 12 percent for most of the experimental results presented. This represents a significant improvement in test repeatability compared to previous experimental studies. The increase in test repeatability was attributed to improvements made to the experimental and data reduction methods. Computations performed with the LEWICE-2D and LEWICE-3D computer codes for all test configurations are presented in this report. For the test cases involving median volumetric diameters of 11 and 21 microns, the correlation between the analytical and experimental impingement efficiency distributions was good. For the median volumetric diameters of 92 and 94-micron cases, however, the analysis produced higher impingement efficiencies and larger impingement limits than the experiment. It is speculated that this discrepancy is due to droplet splashing and breakup experienced by large droplets during impingement.

Papadakis, Michael↗
NASA NTRS2003

Investigation of Dynamic Oxygen Adsorption in Molten Solder Jetting Technology

Surface tension forces play a critical role in fluid dynamic phenomena that are important in materials processing. The surface tension of liquid metals has been shown to be very susceptible to small amounts of adsorbed oxygen. Consequently, the kinetics of oxygen adsorption can influence the capillary breakup of liquid-metal jets targeted for use in electronics assembly applications, where low-melting-point metals (such as tin-containing solders) are utilized as an attachment material for mounting of electronic components to substrates. By interpreting values of surface tension measured at various surface ages, adsorption and diffusion rates of oxygen on the surface of the melt can be estimated. This research program investigates the adsorption kinetics of oxygen on the surface of an atomizing molten-metal jet. A novel oscillating capillary jet method has been developed for the measurement of dynamic surface tension of liquids, and in particular, metal melts which are susceptible to rapid surface degradation caused by oxygen adsorption. The experimental technique captures the evolution of jet swells and necks continuously along the jet propagation axis and is used in conjunction with an existing linear, axisymmetric, constant-property model to determine the variation of the instability growth rate, and, in turn, surface tension of the liquid as a function of surface age measured from the exit orifice. The conditions investigated so far focus on a time window of 2-4ms from the jet orifice. The surface properties of the eutectic 63%Sn-37%Pb solder alloy have been investigated in terms of their variation due to O2 adsorption from a N2 atmosphere containing controlled amounts of oxygen (from 8 ppm to 1000 ppm). The method performed well for situations where the oxygen adsorption was low in that time window. The value of surface tension for the 63Sn-37Pb solder in pure nitrogen was found to be 0.49 N/m, in good agreement with previously published work. A characteristic time of O(1ms) or less was determined for the molten-metal surface to be saturated by oxygen at 1000 ppm concentration in N2.

Megaridis, Constantine M.↗
NASA NTRS2003

Development of A Monte Carlo Radiation Transport Code System For HEDS: Status Update

Modifications of the Monte Carlo radiation transport code HETC are underway to extend the code to include transport of energetic heavy ions, such as are found in the galactic cosmic ray spectrum in space. The new HETC code will be available for use in radiation shielding applications associated with missions, such as the proposed manned mission to Mars. In this work the current status of code modification is described. Methods used to develop the required nuclear reaction models, including total, elastic and nuclear breakup processes, and their associated databases are also presented. Finally, plans for future work on the extended HETC code system and for its validation are described.

Townsend, Lawrence W.↗
NASA NTRS2003

Developing of a New Atmospheric Ionizing Radiation (AIR) Model

As a result of the research leading to the 1998 AIR workshop and the subsequent analysis, the neutron issues posed by Foelsche et al. and further analyzed by Hajnal have been adequately resolved. We are now engaged in developing a new atmospheric ionizing radiation (AIR) model for use in epidemiological studies and air transportation safety assessment. A team was formed to examine a promising code using the basic FLUKA software but with modifications to allow multiple charged ion breakup effects. A limited dataset of the ER-2 measurements and other cosmic ray data will be used to evaluate the use of this code.

Clem, John M.↗
NASA NTRS2003

A Nonlinear Model for Fuel Atomization in Spray Combustion

Most gas turbine combustion codes rely on ad-hoc statistical assumptions regarding the outcome of fuel atomization processes. The modeling effort proposed in this project is aimed at developing a realistic model to produce accurate predictions of fuel atomization parameters. The model involves application of the nonlinear stability theory to analyze the instability and subsequent disintegration of the liquid fuel sheet that is produced by fuel injection nozzles in gas turbine combustors. The fuel sheet is atomized into a multiplicity of small drops of large surface area to volume ratio to enhance the evaporation rate and combustion performance. The proposed model will effect predictions of fuel sheet atomization parameters such as drop size, velocity, and orientation as well as sheet penetration depth, breakup time and thickness. These parameters are essential for combustion simulation codes to perform a controlled and optimized design of gas turbine fuel injectors. Optimizing fuel injection processes is crucial to improving combustion efficiency and hence reducing fuel consumption and pollutants emissions.

Liu, Nan-Suey↗
NASA NTRS2003

A Cartesian Adaptive Level Set Method for Two-Phase Flows

In the present contribution we develop a level set method based on local anisotropic Cartesian adaptation as described in Ham et al. (2002). Such an approach should allow for the smallest possible Cartesian grid capable of resolving a given flow. The remainder of the paper is organized as follows. In section 2 the level set formulation for free surface calculations is presented and its strengths and weaknesses relative to the other free surface methods reviewed. In section 3 the collocated numerical method is described. In section 4 the method is validated by solving the 2D and 3D drop oscilation problem. In section 5 we present some results from more complex cases including the 3D drop breakup in an impulsively accelerated free stream, and the 3D immiscible Rayleigh-Taylor instability. Conclusions are given in section 6.

Ham, F.↗
NASA NTRS2003

Large-eddy Simulation of Stratocumulus-topped Atmospheric Boundary Layers with Dynamic Subgrid-scale Models

The objective of the present study is to evaluate the dynamic procedure in LES of stratocumulus topped atmospheric boundary layer and assess the relative importance of subgrid-scale modeling, cloud microphysics and radiation modeling on the predictions. The simulations will also be used to gain insight into the processes leading to cloud top entrainment instability and cloud breakup. In this report we document the governing equations, numerical schemes and physical models that are employed in the Goddard Cumulus Ensemble model (GCEM3D). We also present the subgrid-scale dynamic procedures that have been implemented in the GCEM3D code for the purpose of the present study.

Senocak, Inane↗
NASA NTRS2003

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗