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 667 records · Page 37

An adaptive embedded mesh procedure for leading-edge vortex flows

A procedure for solving the conical Euler equations on an adaptively refined mesh is presented, along with a method for determining which cells to refine. The solution procedure is a central-difference cell-vertex scheme. The adaptation procedure is made up of a parameter on which the refinement decision is based, and a method for choosing a threshold value of the parameter. The refinement parameter is a measure of mesh-convergence, constructed by comparison of locally coarse- and fine-grid solutions. The threshold for the refinement parameter is based on the curvature of the curve relating the number of cells flagged for refinement to the value of the refinement threshold. Results for three test cases are presented. The test problem is that of a delta wing at angle of attack in a supersonic free-stream. The resulting vortices and shocks are captured efficiently by the adaptive code.

Powell, Kenneth G.↗

A Solar Radiation Parameterization for Atmospheric Studies

The solar radiation parameterization (CLIRAD-SW) developed at the Goddard Climate and Radiation Branch for application to atmospheric models are described. It includes the absorption by water vapor, O3, O2, CO2, clouds, and aerosols and the scattering by clouds, aerosols, and gases. Depending upon the nature of absorption, different approaches are applied to different absorbers. In the ultraviolet and visible regions, the spectrum is divided into 8 bands, and single O3 absorption coefficient and Rayleigh scattering coefficient are used for each band. In the infrared, the spectrum is divided into 3 bands, and the k-distribution method is applied for water vapor absorption. The flux reduction due to O2 is derived from a simple function, while the flux reduction due to CO2 is derived from precomputed tables. Cloud single-scattering properties are parameterized, separately for liquid drops and ice, as functions of water amount and effective particle size. A maximum-random approximation is adopted for the overlapping of clouds at different heights. Fluxes are computed using the Delta-Eddington approximation.

Chou, Ming-Dah↗

Assessment of pancreas cells

Pancreatic islets were obtained from guinea pig pancreas by the collagenase method and kept alive in tissue culture prior to further studies. Pancreas cell morphology was studied by standard histochemical techniques using light microscopy. Preparative vertical electrophoresis-levitation of dispersed fetal guinea pig pancreas cells was conducted in phosphate buffer containing a heavy water (D20) gradient which does not cause clumping of cells or alter the osmolarity of the buffers. The faster migrating fractions tended to be enriched in beta-cell content. Alpha and delta cells were found to some degree in most fractions. A histogram showing the cell count distribution is included.

Vanoss, C. J.↗

Numerical simulation of a powered-lift landing, tracking flow features using overset grids, and simulation of high lift devices on a fighter-lift-and-control wing

Attached as appendices to this report are documents describing work performed on the simulation of a landing powered-lift delta wing, the tracking of flow features using overset grids, and the simulation of flaps on the Wright Patterson Lab's fighter-lift-and-control (FLAC) wing. Numerical simulation of a powered-lift landing includes the computation of flow about a delta wing at four fixed heights as well as a simulated landing, in which the delta wing descends toward the ground. Comparison of computed and experimental lift coefficients indicates that the simulations capture the qualitative trends in lift-loss encountered by thrust-vectoring aircraft operating in ground effect. Power spectra of temporal variations of pressure indicate computed vortex shedding frequencies close to the jet exit are in the experimentally observed frequency range; the power spectra of pressure also provide insights into the mechanisms of lift oscillations. Also, a method for using overset grids to track dynamic flow features is described and the method is validated by tracking a moving shock and vortices shed behind a circular cylinder. Finally, Chimera gridding strategies were used to develop pressure coefficient contours for the FLAC wing for a Mach no. of 0.18 and Reynolds no. of 2.5 million.

Chawla, Kalpana↗

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.↗

A method for using a time interval counter to measure frequency stability

It is shown how a commercial time interval counter can be used to measure the relative stability of two signals that are offset in frequency and mixed down to a beat note of about 1 Hz. To avoid the dead-time problem, the counter is set up to read the time interval between each beat note upcrossing and the next pulse of a 10 Hz reference pulse train. The actual upcrossing times are recovered by a simple algorithm whose outputs can be used for computing residuals and Allan variance. A noise floor-test yielded a delta f/f Allan deviation of 1.3 times 10 to the minus 9th power/tau relative to the beat frequency.

Greenhall, C. A.↗

A method for using a time interval counter to measure frequency stability

It is shown how a commercial time interval counter can be used to measure the relative stability of two signals that are offset in frequency and mixed down to a beat note of about 1 Hz. To avoid the dead-time problem, the counter is set up to read the time interval between each beat note upcrossing and the next pulse of a 10 Hz reference pulse train. The actual upcrossing times are recovered by a simple algorithm whose outputs can be used for computing residuals and Allan variance. A noise floor-test yielded a delta f-f Allan deviation of 1.3 times 10 to the minus 9th power/tau relative to the beat frequency.

Greenhall, Charles A.↗

Fast solutions to the steady state compressible and incompressible fluid dynamic equations

For low speed flows the use of the compressible fluid dynamic equations is inefficient. The use of an explicit scheme requires delta t to be bounded by 1/c. However, the physical parameters change over time scales of order 1/u which is much larger. Hence, it is not appropriate to use explicit schemes for very subsonic flows. Implicit schemes are hard to vectorize and frequently do not converge quickly for very subsonic flows. If one is only interested in the steady state then a minor change to an existing code can greatly increase the efficiency of an explicit method. Even when using an implicit method the proposed changes increase the efficiency of the scheme. The Euler equations for low speed flows will be considered first and then incompressible flows. The method is generalized to include viscous effects. Supersonic flow is accelerated by essentially decoupling the equations.

Turkel, E.↗

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models.

Bayliss, A.↗

On accuracy conditions for the numerical computation of waves

The Helmholtz equation (Delta + K(2)n(2))u = f with a variable index of refraction n, and a suitable radiation condition at infinity serves as a model for a wide variety of wave propagation problems. Such problems can be solved numerically by first truncating the given unbounded domain and imposing a suitable outgoing radiation condition on an artificial boundary and then solving the resulting problem on the bounded domain by direct discretization (for example, using a finite element method). In practical applications, the mesh size h and the wave number K, are not independent but are constrained by the accuracy of the desired computation. It will be shown that the number of points per wavelength, measured by (Kh)(-1), is not sufficient to determine the accuracy of a given discretization. For example, the quantity K(3)h(2) is shown to determine the accuracy in the L(2) norm for a second-order discretization method applied to several propagation models.

Bayliss, A.↗

Very long baseline IPS observations of the solar wind speed in the fast polar streams

Observations of intensity scintillation (IPS) with two or more spaced antennas have been widely used to measure the solar wind velocity. Such methods are particularly valuable in regions which spacecraft have not yet penetrated, but they are also very useful in improving the spatial temporal sampling of the solar wind, even in regions where spacecraft data are available. The principle of the measurement is to measure the time delay tau(sub d) between the scintillations observed with an antenna baseline b. The velocity estimate is just V = b/tau(sub d). The error in estimation of the time delay delta tau(sub d) is independent of the baseline length, thus the error in the velocity estimate delta V given by delta(V)/V approximately equals to (delta tau(sub d))/tau(sub d) is inversely proportional to tau(sub d) and hence to b. However the use of a long baseline b has a less obvious advantage; it provides a means for separating fast and slow contributions when both are present in the scattering region. Here we will present recent observations made using the large cylinder antenna at Ooty in the Nilgiri Hills of South India, and one of the 45 m dishes of GMRT near Pune in West-Central India. The baseline of 900 km is, by a considerable margin, the longest ever used for IPS and provides excellent velocity resolution. These results compared with the ULYSSES observations and other IPS measurements made closer to the sun with higher frequency instruments such as EISCAT and the VLBA will provide a precise measure of the velocity profile of the fast north-polar stream.

Rao, A. Pramesh↗

Land Surface Temperature Measurements form EOS MODIS Data

We have developed a physics-based land-surface temperature (LST) algorithm for simultaneously retrieving surface band-averaged emissivities and temperatures from day/night pairs of MODIS (Moderate Resolution Imaging Spectroradiometer) data in seven thermal infrared bands. The set of 14 nonlinear equations in the algorithm is solved with the statistical regression method and the least-squares fit method. This new LST algorithm was tested with simulated MODIS data for 80 sets of band-averaged emissivities calculated from published spectral data of terrestrial materials in wide ranges of atmospheric and surface temperature conditions. Comprehensive sensitivity and error analysis has been made to evaluate the performance of the new LST algorithm and its dependence on variations in surface emissivity and temperature, upon atmospheric conditions, as well as the noise-equivalent temperature difference (NE(Delta)T) and calibration accuracy specifications of the MODIS instrument. In cases with a systematic calibration error of 0.5%, the standard deviations of errors in retrieved surface daytime and nighttime temperatures fall between 0.4-0.5 K over a wide range of surface temperatures for mid-latitude summer conditions. The standard deviations of errors in retrieved emissivities in bands 31 and 32 (in the 10-12.5 micrometer IR spectral window region) are 0.009, and the maximum error in retrieved LST values falls between 2-3 K. Several issues related to the day/night LST algorithm (uncertainties in the day/night registration and in surface emissivity changes caused by dew occurrence, and the cloud cover) have been investigated. The LST algorithms have been validated with MODIS Airborne Simulator (MAS) dada and ground-based measurement data in two field campaigns conducted in Railroad Valley playa, NV in 1995 and 1996. The MODIS LST version 1 software has been delivered.

Wan, Zhengming↗

Preliminary measurement of the airframe noise from an F-106B delta wing aircraft at low flyover speeds

To establish a realistic lower limit for the noise level of advanced supersonic transport aircraft will require knowledge concerning the amount of noise generated by the airframe itself as it moves through the air. The airframe noise level of an F-106B aircraft was determined and was compared to that predicted from an existing empirical relationship. The data were obtained from flyover and static tests conducted to determine the background noise level of the F-106B aircraft. Preliminary results indicate that the spectrum associated with airframe noise was broadband and peaked at a frequency of about 570 hertz. An existing empirical method successfully predicted the frequency where the spectrum peaked. However, the predicted OASPL value of 105 db was considerably greater than the measures value of 83 db.

Burley, R. R.↗

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.↗

Program of Research in Flight Dynamics in the JIAFS, George Washington University at NASA Langley Research Center

The program objectives are fully defined in the original proposal entitled 'Program of Research in Flight Dynamics in GW at NASA Langley Research Center,' which was originated March 20, 1975, and in the renewals of the research program from December 1, 2000 to November 30, 2001. The program in its present form includes three major topics: 1) the improvement of existing methods and development of new methods for wind tunnel and flight test data analysis, 2) the application of these methods to wind tunnel and flight test data obtained from advanced airplanes, 3) the correlation of flight results with wind tunnel measurements, and theoretical predictions. The Principal Investigator of the program is Dr. Vladislav Klein. Three Graduate Research Scholar Assistants (K. G. Mas, M. M. Eissa and N. M. Szyba) also participated in the program. Specific developments in the program during the period Dec. 1, 2001 through Nov. 30, 2002 included: 1) Data analysis of highly swept delta wing aircraft from wind and water tunnel data, and 2) Aerodynamic characteristics of the radio control aircraft from flight test.

Klein, Vladislav↗

Population Pharmacokinetics of Intranasal Scopolamine

Introduction: An intranasal gel dosage formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness (SMS).The bioavailability and pharmacokinetics (PK) was evaluated using data collected in Phase II IND protocols. We reported earlier statistically significant gender differences in PK parameters of INSCOP at a dose level of 0.4 mg. To identify covariates that influence PK parameters of INSCOP, we examined population covariates of INSCOP PK model for 0.4 mg dose. Methods: Plasma scopolamine concentrations versus time data were collected from 20 normal healthy human subjects (11 male/9 female) after a 0.4 mg dose. Phoenix NLME was employed for PK analysis of these data using gender, body weight and age as covariates for model selection. Model selection was based on a likelihood ratio test on the difference of criteria (-2LL). Statistical significance for base model building and individual covariate analysis was set at P less than 0.05{delta(-2LL)=3.84}. Results: A one-compartment pharmacokinetic model with first-order elimination best described INSCOP concentration ]time profiles. Inclusion of gender, body weight and age as covariates individually significantly reduced -2LL by the cut-off value of 3.84(P less than 0.05) when tested against the base model. After the forward stepwise selection and backward elimination steps, gender was selected to add to the final model which had significant influence on absorption rate constant (ka) and the volume of distribution (V) of INSCOP. Conclusion: A population pharmacokinetic model for INSCOP has been identified and gender was a significant contributing covariate for the final model. The volume of distribution and Ka were significantly higher in males than in females which confirm gender-dependent pharmacokinetics of scopolamine after administration of a 0.4 mg dose.

Wu, L.↗

Multi-source Estimates of Land / Ocean Moisture Transport Variability over the Satellite Era

It is widely appreciated that atmospheric transport of water from the world’s oceans is a process key to planetary energy balance as well as Earth's habitability. What is not yet clear is the extent of variability in moisture transports, the relative importance of interdecadal variability versus climate change signals, and importantly, our ability to quantify these changes. This work assesses variations in moisture transport variability during the satellite era (~1980 to present) by comparing several different estimates. (i) The most direct estimate is the vertically integrated flux convergence of moisture from reanalyses which use observed wind and moisture information. (ii) One alternative estimate comes from P-ET over land taken from global hydrologic models constrained with precipitation and near-surface meteorology. Here we use an ensemble of six models. An adjunct to this method is to employ satellite derive ET (e.g., GLEAM or DOLCE). (iii) Complementary to this is E-P over the global oceans derived from satellite estimates of P such as TRMM, GPM and GPCP and SeaFlux V3 or J-OFURO3 estimates of E, all relying heavily upon microwave measurements. Transport between land and oceans must essentially balance at monthly scales, i.e., vanish globally. (iv) a fourth perspective comes from estimate of terrestrial RO + storage rate, delta S. G-RUN Ensemble which uses observed streamflow and P measurements to calibrate a statistical model provides the former while GRACE, GRACE-FO provide total water storage anomalies used to calculate storage rate changes. GRACE REC uses GRACE data to train a precipitation-driven statistical model to extend storage estimates before the GRACE era. (All of these alternatives to reanalysis estimates also consider the small atmospheric column water vapor contribution.) We examine the transport changes from these three different methodologies, their relative accuracies and discuss the origin of their differences. Regional trends in moisture flux divergence and their role in multi-decadal trends are considered. Interannual variability arising in connection with ENSO variability is a dominant signal, driven largely by P changes. Trends since 1980 include reductions in moisture delivery to the western U.S., eastern Brazil, and central Africa with recovery of moisture convergence to the Sahel and parts of eastern North America.

Franklin Robertson↗

Experiments in Aircraft Roll-Yaw Control using Forebody Tangential Blowing

Advantages of flight at high angles of attack include increased maneuverability and lift capabilities. These are beneficial not only for fighter aircraft, but also for future supersonic and hypersonic transport aircraft during take-off and landing. At high angles of attack the aerodynamics of the vehicle are dominated by separation, vortex shedding and possibly vortex breakdown. These phenomena severely compromise the effectiveness of conventional control surfaces. As a result, controlled flight at high angles of attack is not feasible for current aircraft configurations. Alternate means to augment the control of the vehicle at these flight regimes are therefore necessary. The present work investigates the augmentation of an aircraft flight control system by the injection of a thin sheet of air tangentially to the forebody of the vehicle. This method, known as Forebody Tangential Blowing (FTB), has been proposed as an effective means of increasing the controllability of aircraft at high angles of attack. The idea is based on the fact that a small amount of air is sufficient to change the separation lines on the forebody. As a consequence, the strength and position of the vortices are altered causing a change on the aerodynamic loads. Although a very effective actuator, forebody tangential blowing is also highly non-linear which makes its use for aircraft control very difficult. In this work, the feasibility of using FTB to control the roll-yaw motion of a wind tunnel model was demonstrated both through simulations and experimentally. The wind tunnel model used in the experiments consists of a wing-body configuration incorporating a delta wing with 70-degree sweep angle and a cone-cylinder fuselage. The model is equipped with forebody slots through which blowing is applied. There are no movable control surfaces, therefore blowing is the only form of actuation. Experiments were conducted at a nominal angle of attack of 45 degrees. A unique apparatus that constrains the model to two degrees-of-freedom, roll and yaw, was designed and built. The apparatus was used to conduct dynamic experiments which showed that the system was unstable, its natural motion divergent. A model for the unsteady aerodynamic loads was developed based on the basic physics of the flow and results from flow visualization experiments. Parameters of the aerodynamic model were identified from experimental data. The model was validated using data from dynamic experiments. The aerodynamic model completes the equations of motion of the system which were used in the design of control laws using blowing as the only actuator. The unsteady aerodynamic model was implemented as part of the real-time vehicle control system. A control strategy using asymmetric blowing was demonstrated experimentally. A discrete vortex method was developed to help understand the main physics of the flow. The method correctly captures the interactions between forebody and wing vortices. Moreover, the trends in static loads and flow structure are correctly represented. Flow visualization results revealed the vortical structure of the flow to be asymmetric even for symmetric flight conditions. The effects of blowing, and roll and yaw angles on the flow structure were determined. It is shown that superimposing symmetric and asymmetric blowing has a linearizing effect on the actuator characteristics. Transient responses of roll and yaw moments to step input blowing were characterized, and their differences were explained based on the physical mechanisms through which these loads are generated.

Pedreiro, Nelson↗