Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Structure Corrections in Modeling VLBI Delays for RDV Data

Since 1997, bimonthly S- and X-band observing sessions have been carried out employing the VLBA (Very Long Baseline Array) and as many as ten additional antennas. Maps of the extended structures have been generated for the 160 sources observed in ten of these experiments (approximately 200,000 observations) taking place during 1997 and 1998. This paper reports the results of the first massive application of such structure maps to correct the modeled VLBI (Very Long Baseline Interferometry) delay in astrometric data analysis. For high-accuracy celestial reference frame work, proper choice of a reference point within each extended source is crucial. Here the reference point is taken at the point of maximum emitted flux. Overall, the weighted delay residuals (approximately equal to 30 ps) are reduced by 8 ps in quadrature upon introducing source maps to model the structure delays of the sources. Residuals of some sources with extended or fast-varying structures improve by as much as 40 ps. Scatter of 'arc positions' about a time-linear model decreases substantially for most sources. Based on our results, it is also concluded that source structure is presently not the dominant error source in astrometric/geodetic VLBI.

Sovers, Ojars J.↗
NASA NTRS1988

Focal region fields of distorted reflectors

The problem of the focal region fields scattered by an arbitrary surface reflector under uniform plane wave illumination is solved. The physical optics (PO) approximation is used to calculate the current induced on the reflector. The surface of the reflector is described by a number of triangular domain-wise 5th degree bivariate polynomials. A 2-dimensional Gaussian quadrature is employed to numerically evaluate the integral expressions of the scattered fields. No Freshnel or Fraunhofer zone approximations are made. The relation of the focal fields problem to surface compensation techniques and other applications are mentioned. Several examples of distorted parabolic reflectors are presented. The computer code developed is included, together with instructions on its usage.

Buris, N. E.↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

A defining feature of the discontinuous Galerkin (DG) method for ODE is that the piecewise polynomial solution can have a jump discontinuity at the beginning of each step. Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠 − 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK-DG, collocation, and CG schemes.

numerical methods↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠− 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK DG, collocation, and CG methods.

Numerical Methods for Ordinary Differential Equati↗
NASA NTRS1990

A real-time approximate optimal guidance law for flight in a plane

A real-time guidance scheme is presented for the problem of maximizing the payload into orbit subject to the equations of motion of a rocket over a nonrotating spherical earth. The flight is constrained to a path in the equatorial plane while reaching an orbital altitude at orbital injection speeds. The dynamics of the problem can be separated into primary and perturbation effects by a small parameter, epsilon, which is the ratio of the atmospheric scale height to the radius of the earth. The Hamilton-Jacobi-Bellman or dynamic programming equation is expanded in an asymptotic series where the zeroth-order term (epsilon = 0) can be obtained in closed form. The neglected perturbation terms are included in the higher-order terms of the expansion, which are determined from the solution of first-order linear partial differential equations requiring only integrations which are quadratures. The quadratures can be performed rapidly with emerging computer capability, so that real-time approximate optimization can be used to construct the launch guidance law. The application of this technique to flight in three-dimensions is made apparent from the solution presented.

Feeley, Timothy S.↗
NASA NTRS1982

Conversion and matched filter approximations for serial minimum-shift keyed modulation

Serial minimum-shift keyed (MSK) modulation, a technique for generating and detecting MSK using series filtering, is ideally suited for high data rate applications provided the required conversion and matched filters can be closely approximated. Low-pass implementations of these filters as parallel inphase- and quadrature-mixer structures are characterized in this paper in terms of signal-to-noise ratio (SNR) degradation from ideal and envelope deviation. Several hardware implementation techniques utilizing microwave devices or lumped elements are presented. Optimization of parameter values results in realizations whose SNR degradation is less than 0.5 dB at error probabilities of .000001.

Ziemer, R. E.↗
NASA NTRS1979

Precession of the epsilon ring of Uranus

It is noted that the outer and inner boundaries of the epsilon ring of Uranus can be fitted by aligned Keplerian ellipses. Four possible mechanisms for maintaining uniform precession in the epsilon ring are considered: the ring's self-gravity, precession due to a satellite, smooth pressure gradients, and shocklike phenomena. It is proposed that apse alignment is maintained by the self-gravity of the ring. In this case, a ring mass of approximately 5 x 10 to the 18th g and a mean surface density at quadrature of about 25 g/sq cm are estimated.

Goldreich, P.↗
NASA NTRS1995

Solution of acoustic workshop problems by a spectral multidomain method

We use a new staggered grid Chebyshev spectral multidomain method to solve three of the Workshop benchmark problems. The method defines solution unknowns at the nodes of the Chebyshev Gauss quadrature, and the fluxes at the nodes of the Chebyshev Gauss-Lobatto quadrature. The Chebyshev spectral method gives exponentially convergent phase and dissipation errors. The multidomain approximation gives the method flexibility. Using the method, we solve problems in Categories 1 and 5 of the benchmark problems.

Kopriva, Davis A.↗
NASA NTRS1999

Front end for GPS receivers

The front end in GPS receivers has the functions of amplifying, down-converting, filtering and sampling the received signals. In the preferred embodiment, only two operations, A/D conversion and a sum, bring the signal from RF to filtered quadrature baseband samples. After amplification and filtering at RF, the L1 and L2 signals are each sampled at RF at a high selected subharmonic rate. The subharmonic sample rates are approximately 900 MHz for L1 and 982 MHz for L2. With the selected subharmonic sampling, the A/D conversion effectively down-converts the signal from RF to quadrature components at baseband. The resulting sample streams for L1 and L2 are each reduced to a lower rate with a digital filter, which becomes a straight sum in the simplest embodiment. The frequency subsystem can be very simple, only requiring the generation of a single reference frequency (e.g. 20.46 MHz minus a small offset) and the simple multiplication of this reference up to the subharmonic sample rates for L1 and L2. The small offset in the reference frequency serves the dual purpose of providing an advantageous offset in the down-converted carrier frequency and in the final baseband sample rate.

Thomas, Jr., Jess Brooks↗
NASA NTRS1984

The uniform asymptotic swallowtail approximation - Practical methods for oscillating integrals with four coalescing saddle points

Methods that can be used in the numerical implementation of the uniform swallowtail approximation are described. An explicit expression for that approximation is presented to the lowest order, showing that there are three problems which must be overcome in practice before the approximation can be applied to any given problem. It is shown that a recently developed quadrature method can be used for the accurate numerical evaluation of the swallowtail canonical integral and its partial derivatives. Isometric plots of these are presented to illustrate some of their properties. The problem of obtaining the arguments of the swallowtail integral from an analytical function of its argument is considered, describing two methods of solving this problem. The asymptotic evaluation of the butterfly canonical integral is addressed.

Connor, J. N. L.↗
NASA NTRS1992

Design of an essentially non-oscillatory reconstruction procedure in finite-element type meshes

An essentially non oscillatory reconstruction for functions defined on finite element type meshes is designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitary meshes and the reconstruction of a function from its averages in the control volumes surrounding the nodes of the mesh. Concerning the first problem, the behavior of the highest coefficients of two polynomial interpolations of a function that may admit discontinuities of locally regular curves is studied: the Lagrange interpolation and an approximation such that the mean of the polynomial on any control volume is equal to that of the function to be approximated. This enables the best stencil for the approximation to be chosen. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, two methods were studied: one based on an adaptation of the so called reconstruction via deconvolution method to irregular meshes and one that lies on the approximation on the mean as defined above. The first method is conservative up to a quadrature formula and the second one is exactly conservative. The two methods have the expected order of accuracy, but the second one is much less expensive than the first one. Some numerical examples are given which demonstrate the efficiency of the reconstruction.

Abgrall, Remi↗
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.↗