Search NASA⌕ Search

SEARCH · Search NASA

Results for “POINT MATCHING METHOD”

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

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

At least 235 records · Page 13

An Optimal Orthogonal Decomposition Method for Kalman Filter-Based Turbofan Engine Thrust Estimation

A new linear point design technique is presented for the determination of tuning parameters that enable the optimal estimation of unmeasured engine outputs, such as thrust. The engine s performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters related to each major engine component. Accurate thrust reconstruction depends on knowledge of these health parameters, but there are usually too few sensors to be able to estimate their values. In this new technique, a set of tuning parameters is determined that accounts for degradation by representing the overall effect of the larger set of health parameters as closely as possible in a least-squares sense. The technique takes advantage of the properties of the singular value decomposition of a matrix to generate a tuning parameter vector of low enough dimension that it can be estimated by a Kalman filter. A concise design procedure to generate a tuning vector that specifically takes into account the variables of interest is presented. An example demonstrates the tuning parameters ability to facilitate matching of both measured and unmeasured engine outputs, as well as state variables. Additional properties of the formulation are shown to lend themselves well to diagnostics.

Litt, Jonathan S.↗

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.

Su, Tzu-Jeng↗

Integrating Wind Profiling Radars and Radiosonde Observations with Model Point Data to Develop a Decision Support Tool to Assess Upper-Level Winds for Space Launch

On the day-of-launch, the 45th Weather Squadron (45 WS) Launch Weather Officers (LWOs) monitor the upper-level winds for their launch customers to include NASA's Launch Services Program and NASA's Ground Systems Development and Operations Program. They currently do not have the capability to display and overlay profiles of upper-level observations and numerical weather prediction model forecasts. The LWOs requested the Applied Meteorology Unit (AMU) develop a tool in the form of a graphical user interface (GUI) that will allow them to plot upper-level wind speed and direction observations from the Kennedy Space Center (KSC) 50 MHz tropospheric wind profiling radar, KSC Shuttle Landing Facility 915 MHz boundary layer wind profiling radar and Cape Canaveral Air Force Station (CCAFS) Automated Meteorological Processing System (AMPS) radiosondes, and then overlay forecast wind profiles from the model point data including the North American Mesoscale (NAM) model, Rapid Refresh (RAP) model and Global Forecast System (GFS) model to assess the performance of these models. The AMU developed an Excel-based tool that provides an objective method for the LWOs to compare the model-forecast upper-level winds to the KSC wind profiling radars and CCAFS AMPS observations to assess the model potential to accurately forecast changes in the upperlevel profile through the launch count. The AMU wrote Excel Visual Basic for Applications (VBA) scripts to automatically retrieve model point data for CCAFS (XMR) from the Iowa State University Archive Data Server (http://mtarchive.qeol.iastate.edu) and the 50 MHz, 915 MHz and AMPS observations from the NASA/KSC Spaceport Weather Data Archive web site (http://trmm.ksc.nasa.gov). The AMU then developed code in Excel VBA to automatically ingest and format the observations and model point data in Excel to ready the data for generating Excel charts for the LWO's. The resulting charts allow the LWOs to independently initialize the three models 0-hour forecasts against the observations to determine which is the best performing model and then overlay the model forecasts on time-matched observations during the launch countdown to further assess the model performance and forecasts. This paper will demonstrate integration of observed and predicted atmospheric conditions into a decision support tool and demonstrate how the GUI is implemented in operations.

Bauman, William H., III↗

Algorithms for parallel flow solvers on message passing architectures

The purpose of this project has been to identify and test suitable technologies for implementation of fluid flow solvers -- possibly coupled with structures and heat equation solvers -- on MIMD parallel computers. In the course of this investigation much attention has been paid to efficient domain decomposition strategies for ADI-type algorithms. Multi-partitioning derives its efficiency from the assignment of several blocks of grid points to each processor in the parallel computer. A coarse-grain parallelism is obtained, and a near-perfect load balance results. In uni-partitioning every processor receives responsibility for exactly one block of grid points instead of several. This necessitates fine-grain pipelined program execution in order to obtain a reasonable load balance. Although fine-grain parallelism is less desirable on many systems, especially high-latency networks of workstations, uni-partition methods are still in wide use in production codes for flow problems. Consequently, it remains important to achieve good efficiency with this technique that has essentially been superseded by multi-partitioning for parallel ADI-type algorithms. Another reason for the concentration on improving the performance of pipeline methods is their applicability in other types of flow solver kernels with stronger implied data dependence. Analytical expressions can be derived for the size of the dynamic load imbalance incurred in traditional pipelines. From these it can be determined what is the optimal first-processor retardation that leads to the shortest total completion time for the pipeline process. Theoretical predictions of pipeline performance with and without optimization match experimental observations on the iPSC/860 very well. Analysis of pipeline performance also highlights the effect of uncareful grid partitioning in flow solvers that employ pipeline algorithms. If grid blocks at boundaries are not at least as large in the wall-normal direction as those immediately adjacent to them, then the first processor in the pipeline will receive a computational load that is less than that of subsequent processors, magnifying the pipeline slowdown effect. Extra compensation is needed for grid boundary effects, even if all grid blocks are equally sized.

Vanderwijngaart, Rob F.↗

Planning Satellite Swarm Measurements for Climate Models: Comparing Dynamic Constraint Processing and MILP Methods

We present D-SHIELD, a challenging climate science application to plan coordinated measurements (observations) for a constellation of satellites, each containing two different sensors, each with 61 pointing angle options. The L-band and P-band radar sensors collect data fed into a soil moisture model which tracks and predicts soil moisture across 1.67 million Ground Positions (GP). Soil moisture is an important predictor of wildfires, and then a predictor of floods, landslides and debris flow after a fire. Each measurement covers multiple GP due to the sensor footprint. Each GP has a "model error" which represents the uncertainty of the the soil moisture state prediction. Model error changes at different rates for each GP as the time since last observation increases and after significant events like rain. The planner's goal is to select measurements which maximize soil moisture model improvement (reduce model uncertainty). This problem is combinatorically explosive, involving many degrees of freedom for planner choices. Good domain heuristics can find solutions within a reasonable time for our application needs but cannot be proven optimal. In this paper we compare two different planning approaches to this problem: Dynamic Constraint Processing (DCP) and Mixed Integer Linear Programming (MILP). We match inputs and metrics for both DCP and MILP algorithms to enable a direct apples-to-apples comparison. We demonstrate and discuss the trades between DCP flexibility and performance vs. MILP's promise of provable optimality.

Rich Levinson↗

MANN: A program to transfer designs for diffractive optical elements to a MANN photolithographic mask generator

There are two basic areas of interest for diffractive optics. In the first, the property of wavefront division is exploited for achieving optical fanout, analogous to the more familiar electrical fanout of electronic circuitry. The basic problem here is that when using a simple uniform diffraction grating the energy input is divided unevenly among the output beams. The other area of interest is the use of diffractive elements to replace or supplement standard refractive elements such as lenses. Again, local grating variations can be used to control the amount of bending imparted to optical rays, and the efficiency of the diffractive element will depend on how closely the element can be matched to the design requirements. In general, production restrictions limit how closely the element approaches the design, and for the common case of photolithographic production, a series of binary masks is required to achieve high efficiency. The actual design process is much more involved than in the case of elements for optical fanout, as the desired phase of the optical wavefront over some reference plane must be specified and the phase alteration to be introduced at each point by the diffraction element must be known. This generally requires the utilization of a standard optical design program. Two approaches are possible. In the first approach, the diffractive element is treated as a special type of lens and the ordinary optical design equations are used. Optical design programs tend to follow a second approach, namely, using the equations of optical interference derived from holographic theory and then allowing the introduction of phase front corrections in the form of polynomial equations. By using either of these two methods, diffractive elements can be used not only to compensate for distortions such as chromatic or spherical aberration, but also to perform the work of a variety of other optical elements such as null correctors, beam shapers, etc. The main focus of the project described in this report is how the design information from the lens design program is incorporated into the photolithographic process. It is shown that the MANN program, a photolithographic mask generator, fills the need for a link between lens design programs and mask generation controllers.The generated masks can be used to expose a resist-coated substrate which is etched and then must be re-coated, re-exposed, and re-etched for making copies, just as in the electronics industry.

Matthys, Donald R.↗

A Match Made in Space

Just before the space shuttle reaches orbit, its three main engines shut down so that it can achieve separation from the massive external tank that provided the fuel required for liftoff and ascent. In jettisoning the external tank, which is completely devoid of fuel at this point in the flight, the space shuttle fires a series of thrusters, separate from its main engines, that gives the orbiter the maneuvering ability necessary to safely steer clear of the descending tank and maintain its intended flight path. These thrusters make up the space shuttle s Reaction Control System. While the space shuttle s main engines only provide thrust in one direction (albeit a very powerful thrust), the Reaction Control System engines allow the vehicle to maneuver in any desired direction (via small amounts of thrust). The resulting rotational maneuvers are known as pitch, roll, and yaw, and are very important in ensuring that the shuttle docks properly when it arrives at the International Space Station and safely reenters the Earth s atmosphere upon leaving. To prevent the highly complex Reaction Control System from malfunctioning during space shuttle flights, and to provide a diagnosis if such a mishap were to occur, NASA turned to a method of artificial intelligence that truly defied the traditional laws of computer science.

Source record↗

Measurement of Residual Flexibility for Substructures Having Prominent Flexible Interfaces

Verification of a dynamic model of a constrained structure requires a modal survey test of the physical structure and subsequent modification of the model to obtain the best agreement possible with test data. Constrained-boundary or fixed-base testing has historically been the most common approach for verifying constrained mathematical models, since the boundary conditions of the test article are designed to match the actual constraints in service. However, there are difficulties involved with fixed-base testing, in some cases making the approach impractical. It is not possible to conduct a truly fixed-base test due to coupling between the test article and the fixture. In addition, it is often difficult to accurately simulate the actual boundary constraints, and the cost of designing and constructing the fixture may be prohibitive. For use when fixed-base testing proves impractical or undesirable, alternate free-boundary test methods have been investigated, including the residual flexibility technique. The residual flexibility approach has been treated analytically in considerable detail and has had limited frequency response measurements for the method. This concern is well-justified for a number of reasons. First, residual flexibilities are very small numbers, typically on the order of 1.0E-6 in/lb for translational diagonal terms, and orders of magnitude smaller for off-diagonal values. This poses difficulty in obtaining accurate and noise-free measurements, especially for points removed from the excitation source. A second difficulty encountered in residual measurements lies in obtaining a clean residual function in the process of subtracting synthesized modal data from a measured response function. Inaccuracies occur since modes are not subtracted exactly, but only to the accuracy of the curve fits for each mode; these errors are compounded with increasing distance from the excitation point. In this paper, the residual flexibility method is applied to a simple structure in both test and analysis. Measured and predicted residual functions are compared, and regions of poor data in the measured curves are described. It is found that for accurate residual measurements, frequency response functions having prominent stiffness lines in the acceleration/force format are needed. The lack of such stiffness lines increases measurement errors. Interface drive point frequency respose functions for shuttle orbiter payloads exhibit dominant stiffness lines, making the residual test approach a good candidate for payload modal tests when constrained tests are inappropriate. Difficulties in extracting a residual flexibility value from noisy test data are discussed. It is shown that use of a weighted second order least-squares curve fit of the measured residual function allows identification of residual flexibility that compares very well with predictions for the simple structure. This approach also provides an estimate of second order residual mass effects.

Tinker, Michael L.↗

Nanoflare Heating Frequency of an X-Ray Bright Point Observed By MaGIXS

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multimillion kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with concurrent observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the potential field extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the loop morphology and emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 400 s to 800 s, which strongly suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

coronal heating↗

Determining The Nanoflare Heating Frequency of an X-Ray Bright Point Observed by MaGIXS

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with nearest-available observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the linear-force-free extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the observed emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 1500 s to 3000 s, which suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

coronal heating↗

The Spatial Vision Tree: A Generic Pattern Recognition Engine- Scientific Foundations, Design Principles, and Preliminary Tree Design

New foundational ideas are used to define a novel approach to generic visual pattern recognition. These ideas proceed from the starting point of the intrinsic equivalence of noise reduction and pattern recognition when noise reduction is taken to its theoretical limit of explicit matched filtering. This led us to think of the logical extension of sparse coding using basis function transforms for both de-noising and pattern recognition to the full pattern specificity of a lexicon of matched filter pattern templates. A key hypothesis is that such a lexicon can be constructed and is, in fact, a generic visual alphabet of spatial vision. Hence it provides a tractable solution for the design of a generic pattern recognition engine. Here we present the key scientific ideas, the basic design principles which emerge from these ideas, and a preliminary design of the Spatial Vision Tree (SVT). The latter is based upon a cryptographic approach whereby we measure a large aggregate estimate of the frequency of occurrence (FOO) for each pattern. These distributions are employed together with Hamming distance criteria to design a two-tier tree. Then using information theory, these same FOO distributions are used to define a precise method for pattern representation. Finally the experimental performance of the preliminary SVT on computer generated test images and complex natural images is assessed.

Rahman, Zia-ur↗

Development and Applications of a Stage Stacking Procedure

The preliminary design of multistage axial compressors in gas turbine engines is typically accomplished with mean-line methods. These methods, which rely on empirical correlations, estimate compressor performance well near the design point, but may become less reliable off-design. For land-based applications of gas turbine engines, off-design performance estimates are becoming increasingly important, as turbine plant operators desire peaking or load-following capabilities and hot-day operability. The current work develops a one-dimensional stage stacking procedure, including a newly defined blockage term, which is used to estimate the off-design performance and operability range of a 13-stage axial compressor used in a power generating gas turbine engine. The new blockage term is defined to give mathematical closure on static pressure, and values of blockage are shown to collapse to curves as a function of stage inlet flow coefficient and corrected shaft speed. In addition to these blockage curves, the stage stacking procedure utilizes stage characteristics of ideal work coefficient and adiabatic efficiency. These curves are constructed using flow information extracted from computational fluid dynamics (CFD) simulations of groups of stages within the compressor. Performance estimates resulting from the stage stacking procedure are shown to match the results of CFD simulations of the entire compressor to within 1.6% in overall total pressure ratio and within 0.3 points in overall adiabatic efficiency. Utility of the stage stacking procedure is demonstrated by estimation of the minimum corrected speed which allows stable operation of the compressor. Further utility of the stage stacking procedure is demonstrated with a bleed sensitivity study, which estimates a bleed schedule to expand the compressors operating range.

Kulkarni, Sameer↗

Inner Front Pole Cover Erosion in the 12.5 kW HERMeS Hall Thruster Over a Range of Operating Conditions

Erosion characteristics on the cover of the inner front pole in a 12.5 kW Hall thruster were measured over a wide range of operating conditions in tests of 6 to 14 hours duration using an accelerated test method and a very sensitive radioactive tracer-based erosion diagnostic. The operating points included the nominal 300 - 600 V conditions on a constant 20.8 A throttle curve, but included additional conditions at other currents spanning the throttling envelope and measurements at varying magnetic eld strength, facility pressure, and discharge voltage oscillation amplitude. The results show that the 300 V condition produces the highest wear rates on the 20.8 A throttle curve, but that rates actually increase with decreasing current. The wear rate was insensitive to discharge voltage ripple, but increased monotonically with magnetic eld strength, particularly near the inner radius of the pole cover. The inner region was also sensitive to facility pressure, showing lower rates at a higher pressure level. Separate experiments in which the energy distributions of ions generated by the hollow cathode were measured suggest that the cathode plume may be a source of energetic ions responsible for some of the erosion trends, in addition to ions originating in the thruster plume. The Hall thruster simulation code Hall2De is able to reproduce the erosion characteristics observed at 600 V, 20.8 A, but cannot currently match the rates measured at lower voltages and currents.

Lopez Ortega, Alejandro↗

Handheld White Light Interferometer for Measuring Defect Depth in Windows

Accurate quantification of defects (scratches and impacts) is vital to the certification of flight hardware and other critical components. The amount of damage to a particular component contributes to the performance, reliability, and safety of a system, which ultimately affects the success or failure of a mission or test. The launch-commit criteria on a Space Shuttle Orbiter window are governed by the depth of the defects that are identified by a visual inspection. This measurement of a defect is not easy to obtain given the environment, size of the defect, and location of the window(s). The determination of depth has typically been performed by taking a mold impression and measuring the impression with an optical profiling instrument. Another method of obtaining an estimate of the depth is by using a refocus microscope. To use a refocus microscope, the surface of the glass and bottom of the defect are, in turn, brought into focus by the operator. The amount of movement between the two points corresponds to the depth of the defect. The refocus microscope requires a skilled operator and has been proven to be unreliable when used on Orbiter windows. White light interferometry was chosen as a candidate to replace the refocus microscope. The White Light Interferometer (WLI) was developed to replace the refocus microscope as the instrument used for measuring the depth of defects in Orbiter windows. The WLI consists of a broadband illumination source, interferometer, detector, motion control, displacement sensor, mechanical housing, and support electronics. The illumination source for the WLI is typically a visible light emitting diode (LED) or a near-infrared superluminescent diode (SLD) with power levels of less than a milliwatt. The interferometer is a Michelson configuration consisting of a 1-in. (2.5-cm) cube beam splitter, a 0.5-in. (1.3-cm) optical window as a movable leg (used to closely match the return intensity of the fixed leg from the window), and a mirrored prism to fold the optics into the mechanical housing. The detector may be one of many C-mount CCD (charge-coupled device) cameras. Motion is provided by a commercial nanostepping motor with a serial interface. The displacement sensor is a custom device specifically designed for this application. The mechanical housing and support electronics were designed to integrate the various components into an instrument that could be physically handled by a technician and easily transported.

Youngquist, Robert↗

TPSAS-NF1676L-16988-DND

There has been renewed interest to uniformly recalibrate historical geostationary (GEO) data records to aid in climate monitoring. GEO sensors have annual repeatable angular sampling over a given location. The view angle is fixed and the imaging schedule is usually constant through its lifetime. Given the fact that colocated GEOs are always share the same sub-satellite point and maintain their imaging schedules provides repeatable angular sampling over decades. One of the biggest challenges in transferring a reference sensor calibration using invariant desert targets to another sensor is the accuracy of the bidirectional reflectance distribution function (BDRF). A well-calibrated GEO can be used to predict the daily exoatmospheric radiance model (DERM) over a desert target for a given GMT that is valid for any GEO sensor at the same location. The advantage of this method is that a BRDF is not needed. Another challenge of invariant target calibration is the unique spectra signature of the desert. However, since most GEOs are built in batches, the spectral response functions (SRF) are very similar for most historical GEOs, the spectral band adjustment factor (SBAF) between GEO sensors is much smaller than for MODIS and GEO sensors. Since the water vapor burden over the desert is seasonal, both the TOA and desert surface can be considered invariant for a given day of the year. The reference GEO can be inter-calibrated with MODIS or VIIRS, which have onboard visible calibration using solar diffusers, using other methods, such as ray-matching or deep convective clouds. Also the next generation GEOs will have onboard visible calibration, which will increase the accuracy of this method. Three Meteosats over the Libyan desert will be used to illustrate the DERM method. The reference Meteosat will be inter-calibrated against Aqua-MODIS. The reference GEO DERM will be constructed and used to calibrate the remaining Meteosats. The DERM calibration will be validated by comparing the calibration using Aqua-MODIS ray-matching. Similarly, two GOES sensors using the Sonoran desert will also be highlighted. An uncertainty analysis will also be performed with emphasis on the SBAF, derived over the desert targets using both SCIAMACHY and Hyperion hyper-spectral radiances.

David Doelling↗

Non-Contact Thrust Stand Calibration Method for Repetitively-Pulsed Electric Thrusters

A thrust stand calibration technique for use in testing repetitively-pulsed electric thrusters for in-space propulsion has been developed and tested using a modified hanging pendulum thrust stand. In the implementation of this technique, current pulses are applied to a solenoidal coil to produce a pulsed magnetic field that acts against the magnetic field produced by a permanent magnet mounted to the thrust stand pendulum arm. The force on the magnet is applied in this non-contact manner, with the entire pulsed force transferred to the pendulum arm through a piezoelectric force transducer to provide a time-accurate force measurement. Modeling of the pendulum arm dynamics reveals that after an initial transient in thrust stand motion the quasisteady average deflection of the thrust stand arm away from the unforced or zero position can be related to the average applied force through a simple linear Hooke s law relationship. Modeling demonstrates that this technique is universally applicable except when the pulsing period is increased to the point where it approaches the period of natural thrust stand motion. Calibration data were obtained using a modified hanging pendulum thrust stand previously used for steady-state thrust measurements. Data were obtained for varying impulse bit at constant pulse frequency and for varying pulse frequency. The two data sets exhibit excellent quantitative agreement with each other as the constant relating average deflection and average thrust match within the errors on the linear regression curve fit of the data. Quantitatively, the error on the calibration coefficient is roughly 1% of the coefficient value.

Wong, Andrea R.↗

Status of the NASA GMAO Observing System Simulation Experiment

An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.

OSSE↗

Development And Test of A Digitally Steered Antenna Array for The Navigator GPS Receiver

Global Positioning System (GPS)-based navigation has become common for low-Earth orbit spacecraft as the signal environment is similar to that on the Earth s surface. The situation changes abruptly, however, for spacecraft whose orbital altitudes exceed that of the GPS constellation. Visibility is dramatically reduced and signals that are present may be very weak and more susceptible to interference. GPS receivers effective at these altitudes require increased sensitivity, which often requires a high-gain antenna. Pointing such an antenna can pose a challenge. One efficient approach to mitigate these problems is the use of a digitally steered antenna array. Such an antenna can optimally allocate gain toward desired signal sources and away from interferers. This paper presents preliminary results in the development and test of a digitally steered antenna array for the Navigator GPS research program at NASA s Goddard Space Flight Center. In particular, this paper highlights the development of an array and front-end electronics, the development and test of a real-time software GPS receiver, and implementation of three beamforming methods for combining the signals from the array. Additionally, this paper discusses the development of a GPS signal simulator which produces digital samples of the GPS L1C/A signals as they would be received by an arbitrary antenna array configuration. The simulator models transmitter and receiver dynamics, near-far and multipath interference, and has been a critical component in both the development and test of the GPS receiver. The GPS receiver system was tested with real and simulated GPS signals. Preliminary results show that performance improvement was achieved in both the weak signal and interference environments, matching analytical predictions. This paper summarizes our initial findings and discusses the advantages and limitations of the antenna array and the various beamforming methods.

Pinto, Heitor David↗