Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic resolution adjustment”

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 55 records · Page 3

Version 4 of the SMAP Level-4 Soil Moisture Algorithm and Data Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with a mean latency of ~2.5 days. The underlying L4_SM algorithm assimilates SMAP radiometer brightness temperature (Tb) observations into the NASA Catchment land surface model using a spatially-distributed ensemble Kalman filter. Version 4 of the L4_SM modeling system includes a reduction in the upward recharge of surface soil moisture from below under non-equilibrium conditions, resulting in reduced bias and improved dynamic range of L4_SM surface soil moisture compared to earlier versions. This change and additional technical modifications to the system reduce the mean and standard deviation of the observation-minus-forecast Tb residuals and overall soil moisture analysis increments while maintaining the skill of the L4_SM soil moisture estimates versus independent in situ measurements; the average, bias-adjusted RMSE in Version 4 is 0.039 m(exp 3) m(exp -3) for surface and 0.026 m(exp 3) m(exp -3) for root-zone soil moisture. Moreover, the coverage of assimilated SMAP observations in Version 4 is near-global owing to the use of additional satellite Tb records for algorithm calibration. L4_SM soil moisture uncertainty estimates are biased low (by 0.01-0.02 m(exp 3) m(exp -3)) against actual errors (computed versus in situ measurements). L4_SM runoff estimates, an additional product of the L4_SM algorithm, are biased low (by 35 mm year (exp -1)) against streamflow measurements. Compared to Version 3, bias in Version 4 is reduced by 46% for surface soil moisture uncertainty estimates and by 33% for runoff estimates.

RMSE↗

Changes in the Mass Balance of the Greenland Ice Sheet in a Warming Climate During 2003-2009

Mass changes of the Greenland ice sheet (GIS) derived from ICESat and GRACE data both show that the net mass loss from GIS during 2003-2009 is about 175 Gt/year, which contributes 0.5mm/yr global sea-level rise. The rate of mass loss has increased significantly since the 1990's when the GIS was close to mass balance. Even though the GIS was close to mass balance during the 1990's, it was already showing characteristics of responding to8 warmer climate, specifically thinning at the margins and thickening inland at higher elevations. During 2003-2009, increased ice thinning due to increases in melting and acceleration of outlet glaciers began to strongly exceed the inland thickening from increases in accumulation. Over the entire GIS, the mass loss between the two periods, from increased melting and ice dynamics, increased by about 190 Gt/year while the mass gain, from increased precipitation and accumulation, increased by only about 15Gt/year. These ice changes occurred during a time when the temperature on GIS changed at rate of about 2K/decade. The distribution of elevation and mass changes derived from ICESat have high spatial resolution showing details over outlet glaciers, by drainage systems, and by elevation. However, information on the seasonal cycle of changes from ICESat data is limited, because the ICESat lasers were only operated during two to three campaigns per year of about 35 days duration each. In contrast, the temporal resolution of GRACE data, provided by the continuous data collection, is much better showing details of the seasonal cycle and the inter-annual variability. The differing sensitivity of the ICESat altimetry and the GRACE gravity methods to motion of the underlying bedrock from glacial isostatic adjustment (GIA) is used to evaluate the GIA corrections provided by models. The two data types are also combined to make estimates of the partitioning of the mass gains and losses among accumulation, melting, and ice discharge from outlet glaciers.

Zwally, H. Jay↗

Trajectory design for the Deep Space Program Science Experiment (DSPSE) mission

In 1994, the Deep Space Program Science Experiment (DSPSE) spacecraft will become the first spacecraft to perform, in succession, both a lunar orbiting mission and a deep-space asteroid encounter mission. The primary mission objective is to perform a long-duration flight-test of various new-technology lightweight components, such as sensors, in a deep-space environment. The mission has two secondary science objectives: to provide high-resolution imaging of the entire lunar surface for mapping purposes and flyby imaging of the asteroid 1620 Geographos. The DSPSE mission is sponsored by the Strategic Defense Initiative Organization (SDIO). As prime contractor, the Naval Research Laboratory (NRL) is building the spacecraft and will conduct mission operations. The Goddard Space Flight Center's (GSFC) Flight Dynamics Division is supporting NRL in the areas of The Deep Space Network (DSN) will provide tracking support. The DSPSE mission will begin with a launch from the Western Test Range in late January 1994. Following a minimum 1.5-day stay in a low-Earth parking orbit, a solid kick motor burn will boost DSPSE into an 18-day, 2.5-revolution phasing orbit transfer trajectory to the Moon. Two burns to insert DSPSE into a lunar polar orbit suitable for the mapping mission will be followed by mapping orbit maintenance and adjustment operations over a period of 2 sidereal months. In May 1994, a lunar orbit departure maneuver, in conjunction with a lunar swingby 26 days later, will propel DSPSE onto a heliocentric transfer that will intercept Geographos on September 1, 1994. This paper presents the characteristics, deterministic delta-Vs, and design details of each trajectory phase of this unique mission, together with the requirements, constraints, and design considerations to which each phase is subject. Numerous trajectory plots and tables of significant trajectory events are included. Following a discussion of the results of a preliminary launch window analysis, a summary of the deterministic impulsive delta-V budget required to establish the baseline mission trajectory design is presented.

Carrington, D.↗

(abstract) Overview of NASA's Adaptive Structures Program

NASA's research program in Adaptive Structures was initiated seven years ago to provide the technology required for large (20-50 meters in dimension) precision (submicron) structures for observations from space. The current approach of designing thermally stable passive structures cannot meet the requirements and moreover the systems cannot be validated by ground tests. Adaptive Structures provides the capability to adjust the quasistatic dimensions of the structures, preload joints, add active damping, and provide the actuation forces necessary to attenuate the undesired dynamic motions. The technology is applicable to space platforms with various pointing instruments as well as for microspacecraft. This technology has been incorporated as part of a small Active Fold Mirror on the Wide Field Planetary Camera to be flown within the year to correct the optical errors in the Hubble Telescope. Highlights of NASA's Adaptive Structures research and applications will be presented along with future technology and research requirements for both small inexpensive spacecraft and large precision systems.

adaptive↗

Thermal Structure and Dynamics of Saturn's Northern Springtime Disturbance

This article combined several infrared datasets to study the vertical properties of Saturn's northern springtime storm. Spectroscopic observations of Saturn's northern hemisphere at 0.5 and 2.5 / cm spectral resolution were provided by the Cassini Composite Infrared Spectrometer (CIRS, 17). These were supplemented with narrow-band filtered imaging from the ESO Very Large Telescope VISIR instrument (16) to provide a global spatial context for the Cassini spectroscopy. Finally, nightside imaging from the Cassini Visual and Infrared Mapping Spectrometer (VIMS, 22) provided a glimpse of the undulating cloud activity in the eastern branch of the disturbance. Each of these datasets, and the methods used to reduce and analyse them, will be described in detail below. Spatial maps of atmospheric temperatures, aerosol opacity and gaseous distributions are derived from infrared spectroscopy using a suite of radiative transfer and optimal estimation retrieval tools developed at the University of Oxford, known collectively as Nemesis (23). Synthetic spectra created from a reference atmospheric model for Saturn and appropriate sources of spectroscopic line data (6, 24) are convolved with the instrument function for each dataset. Atmospheric properties are then iteratively adjusted until the measurements are accurately reproduced with physically-realistic temperatures, compositions and cloud opacities.

Fletcher, Leigh N.↗

Towards an "All Speed" Unstructured Upwind Scheme

In the authors previous studies [1], a time-accurate, upwind finite volume method (ETAU scheme) for computing compressible flows on unstructured grids was proposed. The scheme is second order accurate in space and time and yields high resolution in the presence of discontinuities. The scheme features a multidimensional limiter and multidimensional numerical dissipation. These help to stabilize the numerical process and to overcome the annoying pathological behaviors of upwind schemes. In the present paper, it will be further shown that such multidimensional treatments also lead to a nearly all-speed or Mach number insensitive upwind scheme. For flows at very high Mach number, e.g., 10, local numerical instabilities or the pathological behaviors are suppressed, while for flows at very low Mach number, e.g., 0.02, computation can be directly carried out without invoking preconditioning. For flows in different Mach number regimes, i.e., low, medium, and high Mach numbers, one only needs to adjust one or two parameters in the scheme. Several examples with low and high Mach numbers are demonstrated in this paper. Thus, the ETAU scheme is applicable to a broad spectrum of flow regimes ranging from high supersonic to low subsonic, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics).

Loh, Ching Y.↗

Dynamical Studies of the Middle Atmosphere Using High Resolution Doppler Imager Observations

This report summarizes the activities of NASA grant NAG5-11068, "Dynamicai Studies of the Middle Atmosphere Using High Resolution Doppler Imager Observations." The High Resolution Doppler Imager (HRDI) on the Upper Atmosphere Research Satellite (UARS) has been providing direct measurements of the Earth's horizontal wind field in the stratosphere, mesosphere and lower thermosphere. Mesospheric temperatures, ozone, and O((sup 1) D) densities, and stratospheric aerosol extinctions coefficients, are also retrieved. The goal of HRDI is to measure the vector winds in the stratosphere (10-40 km), mesosphere, and lower thermosphere (approximately 50-120 km) during the day, and the lower thermosphere at night (approximately 95 km) to an accuracy of 5 m/s. The horizontal wind vector is measured by observing the Doppler shift of rotational lines of molecular oxygen along two lines of sight. In addition to winds, temperatures and volume emission rates are determined in the mesosphere and lower thermosphere, from which ozone and O((sup 1) D) concentrations can be derived, and aerosol scattering coefficients are determined in the stratosphere. UARS was launched on September 12, 1991, into a 585-km circular orbit inclined 57 degrees to the equator HRDI was activated September 28, 1991 and following a period of checkout and adjustment of the instrument parameters, scientific observations began November 2, 199 1. HRDI operated nearly continuously from launch until April 1995. At that time the UARS solar array drive failed, forcing the instruments to time-share the available power. From July 1995 to July 1996 HRDI operated approximately 50% of the time. At that point, one of the three spacecraft batteries failed and from then until September 1998 the duty cycle was less than 20% per month, At that time it was determined that HRDI could operate during each daytime pass, which increased the daytime duty cycle to close to l00%, while nighttime operations were limited to about a week per month. In the fall of 1999, the second tape recorder failed requiring a real time contact with a TRDSS satellite to retrieve that data. This resulted in about 60% data collection efficiency. Finally, in the summer of 2000, the second star sensor failed requiring the spacecraft attitude to be controlled by a three axis magnetometer and sun sensor. This resulted in a loss of attitude knowledge but operations continue with the anticipation of correcting the attitude. A new method for determining the tide and mean structure from satellite data in conjunction with a new tidal model has been devised. For brevity, it shall be referred to as the TMAT or Tide-Mean Assimilation Technique. Most previous methods of tidal analysis are based on various ways of slicing the data set.

Skinner, Wilbert↗

A conjecture on global flow dynamics

Recent observations by the Yohkoh and Ulysses spacecraft have instilled a new appreciation for the role of the ever-evolving geometry and topology of the coronal magnetic structure in shaping and influencing the global solar wind outflow. While many issues concerning the global evolution of large-scale interplanetary flow structures await resolution, we focus here on one specific problem. Namely, we argue that when considering the impact of solar activity on heliosphere scales, rather more attention should be directed toward ascertaining the nature and implications of long-term. Large-scale changes in the corona than toward the more dramatic and highly visible manifestations of solar activity, such as flares and coronal mass ejections (CMEs). That is, while energetic CMEs unquestionably disrupt the interplanetary medium locally and for a time by the broader measure they prompt lasting, wide-spread heliospheric response only when they occur as the agent of substantial, long-term change in the global configuration of the corona. In support, we point to various evidences of coronal change associated with major episodes of heliospheric disruption, and we discuss the limitations imposed by current observation techniques. In addition, we offer dynamical arguments (illustrated by a simple dynamical simulation) that when the coronal configuration changes quickly enough over sufficient spatial scale, a re-adjustment wave linking the initial and final quasi-steady flow states must propagate out from the Sun. This wave is transitory and global, and its amplitude can grow to substantial proportions with heliocentric distance. Moreover, the readjustment wave occurs whether or not the corona needs to expel CMEs to reconfigure itself. We suggest that it is this wave (and not so much the remnants of the initiating CMEs) that constitutes global merged interaction regions in the outer heliosphere.

Pizzo, V. J.↗

Geologic Mapping of Athabasca Valles

Two factors drive us to map the Athabasca Valles area in unusual detail: (1) the extremely well-preserved and exposed surface morphologies and (2) the extensive high resolution imaging. In particular, the near-complete CTX coverage of Athabasca Valles proper and the extensive coverage of its surroundings have been invaluable. The mapping has been done exclusively in ArcGIS, using individual CTX, THEMIS VIS, and MOC frames overlying the THEMIS IR daytime basemap. MOLA shot points and gridded DTMs are also included. It was found that CTX images processed through ISIS are almost always within 300 m of the MOLA derived locations, and usually within tens of meters, with no adjustments to camera pointing. THEMIS VIS images appear to be systematically shifted to the southwest of their correct positions and MOC images are often kilometers off. The good SNR and minimal artifacts make the CTX images vastly more useful than the THEMIS VIS or MOC images. The bulk of the mapping was done at 1:50,000 scale on CTX images. In more complex areas, mapping at 1:24,000 proved necessary. The CTX images were usually simultaneously viewed on a second monitor using the ISIS3 qview program to display the full dynamic range of the CTX data. Where CTX data was not available, mapping was often done at 1:100,000 and most contacts are mapped as approximate.

Keszthelyi, L. P.↗

JPSS-2 VIIRS Pre-Launch Reflective Solar Band Testing and Performance

The Visible Infrared Imaging Radiometer Suite (VIIRS) instruments on-board the Suomi National Polar-orbiting Partnership (S-NPP) and Joint Polar Satellite System (JPSS) spacecrafts 1 and 2 provides calibrated sensor data record (SDR) reflectance, radiance, and brightness temperatures for use in environment data record (EDR) products. The SDRs and EDRs are used in weather forecasting models, weather imagery and climate applications such as ocean color, sea surface temperature and active fires. The VIIRS has 22 bands covering a spectral range 0.4–12.4 m with resolutions of 375 m and 750 m for imaging and moderate bands respectively on four focal planes. The bands are stratified into three different types based on the source of energy sensed by the bands. The reflective solar bands (RSBs) detect sunlight reflected from the Earth, thermal emissive bands (TEBs) sense emitted energy from the Earth and the day/night band (DNB) detects both solar and lunar reflected energy from the Earth. The SDR calibration uses a combination of pre-launch testing and the solar diffuser (SD), on-board calibrator blackbody (OBCBB) and space view (SV) on-orbit calibrator sources. The pre-launch testing transfers the National Institute of Standards and Technology (NIST) traceable calibration to the SD, for the RSB, and the OBCBB, for the TEB. Post-launch, the on-board calibrators track the changes in instrument response and adjust the SDR product as necessary to maintain the calibration. This paper will discuss the pre-launch radiometric calibration portion of the SDR calibration for the RSBs that includes the dynamic range, detector noise, calibration coefficients and radiometric uncertainties for JPSS-2 VIIRS.

JPSS↗

Ensemble Forecasting/Assimilation/Emissions Estimation with WRF-Chem/DART: Accomplishments, Lessons Learned, and Future Plans

Over the past ten years, we have been conducting research on regional ensemble atmospheric composition forecasting/data assimilation/emissions estimation with WRF-Chem/DART. WRF- Chem/DART integrates the Weather Research and Forecasting model (WRF) with online chemistry (WRF-Chem) into the Data Assimilation Research Testbed (DART). DART is an ensemble data assimilation system based on the ensemble adjustment Kalman filter (EAKF) with adaptive inflation, localization (physical and state space), and an optional non-Gaussian formulation of the EAKF. DART includes assimilation of meteorological and limited chemical observations. WRF-Chem/DART extends DART to include assimilation of: MOPITT CO; IASI CO and O3; MODIS AOD; OMI O3, NO2, and SO2; TROPOMI CO, O3, NO2, and SO2, TES CO, CO2 (research mode), O3, NH3, and CH4 (research mode); CrIS CO, O3, NH3, CH4 (research mode), and PAN; SCIAMACHY NO2; GOME2a NO2; MLS O3 and HNO3; and proxy TEMPO O3, and NO2 satellite retrievals as raw retrievals or as ‘compact phase space retrievals’ (CPSRs) for profile retrievals. WRF-Chem/DART also assimilates in situ atmospheric composition measurements and uses the ‘state augmentation method’ for emissions estimation. In our presentation, we will provide an overview of WRF-Chem/DART: • applications and results; • lessons learned from: (i) independent versus joint assimilation; (ii) total/partial column versus profile retrieval assimilation; (iii) joint in situ and retrieval assimilation; (iv) assimilation at grid resolutions ranging from 100 km to 4 km; (v) dynamic emissions estimation; and • future work related to wildfire emissions estimation and intercomparison of CMAQ (online)/JEDI and CMAQ (online and offline)/DART.

WRF-Chem/DART↗

Probabilistic Risk Assessment for Decision Making During Spacecraft Operations

Decisions made during the operational phase of a space mission often have significant and immediate consequences. Without the explicit consideration of the risks involved and their representation in a solid model, it is very likely that these risks are not considered systematically in trade studies. Wrong decisions during the operational phase of a space mission can lead to immediate system failure whereas correct decisions can help recover the system even from faulty conditions. A problem of special interest is the determination of the system fault protection strategies upon the occurrence of faults within the system. Decisions regarding the fault protection strategy also heavily rely on a correct understanding of the state of the system and an integrated risk model that represents the various possible scenarios and their respective likelihoods. Probabilistic Risk Assessment (PRA) modeling is applicable to the full lifecycle of a space mission project, from concept development to preliminary design, detailed design, development and operations. The benefits and utilities of the model, however, depend on the phase of the mission for which it is used. This is because of the difference in the key strategic decisions that support each mission phase. The focus of this paper is on describing the particular methods used for PRA modeling during the operational phase of a spacecraft by gleaning insight from recently conducted case studies on two operational Mars orbiters. During operations, the key decisions relate to the commands sent to the spacecraft for any kind of diagnostics, anomaly resolution, trajectory changes, or planning. Often, faults and failures occur in the parts of the spacecraft but are contained or mitigated before they can cause serious damage. The failure behavior of the system during operations provides valuable data for updating and adjusting the related PRA models that are built primarily based on historical failure data. The PRA models, in turn, provide insight into the effect of various faults or failures on the risk and failure drivers of the system and the likelihood of possible end case scenarios, thereby facilitating the decision making process during operations. This paper describes the process of adjusting PRA models based on observed spacecraft data, on one hand, and utilizing the models for insight into the future system behavior on the other hand. While PRA models are typically used as a decision aid during the design phase of a space mission, we advocate adjusting them based on the observed behavior of the spacecraft and utilizing them for decision support during the operations phase.

dynamic fault trees↗

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗

Merging Intelligent Systems Technologies with CFD Analysis Strategies: Prototype Development at NASA Ames

A collaborative team of researchers from fields of Computational Fluid Dynamics (CFD), fluid physics, computer architectures, and computer science and knowledge engineering have begun work on a prototype system that addresses several of industry's concerns in using NASA-developed CFD codes as part of the design cycle. A major problem exists in the application of CFD technologies within the aeronautics design cycle due primarily to misunderstandings in the ranges of applicability of the various solver codes or turbulence models. Features that arise during the CFD solution process need to be discriminated and recognized as actual flow features with physical support in the geometry and flow conditions of the problem being solved, or as numerical or non-physical errors arising from mis-application of solver code and its parameters, gridding strategies, or discretization. interpolations. The fundamental concept is to develop an intelligent computational system that can accept the engineer's definition of the problem and construct an optimal CFD solution. To do this requires capturing both the knowledge of how to apply the various CFD tools and how to adapt the application of those tools to flow structures as they evolve during the flow simulation. Embedded within this adaptive system approach is the additional desire to automatically identify and quantify the quality of resolution of the pertinent flow structures, be they genuine or error-induced, and then to adjust the solution strategy accordingly. This paper discusses the status of that prototyping effort.

Thompson, David E.↗

Dynamics and Properties of Global Aerosol using MODIS, AERONET and GOCART Model

Recently produced daily Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol data for the whole year of 2001 are used to show the concentration and dynamics of aerosol over ocean and large parts of the continents. The data were validated against the Aerosol Robotic Network (AERONET) measurements over land and ocean in a special issue in GRL now in press. Monthly averages and a movie based on the daily data are produced and used to demonstrate the spatial and temporal evolution of aerosol. The MODIS wide spectral range is used to distinguish fine smoke and pollution aerosol from coarse dust and salt. The aerosol is observed above ocean and land. The movie produced from the MODIS data provides a new dimension to aerosol observations by showing the dynamics of the system. For example in February smoke and dust emitted from the Sahel and West Africa is shown to travel to the North-East Atlantic. In April heavy dust and pollution from East Asia is shown to travel to North America. In May-June pollution and dust play a dynamical dance in the Arabian Sea and Bay of Bengal. In Aug-September smoke from South Africa and South America is shown to pulsate in tandem and to periodically to be transported to the otherwise pristine Southern part of the Southern Hemisphere. The MODIS data are compared with the Georgia Tech/Goddard Global Ozone Chemistry Aerosol Radiation Transport (GOCART) model to test and adjust source and sink strengths in the model and to study the effect of clouds on the representation of the satellite data.

Kaufman, Yoram↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Fault Detection and Correction for the Solar Dynamics Observatory Attitude Control System

The Solar Dynamics Observatory is an Explorer-class mission that will launch in early 2009. The spacecraft will operate in a geosynchronous orbit, sending data 24 hours a day to a devoted ground station in White Sands, New Mexico. It will carry a suite of instruments designed to observe the Sun in multiple wavelengths at unprecedented resolution. The Atmospheric Imaging Assembly includes four telescopes with focal plane CCDs that can image the full solar disk in four different visible wavelengths. The Extreme-ultraviolet Variability Experiment will collect time-correlated data on the activity of the Sun's corona. The Helioseismic and Magnetic Imager will enable study of pressure waves moving through the body of the Sun. The attitude control system on Solar Dynamics Observatory is responsible for four main phases of activity. The physical safety of the spacecraft after separation must be guaranteed. Fine attitude determination and control must be sufficient for instrument calibration maneuvers. The mission science mode requires 2-arcsecond control according to error signals provided by guide telescopes on the Atmospheric Imaging Assembly, one of the three instruments to be carried. Lastly, accurate execution of linear and angular momentum changes to the spacecraft must be provided for momentum management and orbit maintenance. In th~sp aper, single-fault tolerant fault detection and correction of the Solar Dynamics Observatory attitude control system is described. The attitude control hardware suite for the mission is catalogued, with special attention to redundancy at the hardware level. Four reaction wheels are used where any three are satisfactory. Four pairs of redundant thrusters are employed for orbit change maneuvers and momentum management. Three two-axis gyroscopes provide full redundancy for rate sensing. A digital Sun sensor and two autonomous star trackers provide two-out-of-three redundancy for fine attitude determination. The use of software to maximize chances of recovery from any hardware or software fault is detailed. A generic fault detection and correction software structure is used, allowing additions, deletions, and adjustments to fault detection and correction rules. This software structure is fed by in-line fault tests that are also able to take appropriate actions to avoid corruption of the data stream.

Starin, Scott R.↗

The QBO in Two GISS Global Climate Models: 1. Generation of the QBO

The adjustment of parameterized gravity waves associated with model convection and finer vertical resolution has made possible the generation of the quasi-biennial oscillation (QBO) in two Goddard Institute for Space Studies (GISS) models, GISS Middle Atmosphere Global Climate Model III and a climate/middle atmosphere version of Model E2. Both extend from the surface to 0.002 hPa, with 2deg × 2.5deg resolution and 102 layers. Many realistic features of the QBO are simulated, including magnitude and variability of its period and amplitude. The period itself is affected by the magnitude of parameterized convective gravity wave momentum fluxes and interactive ozone (which also affects the QBO amplitude and variability), among other forcings. Although varying sea surface temperatures affect the parameterized momentum fluxes, neither aspect is responsible for the modeled variation in QBO period. Both the parameterized and resolved waves act to produce the respective easterly and westerly wind descent, although their effect is offset in altitude at each level. The modeled and observed QBO influences on tracers in the stratosphere, such as ozone, methane, and water vapor are also discussed. Due to the link between the gravity wave parameterization and the models' convection, and the dependence on the ozone field, the models may also be used to investigate how the QBO may vary with climate change.

middle atmosphere dynamics↗