SEARCH · Engineering Papers

Results for “error propagation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

NASA NTRS2001

The Linear Bicharacteristic Scheme for Electromagnetics

The upwind leapfrog or Linear Bicharacteristic Scheme (LBS) has previously been implemented and demonstrated on electromagnetic wave propagation problems. This paper extends the Linear Bicharacteristic Scheme for computational electromagnetics to model lossy dielectric and magnetic materials and perfect electrical conductors. This is accomplished by proper implementation of the LBS for homogeneous lossy dielectric and magnetic media and for perfect electrical conductors. Heterogeneous media are modeled through implementation of surface boundary conditions and no special extrapolations or interpolations at dielectric material boundaries are required. Results are presented for one-dimensional model problems on both uniform and nonuniform grids, and the FDTD algorithm is chosen as a convenient reference algorithm for comparison. The results demonstrate that the explicit LBS is a dissipation-free, second-order accurate algorithm which uses a smaller stencil than the FDTD algorithm, yet it has approximately one-third the phase velocity error. The LBS is also more accurate on nonuniform grids.

Beggs, John H.↗

NASA NTRS2011

Orbit Determination of Spacecraft in Earth-Moon L1 and L2 Libration Point Orbits

The ARTEMIS mission, part of the THEMIS extended mission, is the first to fly spacecraft in the Earth-Moon Lissajous regions. In 2009, two of the five THEMIS spacecraft were redeployed from Earth-centered orbits to arrive in Earth-Moon Lissajous orbits in late 2010. Starting in August 2010, the ARTEMIS P1 spacecraft executed numerous stationkeeping maneuvers, initially maintaining a lunar L2 Lissajous orbit before transitioning into a lunar L1 orbit. The ARTEMIS P2 spacecraft entered a L1 Lissajous orbit in October 2010. In April 2011, both ARTEMIS spacecraft will suspend Lissajous stationkeeping and will be maneuvered into lunar orbits. The success of the ARTEMIS mission has allowed the science team to gather unprecedented magnetospheric measurements in the lunar Lissajous regions. In order to effectively perform lunar Lissajous stationkeeping maneuvers, the ARTEMIS operations team has provided orbit determination solutions with typical accuracies on the order of 0.1 km in position and 0.1 cm/s in velocity. The ARTEMIS team utilizes the Goddard Trajectory Determination System (GTDS), using a batch least squares method, to process range and Doppler tracking measurements from the NASA Deep Space Network (DSN), Berkeley Ground Station (BGS), Merritt Island (MILA) station, and United Space Network (USN). The team has also investigated processing of the same tracking data measurements using the Orbit Determination Tool Kit (ODTK) software, which uses an extended Kalman filter and recursive smoother to estimate the orbit. The orbit determination results from each of these methods will be presented and we will discuss the advantages and disadvantages associated with using each method in the lunar Lissajous regions. Orbit determination accuracy is dependent on both the quality and quantity of tracking measurements, fidelity of the orbit force models, and the estimation techniques used. Prior to Lissajous operations, the team determined the appropriate quantity of tracking measurements that would be needed to meet the required orbit determination accuracies. Analysts used the Orbit Determination Error Analysis System (ODEAS) to perform covariance analyses using various tracking data schedules. From this analysis, it was determined that 3.5 hours of DSN TRK-2-34 range and Doppler tracking data every other day would suffice to meet the predictive orbit knowledge accuracies in the Lissajous region. The results of this analysis are presented. Both GTDS and ODTK have high-fidelity environmental orbit force models that allow for very accurate orbit estimation in the lunar Lissajous regime. These models include solar radiation pressure, Earth and Moon gravity models, third body gravitational effects from the Sun, and to a lesser extent third body gravitational effects from Jupiter, Venus, Saturn, and Mars. Increased position and velocity uncertainties following each maneuver, due to small execution performance errors, requires that several days of post-maneuver tracking data be processed to converge on an accurate post-maneuver orbit solution. The effects of maneuvers on orbit determination accuracy will be presented, including a comparison of the batch least squares technique to the extended Kalman filter/smoother technique. We will present the maneuver calibration results derived from processing post-maneuver tracking data. A dominant error in the orbit estimation process is the uncertainty in solar radiation pressure and the resultant force on the spacecraft. An estimation of this value can include many related factors, such as the uncertainty in spacecraft reflectivity and surface area which is a function of spacecraft orientation (spin-axis attitude), uncertainty in spacecraft wet mass, and potential seasonal variability due to the changing direction of the Sun line relative to the Earth-Moon Lissajous reference frame. In addition, each spacecraft occasionally enters into Earth or Moon penumbra or umbra and these shadow crossings reduche solar radiation force for several hours. The effects of these events on orbit determination accuracy will be presented. In order to plan for upcoming stationkeeping maneuvers, the maneuver planning team must take the current orbit estimate, propagate it forward to the planned maneuver time, and determine the optimal maneuver to maintain the Lissajous orbit for one or more revolutions. The propagation is performed using a Runge-Kutta 7/8 integrator and typically the position and velocity uncertainty increases with propagation time, increasing the overall uncertainty of the orbit state at the maneuver execution time. The effect of orbit knowledge uncertainty on stationkeeping operations will be presented.

Woodard, Mark↗

NASA NTRS1975

Microwave power transmission system studies. Volume 2: Introduction, organization, environmental and spaceborne systems analyses

Introduction, organization, analyses, conclusions, and recommendations for each of the spaceborne subsystems are presented. Environmental effects - propagation analyses are presented with appendices covering radio wave diffraction by random ionospheric irregularities, self-focusing plasma instabilities and ohmic heating of the D-region. Analyses of dc to rf conversion subsystems and system considerations for both the amplitron and the klystron are included with appendices for the klystron covering cavity circuit calculations, output power of the solenoid-focused klystron, thermal control system, and confined flow focusing of a relativistic beam. The photovoltaic power source characteristics are discussed as they apply to interfacing with the power distribution flow paths, magnetic field interaction, dc to rf converter protection, power distribution including estimates for the power budget, weights, and costs. Analyses for the transmitting antenna consider the aperture illumination and size, with associated efficiencies and ground power distributions. Analyses of subarray types and dimensions, attitude error, flatness, phase error, subarray layout, frequency tolerance, attenuation, waveguide dimensional tolerances, mechanical including thermal considerations are included. Implications associated with transportation, assembly and packaging, attitude control and alignment are discussed. The phase front control subsystem, including both ground based pilot signal driven adaptive and ground command approaches with their associated phase errors, are analyzed.

Maynard, O. E.↗

NASA NTRS2000

The Linear Bicharacteristic Scheme for Computational Electromagnetics

The upwind leapfrog or Linear Bicharacteristic Scheme (LBS) has previously been implemented and demonstrated on electromagnetic wave propagation problems. This paper extends the Linear Bicharacteristic Scheme for computational electromagnetics to treat lossy dielectric and magnetic materials and perfect electrical conductors. This is accomplished by proper implementation of the LBS for homogeneous lossy dielectric and magnetic media, and treatment of perfect electrical conductors (PECs) are shown to follow directly in the limit of high conductivity. Heterogeneous media are treated through implementation of surface boundary conditions and no special extrapolations or interpolations at dielectric material boundaries are required. Results are presented for one-dimensional model problems on both uniform and nonuniform grids, and the FDTD algorithm is chosen as a convenient reference algorithm for comparison. The results demonstrate that the explicit LBS is a dissipation-free, second-order accurate algorithm which uses a smaller stencil than the FDTD algorithm, yet it has approximately one-third the phase velocity error. The LBS is also more accurate on nonuniform grids.

Beggs, John H.↗

NASA NTRS1981

Estimation of the propagation delay through the troposphere from microwave radiometer data

Microwave propagation delay through the troposphere is studied as a means of estimating the adverse effect of atmospheric water vapor on the accuracy of measurements made with very long baseline interferometry. Vertical profiles of temperature and water vapor density were obtained from a total of 240 radiosonde launches taken simultaneously at three New England locations in 1974; all studies were made at the 19 and 22.2 GHz frequency operating range of the radiometers. The wet path length in the zenith direction could be estimated from the surface water vapor density to an accuracy of 5 cm for the summer data and 2 cm for winter data, and with the brightness temperatures, it could be estimated to an accuracy of 0.3 cm. Two dual-frequency radiometers were also used to determine the accuracy of prediction of the path length from real radiometry data. The rms error for the estimation of wet path length from surface meteorological parameters was 3.2 cm, and for the radiometer brightness temperatures the rms error was 1.5 cm.

Moran, J. M.↗

NASA NTRS1995

Radial alignment simulation of solar wind streams observed by Pioneers 10 and 11 in 1974

A particularly favorable radial lineup between spacecraft in Earth orbit and Pioneers 10 and 11 (near the ecliptic plane at approximately 5.5 AU and approximately 4.5 AU, respectively) occurred in mid-1974, when the solar corona was in a stable and well-defined warped-dipole configuration. The radial alignment study reported here differs from previous applications of the technique in two respects: (1) It is the first time a two-dimensional (2-D) MHD model has been tested over such a lengthy propagation interval; the 2-D capability is crucial for treatment of the nonradial shearing motions occurring across the stream interface. (2) The three-dimensional (3-D) structure observed in the white light corona at that time is related to the systematic patterns of nonradial flow deflections appearing at the Pioneer corotating interaction region (CIR) fronts. Comparison of predicted and observed flows for pairs of streams in two successive rotations reveals that when the parent coronal hole projects far across the solar equator in a predominantly north-south orientation (i.e., nearest the 2-D idealization), the mapping is accurate down to details of the flow structures. But where the spacecraft tracks along a latitudinal boundary of a hole or the associated stream front is inclined at a shallow angle to the equator, the numerical projections deviate systematically from the observations. Among the sources of error are 3-D dynamical interactions neglected in the model, differential rotation effects, and slow temporal evolution of the coronal structures. A characteristic pattern of north-south and east-west deflections indicative of the 3-D geometry appears across the CIRs, but these patterns evidently reflect primarily the local, as opposed to global, orientation of the 3-D stream fronts. Such patterns appear common in CIRs observed by Pioneer during this period. These findings thus hold special relevance for the analysis of Ulysses observations, since the present coronal configuration is similar to that of 1974.

Pizzo, V. J.↗

NASA NTRS2010

Large Telescope Segmented Primary Mirror Alignment

A document discusses a broadband (white light) point source, located at the telescope Cassegrain focus, which generates a cone of light limited by the hole in the secondary mirror (SM). It propagates to the aspheric null-mirror, which is optimized to make all the reflected rays to be normal to the primary mirror (PM) upon reflection. PM retro-reflects the rays back through the system for wavefront analysis. The point source and the wavefront analysis subsystems are all located behind the PM. The PM phasing is absolute (white light) and does not involve the SM. A relatively small, aspheric null-mirror located near the PM center of curvature has been designed to deliver the high level of optical wavefront correction. The phasing of the segments is absolute due to the use of a broadband source. The segmented PM is optically aligned independently and separately from the SM alignment. The separation of the PM segments alignment from the PM to the SM, and other telescope optics alignments, may be a significant advantage, eliminating the errors coupling. The point source of this concept is fully cooperative, unlike a star or laser-generated guide-star, providing the necessary brightness for the optimal S/N ratio, the spectral content, and the stable on-axis position. This concept can be implemented in the lab for the PM initial alignment, or made to be a permanent feature of the space-based or groundbased telescope.

Rud, Mayer↗

NASA NTRS1999

Composite Fuselage Technology

Work was conducted over a ten-year period to address the central issue of damage in primary load-bearing aircraft composite structure, specifically fuselage structure. This included the three facets of damage resistance, damage tolerance, and damage arrest. Experimental, analytical, and numerical work was conducted in order to identify and better understand the mechanisms that control the structural behavior of fuselage structures in their response to the three aspects of damage. Furthermore, work was done to develop straightforward design methodologies that can be employed by structural designers in preliminary design stages to make intelligent choices concerning the material, layup, and structural configurations so that a more efficient structure with structural integrity can be designed and built. Considerable progress was made towards achieving these goals via this work. In regard to damage tolerance considerations, the following were identified as important effects: composite layup and associated orthotropy/structural anisotropy, specifics of initial local damage mechanisms, role of longitudinal versus hoop stress, and large deformation and associated geometric nonlinearity. Means were established to account for effects of radius and for the nonlinear response. In particular, nondimensional parameters were identified to characterize the importance of nonlinearity in the response of pressurized cylinders. This led to the establishment of a iso-nonlinear-error plot for reference in structural design. Finally, in the case of damage tolerance, the general approach of the original methodology to predict the failure pressure involving extending basic plate failure data by accounting for the local stress intensification was accomplished for the general case by accounting for the mechanisms noted by utilizing the capability of the STAGS finite element code and numerically calculating the local stress intensification for the particular configuration to be considered. For the issue of damage arrest, placement of and configuration of stiffeners (including stiffener curvature), and magnitude and orientation of principal strains due to local bending were found to be key considerations. Means were established to account for stiffener effectiveness quantitatively based on radius, slit size, stiffener curvature' and relative bending stifffiesses involved. Geometric nonlinearity was also found to play an - 24 - important role here. Furthermore, it was determined that damage propagation is controlled by different mechanisms (hoop stress versus flapping stress and the associated factors involved in each) depending upon the direction of damage propagation. This latter item results in an inability to scale these phenomena in one test due to the different factors involved. Finally, the importance of shell curvature and associated instability in response to transverse loading including impact were found to be important considerations in damage resistance. A technique, involving asymmetric meshing of a finite element mesh, was developed to predict this behavior and showed excellent correlation with experimental results. Further details of these ten years of work are presented herein with references made to the fourteen documents produced during this work where full details can be found. Implications of this work are discussed and recommendations made. Although it is clear that there is more work to be done to fully understand composite fuselage technology and specifically the overall issue of damage in primary load-bearing composite structures, important understanding and capability has been extended via this work.

Lagace, Paul A.↗

NASA NTRS2005

Spinning Spacecraft Attitude Estimation Using Markley Variables: Filter Implementation And Results

Attitude estimation is often more difficult for spinning spacecraft than for three-axis stabilized platforms due to the need to follow rapidly-varying state vector elements and the lack of three-axis rate measurements from gyros. The estimation problem simplifies when torques are negligible and nutation has damped out, but the general case requires a sequential filter with dynamics propagation. This paper describes the implementation and test results for an extended Kalman filter for spinning spacecraft attitude and rate estimation based on a novel set of variables suggested in a paper by Markley [AAS93-3301 (referred to hereafter as Markley variables). Markley has demonstrated that the new set of variables provides a superior parameterization for numerical integration of the attitude dynamics for spinning or momentum-biased spacecraft. The advantage is that the Markley variables have fewer rapidly-varying elements than other representations such as the attitude quaternion and rate vector. A filter based on these variables was expected to show improved performance due to the more accurate numerical state propagation. However, for a variety of test cases, it has been found that the new filter, as currently implemented, does not perform significantly better than a quaternion-based filter that was developed and tested in parallel. This paper reviews the mathematical background for a filter based on Markley variables. It also describes some features of the implementation and presents test results. The test cases are based on a mission using magnetometer and Sun sensor data and gyro measurements on two axes normal to the spin axis. The orbit and attitude scenarios and spacecraft parameters are modeled after one of the THEMIS (Time History of Events and Macroscale Interactions during Substorms) probes. Several tests are presented that demonstrate the filter accuracy and convergence properties. The tests include torque-free motion with various nutation angles, large constant-torque attitude slews, sensor misalignments, large initial attitude and rate errors, and cases with low data frequency. It is found that the convergence is rapid, the radius of convergence is large, and the results are reasonably accurate even in the presence of unmodeled perturbations.

Sedlak, Joseph E.↗

NASA NTRS2023

An Evaluation of A Physics-Based Firn Model and A Semi-Empirical Firn Model Across the Greenland Ice Sheet (1980–2020)

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (−66.6 km 3 yr −1 ) in SNOWPACK and 1.5 % (−17.4 km 3 yr −1 ) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage.

Megan Thompson-Munson↗

NASA NTRS2020

Performance of NOAA-20 VIIRS Solar Diffuser and Solar Diffuser Stability Monitor

Visible Infrared Imaging Radiometer Suite (VIIRS) radiometrically calibrates its reflective solar bands (RSBs) primarily through a sunlit onboard solar diffuser (SD). The sunlit SD provides a known radiance under the condition that the absolute product of the SD screen transmittance and the bidirectional reflectance distribution function (BRDF) along the SD-to-telescope direction is accurately known. The BRDF changes due to solar exposure. The change, referred to as the H-factor, is monitored by the onboard SD stability monitor (SDSM). The accuracy of the retrieved H-factor propagates to the retrieved F-factor which corrects the scene spectral radiance. High accuracy of the retrieved H-factor relies on high accuracies in the SDSM screen relative effective transmittance and the relative product of the SD screen effective transmittance and the BRDF at the mission start, and a high SDSM detector signal-to-noise ratio (SNR). This article briefly reviews the algorithms used for the NOAA-20 (N20) VIIRS RSB on-orbit radiometric calibration. Additionally, we show the performance of the N20 VIIRS SDSM, giving the SDSM detector SNRs and the SDSM detector gain temporal changes. We develop a model for the SNRs. The model shows that the decreased SNRs in time are due to the detector gain decreases. We also show the N20 VIIRS SD on-orbit performance, measured by the retrieved H-factor and the estimated standard deviation of its error. The H-factor for the telescope view is obtained from the H-factor for the SDSM view, multiplied by an H-factor angular dependence term. We use an innovative method to determine the angular dependence, using the dependence obtained for the Suomi National Polar-orbiting Partnership (SNPP) VIIRS.

VIIRS↗

NASA NTRS2017

A Minimal State Augmentation Algorithm for Vision-Based Navigation without Using Mapped Landmarks

This paper describes MAVeN (Minimal State Augmentation Algorithm for Vision-Based Navigation), which is a new algorithm for vision-based navigation that has only 21 states, yet is able to track features in successive camera images and use them to propagate estimates of the spacecraft position and velocity. The filter dimension drops to 12 if attitude information is already available. The low filter dimension makes MAVeN a very reliable and practical algorithm for real-time flight implementation. The main idea is to project observed features onto a rough shape model of the ground surface, which are then used by the filter as pseudo-landmarks. The shape model is assumed to be known beforehand, as would be obtained from prior surveillance of the landing site from orbit. MAVeN does not require pre-mapped landmarks, so it is able to navigate terrain that has not been previously observed up close. This property is especially important for close proximity operations in small body missions where ground surface features are being seen for the first time at close range. MAVeN is also able to hover motionless above the ground without position error growth, which is unusual for this class of vision-based navigation algorithms.

San Martin, A. Miguel↗

NASA NTRS1994

Measuring rainwater content by radar using propagation differential phase shift

While radars measure several quantities closely coupled to the rainfall rate, for frequencies less than 15 GHz, estimates of the rainwater content W are traditionally computed from the radar reflectivity factor Z or the rate of attenuation A--quantities only weakly related to W. Consequently, instantaneous point estimates of W using Z and A are often erroneous. A more natural, alternative parameter for estimating W at these frequencies is the specific polarization propagation differential phase shift phi(sub DP), which is a measure of the change in the difference between phases of vertically (V) and horizontally (H) polarized waves with increasing distance from a radar. It is now well known that W is nearly linearly related to phi(sub DP) divided by (1 - reversed R), where reversed R is the mass-weighted mean axis ratio of the raindrops. Unfortunately, such relations are not widely used in part because measurements of phi(sub DP) are scarce but also because one must determine reversed R. In this work it is shown that this parameter can be estimated using the differential reflectivity (Z(sub H)/Z(sub V) at 3 GHz. An alternative technique is suggested for higher frequencies when the differential reflectivity becomes degraded by attenuation. While theory indicates that it should be possible using phi(sub DP) to estimate W quite accurately, measurement errors increase the uncertainty to +/- 18%-35% depending on reversed R. While far from ideal, it appears that these estimates are likely to be considerably more accurate than those deduced using currently available methods.

Jameson, A. R.↗

NASA NTRS2000

Hurricane Wind Field Measurements with Scanning Airborne Doppler Lidar During CAMEX-3

During the 1998 Convection and Moisture Experiment (CAMEX-3), the first hurricane wind field measurements with Doppler lidar were achieved. Wind fields were mapped within the eye, along the eyewall, in the central dense overcast, and in the marine boundary layer encompassing the inflow region. Spatial coverage was determined primarily by cloud distribution and opacity. Within optically-thin cirrus slant range of 20- 25 km was achieved, whereas no propagation was obtained during penetration of dense cloud. Measurements were obtained with the Multi-center Airborne Coherent Atmospheric Wind Sensor (MACAWS) on the NASA DC-8 research aircraft. MACAWS was developed and operated cooperatively by the atmospheric lidar remote sensing groups of NOAA Environmental Technology Laboratory, NASA Marshall Space Flight Center, and Jet Propulsion Laboratory. A pseudo-dual Doppler technique ("co-planar scanning") is used to map the horizontal component of the wind at several vertical levels. Pulses from the laser are directed out the left side of the aircraft in the desired directions using computer-controlled rotating prisms. Upon exiting the aircraft, the beam is completely eyesafe. Aircraft attitude and speed are taken into account during real-time signal processing, resulting in determination of the ground-relative wind to an accuracy of about 1 m/s magnitude and about 10 deg direction. Beam pointing angle errors are about 0.1 deg, equivalent to about 17 m at 10 km. Horizontal resolution is about 1 km (along-track) for typical signal processor and scanner settings; vertical resolution varies with range. Results from CAMEX-3 suggest that scanning Doppler wind lidar can complement airborne Doppler radar by providing wind field measurements in regions that are devoid of hydrometeors. At present MACAWS observations are being assimilated into experimental forecast models and satellite Doppler wind lidar simulations to evaluate the relative impact.

Rothermel, Jeffry↗

NASA NTRS2002

A Two-Dimensional Linear Bicharacteristic Scheme for Electromagnetics

The upwind leapfrog or Linear Bicharacteristic Scheme (LBS) has previously been implemented and demonstrated on one-dimensional electromagnetic wave propagation problems. This memorandum extends the Linear Bicharacteristic Scheme for computational electromagnetics to model lossy dielectric and magnetic materials and perfect electrical conductors in two dimensions. This is accomplished by proper implementation of the LBS for homogeneous lossy dielectric and magnetic media and for perfect electrical conductors. Both the Transverse Electric and Transverse Magnetic polarizations are considered. Computational requirements and a Fourier analysis are also discussed. Heterogeneous media are modeled through implementation of surface boundary conditions and no special extrapolations or interpolations at dielectric material boundaries are required. Results are presented for two-dimensional model problems on uniform grids, and the Finite Difference Time Domain (FDTD) algorithm is chosen as a convenient reference algorithm for comparison. The results demonstrate that the two-dimensional explicit LBS is a dissipation-free, second-order accurate algorithm which uses a smaller stencil than the FDTD algorithm, yet it has less phase velocity error.

Beggs, John H.↗

NASA NTRS1999

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

NASA NTRS2005

System for Processing Coded OFDM Under Doppler and Fading

An advanced communication system has been proposed for transmitting and receiving coded digital data conveyed as a form of quadrature amplitude modulation (QAM) on orthogonal frequency-division multiplexing (OFDM) signals in the presence of such adverse propagation-channel effects as large dynamic Doppler shifts and frequency-selective multipath fading. Such adverse channel effects are typical of data communications between mobile units or between mobile and stationary units (e.g., telemetric transmissions from aircraft to ground stations). The proposed system incorporates novel signal processing techniques intended to reduce the losses associated with adverse channel effects while maintaining compatibility with the high-speed physical layer specifications defined for wireless local area networks (LANs) as the standard 802.11a of the Institute of Electrical and Electronics Engineers (IEEE 802.11a). OFDM is a multi-carrier modulation technique that is widely used for wireless transmission of data in LANs and in metropolitan area networks (MANs). OFDM has been adopted in IEEE 802.11a and some other industry standards because it affords robust performance under frequency-selective fading. However, its intrinsic frequency-diversity feature is highly sensitive to synchronization errors; this sensitivity poses a challenge to preserve coherence between the component subcarriers of an OFDM system in order to avoid intercarrier interference in the presence of large dynamic Doppler shifts as well as frequency-selective fading. As a result, heretofore, the use of OFDM has been limited primarily to applications involving small or zero Doppler shifts. The proposed system includes a digital coherent OFDM communication system that would utilize enhanced 802.1la-compatible signal-processing algorithms to overcome effects of frequency-selective fading and large dynamic Doppler shifts. The overall transceiver design would implement a two-frequency-channel architecture (see figure) that would afford frequency diversity for reducing the adverse effects of multipath fading. By using parallel concatenated convolutional codes (also known as Turbo codes) across the dual-channel and advanced OFDM signal processing within each channel, the proposed system is intended to achieve at least an order of magnitude improvement in received signal-to-noise ratio under adverse channel effects while preserving spectral efficiency.

Tsou, Haiping↗

NASA NTRS2012

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements.

Williams, Christopher A.↗