Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized gradient approximations”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Grains in ionized nebulae: Spectral line diagnostics

The depletion of condensable elements onto grains in gaseous nebulae can provide evidence that dust is well mixed with the ionized gas. Al and CA are two of the most depleted elements in the general interstellar medium, and it is therefore important to measure their abundances within the ionized region of nebulae. We compute a large grid of photoionization models and identify sets of line ratios which are relatively insensitive to stellar and nebular parameters, and are thus excellent diagnostics for determining relative abundances. Based on the absence of the ((Ca II) lambda lambda 291, 7324 doublet and the detection of Al II) lambda lambda 2660, 2669 in the ultraviolet, we determine the extent of aluminum and calcium depletion onto grains in NGC 7027 and the Orion Nebula. Our results show a approximately 0.3 dex depletion for Al, but a depletion of more than two and a half orders of magnitude for Ca. A similar calculation based on Mg II lambda 2798 yields roughly a 0.8 dex depletion for Mg. This reaffirms the discrepancy between depletion determined from high and low ionization Mg lines. We also find evidence for a 'depletion gradient' in Ca in NGC 7027, since the calcium depletion we infer for the outer, more neutral regions using (Ca II) is somewhat higher than that inferred for the inner high-ionization region, using (Ca v). This gradient can test current models of the survival of grains within hot ionized gas.

Kingdon, J.↗
NASA NTRS1997

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗
NASA NTRS1993

Advances in Studies of Electrode Kinetics and Mass Transport in AMTEC Cells (abstract)

Previous work reported from JPL has included characterization of electrode kinetics and alkali atom transport from electrodes including Mo, W, WRh(sub x), WPt(sub x)(Mn), in sodium AMTEC cells and vapor exposure cells, and Mo in potassium vapor exposure cells. These studies were generally performed in cells with small area electrodes (about 1 to 5 cm(sup 2)), and device geometry had little effect on transport. Alkali diffusion coefficients through these electrodes have been characterized, and approximate surface diffusion coefficients derived in cases of activated transport. A basic model of electrode kinetic at the alkali metal vapor/porous metal electrode/alkali beta'-alumina solid electrolyte three phase boundary has been proposed which accounts for electrochemical reaction rates with a collision frequency near the three phase boundary and tunneling from the porous electrode partially covered with adsorbed alkali metal atoms. The small electrode effect in AMTEC cells has been discussed in several papers, but quantitative investigations have described only the overall effect and the important contribution of electrolyte resistance. The quantitative characterization of transport losses in cells with large area electrodes has been limited to simulations of large area electrode effects, or characterization of transport losses from large area electrodes with significant longitudinal temperature gradients. This paper describes new investigations of electrochemical kinetics and transport, particularily with WPt(sub 3.5) electrodes, including the influence of electrode size on the mass transport loss in the AMTEC cell. These electrodes possess excellent sodium transport properties making verification of device limitations on transport much more readily attained.

electrode kinetics alkali atom transport transport↗
NASA NTRS1993

Advances in Studies of Electrode Kinetics and Mass Transport in AMTEC Cells

Previous work reported from JPL has included characterization of electrode kinetics and alkali atom transport from electrodes including Mo, W, WRh x (Mn), in sodium AMTEC cells and vapor exposure cells; and Mo in potassium vapor exposure cells. These studies were generally performed in cells with small area electrodes (about 1 to 5 cm 2 ), and device geometry had little effect on transport. Alkali metal diffusion coefficients through these electrodes have been characterized, and approximate surface diffusion coefficients were derived in cases of activated transport. A basic model of electrode kinetics and transport at the alkali metal vapor/pourous metal electrode/alkali beta"-alumina solid electrolyte (BASE) three phase boundary has been proposed which accounts for electrochemical reaction rates with a collision frequency near the three phase boundary and tunneling from the porous electrode partially covered with adsorbed alkali metal atoms. The small electrode effect in AMTEC cells has been discussed in several papers, but quantitative investigations have described only the overall effect and the important contribution of electrolyte resistance. The quantitative characterization of transport losses in cells with large area electrodes has been limited to simulations of large area electrode effects, or characterization of transport losses from large area electrodes with significant longitudinal temperature gradients. This paper describes new investigations of electrochemical kinetics and transport, with four 14.4 cm 2 WPt 3.5 electrodes, including the influence of electrode size on the mass transport loss in the AMTEC cell. These electrodes exhibit very slow sintering, as well as excellent sodium transport properties, making them attractive candidates for AMTEC power conversion use. However, the facile sodium transport in WPt 3.5 electrodes makes characterization of the transport process difficult.

electrode↗
NASA NTRS1994

Large-scale characteristics of interstellar dust from COBE DIRBE observations

Observations from the COBE Diffuse Infrared Background Experiment of the 140 and 240 micrometer emissions from the Galatic plane region (absolute value of b less than 10 deg) are combined with radio surveys that trace the molecular (H2), neutral atomic (H I), and extended low-density (n(sub e) approximately 10 to 100/cm(exp 3)) ionized (H II) gas phases of the interstellar medium to derive physical conditions such as the dust temperature, dust-to-gas mass ratio, and far-infrared emissivity (1) averaged over these gas phases along each line of sight and (2) within each of these three gas phases. This analysis shows large-scale longitudinal and latitudinal gradients in the dust temperature and a decrease in dust temperature with increasing Galactocentric distance. The derived dust temperatures are significantly different from those derived in similar analyses using the Infrared Astronomical Satellite (IRAS) 60 and 100 micrometer data, suggesting that small (5 A approximately less than radius approximately less than 200 A) transiently heated dust particles contribute significantly o the Galactic 60 micrometer emission. It is found that 60% to 75% of the far-infrared luminosity arises from cold (approximately 17 to 22 K) dust associated with diffuse H I clouds, 15% to 30% from cold (approximately 19 K) dust associated with molecular gas, and less than 10% from warm (approximately 29 K) dust in extended low-density H II regions, consistent with the results of the IRAS analyses of the Galactic 60 and 100 micrometer emission. Within 2 deg of longitude of the Galactic center, the derived gas-to-dust mass ratio along the line of sight, G(sub d), reverses its general trend of decreasing G(sub d) toward the inner Galaxy and increases by a factor of approximately 2 to 3 toward the Galactic center. One possible explanation for this result is that the ratio of H2 column density to (12)CO intensity is lower in the Galactic center region than in the Galactic disk.

Sodroski, T. J.↗
NASA NTRS2013

RF Reference Switch for Spaceflight Radiometer Calibration

The goal of this technology is to provide improved calibration and measurement sensitivity to the Soil Moisture Active Passive Mission (SMAP) radiometer. While RF switches have been used in the past to calibrate microwave radiometers, the switch used on SMAP employs several techniques uniquely tailored to the instrument requirements and passive remote-sensing in general to improve radiometer performance. Measurement error and sensitivity are improved by employing techniques to reduce thermal gradients within the device, reduce insertion loss during antenna observations, increase insertion loss temporal stability, and increase rejection of radar and RFI (radio-frequency interference) signals during calibration. The two legs of the single-pole double-throw reference switch employ three PIN diodes per leg in a parallel-shunt configuration to minimize insertion loss and increase stability while exceeding rejection requirements at 1,413 MHz. The high-speed packaged diodes are selected to minimize junction capacitance and resistance while ensuring the parallel devices have very similar I-V curves. Switch rejection is improved by adding high-impedance quarter-wave tapers before and after the diodes, along with replacing the ground via of one diode per leg with an open circuit stub. Errors due to thermal gradients in the switch are reduced by embedding the 50-ohm reference load within the switch, along with using a 0.25-in. (approximately equal to 0.6-cm) aluminum prebacked substrate. Previous spaceflight microwave radiometers did not embed the reference load and thermocouple directly within the calibration switch. In doing so, the SMAP switch reduces error caused by thermal gradients between the load and switch. Thermal issues are further reduced by moving the custom, highspeed regulated driver circuit to a physically separate PWB (printed wiring board). Regarding RF performance, previous spaceflight reference switches have not employed high-impedance tapers to improve rejection. The use of open-circuit stubs instead of a via to provide an improved RF short is unique to this design. The stubs are easily tunable to provide high rejection at specific frequencies while maintaining very low insertion loss in-band.

Knuble, Joseph↗
DOE OSTI2024

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2015

Solar Signals in CMIP-5 Simulations: The Stratospheric Pathway

The 11 year solar-cycle component of climate variability is assessed in historical simulations of models taken from the Coupled Model Intercomparison Project, phase 5 (CMIP-5). Multiple linear regression is applied to estimate the zonal temperature, wind and annular mode responses to a typical solar cycle, with a focus on both the stratosphere and the stratospheric influence on the surface over the period approximately 1850-2005. The analysis is performed on all CMIP-5 models but focuses on the 13 CMIP-5 models that resolve the stratosphere (high-top models) and compares the simulated solar cycle signature with reanalysis data. The 11 year solar cycle component of climate variability is found to be weaker in terms of magnitude and latitudinal gradient around the stratopause in the models than in the reanalysis. The peak in temperature in the lower equatorial stratosphere (approximately 70 hPa) reported in some studies is found in the models to depend on the length of the analysis period, with the last 30 years yielding the strongest response. A modification of the Polar Jet Oscillation (PJO) in response to the 11 year solar cycle is not robust across all models, but is more apparent in models with high spectral resolution in the short-wave region. The PJO evolution is slower in these models, leading to a stronger response during February, whereas observations indicate it to be weaker. In early winter, the magnitude of the modeled response is more consistent with observations when only data from 1979-2005 are considered. The observed North Pacific high-pressure surface response during the solar maximum is only simulated in some models, for which there are no distinguishing model characteristics. The lagged North Atlantic surface response is reproduced in both high- and low-top models, but is more prevalent in the former. In both cases, the magnitude of the response is generally lower than in observations.

stratosphere↗
NASA NTRS1943

Heat Transfer Through Turbulent Friction Layers

The "general Prandtl number" Pr(exp 1) - A(sub q)/A Pr, aside from the Reynolds number determines the ratio of turbulent to molecular heat transfer, and the temperature distribution in turbulent friction layers. A(sub q) = exchange coefficient for heat; A = exchange coefficient for momentum transfer. A formula is derived from the equation defining the general Prandtl number which describes the temperature as a function of the velocity. For fully developed thermal boundary layers all questions relating to heat transfer to and from incompressible fluids can be treated in a simple manner if the ratio of the turbulent shear stress to the total stress T(sub t)/T in the layers near the wall is known, and if the A(sub q)/A can be regarded as independent of the distance from the wall. The velocity distribution across a flat smooth channel and deep into the laminar sublayer was measured for isothermal flow to establish the shear stress ratio T(sub t)/T and to extend the universal wall friction law. The values of T(sub t)/T which resulted from these measurements can be approximately represented by a linear function of the velocity in the laminar-turbulent transition zone. The effect of the temperature relationship of the material values on the flow near the wall is briefly analyzed. It was found that the velocity at the laminar boundary (in contrast to the thickness of the laminar layer) is approximately independent of the temperature distribution. The temperature gradient at the wall and the distribution of temperature and heat flow in the turbulent friction layers were calculated on the basis of the data under two equations. The derived formulas and the figures reveal the effects of the Prandtl number, the Reynolds number, the exchange quantities and the temperature relationship of the material values.

Reichardt, H.↗
NASA NTRS1993

Energy transfer and constrained simulations in isotropic turbulence

The defining characteristic of turbulent flows is their ability to dissipate energy, even in the limit of zero viscosity. The Euler equations, if constrained in such a way that the velocity derivatives remain bounded, conserve energy. But when they arise as the limit of the Navier-Stokes (NS) equations, when the Reynolds number goes to infinity, there is persuasive empirical evidence that the gradients become singular as just the right function of Re for the dissipation to remain non-zero and to approach a well defined limit. It is generally believed that this limiting value of the dissipation is a property of the Euler equations themselves, independent of the particular dissipative mechanism involved, and that it can be normalized with the large scale properties of the turbulent flow (e.g. the kinetic energy per unit volume u'(exp 2)/2, and the integral scale L) without reference to the Reynolds number or to other dissipative quantities. This is usually taken to imply that the low wave number end of the energy spectrum, far from the dissipative range, is also independent of the particular mechanism chosen to dispose of the energy transfer. In the following sections, we present some numerical experiments on the effect of substituting different dissipation models into the truncated Euler equations. We will see that the effect is mainly felt in the 'near dissipation' range of the energy spectrum, but that this range can be quite wide in some cases, contaminating a substantial range of wave numbers. In the process, we will develop a 'practical' approximation to the subgrid energy transfer in isotropic turbulence, and we will gain insight into the structure of the nonlinear interactions among turbulent scales of comparable size, and into the nature of energy backscatter. Some considerations on future research directions are offered at the end.

Jimenez, Javier↗
NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗
NASA NTRS1994

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗
DOE OSTI2024

An extended variational method for the resistive wall mode in toroidal plasma confinement devices

The external-kink stability of a toroidal plasma surrounded by a rigid resistive wall is investigated. The well-known analysis of Haney and Freidberg is rigorously extended to allow for a wall that is sufficiently thick that the thin-shell approximation does not necessarily hold. A generalized Haney–Freidberg formula for the growth-rate of the resistive wall mode is obtained. Thick-wall effects do not change the marginal stability point of the mode but introduce an interesting asymmetry between growing and decaying modes. Growing modes have growth-rates that exceed those predicted by the original Haney–Freidberg formula. On the other hand, decaying modes have decay-rates that are less than those predicted by the original formula. The well-known Hu–Betti formula for the rotational stabilization of the resistive wall mode is also generalized to take thick-wall effects into account. Increasing wall thickness facilitates the rotational stabilization of the mode, because it decreases the critical toroidal electromagnetic torque that the wall must exert on the plasma. On the other hand, the real frequency of the mode at the marginal stability point increases with increasing wall thickness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

'Downward control' of the mean meridional circulation and temperature distribution of the polar winter stratosphere

According to the 'downward control' principle, the extratropical mean vertical velocity on a given pressure level is approximately proportional to the meridional gradient of the vertically integrated zonal force per unit mass exerted by waves above that level. In this paper, a simple numerical model that includes parameterizations of both planetary and gravity wave breaking is used to explore the influence of gravity wave breaking in the mesosphere on the mean meridional circulation and temperature distribution at lower levels in the polar winter stratosphere. The results of these calculations suggest that gravity wave drag in the mesosphere can affect the state of the polar winter stratosphere down to altitudes below 30 km. The effect is most important when planetary wave driving is relatively weak: that is, during southern winter and in early northern winter. In southern winter, downwelling weakens by a factor of 2 near the stratospause and by 20% at 30 km when gravity wave drag is not included in the calculations. As a consequence, temperatures decrease considerably throughout the polar winter stratosphere (over 20 K above 40 km and as much as 8 K at 30 km, where the effect is enhanced by the long radiative relaxation timescale). The polar winter states obtained when gravity wave drag is omitted in this simple model resemble the results of simulations with some general circulation models and suggest that some of the shortcomings of the latter may be due to a deficit in mesospheric momentum deposition by small-scale gravity waves.

Garcia, Rolando R.↗
NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗
NASA NTRS2009

High Temperature VARTM of Phenylethynyl Terminated Imides

Fabrication of composite structures using vacuum assisted resin transfer molding (VARTM) is generally more affordable than conventional autoclave techniques. Recent efforts have focused on adapting VARTM for the fabrication of high temperature composites. Due to their low melt viscosity and long melt stability, certain phenylethynyl terminated imides (PETI) can be processed into composites using high temperature VARTM (HT-VARTM). However, one of the disadvantages of the current HT-VARTM resin systems has been the high porosity of the resultant composites. For aerospace applications, the desired void fraction of less than 2% has not yet been achieved. In the current study, two PETI resins, LaRC PETI-330 and LaRC PETI-8 have been used to make test specimens using HT-VARTM. The resins were infused into ten layers of IM7-6K carbon fiber 5-harness satin fabric at 260 C or 280 C and cured at 371 C. Initial runs yielded composites with high void content, typically greater than 7% by weight. A thermogravimetric-mass spectroscopic study was conducted to determine the source of volatiles leading to high porosity. It was determined that under the thermal cycle used for laminate fabrication, the phenylethynyl endcap was undergoing degradation leading to volatile evolution. By modifying the thermal cycle used in laminate fabrication, the void content was reduced significantly (typically approximately 3%). Densities of the composites were determined using a density gradient column and the glass transition temperatures of the cured composites were measured by dynamic mechanical analysis. Photomicrographs of the panels were taken and void contents were determined by acid digestion. The results of this work are presented herein.

Cano, Roberto J.↗
NASA NTRS1994

On the temperature and velocity through the photosphere of a sunspot penumbra

We investigate the structure in depth of a sunspot penumbra by means of the inversion code of the radiative transfer equation proposed by Ruiz Cobo & del Toro Iniesta (1992), applied to a set of filtergrams of a sunspot, scanning the Fe I line at 5576.1 A, with a sampling interval of 30 mA, from -120 to 120 mA from line center (data previously analyzed by Title et al. 1993). The temperature structure of this penumbra is obtained for each of the 801 pixels selected (0.32 sec x 0.32 sec). On the average, the temperatures seem to decrease as we move inward, but the differences are of the order of the rms values (approximately equal 100-200 K) at a given distance to sunspot center. The outer parts of the penumbra have also a bigger curvature in the T versus log tau(sub 5) relation than the inner parts. We realize, however, that these differences might be influenced by possible stray light effects. Compared to the quiet Sun, penumbral temperatures are cooler at deep layers and hotter at high layers. A mean penumbral model atmosphere is presented. The asymmetries observed in the intensity profile (the line is magnetically insensitive) are deduced to be produced by strong gradients of the line-of-sight velocity that sharply vary spatially along slices of almost constant distance to sunspot center. These variations suggest that such gradients are not only needed to explain the broadband circular polarization observed in sunspots (see Sanchez Almeida & Lites 1992) but are a main characteristic of the fine-scale penumbra. The results are compatible with an Evershed flow present everywhere, but its gradient with depth turns out to vary so that the flow seems to be mainly concentrated in some penumbral fibrils when studied through Dopplergrams. Finally, as by-products of this study, we put constraints to the practical usefulness of the Eddington-Barbier relation, and we explain the values of the Fourier Dopplergrams to be carrying information of layers around the centroid of the generalized response function of Dopplergrams to velocity fluctuations.

Del Toro Iniesta, J. C.↗
NASA NTRS1996

Regulation of Methane Oxidation in a Freshwater Wetland by Water Table Changes and Anoxia

The effects of water table fluctuations and anoxia on methane emission and methane oxidation were studied in a freshwater marsh. Seasonal aerobic methane oxidation rates varied between 15% and 76% of the potential diffusive methane flux (diffusive flux in the absence of aerobic oxidation). On an annual basis, approximately 43% of the methane diffusing into the oxic zone was oxidized before reaching the atmosphere. The highest methane oxidation was observed when the water table was below the peat surface. This was confirmed in laboratory experiments where short-term decreases in water table levels increased methane oxidation but also net methane emission. Although methane emission was generally not observed during the winter, stems of soft rush (Juncus effusus) emitted methane when the marsh was ice covered. Indigenous methanotrophic bacteria from the wetiand studied were relatively anoxia tolerant. Surface peat incubated under anoxic conditions maintained 30% of the initial methane oxidation capacity after 32 days of anoxia. Methanotrophs from anoxic peat initiated aerobic methane oxidation relatively quickly after oxygen addition (1-7 hours). These results were supported by culture experiments with the methanotroph Methylosinus trichosporium OB3b. This organism maintained a greater capacity for aerobic methane oxidation when starved under anoxic compared to oxic conditions. Anoxic incubation of M. trichosporium OB3b in the presence of sulfide (2 mM) and a low redox potential (-110 mV) did not decrease the capacity for methane oxidation relative to anoxic cultures incubated without sulfide. The results suggest that aerobic methane oxidation was a major regulator of seasonal methane emission front the investigated wetland. The observed water table fluctuations affected net methane oxidation presumably due to associated changes in oxygen gradients. However, changes from oxic to anoxic conditions in situ had relatively little effect on survival of the methanotrophic bacteria and thus on methane oxidation potential per se.

Roslev, Peter↗