Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “delta method”

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 487 records · Page 27

Study of the vortex conditions of wings with large sweepback by extrapolation of the Jones method

The pockets of separation originating on the leading edges are surrounded by vortex sheets. Their configuration and intensity were determined by four conditions with the JONES approximation, which is itself corrected by a simple logic. Field pressures and stresses were computed for different cases and are compared with test results (pure deltas, swallow tails, truncations, strakes, ducks, fuselage).

Hirsch, P.↗

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗

Calculation of supersonic viscous flow over delta wings with sharp subsonic leading edges

Two complementary procedures were developed to calculate the viscous supersonic flow over conical shapes at large angles of attack, with application to cones and delta wings. In the first approach the flow is assumed to be conical and the governing equations are solved at a given Reynolds number with a time-marching explicit finite-difference algorithm. In the second method the parabolized Navier-Stokes equations are solved with a space-marching implicit noniterative finite-difference algorithm. This latter approach is not restricted to conical shapes and provides a large improvement in computational efficiency over published methods. Results from the two procedures agree very well with each other and with available experimental data.

Vigneron, Y. C.↗

Calculation of supersonic viscous flow over delta wings with sharp subsonic leading edges

Two complementary procedures have been developed to calculate the viscous supersonic flow over conical shapes at large angles of attack, with application to cones and delta wings. In the first approach the flow is assumed to be conical and the governing equations are solved at a given Reynolds number with a time-marching explicit finite-difference algorithm. In the second method the parabolized Navier-Stokes equations are solved with a space-marching implicit noniterative finite-difference algorithm. This latter approach is not restricted to conical shapes and provides a large improvement in computational efficiency over published methods. Results from the two procedures agree very well with each other and with available experimental data.

Vigneron, Y. C.↗

Heavy doping effects in high efficiency silicon solar cells

The temperature dependence of the emitter saturation current for bipolar devices was studied by varying the surface recombination velocity at the emitter surface. From this dependence, the value was derived for bandgap narrowing that is in better agreement with other determinations that were obtained from the temperature dependence measure on devices with ohmic contacts. Results of the first direct measurement of the minority-carrier transit time in a transparent heavily doped emitter layer were reported. The value was obtained by a high-frequency conductance method recently developed and used for doped Si. Experimental evidence is presented for significantly greater charge storage in highly excited silicon near room temperature than conventional theory would predict. These data are compared with various data for delta E sub G in heavily doped silicon.

Lindholm, F. A.↗

Numerical simulation of the vortical flow over a round-edged double-delta wing

Simulations of three-dimensional vortical flows over a thin double-delta wing with an aspect ratio of 2.05 have been performed. Steady-state solutions to the unsteady incompressible Navier-Stokes equations are obtained using a new implicit upwind-relaxation finite-difference scheme. The method is second-order accurate spatially and naturally dissipative. Numerical results indicate that key features of both vortical interaction and vortex breakdown are successfully simulated. The computed lift coefficients and lateral trajectories of the vortical cores are in good agreement with the experimental data.

Hsu, C.-H.↗

On the Dome Effect of Flux Radiometers to Radiative Forcing

Since the introduction of thermopile, pyranometers (solar, e.g., 0.3-3.0 micrometers) and pyrgeometers (terrestrial, e.g., 4-50 micrometers) have become instruments commonly used for measuring the broadband hemispherical irradiances at the surface in a long-term, monitoring mode for decades. These commercially available radiometers have been manufactured in several countries such as from the United States, Asia, and Europe, and are generally reliable and economical. These worldwide distributions of surface measurements become even more important in the era of Earth remote sensing in studying climate forcing. However, recent studies from field campaigns have pointed out that erroneous factors (e.g., temperature gradients between the filter dome and detector, emissivity of the thermopile) are responsible for the unacceptable level of uncertainty (e.g., 10-20 W m (exp -2)). Using a newly developed instrument of Quantum Well Infrared Photodetector (QWIP), we have characterized the brightness temperature fields of pyranometers and pyrgeometers under various sky conditions. The QWIP is based on the superlattice (GaAs/AlGaAs) technology and has a noise equivalent temperature (NE delta T) less than 0.1 K. The quality of pyranometer and pyrgeometer measurements can be improved largely by applying proper knowledge of the thermal parameters affecting the operation of the thermopile systems. For example, we show a method to determine the "dome factor" (the longwave emission divided by the longwave transmission of a pyrgeometer dome) from field measurements. The results show, and are verified independently by the QWIP, that our dome factors of 0.59 and 0.90 are much smaller than the value of 4.0 assumed by the WMO (World Meteorological Organization). Data correction procedure and algorithm will be presented and discussed.

Tsay, S.-C.↗

On the Dome Effect of Flux Radiometers to Radiative Forcing

Since the introduction of thermopile, pyranometers (solar, e.g., 0.3-3.0 microns) and pyrgeometers (terrestrial, e.g., 4-50 microns) have become instruments commonly used for measuring the broadband hemispherical irradiances at the surface in a long-term, monitoring mode for decades. These commercially available radiometers have been manufactured in several countries such as from the United States, Asia, and Europe, and are generally reliable and economical. These worldwide distributions of surface measurements become even more important in the era of Earth remote sensing in studying climate forcing. However, recent studies from field campaigns have pointed out that erroneous factors (e.g., temperature gradients between the filter dome and detector, emissivity of the thermopile) are responsible for the unacceptable level of uncertainty (e.g., 10-20 W/square Meter). Using a newly developed instrument of Quantum Well Infrared Photodetector (QWIP), we have characterized the brightness temperature fields of pyranometers and pyrgeometers under various sky conditions. The QWIP is based on the superlattice (GaAs/AlGaAs) technology and has a noise equivalent temperature (NE delta T) less than 0.1 K. The quality of pyranometer and pyrgeometer measurements can be improved largely by applying proper knowledge of the thermal parameters affecting the operation of the thermopile systems. For example, we show a method to determine the "dome factor" (the longwave emission divided by the longwave transmission of a pyrgeometer dome) from field measurements. The results show, and are verified independently by the QWIP, that our dome factors of 0.59 and 0.90 are much smaller than the value of 4.0 assumed by the WMO. Data correction procedure and algorithm will be presented and discussed.

Tsay, S.-C.↗

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'propagate-burn-propagate' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons between it and three methods using modified linearized covariance propagation are made. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. In the second method, the maneuver is not modeled but the uncertainty in its magnitude is accounted for by the inclusion of a process noise matrix. In the third method, which is essentially a hybrid of the first two, the nominal portion of the maneuver is included via the state transition matrix while a process noise matrix is used to account for the magnitude uncertainty. Since this method also correctly accounts for the presence of the maneuver in the nominal orbit, it is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, Pc. However, applications for the two other methods exist and are briefly discussed. Despite the fact that the process model ('propagate-burn-propagate') that was studied was very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗

Robustness Analysis of Integrated LPV-FDI Filters and LTI-FTC System for a Transport Aircraft

This paper proposes an analysis framework for robustness analysis of a nonlinear dynamics system that can be represented by a polynomial linear parameter varying (PLPV) system with constant bounded uncertainty. The proposed analysis framework contains three key tools: 1) a function substitution method which can convert a nonlinear system in polynomial form into a PLPV system, 2) a matrix-based linear fractional transformation (LFT) modeling approach, which can convert a PLPV system into an LFT system with the delta block that includes key uncertainty and scheduling parameters, 3) micro-analysis, which is a well known robust analysis tool for linear systems. The proposed analysis framework is applied to evaluating the performance of the LPV-fault detection and isolation (FDI) filters of the closed-loop system of a transport aircraft in the presence of unmodeled actuator dynamics and sensor gain uncertainty. The robustness analysis results are compared with nonlinear time simulations.

Khong, Thuan H.↗

Estimation of End-of-Life Solar Absorptivity for Complex, Multi-Phased Space Missions

Solar absorptivity degradation over time is an important contributor to the thermal performance of a spacecraft at the end of its life. When environments are complex and include many different factors known to impact solar absorptivity degradation, such as UV, charged particles, contamination, or lunar dust coverage, combining data from several different studies can be done to estimate total degradation. While the most conservative way is to stack delta absorptivity from each factor to come up with a total delta, this may be overly conservative for long or complex missions with multiple mission phases. As absorptivity increases the rate of degradation decreases, and synergistic effects need to be taken into account. In this report an alternate method is proposed using Gateway materials as an example. This method determines which environmental factor will be dominant at what time (such as contamination or charged particles) and uses that degradation rate for that mission phase. For the next mission phase, the new absorptivity can be used instead of the original BOL absorptivity to determine degradation rate. This approach prevents overestimating optical degradation when considering degradation from interconnected sources.

Brandon Michael Hoffmann↗

Estimation of End-of-Life Solar Absorptivity for Complex, Multi-Phased Space Missions

Solar absorptivity degradation over time is an important contributor to the thermal performance of a spacecraft at the end of its life. When environments are complex and include many different factors known to impact solar absorptivity degradation, such as UV, charged particles, contamination, or lunar dust coverage, combining data from several different studies can be done to estimate total degradation. While the most conservative way is to stack delta absorptivity from each factor to come up with a total delta, this may be overly conservative for long or complex missions with multiple mission phases. As absorptivity increases the rate of degradation decreases, and synergistic effects need to be taken into account. In this report an alternate method is proposed using Gateway materials as an example. This method determines which environmental factor will be dominant at what time (such as contamination or charged particles) and uses that degradation rate for that mission phase. For the next mission phase, the new absorptivity can be used instead of the original BOL absorptivity to determine degradation rate. This approach prevents overestimating optical degradation when considering degradation from interconnected sources.

thermo-optical properties↗

Flutter of pairs of aerodynamically interfering delta wings.

To examine the effect on flutter of the aerodynamic interference between pairs of closely spaced delta wings, several structurally uncoupled 1/80th-scale models were studied by experiment and analysis. Flutter test boundaries obtained in NASA Langley's 26-in. transonic blowdown wind tunnel were compared with subsonic analytical results generated using the doublet lattice method. Trends for several combinations of vertical and longitudinal wing separation were determined, showing flutter speed significantly affected in the closely spaced configurations. A new flutter mechanism coupling one wing's first bending mode with the other wing's first torsion mode was predicted and observed.

Chipman, R. R.↗

On the application and extension of Harten's high resolution scheme

Extensions of a second order high resolution explicit method for the numerical computation of weak solutions of one dimensonal hyperbolic conservation laws are discussed. The main objectives were (1) to examine the shock resoluton of Harten's method for a two dimensional shock reflection problem, (2) to study the use of a high resolution scheme as a post-processor to an approximate steady state solution, and (3) to construct an implicit in the delta-form using Harten's scheme for the explicit operator and a simplified iteration matrix for the implicit operator.

Yee, H. C.↗

An efficient coordinate transformation technique for unsteady, transonic aerodynamic analysis of low aspect-ratio wings

An efficient coordinate transformation technique is presented for constructing grids for unsteady, transonic aerodynamic computations for delta-type wings. The original shearing transformation yielded computations that were numerically unstable and this paper discusses the sources of those instabilities. The new shearing transformation yields computations that are stable, fast, and accurate. Comparisons of those two methods are shown for the flow over the F5 wing that demonstrate the new stability. Also, comparisons are made with experimental data that demonstrate the accuracy of the new method. The computations were made by using a time-accurate, finite-difference, alternating-direction-implicit (ADI) algorithm for the transonic small-disturbance potential equation.

Guruswamy, G. P.↗

Stabilized determination of geopotential coefficients by the mixed hom-BLUP approach

For the determination of geopotential coefficients, data can be used from rather different sources, e.g., satellite tracking, gravimetry, or altimetry. As each data type is particularly sensitive to certain wavelengths of the spherical harmonic coefficients it is of essential importance how they are treated in a combination solution. For example the longer wavelengths are well described by the coefficients of a model derived by satellite tracking, while other observation types such as gravity anomalies, delta g, and geoid heights, N, from altimetry contain only poor information for these long wavelengths. Therefore, the lower coefficients of the satellite model should be treated as being superior in the combination. In the combination a new method is presented which turns out to be highly suitable for this purpose due to its great flexibility combined with robustness.

Middel, B.↗

Virtual zone Navier-Stokes computations for oscillating control surfaces

A new zoning method called 'virtual zones' has been developed for application to an unsteady finite difference Navier-Stokes code. The virtual zoning method simplifies the zoning and gridding of complex configurations for use with patched multi-zone flow codes. An existing interpolation method has been extensively modified to bring the run time for the interpolation procedure down to the same level as for the flow solver. Unsteady Navier-Stokes computations have been performed for transonic flow over a clipped delta wing with an oscillating control surface. The computed unsteady pressure and response characteristics of the control-surface motion compare well with experimental data.

Klopfer, G. H.↗

Optimization of a GO2/GH2 Impinging Injector Element

An existing injector optimization methodology, method i, is used to investigate optimal design points for a GO2/GH2 impinging injector element. The impinging element, an F-O-F triplet, is optimized in terms of such relevant design variables as fuel pressure drop, DELTA-P(sub f), oxidizer pressure drop, DELTA-P(sub o), combustor length, L(sub comb), and impingement angle, alpha, for a given mixture ratio and chamber pressure.

Tucker, P. Kevin↗