Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “species distribution 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 55 records · Page 3

An analysis of CCD images of the coma of Comet Halley

The analysis of selected CCD images of the coma of comet P/Halley is presented. The images were taken using specially designed filters that isolate regions of a comet's spectrum such that only sunlight which has been scattered by the dust in the coma is recorded. The modeling analysis objective is to make use of the skills developed in the development of Monte Carlo particle trajectory models for the distributions of gas species in cometary comae and to use those models as a basis for a new dust coma model. This model will include a self-consistant picture of the time-dependent dusty-gas dynamics of the inner coma and the three-dimensional time-dependent trajectories of the dust particles under the influence of solar gravity and solar radiation pressure in the outer coma. The model is intended to be used as a tool to analyze selected images from the two sets of CCD images with the hope that it will help the understanding of the effects of a number of important processes on the spatial morphology of the observed dust coma. The processes of importance to the observed dust coma include: (1) the dust particle size distribution function; (2) the terminal velocities of various sized dust particles in the inner coma; (3) the radiation scattering properties of dust particles, which are important both in terms of the observe scattered radiation and the radiation pressure acceleration on dust particles; (4) the fragmentation and/or vaporization of dust particles; and (5) the relative importance of CHON and silicate dust particles as they contribute both to the dusty-gasdynamics in the inner coma (that produce the dust particle terminal velocities) and to the observed spatial morphology on the outer dust coma.

Combi, Michael↗

The OH distribution in cometary atmospheres - A collisional Monte Carlo model for heavy species

The study presents an extension of the cometary atmosphere Monte Carlo particle trajectory model formalism which makes it both physically correct for heavy species and yet computationally reasonable. The derivation accounts for the collision path and scattering redirection of a heavy radical traveling through a fluid coma with a given radial distribution in outflow speed and temperature. The revised model verifies that the earlier fast-H atom approximations used in earlier work are valid, and it is applied to a case where the heavy radical formalism is necessary: the OH distribution. It is found that a steeper variation of water production rate with heliocentric distance is required for a water coma which is consistent with the velocity-resolved observations of Comet P/Halley.

Combi, Michael R.↗

Effects of atmospheric transport on column abundances of nitrogen and chlorine compounds in the arctic stratosphere

In the absence of heterogeneous processes, distributions of NO(y), HNO3, Cl(x) and HCl in the polar winter stratosphere should be determined by transport. A model simulating distributions of these species in the arctic for January and February 1989 is developed using observed fields of potential vorticity and potential temperature. Comparison of model results with column measurements from the DC-8 indicates conversion of NO(x) to HNO3, condensation of HNO3, conversion of HCl to CINO3, and conversion of HCl plus CINO3 to an unmeasured species. Heterogeneous processes strongly affect abundances of NO(y) and Cl(x) species in the winter arctic stratosphere.

Yatteau, John H.↗

Excitation of the ordinary electromagnetic mode in low-beta plasmas

The excitation of the ordinary mode propagating perpendicular to an external magnetic field in a plasma is discussed for the case where the ratio of the kinetic energy associated with electron motion parallel to the external magnetic field to the energy density of the magnetic field is low. The analysis is performed for plasma in which in addition to the thermal electrons there is a hot anisotropic electron component. Relativistic effects are taken into account in the stability analysis. Both a ring distribution and a loss cone distribution are used to model the hot electron species. For the ring distribution typical growth rates are larger than that for the case of the loss cone distribution. Excitation of the ordinary mode is found at numerous harmonics of the cyclotron frequency for the ring distribution, however, the loss cone distribution is found to give instability primarily at the fundamental harmonic.

Freund, H. P.↗

Gas Dynamics and Kinetics in the Cometary Coma: Theory and Observations

Our ability to describe the physical state of the expanding coma affects fundamental areas of cometary study both directly and indirectly. In order to convert measured abundances of gas species in the coma to gas production rates, models for the distribution and kinematics of gas species in the coma are required. Conversely, many different types of observations, together with laboratory data and theory, are still required to determine coma model attributes and parameters. Accurate relative and absolute gas production rates and their variations with time and from comet to comet are crucial to our basic understanding of the composition and structure of cometary nuclei and their place in the solar system. We review the gas dynamics and kinetics of cometary comae from both theoretical and observational perspectives, which are important for understanding the wide variety of physical conditions that are encountered.

Combi, Michael R.↗

Analysis of CCD images of the coma of comet P/Halley

The modeling analysis objective of this project is to make use of the skill acquired in the development of Monte Carlo particle trajectory models for the distributions of gas species in cometary comae as a basis for a new dust coma model. This model will include a self-consistent picture of the time-dependent dusty-gas dynamics of the inner coma and the three-dimensional time-dependent trajectories of the dust particles under the influence of solar gravity and solar radiation pressure in the outer coma. Our purpose is to use this model as a tool to analyze selected images from two sets of data of the comet P/Halley with the hope that we can help to understand the effects of a number of important processes on the spatial morphology of the observed dust coma. The study will proceed much in the same way as our study of the spatially extended hydrogen coma where we were able to understand the spatial morphology of the Lyman-alpha coma in terms of the partial thermalization of the hot H atoms produced by the photodissociation of cometary H2O and OH. The processes of importance to the observed dust coma include: (1) the dust particle size distribution function; (2) the terminal velocities of various sized dust particles in the inner coma; (3) the radiation scattering properties of dust particles, which are important both in terms of the observed scattered radiation and the radiation pressure acceleration on dust particles; (4) the fragmentation and/or vaporization of dust particles; (5) the relative importance of CHON and silicate dust particles as they contribute both to the dusty-gas dynamics in the inner coma (that produce the dust particle terminal velocities) and to the observed spatial morphology of the outer dust coma; and (6) the time and direction dependence of the source of dust.

Combi, Michael R.↗

Hyperspectral Mapping of the Invasive Species Pepperweed and the Development of a Habitat Suitability Model

Mapping and predicting the spatial distribution of invasive plant species is central to habitat management, however difficult to implement at landscape and regional scales. Remote sensing techniques can reduce the impact field campaigns have on these ecologically sensitive areas and can provide a regional and multi-temporal view of invasive species spread. Invasive perennial pepperweed (Lepidium latifolium) is now widespread in fragmented estuaries of the South San Francisco Bay, and is shown to degrade native vegetation in estuaries and adjacent habitats, thereby reducing forage and shelter for wildlife. The purpose of this study is to map the present distribution of pepperweed in estuarine areas of the South San Francisco Bay Salt Pond Restoration Project (Alviso, CA), and create a habitat suitability model to predict future spread. Pepperweed reflectance data were collected in-situ with a GER 1500 spectroradiometer along with 88 corresponding pepperweed presence and absence points used for building the statistical models. The spectral angle mapper (SAM) classification algorithm was used to distinguish the reflectance spectrum of pepperweed and map its distribution using an image from EO-1 Hyperion. To map pepperweed, we performed a supervised classification on an ASTER image with a resulting classification accuracy of 71.8%. We generated a weighted overlay analysis model within a geographic information system (GIS) framework to predict areas in the study site most susceptible to pepperweed colonization. Variables for the model included propensity for disturbance, status of pond restoration, proximity to water channels, and terrain curvature. A Generalized Additive Model (GAM) was also used to generate a probability map and investigate the statistical probability that each variable contributed to predict pepperweed spread. Results from the GAM revealed distance to channels, distance to ponds and curvature were statistically significant (p < 0.01) in determining the locations of suitable pepperweed habitats.

Pepperweed↗

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗

Effect of non-equilibrium flow chemistry on the heating distribution over the MESUR forebody during a Martian entry

Effect of flow field properties on the heating distribution over a 140 deg blunt cone was determined for a Martian atmosphere using Euler, Navier-Stokes (NS), viscous shock layer (VSL), and reacting boundary layer (BLIMPK) equations. The effect of gas kinetics on the flow field and the surface heating distribution were investigated. Gas models with nine species and nine reactions were implemented into the codes. Effects of surface catalysis on the heating distribution were studied using a surface kinetics model having five reactions.

Chen, Yih-Kang↗

Calcium in Mercury's Exosphere: Modeling MESSENGER Data

Mercury is surrounded by a surface-bounded exosphere comprised of atomic species including hydrogen, sodium, potassium, calcium, magnesium, and likely oxygen. Because it is collisionless. the exosphere's composition represents a balance of the active source and loss processes. The Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface. Space ENvironment. GEochemistry. and Ranging (MESSENGER) spacecraft has made high spatial-resolution observations of sodium, calcium, and magnesium near Mercury's surface and in the extended, anti-sunward direction. The most striking feature of these data has been the substantial differences in the spatial distribution of each species, Our modeling demonstrates that these differences cannot be due to post-ejection dynamics such as differences in photo-ionization rate and radiation pressure. but instead point to differences in the source mechanisms and regions on the surface from which each is ejected. The observations of calcium have revealed a strong dawn/dusk asymmetry. with the abundance over the dawn hemisphere significantly greater than over the dusk. To understand this asymmetry, we use a Monte Carlo model of Mercury's exosphere that we developed to track the motions of exospheric neutrals under the influence of gravity and radiation pressure. Ca atoms can be ejected directly from the surface or produced in a molecular exosphere (e.g., one consisting of CaO). Particles are removed from the system if they stick to the surface or escape from the model region of interest (within 15 Mercury radii). Photoionization reduces the final weighting given to each particle when simulating the Ca radiance. Preliminary results suggest a high temperature ( I-2x 10(exp 4) K) source of atomic Ca concentrated over the dawn hemisphere. The high temperature is consistent with the dissociation of CaO in a near-surface exosphere with scale height <= 100 km, which imparts 2 eV to the freshly produced Ca atom. This source region and energy are consistent with data from the three MESSENGER flybys; whether this holds true for the data obtained in orbit is under investigation.

Burger, Matthew H.↗

Temporal and spatial variations observed in the ionospheric composition of Venus - Implication for empirical modelling

The distributions of several ion species in a data base obtained by in situ measurements of the thermal ion composition of the ionosphere of Venus by the Pioneer Venus Orbiter have been sorted to identify temporal and spatial variations and determine the feasibility of an analytical representation of the experimental results. The first results from the sorting of several prominent ions including O(+), O2(+), and H(+) and several minor ions including CO2(+), C(+), and H2(+) reveal significant diurnal variations which consist of strong day to night contrast in the ion concentrations, with differences of one to two orders of magnitude, depending upon ion mass and altitude. It is suggested that repeatable day to night gradients in the ion distribution are adaptable to parametric modeling.

Taylor, H. A., Jr.↗

Range wide genetic differentiation in the bull kelp Nereocystis luetkeana with a seascape genetic focus on the Salish Sea

Introduction: In temperate regions, one of the most critical determinants of present range-wide genetic diversity was the Pleistocene climate oscillations, the most recent one created by the last glacial maximum (LGM). This study aimed to describe N. luetkeana genetic structure across its entire range (Alaska to California) and test different models of population connectivity within the Salish Sea. This region was colonized after the LGM and has been under increased disturbance in recent decades. Methods: We utilized microsatellite markers to study N. luetkeana genetic diversity at 53 sites across its range. Using higher sampling density in the Salish Sea, we employed a seascape genetics approach and tested isolation by hydrodynamic transport and environment models. Results: At the species distribution scale, we found four main groups of genetic co-ancestry, Alaska; Washington with Vancouver Island’s outer coast and Juan de Fuca Strait; Washington’s inner Salish Sea; and Oregon with California. The highest allelic richness (AR) levels were found in California, near the trailing range edge, although AR was also high in Alaska. The inner Salish Sea region had the poorest diversity across the species distribution. Nevertheless, a pattern of isolation by hydrodynamic transport and environment was supported in this region. Discussion: The levels of allelic, private allele richness and genetic differentiation suggest that during the LGM, bull kelp had both northern and southern glacial refugia in the Prince of Wales Island-Haida Gwaii region and Central California, respectively. Genetic diversity in Northern California sites seems resilient to recent disturbances, whereas the low levels of genetic diversity in the inner Salish Sea are concerning.

59 BASIC BIOLOGICAL SCIENCES↗

Space-based Remote Sensing: A Tool for Studying Bird Migration Across Multiple Scales

The study of bird migration on a global scale is one of the compelling and challenging problems of modern biology with major implications for human health and conservation biology. Migration and conservation efforts cross national boundaries and are subject to numerous international agreements and treaties. Space based technology offers new opportunities to shed understanding on the distribution and migration of organisms on the planet and their sensitivity to human disturbances and environmental changes. Our working hypothesis is that individual organism biophysical models of energy and water balance, driven by satellite measurements of spatio-temporal gradients in climate and habitat, will help us to explain the variability in avian species richness and distribution. Further, these models provide an ecological forecasting tool for science and application users to visualize the possible consequences of loss of wetlands, flooding, or other natural disasters such as hurricanes on avian biodiversity and bird migration.

Smith, James A.↗

Imaging and spectroscopy of Comet P/Halley

The goals of this investigation are the analysis of a large set of high-resolution echelle/reticon spectra, and the reduction and analysis of a set of narrow-band-filtered charge coupled device (CCD) images of Comet Halley taken during the preperihelion period at Oak Ridge Observatory by Dr. R. E. McCrosky. The scientific objectives associated with these goals are the determination of the spatial distributions of several important radicals, atoms and ions in the coma. These include C2, CN, C3, H2O(+) and CO(+) from the image data. The analysis of the neutral species distribution with Monte Carlo models will aid in the understanding of their production and decay mechanisms as well as serve as an important indicator of the physical conditions in the inner coma. The spatial distributions of the ions will serve as a guide to constrain the complex model necessary for understanding the interaction of the solar wind and the cometary ions. Work during this past year has been devoted largely to the reduction of the standard star photometry for the CCD image data set, as well as the re-flat-fielding of a number of the comet images. We are pleased to report that despite a number of setbacks and the small effort devoted to this work (2 1/2 months for the PI and a generous share of completely unsupported time by Dr. McCrosky) that this portion of the work has been successfully completed. The goals for the upcoming final year of this project (under a new project number) are to complete the calibration of the CCD image data for inclusion in the IHW archive, to analyze a select portion of the neutral radical images with our Monte Carlo models, and to present the results of the 6300/region spectra as a guide to low-resolution spectral observers in order to yield the unambiguous separation of the contributions of cometary O(1D), airglow O(1D), and the numerous NH2 lines in that region of the spectrum.

Combi, Michael R.↗

New, Improved Goddard Bulk-Microphysical Schemes for Studying Precipitation Processes in WRF

An improved bulk microphysical parameterization is implemented into the Weather Research and Forecasting ()VRF) model. This bulk microphysical scheme has three different options, 2ICE (cloud ice & snow), 3ICE-graupel (cloud ice, snow & graupel) and 3ICE-hail (cloud ice, snow & hail). High-resolution model simulations are conducted to examine the impact of microphysical schemes on two different weather events (a midlatitude linear convective system and an Atlantic hurricane). The results suggest that microphysics has a major impact on the organization and precipitation processes associated with a summer midlatitude convective line system. The Goddard 3ICE scheme with a cloud ice-snow-hail configuration agreed better with observations in terms of rainfall intensity and a narrow convective line than did simulations with a cloud ice-snow-graupel or cloud ice-snow (i.e., 2ICE) configuration. This is because the 3ICE-hail scheme includes dense ice precipitating (hail) particle with very fast fall speed (over 10 in For an Atlantic hurricane case, the Goddard microphysical schemes had no significant impact on the track forecast but did affect the intensity slightly. The improved Goddard schemes are also compared with WRF's three other 3ICE bulk microphysical schemes: WSM6, Purdue-Lin and Thompson. For the summer midlatitude convective line system, all of the schemes resulted in simulated precipitation events that were elongated in the southwest-northeast direction in qualitative agreement with the observed feature. However, the Goddard 3ICE scheme with the hail option and the Thompson scheme agree better with observations in terms of rainfall intensity, expect that the Goddard scheme simulated more heavy rainfall (over 48 mm/h). For the Atlantic hurricane case, none of the schemes had a significant impact on the track forecast; however, the simulated intensity using the Purdue-Lin scheme was much stronger than the other schemes. The vertical distributions of model simulated cloud species (i.e., snow) are quite sensitive to microphysical schemes, which is an important issue for future verification against satellite retrievals. Both the Purdue-Lin and WSM6 schemes simulated very little snow compared to the other schemes for both the midlatitude convective line and hurricane cases. Sensitivity tests are performed for these two WRF schemes to identify that snow productions could be increased by increasing the snow intercept, turning off the auto-conversion from snow to graupel and reducing the transfer processes from cloud-sized particles to precipitation-sized ice.

Tao, Wei-Kuo↗

Chemistry and transport in a three-dimensional stratospheric model - Chlorine species during a simulated stratospheric warming

The distributions in the stratosphere of a variety of chemical species were calculated for a 6-day period during the February 1979 stratospheric major warming, using winds derived from a spectral forecast model which included O(x), NO(x), HO(x), and ClO(x) chemistries as well as longitudinally varying reaction rate coefficients and photolysis rates for these molecules. The results obtained indicate a particular importance of chemistry and transport for the Cl-containing species ClO, ClONO2, HCl, and HOCl. Dynamical effects dominate the variability of HCl, while diurnal effects dominate that of ClONO2 and ClO. The effects of strong planetary wave activity may be seen in terms of large longitudinal variability of the total HCl and ClONO columns in the stratosphere; in the middle and high northern latitudes, it is sufficiently large to exceed the diurnal variability of the column.

Kaye, Jack A.↗

Multiphysics Coupling Methods for Molten Salt Reactor Modeling and Simulation in VERA

New multiphysics coupling capabilities for molten salt reactor (MSR) analysis have been developed in the Virtual Environment for Reactor Applications (VERA). This development consisted of two main efforts. First, a generic species transport module was added in the CTF code, which is the thermal-hydraulics (TH) code for VERA. This module uses the velocity fields for which CTF solves during the TH calculation to transport species through the core and around the primary loop. Additionally, a gas sparging model has been added to CTF to model the movement of certain species, namely, fission products such as xenon gas, to transport between the molten salt and gas bubbles present in the salt. The second effort in this development was coupling this capability to VERA’s neutron transport code MPACT. This effort focused on coupling the detailed TH transport models in CTF to MPACT to account for feedback effects in the neutron transport calculations. Finally, the thermochemistry code Thermochimica has also been coupled to VERA. Thermochimica performs pointwise calculations for chemical potential and Gibbs free energy and determines what phases are produced by the temperature, pressure, and elemental concentrations at different locations in the primary loop. These capabilities are demonstrated using a model of the Molten Salt Reactor Experiment (MSRE). Furthermore, this reactor operated at Oak Ridge National Laboratory in the 1960s, providing sources of experimental data that were used to develop the model. Various combinations of species were modeled using VERA’s new multiphysics coupling capabilities. Species distributions and reactivity effects behaved as expected for the MSRE model, demonstrating that the coupling is behaving correctly and causing appropriate feedback. The results of these calculations show the potential for VERA to be used for a wide variety of MSR analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The stratopause semiannual oscillation in the NCAR Community Climate Model

The middle atmospheric version of the NCAR Community Climate Model (CCM2) has been used to study the development of the equatorial semiannual oscillation (SAO) in the stratosphere. The model domain extends from the ground to about 80 km, with a vertical resolution of 1 km. Transport of nitrous oxide (N2O) with simplified photochemistry is included in the calculation to illustrate the influence of tropical circulations on the distribution of trace species. Diagnosis of model output reveals two distinct phases in the evolution of the zonal mean state on the equator. In early December, a strong and broad easterly jet appears near the stratopause in connection with a midlatitude wave event (sudden stratospheric warming) that reverses the winter westerlies of the Northern Hemisphere throughout the upper stratosphere. When the wave forcing dies out, the radiative drive allows the westerlies to recover at midlatitudes, while easterlies persist in the tropics. The resulting strong meridional gradient of the zonal mean wind provides favorable conditions for the development of inertial instability at lower latitudes. The meridional circulation associated with the instability shapes the 'nose' of the easterly jet, reducing the extension of the unstable region. In equinoctial conditions, a jet of westerlies appears in the lower equatorial mesosphere and descends to lower altitudes; positive accelerations associated with the descending westerlies are due primarily to Kelvin waves. The descent of the westerly jet does not reproduce well the observed behavior of the SAO westerly phase, either in amplitude or in the extent of downward propagation. As a consequence, the model does not simulate the 'double peak' observed in the tropical distribution of N2O. Comparison of wave amplitudes in the model with those derived from satellite observations shows that the calculated amplitudes are larger than observed in the upper stratosphere. It follows that inadequate Kelvin wave forcing is not the cause of the weak westerly phase in the model, and that some other mechanism must be responsible for the generation of the strong westerly phase observed.

Sassi, Fabrizio↗