Search NASA⌕ Search

SEARCH · Search NASA

Results for “Plot”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Elevated temperature viscous remanent magnetization of natural and synthetic multidomain magnetite

The time-temperature relationship of magnetization is a subject of much interest and debate by paleomagnetists, rock magnetists, and magnetic anomaly modellers. We have investigated this relationship by studying the viscous remanent magnetization (VRM) of coarse-grained multidomain (MD) magnetite. Our experiments covered the temperature range from 22 to 400 C, times from minutes to months, and included both Australian granulite samples and multidomain magnetite samples synthesized by the glass ceramic method. VRM acquisition was found to generally increase with temperature but not always at the rate predicted from classical thermal fluctuation theories. Thermal cycling between room temperature (at which the measurements were made) and the VRM acquisition temperature sharply decreased the temperature dependence of the VRM. Room temperature VRM acquisition accelerates with time when plotted on a semilog plot, whereas at elevated temperature the curves are quasilinear against log(time) for both the natural and synthetic samples. This change in behavior may suggest a variation in the VRM acquisition mechanism as a function of temperature for MD magnetite. The granulites have a nearly linear increase in VRM acquisition rate with temperature whereas the glass ceramics display little change in the acquisition rate between 22 and 200 C, but increase by nearly a factor of 3 by 400 C. The increase in VRM of the glass ceramics between 200 and 400 C is in general qualitative agreement with thermal fluctuation theory. There was no systematic change in the rate of VRM acquisition with grain size for the multidomain magnetites used in this study. Elevated temperature (e.g., 400 C) VRM acquisition by the deep crustal granulites, if extrapolated over the Brunhes chron, would produce a magnetization of several A/m which, if true, is of the order required by models for the source of long-wavelength magnetic anomalies.

Kelso, Paul R.↗
NASA NTRS1994

Basalt generation at the Apollo 12 site. Part 1: New data, classification, and re-evaluation

New data are reported from five previously unanalyzed Apollo 12 mare basalts that are incorporated into an evaluation of previous petrogenetic models and classification schemes for these basalts. This paper proposes a classification for Apollo 12 mare basalts on the basis of whole-rock Mg# (molar 100*(Mg/(Mg+Fe))) and Rb/Sr ratio (analyzed by isotope dilution), whereby the ilmenite, olivine, and pigeonite basalt groups are readily distinguished from each other. Scrutiny of the Apollo 12 feldspathic 'suite' demonstrates that two of the three basalts previously assigned to this group (12031, 12038, 12072) can be reclassified: 12031 is a plagioclase-rich pigeonite basalt; and 12072 is an olivine basalt. Only basalt 12038 stands out as a unique sample to the Apollo 12 site, but whether this represents a single sample from another flow at the Apollo 12 site or is exotic to this site is equivocal. The question of whether the olivine and pigeonite basalt suites are co-magmatic is addressed by incompatible trace-element chemistry: the trends defined by these two suites when Co/Sm and Sm/Eu ratios are plotted against Rb/Sr ratio demonstrate that these two basaltic types cannot be co-magmatic. Crystal fractionation/accumulation paths have been calculated and show that neither the pigeonite, olivine, or ilmenite basalts are related by this process. Each suite requires a distinct and separate source region. This study also examines sample heterogeneity and the degree to which whole-rock analyses are representative, which is critical when petrogenetic interpretation is undertaken. Sample heterogeneity has been investigated petrographically (inhomogeneous mineral distribution) with consideration of duplicate analyses, and whether a specific sample (using average data) plots consistently upon a fractionation trend when a number of different compostional parameters are considered. Using these criteria, four basalts have been identified where reported analyses are not representative of the whole-rock composition: 12005, an ilmenite basalt; 12006 and 12036, olivine basalts; and 12031 previously classified as a feldspathic basalt, but reclassified as part of the pigeonite suite.

Neal, Clive R.↗
NASA NTRS1994

Sensible heat flux - Radiometric surface temperature relationship for eight semiarid areas

Measurements of sensible heat flux, radiometric surface temperature, air temperature, and wind speed made at eight semiarid rangeland sites were used to investigate the sensible heat flux-aerodynamic resistance relationship. The individual sites covered a wide range of vegetation (0.1-4 m tall) and cover (3%-95% bare soil) conditions. Mean values of k/B, a quantity related to the resistance of heat versus momentum transfer at the surface, for the individual sites varied between 3.5 and 12.5. A preliminary test of the utility of an excess resistance based on the mean value of k/B showed that the difference between the mean estimated and measured sensible heat fluxes varied +/- 60 W/sq m for the eight semiarid sites. For the eight sites the values of k/B were plotted against the roughness Reynolds number. The plot showed considerable scatter with values ranging between and beyond the theoretical curves for bluff rough and permeable rough surfaces.

Stewart, J. B.↗
NASA NTRS1994

Troilite formation kinetics and growth mechanism in the solar nebula

Troilite formation via the reaction Fe(s) + H2S(g) + H2(g) is the major mechanism for S retention in grains in the solar nebula. Thermodynamic calculations predict that troilite condenses from a solar composition gas. We present experimental results on the kinetics and growth of troilite crystals on Fe metal at temperature (450-650 C) and composition (50-1000 ppm H2S in H2) conditions similar to those in the solar nebula. The fraction of Fe reacted (based on gravimetric data) is plotted at 450, 505, 575, and 650 C. The thickness change of unreacted iron (measured by optical microscopy) is plotted at 575 and 650 C vs. time. the weight change per unit area varies as the square root of time at the lower temperatures and varies linearly with time at the highest temperature. The growth behavior along the lower isotherms is due to diffusion. This behavior suggests sulfide diffusion to the metal-sulfide interface and suggests Fe(2+) diffusion to the sulfide-gas interface. The reaction along the highest isotherm appears to be interface controlled. The formation of troilite crystals is a rapid process forming measurable layers in a few hours. The crystalgrowth is complicated. Initially there are intergrowths of troilite into the pure Fe metal. As the reaction progresses two distinct layers of troilite crystals form. One is in contact with the Fe metal and consists of small randomly oriented crystals with pore space between them. The outermost layer contains large crystals that are all oriented in the same direction. The intergrowth layer is much smaller at 650 C than at 575 C. This suggest that FeS nucleation is inhibited at the higher temperature, accounting for the initially slower reaction rate. Once nucleated, the reaction kinetics are apparently controlled by the growth of the crystals at the interface.

Lauretta, D. S.↗
NASA NTRS1995

Abundance ratios of oxygen, neon, and magnesium in solar active regions and flares: The FIP effect

Relative abundances of oxygen, neon, and magnesium have been derived for a sample of nine solar active regions, flares, and an erupting prominance by combining plots of the ion differential emission measures. The observations were photographed in the 300-600 A range by the Naval Research Laboratory (NRL) spectroheliograph on Skylab. Methods for deriving the Mg/Ne abundance ratio-which measures the separation between the low- first ionization potential (FIP) and high-FIP abundnace plateaus-have been described in previous papers. In this paper we describe the spectroscopic methods for deriving the O/Ne abundance ratio, which gives the ratio between two high-FIP elements. The plot of the O/Ne ratio versus the Mg/Ne ratio in the sample of nine Skylab events is shown. The variation in the Mg/Ne ratio by a factor of 6 is associated with a much smaller range in the O/Ne ratio. This is broadly consistent with the presence of the standard FIP pattern of abundances in the outer atmosphere of the Sun. However, a real change in the relative abundances of oxygen and neon by a factor of 1.5 cannot be excluded.

Widing, K. G.↗
NASA NTRS1995

Deep ROSAT HRI observations of the Orion nebula region

We present results from three deep ROSAT high-resolution imager observations of the Orion Nebula star-forming region. The X-ray images contain over 1500 cataloged stars in a roughly 0.8 sq deg region centered on the Trapezium. In all, 389 distinct X-ray sources have been detected, at least two-thirds of which are associated with a single proper-motion cluster member. X-ray emission is detected from stars of all spectral types, from massive O- and B-type components of the Trapezium to the coolest, low-mass pre-main-sequence (PMS) stars. In this paper we focus primarily on X-ray emission from the late-type PMS stars. Of the approximately 100 late-type cluster members with measured spectral types, approximately three-fourths have been detected; we have derived X-ray luminosity upper limits for the remaining stars. We found coronal X-ray emission turns on around spectral type F6, with the upper envelope of activity increasing with deceasing effective temperature. When plotted in an X-ray luminosity versus bolometric luminosity diagram, late-type PMS stars lie below a 'saturation' line corresponding to L(sub x)/L(sub bol) approximately 10(exp -3). For approximately solar-mass PMS stars, we find a median X-ray luminosity approximately 1 x 10(exp 30) ergs/s. The late G, K, and M stars exhibit nearly a two order of magnitude spread in X-ray luminosity and in L(sub x)/L(sub bol) at a given effective temperature. Plots of X-ray activity versus v sin i rotational velocity and rotational period appear to show no clear dependence of activity on rotation. However, because only a small fraction of late-type PMS stars in the Orion Nebula have measured v sin i or P(sub rot) and because of uncertainties in L(sub x) and L(sub x)/L(sub bol), we believe the data are not conclusive on this point. Light curves of the detected X-ray sources have revealed at least 10 strong X-ray flares with characteristic rise times greater than or approximately equal to 1 hr and decay times ranging from approximately 2 to 12 hr. All the flares have X-ray energies in excess of 3 x 10(exp 35) ergs. Many of the X-ray sources associated with late-type cluster members have been observed in more than one ROSAT HRI exposure, allowing us to look for variability on the approximately 1 yr timescale between HRI observations. A statistical analysis of the resulting count rate ratios suggests that at least one-fourth of the late-type cluster stars are significantly variable on this timescale.

Gagne, Marc↗
NASA NTRS1995

Software Aids In Graphical Depiction Of Flow Data

Interactive Data Display System (IDDS) computer program is graphical-display program designed to assist in visualization of three-dimensional flow in turbomachinery. Grid and simulation data files in PLOT3D format required for input. Able to unwrap volumetric data cone associated with centrifugal compressor and display results in easy-to-understand two- or three-dimensional plots. IDDS provides majority of visualization and analysis capability for Integrated Computational Fluid Dynamics and Experiment (ICE) system. IDDS invoked from any subsystem, or used as stand-alone package of display software. Generates contour, vector, shaded, x-y, and carpet plots. Written in C language. Input file format used by IDDS is that of PLOT3D (COSMIC item ARC-12782).

Stegeman, J. D.↗
NASA NTRS1995

Lidar Analyses

A brief description of enhancements made to the NASA MSFC coherent lidar model is provided. Notable improvements are the addition of routines to automatically determine the 3 dB misalignment loss angle and the backscatter value at which the probability of a good estimate (for a maximum likelihood estimator) falls to 50%. The ability to automatically generate energy/aperture parametrization (EAP) plots which include the effects of angular misalignment has been added. These EAP plots make it very easy to see that for any practical system where there is some degree of misalignment then there is an optimum telescope diameter for which the laser pulse energy required to achieve a particular sensitivity is minimized. Increasing the telescope diameter above this will result in a reduction of sensitivity. These parameterizations also clearly show that the alignment tolerances at shorter wavelengths are much stricter than those at longer wavelengths. A brief outline of the NASA MSFC AEOLUS program is given and a summary of the lidar designs considered during the program is presented. A discussion of some of the design trades is performed both in the text and in a conference publication attached as an appendix.

Spiers, Gary D.↗
NASA NTRS1995

Solar cell angle of incidence corrections

Literature on solar array angle of incidence corrections was found to be sparse and contained no tabular data for support. This lack along with recent data on 27 GaAs/Ge 4 cm by 4 cm cells initiated the analysis presented in this paper. The literature cites seven possible contributors to angle of incidence effects: cosine, optical front surface, edge, shadowing, UV degradation, particulate soiling, and background color. Only the first three are covered in this paper due to lack of sufficient data. The cosine correction is commonly used but is not sufficient when the incident angle is large. Fresnel reflection calculations require knowledge of the index of refraction of the coverglass front surface. The absolute index of refraction for the coverglass front surface was not known nor was it measured due to lack of funds. However, a value for the index of refraction was obtained by examining how the prediction errors varied with different assumed indices and selecting the best fit to the set of measured values. Corrections using front surface Fresnel reflection along with the cosine correction give very good predictive results when compared to measured data, except there is a definite trend away from predicted values at the larger incident angles. This trend could be related to edge effects and is illustrated by a use of a box plot of the errors and by plotting the deviation of the mean against incidence angle. The trend is for larger deviations at larger incidence angles and there may be a fourth order effect involved in the trend. A chi-squared test was used to determine if the measurement errors were normally distributed. At 10 degrees the chi-squared test failed, probably due to the very small numbers involved or a bias from the measurement procedure. All other angles showed a good fit to the normal distribution with increasing goodness-of-fit as the angles increased which reinforces the very small numbers hypothesis. The contributed data only went to 65 degrees from normal which prevented any firm conclusions about extreme angle effects although a trend in the right direction was seen. Measurement errors were estimated and found to be consistent with the conclusions that were drawn. A controlled experiment using coverglasses and cells from the same lots and extending to larger incidence angles would probably lead to further insight into the subject area.

Burger, Dale R.↗
NASA NTRS1996

Electrical Conductivity, Relaxation and the Glass Transition: A New Look at a Familiar Phenomenon

Annealed samples from a single melt of a 10 mol% K2O-90SiO2 glass were reheated to temperatures ranging from 450 to 800 C, held isothermally for 20 min, and then quenched in either air or a silicon oil bath. The complex impedance of both the annealed and quenched samples was measured as a function of temperature from 120 to 250 C using ac impedance spectroscopy from 1 Hz to 1 MHz. The dc conductivity, sigma(sub dc), was measured from the low frequency intercept of depressed semicircle fits to the complex impedance data. When the sigma(sub dc) at 150 C was plotted against soak temperature, the results fell into three separate regions that are explained in terms of the glass structural relaxation time, tau(sub S). This sigma(sub dc) plot provides a new way to look the glass transition range, Delta T(sub r). In addition, sigma(sub dc) was measured for different soak times at 550 C, from which an average relaxation time of 7.3 min was calculated. It was found that the size and position of the Delta T(sub r) is controlled by both the soak time and cooling rate.

Angel, Paul W.↗
NASA NTRS1995

Interactive flutter analysis and parametric study for conceptual wing design

An interactive computer program was developed for wing flutter analysis in the conceptual design stage. The objective was to estimate the flutter instability boundary of a flexible cantilever wing, when well defined structural and aerodynamic data are not available, and then study the effect of change in Mach number, dynamic pressure, torsional frequency, sweep, mass ratio, aspect ratio, taper ratio, center of gravity, and pitch inertia, to guide the development of the concept. The software was developed on MathCad (trademark) platform for Macintosh, with integrated documentation, graphics, database and symbolic mathematics. The analysis method was based on nondimensional parametric plots of two primary flutter parameters, namely Regier number and Flutter number, with normalization factors based on torsional stiffness, sweep, mass ratio, aspect ratio, center of gravity location and pitch inertia radius of gyration. The plots were compiled in a Vaught Corporation report from a vast database of past experiments and wind tunnel tests. The computer program was utilized for flutter analysis of the outer wing of a Blended Wing Body concept, proposed by McDonnell Douglas Corporation. Using a set of assumed data, preliminary flutter boundary and flutter dynamic pressure variation with altitude, Mach number and torsional stiffness were determined.

Mukhopadhyay, Vivek↗
NASA NTRS1995

Study of the rotational-level and temperature dependence of the quenching rate of OH fluorescence due to collisions with water molecules

The importance of the OH radical as an intermediate in many combustion reactions and in atmospheric photochemistry has led many researchers to use it as a diagnostic tool in these processes. The amount of data that has been acquired over the years for this radical is quite considerable. However, the quenching rate of OH with water molecules as a function of temperature and the rotational level of the excited state is not very well understood. The motivation of the studies undertaken is to bridge the gap between the low temperature measurements and the high temperature ones reported in the literature. The technique generally employed in these diagnostics is laser-induced fluorescence (LIF), through which rotational state selective excitation of the radical is possible. Furthermore, in a combustion medium, water is produced in abundance so that knowledge of the quenching rate of OH due to water molecules plays a crucial role in interpreting the data. In general, the precursor to an understanding of the collisional quenching rates of OH involves a characterization of the mode in which the radical is produced; the resulting rotational and translational distribution, followed by a measurement of the OH temperature; and ultimately obtaining the rate constants from the pressure dependence of the fluorescence signal. The experimental implementation of these measurements therefore involved, as a first step, the production of the OH radicals in a microwave discharge cell using water vapor as the source, wherein a hydrogen atom is abstracted from H2O. The second step involved the absorption of photons from the frequency-doubled output of a pulsed amplified, single-frequency cw ring dye laser. By tuning the laser to the peak of the transition and observing the fluorescence decay after the laser pulse, the lifetime of the OH in a particular rotational electronic state was determined (tau = 1.4 microseconds for Q(sub 1)(3)). Knowledge of this parameter led to a determination of the quenching rate. By varying the water vapor pressure in the cell and measuring the lifetime as a function of pressure a linear plot of the quenching rate as a function of pressure was obtained. Using this plot, the quenching cross section was deduced. It has therefore been possible to measure the local translational temperature and the quenching cross section with one laser system.

Koker, Edmond B.↗
NASA NTRS1996

Surface Tension Driven Convection Experiment (STDCE)

This document reports the results obtained from the Surface Tension Driven Convection Experiment (STDCE) conducted aboard the USML-1 Spacelab in 1992. The experiments used 10 cSt silicone oil placed in an open circular container that was 10 cm wide and 5 cm deep. Thermocapillary flow was induced by using either a cylindrical heater placed along the container centerline or by a CO2 laser. The tests were conducted under various power settings, laser beam diameters, and free surface shapes. Thermistors located at various positions in the test section recorded the temperature of the fluid, heater, walls, and air. An infrared imager was used to measure the free surface temperature. The flow field was studied by flow visualization and the data was analyzed by a PTV technique. The results from the flow visualization and the temperature measurements are compared with the numerical analysis that was conducted in conjunction with the experiment. The compared results include the experimental and numerical velocity vector plots, the streamline plots, the fluid temperature, and the surface temperature distribution.

Ostrach, S.↗
NASA NTRS1996

Effects of Core Thickness and Fiber Orientation on Composite Beam Stiffness

This laboratory experiment is intended for students in an introductory polymer materials and processes course or engineering materials course. It can be conducted as an introduction to the hand lay-up process, with additional observations regarding the stiffness of the completed composite beams based on core thickness and fiber orientation. Students gain hands-on experience with the hand lay-up process by constructing glass/epoxy composite panels. Each lab group produces a panel with different core thickness or fiber orientation. The panels are then cut into strips and tested for flexural stiffness in a three-point bending fixture. Students plot deflection versus load data for composite beams with two different fiber orientations, two core thicknesses and one beam with laminate plies only (no core). The deflection plots highlight the effects of core thickness and fiber orientation on composite beam stiffness.

Hillsman, Vernon S.↗
NASA NTRS1996

Computer-Aided System Engineering and Analysis (CASE/A) Programmer's Manual, Version 5.0

The Computer Aided System Engineering and Analysis (CASE/A) Version 5.0 Programmer's Manual provides the programmer and user with information regarding the internal structure of the CASE/A 5.0 software system. CASE/A 5.0 is a trade study tool that provides modeling/simulation capabilities for analyzing environmental control and life support systems and active thermal control systems. CASE/A has been successfully used in studies such as the evaluation of carbon dioxide removal in the space station. CASE/A modeling provides a graphical and command-driven interface for the user. This interface allows the user to construct a model by placing equipment components in a graphical layout of the system hardware, then connect the components via flow streams and define their operating parameters. Once the equipment is placed, the simulation time and other control parameters can be set to run the simulation based on the model constructed. After completion of the simulation, graphical plots or text files can be obtained for evaluation of the simulation results over time. Additionally, users have the capability to control the simulation and extract information at various times in the simulation (e.g., control equipment operating parameters over the simulation time or extract plot data) by using "User Operations (OPS) Code." This OPS code is written in FORTRAN with a canned set of utility subroutines for performing common tasks. CASE/A version 5.0 software runs under the VAX VMS(Trademark) environment. It utilizes the Tektronics 4014(Trademark) graphics display system and the VTIOO(Trademark) text manipulation/display system.

Knox, J. C.↗
NASA NTRS1996

An Interactive Software for Conceptual Wing Flutter Analysis and Parametric Study

An interactive computer program was developed for wing flutter analysis in the conceptual design stage. The objective was to estimate the flutter instability boundary of a flexible cantilever wing, when well-defined structural and aerodynamic data are not available, and then study the effect of change in Mach number, dynamic pressure, torsional frequency, sweep, mass ratio, aspect ratio, taper ratio, center of gravity, and pitch inertia, to guide the development of the concept. The software was developed for Macintosh or IBM compatible personal computers, on MathCad application software with integrated documentation, graphics, data base and symbolic mathematics. The analysis method was based on non-dimensional parametric plots of two primary flutter parameters, namely Regier number and Flutter number, with normalization factors based on torsional stiffness, sweep, mass ratio, taper ratio, aspect ratio, center of gravity location and pitch inertia radius of gyration. The parametric plots were compiled in a Vought Corporation report from a vast data base of past experiments and wind-tunnel tests. The computer program was utilized for flutter analysis of the outer wing of a Blended-Wing-Body concept, proposed by McDonnell Douglas Corp. Using a set of assumed data, preliminary flutter boundary and flutter dynamic pressure variation with altitude, Mach number and torsional stiffness were determined.

Mukhopadhyay, Vivek↗
NASA NTRS1996

Electrochemical Impedance Spectroscopy of Conductive Polymer Coatings

Electrochemical impedance spectroscopy (EIS) was used to investigate the corrosion protection performance of twenty nine proprietary conductive polymer coatings for cold rolled steel under immersion in 3.55 percent NaCl. Corrosion potential as well as Bode plots of the data were obtained for each coating after one hour immersion, All coatings, with the exception of one, have a corrosion potential that is higher in the positive direction than the corrosion potential of bare steel under the same conditions. Group A consisted of twenty one coatings with Bode plots indicative of the capacitive behavior characteristic of barrier coatings. An equivalent circuit consisting of a capacitor in series with a resistor simulated the experimental EIS data for these coatings very well. Group B consisted of eight coatings that exhibited EIS spectra showing an inflection point which indicates that two time constants are present. This may be caused by an electrochemical process taking place which could be indicitive of coating failing. These coatings have a lower impedance that those in Group A.

Calle, Luz Marina↗
NASA NTRS1995

Analytic MHD Theory for Earth's Bow Shock at Low Mach Numbers

A previous MHD theory for the density jump at the Earth's bow shock, which assumed the Alfven M(A) and sonic M(s) Mach numbers are both much greater than 1, is reanalyzed and generalized. It is shown that the MHD jump equation can be analytically solved much more directly using perturbation theory, with the ordering determined by M(A) and M(s), and that the first-order perturbation solution is identical to the solution found in the earlier theory. The second-order perturbation solution is calculated, whereas the earlier approach cannot be used to obtain it. The second-order terms generally are important over most of the range of M(A) and M(s) in the solar wind when the angle theta between the normal to the bow shock and magnetic field is not close to 0 deg or 180 deg (the solutions are symmetric about 90 deg). This new perturbation solution is generally accurate under most solar wind conditions at 1 AU, with the exception of low Mach numbers when theta is close to 90 deg. In this exceptional case the new solution does not improve on the first-order solutions obtained earlier, and the predicted density ratio can vary by 10-20% from the exact numerical MHD solutions. For theta approx. = 90 deg another perturbation solution is derived that predicts the density ratio much more accurately. This second solution is typically accurate for quasi-perpendicular conditions. Taken together, these two analytical solutions are generally accurate for the Earth's bow shock, except in the rare circumstance that M(A) is less than or = 2. MHD and gasdynamic simulations have produced empirical models in which the shock's standoff distance a(s) is linearly related to the density jump ratio X at the subsolar point. Using an empirical relationship between a(s) and X obtained from MHD simulations, a(s) values predicted using the MHD solutions for X are compared with the predictions of phenomenological models commonly used for modeling observational data, and with the predictions of a modified phenomenological model proposed recently. The similarities and differences between these results are illustrated using plots of X and a(s) predicted for the Earth's bow shock. The plots show that the new analytic solutions agree very well with the exact numerical MHD solutions and that these MHD solutions should replace the corresponding phenomenological relations in comparisons with data. Furthermore, significant differences exist between the standoff distances predicted at low M(A) using the MHD models versus those predicted by the new modified phenomenological model. These differences should be amenable to observational testing.

Grabbe, Crockett L.↗