Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spreading models”

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 253 records · Page 14

Modeling the QBO and SAO Driven by Gravity Waves

Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW) is applied in a global scale numerical spectral model (NSM) to describe the semi-annual and quasi-biennial oscillations (SAO and QBO) as well as the long term interannual variations that are driven by wave mean flow interactions. This model has been successful in simulating the salient features observed near the equator at altitudes above 20 km, including the QBO extension into the upper mesosphere inferred from UARS measurements. The model has now been extended to describe also the mean zonal and meridional circulations of the upper troposphere and lower stratosphere that affect the equatorial QBO and its global scale extension. This is accomplished in part through tuning of the GW parameterization, and preliminary results lead to the following conclusions: (1) To reproduce the upwelling at equatorial latitudes associated with the Brewer/Dobson circulation that in part is modulated in the model by the vertical component of the Coriolis force, the eddy diffusivity in the lower stratosphere had to be enhanced and the related GW spectrum modified to bring it in closer agreement with the form recommended for the DSP. (2) To compensate for the required increase in the diffusivity, the observed QBO requires a larger GW source that is closer to the middle of the range recommended for the DSP. (3) Through global scale momentum redistribution, the above developments are conducive to extending the QBO and SAO oscillations to higher latitudes. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. (4) In a 3D version of the model, wave momentum is absorbed and dissipated by tides and planetary waves. Thus, a somewhat larger GW source is required to generate realistic amplitudes for the QBO and SAO.

Mayr, H. G.↗

Dynamic range considerations for EUV MAMA detectors

The multianode microchannel array (MAMA) has been chosen as the detector for two instruments on the ESA/NASA Solar Heliospheric Observatory. The response of the MAMA to the two extreme types of solar spectra, disk and corona, have been modeled with a view toward evaluating dynamic range effects present. The method of MAMA operation is discussed, with emphasis given to modeling the effect of electron cloud charge spreading to several detector anodes and amplifiers (n-fold events). Representative synthetic EUV spectra have been created. The detector response to these spectra is modeled by dissecting the input photon radiation field across the detector array into contributions to the various amplifier channels. The results of this dissection are shown for spectral regions across the entire wavelength region of interest. These results are used to identify regions in which total array photon counting rate or individual amplifier rate may exceed the design limits. This allows the design or operational modes to be tailored to eliminate the problem areas.

Illing, Rainer M. E.↗

The Long-term Transport and Radiative Impacts of the 2017 British Columbia Pyrocumulonimbus Smoke Aerosols in the Stratosphere

Interactions of meteorology with wildfires in British Columbia, Canada during August 2017 led to three major pyrocumulonimbus (pyroCb) events that resulted in the injection of large amounts of smoke aerosols and other combustion products at the local upper troposphere and lower stratosphere (UTLS). These plumes of UTLS smoke with elevated values of aerosol extinction and backscatter compared to the background state were readily tracked by multiple satellite-based instruments as they spread across the Northern Hemisphere (NH). The plumes were observed in the lower stratosphere for about 8-10 months following the fire injections, with a stratospheric aerosol e-folding time of about 5 months. To investigate the radiative impacts of these events on the Earth system, we performed a number of simulations with the Goddard Earth Observing System (GEOS) atmospheric general circulation model (AGCM). Observations from multiple remote-sensing instruments were used to calibrate the injection parameters (location, amount, composition and heights) and optical properties of the smoke aerosols in the model. The resulting simulations of three-dimensional smoke transport were evaluated for a year from the day of injections using daily observations from OMPS-LP (Ozone Mapping Profiler Suite Limb Profiler). The model simulated rate of ascent, hemispheric spread and residence time (or e-folding time) of the smoke aerosols in the stratosphere are in close agreement with OMPS-LP observations. We found that both aerosol self-lofting and the large-scale atmospheric motion play important roles in lifting the smoke plumes from near the tropopause altitudes (~12 km) to about 22-23 km into the atmosphere. Further, our estimations of the radiative impacts of the pyroCb-emitted smoke aerosols showed that the smoke caused an additional warming of the atmosphere by about 0.6-1 W/m2 (zonal mean) that persisted for about 2-3 months after the injections in regions north of 40oN. The surface experienced a comparable magnitude of cooling. The atmospheric warming is mainly located in the stratosphere, coincident with the location of the smoke plumes, leading to an increase in zonal mean shortwave (SW) heating rates of 0.02-0.04 K/day during September 2017.

pyrocumulonimbus (pyroCb) events↗

Non-Migrating Tides, with Zonally Symmetric Component, Generated in the Mesosphere

For comparison with measurements from the TIMED satellite and coordinated ground based observations, we discuss results from our Numerical Spectral Model (NSM) that incorporates the Doppler Spread Parameterization (Hines, 1997) for small-scale gravity waves (GWs). The NSM extends from the ground into the thermosphere and describes the major dynamical features of the atmosphere including the wave driven equatorial oscillations (QBO and SAO), and the seasonal variations of tides and planetary waves. With emphasis on the non-migrating tides, having periods of 24 and 12 hours, we discuss our modeling results that account for the classical migrating solar excitation sources only. As reported earlier, the NSM reproduces the observed seasonal variations and in particular the large equinoctial maxima in the amplitude of the migrating diurnal tide at altitudes around 90 km. Filtering of the tide by the zonal circulation and GW momentum deposition was identified as the cause. The GWs were also shown to produce a strong non-linear interaction between the diurnal and semi-diurnal tides. Confined largely to the mesosphere, the NSM produces through dynamical interactions a relatively large contribution of non-migrating tides. A striking feature is seen in the diurnal and semi-diurnal oscillations of the zonal mean (m = 0). Eastward propagating tides are also generated for zonal wave numbers m = 1 to 4. When the NSM is run without GWs, the amplitudes for the non-migrating tides, including m = 0, are generally small. Planetary wave interaction and non-linear coupling that involves the filtering of GWs and related height integration of dynamical features are discussed as possible mechanisms for generating these non-migrating tides in the NSM. As is the case for the solar migrating tides, the non-migrating tides reveal persistent seasonal variations. Under the influence of the QBO and SAO, interannual variations are produced.

Mayr, H. G.↗

A NON-PRE DOUBLE-PEAKED BURST FROM 4U 1636-536: EVIDENCE FOR BURNING FRONT PROPAGATION

We analyse Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) data of a double-peaked burst from the low mass X-ray binary (LMXB) 4U 1636-536 that shows no evidence for photospheric radius expansion (PRE). We find that the X-ray emitting area on the star increases with time as the burst progresses, even though the photosphere does not expand. We argue that this is a strong indication of thermonuclear flame spreading on the stellar surface during such bursts. We propose a model for such double-peaked bursts, based on thermonuclear flame spreading, that can qualitatively explain their essential features, as well as the rarity of these bursts.

Bhattacharyya, Sudip↗

Remote estimation of the surface characteristics and energy balance over an urban-rural area and the effects of surface heat flux on plume spread and concentration

A one dimensional boundary layer model was used in conjunction with satellite derived infrared surface temperatures to deduce values of moisture availability, thermal inertia, heat and evaporative fluxes. The Penn State satellite image display system, a sophisticated image display facility, was used to remotely sense these various parameters for three cases: St. Louis, Missouri; the Land Between the Lakes, Kentucky; and Clarksville, Tennessee. The urban centers displayed the maximum daytime surface temperatures which correspond to the minimum values of moisture availability. The urban center of St. Louis and the bodies of water displayed the maximum nighttime surface temperatures which correspond to the maximum thermal inertia values. It is shown that moisture availability and thermal inertia are very much responsible for the formation of important temperature variations over the urban rural complex.

Dicristofaro, D. C.↗

Modeling Study of Planetary Waves in the Mesosphere Lower Thermosphere (MLT)

For comparison with measurements from the TIMED satellite and coordinated ground based observations, we present results from our Numerical Spectral Model (NSM) that incorporates the Doppler Spread Parameterization (Hines, 1997) for small-scale gravity waves (GWs). We discuss the planetary waves (PWs) that are purely generated by dynamical interactions, i.e., without explicitly specifying excitation sources related for example to tropospheric convection or topography. With tropospheric heating that reproduces the observed zonal jets near the tropopause and the accompanying reversal in the latitudinal temperature variation, which is conducive to baroclinic instability, long period PWs are produced that propagate up into the stratosphere to affect the wave driven equatorial oscillations (QBO and SAO) extending into the upper mesosphere. The PWs in the model that dominate higher up in the MLT region, however, are to a large extent produced by instabilities under the influence of the zonal circulation and temperature variations in the middle atmosphere and they are amplified by GW interactions. Three classes of PWs are generated there. (1) Rossby waves that slowly propagate westward but are carried by the zonal mean (m = 0) winds to produce eastward and westward propagating PWs respectively in the winter and summer hemispheres below 80 km. Depending on the zonal wave number and magnitudes of the zonal winds under the influence of the equatorial oscillations, the PWs typically have periods between 2 and 20 days and their horizontal wind amplitudes can exceed 40 m/s in the lower mesosphere. (2) Rossby gravity waves that propagate westward at low latitudes, having periods around 2 days for zonal wave numbers m = 2 to 4. (3) Eastward propagating equatorial Kelvin waves generated in the upper mesosphere with periods between 2 and 3 days for m = 1 & 2. The seasonal variations of the PWs reveal that the largest wind amplitudes tend to occur below 80 km in the winter hemisphere, but above that altitude in the summer hemisphere to approach magnitudes as large as 50 m/s.

Mengel, J. G.↗

Computational predictions of flame spread over alcohol pools

The effects of buoyancy and thermocapillarity on pulsating and uniform flame spread above n-propanol fuel pools have been studied using a numerical model. Data obtained indicate that the existence of pulsating flame spread is dependent upon the formation of a gas-phase recirculation cell which entrains evaporating fuel vapor in front of the leading edge of the flame. The size of the recirculation cell which is affected by the extent of liquid motion ahead of the flame, is shown to dictate whether flame spread is uniform or pulsating. The amplitude and period of the flame pulsations are found to be proportional to the maximum extent of the flow head. Under conditions considered, liquid motion was not affected appreciably by buoyancy. Horizontal convection in the liquid is the dominant mechanism for transporting heat ahead of the flame for both the pulsating and uniform regimes.

Schiller, D. N.↗

Deterrent Concentration Measurement with FTIR and Subsequent Ballistic Performance in Medium Caliber Propellant

To meet ballistic requirements, medium and small caliber propellants use deterrent coatings to obtain burn rate progressivity. The required amount and distribution of deterrent varies between gun systems, propellant types, and often between lots. Micro Fourier Transform Infrared (FTIR) spectroscopy was used to measure deterrent gradients in RP36 propellants coated with methyl centralite (MC) at different deterrent levels and different processing conditions. The aromatic C-C bonds at 1496 cm(exp -1) wavenumber were used to monitor the deterrent profiles through the grain. Deterrent gradients measured with FTIR spectroscopy were then used to estimate burn rate gradients in the deterred grains. Burn rates were calculated from literature models and from closed bomb data of RP36 containing uniform deterrent concentration. Finally, the burn rate gradients were input into an IBHFG2 model of a 200 cc-closed bomb. The early flame spreading portion of the closed bomb ballistic cycle (0 to 0.2 P/Pmax) was roughly modeled by dividing the charge up into five propellant decks and igniting them at different times in the ballistic cycle. Pressure traces and vivacity curves from closed bomb shots were compared to predictions. In addition to the burn rate gradient, the closed bomb pressure trace was heavily dependent on ignition and flame spread. These two phenomena were not readily distinguishable from one another in deterred grains. The same RP-36 propellant was shot in a 25 mm M793TP round which was again modeled with IBHVG2. Peak pressure and muzzle velocity were accurately modeled when erosive burning effects were empirically factored into the model.

Furrow, Keith W.↗

Modeling Planetary Waves in the Middle Atmosphere

In our 3D Numerical Spectral Model (NSM) that incorporates the Doppler Spread Parameterization of Hines, planetary waves are generated with zonal wavenumbers 1-4 that have amplitudes as large as 25 m/s and periods between 50 and 2 days. The waves are generated spontaneously, i.e., without external source, and presumably represent normal modes. Numerical experiments are discussed to shed light on the generation mechanisms. For the wavenumbers 1 and 2 with periods primarily between 4 and 15 days, both eastward and westward propagating waves are generated and their amplitudes are affected significantly by the gravity wave source. For the wavenumbers, 3 and 4 in the upper mesosphere, the waves are virtually monochromatic, with periods around 2.5 days and they propagate westward. These waves behave like Rossby gravity waves; their amplitudes are largest during solstice where the baroclinic instability is playing an important role, but the waves decay when the gravity wave momentum source is turned off. The planetary waves generated in the model produce large amplitude modulations of the diurnal tides at altitudes above 80 km, but contrary to expectations they are not playing an important role for the seasonal variations of the tides. Preliminary results show that such Rossby gravity waves are also generated near the tropopause and are influencing significantly the QBO that is generated in the model.

Mengel, J. G.↗

Working group on chromospheric fields - Canopies

Although there are many points of uncertainty and controversy, the working group on chromospheric fields focussed its discussion on the concept of canopies; i.e., no one disagreed that a central issue relating to magnetic fields and chromospheric models is to learn how the photospheric field spreads with height. However, it quickly became apparent that in the time available, there was little prospect of building new unified models of magnetic field phenomena in the chromosphere beyond the scope of the formal presentations. Thus, the discussion was devoted to formulating questions which seemed both possible to address in future work and important for advancing understanding of the chromosphere. It began by discussing unresolved physical issues (almost everything) and then proceeded to consider means, both observational and synthetic, to address them.

Jones, H. P.↗

Supersonic Jet Noise Reductions Predicted With Increased Jet Spreading Rate

In this paper, predictions are made of noise radiation from single, supersonic, axisymmetric jets. We examine the effects of changes in operating conditions and the effects of simulated enhanced mixing that would increase the spreading rate of jet shear layer on radiated noise levels. The radiated noise in the downstream direction is dominated by mixing noise and, at higher speeds, it is well described by the instability wave noise radiation model. Further analysis with the model shows a relationship between changes in spreading rate due to enhanced mixing and changes in the far field radiated peak noise levels. The calculations predict that enhanced jet spreading results in a reduction of the radiated peak noise level.

Dahl, Milo D.↗

A source-sink model of the generation of plate tectonics from non-Newtonian mantle flow

A model of mantle convection which generates plate tectonics requires strain rate- or stress-dependent rheology in order to produce strong platelike flows with weak margins as well as strike-slip deformation and plate spin (i.e., toroidal motion). Here, we employ a simple model of source-sink driven surface flow to determine the form of such a rheology that is appropriate for Earth's present-day plate motions. In this model, lithospheric motion is treated as shallow layer flow driven by sources and sinks which correspond to spreading centers and subduction zones, respectively. Two plate motion models are used to derive the source sink field. As originally implied in the simpler Cartesian version of this model, the classical power law rheologies do not generate platelike flows as well as the hypothetical Whitehead-Gans stick-slip rheology (which incorporates a simple self-lubrication mechanism). None of the fluid rheologies examined, however, produce more than approximately 60% of the original maximum shear. For either plate model, the viscosity fields produced by the power law rheologies are diffuse, and the viscosity lows over strike-slip shear zones or pseudo-margins are not as small as over the prescribed convergent-divergent margins. In contrast, the stick-slip rheology generates very platelike viscosity fields, with sharp gradients at the plate boundaries, and margins with almost uniformly low viscosity. Power law rheologies with high viscosity contrasts, however, lead to almost equally favorable comparisons, though these also yield the least platelike viscosity fields. This implies that the magnitude of toroidal flow and platelike strength distributions are not necessarily related and thus may present independent constraints on the determination of a self-consistent plate-mantle rheology.

Bercovici, David↗

Applicability of the single equivalent point dipole model to represent a spatially distributed bio-electrical source

Although the single equivalent point dipole model has been used to represent well-localised bio-electrical sources, in realistic situations the source is distributed. Consequently, position estimates of point dipoles determined by inverse algorithms suffer from systematic error due to the non-exact applicability of the inverse model. In realistic situations, this systematic error cannot be avoided, a limitation that is independent of the complexity of the torso model used. This study quantitatively investigates the intrinsic limitations in the assignment of a location to the equivalent dipole due to distributed electrical source. To simulate arrhythmic activity in the heart, a model of a wave of depolarisation spreading from a focal source over the surface of a spherical shell is used. The activity is represented by a sequence of concentric belt sources (obtained by slicing the shell with a sequence of parallel plane pairs), with constant dipole moment per unit length (circumferentially) directed parallel to the propagation direction. The distributed source is represented by N dipoles at equal arc lengths along the belt. The sum of the dipole potentials is calculated at predefined electrode locations. The inverse problem involves finding a single equivalent point dipole that best reproduces the electrode potentials due to the distributed source. The inverse problem is implemented by minimising the chi2 per degree of freedom. It is found that the trajectory traced by the equivalent dipole is sensitive to the location of the spherical shell relative to the fixed electrodes. It is shown that this trajectory does not coincide with the sequence of geometrical centres of the consecutive belt sources. For distributed sources within a bounded spherical medium, displaced from the sphere's centre by 40% of the sphere's radius, it is found that the error in the equivalent dipole location varies from 3 to 20% for sources with size between 5 and 50% of the sphere's radius. Finally, a method is devised to obtain the size of the distributed source during the cardiac cycle.

NASA Discipline Cardiopulmonary↗

Numerical modelling of instantaneous plate tectonics

Assuming lithospheric plates to be rigid, 68 spreading rates, 62 fracture zones trends, and 106 earthquake slip vectors are systematically inverted to obtain a self-consistent model of instantaneous relative motions for eleven major plates. The inverse problem is linearized and solved iteratively by a maximum-likelihood procedure. Because the uncertainties in the data are small, Gaussian statistics are shown to be adequate. The use of a linear theory permits (1) the calculation of the uncertainties in the various angular velocity vectors caused by uncertainties in the data, and (2) quantitative examination of the distribution of information within the data set. The existence of a self-consistent model satisfying all the data is strong justification of the rigid plate assumption. Slow movement between North and South America is shown to be resolvable.

Minster, J. B.↗

Comparison of two gas chromatograph models and analysis of binary data

The overall objective of the gas chromatograph system studies is to generate fundamental design criteria and techniques to be used in the optimum design of the system. The particular tasks currently being undertaken are the comparison of two mathematical models of the chromatograph and the analysis of binary system data. The predictions of two mathematical models, an equilibrium absorption model and a non-equilibrium absorption model exhibit the same weaknesses in their inability to predict chromatogram spreading for certain systems. The analysis of binary data using the equilibrium absorption model confirms that, for the systems considered, superposition of predicted single component behaviors is a first order representation of actual binary data. Composition effects produce non-idealities which limit the rigorous validity of superposition.

Keba, P. S.↗

Steepened structures in equatorial spread F. 2: Theory

This paper investigates the properties of a one-dimensional fluid model of plasma convection in the equatorial F region ionosphere. The model equations are similar in form to Burgers equation except for additional higher-order spatial derivatives. Like Burgers equation, solution to the model have the form of propagating, shocklike structures. Numerical simulations of the model closely resemble the steepened structures observed by sounding rocket plasma density probes within equatorial spread F. Simulated denstiy power spectra, like the spectra computed from in situ data, seem to possess power law forms with a break at wavelengths of about 100 m. The precise wavenumber of the spectral break is determined by the ambipolar diffusion coefficient. The model predicts that electric field fluctuations perpendicular to the direction of plasma steepening should be proportional to the plasma density fluctuations. Electric field fluctuations parallel to the steepening will be due primarily to the ambipolar field and have a Boltzmann relationship with density (square of the absolute value of delta E) approximately equal to (K(exp 2))(square of the absolute value of (delta n/n)). At wavelengths less than about 300 m, the ambipolar field should be the dominant component of the total field intensity.

Hysell, D. L↗

Observations of the Saturn D ring

Observations of Saturn at both large and small phase angles were obtained in the 890-nm methane band with a solid state charge coupled device (CCD). When model intensity profiles are convolved with appropriate point spread functions, intensity levels ranging from 5 to 3% of the maximum of ring B are required in the ring D region of the models to match the data, in substantial agreement with previous results obtained through photographic photometry. The normal optical thickness of the D ring is about 0.02 and seems to show the same phase function as the other rings. The division between the C and D ring was not observed, possibly due to the seeing limited image resolution.

Larson, S. M.↗