Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time dependent fluctuations”

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 271 records · Page 15

Frequency-Dependent Viscosity of Xenon Near the Critical Point

We used a novel, overdamped oscillator aboard the Space Shuttle to measure the viscosity eta of xenon near its critical density rho(sub c), and temperature T(sub c). In microgravity, useful data were obtained within 0.1 mK of T(sub c), corresponding to a reduced temperature t = (T -T(sub c))/T(sub c) = 3 x 10(exp -7). The data extend two decades closer to T(sub c) than the best ground measurements, and they directly reveal the expected power-law behavior eta proportional to t(sup -(nu)z(sub eta)). Here nu is the correlation length exponent, and our result for the small viscosity exponent is z(sub eta) = 0.0690 +/- 0.0006. (All uncertainties are one standard uncertainty.) Our value for z(sub eta) depends only weakly on the form of the viscosity crossover function, and it agrees with the value 0.067 +/- 0.002 obtained from a recent two-loop perturbation expansion. The measurements spanned the frequency range 2 Hz less than or equal to f less than or equal to 12 Hz and revealed viscoelasticity when t less than or equal to 10(exp -1), further from T(sub c) than predicted. The viscoelasticity scales as Af(tau), where tau is the fluctuation-decay time. The fitted value of the viscoelastic time-scale parameter A is 2.0 +/- 0.3 times the result of a one-loop perturbation calculation. Near T(sub c), the xenon's calculated time constant for thermal diffusion exceeded days. Nevertheless, the viscosity results were independent of the xenon's temperature history, indicating that the density was kept near rho(sub c), by judicious choices of the temperature vs. time program. Deliberately bad choices led to large density inhomogeneities. At t greater than 10(exp -5), the xenon approached equilibrium much faster than expected, suggesting that convection driven by microgravity and by electric fields slowly stirred the sample.

Berg, Robert F.↗

Towards Full Aircraft Airframe Noise Prediction: Lattice Boltzmann Simulations

Computational results for an 18%-scale, semi-span Gulfstream aircraft model are presented. Exa Corporation's lattice Boltzmann PowerFLOW(trademark) solver was used to perform time-dependent simulations of the flow field associated with this high-fidelity aircraft model. The simulations were obtained for free-air at a Mach number of 0.2 with the flap deflected at 39 deg (landing configuration). We focused on accurately predicting the prominent noise sources at the flap tips and main landing gear for the two baseline configurations, namely, landing flap setting without and with gear deployed. Capitalizing on the inherently transient nature of the lattice Boltzmann formulation, the complex time-dependent flow features associated with the flap were resolved very accurately and efficiently. To properly simulate the noise sources over a broad frequency range, the tailored grid was very dense near the flap inboard and outboard tips. Extensive comparison of the computed time-averaged and unsteady surface pressures with wind tunnel measurements showed excellent agreement for the global aerodynamic characteristics and the local flow field at the flap inboard and outboard tips and the main landing gear. In particular, the computed fluctuating surface pressure field for the flap agreed well with the measurements in both amplitude and frequency content, indicating that the prominent airframe noise sources at the tips were captured successfully. Gear-flap interaction effects were remarkably well predicted and were shown to affect only the inboard flap tip, altering the steady and unsteady pressure fields in that region. The simulated farfield noise spectra for both baseline configurations, obtained using a Ffowcs-Williams and Hawkings acoustic analogy approach, were shown to be in close agreement with measured values.

Khorrami, Mehdi R.↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

JWST Noise Floor. II. Systematic Error Sources in JWST NIRCam Time Series

The James Webb Space Telescope (JWST) holds great promise for characterizing atmospheres of transiting exoplanets, potentially providing insights into Earth-sized planets within the habitable zones of M-dwarf host stars if photon-limited performance can be achieved. Here, we discuss the systematic error sources that are expected to be present in grism time-series observations with the NIRCam instrument. We find that pointing jitter and high-gain antenna moves in addition to the detectors' subpixel crosshatch patterns will produce relatively small variations (less than 6 parts per million, ppm). The time-dependent aperture losses due to thermal instabilities in the optics can also be kept to below 2 ppm. To achieve these low noise values, it is important to employ a sufficiently large (more than 1farcs1) extraction aperture. Persistence due to charge-trapping will have a minor (less than 3 ppm) effect on the time series 20 minutes into an exposure and is expected to play a much smaller role than it does for the Hubble Space Telescope WFC3 detectors. We expect detector temperature fluctuations to contribute less than 3 ppm. In total, our estimated noise floor from known systematic error sources is only 9 ppm per visit. We urge caution, however, because unknown systematic error sources could be present in flight and will only be measurable on astrophysical sources such as quiescent stars. We find that reciprocity failure may introduce a perennial instrument offset at the 40 ppm level, so corrections may be needed when a multi-instrument multi-observatory spectrum is stitched together over wide wavelength ranges.

Everett Schlawin↗

Aspects of the 40-50 day oscillation during the northern winter as inferred from outgoing longwave radiation

One-point correlation statistics and extended empirical orthogonal function analysis are used to show the existence of a principal mode of oscillation in the intraseasonal variability (ISV) of tropical convection. The data analyzed were derived from NOAA Outgoing Longwave Radiation data covering seven 6 mo periods from 1974-83. The data is filtered to remove fluctuations on time scales under 5 days and by subtracting a least-squares parabola to eliminate seasonal means. The previously-known 40-50 day oscillation is found to be the most prevalent component of the ISV of tropical convection. An E-W dipole forms and propagates over the equatorial Indian Ocean, then stalls and intensifies over Indonesia and the equatorial Indian Ocean, then stalls and intensifies over Indonesia and the equatorial central Pacific. Various anomalies generated in association with the dipole and dependent on its exact location when stalled are discussed.

Lau, K.-M.↗

Flux Transfer Event in the Subsolar Region and Near the Cusp: Simultaneous Polar and Cluster Observations

The phenomenon called flux transfer events (FTEs) is widely accepted as the manifestation of time-dependent reconnection. In this paper, we present an observational evidence of a flux transfer event observed simultaneously at low-latitude by Polar and high-latitude by Cluster. This event occurred on March 21, 2002, when both Cluster and Polar were located near the local noon but with large latitudinal distance. Cluster was moving outbound from polar cusp to the magnetosheath, and Polar was in the magnetosheath near the equatorial magnetopause. The observations show that a flux transfer event was formed between the equator and the northern cusp. Polar and Cluster observed the FTE's two open flux tubes: Polar saw the southward moving flux tube near the equator; and Cluster the , northward moving flux tube at high latitude. Unlike low-latitude FTEs, the high-latitude FTE did not exhibit the characteristic bi-polar BN signature. But the plasma data clearly showed its open flux tube configuration. Enhanced electric field fluctuations were observed within the FTE core, both at low- and high-attitudes. This event provides us a unique opportunity to understand high-latitude FTE signatures and the nature of time-varying reconnection.

Le, G.↗

Time dependent heat transfer rates in high Reynolds number hypersonic flowfields

Time dependent heat transfer rates have been calculated from time dependent temperature measurements in the vicinity of shock-wave boundary-layer interactions due to conical compression ramps on an axisymmetric body. The basic model is a cylindrical body with a 10 degree conical nose. Four conical ramps, 20, 25, 30, and 35 degrees serve as shock wave generators. Flowfield surveys have been made in the vicinity of the conical ramp vertex, the separation point, and the reattachment point. A significant effort was made to characterize the natural frequencies and relative powers of the resulting fluctuations in heat transfer rates. This research effort, sponsored jointly by NASA and the Air Force, was conducted in the Air Force Flight Dynamics Directorate High Reynolds Facility. The nominal freestream Mach number was 6, and the freestream Reynolds numbers ranged from 2.2 million/ft to 30.0 million/ft. Experimental results quantify temperature response and the resulting heat transfer rates as a function of ramp angle and Reynolds number. The temperature response within the flowfield appears to be steady-state for all compression ramp angles and all Reynolds numbers, and hence, the heat transfer rates appear to be steady-state.

Flanagan, Michael J.↗

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗

Turbulent correlation measurements in a two-stream mixing layer.

Two point space-time measurements of the axial component of fluctuating velocity were made using linearized hot-wire anemometry. Space scales, convected frame correlations and time scales, and convection velocities in the shear layer were evaluated. Both filtered narrow and broad band convected frame autocorrelations were determined. Differences between axial broad band convection velocities and both mean turbulence and mean entrained fluid velocities were observed. Scaled broad band convection velocities for the mixing layer and simple round jets were found to collapse to a common curve. Axial narrow band convection velocities showed strong frequency and transverse position dependence.

Jones, B. G.↗

Rapid fluctuations in solar flares

Study of rapid fluctuations in the emission of radiation from solar flares provides a promising approach for probing the magneto-plasma structure and plasma processes that are responsible for a flare. It is proposed that elementary flare bursts in X-ray and microwave emission may be attributed to fine structure of the coronal magnetic field, related to the aggregation of photospheric magnetic field into magnetic knots. Fluctuations that occur on a subsecond time-scale may be due to magnetic islands that develop in current sheets during magnetic reconnection. The impulsive phase may sometimes represent the superposition of a large number of the elementary energy-release processes responsible for elementary flare bursts. If so, the challenge of trying to explain the properties of the impulsive phase in terms of the properties of the elementary processes must be faced. Magnetic field configurations that might produce solar flares are divided into a number of categories, depending on: whether or not there is a filament; whether there is no current sheet, a closed current sheet, or an open current sheet; and whether the filament erupts into the corona, or is ejected completely from the Sun's atmosphere. Analysis of the properties of these possible configurations is compared with different types of flares, and to Bai's subdivision of gamma-ray/proton events.

Sturrock, Peter A.↗

Fluid Physics in a Fluctuating Acceleration Environment

Our program of research aims at developing a stochastic description of the residual acceleration field onboard spacecraft (g-jitter) to describe in quantitative detail its effect on fluid motion. Our main premise is that such a statistical description is necessary in those cases in which the characteristic time scales of the process under investigation are long compared with the correlation time of g-jitter. Although a clear separation between time scales makes this approach feasible, there remain several difficulties of practical nature: (i), g-jitter time series are not statistically stationary but rather show definite dependences on factors such as active or rest crew periods; (ii), it is very difficult to extract reliably the low frequency range of the power spectrum of the acceleration field. This range controls the magnitude of diffusive processes; and (iii), models used to date are Gaussian, but there is evidence that large amplitude disturbances occur much more frequently than a Gaussian distribution would predict. The lack of stationarity does not constitute a severe limitation in practice, since the intensity of the stochastic components changes very slowly during space missions (perhaps over times of the order of hours). A separate analysis of large amplitude disturbances has not been undertaken yet, but it does not seem difficult a priori to devise models that may describe this range better than a Gaussian distribution. The effect of low frequency components, on the other hand, is more difficult to ascertain, partly due to the difficulty associated with measuring them, and partly because they may be indistinguishable from slowly changing averages. This latter effect is further complicated by the lack of statistical stationarity of the time series. Recent work has focused on the effect of stochastic modulation on the onset of oscillatory instabilities as an example of resonant interaction between the driving acceleration and normal modes of the system, and on cavity flow as an example of how an oscillatory response under periodic driving becomes diffusive if the forcing is random instead. This paper describes three different topics that illustrate behavior that is peculiar to a stochastic acceleration field. In the first case, we show that g-jitter can induce effective attractive or repulsive forces between a pair of spherical particles that are suspended in an incompressible fluid of different density provided that the momentum diffusion length is larger than the interparticle separation (as in the case in most colloidal suspensions). Second, a stochastic modulation of the control parameter in the vicinity of a pitchfork or supercritical bifurcation is known not to affect the location of the threshold. We show, however, that resonance between the modulation and linearly stable modes close to onset can lead to a shift in threshold. Finally, we discuss the classical problem of vorticity diffusion away from a plane boundary that is being vibrated along its own plane. Periodic motion with zero average vorticity production results in an exponential decay of the vorticity away from the boundary. Random vibration, on the other hand, results in power law decay away from the boundary even if vorticity production averages to zero.

Drolet, Francois↗

Effect of magnetic field geometry on the wave signature of the pickup of interstellar neutrals

By using the technique of Johnstone et al. (1991), the asymptotic wave spectrum induced by the pickup of interstellar ions in the heliosphere is derived. This approach is independent of quasi-linear theory. The asymptotic wave spectrum is calculated from the standard resonance condition by assuming that the asymptotic ion distribution is bispherical and by using conservation of energy. The calculation places no restriction on either the ratio of Alfven speed to solar wind flow speed or the particle pitch angle. Spectra for various geometries (relative orientations of the magnetic field and solar wind velocity) are calculated. Absolute upper limits to the expected wave enhancements are calculated. These upper limits are compared to the background fluctuation spectrum extrapolated from 0.87 AU, assuming a spectral index of -5/3 and radial dependences of -2 and -3. The heliocentric radial distances at which the peaks in the various geometry-dependent wave spectra are expected to be observable above the background are identified. Peaks due to pickup hydrogen at the proton gyrofrequency are found to be always stronger relative to the background in parallel than in perpendicular geometries. The parallel-geometry excitations consist only of sunward propagating waves. In perpendicular geometries a peak also occurs, this time near 0.1 of the proton gyrofrequency (corresponding to the peak calculated by Lee and Ip (1987)). However, the steep spectral dependence of the ambient wave field makes it unlikely that this peak will be observed inside 20 AU. A similar result holds for the case of pickup-helium induced excitation at the helium gyrofrequency. Inside 10 AU, peaks at these lower frequencies tend to be obscured by the steep frequency dependence (proportional to -5/3) of the ambient spectrum and then are only observable in parallel geometries.

Williams, L. L.↗

Propagation velocity and space-time correlation of perturbations in turbulent channel flow

A database obtained from direct numerical simulation of a turbulent channel flow is analyzed to extract the propagation velocity V of velocity, vorticity, and pressure fluctuations from their space-time correlations. A surprising result is that V is approximately the same as the local mean velocity for most of the channel, except for the near-wall region. For y(+) is less than or equal to 15, V is virtually constant, implying that perturbations of all flow variables propagate like waves near the wall. In this region V is 55 percent of the centerline velocity U(sub c) for velocity and vorticity perturbations and 75 percent of U(sub c) for pressure perturbations. Scale-dependence of V is also examined by analyzing the bandpass filtered flow fields. Comprehensive documentation of the propagation velocities and space-time correlation data, which should prove useful in the evaluation of Taylor's hypothesis is presented. An attempt was made to explain some of the data in terms of our current understanding of organized structures, although not all of the data can be explained this way.

Kim, John↗

Is walking a random walk? Evidence for long-range correlations in stride interval of human gait

Complex fluctuation of unknown origin appear in the normal gait pattern. These fluctuations might be described as being (1) uncorrelated white noise, (2) short-range correlations, or (3) long-range correlations with power-law scaling. To test these possibilities, the stride interval of 10 healthy young men was measured as they walked for 9 min at their usual rate. From these time series we calculated scaling indexes by using a modified random walk analysis and power spectral analysis. Both indexes indicated the presence of long-range self-similar correlations extending over hundreds of steps; the stride interval at any time depended on the stride interval at remote previous times, and this dependence decayed in a scale-free (fractallike) power-law fashion. These scaling indexes were significantly different from those obtained after random shuffling of the original time series, indicating the importance of the sequential ordering of the stride interval. We demonstrate that conventional models of gait generation fail to reproduce the observed scaling behavior and introduce a new type of central pattern generator model that sucessfully accounts for the experimentally observed long-range correlations.

Hausdorff, Jeffrey M.↗

Detection of atmospheric pressure loading using very long baseline interferometry measurements

Loading of the Earth by the temporal redistribution of global atmospheric mass is likely to displace the positions of geodetic monuments by tens of millimeters both vertically and horizontally. Estimates of these displacements are determined by convolving National Meteorological Center (NMC) global values of atmospheric surface pressure with Farrell's elastic Green's functions. An analysis of the distances between radio telescopes determined by very long baseline interferometry (VLBI) between 1984 and 1992 reveals that in many of the cases studied there is a significant contribution to baseline length change due to atmospheric pressure loading. Our analysis covers intersite distances of between 1000 and 10,000 km and is restricted to those baselines measured more than 100 times. Accounting for the load effects (after first removing a best fit slope) reduces the weighted root-mean-square (WRMS) scatter of the baseline length residuals on 11 of the 22 baselines investigated. The slight degradation observed in the WRMS scatter on the remaining baselines is largely consistent with the expected statistical fluctuations when a small correction is applied to a data set having a much larger random noise. The results from all baselines are consistent with approximately 60% of the computed pressure contribution being present in the VLBI length determinations. Site dependent coefficients determined by fitting local pressure to the theoretical radial displacement are found to reproduce the deformation caused by the regional pressure to within 25% for most inland sites. The coefficients are less reliable at near coastal and island stations.

Vandam, T. M.↗

Variances of the components and magnitude of the polar heliospheric magnetic field

The heliolatitude dependences of the variances in the components and the magnitude of the heliospheric magnetic field have been analysed, using the Ulysses magnetic field observations from close to the ecliptic plane to 80 southern solar latitude. The normalized variances in the components of the field increased significantly (by a factor about 5) as Ulysses entered the purely polar flows from the southern coronal hole. At the same time, there was at most a small increase in the variance of the field magnitude. The analysis of the different components indicates that the power in the fluctuations is not isotropically distributed: most of the power is in the components of the field transverse to the radial direction. Examining the variances calculated over different time scales from minutes to hours shows that the anisotropy of the field variances is different on different scales, indicating the influence of the two distinct populations of fluctuations in the polar solar wind which have been previously identified. We discuss these results in terms of evolutionary, dynamic processes as a function of heliocentric distance and as a function of the large scale geometry of the magnetic field associated with the polar coronal hole.

Balogh, A.↗

Study of the Unsteady Flow Features on a Stalled Wing

The occurrence of large scale structures in the post stall flow over a rectangular wing at high angles of attack was investigated in a small-scale subsonic wind tunnel. Mean and time dependent measurements within the separated flow field suggest the existence of two distinct angle of attack regimes beyond wing stall. The shallow stall regime occurs over a narrow range of incidence angles (2-3 deg.) immediately following the inception of leading edge separation. In this regime, the principal mean flow structures, termed stall cells, are manifested as a distinct spanwise periodicity in the chordwise extent of the separated region on the model surface with possible lateral mobility not previously reported. Within the stall cells and on the wing surface, large amplitude pressure fluctuations occur with a frequency much lower than anticipated for bluff body shedding, and with minimum effect in the far wake. In the deep stall regime, stall cells are not observed and the separated region near the model is relatively free of large amplitude pressure disturbances.

Yon, Steven A.↗

Dissipation of Turbulence in the Solar Wind as Measured by Cluster

Turbulence in fluids and plasmas is a scale-dependent process that generates fluctuations towards ever-smaller scales until dissipation occurs. Recent Cluster observations in the solar wind demonstrate the existence of a cascade of magnetic energy from the scale of the proton Larmor radius, where kinetic properties of ions invalidate fluid approximations, down to the electron Larmor radius, where electrons become demagnetized. The cascade is quasi-two-dimensional and has been interpreted as consisting of highly oblique kinetic Alfvenic fluctuations that dissipate near at the electron gyroradius scale via proton and electron Landau damping. Here we investigate for the first time the spatial properties of the turbulence at these scales. We report the presence of thin current sheets and discontinuities with spatial sizes greater than or approximately equal to the proton Larmor radius. These isolated structures may be manifestations of intermittency, and such would localize sites of turbulent dissipation. Studying the relationship between turbulent dissipation, reconnection and intermittency is crucial for understanding the dynamics of laboratory and astrophysical plasmas.

Goldstein, Melvyn↗