Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “successive linear approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Mass flow meter using the triboelectric effect for measurement in cryogenics

The use of triboelectric charge to measure the mass flow rate of cryogens for the Space Shuttle Main Engine was investigated. Cross correlation of the triboelectric charge signals was used to determine the transit time of the cryogen between two sensor locations in a .75-in tube. The ring electrode sensors were mounted in a removable spool piece. Three spool pieces were constructed for delivery, each with a different design. One set of electronics for implementation of the cross correlation and flow calculation was constructed for delivery. Tests were made using a laboratory flow loop using liquid freon and transformer oil. The measured flow precision was 1 percent and the response was linear. The natural frequency distribution of the triboelectric signal was approximately 1/f. The sensor electrodes should have an axial length less than approximately one/tenth pipe diameter. The electrode spacing should be less than approximately one pipe diameter. Tests using liquid nitrogen demonstrated poor tribo-signal to noise ratio. Most of the noise was microphonic and common to both electrode systems. The common noise rejection facility of the correlator was successful in compensating for this noise but the signal was too small to enable reliable demonstration of the technique in liquid nitrogen.

Bernatowicz, Henry↗

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗

A comparative look at sunspot cycles

On the basis of cycles 8 through 20, spanning about 143 years, observations of sunspot number, smoothed sunspot number, and their temporal properties were used to compute means, standard deviations, ranges, and frequency of occurrence histograms for a number of sunspot cycle parameters. The resultant schematic sunspot cycle was contrasted with the mean sunspot cycle, obtained by averaging smoothed sunspot number as a function of time, tying all cycles (8 through 20) to their minimum occurence date. A relatively good approximation of the time variation of smoothed sunspot number for a given cycle is possible if sunspot cycles are regarded in terms of being either HIGH- or LOW-R(MAX) cycles or LONG- or SHORT-PERIOD cycles, especially the latter. Linear regression analyses were performed comparing late cycle parameters with early cycle parameters and solar cycle number. The early occurring cycle parameters can be used to estimate later occurring cycle parameters with relatively good success, based on cycle 21 as an example. The sunspot cycle record clearly shows that the trend for both R(MIN) and R(MAX) was toward decreasing value between cycles 8 through 14 and toward increasing value between cycles 14 through 20. Linear regression equations were also obtained for several measures of solar activity.

Wilson, R. M.↗

An Update on the Effectiveness of Postmission Disposal in LEO

The commonly‐adopted orbital debris mitigation measures were developed to reduce the growth of the future debris population. A major component in debris mitigation is post-mission disposal (PMD). The key PMD element for LEO satellites is the 25‐year rule. It is intended to limit the long‐term presence of rocket bodies (R/Bs) and spacecraft (S/C), as well as mission-related debris, in the environment. The effectiveness of PMD has been demonstrated and documented since the development of mitigation measures began in the 1990s. This paper summarizes an updated study, based on the current environment, using the NASA LEGEND model. The study focused on the > or = 10 cm population in LEO. The historical simulation covered 1957 through 2011 and followed the recorded launches and known breakup events. The future projection was carried out for 200 years. An eight‐year launch traffic, 2004 - 2011, was repeated during the projection period. An eight‐year mission lifetime was assumed for future S/C. No stationkeeping and no collision avoidance maneuver were implemented. Only objects 10 cm and larger were included in collision consideration. No explosion was allowed for R/Bs and S/C launched after 2011. The 25‐year PMD rule success rates were set at 0%, 10%, 50%, 75%, and 95%, respectively, for the 5 study scenarios. Results of the simulations were analyzed to quantify the differences among the different compliance rates. As expected, the 0% PMD projection followed a rapid and non‐linear increase in the next 200 years. The LEO population, on average, more than tripled at the end of the simulations. With a 50% compliance of the 25‐year rule, the population growth was reduced approximately by half. However, even with a 95% compliance of the 25‐year rule, the LEO debris population would still increase by an average of more than 50% in 200 years. These simulation results provide an updated assessment of the effectiveness of the 25‐year rule. It is the first and the most cost‐effective defense against future population growth. In addition, the results also confirm the instability of the LEO population and lay the foundation for the need to consider environment remediation in the future.

Liou, J.-C.↗

Geometrical acoustics and transonic helicopter sound

A new method is presented for predicting the impulsive noise generated by a transonic rotor blade. The method is a combined approach involving computational fluid dynamics and geometrical acoustics. A full-potential finite-difference method is used to obtain the pressure field close to the blade. A Kirchhoff integral formulation is then used to extend these finite-difference results into the far field. This Kirchhoff formula is based on geometrical acoustics approximations. It requires initial data across a plane at the sonic radius in a blade-fixed coordinate system. This data is provided by the finite-difference solution. Acoustic pressure predictions show good agreement with hover experimental data for cases with hover tip Mach numbers of 0.88 through 0.96. The cases above 0.92 tip Mach number are dominated by non-linear transonic effects seen as strong shocks on and off the blade tip. This paper gives the first successful predictions of far-field acoustic pressures for high-speed impulsive noise over a range of Mach numbers after delocalization.

Isom, Morris↗

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗

High mass resolution, high angular acceptance time-of-flight mass spectroscopy for planetary missions under the Planetary Instrument Definition and Development Program (PIDDP)

This final report covers three years and several phases of work in which instrumentation for the Planetary Instrument Definition and Development Program (PIDDP) were successfully developed. There were two main thrusts to this research: (1) to develop and test methods for electrostatically scanning detector field-of-views, and (2) to improve the mass resolution of plasma mass spectrometers to M/delta M approximately 25, their field-of-view (FOV) to 360 degrees, and their E-range to cover approximately 1 eV to 50 keV. Prototypes of two different approaches to electrostatic scanning were built and tested. The Isochronous time-of-flight (TOF) and the linear electric field 3D TOF devices were examined.

Young, David T.↗

LSA Large Area Silicon Sheet Task Continuous Czochralski Process Development

A commercial Czochralski crystal growing furnace was converted to a continuous growth facility by installation of a small, in-situ premelter with attendant silicon storage and transport mechanisms. Using a vertical, cylindrical graphite heater containing a small fused quartz test tube linear from which the molten silicon flowed out the bottom, approximately 83 cm of nominal 5 cm diamter crystal was grown with continuous melt addition furnished by the test tube premelter. High perfection crystal was not obtained, however, due primarily to particulate contamination of the melt. A major contributor to the particulate problem was severe silicon oxide buildup on the premelter which would ultimately drop into the primary melt. Elimination of this oxide buildup will require extensive study and experimentation and the ultimate success of continuous Czochralski depends on a successful solution to this problem. Economically, the continuous Czochralski meets near-term cost goals for silicon sheet material.

Rea, S. N.↗

Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids II: Extension to Two Dimensional Scalar Equation

The framework for constructing a high-order, conservative Spectral (Finite) Volume (SV) method is presented for two-dimensional scalar hyperbolic conservation laws on unstructured triangular grids. Each triangular grid cell forms a spectral volume (SV), and the SV is further subdivided into polygonal control volumes (CVs) to supported high-order data reconstructions. Cell-averaged solutions from these CVs are used to reconstruct a high order polynomial approximation in the SV. Each CV is then updated independently with a Godunov-type finite volume method and a high-order Runge-Kutta time integration scheme. A universal reconstruction is obtained by partitioning all SVs in a geometrically similar manner. The convergence of the SV method is shown to depend on how a SV is partitioned. A criterion based on the Lebesgue constant has been developed and used successfully to determine the quality of various partitions. Symmetric, stable, and convergent linear, quadratic, and cubic SVs have been obtained, and many different types of partitions have been evaluated. The SV method is tested for both linear and non-linear model problems with and without discontinuities.

Wang, Z. J.↗

Additional Results on the Radiating Properties of Ltsas and Cwsas

Results are presented on the theoretical experimental investigations of the radiating properties of Linearly Tapered Slot Antennas (LTSAs) and Constant Width Slot Antennas (CWSAs). The antenna is treated as a flared slot radiating in the presence of a conducting half plane. The aperture distribution in the flared slot is determined by approximating the tapered structure as a series of short sections of parallel slots. The slot field of a uniform slotline is determined by employing the spectral Galerkin's technique and a power conservation criterion is employed to match the fields at the step junction between two successive sections. Comparison is made between theory and experiment for several LTSAs and CWSAs etched on epsilon sub r = 2.22 and 2.55 substrates, and or flare angles 4 deg or = gamma or = 10 deg. It is shown that the theory adequately models the physics of the problem when the lateral dimension D of the antenna is large and when the propagation constant of a parallel slot line can be estimated very accurately. Newly observed effects of the lateral truncation on the radiation pattern of the antenna are discussed.

Source record↗

Non-Linear Cosmological Power Spectra in Real and Redshift Space

We present an expression for the non-linear evolution of the cosmological power spectrum based on Lagrangian trajectories. This is simplified using the Zel'dovich approximation to trace particle displacements, assuming Gaussian initial conditions. The model is found to exhibit the transfer of power from large to small scales expected in self-gravitating fields. Some exact solutions are found for power-law initial spectra. We have extended this analysis into red-shift space and found a solution for the non-linear, anisotropic redshift-space power spectrum in the limit of plane-parallel redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of power-law initial spectra. We find that the shape of this ratio depends on the shape of the initial spectrum, but when scaled to linear theory depends only weakly on the redshift-space distortion parameter, beta. The point of zero-crossing of the quadrupole, kappa(sub o), is found to obey a simple scaling relation and we calculate this scale in the Zel'dovich approximation. This model is found to be in good agreement with a series of N-body simulations on scales down to the zero-crossing of the quadrupole, although the wavenumber at zero-crossing is underestimated. These results are applied to the quadrupole-to-monopole ratio found in the merged QDOT plus 1.2-Jy-IRAS redshift survey. Using a likelihood technique we have estimated that the distortion parameter is constrained to be beta greater than 0.5 at the 95 percent level. Our results are fairly insensitive to the local primordial spectral slope, but the likelihood analysis suggests n = -2 un the translinear regime. The zero-crossing scale of the quadrupole is k(sub 0) = 0.5 +/- 0.1 h Mpc(exp -1) and from this we infer that the amplitude of clustering is sigma(sub 8) = 0.7 +/- 0.05. We suggest that the success of this model is due to non-linear redshift-space effects arising from infall on to caustic and is not dominated by virialized cluster cores. The latter should start to dominate on scales below the zero-crossing of the quadrupole, where our model breaks down.

Taylor, A. N.↗

Multiscale Morphological Filtering for Analysis of Noisy and Complex Images

Images acquired with passive sensing techniques suffer from illumination variations and poor local contrasts that create major difficulties in interpretation and identification tasks. On the other hand, images acquired with active sensing techniques based on monochromatic illumination are degraded with speckle noise. Mathematical morphology offers elegant techniques to handle a wide range of image degradation problems. Unlike linear filters, morphological filters do not blur the edges and hence maintain higher image resolution. Their rich mathematical framework facilitates the design and analysis of these filters as well as their hardware implementation. Morphological filters are easier to implement and are more cost effective and efficient than several conventional linear filters. Morphological filters to remove speckle noise while maintaining high resolution and preserving thin image regions that are particularly vulnerable to speckle noise were developed and applied to SAR imagery. These filters used combination of linear (one-dimensional) structuring elements in different (typically four) orientations. Although this approach preserves more details than the simple morphological filters using two-dimensional structuring elements, the limited orientations of one-dimensional elements approximate the fine details of the region boundaries. A more robust filter designed recently overcomes the limitation of the fixed orientations. This filter uses a combination of concave and convex structuring elements. Morphological operators are also useful in extracting features from visible and infrared imagery. A multiresolution image pyramid obtained with successive filtering and a subsampling process aids in the removal of the illumination variations and enhances local contrasts. A morphology-based interpolation scheme was also introduced to reduce intensity discontinuities created in any morphological filtering task. The generality of morphological filtering techniques in extracting information from a wide variety of images obtained with active and passive sensing techniques is discussed. Such techniques are particularly useful in obtaining more information from fusion of complex images by different sensors such as SAR, visible, and infrared.

Kher, A.↗

Characterization of the Nimbus-7 SBUV radiometer for the long-term monitoring of stratospheric ozone

Precise knowledge of in-orbit sensitivity change is critical for the successful monitoring of stratospheric ozone by satellite-based remote sensors. This paper evaluates those aspects of the in-flight operation that influence the long-term stability of the upper stratospheric ozone measurements made by the Nimbus-7 SBUV spectroradiometer and chronicles methods used to maintain the long-term albedo calibration of this UV sensor. It is shown that the instrument's calibration for the ozone measurement, the albedo calibration, has been maintained over the first 6 yr of operation to an accuracy of approximately + or - 2 percent. The instrument's wavelength calibration is shown to drift linearly with time. The knowledge of the SBUV wavelength assignment is maintained to a 0.02-nm precision.

Cebula, Richard P.↗

Parallel Preconditioning for CFD Problems on the CM-5

Up to today, preconditioning methods on massively parallel systems have faced a major difficulty. The most successful preconditioning methods in terms of accelerating the convergence of the iterative solver such as incomplete LU factorizations are notoriously difficult to implement on parallel machines for two reasons: (1) the actual computation of the preconditioner is not very floating-point intensive, but requires a large amount of unstructured communication, and (2) the application of the preconditioning matrix in the iteration phase (i.e. triangular solves) are difficult to parallelize because of the recursive nature of the computation. Here we present a new approach to preconditioning for very large, sparse, unsymmetric, linear systems, which avoids both difficulties. We explicitly compute an approximate inverse to our original matrix. This new preconditioning matrix can be applied most efficiently for iterative methods on massively parallel machines, since the preconditioning phase involves only a matrix-vector multiplication, with possibly a dense matrix. Furthermore the actual computation of the preconditioning matrix has natural parallelism. For a problem of size n, the preconditioning matrix can be computed by solving n independent small least squares problems. The algorithm and its implementation on the Connection Machine CM-5 are discussed in detail and supported by extensive timings obtained from real problem data.

Simon, Horst D.↗

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗

Study of modal coupling procedures for the shuttle: A matrix method for damping synthesis

The damping method was applied successfully to real structures as well as analytical models. It depends on the ability to determine an appropriate modal damping matrix for each substructure. In the past, modal damping matrices were assumed diagonal for lack of being able to determine the coupling terms which are significant in the general case of nonproportional damping. This problem was overcome by formulating the damped equations of motion as a linear perturbation of the undamped equations for light structural damping. Damped modes are defined as complex vectors derived from the complex frequency response vectors of each substructure and are obtained directly from sinusoidal vibration tests. The damped modes are used to compute first order approximations to the modal damping matrices. The perturbation approach avoids ever having to solve a complex eigenvalue problem.

Hasselman, T. K.↗

RAISHIN: A High-Resolution Three-Dimensional General Relativistic Magnetohydrodynamics Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.

Mizuno, Yosuke↗

Sensor Calibration Impacts on Dust Detection Based On MODIS And VIIRS Thermal Emissive Bands

Dust detection using remotely sensed measurements has been one of the challenging problems encountered by atmospheric scientists. MODerate Resolution Imaging Spectroradiometer (MODIS) on the Terra (T)and Aqua (A) platforms have been a versatile sensor for well over 21 and 18 years respectively and have been extremely useful in the retrieval of aerosol information over the entire globe. The MODIS radiances from the Level1B in general are expected to be within 5% accuracy in the reflective wavelengths and within 1% in the thermal emissive wavelengths. In this paper, we evaluate the sensitivity of previously developed dust detection technique based on thermal emissive wavelengths, which correspond to MODIS bands 20, 29, 31, and 32 respectively. The Thermal Emissive Dust Index (TEDI) performed very comparably to the traditional Aerosol Optical Thickness (AOT) retrievals by MODIS reflective channels. Since the MODIS Thermal Emissive Bands (TEB) are well calibrated on-orbit using a BlackBody (BB) source, the calibration of these long wave infrared bands is quite robust. As A-MODIS continues to perform well beyond its designed lifetime of 6 years, the instrument has undergone various levels of degradation during its mission time. As a consequence, it is imperative to check the impacts of calibration on the higher-level retrievals. In this paper, we rigorously analyze the sensitivity of TEDI due to the impact of calibration by the afore-mentioned TEB. The perturbation of the dominant (linear) calibration term demonstrated the following: first, there was a correlation in the sensitivity of the TEDI due to the uncertainty in the linear calibration term. Based on a perturbation in the linear calibration term for all aforementioned bands over a range of ± 5 % yielded the TEDI sensitivity to vary from approximately -3.2 % to about -3.6 %. When considering the uncertainty in each individual band significant changes were observed. The least change was observed for the perturbation in the calibration of band 20 with the TEDI sensitivity and the largest sensitivity in TEDI was observed in the perturbation of band 31 calibration. Thus, in the case of TEDI, noticeable sensitivity due to calibration uncertainty was observed in bands 29, 31, and 32, reiterating the importance of the TEB calibration in these bands. Also, the dust detection scheme based on A-MODIS was successfully transferred to the follow-on sensors such as Suomi (SNPP) and NOAA 20 (N20) VIIRS. The results presented in this paper would be extremely helpful in understanding impacts of calibration on the higher-level products for both current and future missions based on the MODIS heritage. Finally, the work also identifies the importance of radiometric fidelity in maintaining the accuracy of the dust detection. Results presented will show drastic improvement of the Saharan dust detection after the reduction of the electronic crosstalk in the 8.5 μm channel of T-MODIS.

Sriharsha Madhavan↗