Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smoothing methods”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
NASA NTRS2001

Rocket Engine Thrust Chamber Assembly

A thrust chamber assembly for liquid fueled rocket engines and the method of making it wherein a two-piece mandrel having the configuration of an assembly having a combustion chamber portion connected to a nozzle portion through a throat portion is wrapped with a silica tape saturated with a phenolic resin, the tape extending along the mandrel and covering the combustion chamber portion of the mandrel to the throat portion. The width of the tape is positioned at an angle of 30 to 50 deg. to the axis of the mandrel such that one edge of the tape contacts the mandrel while the other edge is spaced from the mandrel. The phenolic in the tape is cured and the end of the wrap is machined to provide a frusto-conical surface extending at an angle of 15 to 30 deg. with respect to the axis of the mandrel for starting a second wrap on the mandrel to cover the throat portion. The remainder of the mandrel is wrapped with a third silica tape having its width positioned at a angle of 5 to 20 deg. from the axis of the mandrel. The resin in the third tape is cured and the assembly is machined to provide a smooth outer surface. The entire assembly is then wrapped with a tow of graphite fibers wetted with an epoxy resin and, after the epoxy resin is cured, the graphite is machined to final dimensions.

Cornelius, Charles S.↗
NASA NTRS2009

Gender-related Changes in Dorsal Hand and Foot Vein Function Following 60 Days of Head Down Bed Rest

It is well known that female astronauts are more likely to experience post-flight orthostatic hypotension and presyncope compared to male astronauts. It has been suggested that the disproportionally higher incidence of presyncope (83% of female vs. 20% male crewmembers) may be due to sex-related differences in vascular function between the upper and lower limbs. However, much of this evidence is specific to changes in resistance vessels. Given that more than 70% of the circulating blood volume resides in compliance vessels, it is conceivable that even small changes in venous function may contribute to post-flight orthostatic hypotension. In spite of this, little is currently known regarding the influence of microgravity exposure on venous function between males and females. PURPOSE: To determine the influence of 60 days of HDBR on dorsal foot and hand vein function between healthy males (M) and females (F). METHODS: Using 2-D ultrasound, dorsal hand and foot vein diameter responses to intravenous infusions phenylephrine (PE), acetylcholine (ACh), and nitroglycerine (NTG) were determined in 26 adults; 10 females (age:37 +/- 2 yr ) and 16 males (age:34 +/- 2 yr ). Changes in venous function were calculated as the difference between diameter at baseline and following each venoactive drug. Differences in venous function between limb and sexes across HDBR were determined using mixed-effects linear regression. RESULTS: In response to 60 days of HDBR, the change in venousconstrictor response to PE in the dorsal hand veins was not significantly different between M and F. Interestingly, the change in constrictor response in the dorsal foot veins (compared to pre HDBR) was approximately 30% greater in the F, whereas the constrictor response was approximately 45% less in the M (p=0.026). HDBR had no influence on the change in dilator response to ACh, or NTG between M and F and between vascular beds. CONCLUSION: These results demonstrate that 60 days of HDBR contributes to sex-related differences in venous function. Importantly, the changes are specific to dorsal foot veins and do not appear to be related to changes in endothelial or smooth muscle function. Previous data from our laboratory indicate that females have an impaired sympathetic response following HDBR. However, the results from the current investigation imply that F maintain the ability to respond to sympathetic stimulation but may not be able to adequately activate a proper response following HDBR.

Westby, Christian M.↗
NASA NTRS2002

Spectral Re-Growth Reduction for CCSDS 8-D 8-PSK TCM

This report presents a study on the CCSDS recommended 8-dimensional 8 PSK Trellis Coded Modulation (TCM) scheme. The important steps of the CCSDS scheme include: conversion of serial data into parallel form, differential encoding, convolutional encoding, constellation mapping, and filtering the 8-PSK symbols using the square root raised cosine (SRRC) pulses. The last step, namely the filtering of the 8 PSK symbols using SRRC pulses, significantly affects the bandwidth of the signal. If a nonlinear power amplifier is used, the SRRC filtered signal creates spectral regrowth. The purpose of this report is to investigate a technique, called the smooth phase interpolated keying (SPIK), that can provide an alternative to SRRC filtering so that good spectral as well as power efficiencies can be obtained with the CCSDS encoder. The results of this study show that the CCSDS encoder does not affect the spectral shape of the SRRC filtered signal or the SPIK signal. When a nonlinear traveling wave tube amplifier (TWTA) is used, the spectral performance of the SRRC signal degrades significantly while the spectral performance of SPIK remains unaffected. The degrading effect of a nonlinear solid state power amplifier (SSPA) on SRRC is found to be less than that due to a nonlinear TWTA. However, in both cases, the spectral performance of the SRRC modulated signal is worse than that of the SPIK signal. The bit error rate (BER) performance of the SRRC signal in a linear amplifier environment is about 2.5 dB better than that of the SPIK signal when both the receivers use algorithms of similar complexity. In a nonlinear TWTA environment, the SRRC signal requires accurate phase tracking since the TWTA introduces additional phase distortion. This problem does not arise with SPIK signal due to its constant envelope property. When a nonlinear amplifier is used, the SRRC method loses nearly 1 dB in the bit error rate performance. The SPIK signal does not lose any performance. Thus the performance gap between SRRC and SPIK reduces. The BER performance of SPIK can be improved even further by using a more optimal receiver. A similar optimal receiver for SRRC is quite complex since the amplifier distorts the pulse shape. However, this requires further investigation and is not covered in this report.

Borah, Deva K.↗
NASA NTRS2007

Fan Rig Noise Spectral Correction for NASA 9'x 15' Low-Speed Wind Tunnel

Aircraft engine noise research and development depends on the ability to study and predict the noise created by each engine component in isolation. Fan noise testing, however, requires a significant support system including a drive mechanism to turn the fan, a device to smooth the flow into the fan, and a stand to raise the fan off the ground each of which has the potential to create its own noise. A methodology was therefore developed to improve the data quality for the 9x15 Low Speed Wind Tunnel (LSWT) at the NASA Glenn Research Center that identifies three noise sources: fan noise, jet noise, and rig noise. The jet noise and rig noise was then measured by mounting a scale model of the 9x15 LSWT setup in a jet rig to simulate everything except the rotating machinery that characterizes fan noise. The data showed that the spectra measured in the LSWT has a strong rig noise component at frequencies as high as 3 kHz depending on the fan and speed. The jet noise was determined to be significantly lower than the rig noise. A mathematical model for the rig noise was then developed using a multi-dimensional least squares fit to the rig noise data. This allows the rig noise to be subtracted or removed, depending on the amplitude of the rig noise relative to the fan noise, at any given frequency, observer angle, or nozzle pressure ratio. The impact of isolating the fan noise with this method on spectra, overall power level (OAPWL), and Effective Perceived Noise Level (EPNL) is studied.

Schifer, Nick↗
NASA NTRS1994

Experimental studies of transpiration cooling with shock interaction in hypersonic flow, part B

This report describes the result of experimental studies conducted to examine the effects of the impingement of an oblique shock on the flowfield and surface characteristics of a transpiration-cooled wall in turbulent hypersonic flow. The principal objective of this work was to determine whether the interaction between the oblique shock and the low-momentum region of the transpiration-cooled boundary layer created a highly distorted flowfield and resulted in a significant reduction in the cooling effectiveness of the transpiration-cooled surface. As a part of this program, we also sought to determine the effectiveness of transpiration cooling with nitrogen and helium injectants for a wide range of blowing rates under constant-pressure conditions in the absence of shock interaction. This experimental program was conducted in the Calspan 48-Inch Shock Tunnel at nominal Mach numbers of 6 and 8, for a Reynolds number of 7.5 x 10(exp 6). For these test conditions, we obtained fully turbulent boundary layers upstream of the interaction regions over the transpiration-cooled segment of the flat plate. The experimental program was conducted in two phases. In the first phase, we examined the effects of mass-addition level and coolant properties on the cooling effectiveness of transpiration-cooled surfaces in the absence of shock interaction. In the second phase of the program, we examined the effects of oblique shock impingement on the flowfield and surface characteristics of a transpiration-cooled surface. The studies were conducted for a range of shock strengths with nitrogen and helium coolants to examine how the distribution of heat transfer and pressure and the characteristics of the flowfield in the interaction region varied with shock strength and the level of mass addition from the transpiration-cooled section of the model. The effects of the distribution of the blowing rate along the interaction regions were also examined for a range of blowing rates through the transpiration-cooled panels. The regions of shockwave/boundary layer interaction examined in these studies were induced by oblique shocks generated with a sharp, flat plate, inclined to the freestream at angles of 5 degrees, 7.5 degrees, and 10 degrees. It was found that, in the absence of an incident shock, transpiration cooling was a very effective method for reducing both the heat transfer and the skin friction loads on the surface. The helium coolant was found to be significantly more effective than nitrogen, because of its low molecular weight and high specific heat. The studies of shock-wave/transpiration-cooled surface interaction demonstrated that the interaction region between the incident shock and the low-momentum transpiration-cooled boundary layer did not result in a significant increase in the size of attached or separated interaction regions, and did not result in significant flowfield distortions above the interaction region. The increase in heating downstream of the shock-impingement point could easily be reduced to the values without shock impingement by a relatively small increase in the transpiration cooling in this region. Surprisingly, this increase in cooling rate did not result in a significant increase in size of the region ahead of the incident shock or create a significantly enlarged interaction region with a resultant increase in the distortion level in the inviscid flow. Thus, transpiration cooling appears to be a very effective technique to cool the internal surfaces of scramjet engines, where shocks in the engine would induce large local increases in wall heating and create viscous/inviscid interactions that could significantly disturb the smooth flow through the combustor. However, if hydrogen is used as the coolant, burning upstream of shock impingement might result in localized hot spots. Clearly, further research is needed in this area.

Holden, Michael S.↗
NASA NTRS2004

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗
NASA NTRS1997

Teamwork in the Terminal Area: Organizational Issues and Solutions

Dynamic growth and technology advances in commercial aviation have turned the terminal area into a complex, multi-organization workplace which requires the smooth coordination of many operational teams. In addition to pilots, cabin crew, air traffic controllers, and dispatch (who nominally work together throughout a flight), surface operations additionally involve local, ground and ramp controllers, ramp agents, maintenance, dozens of service contractors, and any number of teams who are responsible for airport operations. Under abnormal or emergency conditions, even more teams become actively involved. In order to accommodate growth and to meet productivity and safety challenges, numerous changes are being made in surface operations. Unfortunately, it is often the case that changes in technologies, organizational roles, procedures, and training are developed and implemented in isolated and piecemeal fashion without regard to cross organizational impact. Thus, there is a need for evaluation methodologies which assure integrated system safety for all organizations. Such methodologies should aid the understanding of how organizations work together and how changes in one domain affects the next. In this study, we develop one approach toward addressing these organizational issues. Examples of surface operations in abnormal situations are examined in regard to their impact on personnel in the terminal area. Timelines are given for the responses to incidents, along with the necessary communication links, the specific roles that members of terminal teams have, and any overlapping responsibilities. Suggestions to improve cross-operational teamwork are given. Methods of graphic representation are explored, both in regards to human links and access to information. The outcome of such an approach should enhance the understanding which is critical for resolving organizational conflicts and maximizing system effectiveness.

Parke, Bonny K.↗
NASA NTRS1999

Lidar Altimeter Measurements of Canopy Structure: Methods and Validation for Closed Canopy, Broadleaf Forests

Lidar altimeter observations of vegetated landscapes provide a time-resolved measure of laser pulse backscatter energy from canopy surfaces and the underlying ground. Airborne lidar altimeter data was acquired using the Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) for a successional sequence of four, closed-canopy, deciduous forest stands in eastern Maryland. The four stands were selected so as to include a range of canopy structures of importance to forest ecosystem function, including variation in the height and roughness of the outer-most canopy surface and the vertical organization of canopy stories and gaps. The character of the SLICER backscatter signal is described and a method is developed that accounts for occlusion of the laser energy by canopy surfaces, transforming the backscatter signal to a canopy height profile (CHP) that quantitatively represents the relative vertical distribution of canopy surface area. The transformation applies an increased weighting to the backscatter amplitude as a function of closure through the canopy and assumes a horizontally random distribution of the canopy components. SLICER CHPs, averaged over areas of overlap where lidar ground tracks intersect, are shown to be highly reproducible. CHP transects across the four stands reveal spatial variations in vegetation, at the scale of the individual 10 m diameter laser footprints, within and between stands. Averaged SLICER CHPs are compared to analogous height profile results derived from ground-based sightings to plant intercepts measured on plots within the four stands. Tbe plots were located on the segments of the lidar ground tracks from which averaged SLICER CHPs were derived, and the ground observations were acquired within two weeks of the SLICER data acquisition to minimize temporal change. The differences in canopy structure between the four stands is similarly described by the SLICER and ground-based CHP results, however a Chi-square test of similarity documents differences that are statistically significant. The differences are discussed in terms of measurement properties that define the smoothness of the resulting CHPs and Lidar Altimeter Measurements of Canopy Structure - Harding et al. canopy properties that may vertically bias the CHP representations of canopy structure. The statistical differences are most likely due to the more noisy character of the ground-based CHPs, especially high in the canopy where ground-based sightings are rare resulting in an underestimate of canopy surface area and height, and to departures from the assumption of horizontal randomness which bias the CHPs toward the observer (upward for SLICER and downward for ground-based CHPs). The results demonstrate that the SLICER observations reliably provide a measure of canopy structure that reveals ecologically interesting structural variations such as those characterizing a successional sequence of closed-canopy, broadleaf forest stands.

Harding, D. J.↗
NASA NTRS1999

Carbon Fiber Reinforced/Silicon Carbide Turbine Blisk Testing in the SIMPLEX Turbopump

A program designed to implement a ceramic matrix composite integrally bladed disk (blisk) into rocket engine style turbomachinery has successfully completed testing. The Marshall Space Flight Center (MSFC) program, utilizing the MSFC turbomachinery design, analysis, and testing capabilities along with materials development capabilities from both Glenn Research Center (GRC) and MSFC, has tested two carbon fiber reinforced silicon carbide blisks in the Simplex Turbopump at MSFC's Test Stand 500. One blisk contained a polar woven fiber preform, while the second blisk tested utilized a quasi-isotropic preform. Vhile earlier papers have chronicled the program's design, material testing, and torque testing efforts, this paper focuses on the testing of the blisks in the Simplex turbopump. Emphasis will be placed on the actual condition of the blisks before and after the testing test program design methodology, and conclusions that can be drawn from the test data and blisk final conditions. The program performed three separate test series. The first series was needed to validate that the Simplex turbopump was correctly re-built following a major incident to the turbopump. The turbopump had two major differences from the original design. The most obvious difference was the sleeve required throughout the bore of the main housing. The second major difference was modifications to the pump diffuser to improve performance. Several areas were burnt during the incident and were either repaired by weld repair (pump inlet housing) or simply smoothed out (turbine nozzle discharge). The test series was designed to weed out any turbopump design and manufacturing flaws or fatigue issues prior to putting the C/SiC blisks into it. The second and third series were the C/SiC blisk test series. The primary goal of these series was to expose the blisks to as much fatigue causing dynamic stress as possible to examine the material's capability. Initially, the test plan was to put equal time on the two blisks, however, as the test series progressed, the funding allowed additional testing to occur. The additional test time was placed on the polar weave blisk. The total test time accrued on the polar blisk was 2550 seconds with 860 seconds near the turbopump design speed of 25,000 rpm. This testing included 6 tests / 775 seconds pumping liquid nitrogen and 7 tests / 1775 seconds pumping liquid oxygen. The drive gas for all of the tests was gaseous nitrogen due to the lack of hot gas source for the Simplex turbopump. The quasi-isotropic blisk was tested for XX total tests and XXXX seconds with X tests/XXXX seconds pumping liquid nitrogen and X tests/XXXX seconds pumping liquid \oxygen.During the test series, the blisks were inspected following each test. Inspections initially were viewed from the downstream side of the blisks only. Midway through the testing, a method of borescoping the leading edges of the blades was devised, and subsequently, both sides of the blades were inspected following each test. The leading and trailing edges of the polar blisk held up better than the quasi-isotropic blisk. This was a known possibility due to the varying fiber direction in the blades as the rectangular preform weave is cut in a circular pattern. The surprising fact about the testing was that there was no measurable performance loss due to the inaccuracies in the blade manufacturing in the C/SiC blisks, the surface roughness C/SiC of the blades, or the loss of the material in the polar blisk. A performance shift was seen in the quasi-isotropic blisk as portions of the leading and trailing edges were lost. After the testing was completed, detailed inspections of the blisks were performed. The largest surprise was the polar blisk had a obvious crack in a single blade that was located nearly midspan which was not detected in test. The crack ran completely through the blade circumferenciary and through the radial length of the blade. However, the crack does not appear to extend into the blisk hub. Although the cause of the crack is still under investigation, the material appears to be tolerant of this crack, and other hairline cracks discovered under higher magnification. This bodies well for eventual use of this material in actual flight turbopumps where monolithic fracture toughness issues limit its use.

Genge, Gary G.↗
NASA NTRS1995

An experimental study of light scattering by large, irregular particles

The intensity and polarization of light scattered by a variety of types of artificial partices large compared to the wavelength were measured as a function of phase angle. Shape, surface roughness, absorption coefficient, and internal scattering coefficient were varied systematically and their effects studied. Scattering by clear, smooth-surfaced spheres is in quantitative agreement with the predictions of the geometrical optics (ray theory) approximation to physical optics (Mie theory). The phase functions of almost all of the particles measured have both forward and backward scattering lobes. A two-parameter, double Henyey-Greenstein function generally provides reasonably good descriptions of the data, while keeping the number of free parameters to the minimum necessary. On a double Henyey- Greenstein parameter plot all of the particles fall into an L-shaped area of restricted size in which the location is characteristic of the particle type. Formalisms based on the equivalent slab model are also given for estimating the scattering efficiency of a large, irregular particle. For most dielectric particles the transmitted, forward scattered light is partially negatively polarized. It is this component that is respopnsible for the well-known maximum in the polarization curves of planetary regoliths at phase angles around 100 deg. For phase angles between about 30 deg and 70 deg the internally scattered light is found to be randomly polarized in the particles studied here, so that the only contribution to the second component of the Stokes vector is by Fresnel reflection from the particle surface. If this empirical result is general, measurement of the second Stokes vector of the light scattered from a regolith at these angles may provide a method of remotely measuring the mean refractive index.

Mcguire, Audrey F.↗
NASA NTRS2003

Effects of heat stress on baroreflex function in humans

INTRODUCTION: Heat stress significantly reduces orthostatic tolerance in humans. The mechanism(s) causing this response remain unknown. The purpose of this review article is to present data pertaining to the hypothesis that reduced orthostatic tolerance in heat stressed individuals is a result of heat stress induced alterations in baroflex function. METHODS: In both normothermic and heat stressed conditions baroreflex responsiveness was assessed via pharmacological and non-pharmacological methods. In addition, the effects of heat stress on post-synaptic vasoconstrictor responsiveness were assessed. RESULTS: Generally, whole body heating did not alter baroreflex sensitivity defined as the gain of the linear portion of the baroreflex curve around the operating point. However, whole body heating shifted the baroreflex curve to the prevailing (i.e. elevated) heart rate and muscle sympathetic nerve activity. Finally, the heat stress impaired vasoconstrictor responses to exogenous administration of adrenergic agonists. CONCLUSION: Current data do not support the hypothesis that reduced orthostatic tolerance associated with heat stress in humans is due to impaired baroreflex responsiveness. This phenomenon may be partially due to the effects of heat stress on reducing vasoconstrictor responsiveness.

Review, Tutorial↗
NASA NTRS1996

Independent Pixel and Two Dimensional Estimates of LANDSAT-Derived Cloud Field Albedo

A theoretical study has been conducted on the effects of cloud horizontal inhomogeneity on cloud albedo bias. A two-dimensional (2D) version of the Spherical Harmonic Discrete Ordinate Method (SHDOM) is used to estimate the albedo bias of the plane parallel (PP-IPA) and independent pixel (IPA-2D) approximations for a wide range of 2D cloud fields obtained from LANDSAT. They include single layer trade cumulus, open and closed cell broken stratocumulus, and solid stratocumulus boundary layer cloud fields over ocean. Findings are presented on a variety of averaging scales and are summarized as a function of cloud fraction, mean cloud optical depth, cloud aspect ratio, standard deviation of optical depth, and the gamma function parameter Y (a measure of the width of the optical depth distribution). Biases are found to be small for small cloud fraction or mean optical depth, where the cloud fields under study behave linearly. They are large (up to 0.20 for PP-IPA bias, -0.12 for IPA-2D bias) for large v. On a scene average basis PP-IPA bias can reach 0.30, while IPA-2D bias reaches its largest magnitude at -0.07. Biases due to horizontal transport (IPA-2D) are much smaller than PP-IPA biases but account for 20% RMS of the bias overall. Limitations of this work include the particular cloud field set used, assumptions of conservative scattering, constant cloud droplet size, no gas absorption or surface reflectance, and restriction to 2D radiative transport. The LANDSAT data used may also be affected by radiative smoothing.

Chambers, L. H.↗
NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗
NASA NTRS2007

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗
NASA NTRS2010

Direct Relationship Between Perceptual and Motor Variability

The time that elapses between stimulus onset and the onset of a saccadic eye movement is longer and more variable than can be explained by neural transmission times and synaptic delays (Carpenter, 1981, in: Eye Movements: Cognition & Visual Perception, Earlbaum). In theory, noise underlying response-time (RT) variability could arise at any point along the sensorimotor cascade, from sensory noise arising Vvithin the early visual processing shared Vvith perception to noise in the motor criterion or commands necessary to trigger movements. These two loci for internal noise can be distinguished empirically; sensory internal noise predicts that response time Vvill correlate Vvith perceived stimulus magnitude whereas motor internal noise predicts no such correlation. Methods. We used the data described by Liston and Stone (2008, JNS 28:13866-13875), in which subjects performed a 2AFC saccadic brightness discrimination task and the perceived brightness of the chosen stimulus was then quantified in a second 21FC perceptual task. Results. We binned each subject's data into quartiles for both signal strength (from dimmest to brightest) and RT (from slowest to fastest) and analyzed the trends in perceived brightness. We found significant effects of both signal strength (as expected) and RT on normalized perceived brightness (both p less than 0.0001, 2-way ANOVA), without significant interaction (p = 0.95, 2-way ANOVA). A plot of normalized perceived brightness versus normalized RT show's that more than half of the variance was shared (r2 = 0.56, P less than 0.0001). To rule out any possibility that some signal-strength related artifact was generating this effect, we ran a control analysis on pairs of trials with repeated presentations of identical stimuli and found that stimuli are perceived to be brighter on trials with faster saccades (p less than 0.001, paired t-test across subjects). Conclusion. These data show that shared early visual internal noise jitters perceived brightness and the saccadic motor output in parallel. While the present correlation could theoretically result, either directly or indirectly, from some low-level brainstem or retinal mechanism (e.g., arousal, pupil size, photoreceptor noise) that influences both visual and oculomotor circuits, this is unlikely given the earlier fin ding that the variability in perceived motion direction and smooth-pursuit motor output is highly correlated (Stone and Krauzlis, 2003, JOV 3:725-736), suggesting that cortical circuits contribute to the shared internal noise.

Liston, Dorion B.↗
NASA NTRS1994

Optimized design of a hypersonic nozzle

Conventional procedures for designing nozzles involve the design of an inviscid contour (using the method of characteristics) that is corrected with a displacement thickness calculated from boundary-layer theory. However, nozzles designed using this classical procedure have been shown to exhibit poor flow quality at Mach numbers characteristic of hypersonic applications. The nozzle to be designed will be a part of the NASA HYPULSE facility which is being used for hypervelocity flight research. Thus, the flow quality of the nozzle is a critical question that needs to be addressed. Design of nozzles for hypersonic applications requires a proper assessment of the effects of the thick boundary layer on the inviscid flowfield. Since the flow field is largely supersonic, the parabolized form of the Navier-Stokes (PNS) equations can be used. The requirement of a uniform flow at the exit plane of the nozzle can be used to define an objective function as part of an optimization procedure. The design procedure used in this study involves the coupling of a nonlinear (least-squares) optimization algorithm with an efficient, explicit PNS solver. The thick boundary layers growing on the walls of the nozzle limit the extent of the usable core region (region with uniform flow) for testing models (especially rectangular). In order to maximize the region of uniform flow, it was decided to have the exit plane of this nozzle to be (nearly) rectangular. Thus, an additional constraint on the nozzle shape resulted, namely the nozzle will have a shape transitioning from a circular one at the inlet to that of a rectangle at the exit. In order to provide for a smooth shape transition, the cross sectional contour of the nozzle is defined by a superellipse. The nozzle is taken to be a meter in length. The axial variations of the major and minor radii of the superellipse are governed by cubic splines. The design parameters are the coefficients of the splines associated with the local nozzle wall slopes. Extensive calculations have been made (with a three-dimensional Euler code) to understand the effects of various parameters such as location of the knot points of the spline function, different ways of characterizing the uniformity of the flow in the exit plane, as well as the effect of constraining the area of the nozzle to be invariant. Turbulent flow (measurements indicate that the flow at the nozzle inlet is turbulent) calculations are now being performed (with the inviscidly designed nozzle contours) to assess the flow quality.

Krishnamurthy, Ramesh↗
NASA NTRS1994

Automation of reverse engineering process in aircraft modeling and related optimization problems

During the year of 1994, the engineering problems in aircraft modeling were studied. The initial concern was to obtain a surface model with desirable geometric characteristics. Much of the effort during the first half of the year was to find an efficient way of solving a computationally difficult optimization model. Since the smoothing technique in the proposal 'Surface Modeling and Optimization Studies of Aerodynamic Configurations' requires solutions of a sequence of large-scale quadratic programming problems, it is important to design algorithms that can solve each quadratic program in a few interactions. This research led to three papers by Dr. W. Li, which were submitted to SIAM Journal on Optimization and Mathematical Programming. Two of these papers have been accepted for publication. Even though significant progress has been made during this phase of research and computation times was reduced from 30 min. to 2 min. for a sample problem, it was not good enough for on-line processing of digitized data points. After discussion with Dr. Robert E. Smith Jr., it was decided not to enforce shape constraints in order in order to simplify the model. As a consequence, P. Dierckx's nonparametric spline fitting approach was adopted, where one has only one control parameter for the fitting process - the error tolerance. At the same time the surface modeling software developed by Imageware was tested. Research indicated a substantially improved fitting of digitalized data points can be achieved if a proper parameterization of the spline surface is chosen. A winning strategy is to incorporate Dierckx's surface fitting with a natural parameterization for aircraft parts. The report consists of 4 chapters. Chapter 1 provides an overview of reverse engineering related to aircraft modeling and some preliminary findings of the effort in the second half of the year. Chapters 2-4 are the research results by Dr. W. Li on penalty functions and conjugate gradient methods for quadratic programming problems.

Li, W.↗