Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic theory”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Effect of Porosity on Surface Catalytic Efficiency

This paper describes the effect of surface porosity of thermal protection materials on surface catalytic efficiency using test data taken from both arc-jet and side-arm reactor facilities. Relative surface porosity of the samples varied from 6% to 50%. Surface porosity was measured using a flow apparatus and Bernoulli equation. The surface catalytic efficiency of the materials was calculated using aerothermodynamic, and kinetic theories. The catalytic efficiency of the materials are compared at surface temperatures between room temperature and 2500 F. The data are presented in the form of graphs and tables.

Stewart, David A.↗
NASA NTRS2002

Structural Fluctuations and Thermophysical Properties of Molten II-VI Compounds

The objectives of the project are to conduct ground-based experimental and theoretical research on the structural fluctuations and thermophysical properties of molten II-VI compounds to enhance the basic understanding of the existing flight experiments in microgravity materials science programs as well as to study the fundamental heterophase fluctuation phenomena in these melts by: 1) conducting neutron scattering analysis and measuring quantitatively the relevant thermophysical properties of the II-VI melts (such as viscosity, electrical conductivity, thermal diffusivity and density) as well as the relaxation characteristics of these properties to advance the understanding of the structural properties and the relaxation phenomena in these melts and 2) performing theoretical analyses on the melt systems to interpret the experimental results. All the facilities required for the experimental measurements have been procured, installed and tested. Thermal diffusivity of molten tellurium has been measured by a laser flash method in the temperature range of 500 C to 900 C. The measured diffusivity as a function of temperature agrees fairly well with published data. However, a relaxation phenomenon, which shows a slow drift of the measured thermal conductivity toward the equilibrium value after cooling of the melt, was observed for the first time. An apparatus based on the transient torque induced by a rotating magnetic field has been developed to determine the viscosity and electrical conductivity of semiconducting liquids. Viscosity measurements on molten tellurium showed a similar relaxation behavior to the measured diffusivity. The density and volume expansion coefficients for pure Te and HgTe melts were measured as a function of temperature using a pycnometric method. A density maximum was found for both melts but no relaxation behavior was observed. Neutron scattering experiments were performed on the HgTe and HgZnTe melts and the results on pair distribution showed better resolution than previously reported. A simple kinetic theory, which contains the formation reactions of Te polymerization that increases the molecular weight, was proposed to interpret the relaxation behavior of HgZnTe viscosity data.

Su, Ching-Hua↗
NASA NTRS2003

Microscopic Modeling of Intersubband Optical Processes in Type II Semiconductor Quantum Wells: Linear Absorption

Linear absorption spectra arising from intersubband transitions in semiconductor quantum well heterostructures are analyzed using quantum kinetic theory by treating correlations to the first order within Hartree-Fock approximation. The resulting intersubband semiconductor Bloch equations take into account extrinsic dephasing contributions, carrier-longitudinal optical phonon interaction and carrier-interface roughness interaction which is considered with Ando s theory. As input for resonance lineshape calculation, a spurious-states-free 8-band kp Hamiltonian is used, in conjunction with the envelop function approximation, to compute self-consistently the energy subband structure of electrons in type II InAs/AlSb single quantum well structures. We demonstrate the interplay of nonparabolicity and many-body effects in the mid-infrared frequency range for such heterostructures.

Li, Jian-Zhong↗
NASA NTRS1997

A quantitative assay for intercellular aggregation

In an earlier communication (Munn et al., J Immunol. Methods 166: 11-25, 1993), we presented the initial development of a quantitative assay for monitoring the rates of cellular aggregation based on digital image processing and video microscopy. This study describes some important enhancements and modifications to the procedure. A new index is introduced to characterize the three-dimensional morphology of the aggregates. This index is based on temporal changes in the projected area of the cells and cell aggregates during the course of the experiment. By drawing an analogy with the kinetic theory of gases, we have also introduced a procedure to normalize for variations in cell seeding density among different experiments. In addition, the image analysis technique has been improved by introducing a background subtraction algorithm to remove illumination defects and an adaptive segmentation procedure. These improvements allowed us to completely automate the image analysis procedure, thus minimizing user intervention and improving the reproducibility of the measurements. The enhanced visual assay is evaluated using some recent results from our studies on homotypic lymphocyte aggregation.

Non-NASA Center↗
NASA NTRS2004

Two Formal Gas Models For Multi-Agent Sweeping and Obstacle Avoidance

The task addressed here is a dynamic search through a bounded region, while avoiding multiple large obstacles, such as buildings. In the case of limited sensors and communication, maintaining spatial coverage - especially after passing the obstacles - is a challenging problem. Here, we investigate two physics-based approaches to solving this task with multiple simulated mobile robots, one based on artificial forces and the other based on the kinetic theory of gases. The desired behavior is achieved with both methods, and a comparison is made between them. Because both approaches are physics-based, formal assurances about the multi-robot behavior are straightforward, and are included in the paper.

Kerr, Wesley↗
NASA NTRS2004

Time-Average Measurement of Velocity, Density, Temperature, and Turbulence Using Molecular Rayleigh Scattering

Measurement of time-averaged velocity, density, temperature, and turbulence in gas flows using a nonintrusive, point-wise measurement technique based on molecular Rayleigh scattering is discussed. Subsonic and supersonic flows in a 25.4-mm diameter free jet facility were studied. The developed instrumentation utilizes a Fabry-Perot interferometer to spectrally resolve molecularly scattered light from a laser beam passed through a gas flow. The spectrum of the scattered light contains information about velocity, density, and temperature of the gas. The technique uses a slow scan, low noise 16-bit depth CCD camera to record images of the fringes formed by Rayleigh scattered light passing through the interferometer. A kinetic theory model of the Rayleigh scattered light is used in a nonlinear least squares fitting routine to estimate the unknown parameters from the fringe images. The ability to extract turbulence information from the fringe image data proved to be a challenge since the fringe is broadened by not only turbulence, but also thermal fluctuations and aperture effects from collecting light over a range of scattering angles. Figure 1 illustrates broadening of a Rayleigh spectrum typical of flow conditions observed in this work due to aperture effects and turbulence for a scattering angle, chi(sub s), of 90 degrees, f/3.67 collection optics, mean flow velocity, u(sub k), of 300 m/s, and turbulent velocity fluctuations, sigma (sub uk), of 55 m/s. The greatest difficulty in processing the image data was decoupling the thermal and turbulence broadening in the spectrum. To aid in this endeavor, it was necessary to seed the ambient air with smoke and dust particulates; taking advantage of the turbulence broadening in the Mie scattering component of the spectrum of the collected light (not shown in the figure). The primary jet flow was not seeded due to the difficulty of the task. For measurement points lacking particles, velocity, density, and temperature information could reliably be recovered, however the turbulence estimates contained significant uncertainty. Resulting flow parameter estimates are presented for surveys of Mach 0.6, 0.95, and 1.4 jet flows. Velocity, density, and temperature were determined with accuracies of 5 m/s, 1.5%, and 1%, respectively, in flows with no particles present, and with accuracies of 5 m/s, 1-4%, and 2% in flows with particles. Comparison with hotwire data for the Mach 0.6 condition demonstrated turbulence estimates with accuracies of about 5 m/s outside the jet core where Mie scattering from dust/smoke particulates aided in the estimation of turbulence. Turbulence estimates could not be recovered with any significant accuracy for measurement points where no particles were present.

Mielke, Amy F.↗
NASA NTRS2004

Time-Average Molecular Rayleigh Scattering Technique for Measurement of Velocity, Denisty, Temperature, and Turbulence Intensity in High Speed Nozzle Flows

A molecular Rayleigh scattering based flow diagnostic is developed to measure time average velocity, density, temperature, and turbulence intensity in a 25.4-mm diameter nozzle free jet facility. The spectrum of the Rayleigh scattered light is analyzed using a Fabry-Perot interferometer operated in the static imaging mode. The resulting fringe pattern containing spectral information of the scattered light is recorded using a low noise CCD camera. Nonlinear least squares analysis of the fringe pattern using a kinetic theory model of the Rayleigh scattered light provides estimates of density, velocity, temperature, and turbulence intensity of the gas flow. Resulting flow parameter estimates are presented for an axial scan of subsonic flow at Mach 0.95 for comparison with previously acquired pitot tube data, and axial scans of supersonic flow in an underexpanded screeching jet. The issues related to obtaining accurate turbulence intensity measurements using this technique are discussed.

Mielke, Amy F.↗
NASA NTRS2013

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species With Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for nearequilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion of the heating and is then compared to the total heating measured in flight.

Liechty, Derek S.↗
NASA NTRS2014

Extension of a Kinetic Approach to Chemical Reactions to Electronic Energy Levels and Reactions Involving Charged Species with Application to DSMC Simulations

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. Recently introduced molecular-level chemistry models that predict equilibrium and nonequilibrium reaction rates using only kinetic theory and fundamental molecular properties are extended in the current work to include electronic energy level transitions and reactions involving charged particles. These extensions are shown to agree favorably with reported transition and reaction rates from the literature for near-equilibrium conditions. Also, the extensions are applied to the second flight of the Project FIRE flight experiment at 1634 seconds with a Knudsen number of 0.001 at an altitude of 76.4 km. In order to accomplish this, NASA's direct simulation Monte Carlo code DAC was rewritten to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced chemistry model, and to include the extensions presented in this work. The 1634 second data point was chosen for comparisons to be made in order to include a CFD solution. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid because, although near-transitional, the flow is still considered to be continuum. It is shown that the inclusion of electronic energy levels in the DSMC simulation is necessary for flows of this nature and is required for comparison to the CFD solution. The flow field solutions are also post-processed by the nonequilibrium radiation code HARA to compute the radiative portion.

Liechty, Derek S.↗
NASA NTRS2014

Mass Dependency of Isotope Fractionation of Gases Under Thermal Gradient and Its Possible Implications for Planetary Atmosphere Escaping Process

Physical processes that unmix elements/isotopes of gas molecules involve phase changes, diffusion (chemical or thermal), effusion and gravitational settling. Some of those play significant roles for the evolution of chemical and isotopic compositions of gases in planetary bodies which lead to better understanding of surface paleoclimatic conditions, e.g. gas bubbles in Antarctic ice, and planetary evolution, e.g. the solar-wind erosion induced gas escaping from exosphere on terrestrial planets.. A mass dependent relationship is always expected for the kinetic isotope fractionations during these simple physical processes, according to the kinetic theory of gases by Chapman, Enskog and others [3-5]. For O-bearing (O16, -O17, -O18) molecules the alpha O-17/ alpha O-18 is expected at 0.5 to 0.515, and for S-bearing (S32,-S33. -S34, -S36) molecules, the alpha S-33/ αlpha S-34 is expected at 0.5 to 0.508, where alpha is the isotope fractionation factor associated with unmixing processes. Thus, one isotope pair is generally proxied to yield all the information for the physical history of the gases. However, we recently] reported the violation of mass law for isotope fractionation among isotope pairs of multiple isotope system during gas diffusion or convection under thermal gradient (Thermal Gradient Induced Non-Mass Dependent effect, TGI-NMD). The mechanism(s) that is responsible to such striking observation remains unanswered. In our past studies, we investigated polyatomic molecules, O2 and SF6, and we suggested that nuclear spin effect could be responsible to the observed NMD effect in a way of changing diffusion coefficients of certain molecules, owing to the fact of negligible delta S-36 anomaly for SF6.. On the other hand, our results also showed that for both diffusion and convection under thermal gradient, this NMD effect is increased by lower gas pressure, bigger temperature gradient and lower average temperature, which indicate that the nuclear spin effect may not be the significant contributor as the energies involved in the hyperfine effect are much smaller than those with molecular collisions, especially under convective conditions.

Sun, Tao↗
NASA NTRS2015

A Unified Computational Model for Solar and Stellar Flares

We present a unified computational framework that can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker-Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

radiative transfer↗
NASA NTRS2017

Preliminary Analysis of Two Years of the Massive Collision Monitoring Activity

It is hypothesized that the interactions between many of the most massive derelicts in low Earth orbit are more frequent than modeled by the traditional combination of kinetic theory of gases and Poisson probability distribution function. This is suggested by the fact that there are clusters of derelicts where members' inclinations are nearly identical and their apogees/perigees overlap significantly resulting in periodic synchronization of the objects' orbits. In order to address this proposition, an experiment was designed and conducted over the last two years. Results from this monitoring and characterization experiment are presented with implications for proposed debris remediation strategies. Four separate clusters of massive derelicts were examined that are centered around 775km, 850km, 975km, and 1500km, respectively. In aggregate, the constituents of these clusters contain around 500 objects and about 800,000kg of mass; this equates to a third of all derelict mass in LEO. Preliminary analysis indicates that encounter rates over this time period for these objects are greater than is estimated by traditional techniques. Hypothesized dependencies between latitude of encounter, relative velocity, frequency of encounters, inclination, and differential semi-major axis were established and verified. This experiment also identified specific repeatable cluster dynamics that may reduce the cost/risk and enhance the effectiveness of debris remediation activities and also enable new operational debris remediation options.

McKnight, Darren↗
NASA NTRSDate not supplied

CFD Modeling of Phase Change and Pressure Drop during Violent Sloshing of Cryogenic Fluid in a Small-Scale Tank

This paper presents the development of a two-phase CFD model that is used to study sloshing of a cryogenic fluid (LN2) undergoing phase change generated by lateral acceleration and its effects on the ensuing heat and mass transfer and pressure drop in a small-scale tank. In this model the interface is captured using the volume-of-fluid (VOF) method [1]. Kinetic-theory based Schrage relation [2] is used for calculating phase change mass transfer at the liquid-vapor interface. Computational results are compared to the data provided by a non-isothermal sloshing experiment with phase change conducted by the University of Tokyo and JAXA in 2018 using liquid Nitrogen in a transparent tank. The effect of different thicknesses of a liquid temperature stratification layer created at the interface on the pressure drop during sloshing was studied experimentally. CFD results for interface movement and tank pressure are presented and compared in this paper to experimental data for the cases with the smallest and largest liquid temperature stratification thickness.

phase change↗
NASA NTRS2021

Physics and Chemistry of Te and HgTe-based Ternary Melts

Historically, the theoretical treatment of the liquid phases has always been more difficult and complicated than that for solid and gas phases. A liquid has no lattice structure as crystalline solids and the atoms/molecules in the liquid can migrate through it relatively rapidly. On the other hand, it is also interacting with many other atoms/molecules so that the simplifications of the kinetic theory of gases cannot be employed. For more complicated liquids, such as the liquids of high ionicity and those containing hydrogen bonds and electric dipoles, the understanding is far from complete. At the same time, accurate information on the physics and chemistry of semiconductor melts is needed for the quantitative descriptions of the process of crystal growth from melt. The pre-crystallization phenomena in the liquid phase are critical because the properties of the grown crystals depend on the state and structure of the melt as well as the thermal history of the melt during solidification process. However, the data on the liquid phase, such as thermophysical properties of semiconductor melts are scarce, especially for the HgTe-based II-VI ternary compound semiconductors because of their high vapor pressure and extreme toxicity. Analysis of the thermophysical properties of the melt can provide information about structural transitions of the melt during the solidification process. From a broader point of view, the structure of liquids is much more complicated than the crystalline solids, especially the relaxation behavior through different thermal histories. The theory of hetero-phase fluctuations of liquids is applicable to any many-body systems including condensed-matter physics, field theory, physics of nuclear-matter, cosmology, biology and even sociology. This book summarizes the physics and chemistry from the experimental measurements and theoretical analyses of phase diagram, thermodynamic properties, density, thermal conductivity, viscosity, and electrical conductivity on the binary, pseudo-binary and ternary melts of the most advanced IR-detector material systems of HgCdTe and HgZnTe as well as the analyses of these results. The main objectives of this study are: (1) to provide the phase diagrams and thermodynamic properties of Hg-Cd-Te and Hg-Zn-Te systems through quantitatively fitting the experimental data by assuming an associated solution model for the liquid phase, (2) to experimentally measure the thermophysical properties of the Hg-Cd-Te and Hg-Zn-Te melts, including density, viscosity, electrical conductivity and thermal conductivity as functions of temperature and composition and (3) to enhance the fundamental knowledge of hetero-phase fluctuations and relaxation phenomena in the melts and extend our understanding of the solidification process in order to interpret the experimental results of crystal growth so as to improve the melt growth processes of the compound semiconductor. The physics and chemistry of Te and HgTe-based ternary melts were explored through the studies of the structural transformation during melting, the supercooling during solidification, the relaxation phenomena after rapid cooling of the melts and the metal-semiconductor transition in the melts through the analyses of electrical conductivity and Lorenz number. An in-depth study on the thermophysical properties and their time-dependent structural dynamic processes taking place in the vicinity of the solid-liquid phase transition of the narrow homogeneity range HgTe-based ternary semiconductors as well as the analysis of the homogenization process in the melt will also be presented.

Thermophysical properties↗
NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗
NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗