Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated emission”

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 829 records · Page 46

Feasibility of Observing Gamma-Ray Polarization from Cygnus X-1 Using a CubeSat

Instruments flown on CubeSats are small. Meaningful applications of CubeSats in astronomical observations rely on the choice of a particular subject that is feasible for CubeSats. Here we report the result of a feasibility study for observing gamma-ray polarization from Cygnus X-1 using a small Compton polarimeter on board a 3U CubeSat. Silicon detectors and cerium bromide scintillators were employed in the instrument models that we discuss in this study. Through Monte Carlo simulations with a Geant4-based MEGAlib package, we found that, with a 10 Ms on-axis, zenith-direction observation in a low-inclination, low-altitude, Earth-orbit radiation background environment, the minimum detectable polarization degree can be down to about 10% in 160–250 keV, 20% in 250–400 keV, and 65% in 400–2000 keV. A 3U CubeSat dedicated to observing Cygnus X-1 can therefore yield useful information on the polarization state of gamma-ray emissions from the brightest persistent X-ray black hole binary in the sky.

79 ASTRONOMY AND ASTROPHYSICS↗

In silico evaluation of plant nitrification suppression effects on agroecosystem nitrogen loss

Abstract Nitrification regulates the potential for nitrogen (N) loss from ecosystems because it converts ammonium to nitrate, which is susceptible to leaching and gaseous emissions. Crops can suppress the microbes that perform nitrification by exuding nitrification‐inhibiting compounds from their roots and taking up available ammonium. However, the effect of nitrification suppression on agroecosystem N losses remains poorly characterized, and a lack of temporal synchrony between nitrification, N losses, and nitrification suppression by plants could limit the effect of nitrification suppression. We used the DayCent‐CABBI model to evaluate the effectiveness of the suppression of nitrification by sorghum to reduce N 2 O emissions and nitrate leaching in an energy sorghum/soybean rotation at the Energy Farm in Urbana‐Champaign, IL. We simulated nitrification suppression at the measured levels (MeasNS) and at the maximum measured level applied to the entire growing season (MaxNS), and we also explored ways to better utilize nitrification suppression by altering the timing of urea–ammonium–nitrate fertilizer applications. Model experiments showed that most nitrification occurred immediately after fertilizer was applied, whereas nitrification suppression began to ramp up more than a month after planting. On an annual basis, MeasNS experiments showed a 1%–2% reduction in annual N 2 O emissions relative to no nitrification suppression (NoNS), and MaxNS experiments showed a 4%–9% reduction in annual N 2 O emissions relative to NoNS. Both nitrification suppression levels showed <1% reduction in nitrate leaching. Altering the timing of fertilizer applications to better synchronize nitrification suppression with high soil ammonium levels had mixed effects on annual N 2 O emissions and nitrate leaching and sometimes resulted in increased N losses. The timing of simulated N 2 O emissions shifted with the timing of fertilization, and N 2 O emissions from denitrification increased when N 2 O emissions from nitrification decreased. Increasing N retention during the non‐growing season may be more effective than growing‐season nitrification suppression for reducing annual N losses in the rainfed Midwest.

60 APPLIED LIFE SCIENCES↗

Theoretical and experimental studies relevant to interpretation of auroral emissions

This report describes the accomplishments of a program designed to develop the tools necessary to interpret auroral emissions measured from a space-based platform. The research was divided into two major areas. The first area was a laboratory study designed to improve our understanding of the space vehicle external environment and how it will affect the space-based measurement of auroral emissions. Facilities have been setup and measurements taken to simulate the gas phase environment around a space vehicle; the radiation environment encountered by an orbiting vehicle that passes through the Earth's radiation belts; and the thermal environment of a vehicle in Earth orbit. The second major area of study was a modeling program to develop the capability of using auroral images at various wavelengths to infer the total energy influx and characteristic energy of the incident auroral particles. An ab initio auroral calculation has been added to the extant ionospheric/thermospheric global modeling capabilities within our group. Once the addition of the code was complete, the combined model was used to compare the relative intensities and behavior of various emission sources (dayglow, aurora, etc.). Attached papers included are: 'Laboratory Facility for Simulation of Vehicle-Environment Interactions'; 'Workshop on the Induced Environment of Space Station Freedom'; 'Radiation Damage Effects in Far Ultraviolet Filters and Substrates'; 'Radiation Damage Effects in Far Ultraviolet Filters, Thin Films, and Substrates'; 'Use of FUV Auroral Emissions as Diagnostic Indicators'; and 'Determination of Ionospheric Conductivities from FUV Auroral Emissions'.

Keffer, Charles E.↗

Reducing Grid Costs while Abating Emissions: Opportunities for Flexible Building Loads

This study evaluates the value of technology-agnostic, shiftable flexible building loads in modeled 2030 and 2040 U.S. grids for four types of customers under three potential aggregated distributed energy resources programs. The value examined includes monetary value from providing grid services (e.g., energy, capacity, flexibility reserve, regulation reserve, contingency reserve) and from reducing greenhouse gas emissions. By comparing more than 845,164,800 simulated shifting opportunities, the study finds that the timing of consumption is critical for profit-driven customers. A program that is activated for 30 critical hours of system operation can result in to $73 per year in revenue for 1 kWh of shiftable load. Emission reduction, on the other hand, is best accumulated through a year-round program: 1 kWh of shiftable load can lead to up to 487.9 kg CO2e carbon reduction per year. The report also provides detailed insights on the trade-off between revenue and emissions, regional variation, short- versus medium- term value, and impacts from various building flexibility parameters, such as shifting window and dissipation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

MHD Modelling of Coronal Loops: Injection of High-Speed Chromospheric Flows

Context. Observations reveal a correspondence between chromospheric type II spicules and bright upward-moving fronts in the corona observed in the extreme-ultraviolet (EUV) band. However, theoretical considerations suggest that these flows are probably not the main source of heating in coronal magnetic loops. Aims. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes and the possible production of emission in the EUV band. Methods. We simulated the propagation of a dense 104 K chromospheric jet upward along a coronal loop by means of a 2D cylindrical MHD model that includes gravity, radiative losses, thermal conduction, and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool (approximately 0.8 MK) corona, linked through a steep transition region. In our reference model, the jet initial speed is 70 km per second, its initial density is 10(exp 11) per cubic centimeter, and the ambient uniform magnetic field is 10 G. We also explored other values of jet speed and density in 1D and different magnetic field values in 2D, as well as the jet propagation in a hotter (approximately 1.5 MK) background loop. Results. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock-front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about 10(exp 6) K. As a result, a bright moving front becomes visible in the 171 Angstrom channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. Conclusions. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

Sun: chromosphere↗

A Test of Sensitivity to Convective Transport in a Global Atmospheric CO2 Simulation

Two approximations to convective transport have been implemented in an offline chemistry transport model (CTM) to explore the impact on calculated atmospheric CO2 distributions. GlobalCO2 in the year 2000 is simulated using theCTM driven by assimilated meteorological fields from the NASA s Goddard Earth Observation System Data Assimilation System, Version 4 (GEOS-4). The model simulates atmospheric CO2 by adopting the same CO2 emission inventory and dynamical modules as described in Kawa et al. (convective transport scheme denoted as Conv1). Conv1 approximates the convective transport by using the bulk convective mass fluxes to redistribute trace gases. The alternate approximation, Conv2, partitions fluxes into updraft and downdraft, as well as into entrainment and detrainment, and has potential to yield a more realistic simulation of vertical redistribution through deep convection. Replacing Conv1 by Conv2 results in an overestimate of CO2 over biospheric sink regions. The largest discrepancies result in a CO2 difference of about 7.8 ppm in the July NH boreal forest, which is about 30% of the CO2 seasonality for that area. These differences are compared to those produced by emission scenario variations constrained by the framework of Intergovernmental Panel on Climate Change (IPCC) to account for possible land use change and residual terrestrial CO2 sink. It is shown that the overestimated CO2 driven by Conv2 can be offset by introducing these supplemental emissions.

Bian, H.↗

Role of Surface Wind and Vegetation Cover in Multi-decadal Variations of Dust Emission in the Sahara and Sahel

North Africa, the world's largest dust source, is non-uniform, consisting of a permanently arid region (Sahara), a semi-arid region (Sahel), and a relatively moist vegetated region (Savanna), each with very different rainfall patterns and surface conditions. This study aims to better understand the controlling factors that determine the variation of dust emission in North Africa over a 27-year period from 1982 to 2008, using observational data and model simulations. The results show that the model-derived Saharan dust emission is only correlated with the 10-m winds (W10m) obtained from reanalysis data, but the model-derived Sahel dust emission is correlated with both W10m and the Normalized Difference Vegetation Index (NDVI) that is obtained from satellite. While the Saharan dust accounts for 82 of the continental North Africa dust emission (1340-1570 Tg year(exp -1) in the 27-year average, the Sahel accounts for 17 with a larger seasonal and inter-annual variation (230-380 Tg year(exp -1), contributing about a quarter of the transatlantic dust transported to the northern part of South America. The decreasing dust emission trend over the 27-year period is highly correlated with W10m over the Sahara (R equals 0.92). Over the Sahel, the dust emission is correlated with W10m (R 0.69) but is also anti-correlated with the trend of NDVI (R equals 0.65). W10m is decreasing over both the Sahara and the Sahel between 1982 and 2008, and the trends are correlated (R equals 0.53), suggesting that Saharan Sahelian surface winds are a coupled system, driving the inter-annual variation of dust emission.

Inter-annual variation↗

Scoping study of detecting high harmonic fast waves in NSTX-U hot core plasma directly using beam emission spectroscopy

In this paper, we used radio frequency (RF) full wave simulation to explore the feasibility of generating a density beat oscillation formed by high harmonic fast waves (HHFWs) in the National Spherical Tokamak Experiment-Upgrade (NSTX-U) and of measuring it numerically with 2D beam emission spectroscopy (BES). This oscillation is of interest as it can be used to determine the injected RF fields. We began by computing HHFW fields in NSTX-U under a variety of experimental conditions and antenna configurations with Petra-M, a 3D RF wavefield solver. These fields were then used to evaluate the amplitude of δne,beat/ne0 in the volume of measurement, which we predict to be between 5 × 10−10 and 10−9, although it is highly dependent on plasma parameters. Trends in magnetic field, core temperature, core density, and antenna phasing are presented, and the parameters most suitable to this measurement are identified. Finally, we developed a synthetic BES diagnostic to compare the experimental BES signal against.

Rutherford, Grant (ORCID:0000000163620037)↗

Two-dimensional parametric design of short pulse laser driven iron opacity experiments

High intensity short pulse lasers are commonly used to create high temperature and high density plasmas. These plasmas are used to study fundamental properties such as the equation of state and opacity. The targets often use small-diameter, thin buried layers of the materials of interest, surrounded by a tamper of low atomic number material such as plastic. Computational modeling is used to design and interpret experiments with short pulse lasers. Most of the modeling to date has assumed one-dimensional plane-parallel geometry. In this paper, the effects of radial gradients in the irradiation of thin planar targets are studied with one- and two-dimensional radiation/hydrodynamic simulations. It is found that the main effect of radial gradients is the averaging of plasma conditions and x-ray emission over the pattern of irradiation. Additionally, differences between one- and two-dimensional simulations arise because the plasma conditions and x-ray emission are, in general, nonlinear functions of the temperature of the plasma and thereby also nonlinear functions of the irradiation intensity. The differences increase with the ratio of the buried layer radius to the laser spot radius. The root mean square difference in the inferred iron L-shell opacity is less than 30% when the ratio is less than 1.0 and about 70% when the ratio is 1.5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Meteorological Controls on Biomass Burning During Santa Ana Events in Southern California

Fires occurring during Santa Ana (SA) events in southern California are driven by extreme fire weather characterized by high temperatures, low humidities, and high wind speeds. We studied the controls on burned area and carbon emissions during two intensive SA burning periods in 2003 and 2007. We therefore used remote sensing data in parallel with fire weather simulations of the Weather and Regional Forecast model. Total carbon emissions were approximately 1800 gigagrams in 2003 and 900 gigagrams in 2007, based on a daily burned area and a fire emission model that accounted for spatial variability in fuel loads and combustion completeness. On a regional scale, relatively strong positive correlations were found between the daily Fosberg fire weather index and burned area/emissions (probability is less than 0.01). Our analysis provides a quantitative assessment of relationships between fire activity and weather during severe SA fires in southern California.

Santa Ana (SA)↗

Update on AeroCom Biomass Burning Emission Injection Height experiment (BBEIH)

In BBEIH, we introduced the biomass burning injection heights based on MISR stereo-derived plume-height retrievals (Val Martin et al., 2010; 2018) in four participating CTMs. Specifically, we proposed 4 simulations for the year 2008: 1) BASE: using the biomass burning emissions from Global Fire Emissions Database version 4 with small fires (GFED4s) and the model default fire emission injection height; 2) BBIH: Same as BASE, but using the injection height derived from the seasonally and regionally varying MISR-retrieved plume heights; 3) NOBB: no fire emission; 4) Same as BASE, but using Fire Energetics and Emissions Research version 1.0 (FEER 1.0). We will address the following scientific questions: 1) To what extent are the model simulations sensitive to the assumed biomass burning injection height? 2) In which regions/seasons/surface-types are the aforementioned sensitivities most important? We will update the multi-model inter-comparison with a focus on the vertical aerosol distribution in near-source characteristics and downwind plume evolution." In order to test the sensitivity of model results to smoke injection height, we proposed the biomass burning emission injection height experiment (BBEIH) in the AeroCom 2019 workshop. The description of implementation methods and input data for BBEIH can be found in https://wiki.met.no/aerocom/phase3-experiments#biomass_burning_emission_injection_height_experiment_bbeih.

Xiaohua Pan↗

Neutron capture measurements at UMass lowell research reactor

A new facility was designed around the thermal column beam port of the 1 MW Research Reactor at the University of Massachusetts Lowell. Thermal neutrons are collimated to a 1-inch diameter beam and incident on samples to induce the radiative neutron capture. New measurements of capture gamma rays are planned for Mn, Cu, Ni, Cr, and Gd samples in the next few years. The experiments will be carried out in close collaboration with the ENSDF (Evaluated Nuclear Structure Data File) evaluation group at Brookhaven National Laboratory. The gamma rays are measured using an array of high-resolution HPGe detectors. Some of the HPGe detectors use active Compton shields that were designed from scintillation detectors and improve significantly the signal-to-background in the measured gamma-ray spectra with HPGe. The experimental spectra will be validated by the GEANT4 simulations of the array and theoretical models of the emission of gamma rays from the compound nucleus. In conclusion, the new experimental results on the capture gamma ray intensities will be incorporated in future ENSDF evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost: Model Inputs and Outputs

This dataset contains model inputs and simulation outputs associated with the study “Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost.” The ecosys model was extended to represent thermokarst processes in laterally coupled two-dimensional transects and applied to the Scotty Creek basin (Northwest Territories, Canada). The ecosys model and documentation is available for download at https://github.com/jinyun1tang/ECOSYS. The dataset contains the inputs used to run the model for this study (meteorological forcing data (1950–2018) and site-specific soil and vegetation parameters) and model outputs used to generate Figure 3 and Figure 4 in the paper (simulated ground thermal and hydrologic states, vegetation dynamics, and carbon fluxes).

54 ENVIRONMENTAL SCIENCES↗

CO2 Capture Strategies via Mineralization with Industrial Waste Brines

Large coal-fired power plants (>500 MW) account for 30% of global CO2 emissions, and long-term management of this CO2 to is urgently needed mitigate global temperature increases. Sequestration of CO2 within stable mineral carbonates (e.g., CaCO3) represents an attractive emission reduction strategy because it offers a leakage-free alternative to geological storage of CO2 in an environmentally friendly form. We have previously described a mineralization process in which divalent cations are sourced from various waste streams (e.g., produced water and brackish water) and alkalinity is induced via regenerable ion-exchange materials (Bustillos et. al. Frontiers in Energy Research. 2020, 8, 352). In our process, aqueous carbonate-bearing streams with pH > 8 are produced by contacting fresh water and carbon dioxide with various ion-exchange materials (e.g., Na form zeolites or ion exchange resins). These streams are mixed with produced water containing varying concentrations (~0.01 – 1.0 M) of Ca2+ leading to the precipitation of solid calcium carbonate (PCC). This process has the advantages of using regenerable solids in a simple and continuous process to increase the pH of water by ion exchange instead of relying on the consumption of costly and unsustainable sources of alkalinity (e.g., sodium hydroxide). While once-through column experiments showed the above benefits, the same were yet to established in a steady-state process with recycle streams. In this work, we set up a process simulation to quantify the energy requirements and CO2 emissions associated with the process and seek optimal produced water compositions and CO2 concentrations (5 – 20 vol%). The process simulation was set up in ASPEN Plus using eRNTL as the thermodynamic property method and sequential modular strategy. Ion exchange alkaline solution was simulated using sodium hydroxide and validated against the experimental data obtained from once-through kinetic experiments. Nanofiltration and reverse osmosis membrane steps were also implemented for the separation of divalent cations and production of fresh water and a regeneration stream following mineralization. Sensitivity analysis was carried out using a range of produced water compositions (0.01 – 1.0 M Ca2+, 0.001 – 0.15 M Mg2+, 0.5 – 3.5 M Na+ and 0.0004 – 0.002 M Fe2+) according to the United States Geological Survey (USGS) database. Calcium carbonate yields increased with increasing CO2 concentrations and were maximized using produced water compositions with larger Ca2+ concentrations. Maximum calcium carbonate yields produced at 5 vol%, 12 vol% and 20 vol% CO2 were 2.3 mmol/L, 5.5 mmol/L, and 9.3 mmol/L, respectively, with the formation of brucite (a magnesium hydroxide phase, Mg(OH)2) and goethite (an iron hydroxide phase, FeOOH) as the primary contaminant phases (99% calcite, 0.6% brucite, 0.4% goethite), which agree with phases detected by XRD experimentally. These results indicate high purity calcium carbonate can be precipitated using industrial waste streams. Consequentially, energy consumption and net CO2 emissions were minimized where precipitated calcium carbonate was maximized for all produced water compositions and CO2 concentrations. Minimum energy consumptions were 0.21 kWh/ton CO2 processed, with 98% of the energy input required coming from the membrane filtration steps. Produced water compositions with large Na+ concentrations (> 0.5 M) were effective at reducing energy consumptions due to faster regeneration time of ion exchange materials. Additionally, calculated net CO2 emissions were negative for the process and ranged from -0.02 kg/ton CO2 to -0.15 kg/ton CO2 processed, indicating a low emission process. We will also present techno-economic assessment showing the economic benefits of the current process as an alternative to the addition of stoichiometric bases to induce alkalinity for the precipitation of CaCO3.

Simonetti, Dante↗

Evaluation of Stereology for Snow Microstructure Measurement and Microwave Emission Modeling: A Case Study

Reliable microstructure measurement of snow is a requirement for microwave radiative transfer model validation. Snow specific surface area (SSA) can be measured using stereological methods, in which snow samples are cast in the field and photographed in the laboratory. Processing stereology photographs manually by counting intersections of test cycloids with air–ice boundaries reduces the problems in binary segmentation. This paper is a case study to evaluate the repeatability of the manually stereology interpretation by two independent research groups. We further assessed how uncertainty in snow SSA influences simulated brightness temperature(TB) driven by the Microwave Emission Model of Layered Snowpacks (MEMLS), and how stereology compares to Near Infrared (NIR) camera and hand lens. Data was obtained from two alpine snow profiles from Steamboat Springs, Colorado. Results showed that stereological SSA values measured by two groups are highly consistent, and the ground radiometer measured TBat 19 and 37 GHz was successfully predicted (RMSE<3.8 K);simulations using NIR SSA and hand-lens geometric grain size (Dg)measurements have larger errors. This conclusion was not sensitive to uncertainty in the free parameters of TB modeling.

Microwave radiometry↗

Jetted and Turbulent Stellar Deaths: New LVK-detectable Gravitational-wave Sources

Upcoming LIGO–Virgo–KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events that come from black hole and neutron star binary mergers. Detection of noninspiral GW sources is also anticipated. We report the discovery of a new class of noninspiral GW sources—the end states of massive stars—that can produce the brightest simulated stochastic GW burst signal in the LVK bands known to date, and could be detectable in LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble—cocoon—inside of the star. We simulate such a system using state-of-the-art 3D general relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, ~10–100 Hz, over a characteristic jet activity timescale ~10–100 s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third-generation GW detectors with an estimated rate of ~10 events yr –1 . The detection rate drops to ~1% of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multimessenger events.

79 ASTRONOMY AND ASTROPHYSICS↗

Changes in Anthropogenic Precursor Emissions Drive Shifts in the Ozone Seasonal Cycle Throughout the Northern Midlatitude Troposphere

Simulations by six Coupled Model Intercomparison Project Phase 6 (CMIP6) Earth system models indicate that the seasonal cycle of baseline tropospheric ozone at northern midlatitudes has been shifting since the mid-20th century. Beginning in ∼ 1940, the magnitude of the seasonal cycle increased by ∼10 ppb (measured from seasonal minimum to maximum), and the seasonal maximum shifted to later in the year by about 3 weeks. This shift maximized in the mid-1980s, followed by a reversal – the seasonal cycle decreased in amplitude and the maximum shifted back to earlier in the year. Similar changes are seen in measurements collected from the 1970s to the present. The timing of the seasonal cycle changes is generally concurrent with the rise and fall of anthropogenic emissions that followed industrialization and the subsequent implementation of air quality emission controls. A quantitative comparison of the temporal changes in the ozone seasonal cycle at sites in both Europe and North America with the temporal changes in ozone precursor emissions across the northern midlatitudes found a high degree of similarity between these two temporal patterns. We hypothesize that changing precursor emissions are responsible for the shift in the ozone seasonal cycle; this is supported by the absence of such seasonal shifts in southern midlatitudes where anthropogenic emissions are much smaller. We also suggest a mechanism by which changing emissions drive the changing seasonal cycle: increasing emissions of NOx allow summertime photochemical production of ozone to become more important than ozone transported from the stratosphere, and increasing volatile organic compounds (VOCs) lead to progressively greater photochemical ozone production in the summer months, thereby increasing the amplitude of the seasonal ozone cycle. Decreasing emissions of both precursor classes then reverse these changes. The quantitative parameter values that characterize the seasonal shifts provide useful benchmarks for evaluating model simulations, both against observations and between models.

Ozone↗

Discrete electromagnetic emissions in planetary magnetospheres

The present state of knowledge and observations of discrete electromagnetic emissions detected in the magnetospheres of earth and the outer planets are reviewed. The history of ground observations and satellite observations is included. Active experiments studying wave propagation, wave-wave and wave-particle interactions, and the generation and triggering of discrete emissions are addressed. The various theories proposed to explain the phenomena are examined, and computer simulations of wave-particle interactions and the generation and triggering of discrete emissions are reviewed. Related phenomena observed at frequencies above both the electron cyclotron frequency and the plasma frequency are considered. Areas of future investigation are discussed.

Anderson, Roger R.↗