Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “species dynamics”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

CFD analyses of combustor and nozzle flowfields

The objectives of the research are to improve design capabilities for low thrust rocket engines through understanding of the detailed mixing and combustion processes. A Computational Fluid Dynamic (CFD) technique is employed to model the flowfields within the combustor, nozzle, and near plume field. The computational modeling of the rocket engine flowfields requires the application of the complete Navier-Stokes equations, coupled with species diffusion equations. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an ongoing experimental program at NASA LeRC. The numerical procedure is performed on both time-marching and time-accurate algorithms, using an LU approximate factorization in time, flux split upwinding differencing in space. The integrity of fuel film cooling along the wall, its effectiveness in the mixing with the core flow including unsteady large scale effects, the resultant impact on performance and the assessment of the near plume flow expansion to finite pressure altitude chamber are addressed.

Tsuei, Hsin-Hua↗
NASA NTRS1995

Aerothermodynamic Facilities And Measurement: Flow Characterization in Shock Tunnels

This presentation will examine the key performance aspects of shock tunnels as they relate to their use as aerothermodynamic flow simulation facilities. Assessment of shock tube reservoir conditions and flow contaminants generated in the shock tube will be presented along with their limiting impact on viable test envelopes, Facility nozzle performance as it pertains to test time assessment and nozzle exit flow quality (survey of pressure, temperature, and species) will be addressed. Also included will be a discussion of free stream flow diagnostics, both intrusive and nonintrusive, for measurement of critical flow properties not directly inferred from surface mounted transducers. The use of computational fluid dynamics for purposes of validating experimental measurements as well as predicting performance in regimes where measurements are not feasible or possible will be discussed. The use of CFD for facility research and design will also be presented.

Cavolowsky, John A.↗
NASA NTRSDate not supplied

TPSAS-NF1676L-27509-DND

The information about composition of aerosols has a great importance for monitoring and understanding of climate and environment dynamics. Such information can be obtained using situ measurements or chemical transport models. However, in situ sampling has limited spatial and temporal coverage, while estimations have large uncertainties. The present work enables the monitoring of aerosol chemical species from space-borne observations, providing observationally-based results with spatial and temporal coverage. Following the ideas of Schuster et al.(2005, 2009, 2016), we retrieve chemical composition directly from satellite observations without intermediate retrieval of the refractive index (in contrast with Schuster's approach). This approach is expected to reduce the influence of modeling uncertainties, and to provide additional constraints in situations where remote sensing observations do not have enough spectral sensitivity to refractive index. One of principal difficulties is the identification of an adequate model for linking refractive index to chemical composition. Therefore, the initial effort of this work has focused on identifying an optimal physico-chemistry to refractive index conversion model. With that purpose, we first tested the retrieval approach using a simplified volume-weighting model and updated by the Maxwell Garnett mixing model. This concept was incorporated into the GRASP algorithm designed to retrieve an extended set of atmospheric parameters from remote sensing observations. Then a series of sensitivity tests using synthetic data of POLDER/PARASOL were conducted, and followed by inversion of real PARASOL observations over Africa. The sensitivity tests showed that these two models allow the retrieval to distinguish amongst the assumed chemical species. Results obtained from real PARASOL data demonstrated good agreement with the optical characteristics provided by AERONET. In addition, the obtained patterns of chemical component distribution agreed with known physical expectations. We also compared the chemical composition retrievals to GOCART model simulations and obtained good consistencies.

Lei Li↗
DOE OSTI2024

Monitoring the Reaction Dynamics of UF6 by Cryogenic Layering and FTIR Spectroscopy

Uranium hexafluoride (UF6) is a commonly used material feedstock for uranium enrichment processes. When introduced to water in the atmosphere, it reacts rapidly to form uranyl fluoride (UO2F2). Here, we investigate the UF6 hydrolysis reaction by cryogenically trapping reaction intermediates and characterizing the trapped species by FTIR. The reactant species are sequentially layered onto a diamond substrate held at 10K by a closed cycle liquid helium cryostat. At this temperature, the hydrolysis reaction is not spontaneous and can be catalyzed by the introduction of heat. Upon heating, the reaction moves through several intermediate compounds before proceeding to the final UO2F2 product. Several previously unobserved bands appear while the reaction progresses which may help to elucidate the mechanism behind UF6 hydrolysis.

McNamara, III, Louis E.↗
DOE OSTI2024

Chloroplast Stress Signals: Control of Retrograde Signaling, Chloroplast Turn-Over, and Cell Fate Decisions

Chloroplasts (photosynthetic plastids) are semiautonomous organelles that contain their own small genomes. The proteomes of chloroplasts, however, are a mixture of plastid and nuclear-encoded proteins. Chloroplasts perform photosynthesis, which is prone to damaging the organelles, leading to the production of reactive oxygen species (ROS) that damage the cell under environmental stresses. Thus, for the cell to maintain proper chloroplast function, efficient photosynthesis, and avoid ROS damage, it relies on complex crosstalk between the chloroplast, the nucleus, other organelles within the cell, and the cytoplasm in between. This communication involves retrograde signals from chloroplasts to control nuclear gene expression, programmed cell death (PCD), and chloroplast degradation. Here we review these signals with an emphasis on the roles of the ROS singlet oxygen ( 1 O 2 ) and plastid gene expression. We cover (1) recent work on understanding how multiple 1 O 2 signaling pathways can be initiated within stressed chloroplasts, (2) how individualized post-translational regulatory systems allow chloroplasts to control their proteomes and degradation, and (3) how chloroplast signals ultimately control cell fate decisions, such as PCD, senescence, and vacuole-mediated degradation of chloroplasts (chloroplast quality control). Overall, this chapter discusses how chloroplasts can act as environmental sensors for the cell and allow plants to acclimate to stress and thrive in dynamic environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Designed 2D protein crystals as dynamic molecular gatekeepers for a solid-state device

The sensitivity and responsiveness of living cells to environmental changes are enabled by dynamic protein structures, inspiring efforts to construct artificial supramolecular protein assemblies. However, despite their sophisticated structures, designed protein assemblies have yet to be incorporated into macroscale devices for real-life applications. We report a 2D crystalline protein assembly of C98/E57/E66 L-rhamnulose-1-phosphate aldolase ( CEE RhuA) that selectively blocks or passes molecular species when exposed to a chemical trigger. CEE RhuA crystals are engineered via cobalt(II) coordination bonds to undergo a coherent conformational change from a closed state (pore dimensions <1 nm) to an ajar state (pore dimensions ~4 nm) when exposed to an HCN(g) trigger. When layered onto a mesoporous silicon (pSi) photonic crystal optical sensor configured to detect HCN (g) , the 2D CEE RhuA crystal layer effectively blocks interferents that would otherwise result in a false positive signal. The 2D CEE RhuA crystal layer opens in selective response to low-ppm levels of HCN (g) , allowing analyte penetration into the pSi sensor layer for detection. These findings illustrate that designed protein assemblies can function as dynamic components of solid-state devices in non-aqueous environments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗
DOE OSTI2023

The impact of physiologically relevant temperatures on physical properties of thylakoid membranes: a molecular dynamics study

Thylakoid membranes are energy-converting membranes with a unique lipid composition. Though the membranes are primarily composed of proteins, their photosynthetic function is strongly influenced by the lipid constituents. Here we characterize, with molecular dynamics (MD) simulations, lipid bilayers with compositions representative of plant thylakoid membranes. We determine, in a wide range of temperatures, the physical parameters of the model membranes which are relevant for the photosynthetic function. We found a marked impact of temperature on membrane permeability due to a combination of increased compressibility and curvature of the membrane at elevated temperatures. With increasing temperatures, we observed increasingly smeared transmembrane density profiles of the membrane forming lipid headgroups predicting increased membrane flexibility. The diffusion coefficient of the lipids increased with temperature without apparent specificity for lipid species. Instead of a comprehensive experimental dataset in the relevant temperature range, we quantitatively compared and validated our MD results with MD simulations on a dipalmitoylphosphatidylcholine model system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Arctic Shrub Expansion, Plant Functional Trait Variation, and Effects on Belowground Carbon Cycling (Final Technical Report)

Terrestrial ecosystems are undergoing dramatic changes in response to climate warming, and these changes are expected to feedback to the atmosphere, potentially altering the trajectory of future climate change. Feedbacks from Arctic ecosystems are a major concern because the Arctic is projected to warm significantly in the 21 st century and because >50% of global belowground organic carbon is stored in permafrost and overlying soils. Warming-driven release of this carbon could drastically increase atmospheric greenhouse gas concentrations and accelerate climate warming. Plant communities are also responding to warming, as evidenced by the widely documented increase in woody-shrub growth and “greening” across much of the Arctic tundra biome. This vegetation shift may offset or amplify warming by altering carbon cycling. The direction and magnitude of shrub effects remain highly uncertain, however, due to limited understanding of the consequences of shrub expansion for belowground carbon cycling and simplification of these relationships in models. The major shrubs expanding in the Arctic (Betula, Salix, and Alnus) vary widely with respect to aboveground and belowground traits (e.g., tissue production and chemistry, rooting depth, microbial symbionts), and may also exhibit substantial intraspecific variation in these traits in response to environmental conditions. Such variation is likely to have profound implications for soil carbon cycling. The overarching goal of this project was to improve process-based understanding of the influence of shrub expansion on carbon cycling to enable improved representation of carbon dynamics in ecosystem and Earth system models. We investigated how plant functional traits vary among shrub genera, respond to environmental conditions, and affect belowground carbon and nutrient cycling by quantifying relationships among functional traits and biogeochemical cycling along edaphic gradients nested within a climate gradient in the Alaskan tundra. We found consistent differences in leaf and root traits among shrub genera and between shrubs and a widespread sedge species, indicating diverse nutrient acquisition strategies and belowground impacts among different arctic shrubs. We also found striking differences in trait values among individuals within the same species or genera within sites. Soil parameters were more important than climate parameters for predicting size and leaf trait variation, and root trait responses were less dependent on climate overall. For all but one root trait, including parameters representing aboveground traits improved the predictive ability of models. These results demonstrate that tundra shrub traits vary considerably at local scales and soil factors drive this variation, especially belowground. Furthermore, leveraging information about aboveground traits and soil conditions can improve predictions of how belowground traits will respond to climate change. Despite these differences, soil carbon and nitrogen pools in the active layer did not vary among plots dominated by different shrub or sedge genera. Instead, pool sizes generally decreased from warmer to colder sites, consistent with a productivity gradient. Patterns of isotopic N composition indicate that shrubs tighten nitrogen cycling via nitrogen resorption or immobilization of shrub litter. Overall, these results suggest that further identifying the specific shrub genera in the tundra landscape will ultimately provide better predictions of belowground dynamics across the changing arctic. We also performed simulation experiments with the Terrestrial Ecosystem Model (TEM) incorporated in the Predictive Ecosystem Analyzer (PEcAn) framework, treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We performed simulations across different types of tundra, including shrub tundra. One key finding was that both model sensitivity and uncertainty to a given parameter could vary within the same type of tundra, but in a different geographical location, such as over the climate gradient of shrub tundra described above. We organized a special session at the annual meeting of the Ecological Society of America in August 2019 to disseminate our results, refine recommendations for model improvement, and initiate collaborations to implement these recommendations in existing models of tundra carbon dynamics at ecosystem to Earth system scales. Our results support DOE near-term priorities by providing mechanistic insights into the role of vegetation change in the terrestrial carbon cycle in a region that is inadequately represented in Earth system models. Current models reduce the complexity of Arctic vegetation to a small number of plant functional types (PFTs). This approach implicitly assumes that each PFT represents the average ecological function of its constituent species, thus ignoring the effects of trait variation on biogeochemical cycling and potentially leading to large uncertainty in the sign and magnitude of ecosystem feedbacks to climate. By quantifying variation of plant functional traits across broad gradients of climatic and edaphic conditions and elucidating the linkages of such variation with carbon and nutrient cycling, our results illustrate the need and create a foundation for further developing trait-based modeling approaches that allow the traits of PFTs to vary as a function of environmental conditions. These approaches should improve the capacity of simulation models to offer insights into ecosystem carbon dynamics associated with novel plant communities in a rapidly changing Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The effects of allelochemical transfer on the dynamics of hosts, parasitoids, and competing hyperparasitoids

Allelochemicals produced by plants may be ingested by herbivorous insects and transferred to higher trophic levels with potentially deleterious effects. We develop a system of differential equations to investigate the effect of the transfer of allelochemicals, such as nicotine, on the population dynamics of a system of hosts, parasitoids, and two competing hyperparasitoids that attack different life stages of the parasitoids. We find both somewhat deleterious effects of nicotine on the larvae-attacking hyperparasitoids and increased attack rates for the pupae-attacking hyperparasitoids can promote coexistence. We also use an evolutionary game-theoretic approach to determine the optimal distribution of hyperparasitoid attacks among nicotine-producing and nicotine-free plants. With strong deleterious effects of nicotine and increased attack rates for the pupae-attacking hyperparasitoid, we find both species attack parasitoids on the nicotine-free plant but only pupae-attacking hyperparasitoids attack parasitoids on the nicotine-producing plant.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Oxygen Electrocatalysis on Mixed-Metal Oxides/Oxyhydroxides: From Fundamentals to Membrane Electrolyzer Technology

Catalyzing the oxygen evolution reaction (OER) is important for key energy-storage technologies, particularly water electrolysis and photoelectrolysis for hydrogen fuel production. Under neutral-to-alkaline conditions, first-row transition-metal oxides/(oxy)hydroxides are the fastest-known OER catalysts and have been the subject of intense study for the past decade. Critical to their high performance is the intentional or accidental addition of Fe to Ni/Co oxides that convert to layered (oxy)hydroxide structures during the OER. Unraveling the role that Fe plays in the catalysis and the molecular identity of the true “active site” has proved challenging, however, due to the dynamics of the host structure and absorbed Fe sites as well as the diversity of local structures in these disordered active phases. In this Account, we highlight our work to understand the role of Fe in Ni/Co (oxy)hydroxide OER catalysts. We first discuss how we characterize the intrinsic activity of the first-row transition-metal (oxy)hydroxide catalysts as thin films by accounting for the contributions of the catalyst-layer thickness (mass loading) and electrical conductivity as well as the underlying substrate’s chemical interactions with the catalyst and the presence of Fe species in the electrolyte. We show how Fe-doped Ni/Co (oxy)hydroxides restructure during catalysis, absorb/desorb Fe, and in some cases degrade or regenerate their activity during electrochemical testing. We highlight the relevant techniques and procedures that allowed us to better understand the role of Fe in activating other first-row transition metals for OER. We find several modes of Fe incorporation in Ni/Co (oxy)hydroxides and show how those modes correlate with activity and durability. We also discuss how this understanding informs the incorporation of earth-abundant transition-metal OER catalysts in anion-exchange-membrane water electrolyzers (AEMWE) that provide a locally basic anode environment but run on pure water and have advantages over the more-developed proton-exchange-membrane water electrolyzers (PEMWE) that use platinum-group-metal (PGM) catalysts. We outline the key issues of introducing Fe-doped Ni/Co (oxy)hydroxide catalysts at the anode of the AEMWE, such as the oxidative processes triggered by Fe species traveling through the polymer membrane, pH-gradient effects on the catalyst stability, and possibly limited catalyst utilization in the compressed stack configuration. We also suggest possible mitigation strategies for these issues. Lastly, we summarize remaining challenges including the long-term stability of Fe-doped Ni/Co (oxy)hydroxides under OER conditions and the lack of accurate models of the dynamic active surface that hinder our understanding of, and thus ability to design, these catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Observation of Slow Eigen-Zundel Interconversion in H + (H 2 O) 6 Clusters upon Isomer-Selective Vibrational Excitation and Buffer Gas Cooling in a Cryogenic Ion Trap

The formation of isomers when trapping floppy cluster ions in a temperature-controlled ion trap is a generally observed phenomenon. This involves collisional quenching of the ions initially formed at high temperature by buffer gas cooling until their internal energies fall below the barriers in the potential energy surface that separate them. Here we explore the kinetics at play in the case of the two isomers adopted by the H + (H 2 O) 6 cluster ion that differ in the proton accommodation motif. One of these is most like the Eigen cation with a tricoordinated hydronium motif (denoted E), and the other is most like the Zundel ion with the proton equally shared between two water molecules (denoted Z). After initial cooling to about 20 K in the radiofrequency (Paul) trap, the relative populations of these two spectroscopically distinct isomers are abruptly changed through isomer-selective photoexcitation of bands in the OH stretching region with a pulsed (~6 ns) infrared laser while the ions are in the trap. We then monitor the relaxation of the vibrationally excited clusters and reformation of the two cold isomers by recording infrared photodissociation spectra with a second IR laser as a function of delay time from the initial excitation. The latter spectra are obtained after ejecting the trapped ions into a time-of-flight photofragmentation mass spectrometer, thus enabling long (~0.1 s) delay times. Excitation of the Z isomer is observed to display long-lived vibrationally excited states that are collisionally cooled on a ms time scale, some of which quench into the E isomer. These excited E species then display spontaneous interconversion to the Z form on a ~10 ms time scale. Furthermore, these qualitative observations set the stage for a series of experimental measurements that can provide quantitative benchmarks for theoretical simulations of cluster dynamics and the potential energy surfaces that underlie them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Optical O 2 Sensors Also Respond to Redox Active Molecules Commonly Secreted by Bacteria

From a metabolic perspective, molecular oxygen (O 2 ) is arguably the most significant constituent of Earth’s atmosphere. Nearly every facet of microbial physiology is sensitive to the presence and concentration of O 2 , which is the most favorable terminal electron acceptor used by organisms and also a dangerously reactive oxidant. As O 2 has such sweeping implications for physiology, researchers have developed diverse approaches to measure O 2 concentrations in natural and laboratory settings. Recent improvements to phosphorescent O 2 sensors piqued our interest due to the promise of optical measurement of spatiotemporal O 2 dynamics. However, we found that our preferred bacterial model, Pseudomonas aeruginosa PA14, secretes more than one molecule that quenches such sen sors, complicating O 2 measurements in PA14 cultures and biofilms. Assaying supernatants from cultures of 9 bacterial species demonstrated that this phenotype is common: all super natants quenched a soluble O 2 probe substantially. Phosphorescent O 2 probes are often embedded in solid support for protection, but an embedded probe called O 2 NS was quenched by most supernatants as well. Measurements using pure compounds indicated that quenching is due to interactions with redox-active small molecules, including phena zines and flavins. Uncharged and weakly polar molecules like pyocyanin were especially potent quenchers of O 2 NS. These findings underscore that optical O 2 measurements made in the presence of bacteria should be carefully controlled to ensure that O 2 , and not bacterial secretions, is measured, and motivate the design of custom O 2 probes for specific organisms to circumvent sensitivity to redox-active metabolites.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1979

A study of the feasibility of ultraviolet spectrometry for cometary missions

Ultraviolet spectra fo the comet West obtained by sounding rocket experiments in 1976 are reproduced and interpreted in order to estimate the expected brightness of the emission features and determine the spatial extent of these features for the proposed Halley Flyby/Tempel 2 rendezvous and the possible Halley or Encke flybys close to perihelion. A coma model was constructed and evaluated for the physical condition of candidate targets such as heliocentric distance, gas production, and composition. In addition to brightness profiles, the neutral and ion densities of the principal species are also dervied. The brightness profiles can be used to determine the feasibility of utilizing the space telescope to provide supporting observations during the mission. Basic parameters identified are spectral range, wavelength resolution, spatial resolution, sensitivity and dynamic range, rejection of scattered light, and integration or accumulation time.

Feldman, P. D.↗
NASA NTRS1992

Linkages between terrestrial ecosystems and the atmosphere

The primary research issue in understanding the role of terrestrial ecosystems in global change is analyzing the coupling between processes with vastly differing rates of change, from photosynthesis to community change. Representing this coupling in models is the central challenge to modeling the terrestrial biosphere as part of the earth system. Terrestrial ecosystems participate in climate and in the biogeochemical cycles on several temporal scales. Some of the carbon fixed by photosynthesis is incorporated into plant tissue and is delayed from returning to the atmosphere until it is oxidized by decomposition or fire. This slower (i.e., days to months) carbon loop through the terrestrial component of the carbon cycle, which is matched by cycles of nutrients required by plants and decomposers, affects the increasing trend in atmospheric CO2 concentration and imposes a seasonal cycle on that trend. Moreover, this cycle includes key controls over biogenic trace gas production. The structure of terrestrial ecosystems, which responds on even longer time scales (annual to century), is the integrated response to the biogeochemical and environmental constraints that develop over the intermediate time scale. The loop is closed back to the climate system since it is the structure of ecosystems, including species composition, that sets the terrestrial boundary condition in the climate system through modification of surface roughness, albedo, and, to a great extent, latent heat exchange. These separate temporal scales contain explicit feedback loops which may modify ecosystem dynamics and linkages between ecosystems and the atmosphere. The long-term change in climate, resulting from increased atmospheric concentrations of greenhouse gases (e.g., CO2, CH4, and nitrous oxide (N2O)) will further modify the global environment and potentially induce further ecosystem change. Modeling these interactions requires coupling successional models to biogeochemical models to physiological models that describe the exchange of water, energy, and biogenic trace gases between the vegetation and the atmosphere at fine time scales. There does not appear to be any obvious way to allow direct reciprocal coupling of atmospheric general circulation models (GCM's), which inherently run with fine time steps, to ecosystem or successional models, which have coarse temporal resolution, without the interposition of physiological canopy models. This is equally true for biogeochemical models of the exchange of carbon dioxide and trace gases. This coupling across time scales is nontrivial and sets the focus for the modeling strategy.

Bretherton, Francis↗
NASA NTRS2000

The Turbulent Life of Phytoplankton

Phytoplankton is a generic name for photosynthesizing microscopic organisms that inhabit the upper sunlit layer (euphotic zone) of almost all oceans and bodies of freshwater. They are agents for "primary production," the incorporation of carbon from the environment into living organisms, a process that, sustains the aquatic food web. It is estimated that phytoplankton contribute about half of the global primary production, the other half being due to terrestrial plants. By sustaining the aquatic food web and controlling the biogeochemical cycles through primary production, phytoplankton exert a dominant influence on life on earth. Turbulence influences this process in three very important ways. First, essential mineral nutrients are transported from the deeper layers to the euphotic zone through turbulence. Second, turbulence helps to suspend phytoplankton in the euphotic zone since in still water, the phytoplankton, especially the larger species, tend to settle out of the sunlit layers. Third, turbulence transports phytoplankton from the surface to the dark sterile waters, and this is an important mechanism of loss. Thus, stable phytoplankton populations are maintained through a delicate dynamic balance between the processes of turbulence, reproduction, and sinking. The first quantitative model for this was introduced by Riley, Stommel and Bumpus in 1949. This is an attempt to extend their efforts through a combination of analysis and computer simulation in order to better understand the principal qualitative aspects of the physical/biological coupling of this natural system.

S Ghosal↗