Search NASA⌕ Search

SEARCH · Search NASA

Results for “QUADRATURE APPROXIMATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗
NASA NTRS2000

Flight Dynamics Analysis for Leonardo-BRDF

Leonardo-BRDF (Bidirectional Reflectance Distribution Function) is a new NASA mission concept proposed to allow the investigation of radiative transfer and its effect on the Earth's climate and atmospheric phenomenon. Enabled by the recent developments in small-satellite and formation flying technology, the mission is envisioned to be composed of an array of spacecraft in carefully designed orbits. The different perspectives provided by a distributed array of spacecraft offer a unique advantage to study the Earth's albedo. This paper presents the flight dynamics analysis performed in the context of the Leonardo-BRDF science requirements. First, the albedo integral is investigated and the effect of viewing geometry on science return is studied. The method used in this paper, based on Gauss quadrature, provides the optimal formation geometry to ensure that the value of the integral is accurately approximated. An orbit design approach is presented to achieve specific relative orbit geometries while simultaneously satisfying orbit dynamics constraints to reduce formation-keeping fuel expenditure. The relative geometry afforded by the design is discussed in terms of mission requirements. An optimal Lambert initialization scheme is presented with the required Delta-V to distribute all spacecraft from a common parking orbit into their appropriate orbits in the formation. Finally, formation-keeping strategies are developed and the associated Delta-V's are calculated to maintain the formation in the presence of perturbations.

Hughes, Steven P.↗
NASA NTRS2001

A Preliminary Formation Flying Orbit Dynamics Analysis for Leonardo-BRDF

Leonardo-BRDF is a NASA mission concept proposed to allow the investigation of radiative transfer and its effect on the Earth's climate and atmospheric phenomenon. Enabled by the recent developments in small-satellite and formation flying technology, the mission is envisioned to be composed of an array of spacecraft in carefully designed orbits. The different perspectives provided by a distributed array of spacecraft offer a unique advantage to study the Earth's albedo. This paper presents the orbit dynamics analysis performed in the context of the Leonardo-BRDF science requirements. First, the albedo integral is investigated and the effect of viewing geometry on science return is studied. The method used in this paper, based on Gauss quadrature, provides the optimal formation geometry to ensure that the value of the integral is accurately approximated. An orbit design approach is presented to achieve specific relative orbit geometries while simultaneously satisfying orbit dynamics constraints to reduce formation-keeping fuel expenditure. The relative geometry afforded by the design is discussed in terms of mission requirements. An optimal two-burn initialization scheme is presented with the required delta-V to distribute all spacecraft from a common parking orbit into their appropriate orbits in the formation. Finally, formation-keeping strategies are developed and the associated delta-V's are calculated to maintain the formation in the presence of perturbations.

Hughes, Steven P.↗
NASA NTRS2001

A Preliminary Formation Flying Orbit Dynamics Analysis for Leonardo-BRDF

Leonardo-BRDF is a new NASA mission concept proposed to allow the investigation of radiative transfer and its effect on the Earth's climate and atmospheric phenomenon. Enabled by the recent developments in small-satellite and formation flying technology, the mission is envisioned to be composed of an array of spacecraft in carefully designed orbits. The different perspectives provided by a distributed array of spacecraft offer a unique advantage to study the Earth's albedo. This paper presents the flight dynamics analysis performed in the context of the Leonardo-BRDF science requirements. First, the albedo integral is investigated and the effect of viewing geometry on science return is studied. The method used in this paper, based on Gauss quadrature, provides the optimal formation geometry to ensure that the value of the integral is accurately approximated. An orbit design approach is presented to achieve specific relative orbit geometries while simultaneously satisfying orbit dynamics constraints to reduce formation-keeping fuel expenditure. The relative geometry afforded by the design is discussed in terms of mission requirements. An optimal Lambert initialization scheme is presented with the required DeltaV to distribute all spacecraft from a common parking orbit into their appropriate orbits in the formation. Finally, formation-keeping strategies are developed and the associated DeltaV's are calculated to maintain the formation in the presence of perturbations.

Hughes, Steven P.↗
NASA NTRS2020

Examining Biases in Diurnally-Integrated Shortwave Irradiances due to Two- and Four-Stream Approximations in Cloudy Atmosphere

Shortwave irradiance biases due to two- and four-stream approximations have been studied for the last couple of decades, but biases in estimating Earth’s radiation budget have not been examined in earlier studies. In order to quantify biases in diurnally-averaged irradiances, we integrate the two- and four-stream biases using realistic diurnal variations of cloud properties from Clouds and the Earth’s Radiant Energy System (CERES) synoptic (SYN) hourly product. Three approximations are examined in this study, delta-two-stream-Eddington (D2strEdd), delta- two-stream-quadrature (D2strQuad), and delta-four-stream-quadrature (D4strQuad). Irradiances computed by the Discrete Ordinates Radiative Transfer (DISORT) and Monte Carlo (MC) methods are used as references. The MC noises are further examined by comparing with DISORT results. When the biases are integrated with a one-day of solar zenith angle variation, regional biases of D2strEdd and D2strQuad reach up to 8 W/sq.m, while biases of D4strQuad reach up to 2 W/sq.m. When the biases are further averaged monthly or annually, regional biases of D2strEdd and D2strQuad can reach –1.5 W/sq.m in SW top-of-atmosphere (TOA) upward irradiances and +3 W/sq.m in surface downward irradiances. In contrast, regional biases of D4strQuad are within +0.9 for TOA irradiances and –1.2 W/sq.m for surface irradiances. Except for polar regions, monthly and annual global mean biases are similar, suggesting that the biases are nearly independent to season. Biases in SW heating rate profiles are up to –0.008 K/d for D2strEdd and –0.016 K/d for D2strQuad, while the biases of the D4strQuad method are negligible.

Ham, Seung-Hee↗
NASA NTRS1971

Effect of angular quadrature on results of two-dimensional space power reactor shield calculations.

Study of the effect of angular quadrature on the results of two-dimensional transport calculations made to evaluate radiation environment in the vicinity of a shielded space power reactor. The application considered is for a manned space-station; crew shielding must be provided within a 25 deg cone angle. The asymmetric shield layout is evaluated using the two-dimensional transport code DOT. Within the shield, an S sub 6 angular quadrature yields satisfactory results. For regions where the ray effect may occur it is not apparent that an S sub 10 quadrature is adequate. Anomalies in results due to transport of particles in a void and necessary geometric approximations were observed.

Connolley, D. J.↗
NASA NTRS1990

Two-dimensional radiation and scattering at short wave length

In the present investigation of radiation and scattering by objects when the wavelengths are much smaller than the characteristic dimensions of the radiator or scatterer, the boundary-element method is used to obtain computational accuracy and efficiency. The approach employed for wavelengths less than 5 percent of object characteristic dimensions involves cubic elements, approximate polynomial and asymptotic evaluations of the fundamental solution, and a tailoring of the order of the Gaussian quadrature according to the local demands dictated by the distance between sending and receiving points. The method addresses the propagation of low-frequency sound over large terrain features.

Yoon, W. S.↗
NASA NTRS1988

Single-grid spectral collocation for the Navier-Stokes equations

The aim of the paper is to study a collocation spectral method to approximate the Navier-Stokes equations: only one grid is used, which is built from the nodes of a Gauss-Lobatto quadrature formula, either of Legendre or of Chebyshev type. The convergence is proven for the Stokes problem provided with inhomogeneous Dirichlet conditions, then thoroughly analyzed for the Navier-Stokes equations. The practical implementation algorithm is presented, together with numerical results.

Bernardi, Christine↗
NASA NTRS2004

A Slow Streamer Blowout at the Sun and Ulysses

On 10 June 2000 a streamer on the southeast limb slowly disappeared from LASCO/C2 over approximately 10 hours. A small CME was reported in C2. A substantial interplanetary CME (ICME) was later detected at Ulysses, which was at quadrature with the Sun and SOHO at the time. This detection illustrates the properties of an ICME for a known solar source and demonstrates that the identification can be done even beyond 3 AU. Slow streamer blowouts such as this have long been known but are little studied. We report on the SOHO observation of a coronal mass ejection (CME) on the solar limb and the subsequent in situ detection at Ulysses, which was near quadrature at the time, above the location of the CME. SOHO-Ulysses quadrature was 13 June, when Ulysses was 3.36 AU from the Sun and 58.2 degrees south of the equator off the east limb. The slow streamer blowout was on 10 June, when the SOHO-Sun-Ulysses angle was 87 degrees.

Seuss, S. T.↗
NASA NTRS2018

Combined Error and Uncertainty Estimates for CFD Problems

Given input sources of uncertainty, non-intrusive uncertainty propagation methods quantify the uncertainty in output quantities of interest (QoI) by performing a nite number of CFD (Computational Fluid Dynamics) instance realizations needed in the calculation of output statistics. It is well known that this introduces multiple sources of error. CFD codes often utilize finite-dimensional approximation (grids, basis functions, etc.) thus incurring CFD numerical errors often approximately reinterpreted as a statistical bias. Uncertainty propagation methods calculate uncertainty statistics for output quantities of interest using a numerical method (e.g. deterministic quadrature, sampling, etc.) thus incurring UQ (Uncertainty Quantification) numerical errors. Importance of quantifying these errors in large scale scientific computing: How accurate is an output statistic?; How should additional computational resources be invested to further reduce the error in a statistic?

Posteriori↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
NASA NTRS2002

Digital Distortion Caused by Traveling- Wave-Tube Amplifiers Simulated

Future NASA missions demand increased data rates in satellite communications for near real-time transmission of large volumes of remote data. Increased data rates necessitate higher order digital modulation schemes and larger system bandwidth, which place stricter requirements on the allowable distortion caused by the high-power amplifier, or the traveling-wave-tube amplifier (TWTA). In particular, intersymbol interference caused by the TWTA becomes a major consideration for accurate data detection at the receiver. Experimentally investigating the effects of the physical TWTA on intersymbol interference would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to acquiring the required digital hardware. Thus, an accurate computational model is essential to predict the effects of the TWTA on system-level performance when a communication system is being designed with adequate digital integrity for high data rates. A fully three-dimensional, time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell code MAFIA (Solution of Maxwell's equations by the Finite-Integration-Algorithm). It comprehensively takes into account the effects of frequency-dependent AM (amplitude modulation)/AM and AM/PM (phase modulation) conversion, gain and phase ripple due to reflections, drive-induced oscillations, harmonic generation, intermodulation products, and backward waves. This physics-based TWT model can be used to give a direct description of the effects of the nonlinear TWT on the operational signal as a function of the physical device. Users can define arbitrary excitation functions so that higher order modulated digital signals can be used as input and that computations can directly correlate intersymbol interference with TWT parameters. Standard practice involves using communication-system-level software packages, such as SPW, to predict if adequate signal detection will be achieved. These models use a nonlinear, black-box model to represent the TWTA. The models vary in complexity, but most make several assumptions regarding the operation of the high-power amplifier. When the MAFIA TWT interaction model was used, these assumptions were found to be in significant error. In addition, digital signal performance, including intersymbol interference, was compared using direct data input into the MAFIA model and using the system-level analysis tool SPW for several higher order modulation schemes. Results show significant differences in predicted degradation between SPW and MAFIA simulations, demonstrating the significance of the TWTA approximations made in the SPW model on digital signal performance. For example, a comparison of the SPW and MAFIA output constellation diagrams for a 16-ary quadrature amplitude modulation (16-QAM) signal (data shown only for second and fourth quadrants) is shown. The upper-bound degradation was calculated from the corresponding eye diagrams. In comparison to SPW simulations, the MAFIA data resulted in a 3.6-dB larger degradation.

Kory, Carol L.↗
NASA NTRS1976

Some new results on the statistics of radio wave scintillation. I - Empirical evidence for Gaussian statistics

This paper presents an analysis of ionospheric scintillation data which shows that the underlying statistical structure of the signal can be accurately modeled by the additive complex Gaussian perturbation predicted by the Born approximation in conjunction with an application of the central limit theorem. By making use of this fact, it is possible to estimate the in-phase, phase quadrature, and cophased scattered power by curve fitting to measured intensity histograms. By using this procedure, it is found that typically more than 80% of the scattered power is in phase quadrature with the undeviated signal component. Thus, the signal is modeled by a Gaussian, but highly non-Rician process. From simultaneous UHF and VHF data, only a weak dependence of this statistical structure on changes in the Fresnel radius is deduced. The signal variance is found to have a nonquadratic wavelength dependence. It is hypothesized that this latter effect is a subtle manifestation of locally homogeneous irregularity structures, a mathematical model proposed by Kolmogorov (1941) in his early studies of incompressible fluid turbulence.

Rino, C. L.↗
NASA NTRS2016

On the Calculation of Uncertainty Statistics with Error Bounds for CFD Calculations Containing Random Parameters and Fields

This chapter discusses the ongoing development of combined uncertainty and error bound estimates for computational fluid dynamics (CFD) calculations subject to imposed random parameters and random fields. An objective of this work is the construction of computable error bound formulas for output uncertainty statistics that guide CFD practitioners in systematically determining how accurately CFD realizations should be approximated and how accurately uncertainty statistics should be approximated for output quantities of interest. Formal error bounds formulas for moment statistics that properly account for the presence of numerical errors in CFD calculations and numerical quadrature errors in the calculation of moment statistics have been previously presented in [8]. In this past work, hierarchical node-nested dense and sparse tensor product quadratures are used to calculate moment statistics integrals. In the present work, a framework has been developed that exploits the hierarchical structure of these quadratures in order to simplify the calculation of an estimate of the quadrature error needed in error bound formulas. When signed estimates of realization error are available, this signed error may also be used to estimate output quantity of interest probability densities as a means to assess the impact of realization error on these density estimates. Numerical results are presented for CFD problems with uncertainty to demonstrate the capabilities of this framework.

Uncertainty↗
NASA NTRS1986

A method to dramatically improve subcarrier tracking

A method is presented for achieving a dramatic improvement in phase tracking of square wave subcarriers or other square waves. The method is to set the amplitude of the phase quadrature reference signal to zero except near the zero crossings of the input signal. Without changing the loop bandwidth, the variance of the phase error can be reduced to approximately W sigma(sub 0)(2), were sigma (sub 0)(2) is the phase error variance without windowing, and W is the fraction of cycle in which the reference signal has a nonzero value. Simulation results confirm the analysis and establish minimum W versus signal-to-noise ratio. Typically, the window can be made so narrow as to achieve a phase error variance of 1.5 sigma(sub 0)(4).

Hurd, W. J.↗
NASA NTRS1988

A method to dramatically improve subcarrier tracking

A method is presented for achieving a dramatic improvement in phase tracking of square wave subcarriers or other square waves. The method is to set the amplitude of the phase quadrature reference signal to zero except near the zero crossings of the input signal. Without changing the loop bandwidth, the variance of the phase error can be reduced to approximately W sigma(sub 0)(2), where sigma (sub 0)(2) is the phase error variance without windowing, and W is the fraction of cycle in which the reference signal has a nonzero value. Simulation results confirm the analysis and establish minimum W versus signal-to-noise ratio. Typically, the window can be made so narrow as to achieve a phase error variance of 1.5 sigma(sub 0)(4).

Hurd, William J.↗
NASA NTRS1990

An approach to optimal guidance of an advanced launch vehicle concept

An approximate solution for the maximum payload trajectory of a two-stage launch vehicle using a regular perturbation technique is presented. A zero-order solution for a two-stage vehicle based on a flat-earth approximation and negligible atmospheric effects is obtained in closed form. High-order correction terms are obtained from the solution of nonhomogeneous, first-order linear differential equations by quadrature. This promises the capability for an onboard optimal guidance law implementation.

Leung, Martin S.↗
NASA NTRS1994

Space proton transport in one dimension

An approximate evaluation procedure is derived for a second-order theory of coupled nucleon transport in one dimension. An analytical solution with a simplified interaction model is used to determine quadrature parameters to minimize truncation error. Effects of the improved method on transport solutions with the BRYNTRN data base are evaluated. Comparisons with Monte Carlo benchmarks are given. Using different shield materials, the computational procedure is used to study the physics of space protons. A transition effect occurs in tissue near the shield interface and is most important in shields of high atomic number.

NASA Discipline Number 45-10↗