Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 901 records · Page 50

The Benefits of Global High Resolution for Climate Simulation: Process Understanding and the Enabling of Stakeholder Decisions at the Regional Scale

Many think that the time scales of the Paris Climate Agreement indicate urgent action is required on climate policies over the next few decades, in order to avoid the worst risks posed by climate change. On these relatively short time scales the combined effect of climate variability and change are both key drivers of extreme events, with decadal time scales also important for infrastructure planning. Hence, in order to assess climate risk on such time scales, we require climate models to be able to represent key aspects of both internally driven climate variability and the response to changing forcings. In this paper we argue that we now have the modeling capability to address these requirements—specifically with global models having horizontal resolutions considerably enhanced from those typically used in previous Intergovernmental Panel on Climate Change (IPCC) and Coupled Model Intercomparison Project (CMIP) exercises. The improved representation of weather and climate processes in such models underpins our enhanced confidence in predictions and projections, as well as providing improved forcing to regional models, which are better able to represent local-scale extremes (such as convective precipitation). We choose the global water cycle as an illustrative example because it is governed by a chain of processes for which there is growing evidence of the benefits of higher resolution. At the same time it comprises key processes involved in many of the expected future climate extremes (e.g., flooding, drought, tropical and midlatitude storms).

54 ENVIRONMENTAL SCIENCES↗

The earth resources technology satellite - Photographic facility quality control program

A detailed description is given of the quality control program used in the photographic laboratory of the NASA-ERTS Ground Data Handling System. The product response variables measured include tone reproduction, resolution, and low spatial frequency noise. In addition to product response variables, certain performance parameters of the laboratory printers and processors are frequently measured in order to produce consistent duplications of archival photography. A description is given of the operation and use of a densitometer/computer interface which is used to calculate three tone reproduction response variables - film speed, average gradient, and base plus fog density. This procedure eliminates the need for any hand plotting of D log E curves to manually determine response variables.

Shaffer, R. M.↗

Quantifying Contributions of Land Carbon Fluxes Variability and Atmospheric Transport Variability to Atmospheric CO2 Variability

Much of the uncertainty in predicting variations in greenhouse gases originates in the complex dynamics of the land ecosystem and the atmosphere. To reduce the uncertainties, it is useful to decouple the contributions of land carbon fluxes and atmospheric transport to atmospheric carbon variability. Here we isolate these contributions using a version of the NASA GEOS model that couples carbon, energy and water cycles between the land and the atmosphere. Our current study is a follow-on to a preliminary analysis that suggested that an extreme event (e.g., imposed drought) in a free running AGCM simulation affects atmospheric CO2 more through its impact on atmospheric transport than through a modification of land carbon fluxes. In order to more carefully isolate the effects of the land carbon variability and atmospheric transport variability on atmospheric CO2 variability, we conducted two coupled-AGCM simulations in replay mode, a mode that forces the model's evolution of weather to match that of the MERRA-2 reanalysis. In the control simulation, the land carbon fluxes and the atmospheric CO2 concentrations, as well as the meteorology, are simulated over 2001-2015. We compute the climatological seasonal cycles of net land carbon production from this control simulation and then prescribe, in a second simulation, these climatological carbon fluxes to the atmosphere in the same replay mode. By comparing the atmospheric CO2 variability produced in the two simulations, we fully isolate the part of this variability associated with land surface fluxes. Relative contributions of land flux variability and atmospheric transport variability to CO2 variability are quantified both on a regional basis and as a function of height into the atmosphere to support interpretation of both ground-based and satellite observations.

Lee, Eunjee↗

Effects of spatial order of accuracy on the computation of vortical flowfields

The effect of the order-of-accuracy, used for the spatial discretization, on the resolution of the leading edge vortices over sharp-edged delta wings is investigated. The flowfield is computed using a viscous/inviscid zonal approach. The viscous flow in the vicinity of the wing is computed using the conservative formulation of the compressible, thin-layer Navier-Stokes equations. The leeward-side vortical flowfield and the other flow regions away from the surface are computed as inviscid. The time integration is performed with both an explicit fourth-order Runge-Kutta scheme and an implicit, factorized, iterative scheme. High-order-accurate inviscid fluxes are computed using both a conservative and a non-conservative (primitive variable) formulation. The nonlinear, inviscid terms of the primitive variable form of the governing equations are evaluated with a finite-difference numerical scheme based on the sign of the eigenvalues. High-order, upwind-biased, finite difference formulas are used to evaluate the derivatives of the nonlinear convective terms. Computed results are compared with available experimental data, and comparisons of the flowfield in the vicinity of the vortex cores are presented.

Ekaterinaris, J. A.↗

True Airspeed Measurement by Ionization-Tracer Technique

Ion bundles produced in a pulse-excited corona discharge are used as tracers with a radar-like pulse transit-time measuring instrument in order to provide a measurement of airspeed that is independent of all variables except time and distance. The resulting instrumentation need not project into the air stream and, therefore, will not cause any interference in supersonic flow. The instrument was tested at Mach numbers ranging from 0.3 to 3.8. Use of the proper instrumentation and technique results in accuracy of the order of 1 percent.

Boyd, B.↗

Algorithms For Integrating Nonlinear Differential Equations

Improved algorithms developed for use in numerical integration of systems of nonhomogenous, nonlinear, first-order, ordinary differential equations. In comparison with integration algorithms, these algorithms offer greater stability and accuracy. Several asymptotically correct, thereby enabling retention of stability and accuracy when large increments of independent variable used. Accuracies attainable demonstrated by applying them to systems of nonlinear, first-order, differential equations that arise in study of viscoplastic behavior, spread of acquired immune-deficiency syndrome (AIDS) virus and predator/prey populations.

Freed, A. D.↗

Calibration of Atmospherically Induced Delay Fluctuations Due to Water Vapor

We have completed a new generation of water vapor radiometers (WVR), the A- series, in order to support radio science experiments with the Cassini spacecraft. These new instruments sense three frequencies in the vicinity of the 22 GHz emission line of atmospheric water vapor within a 1 degree beamwidth from a clear aperture antenna that is co-pointed with the radio telescope down to 10 degree elevation. The radiometer electronics features almost an order of magnitude improvement in temperature stability compared with earlier WVR designs. For many radio science experiments, the error budget is likely to be dominated by path delay fluctuations due to variable atmospheric water vapor along the line-of-sight to the spacecraft. In order to demonstrate the performance of these new WVRs we are attempting to calibrate the delay fluctuations as seen by a radio interferometer operating over a 21 km baseline with a WVR near each antenna. The characteristics of these new WVRs will be described and the results of our preliminary analysis will be presented indicating an accuracy of 0.2 to 0.5 mm in tracking path delay fluctuations over time scales of 10 to 10,000 seconds.

Resch, George↗

Hermite analogs of the lowest order Raviart–Thomas mixed method for convection–diffusion equations

The Raviart–Thomas mixed finite-element method of the lowest order (Raviart and Thomas in mixed finite-element methods for second-order elliptic problems, Lecture Notes in mathematics, Springer, New York, 1977) commonly known as the RT{sub 0} method, is a well-established and popular numerical tool to solve diffusion-like problems providing flux continuity across inter-element boundaries. Douglas and Roberts extended the method to the case of more general second-order boundary-value problems including the convection–diffusion equations (cf. this journal Douglas in Comput Appl Math 1:91–103; 1982). The main drawback of these methods, however, is the poor representation of the primal variable by piecewise constant functions. The Hermite analog of the RT{sub 0} method for treating pure diffusion phenomena proposed in Ruas (J Comput Appl Math, 246:234–242; 2013) proved to be a valid alternative to attain higher order approximation of the primal variable while keeping intact the matrix structure and the quality of the discrete flux variable of the original RT{sub 0} method. Non-trivial extensions of this method are studied here that can be viewed as Hermite analogs of the two Douglas and Roberts’ versions of the RT{sub 0} method, to solve convection–diffusion equations. A detailed convergence study is carried out for one of the Hermite methods, and numerical results illustrate the performance of both of them, as compared to each other and to the corresponding mixed methods.

97 MATHEMATICS AND COMPUTING↗

Taylor Series Trajectory Calculations Including Oblateness Effects and Variable Atmospheric Density

Taylor series integration is implemented in NASA Glenn's Spacecraft N-body Analysis Program, and compared head-to-head with the code's existing 8th- order Runge-Kutta Fehlberg time integration scheme. This paper focuses on trajectory problems that include oblateness and/or variable atmospheric density. Taylor series is shown to be significantly faster and more accurate for oblateness problems up through a 4x4 field, with speedups ranging from a factor of 2 to 13. For problems with variable atmospheric density, speedups average 24 for atmospheric density alone, and average 1.6 to 8.2 when density and oblateness are combined.

Scott, James R.↗

Pole placement and order reduction in two-time-scale control systems through Riccati iteration

A transformation of variables taken from singular perturbations may be applied to two-time-scale linear systems in state space form to reduce the system to block-diagonal form with slow and fast modes decoupled. The transformation is easily computed by applying the new Riccati iteration. The iteration yields a solution to the nonsymmetric algebraic Riccati equation obtained by partitioning the original system matrix A. The numerical procedure is initiated with the trivial iterate L(0) = 0, and is globally convergent to the desired unique time scale decoupling solution. After transformation, the decoupled system may be used in controller design to achieve exact closed loop pole placement in the slow subsystem without altering the poles of the fast subsystem. The decoupled form may also be used to reduce system order by wetting a small parameter to zero. Provided the fast subsystem is stable, the order reduction can be expected to yield a good approximation to the original system. These methods are demonstrated using the 16th order linear model of a turbofan engine.

Anderson, L. R.↗

Impact of FERC Order 2222 on DER Participation Rules in US Electricity Markets

Electricity markets in the bulk grid are beginning to implement market mechanisms that support the procurement of flexible capabilities from wide range of technologies, including distributed energy resources (DERs). The flexibility of these resources will help counterbalance supply uncertainties from large-scale integration of variable renewable generation. To encourage development of distributed and aggregated market participants, FERC Order 2222 was issued in September 2020 to require each Independent System Operator (ISO) in the US to implement rules that enable broader participation from aggregations of DERs in the bulk market. The following paper first describes the generic design of ISO markets before introducing the new market participation rules that ISOs have proposed for compliance with Order 2222. The paper then describes how software performance issues may continue to affect the eligibility requirements and offer structures for DER aggregations participating in ISOs, noting that continued research on computational methods may help reduce burdens for DER integration. The prospects for transmission and distribution system coordination is second major issue discussed, which will require minor changes to existing processes in the short term. In the longer term, there is more opportunity for more wide-ranging reforms, such as the development of a Distribution System Operator (DSO) framework. Newly proposed market rules may affect how Transactive Energy Systems (TES) will help facilitate efficient formation of DER aggregations and operation of the individual DERs within an aggregation. Within the TES context, the challenge is to fully understand how resource eligibility and operational and planning coordination methods will affect the design and implementation of TES.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The relationship between eddy-transport and second-order closure models for stratified media and for vortices

The question is considered of how complex a model should be used for the calculation of turbulent shear flows. At the present time there are models varying in complexity from very simple eddy-transport models to models in which all the equations for the nonzero second-order correlations are solved simultaneously with the equations for the mean variables. A discussion is presented of the relationship between these two models of turbulent shear flow. Two types of motion are discussed: first, turbulent shear flow in a stratified medium and, second, the motion in a turbulent line vortex. These two cases are instructive because in the first example eddy-transport methods have proven reasonably effective, whereas in the second, they have led to erroneous conclusions. It is not generally appreciated that the simplest form of eddy-transport theory can be derived from second-order closure models of turbulent flow by a suitably limiting process. This limiting process and the suitability of eddy-transport modeling for stratified media and line vortices are discussed.

Donaldson, C. D.↗

Characterizing the Influence of the General Circulation on Marine Boundary Layer Clouds

The seasonal and intraseasonal variability of boundary layer cloud in the subtropical eastern oceans are studied using combined data from the International Satellite Cloud Climatology Project (ISCCP) and the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis. Spectral analysis reveals that most of the time variability of cloud properties occurs on seasonal to annual time scales. The variance decreases one to two orders of magnitude for each decade of time scale decrease, indicating that daily to monthly time scales have smaller, but non-negligible variability. The length of these dominant time scales suggests that the majority of the variability is influenced by the general circulation and its interaction with boundary layer turbulence, rather than a product of boundary layer turbulence alone. Previous datasets have lacked the necessary resolution in either time or in space to properly characterize variability on synoptic scales; this is remedied by using global satellite-retrieved cloud properties. We characterize the intraseasonal subtropical cloud variability in both hemispheres and in different seasons. In addition to cloud fraction, we examine variability of cloud optical thickness - cloud top pressure frequency distributions. Despite the large concentration of research on the variability of Northern Hemisphere (NH) regions during summer, it is noted that the largest amplitude intraseasonal variability in the NH regions occurs during local winter. The effect of intraseasonal variability on the calculation and interpretation of seasonal results is investigated. Decreases in seasonally averaged cloud cover, optical thickness and cloud top pressure from the May-through-September season to the November-through-March season are most apparent in the NH regions. Further analysis indicates that these changes are due to an increase in frequency, but a decrease in the persistence of synoptic events. In addition, changes in cloud top pressure and optical thickness characteristics from the summer to winter seasons indicate that the NH subtropics undergo a change in dynamic regime with season. This change appears in the cloud fields as a shift from the more commonly seen lower-altitude, thicker optical thickness clouds to higher-altitude, thinner clouds. The latter cloud-type is associated with the lower sea level pressure, upward vertical velocity phase of the synoptic wave. Intraseasonal changes in cloud properties in the Southern Hemisphere and NH summer are much smaller in amplitude. Although they also appear to be linked to changes in the large-scale dynamics, similarly to NH winter variations, the relationships are more ambiguous due to the small amplitudes and longer time scales. We attempt to interpret some of these relationships using the results of the Betts and Ridgway (1989) box model. However, these results cannot consistently explain the patterns when results from all regions are considered, implying that this model may not adequately explain all the processes involved in the variability.

Rozendaal, Margaret A.↗

Keplerian Motion of Broad-Line Region Gas as Evidence for Supermassive Black Holes in Active Galactic Nuclei

Emission-line variability data on NGC 5548 argue strongly for the existence of a mass of order 7 x l0(exp 7) solar mass, within the inner few light-days of the nucleus in the Seyfert 1 galaxy NGC 5548. The time-delayed response of the emission lines to continuum variations is used to infer the size of the line-emitting region, and these determinations are combined with measurements of the Doppler widths of the variable line components to estimate a virial mass. The data for several different emission lines spanning an order of magnitude in distance from the central source show the expected V proportional to r(sup -1/2) correlation and are consistent with a single value for the mass.

Peterson, Bradley M.↗

Competition between magnetic order and charge localization in Na 2 IrO 3 thin crystal devices thin crystal devices

Spin orbit assisted Mott insulators such as sodium iridate (Na 2 IrO 3 ) have been an important subject of study in recent years. In these materials, the interplay of electronic correlations, spin-orbit coupling, crystal field effects, and a honeycomb arrangement of ions bring exciting ground states, predicted in the frame of the Kitaev model. The insulating character of Na 2 IrO 3 has hampered its integration to an electronic device, desirable for applications, such as the manipulation of quasiparticles interesting for topological quantum computing. Here we show through electronic transport measurements supported by angle-resolved photoemission spectroscopy (ARPES) experiments, that electronic transport in Na 2 IrO 3 is ruled by variable range hopping and it is strongly dependent on the magnetic ordering transition known for bulk Na 2 IrO 3 , as well as on external electric fields. Furthermore, electronic transport measurements allow us to deduce a value for the localization length and the density of states in our Na 2 IrO 3 thin crystal devices, and offer an alternative approach to study insulating 2D-materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Order of accuracy of QUICK and related convection-diffusion schemes

This report attempts to correct some misunderstandings that have appeared in the literature concerning the order of accuracy of the QUICK scheme for steady-state convective modeling. Other related convection-diffusion schemes are also considered. The original one-dimensional QUICK scheme written in terms of nodal-point values of the convected variable (with a 1/8-factor multiplying the 'curvature' term) is indeed a third-order representation of the finite volume formulation of the convection operator average across the control volume, written naturally in flux-difference form. An alternative single-point upwind difference scheme (SPUDS) using node values (with a 1/6-factor) is a third-order representation of the finite difference single-point formulation; this can be written in a pseudo-flux difference form. These are both third-order convection schemes; however, the QUICK finite volume convection operator is 33 percent more accurate than the single-point implementation of SPUDS. Another finite volume scheme, writing convective fluxes in terms of cell-average values, requires a 1/6-factor for third-order accuracy. For completeness, one can also write a single-point formulation of the convective derivative in terms of cell averages, and then express this in pseudo-flux difference form; for third-order accuracy, this requires a curvature factor of 5/24. Diffusion operators are also considered in both single-point and finite volume formulations. Finite volume formulations are found to be significantly more accurate. For example, classical second-order central differencing for the second derivative is exactly twice as accurate in a finite volume formulation as it is in single-point.

Leonard, B. P.↗

An investigation of a mathematical model of an optically pumped Ti(3+):Al2O3 laser system

During the last several years, solid state lasers were developed that have the potential for meeting rigorous performance requirements for space-based remote sensing of the atmosphere. In order to design a stable and efficient laser and to understand the effect on laser output of changes in the physical and design parameters, an understanding of the development of the dynamical processes of the laser is necessary. Typically, the dynamical processes in a laser system are investigated via rate equations describing the evolution of the occupancy in the electronic levels and of the photon density in the laser cavity. There are two approaches to this type of study. Most often, for the sake of simplicity, the spatial variations of the dynamic variables in the laser system are disregarded and the mathematical model consists of a system of first order nonlinear ordinary differential equations (ODE). The second approach is to take into account both spatial and temporal variations in the dynamic variables in the laser cavity. The resulting model consists of a first order semilinear system of partial differential equations (PDE). The model which was studied was studied was generic in the sense that it was a four-level laser system, but the parameters used in the numerical study were specific to Titanium-doped sapphire. For simplicity, a constant, spatially uniform pumping scheme was considered. In addition, a simplification of the model was made so that it treats a single lasing wavelength with a narrow bandwidth. The purpose was to investigate both versions of the mathematical model and to determine whether the numerical solutions are similar both qualitatively and quantitatively. The systems of ordinary differential equations were solved numerically using a Runge-Kutta-Fehlberg algorithm which was very efficient for typical values of the physical parameters. A numerical scheme, based on the Modified Euler method, for computing solutions to the system of partial differential equations was developed and implemented. The PDE model was solved numerically at the expense of greatly increased computer time.

Roberts, Lila F.↗

Uncertainty and Sensitivity Analysis of Afterbody Radiative Heating Predictions for Earth Entry

The objective of this work was to perform sensitivity analysis and uncertainty quantification for afterbody radiative heating predictions of Stardust capsule during Earth entry at peak afterbody radiation conditions. The radiation environment in the afterbody region poses significant challenges for accurate uncertainty quantification and sensitivity analysis due to the complexity of the flow physics, computational cost, and large number of un-certain variables. In this study, first a sparse collocation non-intrusive polynomial chaos approach along with global non-linear sensitivity analysis was used to identify the most significant uncertain variables and reduce the dimensions of the stochastic problem. Then, a total order stochastic expansion was constructed over only the important parameters for an efficient and accurate estimate of the uncertainty in radiation. Based on previous work, 388 uncertain parameters were considered in the radiation model, which came from the thermodynamics, flow field chemistry, and radiation modeling. The sensitivity analysis showed that only four of these variables contributed significantly to afterbody radiation uncertainty, accounting for almost 95% of the uncertainty. These included the electronic- impact excitation rate for N between level 2 and level 5 and rates of three chemical reactions in uencing N, N(+), O, and O(+) number densities in the flow field.

West, Thomas K., IV↗