Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quadrature by expansion”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Relationship Between the Expansion Speed and Radial Speed of CMEs Confirmed Using Quadrature Observations from SOHO and STEREO

It is difficult to measure the true speed of Earth-directed CMEs from a coronagraph along the Sun-Earth line because of the occulting disk. However, the expansion speed (the speed with which the CME appears to spread in the sky plane) can be measured by such coronagraph. In order to convert the expansion speed to radial speed (which is important for space weather applications) one can use empirical relationship between the two that assumes an average width for all CMEs. If we have the width information from quadrature observations, we can confirm the relationship between expansion and radial speeds derived by Gopalswamy et al. (2009, CEAB, 33, 115,2009). The STEREO spacecraft were in quadrature with SOHO (STEREO-A ahead of Earth by 87 and STEREO-B 94 behind Earth) on 2011 February 15, when a fast Earth-directed CME occurred. The CME was observed as a halo by the Large-Angle and Spectrometric Coronagraph (LASCO) on board SOHO. The sky-plane speed was measured by SOHO/LASCO as the expansion speed, while the radial speed was measured by STEREO-A and STEREO-B. In addition, STEREO-A and STEREO-B images measured the width of the CME, which is unknown from Earth view. From the SOHO and STEREO measurements, we confirm the relationship between the expansion speed (Vexp ) and radial speed (Vrad ) derived previously from geometrical considerations (Gopalswamy et al. 2009): Vrad = 1/2 (1 + cot w) Vexp, where w is the half width of the CME. STEREO-B images of the CME, we found that CME had a full width of 75 degrees, so w = 37.5 degrees. This gives the relation as Vrad = 1.15 Vexp. From LASCO observations, we measured Vexp = 897 km/s, so we get the radial speed as 1033 km/s. Direct measurement of radial speed from STEREO gives 945 km/s (STEREO-A) and 1057 km/s (STEREO-B). These numbers are different only by 2.3% and 8.5% (for STEREO-A and STEREO-B, respectively) from the computed value.

Gopalswamy, Nat↗
NASA NTRSDate not supplied

Development of An Improved BRDF Hotspot Model and its Use in VLIDORT to Study the Impact of Atmospheric Scattering on Hotspot Directional Signatures in the Atmosphere

The term “hotspot” refers to the sharp increase of reflectance occurring when incident (solar) and reflected (viewing) directions almost coincide in the backscatter direction. The accurate simulation of hotspot directional signatures is important for many remote sensing applications. The RossThick-LiSparse-Reciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is widely used in radiative transfer simulations, and the hotspot model mostly used is from Maignan- Bréon but it typically requires large values of numerical quadrature and Fourier expansion terms in order to represent the hotspot accurately. To improve its use in atmospheric radiative transfer (RT) model simulations, in this paper we have developed a modified version based on the Maignan-Bréon’s hotspot BRDF model that converge much faster numerically, making it more practical for use in RT models that require Fourier expansion of BRDF to simulate the top-of-atmosphere (TOA) hotspot signatures. Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), we found that reasonable TOA hotspot accuracy can be obtained with just 23 Fourier terms for clear atmospheres, and 63 Fourier terms for atmospheres with aerosol scattering. One advantage of this modified model is that the new hotspot model agrees very well with the original RossThick model away the hotspot region, making it is very convenient to use in the condition with and without hotspot in applications. This model can calculate the amplitude of hot spot accurately, and has been added in the most recent version of VLIDORT. However, there are some difference of this modified model with the original model for scattering angle close the hot spot, and it may not be appropriate for those who need an exact representation of the hot spot angular signature close to hot spot.

Xiaozhen (Shawn) Xiong↗
NASA NTRSDate not supplied

A New BRDF Hotspot Model and its Use in Improving Radiative Transfer Model Efficiency and Application to VLIDORT-based PCRTM Model

The term “hotspot” refers to the sharp increase of reflectance occurring when incident (solar) and reflected (viewing) directions almost coincide in the backscatter direction. The accurate simulation of hotspot directional signatures is important for many remote sensing applications. The RossThick-LiSparse-Reciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is widely used in radiative transfer simulations, and the hotspot model mostly used is from Maignan- Bréon but it typically requires large values of numerical quadrature and Fourier expansion terms in order to represent the hotspot accurately. To improve its use in atmospheric radiative transfer (RT) model simulations, we have developed a modified version based on the Maignan-Bréon’s hotspot BRDF model that converge much faster numerically, making it more practical for use in RT models that require Fourier expansion of BRDF to simulate the top-of-atmosphere (TOA) hotspot signatures. On basis of the vector linearized discrete ordinate radiative transfer model (VLIDORT), we have built a line-by-line based simulation system that can simulate the TOA radiance for atmospheres with 26 gases, different aerosols, dust and clouds, and surface with different BRDF. Using simulated results under diverse atmospheric conditions, we have developed VLIDORT-PCRTM model.

Xiaozhen Xiong↗
DOE OSTI2024

A New Galerkin Quadrature Method Not Requiring a Matrix Inverse

We derive a new Galerkin quadrature (GQ) method for S 𝑛 calculations that differs from the two methods preceding it in that a matrix inverse for an 𝑁 𝑑 × 𝑁 𝑑 matrix, where 𝑁𝑑 is the number of directions in the quadrature set, is no longer required. Galerkin quadrature methods are designed for calculations with highly anisotropic scattering. Such methods are not simply special angular quadratures but also are methods for representing the S 𝑛 scattering source that offers several advantages relative to the standard scattering source representation when highly truncated Legendre cross-section expansions must be used. Galerkin quadrature methods are also useful when the scattering is moderately anisotropic, but the quadrature being used is not sufficiently accurate for the order of the scattering source expansion that is required. Furthermore, we derive the new method and present computational results showing that its performance for two challenging problems is comparable to those of the two GQ methods that preceded it.

Galerkin quadrature↗
DOE OSTI2023

High-Performance Electron Sources: Numerical Methods and Beam Dynamics at the Precision Frontier

Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem

43 PARTICLE ACCELERATORS↗
NASA NTRS1975

Averaged initial Cartesian coordinates for long lifetime satellite studies

A set of initial Cartesian coordinates, which are free of ambiguities and resonance singularities, is developed to study satellite mission requirements and dispersions over long lifetimes. The method outlined herein possesses two distinct advantages over most other averaging procedures. First, the averaging is carried out numerically using Gaussian quadratures, thus avoiding tedious expansions and the resulting resonances for critical inclinations, etc. Secondly, by using the initial rectangular Cartesian coordinates, conventional, existing acceleration perturbation routines can be absorbed into the program without further modifications, thus making the method easily adaptable to the addition of new perturbation effects. The averaged nonlinear differential equations are integrated by means of a Runge Kutta method. A typical step size of several orbits permits rapid integration of long lifetime orbits in a short computing time.

Pines, S.↗
NASA NTRS1976

Rapid computation of the Voigt profile

Computational procedures for evaluating the Voigt profile function with maximum relative error about one part in ten thousand are discussed. The computational region is split up into four subregions, each of which is worked on with a different computational technique (Chebyshev expansion, continued fraction expansion, 2-point Gauss-Hermite quadrature, and 4-point Gauss-Hermite quadrature). The overall procedure is designed for line-by-line transmittance calculations and similar applications, and an efficient FORTRAN IV subprogram is outlined in an appendix.

Drayson, S. R.↗
NASA NTRS2024

Quantification of Uncertainty and Risk Sensitivity for Safety of Emerging Operations

The growing need to develop and deploy small unmanned aerial vehicles (sUAVs) for various applications in the airspace necessitates reliable tools to accurately predict the flight trajectories of the sUAVs. The knowledge of the predicted trajectories help decision makers anticipate potential conflict, assess the risk, and take appropriate risk mitigation actions. In addition, uncertainties in vehicle models, weather, and controller action further highlights the need for reliable prediction tools. In this project, the application of mixed sparse grid-based quadrature and generalized polynomial chaos(gPC) expansion method for uncertainty quantification and collision assessment in air traffic consisting of fixed-wing small unmanned aerial vehicles (sUAV) was studied. From the results obtained, it can be concluded that this provides a reliable framework to carry out quantitative conflict assessment in an unmanned air traffic, which when employed, can improve the functionalities of the unmanned traffic management system. It was observed that the results from the gPC expansion framework developed in the project can be utilized to conduct rapid probabilistic collision assessment for near real-time unmanned traffic management in the airspace. From the vehicle models, position updates, and wind-field data, a priori gPC based 3-σcon-fidence ellipses can provide estimates of potential conflict at some future instants. The computational costs scaled linearly when the uncertain inputs were fewer. Further, the largest allowable distribution of para-metric uncertainties that leads to the smallest risk of collision in traffic of small unmanned aerial vehicles could be calculated. The time of closest approach between two sUAVs can be established paving way for development of proactive mitigation strategies. The separation between the sUAVs was found to be most significantly affected by uncertainties in the maximum available thrusts, zero-lift drag coefficients, and wing planform areas of the sUAVs. The study of uncertain wind-fields indicated that a heterogeneous traffic mix resulted in an increased probability of conflict. Increased measurement update rate reduced the uncertain-ties in the trajectories of the vehicles, further reducing the probability of conflict but rapid updates of all vehicles in the airspace poses a stringent communication limitation. The gPC framework also provided the means to analyze vehicle impact (crash region) due to loss of control resulting from actuator failure in sUAS traffic, essentially to predict impact and crash zones for representative vehicles. The predicted regions when compared with non-participant density, provides a means to develop an early mitigation strategy, should the sUAV detect an imminent actuator failure.

Rajnish Bhusal↗
NASA NTRS1985

Unsteady lifting-line theory as a singular-perturbation problem

Unsteady lifting-line theory is developed for a wing of large aspect ratio oscillating at low frequency in inviscid incompressible flow. The wing is assumed to have a rigid chord but a flexible span. Use of the method of matched asymptotic expansions reduces the problem from a singular integral equation to quadrature. The pressure field and airloads, for a prescribed wing shape and motion, are obtained in closed form as expansions in inverse aspect ratio. A rigorous definition of unsteady induced downwash is also obtained. Numerical calculations are presented for an elliptic wing in pitch and heave; compared with numerical lifting-surface theory, computation time is reduced significantly. The present work also identifies and resolves errors in the unsteady lifting line theory of James (1975), and points out a limitation in that of Van Holten (1975, 1976, 1977).

Ahmadi, A.↗
NASA NTRS2000

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗
NASA NTRS1986

Numerical quadrature methods for integrals of singular periodic functions and their application to singular and weakly singular integral equations

High accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Such periodic equations are used in the solution of planar elliptic boundary value problems, elasticity, potential theory, conformal mapping, boundary element methods, free surface flows, etc. The use of the quadrature methods is demonstrated with numerical examples.

Sidi, A.↗
NASA NTRS1988

Quadrature methods for periodic singular and weakly singular Fredholm integral equations

High-accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are subsequently used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Throughout the development the periodic nature of the problem plays a crucial role. Such periodic equations are used in the solution of planar elliptic boundary value problems such as those that arise in elasticity, potential theory, conformal mapping, and free surface flows. The use of the quadrature methods is demonstrated with numerical examples.

Sidi, Avram↗
NASA NTRS1992

A doubly averaging method for third body perturbations in planet equator coordinates

The first order doubly averaged potential due to third-body gravity is derived in any arbitrary coordinates. The equations of motion are nonsingular at zero eccentricity. The derivation uses a recursive method which allows easy expansion to higher order terms. Instead of using analytical quadrature to obtain the doubly averaged potential, the method presented in this paper simply eliminates the mean anomaly of the perturbed and perturbing bodies by inspection of the recursive formulation. The derivatives of the orbital elements can be numerically integrated rapidly. When a planet equator coordinate system is used, they can be added directly to the derivatives due to gravity harmonics without any coordinate transformation. The method is applied to various high altitude missions. The results are compared with a high precision numerical integration method and are found to provide excellent agreement.

Hwok, Johnny H.↗
NASA NTRS1999

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations that Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + integral(0 to 1) g(x,t) F(t, y(t)) dt, 0 less than or equal to x less than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integrals equations arise, e.g., when one applies Green's function techniques to nonlinear two-point boundary value problems of the form U''(x) = f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and g(l) = y(sub 1), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trapezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram↗
NASA NTRS1981

An asymptotic unsteady lifting-line theory with energetics and optimum motion of thrust-producing lifting surfaces

A low frequency unsteady lifting-line theory is developed for a harmonically oscillating wing of large aspect ratio. The wing is assumed to be chordwise rigid but completely flexible in the span direction. The theory is developed by use of the method of matched asymptotic expansions which reduces the problem from a singular integral equation to quadrature. The wing displacements are prescribed and the pressure field, airloads, and unsteady induced downwash are obtained in closed form. The influence of reduced frequency, aspect ratio, planform shape, and mode of oscillation on wing aerodynamics is demonstrated through numerical examples. Compared with lifting-surface theory, computation time is reduced significantly. Using the present theory, the energetic quantities associated with the propulsive performance of a finite wing oscillating in combined pitch and heave are obtained in closed form. Numerical examples are presented for an elliptic wing.

Ahmadi, A. R.↗