Engineering PapersSearch

SEARCH · Engineering Papers

Results for “total variation regularization”

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 37 records · Page 2

Solar irradiance variations and solar activity

A mean value for the 1 AU total solar irradiance of 1368.2 W/sq m and a downward trend of 0.05% per year were derived from measurements by the Active Cavity Radiometer Irradiance Monitor (ACRIM) experiment on the Solar Maximum Mission during 1980. Distinct temporary solar irradiance decreases associated with solar activity maxima were observed with a series of nine dips from April to October recurring at fairly regular intervals averaging 24 days. The decreases correlate inversely with sunspot area, 2800-MHz flux, and Zurich sunspot number. Dominant periods common to the irradiance and sunspot area power spectra link the irradiance decreases to sunspot flux deficit in solar active regions. Evidence of significant total irradiance modulation by facular flux excess is cited. A persistent radiative cycle of active regions consistent with the ACRIM irradiance results and the morphology of solar active regions was found. The pattern of regularly recurrent active region maxima between April and October suggests an asymmetry in solar activity generation during this period.

Willson, R. C.

Cyclic Oxidation Modeling and Life Prediction

The cyclic oxidation process can be described as an iterative scale growth and spallation sequence by a number of similar models. Model input variable include oxide scale type and growth parameters, spalling geometry, spall constant, and cycle duration. Outputs include net weight change, the amounts of retained and spalled oxide, the total oxygen and metal consumed, and the terminal rates of weight loss and metal consumption. All models and their variations produce a number of similar characteristic features. In general, spalling and material consumption increase to a steady state rate, at which point the retained scale approaches a constant and the rate of weight loss becomes linear. For one model, this regularity was demonstrated as dimensionless, universal expressions, obtained by normalizing the variables by critical performance factors. These insights were enabled through the use of the COSP for Windows cyclic oxidation spalling program.

Smialek, James L.

Postseismic Deformation after the 1964 Great Alaskan Earthquake: Collaborative Research with Goddard Space Flight Center

The purpose of this project was to carry out GPS observations on the Kenai Peninsula, southern Alaska, in order to study the postseismic and contemporary deformation following the 1964 Alaska earthquake. All of the research supported in this grant was carried out in collaboration with Dr. Steven Cohen of Goddard Space Flight Center. The research funding from this grant primarily supported GPS fieldwork, along with the acquisition of computer equipment to allow analysis and modeling of the GPS data. A minor amount of salary support was provided by the PI, but the great majority of the salary support was provided by the Geophysical Institute. After the expiration of this grant, additional funding was obtained from the National Science Foundation to continue the work. This grant supported GPS field campaigns in August 1995, June 1996, May-June and September 1997, and May-June 1998. We initially began the work by surveying leveling benchmarks on the Kenai peninsula that had been surveyed after the 1964 earthquake. Changes in height from the 1964 leveling data to the 1995+ GPS data, corrected for the geoid-ellipsoid separation, give the total elevation change since the earthquake. Beginning in 1995, we also identified or established sites that were suitable for long-term surveying using GPS. In the subsequent annual GPS campaigns, we made regular measurements at these GPS marks, and steadily enhanced our set of points for which cumulative postseismic uplift data were available. From 4 years of Global Positioning System (GPS) measurements, we find significant spatial variations in present-day deformation between the eastern and western Kenai peninsula, Alaska. Sites in the eastern Kenai peninsula and Prince William Sound move to the NNW relative to North America, in the direction of Pacific-North America relative plate motion. Velocities decrease in magnitude from nearly the full plate rate in southern Prince William Sound to about 30 mm/yr at Seward and to about 5 mm/yr near Anchorage. In contrast, sites in the western Kenai peninsula move to the SW, in a nearly trenchward direction, with a velocity of about 20 mm/yr. The data are consistent with the shallow plate interface offshore and beneath the eastern Kenai and Prince William Sound being completely locked or nearly so, with elastic strain accumulation resulting in rapid motion in the direction of relative plate motion of sites in the overriding plate. The velocities of sites in the western Kenai, along strike to the southwest, are opposite in sign with those predicted from elastic strain accumulation. These data are incompatible with a significant locked region in this segment of the plate boundary. Trenchward velocities are found also for some sites in the Anchorage area. We interpret the trenchward velocities as being caused by a continuing postseismic transient from the 1964 great Alaska earthquake.

Freymueller, Jeffrey T.

Data traffic reduction schemes for Cholesky factorization on asynchronous multiprocessor systems

Communication requirements of Cholesky factorization of dense and sparse symmetric, positive definite matrices are analyzed. The communication requirement is characterized by the data traffic generated on multiprocessor systems with local and shared memory. Lower bound proofs are given to show that when the load is uniformly distributed the data traffic associated with factoring an n x n dense matrix using n to the alpha power (alpha less than or equal 2) processors is omega(n to the 2 + alpha/2 power). For n x n sparse matrices representing a square root of n x square root of n regular grid graph the data traffic is shown to be omega(n to the 1 + alpha/2 power), alpha less than or equal 1. Partitioning schemes that are variations of block assignment scheme are described and it is shown that the data traffic generated by these schemes are asymptotically optimal. The schemes allow efficient use of up to O(n to the 2nd power) processors in the dense case and up to O(n) processors in the sparse case before the total data traffic reaches the maximum value of O(n to the 3rd power) and O(n to the 3/2 power), respectively. It is shown that the block based partitioning schemes allow a better utilization of the data accessed from shared memory and thus reduce the data traffic than those based on column-wise wrap around assignment schemes.

Naik, Vijay K.

Dwarf novae

Dwarf novae are defined on grounds of their semi-regular brightness variations of some two to five magnitudes on time scales of typically 10 to 100 days. Historically several different classification schemes have been used. Today, dwarf novae are divided into three sub-classes: the U Geminorum stars, the SU Ursae Majoris stars, and the Z Camelopardalis stars. Outbursts of dwarf novae occur at semi-periodic intervals of time, typically every 10 to 100 days; amplitudes range from typically 2 to 5 mag. Within certain limits values are characteristic for each object. Relations between the outburst amplitude, or the total energy released during outburst, and the recurrence time have been found, as well as relations between the orbital period and the outburst decay time, the absolute magnitude during outburst maximum, and the widths of long and short outbursts, respectively. Some dwarf novae are known to have suspended their normal outburst activity altogether for a while. They later resumed it without having undergone any observable changes. The optical colors of dwarf novae all are quite similar during outburst, considerably bluer than during the quiescent state. During the outburst cycle, characteristic loops in the two color diagram are performed. At a time resolution on the order of minutes, strictly periodic photometric changes due to orbital motion become visible in the light curves of dwarf novae. These are characteristic for each system. Remarkably little is known about orbital variations during the course of an outburst. On time-scales of minutes and seconds, further more or less periodic types of variability are seen in dwarf novae. Appreciable flux is emitted by dwarf novae at all wavelengths from the X-rays to the longest IR wavelengths, and in some cases even in the radio. Most dwarf novae exhibit strong emission line spectra in the optical and UV during quiescence, although some have only very weak emissions in the optical and/or weak absorptions at UV wavelengths.

Ladous, Constanze

Elimination of secular terms from the differential equations for the elements of perturbed two-body motion

In 1961, Sperling linearized and regularized the differential equations of motion of the two-body problem by changing the independent variable from time to fictitious time by Sundman's transformation (r = dt/ds) and by embedding the two-body energy integral and the Laplace vector. In 1968, Burdet developed a perturbation theory which was uniformly valid for all types of orbits using a variation of parameters approach on the elements which appeared in Sperling's equations for the two-body solution. In 1973, Bond and Hanssen improved Burdet's set of differential equations by embedding the total energy (which is a constant when the potential function is explicitly dependent upon time.) The Jacobian constant was used as an element to replace the total energy in a reformulation of the differential equations of motion. In the process, another element which is proportional to a component of the angular momentum was introduced. Recently trajectories computed during numerical studies of atmospheric entry from circular orbits and low thrust beginning in near-circular orbits exhibited numerical instability when solved by the method of Bond and Gottlieb (1989) for long time intervals. It was found that this instability was due to secular terms which appear on the righthand sides of the differential equations of some of the elements. In this paper, this instability is removed by the introduction of another vector integral called the delta integral (which replaces the Laplace Vector) and another scalar integral which removes the secular terms. The introduction of these integrals requires a new derivation of the differential equations for most of the elements. For this rederivation, the Lagrange method of variation of parameters is used, making the development more concise. Numerical examples of this improvement are presented.

Bond, Victor R.

Empirical correction for earth sensor horizon radiance variation

A major limitation on the use of infrared horizon sensors for attitude determination is the variability of the height of the infrared Earth horizon. This variation includes a climatological component and a stochastic component of approximately equal importance. The climatological component shows regular variation with season and latitude. Models based on historical measurements have been used to compensate for these systematic changes. The stochastic component is analogous to tropospheric weather. It can cause extreme, localized changes that for a period of days, overwhelm the climatological variation. An algorithm has been developed to compensate partially for the climatological variation of horizon height and at least to mitigate the stochastic variation. This method uses attitude and horizon sensor data from spacecraft to update a horizon height history as a function of latitude. For spacecraft that depend on horizon sensors for their attitudes (such as the Total Ozone Mapping Spectrometer-Earth Probe-TOMS-EP) a batch least squares attitude determination system is used. It is assumed that minimizing the average sensor residual throughout a full orbit of data results in attitudes that are nearly independent of local horizon height variations. The method depends on the additional assumption that the mean horizon height over all latitudes is approximately independent of season. Using these assumptions, the method yields the latitude dependent portion of local horizon height variations. This paper describes the algorithm used to generate an empirical horizon height. Ideally, an international horizon height database could be established that would rapidly merge data from various spacecraft to provide timely corrections that could be used by all.

Hashmall, Joseph A.

Measurement of total electron content of midlatitude ionosphere and protonosphere via Faraday rotation and group relay techniques using transmission from geostationary satellites ATS-3 and ATS-6

Measurement of integrated columnar electron content and total electron content for the local ionosphere and the overlying protonosphere via Faraday rotation and group delay techniques has proven very useful. A field station was established having the geographic location of 31.5 deg N latitude and 91.06 deg W longitude to accomplish these objectives. A polarimeter receiving system was set up in the beginning to measure the Faraday rotation of 137.35 MHz radio signal from geostationary satellite ATS 3 to yield the integrated columnar electron content of the local ionosphere. The measurement was continued regularly, and the analysis of the data thus collected provided a synopsis of the statistical variation of the ionosphere along with the transient variations that occurred during the periods of geomagnetic and other disturbances.

Paul, M. P.

Steps toward determination of the size and structure of the broad-line region in active galactic nuclei. 5: Variability of the ultraviolet continuum and emission lines of NGC 3783

We report on the results of intensive ultraviolet spectral monitoring of the Seyfert 1 galaxy NGC 3783. The nucleus of NGC 3783 was observed with the International Ultraviolet Explorer satellite on a regular basis for a total of 7 months, once every 4 days for the first 172 days and once every other day for the final 50 days. Significant variability was observed in both continuum and emission-line fluxes. The light curves for the continuum fluxes exhibited two well-defined local minima or 'dips,' the first lasting is less than or approximately 20 days and the second is less than or approximately 4 days, with additional episodes of relatively rapid flickering of approximately the same amplitude. As in the case of NGC 5548 (the only other Seyfert galaxy that has been the subject of such an intensive, sustained monitoring effort), the largest continuum variations were seen at the shortest wavelengths, so that the continuum became 'harder' when brighter. The variations in the continuum occurred simultaneously at all wavelengths (delta(t) is less than 2 days). Generally, the amplitude of variability of the emission lines was lower than (or comparable to) that of the continuum. Apart from Mg II (which varied little) and N V (which is relatively weak and badly blended with Ly(alpha), the light curves of the emission lines are very similar to the continuum light curves, in each case with a small systematic delay or 'lag.' As for NGC 5548, the highest ionization lines seem to respond with shorter lags than the lower ionization lines. The lags found for NGC 3783 are considerably shorter than those obtained for NGC 5548, with values of (formally) approximately 0 days for He II + O III), and approximately 4 days for Ly(alpha) and C IV. The data further suggest lags of approximately 4 days for Si IV + O IV) and 8-30 days for Si III + C III). Mg II lagged the 1460 A continuum by approximately 9 days, although this result depends on the method of measuring the line flux and may in fact be due to variability of the underlying Fe II lines. Correlation analysis further shows that the power density spectrum contains substantial unresolved power over timescales of is less than or approximately 2 days, and that the character of the continuum variability may change with time.

Reichert, G. A.

IUE observations of Comet Halley: Evolution of the UV spectrum between September 1985 and July 1986

The ultraviolet spectrum of comet P/Halley was monitored with the IUE between 12 September 1985 and 8 July 1986 (r <2.6 AU pre and post-perihelion) at regular time intervals except for a two-month period around the time of perihelion. A complete characterization of the UV spectrum of the comet was obtained to derive coma abundances and to study the light emission mechanisms of the observed species. The Fine Error Sensor (FES) camera of the IUE was used to photometrically investigate the coma brightness variation on time scales of the order of hours. Spectroscopic observations as well as FES measurements show that the activity of the nucleus is highly variable, particularly at the end of December 1985 and during March and April 1986. The production rates of OH, CS and dust are derived for the entire period of the observations. The total water loss rate for this period is estimated to be 150 million metric tons.

Feldman, P. D.

Estimating When the Antarctic Ozone Hole Will Recover

The Antarctic ozone hole develops each year and culminates by early spring (late September - early October). The severity of the hole has been assessed from satellites using the minimum total ozone value from the October monthly mean (depth of the hole) and by calculating the average area coverage during this September-October period. Profile information shows that ozone is completely destroyed in the 14-21 km layer by early October. Ozone is mainly destroyed by halogen (chlorine and bromine) catalytic cycles, and these losses are modulated by temperature variations. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. Using projections of halogen levels combined with age-of-air estimates, we find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. We estimate that the ozone hole will begin to show first signs of size decrease in about 2023, and the hole will fully recover to pre-1980 levels in approximately 2070. Estimates of the ozone hole's recovery from models reveal important differences that will be discussed.

Newman, Paul A.

Estimating when the Antarctic Ozone Hole will Recover

The Antarctic ozone hole develops each year and culminates by early spring (late September - early October). The severity of the hole has been assessed from satellites using the minimum total ozone value from the October monthly mean (depth of the hole) and by calculating the average area coverage during this September-October period. Profile information shows that ozone is completely destroyed in the 14-2 1 km layer by early October. Ozone is mainly destroyed by halogen (chlorine and bromine) catalytic cycles, and these losses are modulated by temperature variations. Because atmospheric halogen levels are responding to international a'greements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. Using projections of halogen levels combined with age-of-air estimates, we find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. We estimate that the ozone hole will begin to show first signs of size decrease in about 2023, and the hole will fully recover to pre-1980 levels in approximately 2070. Estimates of the ozone hole's recovery from models reveal important differences that will be discussed.

Newman, Paul A.

Some Numerical Experiments on Detonation Wave Propagation

In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Cambier, Jean-Luc

Detecting the Recovery of the Antarctic Ozone Hole

The Antarctic ozone hole develops each year and culminates by early Spring. Antarctic ozone values have been monitored since 1979 using satellite observations from the TOMS instrument. The severity of the hole has been assessed from TOMS using the minimum total ozone value from the October monthly mean (depth of the hole) and by calculating the average size during the September-October period. Ozone is mainly destroyed by halogen catalytic cycles, and these losses are modulated by temperature variations in the collar of the polar lower stratospheric vortex. In this presentation, we show the relationships of halogens and temperature to both the size and depth of the hole. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. Using projections of halogen levels combined with age-of-air estimates, we find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. We will show estimates of both when the ozone hole will begin to show first signs of recovery, and when the hole will fully recover to pre-1980 levels.

Newman, Paul A.

When Will the Antarctic Ozone Hole Recover?

The Antarctic ozone hole develops each year and culminates by early Spring. Antarctic ozone values have been monitored since 1979 using satellite observations from the TOMS instrument. The severity of the hole has been assessed from TOMS using the minimum total ozone value from the October monthly mean (depth of the hole) and by calculating the average size during the September-October period. Ozone is mainly destroyed by halogen catalytic cycles, and these losses are modulated by temperature variations in the collar of the polar lower stratospheric vortex. In this presentation, we show the relationships of halogens and temperature to both the size and depth of the hole. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. Using projections of halogen levels combined with age-of-air estimates, we find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. We will show estimates of both when the ozone hole will begin to show first signs of recovery, and when the hole will fully recover to pre-1980 levels.

Newman, Paul A.

When will the Antarctic Ozone Hole Recover?

The Antarctic ozone hole develops each year and culminates by early Spring. Antarctic ozone values have been monitored since 1979 using satellite observations from the .TOMS instrument. The severity of the hole has been assessed from TOMS using the minimum total ozone value from the October monthly mean (depth of the hole) and by calculating the average size during the September-October period. Ozone is mainly destroyed by halogen catalytic cycles, and these losses are modulated by temperature variations in the collar of the polar lower stratospheric vortex. In this presentation, we show the relationships of halogens and temperature to, both the size and depth of the hole. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. Using projections of halogen levels combined with age-of-air estimates, we find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. The ozone hole will begin to show first signs of recovery in about 2023, and the hole will fully recover to pre-1980 levels in approximately 2070. This 2070 recovery is 20 years later than recent projections.

Newman, Paul A.

When Will the Antarctic Ozone Hole Recover?

The Antarctic ozone hole develops each year and culminates by early spring (late September - early October). Antarctic ozone values have been monitored since 1979 using satellite observations from the TOMS instrument. The severity of the hole has been assessed from TOMS using the minimum total ozone value from the October monthly mean (depth of the hole) and by calculating the average area coverage during this September-October period. Ozone is mainly destroyed by halogen (chlorine and bromine) catalytic cycles, and these losses are modulated by temperature variations in the collar of the polar lower stratospheric vortex. In this talk, I will show the relationships of halogens and temperature to both the size and depth of the hole. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. Using projections of halogen levels combined with age-of-air estimates, we find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. The ozone hole will begin to show first signs of recovery in about 2023, and the hole will fully recover to pre-1980 levels in approximately 2070. This 2070 recovery is 20 years later than recent projections. I will also discuss current assessments of mid-latitude ozone recovery.

Newman, Paul A.

The Global Methane Budget 2000–2017

Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric lifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). For the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 Tg CH4/yr (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 Tg CH4/yr or ∼ 60 % is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 Tg CH4/yr or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 Tg CH4/yr larger than our estimate for the previous decade (2000–2009), and 24 Tg CH4/yr larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30 % larger global emissions (737 Tg CH4/yr, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼ 65 % of the global budget, < 30° N) compared to mid-latitudes (∼ 30 %, 30–60° N) and high northern latitudes (∼ 4 %, 60–90° N). The most important source of uncertainty in the methane budget is attributable to natural emissions, especially those from wetlands and other inland waters. Some of our global source estimates are smaller than those in previously published budgets (Saunois et al., 2016; Kirschke et al., 2013). In particular wetland emissions are about 35 Tg CH4/yr lower due to improved partition wetlands and other inland waters. Emissions from geological sources and wild animals are also found to be smaller by 7 Tg CH4/yr by 8 Tg CH4/yr, respectively. However, the overall discrepancy between bottom-up and top-down estimates has been reduced by only 5 % compared to Saunois et al. (2016), due to a higher estimate of emissions from inland waters, highlighting the need for more detailed research on emissions factors. Priorities for improving the methane budget include (i) a global, high-resolution map of water-saturated soils and inundated areas emitting methane based on a robust classification of different types of emitting habitats; (ii) further development of process-based models for inland-water emissions; (iii) intensification of methane observations at local scales (e.g., FLUXNET-CH4 measurements) and urban-scale monitoring to constrain bottom-up land surface models, and at regional scales (surface networks and satellites) to constrain atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) development of a 3D variational inversion system using isotopic and/or co-emitted species such as ethane to improve source partitioning.

methane budget