Search NASA⌕ Search

SEARCH · Search NASA

Results for “front propagation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

The Dynamics of Miscible Interfaces: Simulations

The goal of this experimental/computational investigation (joint with Prof Maxworthy at USC) has been to study the dynamics of miscible interfaces, both from a scientific and a practical point of view, and to prepare a related experiment to be flown on the International Space Station. In order to address these effects, we have focused experimental and computational investigations on miscible displacements in cylindrical capillary tubes, as well as in Hele-Shaw cells. Regarding the flow in a capillary tube, the question was addressed as to whether Korteweg stresses and/or divergence effects can potentially account for discrepancies observed between conventional Stokes flow simulations and experiments for miscible flows in capillary tubes. An estimate of the vorticity and streamfunction fields induced by the Kortewegs stresses was derived, which shows these stresses to result in the formation of a vortex ring structure near the tip of the concentration front. Through this mechanism the propagation velocity of the concentration front is reduced, in agreement with the experimental observations. Divergence effects, on the other hand, were seen to be very small, and they have a negligible influence on the tip velocity. As a result, it can be concluded that they are not responsible for the discrepancies between experiments and conventional Stokes simulations. A further part of our investigation focussed on the development of high-accuracy three-dimensional spectral element simulation techniques for miscible flows in capillary tubes, including the effects of variable density and viscosity. Towards this end, the conservation equations are treated in cylindrical coordinates.

Meiburg, Eckart↗
NASA NTRS1989

Numerical simulation of hyperbolic heat conduction with convection boundary conditions and pulse heating effects

The paper describes the numerical simulation of hyperbolic heat conduction with convection boundary conditions. The effects of a step heat loading, a sudden pulse heat loading, and an internal heat source are considered in conjunction with convection boundary conditions. Two methods of solution are presened for predicting the transient behavior of the propagating thermal disturbances. In the first method, MacCormack's predictor-corrector method is employed for integrating the hyperbolic system of equations. Next, the transfinite element method, which employs specially tailored elements, is used for accurately representing the transient response of the propagating thermal wave fronts. The agreement between the results of various numerical test cases validate the representative behavior of the thermal wave fronts. Both methods represent hyperbolic heat conduction behavior by effectively modeling the sharp discontinuities of the propagating thermal disturbances.

Glass, David E.↗
NASA NTRS1995

Nonlinear Propagating Features in the Stratosphere of Jupiter Generated by the Impact of SL-9

Among the more intriguing features observed by HST after the impacts of the fragments of Shoemaker-Levy 9 were rings, encircling several of the impact sites, that expanded radially at constant velocities of order 450-500 m/s. One credible suggestion has been that the rings are stratospheric manifestations of linear gravity waves seated in the troposphere (Ingersoll and Kanamori, 1995). In this scenario the waves are made visible by a condensible material of unknown composition. Based on simulations of the comet impact by Zahnle and Mac Low (1995), we suggest an alternative hypothesis, namely that the rings are nonlinear stratospheric gravity waves which break and entrain impact debris material in the breaking wave front. The computations produce just such a feature, with the nonlinear breaking wave front in the stratosphere propagating radially away from the impact site with the correct speed. The computed ring of material is located in an altitude region where the static stability profile has a local maximum and therefore forms a waveguide. Because entrained impact debris forms the ring feature in the model, a condensate of some unknown composition is not required in order to view the wave front. However, there is one severe problem with the computed nonlinear wave: in contrast to a linear wave, it slows down. The computed feature stops within 1000 seconds, while the observed feature lasted at least 10,000 seconds. Since the nonlinear breaking wave front is a natural product of our impact simulations, we are now investigating the conditions that will enable the breaking wave front to propagate further, in accord with the observations. The results of these computations will be presented.

Young, R. E.↗
DOE OSTI2022

Transition from sub-Rayleigh anticrack to supershear crack propagation in snow avalanches

Snow slab avalanches, characterized by a distinct, broad fracture line, are released following anticrack propagation in highly porous weak snow layers buried below cohesive slabs. The anticrack mechanism is driven by the volumetric collapse of the weak layer, which leads to the closure of crack faces and to the onset of frictional contact. Here, on the basis of snow fracture experiments, full-scale avalanche measurements and numerical simulations, we report the existence of a transition from sub-Rayleigh anticrack to supershear crack propagation. This transition follows the Burridge–Andrews mechanism, in which a supershear daughter crack nucleates ahead of the main fracture front and eventually propagates faster than the shear wave speed. Furthermore, we show that the supershear propagation regime can exist even if the shear-to-normal stress ratio is lower than the static friction coefficient as a result of the loss of frictional resistance during collapse. This finding shows that snow slab avalanches have fundamental similarities with strike-slip earthquakes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2005

A Multi-Instrument Measurement of a Mesospheric Bore at the Equator

We have made a comprehensive measurement of mesospheric bore phenomenon at the equator at Kototabang, Indonesia (0.2 deg S, 100.3 deg E), using an airglow imager, an airglow temperature photometer, a meteor radar, and the SABER instrument on board the TIMED satellite. The bore was detected in airglow images of both OH-band (peak emission altitude: 87 km) and 557.7-nm (96 km) emissions, as east-west front-like structure propagating northward with a velocity of 52-58 m/s. Wave trains with a horizontal wavelength of 30-70 km are observed behind the bore front. The airglow intensity decreases for all the mesospheric emissions of OI (557.7 nm), OH-band, O2-band (altitude: 94 km), and Na (589.3 nm) (90 km) after the bore passage. The rotational temperatures of both OH-band and O2-band also decrease approximately 10 K after the bore passage. An intense shear in northward wind velocity of 80m/s was observed at altitudes of 84-90 km by the meteor radar. Kinetic temperature profile at altitudes of 20-120 km was observed near Kototabang by TIMED/SABER. On the basis of these observations, we discuss generation and ducting of the observed mesospheric bore.

Shiokawa, K.↗
NASA NTRS1999

Moreton Waves

"Moreton waves," named for the observer who popularized them, are a solar phenomenon also known in scientific literature as "Moreton-Ramsey wave," "flare waves," "flare-associated waves," "MHD blast waves," "chromospheric shock fronts" and various other combinations of terms which connote violently propagating impulsive disturbances. It is unclear whether all of the observations to which these terms have been applied pertain to a single physical phenomenon: there has perhaps been some overlap between the observations and the assumed physical properties of the observed occurrence. Moreton waves are ideally observed in the wings of H alpha, and appear as semi-circular fronts propagating at speeds ranging from several hundred to over a thousand km/sec. They form an arc, or "brow shape" which can span up to 180 degrees. Extrapolating the speed and locations of the arc indicates that the phenomenon's origin intersects well with the impulsive phase of the associated H alpha flare (if the flare exhibits an impulsive phase). However, the arc may not form or may not be observable until it is tens of megameters from the flaring region, and subsequently can propagate to distances exceeding 100 megameters. The high speeds and distances of propagation, plus the associated radio and energetic particle observations, provided strong evidence of a coronal, rather than a chromospheric origin. The H alpha manifestation of the wave is assumed to be the "ground track" or "skirt" of a three-dimensional disturbance.

Thompson, B. J.↗
NASA NTRS1980

Numerical techniques for solving nonlinear instability problems in smokeless tactical solid rocket motors

The selection of a satisfactory numerical method for calculating the propagation of steep fronted shock life waveforms in a solid rocket motor combustion chamber is discussed. A number of different numerical schemes were evaluated by comparing the results obtained for three problems: the shock tube problems; the linear wave equation, and nonlinear wave propagation in a closed tube. The most promising method--a combination of the Lax-Wendroff, Hybrid and Artificial Compression techniques, was incorporated into an existing nonlinear instability program. The capability of the modified program to treat steep fronted wave instabilities in low smoke tactical motors was verified by solving a number of motor test cases with disturbance amplitudes as high as 80% of the mean pressure.

Baum, J. D.↗
DOE OSTI2021

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗
DOE OSTI2020

Configurational stability of a crack propagating in a material with mode-dependent fracture energy – Part II: Drift of fracture facets in mixed-mode I+II+III

In earlier papers (Leblond et al., 2011; 2019), we presented linear stability analyses of the coplanar propagation of a crack loaded in mixed-mode I+III, based on a “double” propagation criterion combining (Griffith, 1920)’s energetic condition and (Goldstein and Salganik, 1974)’s principle of local symmetry. The difference between the two papers was that in the more recent one, the local value of the critical energy-release-rate was no longer considered as a constant, but heuristically allowed to depend upon the ratio of the local mode III to mode I stress intensity factors. This led to a much improved, qualitatively acceptable agreement of theory and experiments, for the “threshold” value of the ratio of the unperturbed mode III to mode I stress intensity factors, above which coplanar propagation becomes unstable. In this paper, the analysis is extended to the case where a small additional mode II loading component is present in the initially planar configuration of the crack, generating a small, general kink of this crack from the moment it is applied. The main new effect resulting from presence of such a loading component is that the instability modes present above the threshold must drift along the crack front during its propagation. This prediction may be useful for future theoretical interpretations of a number of experiments where such a drifting motion was indeed observed.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Erosion of an electron-beam front in a long beam-plasma system

The front of an electron beam injected into a uniform, magnetized plasma is observed to erode rapidly because of scattering of the beam electrons via collective beam-plasma interactions. The propagation of the beam front is measured with temporal and axial resolution via beam excitation light. A broadened light front can result that propagates much slower than the injected beam particle speed.

Chan, L. Y.↗
NASA NTRS1997

The Effect of Gravity on the Combustion Synthesis of Porous Ceramics and Metal Matrix Composites

Combustion synthesis (self propagating, high temperature synthesis-SHS) is a novel technique that is capable of producing many advanced materials. The ignition temperature (Tig) of such combustion synthesis reactions is often coincident with that of the lowest melting point reactant. The resultant liquid metal wets and spreads around the other solid reactant particles of higher melting points, thereby improving the reactant contact and kinetics, followed by formation of the required compounds. This ignition initiates a combustion propagating wave whose narrow reaction front rapidly travels through the reactants. Since this process is highly exothermic, the heat released by combustion often melts the reactant particles ahead of the combustion front and ignites the adjacent reactant layer, resulting in a self-sustaining reaction. Whenever a fluid phase (liquid or gas) is generated by the reaction system, gravity-driven phenomena can occur. Such phenomena include convective flows of fluid by conventional or unstable convection and settling of the higher density phases. A combustion process is often associated with various kinds of fluid flow. For instance, if the SHS reaction is carried out under inert or reactive gas atmospheres, or a volatile, e.g., B2O3, is deliberately introduced as a reactant, convective flows of the gas will occur due to a temperature gradient existing in the atmosphere when a combustion wave is initiated. The increased gas flow will produce a porous (or expanded) SHS product. Owing to the highly exothermic nature of many SHS reactions, liquid phase(s) can also form before, at, or after the combustion front. The huge temperature gradient at the combustion front can induce convective flows (conventional or unstable) of the liquid phase. Each of these types of convective fluid flow can change the combustion behavior of the synthesizing reaction, and, therefore, the resultant product microstructure. In addition, when two or more phases of different density are produced at or ahead of the propagating combustion front settling of the higher density phase will occur resulting in a non-uniform product microstructure and properties.

Moore, J. J.↗
NASA NTRS1986

Simulations of gravity waves in frontal zones

Recent observations of frontal systems passing over the PROFS network in Colorado by Shapiro (1984) indicate the horizontal scale of some fronts can be on the order of a few kilometers or less. Motivated by Shapiros results, esearchers re-ran an earlier numerical simulation of frontogenesis by Williams (1972) (where a very simple stretching deformation forcing of a front is considered) using as high a resolution in both the vertical and horizontal as was possible. Highest resolution that we considered consisted of a grid with a vertical spacing of 32 meters and a horizontal spacing of 260 meters in the immediate vicinity of the surface front. The purpose was to examine development of the frontal structures as the scale of the front became very small. Two sources for banded structures in clouds were found: (1) in the stationary waves above the front and (2) in waves propagating upward from breaking waves under the front.

Gall, R.↗
DOE OSTI2024

DNS of ignition and flame stabilization in a simplified gas turbine premixer

With the increasing need for fuel flexibility, mitigation of auto-ignition (AI) inside gas turbine (GT) premixers becomes crucial. They must be designed to yield a sufficiently homogeneous fuel-air mixture to achieve low emissions while at the same time avoiding the occurrence of AI and subsequent flame stabilization. This challenge requires a detailed understanding of turbulent mixing and chemistry interactions. In the present work, a direct numerical simulation (DNS) of an array of jets in crossflow (JICF), representative of an industrial GT premixer, is reported to shed light on these complex phenomena. It is found that AI kernels form in the aft part of the premixer and coalesce into a flame front that then propagates upstream, mainly through the boundary layer, and successively engulfs the jets. This, therefore, suggests a significant role of the jet array pattern on the flame stabilization. It is noted that AI kernels continue to form independently during the whole time of the simulation. To clarify the contribution of AI and diffusion in the ignition kernels and the main flame, chemical explosive mode analysis (CEMA) is employed jointly with a kernel tracking algorithm. It is found that during the initial formation of the flame, many ignition kernels form in mixtures with low scalar dissipation rate and large contribution from AI mode. As they quickly grow, they merge into a single flame front that becomes increasingly more diffusion-assisted over time, balancing the AI mode. Turbulence is shown to have a significant enhancing effect in lean premixed flames, but further analysis is required to fully characterize it. These findings are relevant for the industrial premixer studied, and also for novel micromix concepts that may be used in the next generation of GT combustion systems.

ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1997

Premixed Gas Combustion: An Excitable System

Rotating spiral and target patterns have been observed experimentally on freely-propagating premixed gas flames in large diameter tubes at normal gravity (1-g). These modes of propagation occur in near-limit mixtures which have a Lewis number (Le, defined as the ratio of the thermal diffusivity of the cold mixture to the mass diffusivity of the scarce component into the mixture) sufficiently greater than one. However, at 1-g, buoyant flows strongly distort the flame curvature, hydrodynamics (thus stretch) and convective transport of species and heat. In turn, these alter the critical Le required for onset of instability. To isolate and better understand the mechanisms which drive the observed patterns and their dynamics, 1-g and microgravity (micro-g) experiments are being conducted to determine: (1) the structure and dynamics of the patterns, (2) a map of the critical Le and heat loss for their occurrence, (3) the relative significance of the chemical kinetics, and (4) the effect of curvature (local wave and global flame front) on wave propagation. With this in hand, we will be better prepared to discuss an additional mode, a state of 'chemical turbulence,' which seems to be the ultimate fate of many of these near-limit flames prior to extinction.

Pearlman, Howard↗
DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗
DOE OSTI2021

Replacement of Calcium Carbonate Polymorphs by Cerussite

Calcium carbonate (CaCO 3 ) polymorphs, calcite, aragonite, and vaterite, serve as a major sink to retain various metal ions in natural and engineered systems. Here, we visualize the systematic trends in reactivities of calcite, vaterite, and aragonite to Pb 2+ dissolved in acidic aqueous solutions using in situ optical microscopy combined with ex situ scanning electron and transmission X-ray microscopies. All three polymorphs undergo pseudomorphic replacement by cerussite (PbCO 3 ) but with distinct differences in the evolution of their morphologies. The replacement of calcite and aragonite occurs through the formation of a pseudomorphic cerussite shell (typically 5–10 μm thick) followed by a slower inward propagation of reaction fronts through a thin solution gap (~0.1 μm wide) between the shell and the CaCO 3 substrate. The replacement of vaterite is characterized by the formation of a thinner cerussite shell (≤1 μm thick) and a larger cavity between the shell and the host mineral. These systematic differences in cerussite morphology for different CaCO 3 polymorphs are explained by the relative dissolution and precipitation rates of the reactant and product minerals, coupled with the role of ion transport through the cerussite shells. We also find that the replacement of calcite by cerussite is the slowest when all three polymorphs coexisted. Furthermore, our results provide mechanistic insights into the growth mode of cerussite on dissolving calcium carbonate and demonstrate these CaCO 3 polymorphs as promising substrate materials for removal and recycling of Pb from acidic polluted water and industrial effluents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Characterizing Greater Houston's Aerosol by Air Mass During TRACER

During the TRacking Aerosol and Convection interaction ExpeRiment (TRACER), a suite of aerosol and cloud measurements were made using the Texas A&M Rapid Onsite Atmospheric Measurement Van (ROAM-V) mobile instrument platform. Joint ROAM-V/radiosonde deployments focused on sampling polluted marine and continental air masses to understand summertime convection along and across the sea-breeze front as it propagated through Houston. The polluted marine air mass at Seawolf Park, Galveston, TX was defined by two characteristic aerosol populations. Although the background total particle concentrations averaged 2,500 cm −3 , the air mass was also influenced by frequent but irregular periods of ship emission with significantly higher aerosol concentrations at times exceeding 34,000 cm −3 . Ship emission influenced periods, typically lasting <10 min, contained smaller, less hygroscopic particles resulting in a 69% relative decline in cloud condensation nuclei (CCN) activated fraction at 1% supersaturation compared to background periods. Measurements in continental air masses northwest of Houston revealed an average aerosol concentration of 5,208 cm −3 with a κ value near 0.1 reasonably describing the CCN population, independent of particle size. In continental air masses, substantial differences in particle size and CCN activation over small distances (<42 miles between sites) suggest considerable site-to-site variability in addition to the expected day-to-day differences. This small-scale variability makes it difficult to generalize continental air mass aerosol properties. Both coastal and inland locations had effective ice nucleating particles, but inland deployments observed the warmest nucleation temperature at −15.6°C compared to −17.8°C close to the coast. These measurements can reduce uncertainties in regional convection allowing models to improve understanding of aerosol-cloud interactions.

54 ENVIRONMENTAL SCIENCES↗