SEARCH · Search NASA

Results for “TRANSIENT OSCILLATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

301 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

NASA NTRS2021

Plumelets: Dynamic Filamentary Structures in Solar Coronal Plumes

Solar coronal plumes long seemed to possess a simple geometry supporting spatially coherent, stable outflow without significant fine structure. Recent high-resolution observations have challenged this picture by revealing numerous transient, small-scale, collimated out-flows (“jetlets”) at the base of plumes. The dynamic filamentary structure of solar plumes above these outflows, and its relationship with the overall plume structure, have remained largely unexplored. We analyzed the statistics of continuously observed fine structure in-side a single representative bright plume within a mid-latitude coronal hole during 2016 July2-3. By applying advanced edge-enhancement and spatiotemporal analysis techniques to ex-tended series of high-resolution images from the Solar Dynamics Observatory’s Atmospheric Imaging Assembly, we determined that the plume was composed of numerous time-evolving filamentary substructures, referred to as “plumelets” in this paper, that accounted for most of the plume emission. The number of simultaneously identifiable plumelets was positively correlated with plume brightness, peaked in the fully formed plume, and remained saturated thereafter. The plumelets had transverse widths of 10 Mm and intermittently supported upwardly propagating periodic disturbances with phase speeds of 190-260 km s−1and longitudinal wavelengths of 55-65 Mm. The characteristic frequency (3.5 mHz) is commensurate with that of solar p-modes. Oscillations in neighboring plumelets are uncorrelated, indicating that the waves could be driven by p-mode flows at spatial scales smaller than the plumelet separation. Multiple independent sources of outflow within a single coronal plume should impart significant fine structure to the solar wind that may be detectable by Parker Solar Probe and Solar Orbiter.

V M Uritsky↗

NASA NTRS1994

Research on Orbital Plasma-Electrodynamics (ROPE)

Since the development of probe theory by Langmuir and Blodgett, the problem of current collection by a charged spherically or cylindrically symmetric body has been investigated by a number of authors. This paper overviews the development of a fully three-dimensional particle simulation code which can be used to understand the physics of current collection in three dimensions and can be used to analyze data resulting from the future tethered satellite system (TSS). According to the TSS configurations, two types of particle simulation models were constructed: a simple particle simulation (SIPS) and a super particle simulation (SUPS). The models study the electron transient response and its asymptotic behavior around a three dimensional, highly biased satellite. The potential distribution surrounding the satellite is determined by solving Laplace's equation in the SIPS model and by solving Poisson's equation in the SUPS model. Thus, the potential distribution in space is independent of the density distribution of the particles in the SUPS model but it does depend on the density distribution of the particles in the SUPS model. The evolution of the potential distribution in the SUPS model is described. When the spherical satellite is charged to a highly positive potential and immersed in a plasma with a uniform magnetic field, the formation of an electron torus in the equatorial plane (the plane in perpendicular to the magnetic field) and elongation of the torus along the magnetic field are found in both the SIPS and the SUPS models but the shape of the torus is different. The areas of high potential that exist in the polar regions in the SUPS model exaggerate the elongation of the electron torus along the magnetic field. The current collected by the satellite for different magentic field strengths is investigated in both models. Due to the nonlinear effects present in SUPS, the oscillating phenomenon of the current collection curve during the first 10 plasma periods can be seen (this does not appear in SIPS). From the parametric studies, it appears that the oscillating phenomenon of the current collection curve occurs only when the magnetic field strength is less than 0.2 gauss for the present model.

Wu, S. T.↗

NASA NTRS2013

Megadroughts in Southwestern North America in ECHO-G Millennial Simulations and Their Comparison to Proxy Drought Reconstructions

Simulated hydroclimate variability in millennium-length forced transient and control simulations from the ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G) coupled atmosphere-ocean general circulation model (AOGCM) is analyzed and compared to 1000 years of reconstructed Palmer drought severity index (PDSI) variability from the North American Drought Atlas (NADA). The ability of the model to simulate megadroughts in the North American southwest is evaluated. (NASW: 25deg42.5degN, 125deg-105degW). Megadroughts in the ECHO-G AOGCM are found to be similar in duration and magnitude to those estimated from the NADA. The droughts in the forced simulation are not, however, temporally synchronous with those in the paleoclimate record, nor are there significant differences between the drought features simulated in the forced and control runs. These results indicate that model-simulated megadroughts can result from internal variability of the modeled climate system rather than as a response to changes in exogenous forcings. Although the ECHO-G AOGCM is capable of simulating megadroughts through persistent La Nina-like conditions in the tropical Pacific, other mechanisms can produce similarly extreme NASW moisture anomalies in the model. In particular, the lack of low-frequency coherence between NASW soil moisture and simulated modes of climate variability like the El Nino-Southern Oscillation, Pacific decadal oscillation, and Atlantic multidecadal oscillation during identified drought periods suggests that stochastic atmospheric variability can contribute significantly to the occurrence of simulated megadroughts in the NASW. These findings indicate that either an expanded paradigm is needed to understand multidecadal hydroclimate variability in the NASW or AOGCMs may incorrectly simulate the strength and/or dynamics of the connection between NASW hydroclimate variability and the tropical Pacific.

variability↗

NASA NTRS2019

Tropopause Laminar Cirrus and Its Role in the Lower Stratosphere Total Water Budget

Laminar cirrus are thin, extensive, isolated layers of ice clouds frequently observed in the tropical tropopause layer. Widespread laminar cirrus significantly affects tropical tropopause layer total water and thermal budget. In this study, we extract laminar cirrus from the Cloud‐Aerosol Lidar with Orthogonal Polarization Level 1 attenuated total backscatter images for January 2009, in order to characterize statistical properties of laminar cirrus cloud length, base, thickness, optical depth, and layer partial ice water path. These characteristics are used to develop an algorithm identifying laminar cirrus automatically from the Cloud‐Aerosol Lidar with Orthogonal Polarization Level 2 layer product for 2008–2017. The nearly 10‐year records reveal that tropopause laminar cirrus occurrence (30–40% of total cirrus) is strongly anticorrelated with the tropopause temperatures in that colder tropopause in frequent (super)saturation during boreal winter favors in situ formation of clouds. Interannually, anomalously warmer troposphere temperature (ΔT), easterly shear of the quasi‐biennial oscillation, and stronger upwelling branch of the Brewer‐Dobson circulation enhance laminar cirrus formation via cooling of the tropopause. The tropopause laminar cirrus carries ~0.05 mg/m(exp 3) (~0.5 ppmv) of ice water content during boreal winter and <0.01 mg/m(exp 3) during summer, which is anticorrelated with the seasonal variations of water vapor (H2O) observed by the Microwave Lime Sounder, indicating a temperature‐regulated partition between vapor and ice. Interannually, in cirrus‐rich region 1 ppmv decrease in H2O corresponds to 0.2–0.3 ppmv increase in ice water content. Frequently situated in (super)saturated air, tropopause laminar cirrus are likely to survive multiple lifecycles of the sublimation‐deposition processes, and may contribute up to 10% to the total water budget in the lower stratosphere. Satellites constantly observe thin, isolated, extensive layer of cirrus around the tropopause. The so‐called “laminar” cirrus occurrence and their ice amount are strongly regulated by temperature, such that colder temperatures favor more frequent (super)saturation, which results in more frequent laminar cirrus with more ice amount and therefore less water vapor. In this study we analyze laminar cirrus and water vapor from the Cloud‐Aerosol Lidar with Orthogonal Polarization and Microwave Lime Sounder observations, and hypothesize that laminar cirrus could act as an important transient water storage and contribute to the total water budget in the lower stratosphere.

Tao Wang↗

NASA NTRS2001

A Decade in the Life of EXO 2030+375: A Multiwavelength Study of an Accreting X-ray Pulsar

Using BATSE and RXTE observations from 1991 April to 2001 August we have detected 71 outbursts from 82 periastron passages of EXO 2030+375, a 42-second transient X-ray pulsar with a Be star companion, including several outbursts from 1993 August to 1996 April when the source was previously believed to be quiescent. Combining BATSE, RXTE, and EXOSAT data we have derived an improved orbital solution. Applying this solution results in a smooth profile for the spin-up rate during the giant outbursts and results in evidence for a correlation between the spin-up rate and observed flux in the brighter BATSE outbursts. Infrared and H(alpha) measurements show a decline in the density of the circumstellar disk around the Be star. This decline is followed by a sudden drop in the X-ray flux and a turn-over from a spin-up trend to spin-down in the frequency history. This is the first Be/X-ray binary which shows an extended interval, about 2.5 years, where the global trend is spin-down, but the outbursts continue. In 1995 the orbital phase of EXO 2030+375's outbursts shifted from peaking about 6 days after periastron to peaking before periastron. The outburst phase slowly recovered to peaking at about 2.5 days after periastron. We interpret this shift in orbital phase followed by a slow recovery as evidence for a global one-armed oscillation propagating in the Be disk. This is further supported by changes in the shape of the H(alpha) profile which are commonly believed to be produced by a reconfiguration of the Be disk. The truncated viscous decretion disk model provides an explanation for the long series of normal outbursts and the evidence for an accretion disk in the brighter normal outbursts. Long-term multi-wavelength observations such as these clearly add considerably to our knowledge of Be/X-ray binaries and the relationship between optical, infrared and X-ray observations.

Wilson, Colleen A.↗

NASA NTRS2020

Large-Scale Influences on Atmospheric River–Induced Extreme Precipitation Events Along the Coast of Washington State

Transient, narrow plumes of strong water vapor transport, referred to as Atmospheric Rivers (ARs) are responsible for much of the precipitation along the west coast of the United States. Along the coast of Oregon and Washington, the most intense cool season precipitation events are almost always induced by an AR and can result in detrimental impacts on society due to mudslides and flooding. It is therefore important to understand the large scale influence on extreme AR events so that they can be accurately predicted on timescales ranging from numerical weather prediction to seasonal forecasts. Here, characteristics of ARs that result in observed extreme precipitation events are compared to typical ARs on the coast of Washington State using data from the Modern Era Retrospective analysis for Research and Applications, Version 2. In addition to more intense water vapor transport, notable differences in the synoptic scale forcing are present during extreme precipitation events that are not present during typical AR events. In particular, an anomalously deep low pressure system is stationed to the west in the Gulf of Alaska, alongside a jet streak overhead. Attention will also be given to subseasonal and seasonal teleconnection patterns that are known to influence the weather in the Pacific Northwest of the United States. While little influence can be seen from the phase of the El Nino Southern Oscillation, Pacific Decadal Oscillation, and Pacific North American Pattern, the Madden Julian Oscillation (MJO) can play a role in determining the strength of precipitation associated with in AR on the Washington Coast. Lastly, interactions between the MJO and other teleconnection patterns will be explored to determine key features that should be investigated when making subseasonal predictions for AR activity and the associated precipitation in the Pacific Northwest.

Collow, Allison↗

NASA NTRS1994

The volcanic signal in Goddard Institute for Space Studies three-dimensional model simulations

Transient calculations of the Goddard Institute for Space Studies general circulation model for the climatic signal of volcanic eruptions are analyzed. By compositing the output for two different volcanoes for scenario A and five different volcanoes for scenario B, the natural variability is suppressed and the volcanic signals are extracted. Significant global mean surface air temperature cooling and precipitation reduction are found for several years following the eruptions, with larger changes in the Northern Hemisphere (NH) than in the Southern Hemisphere. The global-average temperature response lasts for more than four years, but the precipitation response disappears after three years. The largest cooling in the model occurs in the NH summer of the year after spring eruptions. Significant zonal-average temperature reductions begin in the tropics immediately after the eruptions and extend to 45 deg S - 45 deg N in the year after the eruptions. In the second year, cooling is still seen from 30 deg S to 30 deg N. Because of the low variability in this model as compared to the real world, these signals may appear more significant here than they would by attempting to isolate them using real data. The results suggest that volcanoes can enhance the drought in the Sahel. No evidence was found that stratospheric aerosols from the low-atitude volcanic eruptions can trigger El Nino-Southern Oscillation (ENSO) events in this model.

Robock, Alan↗

NASA NTRS2021

Photospheric Radius Expansion and a Double-peaked Type-I X-Ray Burst from GRS 1741.9-2853

We present an analysis of two type-I X-ray bursts observed by NuSTAR originating from the very faint transient neutron star (NS) low-mass X-ray binary GRS 1741.9–2853 during a period of outburst in 2020 May. We show that the persistent emission can be modeled as an absorbed, Comptonized blackbody in addition to Fe Kα emission, which can be attributed to relativistic disk reflection. We measure a persistent bolometric, unabsorbed luminosity of L(bol) = 7.03(+0.04, -0.05) x 10^(36) erg/s, assuming a distance of 7 kpc, corresponding to an Eddington ratio of 4.5%. This persistent luminosity combined with light-curve analysis leads us to infer that the bursts were the result of pure He burning rather than mixed H/He burning. Time-resolved spectroscopy reveals that the bolometric flux of the first burst exhibits a double-peaked structure, placing the source within a small population of accreting NSs that exhibit multiple-peaked type-I X-ray bursts. We find that the second, brighter burst shows evidence for photospheric radius expansion (PRE) and that at its peak, this PRE event had an unabsorbed bolometric flux of F(peak)= 2.94 (+0.28, -0.26) 10^(-8) erg/sq.cm s. This yields a new distance estimate of d = 9.0 ± 0.5 kpc, assuming that this corresponds to the Eddington limit for pure He burning on the surface of a canonical NS. Additionally, we performed a detailed timing analysis that failed to find evidence for quasi-periodic oscillations or burst oscillations, and we place an upper limit of 16% on the rms variability around 589 Hz, the frequency at which oscillations have previously been reported.

Sean N. Pike↗

NASA NTRS2020

The Fast Radio Burst Luminosity Function and Death Line in the Low-Twist Magnetar Model

We explore the burst energy distribution of fast radio bursts(FRBs)in the low-twist magnetar model ofWadiasingh & Timokhin(WT19). Motivated by the power-lawfluence distributions of FRB121102, we proposean elementary model for the FRB luminosity function of individual repeaters with an inversion protocol thatdirectly relates the power-law distribution index of magnetar short burstfluences to that for FRBs. The protocolindicates that the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar shortbursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts(required inWT19)for individual repeaters implies the predicted burstfluence distribution is narrow,3 decades for yieldingstrains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and chargestarvation, we obtain a death line for FRBs, which segregates magnetars from the normal pulsar population,suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individualmagnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-lawluminosity function’s low-energy character depends on the population model, while the high-energy index tracesthat of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalentenergy/luminosity function peaks at∼1037–1040erg with a low-energy cutoff at∼1037erg. Lastly, we considerthe localfluence distribution of FRBs andfind that it can constrain the subset of FRB-producing magnetarprogenitors. Our model suggests that improvements in sensitivity may reveal aflattening of the global FRBfluencedistribution and saturation in FRB rates.

Cosmology↗

NASA NTRS2002

XTE J1908+094

Goddard Space Flight Center reports the serendipitous discovery of a new x-ray transient, XTE J1908+094, in RXTE (Rossi X-ray Timing Explorer) PCA (Proportional Counter Array) observations of the soft-gamma-ray repeater SGR 1900+14, triggered following the burst activity on Feb. 17-18 (GCN 1253). These observations failed to detect the 5.2-s SGR pulsations, pointing towards a possible new source as the origin of the high x-ray flux. An RXTE PCA scan of the region around SGR 1900+14 on Feb. 21 was consistent with emission only from known sources (and no new sources). However, the scans required SGR 1900+14 to be 20 times brighter than its quiescent flux level (GCN 1256). A Director's Discretionary Time Chandra observation on Mar. 11 showed that the SGR was quiescent and did not reveal any new source within the Chandra ACIS (Advanced CCD (charge coupled device) Imaging Spectrometer) field-of-view. A subsequent RXTE PCA scan on Mar. 17, taken in combination with the first scan, required that a new source be included in the fit. The best-fit position is R.A. 19h 08m 50s, Decl. = +9 22 deg .5 (equinox J2000.0; estimated 2 deg systematic error radius), or approximately 24 deg away from the SGR source. The source spectrum (2-30 kev) can be best fit with a power-law function including photoelectric absorption (column density N_h = 2.3 x 10(exp 22), photon index = 1.55). Iron line emission is present, but may be due to the Galactic ridge. Between Feb. 19 and Mar. 17, the source flux (2-10 keV) has risen from 26 to 64 mCrab. The power spectrum is flat between 1 mHz and 0.1 Hz, falling approximately as 1/f**0.5 up to 1 Hz. At 1 Hz is seen a broad quasiperiodic oscillation peak and a break to a 1/f**2 power law, which continues to 4 Hz. The fractional rms (root mean square) amplitude from 1 mHz to 4 Hz is 43 percent. No coherent pulsations are seen between 0.001 and 1024 Hz. The authors conclude that XTE J1908+094 is a new blackhole candidate.

Woods, P. M.↗

NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

NASA NTRS2005

Dynamic Weakening (Extinction) of Simple Hydrocarbon-air Counterflow Diffusion Flames by Oscillatory Inflows

This study of laminar non-premixed HC-air flames used an Oscillatory-input Opposed Jet Burner (OOJB) system developed from a previously well-characterized 7.2-mm Pyrex-nozzle OJB system. Over 600 dynamic Flame Strength (FS) measurements were obtained on unanchored (free-floating) laminar Counterflow Diffusion Flames (CFDF's). Flames were stabilized using plug inflows having steady-plus-sinusoidal axial velocities of varied magnitude, frequency, f, up to 1600 Hz, and phase angle from 0 (most data) to 360 degrees. Dynamic FS is defined as the maximum average air input velocity (U(sub air), at nozzle exit) a CFDF can sustain before strain-induced extinction occurs due to prescribed oscillatory peak-to-peak velocity inputs superimposed on steady inputs. Initially, dynamic flame extinction data were obtained at low f, and were supported by 25-120 Hz Hot-Wire cold-flow velocity data at nozzle exits. Later, expanded extinction data were supported by 4-1600 Hz Probe Microphone (PM) pk/pk P data at nozzle exits. The PM data were first obtained without flows, and later with cold stagnating flows, which better represent speaker-diaphragm dynamics during runs. The PM approach enabled characterizations of Dynamic Flame Weakening (DFW) of CFDF's from 8 to 1600 Hz. DFW was defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = - 100 d(U(sub air) / U(sub air),0Hz) / d(pkpk P). The linear normalization with respect to acoustic pressure magnitude (and steady state (SS) FS) led to a DFW unaffected by strong internal resonances. For the C2H4/N2-air system, from 8 to 20 Hz, DFW is constant at 8.52 plus or minus 0.20 (% weakening)/Pa. This reflects a quasi-steady flame response to an acoustically induced dU(sub air)/dP. Also, it is surprisingly independent of C2H4/N2 mole fraction due to normalization by SS FS. From 20 to approximately 150 Hz, the C2H4/N2 air-flames weakened progressively less, with an inflection at approximately 70 Hz, and became asymptotically insensitive (DFW approximately 0) at approximately 300 Hz, which continued to 1600 Hz. The DFW of CH4-air flames followed a similar pattern, but showed much greater weakening than C2H4/N2-air flames; i.e., the quasi-steady DFW (8 to approximately 15 Hz) was 44.3 %/Pa, or approximately 5x larger, even though the 0 Hz (SS) FS was only 3.0 x smaller. The quasi-steady DFW's of C3H8-air and C2H6-air were intermediate at 34.8 and 20.9 %Pa, respectively. The DFW profiles of all four fuels, at various frequencies, correlated well but non-linearly with respective SS FS's. Notably, the DFW profile for C3H8 air fell more rapidly in the range greater than 15 to 60 Hz, compared with the 1- and 2-carbon fuels. This may indicate a shift in chemical kinetics, and/or O2 transport to a flame that moved closer to the fuel-side. In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and correlate closely, but non-linearly, with Steady-State Flame Strengths at any given frequency. For reasons unknown, the dynamic flames didn't weaken more at intermediate frequencies (e.g., at 20-50 Hz) than they did at low frequencies (less than 15 Hz), where quasi-steady weakening appears to dominate. Quasi-steady flame weakening ostensibly represents a transient input strain rate maximum that just exceeds the steady-state strain-rate-limited extinction limit for a few cycles. Clearly, further detailed mechanistic understanding is needed in the fall-off region.

Pellett, G.↗