Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Time-dependent experimental identification of inter-ELM microtearing modes in the tokamak edge on DIII-D

In a series of discharges on the DIII-D tokamak, fast vertical plasma jogs are used to induce current perturbations in the steep gradient region of the H-mode edge. These current perturbations directly influence the edge q profile, decoupling the resonant location and instability drive of pedestal-localized microtearing modes (MTMs.) By exploiting this effect, we develop and apply a new experimental technique to track the dynamical frequency evolution of MTMs in the pedestal region, providing a compelling validation of the MTM model. The frequency of potential MTMs is calculated as the Doppler-shifted electron diamagnetic frequency at rational q=m/n surfaces, showing remarkable agreement with chirped frequency behavior of n=3, 4 and 5 modes detected with fast magnetics. Data is collected throughout multiple ELM cycles in order to build robust statistics describing the time-dependent frequency evolution of MTMs, which can be explained by examining the recovery of pedestal gradients after an ELM event. MTMs have a dominant transport contribution in the electron thermal channel, so the presented results indicate that reduced models of pedestal transport must be electromagnetic in nature and constructed with accurate calculations of MTM stability; inclusion of this physics is essential for accurate predictions of the electron temperature pedestal profile. Supporting measurements of mode saturation, propagation direction and transport fingerprints are made to support the dynamic frequency determination, unambiguously and experimentally identifying MTMs in the pedestal region of DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Basal Melting and Oceanic Observations Beneath Central Fimbulisen, East Antarctica

Abstract Basal melting of ice shelves is fundamental to Antarctic ice sheet mass loss, yet direct observations remain sparse. We present the first year‐round melt record (2017–2021) from a phase‐sensitive radar on Fimbulisen, one of the fastest flowing ice shelves in Dronning Maud Land, East Antarctica. The observed long‐term mean ablation rate at 350 m depth below the central ice shelf was 1.0 ± 0.5 m yr −1 , marked by substantial sub‐weekly variability ranging from 0.4 to 3.5 m yr −1 . 36‐h filtered basal melt rate fluctuations closely align with ocean velocity. On seasonal time scales, melt rates peak during austral spring to autumn (September–March), driven by both elevated ocean velocities and thermal driving near the base. The combined effect of thermal driving and current speed explains the majority of the melt rate variability ( r = 0.84), highlighting the dominant role of shear‐driven turbulence. This relationship enables parameterization of melt rates for the decade‐long ocean record (2010–2021), although deviations appear under low and high forcing conditions. Both observed and parameterized melt rates show similar yearly mean magnitudes compared to satellite‐derived melt rates but with a tenfold lower seasonal amplitude and a 3‐month delay in seasonality. These detailed concurrent ice–ocean observations provide essential validation data for remote sensing and numerical models that aim to quantify and project ice‐shelf response to a change in ocean forcing. In situ measurements and continued monitoring are crucial for accurately assessing and modeling future basal melt rates, and for understanding the complex dynamics driving ice‐shelf stability and sea‐level change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Achieving designed texture and flows in bulk active nematics using optimal control theory

Being intrinsically nonequilibrium, active materials can potentially perform functions that would be thermodynamically forbidden in passive materials. However, active systems have diverse local attractors that correspond to distinct dynamical states, many of which exhibit chaotic turbulent-like dynamics and thus cannot perform work or useful functions. Designing such a system to choose a specific dynamical state is a formidable challenge. Motivated by recent advances enabling optogenetic control of experimental active materials, we describe an optimal control theory framework that identifies a spatiotemporal sequence of light-generated activity that drives an active nematic system toward a prescribed dynamical steady state. Active nematics are unstable to spontaneous defect proliferation and chaotic streaming dynamics in the absence of control. We demonstrate that optimal control theory can compute activity fields that redirect the dynamics into a variety of alternative dynamical programs and functions. This includes dynamically reconfiguring between states, selecting and stabilizing emergent behaviors that do not correspond to attractors, and are hence unstable in the uncontrolled system. Furthermore, our results provide a roadmap to leverage optical control methods to rationally design structure, dynamics, and function in a wide variety of active materials.

Complex systems theory↗
NASA NTRS1982

The measurement of the synoptic scale wind over the ocean

Mesoscale and microscale features of the turbulent winds over the ocean are related to the synoptic scale winds in terms of published spectral forms for the microscale, a mesoscale valley and published values of U*, VAR u', VAR v' and z/L, as defined in the text and as obtained for moderate to gale force winds. The frequencies involved correspond to periods longer than 1 hour and extend to the microscale, which starts at a period near 2 minutes, or so, and continues to the Kolmogorov inertial range. Nondimensional spectra that span both the mesoscale and the microscale are derived as a function of u, f(= n z/u) and z/L, where z is 10 meters, L is the Monin Obukov stability length and u is evaluated at 10 meters. For the same u, different values of z/L produce a range of values of u which in turn result in variations of the eddy structure of the mesoscale and microscale spectra. Both conventional anemometer averages and remotely sensed winds contain a random component of the mesoscale wind in their values. These components are differnces and not errors when winds are compared, and quantitative values for these differences are given. Ways to improve the measurement of the synoptic scale wind by transient ships, data buoys and scatterometers on future spacecraft are described. These ways are loner averaging times for ships and data buoys, depending on the synoptic conditions, and pooling spacecraft to form super observations. Design considerations for future remote sensing systems are given.

Pierson, W. J.↗
NASA NTRS2020

Assessment of Transition Modeling Capability in OVERFLOW with Emphasis on Swept-Wing Configurations

In preparation for comparisons with data obtained from the recently concluded experiments in the National Transonic Facility at the NASA Langley Research Center on the new common research model with natural laminar flow (CRM-NLF), an assessment of the current transition modeling capability in the NASA OVERFLOW 2.2o code has been carried out. A combination of the available experimental data and linear stability analysis is used to evaluate the accuracy and robustness of these models for selected swept-wing type configurations, with significant crossflow. An additional goal for this work involves providing a comparative assessment of the relevant transition models in the context of a single flow solver and identifying model limitations as well as the potential for future improvements that would help strengthen the physical basis of such transition models. Included in this investigation is an assessment of the sensitivities of the underlying transition models to grid resolution (wall-normal, as well as streamwise and spanwise) and the values of extra input parameters such as the level of surface roughness and freestream turbulence variables. The flow configurations targeted in this assessment include the NASA NLF(2)-0415 swept-wing configuration, the sickle-shaped wing introduced by the Technical University of Braunschweig, and the wing-body configuration of the CRM model from the fourth and fifth AIAA CFD Drag Prediction Workshops.

Freestream velocity↗
NASA NTRS1985

Wind-Tunnel Investigation of a Full-Scale Canard-Configured General Aviation Airplane

An investigation was conducted in the Langley 30- by 60-Foot Tunnel to determine the aerodynamic characteristics of a powered, full-scale model of a general aviation airplane employing a canard. Although primary emphasis of the investigation was placed on evaluating the aerodynamic performance and the stability and control characteristics of the basic configuration, tests were also conducted to study the following effects of varying the basic configuration: effect of Reynolds number; effect of canard; effect of outboard wing leading-edge droop; effect of center-of-gravity location; effect of elevator trim; effect of landing gear; effect of lateral-directional control; effect of power; effect of fixed transition; effect of water spray; effects of canard incidence, canard airfoil section, and canard position; and effects of winglets and upper winglet size. Additional aspects of the study were to determine the boundary-layer transition characteristics of airfoil surfaces and the effect of fixing the boundary layer to be turbulent by means of a transition strip near the leading edge. The tests were conducted at Reynolds numbers from 0.60 x 10 to the 6th power to 2.25x10 to the 6th power, based on the wing mean aerodynamic chord, at angles of attack from -4.5 deg to 41.5 deg, and at angles of sideslip from -15 deg to 15 deg.

Yip, L. P.↗
NASA NTRS1989

The use of direct numerical simulation in the study of turbulent, chemically-reacting flows

Full turbulence simulations are used here to study the effects of chemical heat release on the large-scale structures in turbulent mixing layers. In agreement with laboratory results, it is found that the heat release lowers the rate at which the mixing layer grows and reduces the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer produce changes in the flame vortex structure that act to produce more diffuse vortices than in the constant density case, resulting in lower rotation rates of the large-scale structures. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layer are shown to result from vorticity generation by baroclinic torques. Calculations of the energy in the various wavenumbers shows that the heat release has a stabilizing effect on the growth rate of individual modes. This methodology can be applied for Reynolds numbers less than several hundred and for Damkoehler numbers less than about ten.

Riley, J. J.↗
NASA NTRS2009

Smoke Point in Co-flow Experiment

The Smoke Point In Co-flow Experiment (SPICE) determines the point at which gas-jet flames (similar to a butane-lighter flame) begin to emit soot (dark carbonaceous particulate formed inside the flame) in microgravity. Studying a soot emitting flame is important in understanding the ability of fires to spread and in control of soot in practical combustion systems space. Previous experiments show that soot dominates the heat emitted from flames in normal gravity and microgravity fires. Control of this heat emission is critical for prevention of the spread of fires on Earth and in space for the design of efficient combustion systems (jet engines and power generation boilers). The onset of soot emission from small gas jet flames (similar to a butane-lighter flame) will be studied to provide a database that can be used to assess the interaction between fuel chemistry and flow conditions on soot formation. These results will be used to support combustion theories and to assess fire behavior in microgravity. The Smoke Point In Co-flow Experiment (SPICE) will lead to a o improved design of practical combustors through improved control of soot formation; o improved understanding of and ability to predict heat release, soot production and emission in microgravity fires; o improved flammability criteria for selection of materials for use in the next generation of spacecraft. The Smoke Point In Co-flow Experiment (SPICE) will continue the study of fundamental phenomena related to understanding the mechanisms controlling the stability and extinction of jet diffusion flames begun with the Laminar Soot Processes (LSP) on STS-94. SPICE will stabilize an enclosed laminar flame in a co-flowing oxidizer, measure the overall flame shape to validate the theoretical and numerical predictions, measure the flame stabilization heights, and measure the temperature field to verify flame structure predictions. SPICE will determine the laminar smoke point properties of non-buoyant jet diffusion flames (i.e., the properties of the largest laminar jet diffusion flames that do not emit soot) for several fuels under different nozzle diameter/co-flow velocity configurations. Luminous flame shape measurements would also be made to verify models of the flame shapes under co-flow conditions. The smoke point is a simple measurement that has been found useful to study the influence of flow and fuel properties on the sooting propensity of flames. This information would help support current understanding of soot processes in laminar flames and by analogy in turbulent flames of practical interest.

Urban, David L.↗
DOE OSTI2020

Nek5000 developments in support of industry and the NRC

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics verification and validation (V&V) work has focused in three areas of Nek5000 V&V-driven development. First, in a close collaborative effort with the U. S. Nuclear Regulatory Commission (NRC) staff, we have continued V&V efforts for the HYMERES-2 project using the OECD/NEA sponsored testing in the PSI PANDA facility. This year’s focus of ANL-NRC collaboration involves Nek5000 setups and validation for a range of problems relevant to and including the HYMERES-2 benchmark from PSI. The primary outcome of this year efforts is a more efficient geometry and inlet modeling simplification after a careful sensitivity study of the inlet profiles and pipe geometries. The resulting modeling choice of a short recycling/fully-developed turbulent inlet is within the experimental uncertainty estimate. This finding simplifies the next step of the cross-V&V HYMERES-2 project. In addition, the ANL team continue to provide assistance to the NRC staff in the form of Nek5000 application support in general and on the use of the HPC platforms of ALCF and INL in particular. This supports the NRC’s assessment of Nek5000 for use with the NRC Blue CRAB code suite. Second, we have implemented and tested more robust model of URANS, namely the k – τ model, a variant of the k-ω model, along with other improvements to RANS Nek5000 modeling in general. Because of its demonstrated robustness and stability, the k – τ model is the only RANS model that has been implemented in the new GPU version of the Nek5000 code, nekRS. Lastly, we report the initial implementation of Jacobian-free Newton Krylov approach to the direct Newton method for steady fluid solvers aimed at acceleration of RANS modeling and at IC improvement for LES campaigns. Also leveraging the Exascale Computing Project (ECP) ANL/CEED & SMR team’s software development effort to support NEAMS problems at large scale of the advanced computing architectures, NekRS, a GPU variant of Nek5000, built on top of kernels from libParanumal using OCCA for portability, has been successfully run on the full system of Summit (4608 nodes, 27648 GPUs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1980

Flow measurements in a model swirl combustor

Flow measurements in a model combustor composed of two confined coaxial swirling jets under noncombusting conditions are presented. Mean flow results are obtained for five different flow conditions to determine the effect of outer swirl on the recirculation zone (which is used for flame stabilization under combustion conditions). As the outer swirl magnitude is first decreased from maximum counter-swirl to zero and then increased to give co-swirl conditions, the size and the reverse flow velocity in the recirculation zone diminish. Detailed time mean and fluctuating flow measurements are obtained for a co-swirl and a counter-swirl condition with a directional pitot probe and hot-wire anenometry. For these two cases, recirculation zone occurs only with counter-swirl, near the exit of the inner jet. The recirculation zone is in the form of a one celled toroidal vortex having very low swirl velocity. Axial development of fluctuation levels, energy dissipation rates and turbulence length scales are described for the two flow conditions. Spectral analysis reveals periodic oscillations in both flows. The oscillations originate from the inner jet and up to 4 harmonics are observed. The fundamental frequencies are comparable to the rotational frequencies of the inner jets at the exit which are approximately under solid body rotation. The significance of the present results for the combustion process is discussed.

Vu, B. T.↗
DOE OSTI2020

Flow-induced oscillation patterns for two tandem cylinders with turbulence stimulation and variable stiffness and damping

Herein, the oscillation patterns of two 1-DOF cylinder oscillators, undergoing VIV and galloping, are investigated in a free surface water channel for 3.2 × 104 ≤ Re ≤ 1.2 × 10 5 . The cylinders are arranged in tandem and supported by springs for a range of different spring stiffness and damping parameters. The efficiency of a current energy converter (CEC), based on flow-induced oscillation (FIO) of multiple cylinders in tandem, is critically related to the cylinder oscillation patterns. Due to the limited number of studies on the FIO for multiple cylindrical or prismatic bodies, their oscillation patterns have neither been identified nor classified. The surfaces of the cylinders are modified by turbulence stimulation to enhance FIO. Three different center-to-center spacing, five stiffness, and six damping ratios, for a total of 90 sets of experiments were conducted. The current velocity range is from 0.34 m/s to 1.32 m/s. From more than 2000 tests, five major patterns, nine sub-patterns are identified and classified. The patterns are defined based on the amplitudes, frequencies and the phase angle differences between the two cylinders. The characteristics and mechanics of each oscillation pattern are explored and explained from the perspective of fluid-structure interaction (FSI). By systematically varying the parameters, the underlying hydrodynamic mechanisms, including the coupling level between vortices and cylinders, the shielding effect, the wake effect, and the stability states are revealed. Few preliminary observations on the connection between oscillation pattern and harnessed power by the tandem cylinders are reported.

42 ENGINEERING↗
DOE OSTI2024

Dependence of Cloud Macrophysical Properties and Phase Distributions on Environmental Conditions Over the North Atlantic and Southern Ocean: Results From COMBLE and MARCUS

The accurate representation of Cold Air Outbreaks (CAOs) and affiliated mixed-phase boundary layer (BL) clouds in models is challenging. How BL cloud properties evolve during CAOs and their dependence on meteorological conditions is not well understood but is important for the simulation of Earth's energy budgets. Here the properties of polar BL clouds over the North Atlantic (NA) and Southern Ocean (SO) are compared using observations from the Measurements of Aerosol Radiation and CloUds over the SO (MARCUS) and CAOs in the Marine BL Experiment (COMBLE) conducted over the NA. MARCUS observations show a stronger BL inversion than COMBLE, with a higher mean EIS (estimated inversion strength)/LTS (lower tropospheric stability) of -0.03 K/13 K compared to COMBLE’s -3.2 K/9.3 K. 39% of CAOs observed during COMBLE were intense with M > 5 K, while MARCUS only had 1.3%. 78%/72% of clouds sampled in CAOs during COMBLE/MARCUS had cloud top heights <4 km. The mean BL cloud top height was over 400 m higher, and the BL was over 500 m deeper for M of 10 K compared to 0 K for both regions. MARCUS observed a 27% moister BL structure than COMBLE when M > 5 K due to stronger BL inversion trapping more moisture within the BL. Under the same LTS, EIS, and M conditions, MARCUS observed a 12% drier BL structure, and clouds were 46% more turbulent than COMBLE. During CAOs, 54% of single-layer BL clouds sampled during MARCUS had liquid-dominated bases compared to 39% during COMBLE.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

On the stability of the ionopause of Venus

The stability of the Venus ionopause is examined in light of the importance of gravitation and curvature. Using a one-fluid approximation for the equation of motion of the plasma, and ignoring the effects of neutrals, a dispersion relation is obtained that includes the effects of the magnetic field, sheared plasma flow, buoyancy, centrifugal force and magnetic tension due to boundary curvature. It is found that buoyancy acts to neutralize the flute instability. As expected, the Kelvin-Helmholtz mode is the dominant instability over most of the dayside ionopause. The expected growth times of this mode are short in comparison with the wave-convection time over the boundary; the waves can grow and saturate quickly, producing a turbulent boundary that may affect electrodynamic coupling between the solar wind and ionospheric plasmas.

Elphic, R. C.↗
DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗
DOE OSTI2022

Similarity of the turbulent kinetic energy dissipation rate distribution in the upper mixed layer of the tropical Indian Ocean

Turbulence within the upper ocean mixed layer plays a key role in various physical, biological, and chemical processes. Between September and November 2011, a dataset of 570 vertical profiles of the turbulent kinetic energy (TKE) dissipation rate, as well as conventional hydrological and meteorological data, were collected in the upper layer of the tropical Indian Ocean. These data were used to statistically analyze the vertical distribution of the TKE dissipation rate in the mixed layer. The arithmetic-mean method made the statistical TKE dissipation rate profile more scattered than the median and geometric-mean methods. The statistical TKE dissipation rate were respectively scaled by the surface buoyancy flux and the TKE dissipation rate at the mixed-layer base. It was found that the TKE dissipation rate scaled by that at the mixed-layer base exhibited better similarity characteristics than that scaled by the surface buoyancy flux, whether the stability parameter D/|L MO | was greater or less than 10, indicating that the TKE dissipation rate at the mixed-layer base is a better characteristic scaling parameter for reflecting the intrinsic structure of the TKE dissipation rate in the mixed layer, where D and L MO are respectively the mixed-layer thickness and the Monin-Obukhov length scale. The parameterization of the TKE dissipation rate at the mixed-layer base on the shear-driven dissipation rate and the surface buoyancy flux was further explored. It was found that the TKE dissipation rate at the mixed-layer base could be well fitted by a linear combination of three terms: the wind-shear-driven dissipation rate, the surface buoyancy flux, and a simple nonlinear coupling term of these two.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Analysis of rotordynamic coefficients of helically-grooved turbulent annular seals

An analysis for helically-grooved turbulent annular seals is developed to predict leakage and dynamic coefficients, as related to rotordynamics. The grooved surface pattern is formulated as an inhomogeneous directivity in surface shear stress. The zeroth-order equations define the steady-state leakage and the circumferential velocity development due to wall shear for a centered rotor position. The first-order equations define perturbations in the pressure and axial and circumferential velocity fields due to small motion of the rotor about the centered position. Numerical results are presented for proposed grooved seals in the High Pressure Oxygen Turbopump (HPTOP) of the Space Shuttle Main Engine (SSME) and for a water-pump application. The results show that an optimum helix angle exists from a rotordynamic stability viewpoint. Further, a properly designed helically-grooved stator is predicted to have pronounced stability advantages over other currently used seals.

Kim, C.-H.↗
DOE OSTI2022

A new discovery of edge localized modes suppression using ICRH

Here, the high-confinement mode (H-mode) is very important for the fusion performance of ITER and future fusion reactors. The H-mode pedestal is a region featuring strongly reduced turbulence and transport just inside the limiting flux surface with strong plasma gradients, which drive edge localized modes (ELMs) in tokamaks. The ELMs, however, are quasiperiodic relaxations of the pedestal, resulting in a series of hot plasma eruptions similar to magnetospheric substorms, solar and stellar flares, which could potentially damage the ITER divertor plates and first walls. To mitigate and suppress the ELMs in H-mode plasmas, the fusion community has tried to develop effective actuators over the past two decades. Two separate excellent inventions have been awarded to recognize their achievements. The 2018 APS excellence award in Plasma Physics Research was given to Drs. Todd E. Evans, Max E. Fenstermacher, and Richard Alan Moyer for their first experimental demonstration of the stabilization of ELMs in high-confinement diverted discharges by application of very small edge-resonant magnetic perturbations, which led to the adoption of suppression coils in the ITER design. The 2018 Fusion Technology Award from the IEEE’s Nuclear and Plasma Science Society (NPSS) was given to Dr. Larry Baylor for his work designing the fueling, pumping and disruption mitigation system for the US ITER Project and for the implementation of pellet ELM mitigation for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Reynolds stress closure in jet flows using wave models

Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.

Morris, P. J.↗