Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 1,315 records · Page 73

High Degree and Order Gravity Fields of the Moon Derived from GRAIL Data

The Gravity Recovery and Interior Laboratory (GRAIL) spacecraft conducted the mapping of the gravity field of the Moon from March 1, 2012 to May 29, 2012. The twin spacecraft acquired highly precise K Band range-rate (KBRR) intersatellite ranging data and Deep Space Network (DSN) data during this prime mission phase from altitudes of 15 to 75 km above the lunar surface over three lunar months. We have processed these data using the NASA GSFC GEODYN orbit determination and geodetic parameter estimation program, and we have determined gravity fields up to degree and order 420 in spherical harmonics. The new gravity solutions show improved correlations with LOLA-derived topography to high degree and order and resolve many lunar features in the geopotential with a resolution of less than 30 km, including for example the central peak of the crater Tycho. We discuss the methodology used for the processing of the GRAIL data, the quality of the orbit determination on the GRAIL satellites and the derivation of the solutions, and their evaluation with independent data, including Lunar Prospector. We show that with these new GRAIL gravity solutions, we can now fit the low altitude, extended mission Lunar Prospector tracking data better than with any previous gravity model that included the LP data.

Lemoine, F. G.↗

Maneuver Design and Calibration for the Genesis Spacecraft

Genesis is the fifth mission selected as part of NASA's Discovery Program. The objective of Genesis is to collect solar wind samples for a period of approximately two years while in a halo orbit about the Earth-Sun L I point. At the end of this period, the samples are to be returned to a specific recovery point on the Earth for subsequent analysis. This goal has never been attempted before and presents a formidable challenge in terms of mission design and operations, particularly planning and execution of propulsive maneuvers. To achieve a level of cost-effectiveness consistent with a Discovery-class mission, the Genesis spacecraft design was adapted to the maximum extent possible from designs used on earlier missions, such as Mars Global Surveyor (MGS) and Stardust, another sample collection mission. The spacecraft design for Genesis is shown. Spin stabilization was chosen for attitude control, in lieu of three-axis stabilization, with neither reaction wheels nor accelerometers included. This precludes closed-loop control of propulsive maneuvers and implies that any attitude changes, including spin changes and precessions, will behave like translational propulsive maneuvers and affect the spacecraft trajectory. Moreover, to minimize contamination risk to the samples collected, all thrusters were placed on the side opposite the sample collection canister. The orientation and characteristics of thrusters are indicated. For large maneuvers (>2.5 m/s), two 5 lbf thrusters will be used for delta v, with precession to the burn attitude, followed by spin-up from 1.6 to 10 rpm before the burn and spin down to 1.6 rpm afterwards, then precession back to the original spin attitude. For small maneuvers (<2.5 m/s), no spin change is needed and four 0.2 lbf thrusters are used for Av. Single or double 360 deg. precession changes are required whenever the desired delta v falls inside the two-way turn circle (about 0.4 m/s) based on the mass properties, spin rate and lever arm lengths based on thruster locations. In such instances, delta v resulting from spacecraft precession cannot be used effectively as a component of the desired delta v, and must therefore be removed by precessing at least one complete revolution around the turn circle. To eliminate cross-track execution errors, a second revolution in the opposite direction would also be performed. This paper will address the design of propulsive maneuvers in light of the aforementioned challenges and other constraints. Maneuver design will be performed jointly by the Navigation Team at JPL and the Spacecraft Team at LMA, based on the process indicated . Typical maneuver timelines will be presented which address considerations introduced by attitude changes. These include nutation, which is introduced by precessing or spinning down and must be given sufficient time to damp out prior to execution of subsequent events, as well as sun and earth pointing constraints, which must be considered to ensure sufficient spacecraft power and to minimize telecommunications interruptions, respectively. The paper will include a description of how individual propulsive maneuvers are resolved into components to account for delta v from translational burns and spacecraft attitude changes required to carry out such maneuvers. Contributions to maneuver delta v arising from attitude changes, based on mass properties for the period just after launch, are indicated. Similar curves will be presented spanning all mission phases from launch through return. A set of closed-form equations for resolving maneuver components, base on a specific delta v required for correction or deterministic changes to the spacecraft trajectory will be presented, as well. In addition to nominal maneuvers, special calibration maneuvers are planned to improve open-loop modeling of maneuvers and to reduce execution errors. Uncalibrated execution errors are indicated. Such errors could be reduced by 50% or more over the course of the mission. Special calibrations are of particular importance for the return leg of the mission, since the sample canister must be returned to a specific location within the Utah Test and Training Range (UTTR) for mid-air retrieval. An entry angle tolerance of no less than +/- 0.08 deg. is required to achieve this objective. Biasing of the final return maneuvers coupled with a specific maneuver mode to use a series of well-characterized spin changes to effect these maneuvers is part of the current Genesis baseline mission plan. Another important objective of calibrations is to better characterize precession maneuvers. Such maneuvers are part of most propulsive maneuvers, but are also required periodically to maintain sun-pointing for power or daily during solar-wind pointing during collection periods. Although relatively small, such maneuvers will have a significant cumulative impact on orbit determination, particularly in the halo portion of the mission. The current mission design also calls for three stationkeeping maneuvers during each halo orbit of approximately six months duration. These stationkeeping maneuvers may be sufficiently small that single or double 360 deg. precession changes may be required. Because there are no accelerometers on board the spacecraft, calibration can only be performed with the aid of ground-based radiometric tracking. To establish a high degree of accuracy in characterizing the magnitude of burns, the spacecraft spin axis should be along the line of sight to the Earth, providing Doppler measurements with <1 mm/sec accuracy in S-Band. Emission constraints allow such alignment only during certain portions of the mission when the Earth-spacecraft-sun geometry is favorable. The impact of precessions, or burns at times when geometry is not favorable, can be assessed by reconstruction of the spacecraft trajectory using tracking arcs of several days before and after the event.

Williams, Kenneth E.↗

Models of Angular Momentum Input to a Circumterrestrial Swarm from Encounters with Heliocentric Planetesimals

Models of lunar origin in which the Moon accretes in orbit about the Earth from material approaching the Earth from heliocentric orbits must overcome a fundamental problem: the approach orbits of such material would be, in the simplest approximation, equally likely to be prograde or retrograde about the Earth, with the result that accretion of such material adds mass but not angular momentum to circumterrestrial satellites. Satellite orbits would then decay due to the resulting drag, ultimately impacting onto the Earth. One possibility for adding both material and angular momentum to Earth orbit is investigated: imbalance in the delivered angular momentum between pro and retrograde Earth passing orbits which arises from the three body dynamics of planetesimals approaching the Earth from heliocentric space. In order to study angular momentum delivery to circumterrestrial satellites, the near Earth velocities were numerically computed as a function of distance from the Earth for a large array of orbits systematically spanning heliocentric phase space.

Davis, D. R.↗

Robust Control Algorithm for a Two Cart System and an Inverted Pendulum

The Rectilinear Control System can be used to simulate a launch vehicle during liftoff. Several control schemes have been developed that can control different dynamic models of the rectilinear plant. A robust control algorithm was developed that can control a pendulum to maintain an inverted position. A fluid slosh tank will be attached to the pendulum in order to test robustness in the presence of unknown slosh characteristics. The rectilinear plant consists of a DC motor and three carts mounted in series. Each cart s weight can be adjusted with brass masses and the carts can be coupled with springs. The pendulum is mounted on the first cart and an adjustable air damper can be attached to the third cart if desired. Each cart and the pendulum have a quadrature encoder to determine position. Full state feedback was implemented in order to develop the control algorithm along with a state estimator to determine the velocity states of the system. A MATLAB program was used to convert the state space matrices from continuous time to discrete time. This program also used a desired phase margin and damping ratio to determine the feedback gain matrix that would be used in the LabVIEW program. This experiment will allow engineers to gain a better understanding of liquid propellant slosh dynamics, therefore enabling them to develop more robust control algorithms for launch vehicle systems

Wilson, Chris L.↗

Subgrid Combustion Modeling for the Next Generation National Combustion Code

In the first year of this research, a subgrid turbulent mixing and combustion methodology developed earlier at Georgia Tech has been provided to researchers at NASA/GRC for incorporation into the next generation National Combustion Code (called NCCLES hereafter). A key feature of this approach is that scalar mixing and combustion processes are simulated within the LES grid using a stochastic 1D model. The subgrid simulation approach recovers locally molecular diffusion and reaction kinetics exactly without requiring closure and thus, provides an attractive feature to simulate complex, highly turbulent reacting flows of interest. Data acquisition algorithms and statistical analysis strategies and routines to analyze NCCLES results have also been provided to NASA/GRC. The overall goal of this research is to systematically develop and implement LES capability into the current NCC. For this purpose, issues regarding initialization and running LES are also addressed in the collaborative effort. In parallel to this technology transfer effort (that is continuously on going), research has also been underway at Georgia Tech to enhance the LES capability to tackle more complex flows. In particular, subgrid scalar mixing and combustion method has been evaluated in three distinctly different flow field in order to demonstrate its generality: (a) Flame-Turbulence Interactions using premixed combustion, (b) Spatially evolving supersonic mixing layers, and (c) Temporal single and two-phase mixing layers. The configurations chosen are such that they can be implemented in NCCLES and used to evaluate the ability of the new code. Future development and validation will be in spray combustion in gas turbine engine and supersonic scalar mixing.

Menon, Suresh↗

Permeation of halide anions through phospholipid bilayers occurs by the solubility-diffusion mechanism

Two alternative mechanisms are frequently used to describe ionic permeation of lipid bilayers. In the first, ions partition into the hydrophobic phase and then diffuse across (the solubility-diffusion mechanism). The second mechanism assumes that ions traverse the bilayer through transient hydrophilic defects caused by thermal fluctuations (the pore mechanism). The theoretical predictions made by both models were tested for halide anions by measuring the permeability coefficients for chloride, bromide, and iodide as a function of bilayer thickness, ionic radius, and sign of charge. To vary the bilayer thickness systematically, liposomes were prepared from monounsaturated phosphatidylcholines (PC) with chain lengths between 16 and 24 carbon atoms. The fluorescent dye MQAE (N-(ethoxycarbonylmethyl)-6-methoxyquinolinium bromide) served as an indicator for halide concentration inside the liposomes and was used to follow the kinetics of halide flux across the bilayer membranes. The observed permeability coefficients ranged from 10(-9) to 10(-7) cm/s and increased as the bilayer thickness was reduced. Bromide was found to permeate approximately six times faster than chloride through bilayers of identical thickness, and iodide permeated three to four times faster than bromide. The dependence of the halide permeability coefficients on bilayer thickness and on ionic size were consistent with permeation of hydrated ions by a solubility-diffusion mechanism rather than through transient pores. Halide permeation therefore differs from that of a monovalent cation such as potassium, which has been accounted for by a combination of the two mechanisms depending on bilayer thickness.

NASA Discipline Exobiology↗

Extinction-to-Backscatter Ratios of Lofted Aerosol Layers Observed During the First Three Months of CALIPSO Measurements

Case studies from the first three months of the Cloud and Aerosol Lidar and Infrared Pathfinder Spaceborne Observations (CALIPSO) measurements of lofted aerosol layers are analyzed using transmittance [Young, 1995] and two-wavelength algorithms [Vaughan et al., 2004] to determine the aerosol extinction-to-backscatter ratios at 532 and 1064 nm. The transmittance method requires clear air below the layer so that the transmittance through the layer can be determined. Suitable scenes are selected from the browse images and clear air below features is identified by low 532 nm backscatter signal and confirmed by low depolarization and color ratios. The transmittance and two-wavelength techniques are applied to a number of lofted layers and the extinction-to-backscatter ratios are compared with values obtained from the CALIPSO aerosol models [Omar et al., 2004]. The results obtained from these studies are used to adjust the aerosol models and develop observations based extinction-to-backscatter ratio look-up tables and phase functions. Values obtained by these techniques are compared to Sa determinations using other independent methods with a goal of developing probability distribution functions of aerosol type-specific extinction to backscatter ratios. In particular, the results are compared to values determined directly by the High Spectral Resolution Lidar (HSRL) during the CALIPSO CloudSat Validation Experiments (CCVEX) and Sa determined by the application of the two-wavelength lidar Constrained Ratio Aerosol Model-fit (CRAM) retrieval approach [Cattrall et al., 2005; Reagan et al., 2004] to the HSRL data. The results are also compared to values derived using the empirical relationship between the multiple-scattering fraction and the linear depolarization ratio by using Monte Carlo simulations of water clouds [Hu et al., 2006].

Omar, Ali H.↗

Using Decision Trees to Detect and Isolate Simulated Leaks in the J-2X Rocket Engine

The goal of this work was to use data-driven methods to automatically detect and isolate faults in the J-2X rocket engine. It was decided to use decision trees, since they tend to be easier to interpret than other data-driven methods. The decision tree algorithm automatically "learns" a decision tree by performing a search through the space of possible decision trees to find one that fits the training data. The particular decision tree algorithm used is known as C4.5. Simulated J-2X data from a high-fidelity simulator developed at Pratt & Whitney Rocketdyne and known as the Detailed Real-Time Model (DRTM) was used to "train" and test the decision tree. Fifty-six DRTM simulations were performed for this purpose, with different leak sizes, different leak locations, and different times of leak onset. To make the simulations as realistic as possible, they included simulated sensor noise, and included a gradual degradation in both fuel and oxidizer turbine efficiency. A decision tree was trained using 11 of these simulations, and tested using the remaining 45 simulations. In the training phase, the C4.5 algorithm was provided with labeled examples of data from nominal operation and data including leaks in each leak location. From the data, it "learned" a decision tree that can classify unseen data as having no leak or having a leak in one of the five leak locations. In the test phase, the decision tree produced very low false alarm rates and low missed detection rates on the unseen data. It had very good fault isolation rates for three of the five simulated leak locations, but it tended to confuse the remaining two locations, perhaps because a large leak at one of these two locations can look very similar to a small leak at the other location.

Schwabacher, Mark A.↗

A Model for the Breakup of Comet Linear (C/1999 S4)

We propose a mechanism based on the rubble-pile hypothesis of the cometary nucleus (Weissman 1986) to explain the catastrophic breakup of comet LINEAR (C/1999 S4) observed during July-August 2000. We suggest that a solid nucleus made up of 10-100 m "cometesimals" (Weidenschilling 1997) contains a network of inter-connected voids in the inter-cometesimal regions. The production of super-volatile (i.e., species more volatile than water) gases into these voids occurs due to the thermal wave propagating through the nucleus and associated phase transitions of water ice. The network of voids provides an efficient pathway for rapid propagation of these gases within the nucleus resulting in gas pressure caused stresses over a wide regime of the nucleus. This provides a mechanism for catastrophic breakups of small cometary nuclei such as comet LINEAR (C/1999 S4) as well as for some observed cometary outbursts including those that occur at large heliocentric distances (e.g., West et al. 1991). We emphasize the importance of techniques such as radar reflection tomography and radiowave transmission tomography (e.g., Kofman et al. 1998) aboard cometary missions to determine the three dimensional structure of the nucleus in particular the extent of large scale voids.

Samarasinha, Nalin H.↗

Early differentiation of the Moon: Experimental and modeling studies and experimental and modeling studies of massif anorthosites

NASA grant NAG9-329 was in effect from 3/1/89 to 8/31/94, the last 18 months being a no-cost extension. While the grant was in effect, the P.I., coworkers, and students gave 22 talks and poster sessions at professional meetings, published 12 articles in referred journals (one more is in press, and another is in review), and edited 2 workshop reports relevant to this project. Copies of all the publications are appended to this report. The major accomplishments during the grant period have derived from three quarters: 1) the application of quantitative models of fractional crystallization and partial melting to various problems in planetary science, such as the petrogenesis of picritic glasses and mare basalts and the implications of the SNC meteorites for martian evolution; 2) an experimental study of silicate liquid immiscibility relevant to early lunar differentiation and the petrogenesis of evolved highlands rocks; and 3) experimental studies of massif anorthosites and related rocks that provide terrestrial analogs for the proposed origin of lunar anorthosites by multistage processes. The low-pressure aspects of the quantitative models were developed by the P.I. in the 1980s with NASA support and culminated with a paper comparing the crystallization of terrestrial and lunar lavas. The basis for the high-pressure modifications to the quantitative models is a data set gleaned from high pressure melting experiments done at Lamont and is supplemented by published data from other labs that constrain the baric and compositional dependences of various liquidus phase boundaries such as olivine/orthopyroxene, relevant to the melting of the mantles of the terrestrial planets. With these models it is possible to predict not only the thermal and compositional evolution of magmatic liquids ranging in composition from lumar mare basalt to terrestrial calc-alkaline basalts, but also the small increments of fractional melting that are produced when mantle rises adiabatically. Copies of the crystallization/melting programs have been given to several colleagues in planetary science. Additionally, a series of computer graphics programs, based on the algorithms in the crystallization programs have been developed that display liquidus diagrams appropriate to input compositions.

Longhi, John↗

RXTE Monitoring of the Anomalous X-ray Pulsar 1E 1048.1-5937: Long-Term Variability and the 2007 March Event

After three years of no unusual activity, Anomalous X-ray Pulsar 1E 1048.1-5937 reactivated in 2007 March. We report on the detection of a large glitch (deltav/v = 1.63(2) x 10(exp -5)) on 2007 March 26 (MJD 54185.9), contemporaneous with the onset of a pulsed-flux flare, the third flare observed from this source in 10 years of monitoring with the Rossi X-ray Timing Explorer. Additionally, we report on a detailed study of the evolution of the timing properties, the pulsed flux, and the pulse profile of this source as measured by RXTE from 1996 July to 2008 January. In our timing study, we attempted phase coherent timing of all available observations. We show that in 2001, a timing anomaly of uncertain nature occurred near the rise of the first pulsed flux flare; we show that a likely glitch (deltav/v = 2.91(9) x 10(exp -6)) occurred in 2002, near the rise of the second flare, and we present a detailed description of the variations in the spin-down. In our pulsed flux study, we compare the decays of the three flares and discuss changes in the hardness ratio. In our pulse profile study, we show that the profile exhibited large variations near the peak of the first two flares, and several small short-term profile variations during the most recent flare. Finally, we report on the discovery of a small burst 27 days after the peak of the last flare, the fourth burst discovered from this source. We discuss the relationships between the observed properties in the framework of the magnetar model.

Dib, Rim↗

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

Development of Numerical Extended Hydrodynamics for Transition-Regime Non-Equilibrium Flows Encountered in Semiconductor Manufacturing Processes

Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Groth, Clinton P. T.↗

Main Power Distribution Unit for the Jupiter Icy Moons Orbiter (JIMO)

Around the year 2011, the Jupiter Icy Moons Orbiter (JIMO) will be launched and on its way to orbit three of Jupiter s planet-sized moons. The mission goals for the JIMO project revolve heavily around gathering scientific data concerning ingredients we, as humans, consider essential: water, energy and necessary chemical elements. The JIM0 is an ambitious mission which will implore propulsion from an ION thruster powered by a nuclear fission reactor. Glenn Research Center is responsible for the development of the dynamic power conversion, power management and distribution, heat rejection and ION thrusters. The first test phase for the JIM0 program concerns the High Power AC Power Management and Distribution (PMAD) Test Bed. The goal of this testing is to support electrical performance verification of the power systems. The test bed will incorporate a 2kW Brayton Rotating Unit (BRU) to simulate the nuclear reactor as well as two ION thrusters. The first module of the PMAD Test Bed to be designed is the Main Power Distribution Unit (MPDU) which relays the power input to the various propulsion systems and scientific instruments. The MPDU involves circuitry design as well as mechanical design to determine the placement of the components. The MPDU consists of fourteen relays of four different variations used to convert the input power into the appropriate power output. The three phase system uses 400 Vo1ts(sub L-L) rms at 1000 Hertz. The power is relayed through the circuit and distributed to the scientific instruments, the ION thrusters and other controlled systems. The mechanical design requires the components to be positioned for easy electrical wiring as well as allowing adequate room for the main buss bars, individual circuit boards connected to each component and power supplies. To accomplish creating a suitable design, AutoCAD was used as a drafting tool. By showing a visual layout of the components, it is easy to see where there is extra room or where the components may interfere with one another. By working with the electrical engineer who is designing the circuit, the specific design requirements for the MPDU were determined and used as guidelines. Space is limited due to the size of the mounting plate therefore each component must be strategically placed. Since the MPDU is being designed to fit into a simulated model of the spacecraft systems on the JIMO, components must be positioned where they are easily accessible to be wired to the other onboard systems. Mechanical and electrical requirements provided equally important limits which are combined to produce the best possible design of the MPDU.

Papa, Melissa R.↗

Impact of Growth Phase, Pigment Adaptation, and Climate Change Conditions on the Cellular Pigment and Carbon Content of Fifty‐One Phytoplankton Isolates

Owing to their importance in aquatic ecosystems, the demand for models that estimate phytoplankton biomass and community composition in the global ocean has increased over the last decade. Moreover, the impacts of climate change, including elevated carbon dioxide (CO2), increased stratification, and warmer sea surface temperatures, will likely shape phytoplankton community composition in the global ocean. Chemotaxonomic methods are useful for modeling phytoplankton community composition from marker pigments normalized to chlorophyll a (Chl a). However, photosynthetic pigments, particularly Chl a, are sensitive to nutrient and light conditions. Cellular carbon is less sensitive, so using carbon biomass instead may provide an alternative approach. To this end, cellular pigment and carbon concentrations were measured in 51 strains of globally relevant, cultured phytoplankton. Pigment-to-Chl a and pigment-to-carbon ratios were computed for each strain. For 25 strains, measurements were taken during two growth phases. While some differences between growth phases were observed, they did not exceed within class differences. Multiple strains of Amphidinium carterae, Ditylum brightwellii and Heterosigma akashiwo were measured to determine whether time in culture influenced pigment and carbon composition. No appreciable trends in cellular pigment or carbon content were observed. Lastly, the potential impact of climate change conditions on the pigment ratios was assessed using a multistressor experiment that included increased mean light, temperature, and elevated pCO2 on three species: Thalassiosira oceanica, Ostreococcus lucimarinus, and Synechococcus. The largest differences were observed in the pigment-to-carbon ratios, while the marker pigments largely covaried with Chl a. The implications of these observations to chemotaxonomic applications are discussed.

Pigments↗

Data correlation and analysis of arc tunnel and wind tunnel tests of RSI joints and gaps, phase 2. Volume 1: Technical report

Heat transfer data measured in gaps representative of those being employed for joints in the space shuttle reusable surface insulation (RSI) thermal protection systems (TPS) were assimilated, analyzed, and correlated. Several types of gap were investigated with emphasis on simple butt joints. Gap widths ranged from 0.0 to 0.76 cm and depths ranged from 1 to 6 cm. Laminar, transitional, and turbulent boundary layer flows over the gap opening were investigated. The angle between gap axis and external flow was varied between 0 and pi/2 radians. The contoured cross section gap performed significantly better than all other wide gaps and slightly better than all other narrow gap geometries. Three dimensional heating variations were observed within gaps in the absence of external flow pressure gradients. Interactions between heating within gaps and heating of adjacent top tile surfaces were observed. Gaps aligned with the flow were observed to promote boundary layer transition. Heat transfer correlation equations were obtained for many of the tests. The TPS thickness requirements with and without gaps were computed for a current shuttle entry trajectory. Experimental data employed in the study are summarized. A description of each test facility, run schedule and test conditions, model descriptive information, and heat flux data are included.

Cristensen, H. E.↗

Computational Modeling of Semiconductor Dynamics at Femtosecond Time Scales

The main objective of the Joint-Research Interchange NCC2-5149 was to develop computer codes for accurate simulation of femtosecond pulse propagation in semiconductor lasers and semiconductor amplifiers [I]. The code should take into account all relevant processes such as the interband and intraband carrier relaxation mechanisms and the many-body effects arising from the Coulomb interaction among charge carriers [2]. This objective was fully accomplished. We made use of a previously developed algorithm developed at NASA Ames [3]-[5]. The new algorithm was tested on several problems of practical importance. One such problem was related to the amplification of femtosecond optical pulses in semiconductors. These results were presented in several international conferences over a period of three years. With the help of a postdoctoral fellow, we also investigated the origin of instabilities that can lead to the formation of femtosecond pulses in different kinds of lasers. We analyzed the occurrence of absolute instabilities in lasers that contain a dispersive host material with third-order nonlinearities. Starting from the Maxwell-Bloch equations, we derived general multimode equations to distinguish between convective and absolute instabilities. We find that both self-phase modulation and intensity-dependent absorption can dramatically affect the absolute stability of such lasers. In particular, the self-pulsing threshold (the so-called second laser threshold) can occur at few times the first laser threshold even in good-cavity lasers for which no self-pulsing occurs in the absence of intensity-dependent absorption. These results were presented in an international conference and published in the form of two papers.

Agrawal, Govind P.↗

Retrieval of Soil Moisture and Roughness from the Polarimetric Radar Response

The main objective of this investigation was the characterization of soil moisture using imaging radars. In order to accomplish this task, a number of intermediate steps had to be undertaken. In this proposal, the theoretical, numerical, and experimental aspects of electromagnetic scattering from natural surfaces was considered with emphasis on remote sensing of soil moisture. In the general case, the microwave backscatter from natural surfaces is mainly influenced by three major factors: (1) the roughness statistics of the soil surface, (2) soil moisture content, and (3) soil surface cover. First the scattering problem from bare-soil surfaces was considered and a hybrid model that relates the radar backscattering coefficient to soil moisture and surface roughness was developed. This model is based on extensive experimental measurements of the radar polarimetric backscatter response of bare soil surfaces at microwave frequencies over a wide range of moisture conditions and roughness scales in conjunction with existing theoretical surface scattering models in limiting cases (small perturbation, physical optics, and geometrical optics models). Also a simple inversion algorithm capable of providing accurate estimates of soil moisture content and surface rms height from single-frequency multi-polarization radar observations was developed. The accuracy of the model and its inversion algorithm is demonstrated using independent data sets. Next the hybrid model for bare-soil surfaces is made fully polarimetric by incorporating the parameters of the co- and cross-polarized phase difference into the model. Experimental data in conjunction with numerical simulations are used to relate the soil moisture content and surface roughness to the phase difference statistics. For this purpose, a novel numerical scattering simulation for inhomogeneous dielectric random surfaces was developed. Finally the scattering problem of short vegetation cover above a rough soil surface was considered. A general scattering model for grass-blades of arbitrary cross section was developed and incorporated in a first order random media model. The vegetation model and the bare-soil model are combined and the accuracy of the combined model is evaluated against experimental observations from a wheat field over the entire growing season. A complete set of ground-truth data and polarimetric backscatter data were collected. Also an inversion algorithm for estimating soil moisture and surface roughness from multi-polarized multi-frequency observations of vegetation-covered ground is developed.

Sarabandi, Kamal↗