Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing technique”

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

Monte Carlo simulation of collisionless shocks showing preferential acceleration of high A/Z particles

A collisionless quasi-parallel shock is simulated by Monte Carlo techniques. The scattering of all velocity particles from thermal to high energy is assumed to occur so that the mean free path is directly proportional to velocity times the mass-to-charge-ratio, and inversely proporational to the plasma density. The shock profile and velocity spectra are obtained, showing preferential acceleration of high A/Z particles relative to protons. The inclusion of the back pressure of the scattering particles on the inflowing plasma produces a smoothing of the shock profile, which implies that the spectra are steeper than for a discontinuous shock.

Ellison, D. C.↗

Geosat altimeter observations of the surface circulation of the Southern Ocean

Using Geosat altimeter data for 26 months from November 1986 to December 1988 and a newly developed technique for the analysis of height data, the variability of the sea level and the surface geostrophic currents in the Southern Ocean is investigated. The processed Geosat data are used to examine the relationship between the mesoscale variability and the values of mean circulation, determined from historical hydrographic data. It is shown that the geographical patterns of both the mean flow and the mesoscale variability are correlated. An efficient objective-analysis algorithm for generating smoothed fields from observations randomly distributed in time and two space dimensions is developed and applied to 26 months of Geosat data. The smoothed fields are then used to investigate the large-scale low-frequency variability of the sea level and the surface geostrophic velocity in the Southern Ocean, in order to identify the mode of the observed variations.

Chelton, Dudley B.↗

Investigation of intrinsic variability in one-dimensional parallel shocks using steady state hybrid simulations

We have developed a means of producing a steady state hybrid simulation of a collisionless shock. The shock is stopped in the simulation box by transforming into the shock frame and by modifying the downstream boundary conditions to allow the plasma to flow through the simulation box. Once the shock is stationary in the box frame, the simulation can be run for an arbitrary time with a fixed box size and a fixed number of simulation particles. Using this technique, we have shown that certain gross properties associated with the shock, such as the particle distribution function (including energetic particles produced by Fermi acceleration) and the flow speed profile, are constant (except for statistical variations) over hundreds of gyroperiods when averaged over times short compared to the average residence time of energetic particles. Our results imply that any microphysical processes responsible for particle heating and/or injection into the Fermi mechanism can be viewed as smooth and continuous on timescales longer than a few gyroperiods.

Bennett, Lee↗

An efficient technique for determining apparent temperature distributions from antenna temperature measurements

A method by which the apparent microwave temperature characteristic of a flat scene is estimated from dual polarized measurements is derived and interpreted. Approximate linear relationships between antenna and apparent temperatures are established by weighting emission components in spherical bands under the assumption that the surface is isotropic. The weighting factors are formed by integrating the antenna pattern functions over these bands. The vector aspect of the formulation is retained to account for the difference between the definition of the antenna polarizations and the polarizations of the emitted fields. The method is largely applicable to the measurement of smooth temperature distributions by an antenna having good spatial resolution of the distributions and is considered efficient for inverting large volumes of measurements. Sample cases are presented and the implications of these cases on remote radiometer observations are discussed. It is shown that cross-coupling occurs between the polarizations of the emitted fields and the polarizations of the antenna. For this reason and because practical antennas have cross-polarized patterns associated with them, it is necessary to conduct measurements at both horizontal and vertical polarizations to realize the inversion. It is also made evident that thorough inversions require that the apparent temperatures be sampled at a sufficient number of points between nadir and zenith.

Claassen, J. P.↗

Incremental States for Precise On-Orbit Relative Knowledge in Formation Flight

High precision and close proximity formation flight is an enabling technology for future space missions and requires an on-board relative navigation capability that is accurate to the mm-level and robust to formation parameters. Common estimation techniques linearize the entire formation about one spacecraft’s position, resulting in degraded accuracy due to linearization errors when separation distances become large. Additionally, methods which decouple absolute and relative state estimates usually require ad-hoc methods to incorporate the estimates together. This work discusses an alternative ”incremental” formulation of the relative navigation problem which is in variant to formation size, robust to coupling between absolute and relative dynamics, and can undergo similarity transformations to smoothly incorporate either absolute or relative information without numerical issues. A specific example of this architecture is presented in the context of formation navigation using Carrier-Differential Global Positioning System measurements, and is compared to a traditional leader-linearized filter. In the presence of accurate measurements, the incremental architecture is shown to reduce linearization errors and mean estimate errors by up to three orders of magnitude without the incorporation of new sensor information.

Seubert, Carl↗

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitary time, G(bar a, bar t), and that of the time when a crack reaches an arbitary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. TAking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗

Creep life prediction based on stochastic model of microstructurally short crack growth

A nondimensional model of microstructurally short crack growth in creep is developed based on a detailed observation of the creep fracture process of 304 stainless steel. In order to deal with the scatter of small crack growth rate data caused by microstructural inhomogeneity, a random variable technique is used in the model. A cumulative probability of the crack length at an arbitrary time, G(bar a, bar t), and that of the time when a crack reaches an arbitrary length, F(bar t, bar a), are obtained numerically by means of a Monte Carlo method. G(bar a, bar t), and F(bar t, bar a) are the probabilities for a single crack. However, multiple cracks generally initiate on the surface of a smooth specimen from the early stage of creep life to the final stage. Taking into account the multiple crack initiations, the actual crack length distribution observed on the surface of a specimen is predicted by the combination of probabilities for a single crack. The prediction shows a fairly good agreement with the experimental result for creep of 304 stainless steel at 923 K. The probability of creep life is obtained from an assumption that creep fracture takes place when the longest crack reaches a critical length. The observed and predicted scatter of the life is fairly small for the specimens tested.

Kitamura, Takayuki↗

Approach for Structurally Clearing an Adaptive Compliant Trailing Edge Flap for Flight

The Adaptive Compliant Trailing Edge (ACTE) flap was flown on the National Aeronautics and Space Administration (NASA) Gulfstream GIII testbed at the NASA Armstrong Flight Research Center. This smoothly curving flap replaced the existing Fowler flaps creating a seamless control surface. This compliant structure, developed by FlexSys Inc. in partnership with the Air Force Research Laboratory, supported NASA objectives for airframe structural noise reduction, aerodynamic efficiency, and wing weight reduction through gust load alleviation. A thorough structures airworthiness approach was developed to move this project safely to flight. A combination of industry and NASA standard practice require various structural analyses, ground testing, and health monitoring techniques for showing an airworthy structure. This paper provides an overview of compliant structures design, the structural ground testing leading up to flight, and the flight envelope expansion and monitoring strategy. Flight data will be presented, and lessons learned along the way will be highlighted.

flight test↗

Image inversion analysis of the HST OTA (Hubble Space Telescope Optical Telescope Assembly), phase A

Technical work during September-December 1990 consisted of: (1) analyzing HST point source images obtained from JPL; (2) retrieving phase information from the images by a direct (noniterative) technique; and (3) characterizing the wavefront aberration due to the errors in the Hubble Space Telescope (HST) mirrors, in a preliminary manner. This work was in support of JPL design of compensating optics for the next generation wide-field planetary camera on HST. This digital technique for phase retrieval from pairs of defocused images, is based on the energy transport equation between these image planes. In addition, an end-to-end wave optics routine, based on the JPL Code 5 prescription of the unaberrated HST and WFPC, was derived for output of the reference phase front when mirror error is absent. Also, the Roddier routine unwrapped the retrieved phase by inserting the required jumps of +/- 2(pi) radians for the sake of smoothness. A least-squares fitting routine, insensitive to phase unwrapping, but nonlinear, was used to obtain estimates of the Zernike polynomial coefficients that describe the aberration. The phase results were close to, but higher than, the expected error in conic constant of the primary mirror suggested by the fossil evidence. The analysis of aberration contributed by the camera itself could be responsible for the small discrepancy, but was not verified by analysis.

Litvak, M. M.↗

Empirical convection models for northward IMF

It is clear that polar cap convection during times of northward Interplanetary Magnetic Field (IMF) is more structured and of lower mean speed than at times of southward IMF. This, coupled with the fact that the polar cap is smaller, means that empirical models are more difficult to construct with certainty. It is also clear that sunward flow deep in the polar cap is often observed, but its connection with the rest of the flow pattern is controversial. At present, empirical models are of three types: 'statistical' models wherein data from different days but with similar IMF conditions are averaged together; 'pattern recognition' models, which are built up by examining individually hundreds of passes to derive a 'typical' pattern which embodies features frequently observed; and 'assimilative' models, which use data of different types and from as many locations as possible, but all taken at the same time, in order to derive a snapshot (or series of snapshots) of the entire pattern. Each type of model has its own difficulties. Statistical models, by their very nature, smooth out flow features (e.g. the convection reversal, and the locus of sunward flow deep in the polar cap) which are not found at precisely the same invariant latitudes and magnetic local times on different days. Pattern recognition models are better at reproducing small-scale features, but the large-scale pattern can be a matter of interpretation. Assimilative models (such as AMIE) hold out the best hope for creating instantaneous, global convection patterns; however, the analysis technique tends to be most irregular (and least reliable) in the regions which are not well covered by in situ data. It appears that, at least at times, a four cell model with sunward flow at the highest and lowest latitudes, and antisunward flow in between, is consistent with the observations. At other times, the observations may be consistent with a two-cell convection pattern, but which includes significant meanders within the polar cap.

Moses, Julie J.↗

Turbulent flow between a rotating disk and a parallel wall

A scalar effective viscosity method is used to calculate the incompressible flow in the gap between two infinite, parallel smooth disks, one of them rotating and the other stationary. The flow is calculated for two values of the gap Reynolds number omega(s squared)/nu, namely 296 and 2852, where omega is the disk angular velocity, s the gap width, and nu the kinematic viscosity of the fluid; and over a range of radial Reynolds number omega(r squared)/nu from 0 to 10 million, where r is the radius from the axis of rotation. A numerical shooting technique is used to match the two boundary-layer flows arising on the disk and wall. Gradual transition is assumed for radial Reynolds numbers from 160,000 to 250,000. The results are in good agreement with Daily and Nece's measurements of velocity profiles and disk friction drag for enclosed rotating disks, and at sufficiently large radius the local skin friction approaches that of turbulent Couette flow.

Cooper, P.↗

The reduction correction in North America

An inverse Poisson integral technique was used to determine a gravity field on the geoid which, when continued by analytic free space methods to the topographic surface, agrees with the observed field. The computation is performed in three stages, each stage refining the previous solution using data at progressively increasing resolution (1 x 1 deg, 5 x 5', 5/8 x 5/8') from a decreasing area of integration. Reduction corrections are computed at 5/8 x 5/8' granularity by differencing the geoidal and surface values, smoothed by low-pass filtering and sub-sampled at 5' intervals. The 1 x 1 deg averages of the reduction corrections thus obtained for 172 1 x 1 deg squares in western North America are discussed. The 1 x 1 deg mean reduction corrections are predominantly positive, varying from -3 to +15 mgal, with values in excess of 5 mgal for 26 squares. Their mean and rms values are +2.4 and 3.6 mgal respectively and they correlate well with the mean terrain corrections. The mean and rms contributions from the three stages of computation are: 1 x 1 deg stage +0.15 and 0.7 mgal; 5 x 5' stage + 1.0 and 1.6 mgal; and 5/8 x 5/8' stage +1.3 and 1.8 mgal. These results reflect a tendency for the contributions to become larger and more systematically positive as the wavelengths involved become shorter. The results are discussed in terms of two mechanisms; the first is a tendency for the absolute values of both positive and negative anomalies to become larger when continued downwards and, the second, a non-linear rectification, due to the correlation between gravity anomaly and topographic height, which results in the values continued to a level surface being systematically more positive than those on the topography.

Martzen, P. D.↗

Mercury Crustal and Lithospheric Properties from High-resolution Gravity Field Models

We have analyzed all the available tracking data from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission to determine high-resolution gravity field models. We have applied a technique based on line-of-sight accelerations that results in models with improved correlations with topography in areas where the spacecraft altitude was low. We have used these models in an analysis of admittance to determine properties of the crust and lithosphere. We focus on four areas where correlations between gravity and topography are high. These locations represent different areas on Mercury: the high-Mg region, the Strindberg crater plus some lobate scarps, heavily cratered terrain, and smooth plains (the northern rise). We find crustal densities around 2600 kg m-3 for the first three areas, and a higher density for the northern rise, indicative of volcanic materials. Elastic thickness is generally low and varies between 11 and 30 km.

Sander Johannes Goossens↗

A catalog of jovian decameter-wave radio observations from 1957 - 1978

Data from over 200,000 hours of observation of Jupiter radio emission in the decameter-wavelength band, were collected from 13 observing sites and are available on magnetic tape. Observations were made at 14 fixed frequencies from 5 to 30 MHz. The characteristics of the tape recording technique and the data format are described. The combination of overlapping data from observing sites scattered world-wide lessens the effect of the earth's daily interruption of the ground-received signal. A power spectral analysis of the data shows no evidence of periodicities within the data other than the well-known influences of Jupiter, Io, and the earth. The dependence of the occurrence probability of emission on System 3 longitude and the phase of Io varies smoothly with frequency down to 15 MHz and then appears quite different at 10 MHz. The morphology of the radio sources is both complex and stable for periods of at least months and probably much longer.

Thieman, J. R.↗

A multigrid method for steady Euler equations on unstructured adaptive grids

A flux-difference splitting type algorithm is formulated for the steady Euler equations on unstructured grids. The polynomial flux-difference splitting technique is used. A vertex-centered finite volume method is employed on a triangular mesh. The multigrid method is in defect-correction form. A relaxation procedure with a first order accurate inner iteration and a second-order correction performed only on the finest grid, is used. A multi-stage Jacobi relaxation method is employed as a smoother. Since the grid is unstructured a Jacobi type is chosen. The multi-staging is necessary to provide sufficient smoothing properties. The domain is discretized using a Delaunay triangular mesh generator. Three grids with more or less uniform distribution of nodes but with different resolution are generated by successive refinement of the coarsest grid. Nodes of coarser grids appear in the finer grids. The multigrid method is started on these grids. As soon as the residual drops below a threshold value, an adaptive refinement is started. The solution on the adaptively refined grid is accelerated by a multigrid procedure. The coarser multigrid grids are generated by successive coarsening through point removement. The adaption cycle is repeated a few times. Results are given for the transonic flow over a NACA-0012 airfoil.

Riemslagh, Kris↗

Ditching Investigation of a 1/24-Scale Model of the Boeing B-47 Airplane

An investigation of a 1/24-scale dynamically similar model of the Boeing B-47 airplane was made to determine the ditching characteristics and proper ditching technique for the airplane. Various conditions of damage, landing attitude, flap setting, and speed were investigated. The behavior of the model was determined from visual observations, motion-picture records, and time-history deceleration records. The results of the investigation are presented in table form, photographs, and curves. The airplane should be ditched at the lowest speed and highest attitude consistent with adequate control; the flaps should be full down. The airplane will probably make a deep but fairly smooth run. The fuselage bottom will be damaged and partially filled with water; consequently, crew members should be assigned ditching stations near an exit in the upper or forward part of the fuselage. The nacelles may be expected to be torn away from the wing. In calm water the maximum decelerations will be about 3g and the landing run will be about 6 fuselage lengths.

Fisher, Lloyd J.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗