Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational”

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 523 records · Page 29

Variational Calculations of Ro-Vibrational Energy Levels and Transition Intensities for Tetratomic Molecules

A description is given of an algorithm for computing ro-vibrational energy levels for tetratomic molecules. The expressions required for evaluating transition intensities are also given. The variational principle is used to determine the energy levels and the kinetic energy operator is simple and evaluated exactly. The computational procedure is split up into the determination of one dimensional radial basis functions, the computation of a contracted rotational-bending basis, followed by a final variational step coupling all degrees of freedom. An angular basis is proposed whereby the rotational-bending contraction takes place in three steps. Angular matrix elements of the potential are evaluated by expansion in terms of a suitable basis and the angular integrals are given in a factorized form which simplifies their evaluation. The basis functions in the final variational step have the full permutation symmetries of the identical particles. Sample results are given for HCCH and BH3.

Schwenke, David W.↗

X-ray versus Optical Variations in the Seyfert 1 Nucleus NGC 3516: A Puzzling Disconnectedness

We present optical broadband (B and R) observations of the Seyfert 1 nucleus NGC 3516, obtained at Wise Observatory from March 1997 to March 2002, contemporaneously with X-ray 2-10 keV measurements with RXTE. With these data we increase the temporal baseline of this dataset to 5 years, more than triple to the coverage we have previously presented for this object. Analysis of the new data does not confirm the 100-day lag of X-ray behind optical variations, tentatively reported in our previous work. Indeed, excluding the first year's data, which drive the previous result, there is no significant correlation at any Lag between the X-ray and optical bands. We also find no correlation at any lag between optical flux and various X-ray hardness ratios. We conclude that the close relation observed between the bands during the first year of our program was either a fluke, or perhaps the result of the exceptionally bright state of NGC 3516 in 1997, to which it has yet to return. Reviewing the results of published joint X-ray and UV/optical Seyfert monitoring programs, we speculate that there are at least two components or mechanisms contributing to the X-ray continuum emission up to 10 key: a soft component that is correlated with UV/optical variations on timescales approx. greater than 1 day, and whose presence can be detected when the source is observed at low enough energies (approx. 1 keV), is unabsorbed, or is in a sufficiently bright phase; and a hard component whose variations are uncorrelated with the UV/optical.

Maoz, Dan↗

Toward Failure Modeling In Complex Dynamic Systems: Impact of Design and Manufacturing Variations

When designing vehicle vibration monitoring systems for aerospace devices, it is common to use well-established models of vibration features to determine whether failures or defects exist. Most of the algorithms used for failure detection rely on these models to detect significant changes during a flight environment. In actual practice, however, most vehicle vibration monitoring systems are corrupted by high rates of false alarms and missed detections. Research conducted at the NASA Ames Research Center has determined that a major reason for the high rates of false alarms and missed detections is the numerous sources of statistical variations that are not taken into account in the. modeling assumptions. In this paper, we address one such source of variations, namely, those caused during the design and manufacturing of rotating machinery components that make up aerospace systems. We present a novel way of modeling the vibration response by including design variations via probabilistic methods. The results demonstrate initial feasibility of the method, showing great promise in developing a general methodology for designing more accurate aerospace vehicle vibration monitoring systems.

Tumer, Irem Y.↗

Evaluating Ground-based Proxies for Solar Irradiance Variation

In order to determine what ground-based proxies are best for evaluating solar irradiance variation before the advent of space observations, it is necessary to test these proxies against space observations. We have tested sunspot number, total sunspot area, and sunspot umbral area against the Nimbus-7 measurements of total solar irradiance variation cover the eleven year period 1980-1990. The umbral area yields the best correlation and the total sunspot area yields the poorest. Reasons for expecting the umbral area to yield the best correlation are given, the statistical procedure followed to obtain the results is described, and the value of determining the best proxy is discussed. The latter is based upon the availability of an excellent database from the Greenwich Observatory obtained over the period 1876-1976, which can be used to estimate the total solar irradiance variation before sensitive space observations were available. The ground-based observations used were obtained at the Coimbra Solar Observatory. The analysis was done at Goddard using these data and data from the Nimbus-7 satellite.

Oegerle, William↗

Spatial Variation of Surface Moisture Fluxes in SGP

Analysis of aircraft data in "spatial variations of surface moisture flux from aircraft data" indicates that the impact of small-scale surface heterogeneity on the spatial variation of surface moisture fluxes into the atmosphere is reduced by horizontal mixing. This mixing generally increases with the development of the daytime convective mixed layer, thus reducing the relative influence of surface heterogeneity on the spatial variation of moisture fluxes.

Mahrt, L.↗

Cross-Correlation Analysis of UV Profile Variations and Nonradial Pulsations in Be Stars

This study is an investigation of the role of photospheric nonradial pulsations in the mass loss processes of Be stars. Some 18 Be stars were the focus of intensive, multi-wavelength campaigns between 1985 and 1996 with the NASA/ESA International Ultraviolet Explorer Satellite (IUE) . The campaigns were designed to monitor variations in the UV flux and stellar wind lines with near simultaneous, high signal-to-noise, optical spectroscopy of photospheric lines. The high quality, ground-based spectra were required to record the subtle, blue-to-red moving bump patterns in spectral line profiles that result from nonradial pulsations. We found that the same fine variations in the UV photospheric lines can be extracted from noisy IUE spectra by cross-correlating each target spectrum with a narrow-lined standard spectrum. Time series analyses of the cross-correlation functions then reveal the periodic signals, amplitudes, phases, and modal identifications of the pulsations. Our goals are to compare such results from UV and optical time series in order to model the stellar photospheric temperature and velocity fluctuations, and to search for correlated UV wind line variations to determine whether there is a connection between pulsation and mass loss.

Gies, Douglas R.↗

Controls on Variations of Surface Energy, Water, and Carbon Budgets within Large-Scale Amazon Basin

A key research focus of the LBA Research Program is understanding the space-time variations in interlinked surface energy, water, and carbon budgets, the controls on these variations, and the implications of these controls on the carbon sequestering capacity of the large scale forest-pasture system that dominates the Amaz6nia landscape. Quantification of these variations and controls are investigated by a combination of in situ measurements, remotely sensed measurements from space, and a realistically forced hydrometeorological model coupled to a carbon assimilation model, capable of simulating details within the surface energy and water budgets along with the principle processes of photosynthesis and respiration. Herein we describe the results of an investigation concerning the space-time controls of carbon sources and sinks distributed over the large scale Amazon basin. The results are derived from a carbon-water-energy budget retrieval system for the large scale Amazon basin, which uses a coupled carbon assimilation-hydrometeorological model as an integrating system, forced by both in situ meteorological measurements and remotely sensed radiation and precipitation fluxes obtained from a combination of GOES, SSM/I, TOMS, and TRh4M satellite measurements. Results include validation of (a) retrieved surface radiation and precipitation fluxes based on 30-min averaged surface measurements taken at Ji-Parani in Rondania and Manaus in Amazonas, and (b) modeled sensible, latent, and C02 fluxes based on tower measurements taken at Reserva Jaru, Manaus and Fazenda Nossa Senhora. The space-time controls on carbon sequestration are partitioned into sets of factors classified by: (1) above canopy meteorology, (2) incoming surface radiation, (3) precipitation interception, and (4) indigenous stomatal processes varied over the different land covers of pristine rainforest, partially, and fully logged rainforests, and pasture lands. These are the principle meteorological, thermodynamical, hydrological, and biophysical control paths which perturb net carbon fluxes and sequestration, produce time-space switching of carbon sources and sinks, undergo modulation through atmospheric boundary layer feedbacks, and respond to any discontinuous intervention on the landscape itself such as produced by human intervention in converting rainforest to pasture or conducting selective/clearcut logging operations. The results demonstrate how relative carbon sequestration capacity of the Amazonian ecosystem responds to these controls, and how interpretation of space-time heterogeneities in carbon sequestration depends on a fairly exact quantification of the interacting non-linear properties of photosynthesis in response to incoming solar flux, air-canopy temperatures, and leaf water interception -- and soil respiration in response to upper layer soil temperature and water content. The results also show how the interpretation of the control processes is highly sensitive to the scales at which the surface fluxes are analyzed.

Smith, Eric A.↗

Study of CFD Variation on Transport Configurations from the Second Drag-Prediction Workshop

This paper describes and analyzes a series of nearly 90 CFD test cases performed as a contribution to the second Drag Prediction Workshop, held in association with the AIAA in June 2003. Two configurations are included: DLR-F6 wing-body and wing-body-nacelle-pylon. The ability of CFD to predict the drag, lift, and pitching moment from experiment-including the "delta" arising from the addition of the nacelle and pylon-is assessed. In general, at a fixed angle of attack CFD overpredicts lift, but predicts the delta C (sub L) reasonably well. At low lift levels (C (sub L) less than 0.3)), delta C (sub D) is 20-30 drag counts (30-45%) high. At the target lift coefficient of C(sub L) = 0.5, delta C (sub D) is overpredicted by between 11-16 counts. However, the primary contribution of this paper is mot so much the assessment of CFD against experiment, but rather a detailed assessment and analysis of CFD variation. The series of test cases are designed to determine the sensitivity/variability of CFD to a variety of factors, including grid, turbulence model, transition code, and viscous model. Using medium-level grids (6-11 million points) at the target lift coefficient, the maximum variation in drag due to different grids is 5-11 drag counts, due to code is 5-10 counts, due to turbulence model is 7-15 counts, due to transition is 10-11 counts, and due to viscous model is 4-5 counts. Other specific variations are described in the paper.

Christopher L Rumsey↗

Limitation of Ground-based Estimates of Solar Irradiance Due to Atmospheric Variations

The uncertainty in ground-based estimates of solar irradiance is quantitatively related to the temporal variability of the atmosphere's optical thickness. The upper and lower bounds of the accuracy of estimates using the Langley Plot technique are proportional to the standard deviation of aerosol optical thickness (approx. +/- 13 sigma(delta tau)). The estimates of spectral solar irradiance (SSI) in two Cimel sun photometer channels from the Mauna Loa site of AERONET are compared with satellite observations from SOLSTICE (Solar Stellar Irradiance Comparison Experiment) on UARS (Upper Atmospheric Research Satellite) for almost two years of data. The true solar variations related to the 27-day solar rotation cycle observed from SOLSTICE are about 0.15% at the two sun photometer channels. The variability in ground-based estimates is statistically one order of magnitude larger. Even though about 30% of these estimates from all Level 2.0 Cimel data fall within the 0.4 to approx. 0.5% variation level, ground-based estimates are not able to capture the 27-day solar variation observed from SOLSTICE.

Wen, Guoyong↗

Meridional Variations of C2H2 and C2H6 in Jupiter's Atmosphere from Cassini CIRS Infrared Spectra

The abundances of hydrocarbons such as acetylene (C2H2) and ethane (C2H6) in Jupiter's atmosphere are important physical quantities, constraining our models of the chemical and dynamical processes. However, our knowledge of these quantities and their vertical and latitudinal variations has remained sparse. The flyby of the Cassini spacecraft with Jupiter at the end of 2000 provided an excellent opportunity to observe the infrared spectrum with the Composite Infrared Spectrometer (CIRS) instrument, mapping the spatial variation of emissions from 10-1400 cm-1. CIRS spectra taken at the highest resolution (0.5 cm-1) in early December 2000 have been analysed to infer atmospheric temperatures in the stratosphere at 0.5-20 mbar via the v4 of CH4, and in the troposphere at 100-400 mbar, via the hydrogen collision-induced continuum absorption at 600-800 cm. Simultaneously, we have searched for meridional abundance variations in C2H2 and C2H6 via the v5 and vg bands respectively. Tropospheric absorption and stratospheric emission are highly anti-correlated at the CIM resolution, introducing a non-uniqueness into the retrievals, which means that vertical gradient and column abundance cannot be simultaneously found without additional constraints. If we assume the profile shapes from photochemical model calculations, we show that the column abundance of C2H2 must decrease sharply towards the poles, while C2H6 is constant or slightly increasing. The relevance of these results to current photochemical and dynamical knowledge of Jupiter's atmosphere is discussed.

Nixon, C. A.↗

Neoproterozoic variations in the C-isotopic composition of seawater: stratigraphic and biogeochemical implications

The recent proliferation of stratigraphic studies of delta 13C variation in carbonates and organic C in later Neoproterozoic and basal Cambrian successions (approximately 850-530 Ma) indicates a strong oscillating trend in the C-isotopic composition of surface seawater. Alone, this trend does not adequately characterize discrete intervals in Neoproterozoic time. However, integrated with the vectorial signals provided by fossils and Sr-isotopic variations, C isotope chemostratigraphy facilitates the interbasinal correlation of later Neoproterozoic successions. Results of these studies are evaluated in terms of four stratigraphic intervals: (1) the Precambrian/Cambrian boundary, (2) the post-Varanger terminal Proterozoic, (3) the late Cryogenian, and (4) the early Cryogenian. Where biostratigraphic or radiometric data constrain the age of Neoproterozoic sedimentary sequences, secular variations in C and Sr isotopes can provide a level of stratigraphic resolution exceeding that provided by fossils alone. Isotopic data place strong constraints on the chemical evolution of seawater, linking it to major tectonic and paleoclimatic events. They also provide a biogeochemical framework for the understanding of the initial radiation of macroscopic metazoans, which is associated stratigraphically, and perhaps causally, with a global increase in the burial of organic C and a concomitant rise of atmospheric O2.

NASA Discipline Exobiology↗

Millimeter-scale variations of stable isotope abundances in carbonates from banded iron-formations in the Hamersley Group of Western Australia

Several diamond drill cores from formations within the Hamersley Group of Western Australia have been studied for evidence of short-range variations in the isotopic compositions of the carbonates. For a set of 32 adjacent microbands analyzed in a specimen from the Marra Mamba Iron Formation, carbon isotope compositions of individual microbands ranged from -2.8 to -19.8 per mil compared to PDB and oxygen isotope compositions ranged from 10.2 to 20.8 per mil compared to SMOW. A pattern of alternating abundances was present, with the average isotopic contrasts between adjacent microbands being 3.0 per mil for carbon and 3.1 per mil for oxygen. Similar results were obtained for a suite of 34 microbands (in four groups) from the Bruno's Band unit of the Mount Sylvia Formation. Difficulties were experienced in preparing samples of single microbands from the Dales Gorge Member of the Brockman Iron Formation, but overall isotopic compositions were in good agreement with values reported by previous authors. Chemical analyses showed that isotopically light carbon and oxygen were correlated with increased concentrations of iron. The preservation of these millimeter-scale variations in isotopic abundances is interpreted as inconsistent with a metamorphic origin for the isotopically light carbon in the BIF carbonates. A biological origin is favored for the correlated variations in 13C and Fe, and it is suggested that the 13C-depleted carbonates may derive either from fermentative metabolism or from anaerobic respiration. A model is presented in which these processes occur near the sediment-water interface and are coupled with an initial oxidative precipitation of the iron.

NASA Program Exobiology↗

Isotopic compositions of carbonates and organic carbon from upper Proterozoic successions in Namibia: stratigraphic variation and the effects of diagenesis and metamorphism

The carbon isotope geochemistry of carbonates and organic carbon in the late Proterozoic Damara Supergroup of Namibia, including the Nama, Witvlei, and Gariep groups on the Kalahari Craton and the Mulden and Otavi groups on the Congo Craton, has been investigated as an extension of previous studies of secular variations in the isotopic composition of late Proterozoic seawater. Subsamples of microspar and dolomicrospar were determined, through petrographic and cathodoluminescence examination, to represent the "least-altered" portions of the rock. Carbon-isotopic abundances in these phases are nearly equal to those in total carbonate, suggesting that 13C abundances of late Proterozoic fine-grained carbonates have not been significantly altered by meteoric diagenesis, although 18O abundances often differ significantly. Reduced and variable carbon-isotopic differences between carbonates and organic carbon in these sediments indicate that isotopic compositions of organic carbon have been altered significantly by thermal and deformational processes, likely associated with the Pan-African Orogeny. Distinctive stratigraphic patterns of secular variation, similar to those noted in other, widely separated late Proterozoic basins, are found in carbon-isotopic compositions of carbonates from the Nama and Otavi groups. For example, in Nama Group carbonates delta 13C values rise dramatically from -4 to +5% within a short stratigraphic interval. This excursion suggests correlation with similar excursions noted in Ediacaran-aged successions of Siberia, India, and China. Enrichment of 13C (delta 13C> +5%) in Otavi Group carbonates reflects those in Upper Riphean successions of the Akademikerbreen Group, Svalbard, its correlatives in East Greenland, and the Shaler Group, northwest Canada. The widespread distribution of successions with comparable isotopic signatures supports hypotheses that variations in delta 13C reflect global changes in the isotopic composition of late Proterozoic seawater. Within the Damara basin, carbon-isotopic compositions of carbonates provide a potentially useful tool for the correlation of units between the Kalahari and Congo cratons. Carbonates depleted in 13C were deposited during and immediately following three separate glacial episodes in Namibia. The correspondence between ice ages and negative delta 13C excursions may reflect the effects of lowered sea levels; enhanced circulation of deep, cold, O2-rich seawater; and/or the upwelling of 13C-depleted deep water. Iron-formation is additionally associated with one of the glacial horizons, the Chuos tillite. Carbon-13 enriched isotopic abundances in immediately pre-glacial carbonates suggest that oceanographic conditions favored high rates of organic burial. It is likely that marine waters were stratified, with deep waters anoxic. A prolonged period of ocean stratification would permit the build-up of ferrous iron, probably from hydrothermal sources. At the onset of glaciation, upwelling would have brought 13C-depleted and iron-rich deep water onto shallow shelves where contact with cold, oxygenated surface waters led to the precipitation of ferric iron.

NASA Program Exobiology↗

The effects of UV radiation A and B on diurnal variation in photosynthesis in three taxonomically and ecologically diverse microbial mats

Photosynthetic primary production, the basis of most global food chains, is inhibited by UV radiation. Evaluating UV inhibition is therefore important for assessing the role of natural levels of UV radiation in regulating ecosystem behavior as well as the potential impact of stratospheric ozone depletion on global ecosystems. As both photosynthesis and UV fluxes are subject to diurnal variations, we examined the diurnal variability of the effect of UV radiation on photosynthesis in three diverse algal mats. In one of the mats (Cyanidium caldarium) a small mean decrease in primary productivity over the whole day occurred when both UVA and UVB were screened out. In two of the mats (Lyngbya aestuarii and Zygogonium sp.) we found a mean increase in the total primary productivity over the day when UVB alone was screened and a further increase when UVA and UVB were both screened out. Variations in the effects of UV radiation were found at different times of the day. This diurnal variability may be because even under the same solar radiation flux, there are different factors that may control photosynthetic rate, including nutritional status and other physiological processes in the cell. The results show the importance of assessing the complete diurnal productivity. For some of the time points the increase in the mean was still within the standard deviations in primary productivity, illustrating the difficulty in dissecting UV effects from other natural variations.

Ultraviolet Rays↗

Muscle recruitment variations during wrist flexion exercise: MR evaluation

OBJECTIVE: Many exercise protocols used in physiological studies assume homogeneous and diffuse muscle recruitment. To test this assumption during a "standard" wrist flexion protocol, variations in muscle recruitment were assessed using MRI in eight healthy subjects. MATERIALS AND METHODS: Variations were assessed by comparing the right to the left forearms and the effect of slight (15 degrees) pronation or supination at the wrist. RESULTS: Postexercise imaging showed focal regions of increased signal intensity (SI), indicating relatively strong recruitment, most often in entire muscles, although occasionally only in subvolumes of muscles. In 15 of 26 studies, flexor carpi radialis (FCR) showed more SI than flexor carpi ulnaris, while in 11 studies SI in these muscles increased equivalently. Relatively greater FCR recruitment was seen during pronation and/or use of the nondominant side. Palmaris longus, a wrist flexor, did not appear recruited in 4 of 11 forearms in which it was present. A portion of the superficial finger flexor became hyperintense in 89% of studies, while recruitment of the deep finger flexor was seen only in 43%. CONCLUSION: Inter- and intraindividual variations in forearm muscle recruitment should be anticipated in physiological studies of standard wrist flexion exercise protocols.

NASA Discipline Cardiopulmonary↗

Zonal-Mean Temperature Variations Inferred from SABER Measurements on TIMED Compared with UARS Observations

In the Numerical Spectral Model (NSM, Mayr et al., 2003), small-scale gravity waves propagating in the north/south direction can generate zonal mean (m = 0) meridional wind oscillations with periods between 2 and 4 months. These oscillations tend to be confined to low latitudes and have been interpreted to be the meridional counterpart of the wave-driven Quasi Biennial Oscillation in the zonal circulation. Wave driven meridional winds across the equator should generate, due to dynamical heating and cooling, temperature oscillations with opposite phase in the two hemispheres. We have analyzed SABER temperature measurements in the altitude range between 55 and 95 km to investigate the existence such variations. Because there are also strong tidal signatures (up to approximately 20 K) in the data, our algorithm estimates both mean values and tides together from the data. Based on SABER temperature data, the intra-annual variations with periods between 2 and 4 months can have amplitudes up to 5 K or more, depending on the altitude. Their amplitudes are in qualitative agreement with those inferred Erom UARS data (from different years). The SABER temperature variations also reveal pronounced hemispherical asymmetries, which are qualitatively consistent with wave driven meridional wind oscillations across the equator. Oscillations with similar periods have been seen in the meridional winds based on UARS data (Huang and Reber, 2003).

Huang, Frank T.↗

Single-Transducer, Ultrasonic Imaging Method for High-Temperature Structural Materials Eliminates the Effect of Thickness Variation in the Image

NASA Lewis Research Center's Life Prediction Branch, in partnership with Sonix, Inc., and Cleveland State University, recently advanced the development of, refined, and commercialized an advanced nondestructive evaluation (NDE) inspection method entitled the Single Transducer Thickness-Independent Ultrasonic Imaging Method. Selected by R&D Magazine as one of the 100 most technologically significant new products of 1996, the method uses a single transducer to eliminate the superimposing effects of thickness variation in the ultrasonic images of materials. As a result, any variation seen in the image is due solely to microstructural variation. This nondestructive method precisely and accurately characterizes material gradients (pore fraction, density, or chemical) that affect the uniformity of a material's physical performance (mechanical, thermal, or electrical). Advantages of the method over conventional ultrasonic imaging include (1) elimination of machining costs (for precision thickness control) during the quality control stages of material processing and development and (2) elimination of labor costs and subjectivity involved in further image processing and image interpretation. At NASA Lewis, the method has been used primarily for accurate inspections of high temperature structural materials including monolithic ceramics, metal matrix composites, and polymer matrix composites. Data were published this year for platelike samples, and current research is focusing on applying the method to tubular components. The initial publicity regarding the development of the method generated 150 requests for further information from a wide variety of institutions and individuals including the Federal Bureau of Investigation (FBI), Lockheed Martin Corporation, Rockwell International, Hewlett Packard Company, and Procter & Gamble Company. In addition, NASA has been solicited by the 3M Company and Allison Abrasives to use this method to inspect composite materials that are manufactured by these companies.

Roth, Don J.↗

Onset of Convection Due to Surface Tension Variations in Multicomponent and Binary Fluid Layers

Under certain conditions, such as in thin liquid films or microgravity, surface tension variations along a free surface can induce convection. Convection onset due to surface tension variation is important to many terrestrial technological processes in addition to microgravity materials processing applications. Examples include coating, drying crystallization, solidification, liquid surface contamination, and containerless processing. In double-diffusive and multicomponent systems, the spatial variations of surface tension are associated with two or more stratifying agencies, respectively. For example, both temperature and species (concentration) gradients are associated with convection in the solidification of binary alloys or salt ponds. The direction of the two (or more) gradients has a profound effect on the nature of the flow at or slightly beyond the onset of convection. Our recent work at the NASA Lewis Research Center focused on characterizing surface-tension-induced onset of convection, often referred to as Marangoni-Benard convection. Exact solutions for the stationary neutral stability of multicomponent fluid layers with interfacial deformation were derived. These solutions also permit the computation of a boundary curve that separates the long and finite wavelength instabilities. Computing points along this boundary using the exact solution (when possible) is more efficient than the typical numerical approaches, such as finite difference or spectral methods. Above the curve, a long wavelength instability was predicted, suggesting that convection would occur principally through one large flow cell in the layer, whereas below the curve, finite wavelength instabilities occur which suggest multiple finite-sized circulation cells. For many common liquids with layer depths greater than 100 mm, finite wave instability is predicted under terrestrial conditions; however, with little exception, long wavelength instability is predicted in microgravity for the identical fluid systems.

Skarda, J. Raymond Lee↗