Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “contact loss”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Containerless, Low-Gravity Undercooling of Ti-Ce Alloys in the MSFC Drop Tube

Previous tests of the classical nucleation theory as applied to liquid-liquid gap miscibility systems found a discrepancy between experiment and theory in the ability to undercool one of the liquids before the L1-L2 separation occurs. To model the initial separation process in a two-phase liquid mixture, different theoretical approaches, such as free-energy gradient and density gradient theories, have been put forth. If there is a large enough interaction between the critical liquid and the crucible, both models predict a wetting temperature (T(sub w)) above which the minority liquid perfectly wets and layers the crucible interface, but only on one side of the immiscibility dome. Materials with compositions on the other side of the dome will have simple surface adsorption by the minority liquid before bulk separation occurs when the coexistence (i.e., binoidal) line in reached. If the interaction between the critical liquid and the crucible were to decrease, T(sub w) would increase, eventually approaching the critical consolute temperature (T(sub cc)). If this situation occurs, then there could be large regions of the miscibility gap where non-perfect wetting conditions prevail resulting in droplets of L1 liquid at the surface having a non-zero contact angle. The resulting bulk structure will then depend on what happens on the surface and the subsequent processing conditions. In the past several decades, many experiments in space have been performed on liquid metal binary immiscible systems for the purpose of determining the effects that different crucibles may have on the wetting and separation process of the liquids. Potard performed experiments that showed different crucible materials could cause the majority phase to preferentially wet the container and thus produce a dispersed microstructure of the minority phase. Several other studies have been performed on immiscibles in a semi-container environment using an emulsion technique. Only one previous study was performed using completely containerless processing of immiscible metals and the results of that investigation are similar to some of the emulsion studies. In all the studies, surface wetting was attributed as the cause for the similar microstructures or the asymmetry in the ability to undercool the liquid below the binoidal on one side of the immiscibility dome. By removing the container completely from the separation process, it was proposed that the loss of the crucible/liquid interaction would produce a large shift in T(sub w) and thus change the wetting characteristics at the surface. By investigating various compositions across the miscibility gap, a change in the type and amount of liquid wetting at the surface of a containerless droplet should change the surface nucleating behavior of the droplet - whether it be the liquid-liquid wetting or the liquid-to-solid transition. Undercooling of the liquid into the metastable region should produce significant differences in the separation process and the microstructure upon solidification. In this study, we attempt to measure these transitions by monitoring the temperature of the sample by optical pyrometry. Microstructural analysis will be made to correlate with the degree of undercooling and the separation mechanisms involved.

Robinson, M. B.↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Sputter-Coated TiO 2 Films as Passivation and Hole Transfer Layers for Improved Energy Conversion with Solar Fuel WO 3 /CuWO 4 Photoanodes

Atomic layer deposited (ALD) “leaky” TiO 2 have gained interest as charge-selective protection layers for semiconductor solar fuel electrodes. Here we demonstrate the use of sputter-deposited TiO 2 layers as hole selective contacts for WO 3 /CuWO 4 type 2 heterojunction water oxidation photoanodes for the first time. TiO 2 protection layers with varying thicknesses (2 to 128 nm) were deposited using the RF magnetron sputtering technique. The resulting TiO 2 films are amorphous based on Raman spectroscopy and powder XRD. Photoelectrochemical scans and Vibrating Kelvin probe photovoltage spectroscopy show that 2-8 nm TiO 2 layers nearly double the photocurrent to 0.97 mA cm -2 under AM 1.5 illumination (19% AQE at 350 nm), increase the surface photovoltage signal by 25%, and increase the WO 3 /CuWO 4 bandgap. These effects can be attributed to the selectivity of TiO 2 for photoholes. Additionally, SPV data suggest that TiO 2 overlayers suppress copper-based surface recombination defects. Reduced photocurrent and the photovoltage are seen in thicker TiO 2 films (16 to 128 nm) as a result of an increasing hole transfer resistance and because of light shading effects according to photoaction spectra. The TiO 2 films also improve the stability of the WO 3 /CuWO 4 photoelectrodes, allowing nearly constant O 2 evolution over 3 hours after an initial 20-35% loss. Overall, this work establishes RF magnetron sputtering as a useful method to install amorphous TiO 2 passivation layers for improved WO 3 /CuWO 4 solar fuel photoelectrodes. Furthermore, we show how the combination of PEC with SPV measurements provides insight into the function of the TiO 2 coatings.

CuWO4↗
NASA NTRS1962

Titanium-Oxygen Reactivity Study

A program has been conducted at Astronautics to investigate the likelihood of occurrence of the catastrophic oxidation of titanium alloy sheet under conditions which simulate certain cases of accidental failure of the metal while it is in contact with liquid or gaseous oxygen. Three methods of fracturing the metal were used; they consisted of mechanical puncture, tensile fracture of welded joints, and perforation by very high velocity particles. The results of the tests which have been conducted provide further evidence of the reactivity of titanium with liquid and gaseous oxygen. The evidence indicates that the rapid fracturing of titanium sheet while it is in contact with oxygen initiates the catastrophic oxidation reaction. Initiation occurred when the speed of the fracture was some few feet per second, as in both the drop-weight puncture tests and the static tensile fracture tests of welded joints, as well as when the speed was several thousand feet per second, as in the simulated micrometeoroid penetration tests. The slow propagation of a crack, however, did not initiate the reaction. It may logically be concluded that the localized frictional heat of rapid fracture and/or spontaneous oxidation (exothermic) of minute particles emanating from the fracture cause initiation of the reaction. Under conditions of slow fracture, however, the small heat generated may be adequately dissipated and the reaction is not initiated. A portion of the study conducted consisted of investigating various means by which the reaction might be retarded or prevented. Providing a "barrier" at the titanium-oxygen interface consisting of either aluminum metal or a coating of a petroleum base corrosion inhibitor appeared to be only partially effective in retarding the reaction. The accidental puncturing or similar rupturing of thin-walled pressurized oxygen tanks on missiles and space vehicle will usually constitute loss of function, and may sometimes cause their catastrophic destruction by explosive decompression regardless of the type of material used for their construction. In the case of tanks constructed of titanium alloys the added risk is incurred of catastrophic burning of the tanks. In view of this it is recommended that thin-walled tanks constructed of titanium alloys should not be used to contain liquid or gaseous oxygen.

Chafey, J. E.↗
NASA NTRS2011

NASA Tech Briefs, February 2011

Topics covered include: Multi-Segment Radius Measurement Using an Absolute Distance Meter Through a Null Assembly; Fiber-Optic Magnetic-Field-Strength Measurement System for Lightning Detection; Photocatalytic Active Radiation Measurements and Use; Computer Generated Hologram System for Wavefront Measurement System Calibration; Non-Contact Thermal Properties Measurement with Low-Power Laser and IR Camera System; SpaceCube 2.0: An Advanced Hybrid Onboard Data Processor; CMOS Imager Has Better Cross-Talk and Full-Well Performance; High-Performance Wireless Telemetry; Telemetry-Based Ranging; JWST Wavefront Control Toolbox; Java Image I/O for VICAR, PDS, and ISIS; X-Band Acquisition Aid Software; Antimicrobial-Coated Granules for Disinfecting Water; Range 7 Scanner Integration with PaR Robot Scanning System; Methods of Antimicrobial Coating of Diverse Materials; High-Operating-Temperature Barrier Infrared Detector with Tailorable Cutoff Wavelength; A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane; Thermally Conductive Tape Based on Carbon Nanotube Arrays; Two Catalysts for Selective Oxidation of Contaminant Gases; Nanoscale Metal Oxide Semiconductors for Gas Sensing; Lightweight, Ultra-High-Temperature, CMC-Lined Carbon/Carbon Structures; Sample Acquisition and Handling System from a Remote Platform; Improved Rare-Earth Emitter Hollow Cathode; High-Temperature Smart Structures for Engine Noise Reduction and Performance Enhancement; Cryogenic Scan Mechanism for Fourier Transform Spectrometer; Piezoelectric Rotary Tube Motor; Thermoelectric Energy Conversion Technology for High-Altitude Airships; Combustor Computations for CO2-Neutral Aviation; Use of Dynamic Distortion to Predict and Alleviate Loss of Control; Cycle Time Reduction in Trapped Mercury Ion Atomic Frequency Standards; and A (201)Hg+ Comagnetometer for (199)Hg+ Trapped Ion Space Atomic Clocks.

Source record↗
NASA NTRS2023

Magnetic Anisotropy and Stress Dependent Epoxy Wetting in FeNi-based Metal Amorphous Nanocomposite

Recently developed FeNi-based metal amorphous nanocomposites (MANCs) used in high-speed motors (HSMs) exhibit reduced eddy current losses while maintaining good mechanical properties and glass-forming abilities. Magnetic anisotropy in (Fe 70 Ni 30 ) 80 Nb 4 B 14 Si 2 amorphous magnetic ribbon (AMR) in the as-cast state and upon conventional (T CA ) and strain (T SA ) annealing heat treatment is investigated. From ribbon samples in as-cast condition, quenched in stress from planar flow casting (PFC) induced as-cast curvature derived uniaxial magnetic anisotropy. Stress relief by conventional furnace annealing at T CA > 350 °C achieved isotropic properties in the bulk. Annealing about the primary crystallization temperature, T CA ~ 450 °C, resulted in the formation of both FCC and BCC nanocrystallites and evolution to isotropic bulk magnetic properties confirming the random anisotropy model. In samples strain annealed at T SA = 440 °C at various tensions, relatively large controlled induced uniaxial anisotropy is achieved. The largest magnetic anisotropy occurs in annealing under the stress of 250 MPa yielding an anisotropy field of 7.1 x 10 4 A/m. Surface anisotropy observed by MOKE differs from bulk anisotropy due to image contrast from closure domains. Epoxy coatings are important for improved bonding, mechanical properties, and resistivity in tape-wound MANC cores for HSMs. Using a sessile droplet method, the equilibrium contact angle of an epoxy droplet on a tensile stress annealed MANC exhibits stress dependent surface energies. Anisotropic wetting in FeNi-based MANC heat treated at T CA = 440 °C mimics surface magnetic anisotropy observed by MOKE.

soft magnetic materials↗
NASA NTRS2003

Estimating the Tropospheric Ozone Distribution by the Assimilation of Satellite Data

Tropospheric ozone is important to the environment, because it acts as a strong oxidant to control the concentrations of many reduced gases (methane, carbon monoxide, ... ), its radiative forcing plays a significant role in the greenhouse effect, and direct contact with ozone is harmful to human health. Tropospheric ozone, whose main sources are intrusion from the stratosphere and chemical production from source gases associated with urban pollution or biomass burning, varies on a wide range of spatial and temporal scales. Its transport and chemistry can be influenced by weather, seasonal, or multiannual variability. Despite the importance of tropospheric ozone, it contributes only about 10% of the total ozone loading in the atmosphere. Consequently, satellite instruments lose sensitivity below the stratospheric ozone peak, and provide little information about middle and lower tropospheric ozone. This talk will discuss recent modifications made to the satellite ozone data assimilation system at NASA's Data Assimilation Office (DAO) in order to provide better tropospheric ozone columns and profiles. We use a version of the system that assimilates only the data from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument. The quality of the assimilated ozone in the tropical troposphere is evaluated by comparison with independent observations obtained from the Southern Hemispheric Additional Ozonesondes (SHADOZ) network. It is shown that the quality of ozone fields is sensitive to the winds used in the transport model. Increasing the vertical resolution of the model also has a beneficial impact. The assimilated ozone in the lower troposphere was substantially improved by inclusion of tropospheric ozone production, loss, and dry deposition rates from the Harvard GEOS-CHEM model. The mechanisms behind these results will be examined and the implications for our understanding of tropospheric ozone will be discussed.

Hayashi, Hiroo↗
DOE OSTI2024

Catalytic Reduction of Esters over Zirconia-Supported Metal Catalysts

Esters are often produced as unwanted byproducts during the catalytic upgrading of ethanol to diesel fuel precursors through Guerbet coupling. Removal of esters from the product stream is important to prevent the loss of downstream catalyst activity from ester-derived carboxylic acids. In this work, we studied ester hydrogenolysis to the parent alcohols as a viable route for enhanced diesel fuel production. Specifically, we investigated the reduction of hexyl acetate in butanol over ZrO 2 -supported Ni, Co, Cu, Rh, Pd, and Pt catalysts, where Cu/ZrO 2 was the most selective catalyst for the hydrogenolysis of hexyl acetate into hexanol and ethanol. Thermodynamic analysis reveals that a 90% alcohol yield can be obtained at 200 °C, 30 bar, and a relatively high H 2 :hexyl acetate molar ratio of 480:1. Experimentally, an alcohol yield of 88% yield was obtained with a 10 wt % Cu/ZrO 2 catalyst at these conditions with a residence time of 5.4 h kg cat kmol gas –1 . Catalytic tests on the support revealed that ZrO 2 catalyzes the transesterification reaction between hexyl acetate and butanol. However, only the Cu sites can catalyze the hydrogenolysis of the esters into the final alcohols. We developed a kinetic model for our experimental results, which shows that the transesterification and hydrogenolysis reactions run at two different timescales, the former being 10 times faster than the latter. Data regression has been used to develop a model to predict the mole fraction distribution of ester hydrogenolysis products over a wide range of contact times. Cu/ZrO 2 loses half its catalytic activity after 80 h of time on stream. Modeling of deactivation data reveals that the ZrO 2 support conserves a residual activity due to external active sites, while active sites over the Cu surface deactivate at different rates. Furthermore, the catalytic conversion of esters into their parent alcohols is relevant to the production of surrogate liquid fuels since alcohols can be bimolecularly dehydrated to produce a blend of ethers with diesel fuel-like properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2021

Power Operations of the Mars Exploration Rovers

The rovers of the National Aeronautics and Space Administration’s (NASA) Mars Exploration Rovers (MER) project, Spirit (MER-A) and Opportunity (MER-B), safely landed on the surface of Mars three weeks apart during January 2004. Spirit and Opportunity were built and operated by the Jet Propulsion Laboratory (JPL), which is managed by the California Institute of Technology (Caltech) for NASA. Spirit landed at Gusev Crater, 14.8 degrees south of the equator, and operated continuously on the surface of Mars from January 4, 2004 until last contact from Spirit on March 22, 2010. Opportunity landed at Meridiani Planum, a location 2.5 degrees south of the Martian equator. Opportunity operated continuously on the Martian surface from January 25, 2004 until the last received transmission from Opportunity on June 10, 2018. The goal of the MER project was to determine if Mars ever had a habitable environment, in particular, if it ever had water. During their missions, both Spirit and Opportunity found evidence that liquid water once flowed on the surface of Mars. Both Spirit and Opportunity used a 1.33 m2 triple-junction solar array as their power sources. Based on observations of the original Mars rover, the Sojourner rover of the Mars Pathfinder mission that landed on Mars on July 4, 1997, the expectation was that the Martian dust would rapidly accumulate on the solar arrays of Spirit and Opportunity, and that the rovers would not have enough energy to continue operations after 90 Martian days (sols). Instead, due in part to lower dust accumulation rates than expected and numerous dust cleaning events, the Spirit and Opportunity rovers continued to operate on the Martian surface for over 2000 sols (MER-A) and 5000 sols (MER-B), respectively. During this time, the rovers experienced multiple Martian winters and several dust storms. Because the sources of solar array energy loss were known, the solar array energy output offered a method to scientifically estimate the loading and aeolian removal of dust from the solar arrays each sol. The MER Power operations team called this value that they calculated the solar array Dust Factor (DF). Dust Factor was defined as the fraction of sunlight that penetrates the accumulated dust on the surface of the solar array. A Dust Factor of 1.0 would indicate that the solar array was perfectly clean. A Dust Factor of 0.6 would indicate that only 60% of the available sunlight was able to penetrate the accumulated dust on the solar arrays. The MER Power subsystem operations team used the Multi-Mission Power Analysis Tool (MMPAT) to perform these Dust Factor calculations. The MMPAT software tool modeled the behavior of the solar arrays and the batteries as they interacted with the spacecraft power loads over the mission timeline. MMPAT also had knowledge (through telemetry and user inputs) of telemetered Power subsystem voltages and currents, atmospheric opacity (Tau), rover surface location, rover attitude, the planetary tilt and distance of Mars from the sun based on day of year, the instantaneous elevation of the sun based on time of day, temperatures (internal and external), terrain masking, and shadowing (due to the camera mast and antennas). Once all of the known sources of array energy loss are accounted, the remaining difference between the expected array energy and the actual array energy determines the solar array Dust Factor. The assumptions made while determining the Dust Factor are 1) that the single measured atmospheric opacity value (Tau) is constant over the course of the entire sol, 2) there is no measurable solar cell degradation, 3) there are no shorted strings and therefore, 4) all unexpected solar array energy losses are due to accumulated dust on the solar array. Although it cannot provide an absolute measure of dust loading, the determination of solar array Dust Factor provides a useful way of tracking dust accumulation and aeolian dust removal trends on Martian spacecraft. Any spacecraft on the Martian surface is vulnerable to dust, especially solar-powered spacecraft. The Spirit and Opportunity rovers were operational on the Martian surface for much longer than expected due in part to aeolian removal of dust from their solar arrays. The first few dust removal events were a pleasant surprise to the MER operations teams; however, over time a pattern began to arise. In over three Mars Years on the surface, the Power operations team tracked the solar array Dust Factor at Gusev Crater (location of MER-A) and observed that there were several significant dust removal and deposition events in Mars Year (MY) 27, even in the absence of a large dust storm. In MY 28 at Gusev Crater, the large atmospheric opacity (Tau) increase lagged significant dust removal events. Overall at Gusev Crater, there was a pattern of steady dust accumulation on the solar arrays, with a small number of significant dust cleaning events. At Meridiani Planum, where Opportunity rover operated for over seven Mars Years (late MY 26 to mid MY 34), a clear and consistent pattern of dust movement emerged. Meridiani Planum had a predictable, seasonally dependent pattern of gradual and continuous dust accumulation and removal. In summary, this paper explains the reasons for the development of the solar array Dust Factor and how it was used in mission operations. In particular, this paper describes the MMPAT software package, how it models array energy, including the important assumptions, model inputs and sources of error. And finally, this paper will show how the calculated solar array Dust Factor was used at Gusev Crater and Meridiani Planum to predict dust accumulation rates and weather patterns, and the importance of this generated data set for current and future solar- powered missions to Mars, such as the InSight lander and the planned Mars Sample Return rover.

Chin, Keith B.↗
DOE OSTI2023

Modeling Advanced Strippers for CO2 Capture from Gas-Fired Power Plants using Aqueous Piperazine

Energy performance of 5 m piperazine was evaluated with an advanced stripper at the pilot scale, accounting for heat loss. Modeling studies indicated that heat loss at locations other than the heat source impacted the heat duty. For strippers with excess packing, the column was the most important source of heat loss, and values as low as 0.03 GJ/hr can cause pinch and reabsorption of CO2. Solvents like 5 m PZ are more affected by heat loss in the column due to a top-side temperature pinch at high lean loading. Heat loss impacts NGCC CO2 capture more than coal-based CO2 capture. Surface temperature measurements pilot plants showed that heat loss was 50-75% controlled by natural convection. Measured heat loss was correlated with steam flowrate and wind speed and ranged from 8126 to 219668 Btu/hr, representing 35% of the measured heat rate on average. Net heat duty linearly varied with CO2 removal at low lean loading and was similar for coal and NGCC conditions at fixed removal. At 90% removal, net heat duty was about 2.5 GJ/t with 4 and 12% CO2 in the flue gas. 7 At NGCC conditions, strippers with finite packing can benefit from a flashing feed to the top of the column by using a hot bypass. A flashing feed was linked to a reduction in irreversibility of the stripper from temperature driving forces. It also reduces heat duty by up to 6% at fixed packing height and reduces packing requirement for a fixed steam heater size compared to a warm bypass. Lean loading and lean solvent rate were effective handles to maximize profitability of a fixed stripper design at low gas price. High ambient temperature operation can benefit from low pressure stripping to about 0.18 lean loading at 150 ℃ to maintain the cyclic capacity at reduced rich loading. At low gas price, the capture plant was able to maximize profitability even with a high heat duty of 2.7-3.1 GJ/t at a lean loading of 0.2-0.22 mol/mol. Heat recovery by partial water vapor condensation in a gas-liquid exchanger can be replaced by a direct contact condenser (DCC). The DCC improved gas cooling and reduced heat duty compared to the CO2 exchanger at 0.2 lean loading. The DCC when used with a 150 ℃/5.5 bar stripper can reduce the cost of capture compared to the base case by $3- 4/tonne at 15 ft of packing (optimum) but increases cost of capture by $6/tonne with a 120 ℃/2 bar stripper and 10 ft of packing (optimum). The DCC can work effectively with a solvent with a high heat of absorption and thermal stability.

Suresh Babu, Athreya↗
DOE OSTI2022

Local conformations and heterogeneities in structures and dynamics of isotactic polypropylene adsorbed onto carbon fiber

Carbon fiber (CF) reinforced polymers (CFRPs) have experienced widespread use in various industries. One of the most important parameters that controls the macroscopic property of CFRPs is the interface between a polymer matrix and CF. There is growing evidence to suggest the formation of a bound polymer layer (BPL), i.e., polymer chains that physically adsorb on a filler surface. However, this interface is always in contact with the thicker part of a polymer matrix, rendering its understanding a difficult task. Therein, we use CF-reinforced isotactic polypropylene (iPP) as a rational CFRP. To characterize the BPL on the CF surface, we extracted it from the CFRP using solvent-rinsing with p-xylene. The physical and thermal properties of the BPL were characterized by differential scanning calorimetry and thermogravimetric analysis, while its microscopic structures and dynamics were probed by small-angle neutron scattering and quasi-elastic neutron scattering (QENS) techniques. Subsequently, we employed atomistic molecular dynamics (MD) simulations to complement the QENS results above the bulk melting temperature and reveal details that were experimentally inaccessible. In this work, we observed that the degree of crystallinity of the BPL was quite lower than the bulk, while the melting temperature of the BPL remained the same as the bulks. Within the given length and time scales probed by QENS and MD, we also observed that most of the bound chains were mobile, with the formation of a high-density region (less than 1 nm in thickness) near the CF surface. The segmental dynamics of the bound chains probed by both QENS and MD were also much faster than those of the free chains, possibly due to the presence of a free surface region at the topmost surface of the BPL. Furthermore, the MD results demonstrated that the backbone chains and side groups lie nearly flat on the CF surface, which is the driving force for the flattening process of the iPP BPL to overcome the conformational entropy loss in the total free energy.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗
DOE OSTI2023

CCN data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.5 LPM is directed to the CCN. To calculate particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate particle loss through the inlet during sampling. Second, the corrected SMPS data was used to calculate the average of the total losses per scan down the CCN line. Then, the correction was applied to the CCN data. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. The supersaturation uncertainty is estimated conservatively at +/- 0.03%, where variation in the inlet temperature, pressure, and calibration technique prevents a more accurate measurement. Confidence in the reported supersaturation measurements is based on a pre-campaign calibration (following the methods from our previous work and Deng, 2014 based on Rose, 2008) in addition to inter-comparisons with the DOE for two days (22/08/18 and 22/09/01) where TAMU was co-located with AMF1. The inter-comparisons show good agreement between our instrument and the DOEs instrument on both days at all supersaturations. After the last inter-comparison on 22/09/01, there was no indication of a malfunction by our instrument through the rest of the campaign. Unfortunately, the instrument was dropped during demobilization. A post-campaign calibration was conducted, which showed a substantial departure from the pre-campaign calibration. The drop may have damaged the instrument’s ability to produce the desired supersaturations. Therefore, we do consider the data after 22/09/01 to be correct, but it should be used with caution. The CCN counter sampled for 3 minutes at each supersaturation setpoint (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%). At the end of a cycle, the instrument was set to 0.01% supersaturation for 5 minutes. The data is comprised of the last 60 seconds of each supersaturation set point (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%) to ensure the instrument stabilized and was able to reach thermal equilibrium. We removed the data during the periods where there were operational difficulties, setup, or maintenance. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Rose et. al. Calibration and Measurement Uncertainties of a Continuous-Flow Cloud Condensation Nuclei Counter (DMT-CCNC): CCN Activation of Ammonium Sulfate and Sodium Chloride Aerosol Particles in Theory and Experiment. Atmos. Chem. Phys., 8, 1153-1179, 2008. Deng et. al. Using Raman Microspectroscopy to Determine Chemical Composition and Mixing State of Airborne Marine Aerosols over the Pacific Ocean. Aerosol Science and Technology, Vol 48, Issue 2, 2014. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

The circum-Chryse region as a possible example of a hydrologic cycle on Mars: Geologic observations and theoretical evaluation

The transection and superposition relationships among channels, chaos, surface materials units, and other features in the circum-Chryse region of Mars were used to evaluate relative age relationships and evolution of flood events. Channels and chaos in contact (with one another) were treated as single discrete flood-carved systems. Some outflow channel systems form networks and are inferred to have been created by multiple flood events. Within some outflow channel networks, several separate individual channel systems can be traced to a specific chaos which acted as flood-source area to that specific flood channel. Individual flood-carved systems were related to widespread materials units or other surface features that served as stratigraphic horizons. Chryse outflow channels are inferred to have formed over most of the perceivable history of Mars. Outflow channels are inferred to become younger with increasing proximity to the Chryse basin. In addition, outflow channels closer to the basin show a greater diversity in age. The relationship of subsequent outflow channel sources to the sources of earlier floods is inferred to disfavor episodic flooding due to the progressive tapping of a juvenile near-surface water supply. Instead, we propose the circum-Chryse region as a candidate site of past hydrological recycling. The discharge rates necessary to carve the circum-Chryse outflow channels would have inevitably formed temporary standing bodies of H2O on the Martian surface where the flood-waters stagnated and pooled (the Chryse basin is topographically enclosed). These observations and inferences have led us to formulate and evaluate two hypotheses: (1) large amounts of the sublimated H2O off the Chryse basin flood lakes precipitated (snowed) onto the flood-source highlands and this H2O was incorporated into the near surface, recharging the H2O sources, making possible subsequent deluges; and (2) ponded flood-water in Chryse basin drained back down an anti basinward dipping subsurface layer accessed along the southern edge of the lake, recharging the flood-source aquifers. H2O not redeposited in the flood-source region was largely lost to the hydrologic cycle. This loss progressively lowered the vitality of the cycle, probably by now killing it. Our numerical evaluations indicate that of the two hypotheses we formulated, the groundwater seep cycle seems by far the more viable. Optimally, approx. 3/4 of the original mass of an ice-covered cylindrical lake (albedo 0.5, 1 km deep, 100-km radius, draining along its rim for one quarter of its circumference into substrata with a permeability of 3000 darcies) can be modeled to have moved underground (on timescales of the order of 10(exp 3) years) before the competing mechanisms of sublimation and freeze down choked off further water removal. Once underground, this water can travel distances equal to the separation between Chryse basin and flood-source sites in geologically short (approx. 10(exp 6) year-scale) times. Conversely, we calculate that optimally only approx. 40% of the H2O carried from Chryse can condense at the highlands, and most of the precipitate would either collect at the base of the highlands/lowlands scarp or sublimate at rates greater than it would accumulate over the flood-source sites. Further observations from forthcoming missions may permit the determination of which mechanisms may have operated to recycle the Chryse flood-waters.

Moore, Jeffrey M.↗
NASA NTRS1995

The Circum-Chryse Region as a Possible Example of a Hydrologic Cycle on Mars: Geologic Observations and Theoretical Evaluation

The transection and superposition relationships among channels, chaos, surface materials units, and other features in the circum-Chryse region of Mars were used to evaluate relative age relationships and evolution of flood events. Channels and chaos in contact (with one another) were treated as single discrete flood-carved systems. Some outflow channel systems form networks and are inferred to have been created by multiple flood events. Within some outflow channel networks, several separate individual channel systems can be traced to a specific chaos which acted as flood-source area to that specific flood channel. Individual flood-carved systems were related to widespread materials units or other surface features that served as stratigraphic horizons. Chryse outflow channels are inferred to have formed over most of the perceivable history of Mars. Outflow channels are inferred to become younger with increasing proximity to the Chryse basin. In addition, outflow channels closer to the basin show a greater diversity in age. The relationship of subsequent outflow channel sources to the sources of earlier floods is inferred to disfavor episodic flooding due to the progressive tapping of a juvenile near-surface water supply. Instead, we propose the circum-Chryse region as a candidate site of past hydrological recycling. The discharge rates necessary to carve the circum-Chryse outflow channels would have inevitably formed temporary standing bodies of H2O on the Martian surface where the flood-waters stagnated and pooled (the Chryse basin is topographically enclosed). These observations and inferences have led us to formulate and evaluate two hypotheses: (1) large amounts of the sublimated H2O off the Chryse basin flood lakes precipitated (snowed) onto the flood-source highlands and this H2O was incorporated into the near surface, recharging the H2O sources, making possible subsequent deluges; and (2) ponded flood-water in Chryse basin drained back down an anti basinward dipping subsurface layer accessed along the southern edge of the lake, recharging the flood-source aquifers. H2O not redeposited in the flood-source region was largely lost to the hydrologic cycle. This loss progressively lowered the vitality of the cycle, probably by now killing it. Our numerical evaluations indicate that of the two hypotheses we formulated, the groundwater seep cycle seems by far the more viable. Optimally, approximately 3/4 of the original mass of an ice-covered cylindrical lake (albedo 0.5, 1 km deep, 100-km radius, draining along its rim for one quarter of its circumference into substrata with a permeability of 3000 darcies) can be modeled to have moved underground (on timescales of the order of 10(exp 3) years) before the competing mechanisms of sublimation and freeze down choked off further water removal. Once underground, this water can travel distances equal to the separation between Chryse basin and flood-source sites in geologically short (approximately 10(exp 6) year-scale) times. Conversely, we calculate that optimally only approximately 40% of the H2O carried from Chryse can condense at the highlands, and most of the precipitate would either collect at the base of the highlands/lowlands scrap or sublimate at rates greater than it would accumulate over the flood-source sites. Further observations from forthcoming missions may permit the determination of which mechanisms may have operated to recycle the Chryse flood-waters.

Moore, Jeffrey M.↗
NASA NTRS2009

Assessment of the Short-Term Radiometric Stability between Terra MODIS and Landsat 7 ETM+ Sensors

The Landsat 7 (L7) Enhanced Thematic Mapper (ETM+) sensor was launched on April 15th, 1999 and has been in operation for over nine years. It has six reflective solar spectral bands located in the visible and shortwave infrared part of the electromagnetic spectrum (0.5 - 2.5 micron) at a spatial resolution of 30 m. The on-board calibrators are used to monitor the on-orbit sensor system changes. The ETM+ performs solar calibrations using on-board Full Aperture Solar Calibrator (FASC) and the Partial Aperture Solar Calibrator (PASC). The Internal Calibrator Lamp (IC) lamps, a blackbody and shutter optics constitute the on-orbit calibration mechanism for ETM+. On 31 May 2003, a malfunction of the scan-line corrector (SLC) mirror assembly resulted in the loss of approximately 22% of the normal scene area. The missing data affects most of the image with scan gaps varying in width from one pixel or less near the centre of the image to 14 pixels along the east and west edges of the image, creating a wedge-shaped pattern. However, the SLC failure has no impacts on the radiometric performance of the valid pixels. On December 18, 1999, the Moderate Resolution Imaging Spectroradiometer (MODIS) Proto-Flight Model (PFM) was launched on-board the NASA's EOS Terra spacecraft. Terra MODIS has 36 spectral bands with wavelengths ranging from 0.41 to 14.5 micron and collects data over a wide field of view angle (+/-55 deg) at three nadir spatial resolutions of 250 m, 500 in 1 km for bands 1 to 2, 3 to 7, and 8 to 36, respectively. It has 20 reflective solar bands (RSB) with spectral wavelengths from 0.41 to 2.1 micron. The RSB radiometric calibration is performed by using on-board solar diffuser (SD), solar diffuser stability monitor (SDSM), space-view (SV), and spectro-radiometric calibration assembly (SRCA). Through the SV port, periodic lunar observations are used to track radiometric response changes at different angles of incidence (AOI) of the scan mirror. As a part of the AM Constellation satellites, Terra MODIS flies approximately 30 minutes behind L7 ETM+ in the same orbit. The orbit of L7 is repetitive, circular, sunsynchronous, and near polar at a nominal altitude of 705 km (438 miles) at the Equator. The spacecraft crosses the Equator from north to south on a descending node between 10:00 AM and 10:15 AM. Circling the Earth at 7.5 km/sec, each orbit takes nearly 99 minutes. The spacecraft completes just over 14 orbits per day, covering the entire Earth between 81 degrees north and south latitude every 16 days. The longest continuous imaging swath that L7 sensor can collect is for a 14-minute subinterval contact period which is equivalent to 35 full WRS-2 scenes. On the other hand, Terra can provide the entire corresponding orbit with wider swath at any given ETM+ collection without contact time limitation. There are six spectral matching band pairs between MODIS (bands 3, 4, 1, 2, 6, 7) and ETM+ (bands 1, 2, 3, 4, 5, 7) sensor. MODIS has narrower spectral responses than ETM+ in all the bands. A short-term radiometric stability was evaluated using continuous ETM+ scenes within the contact period and the corresponding half orbit MODIS scenes. The near simultaneous earth observations (SNO) were limited by the smaller swath size of ETM+ (187 km) as compared to MODIS (2330 km). Two sets of continuous granules for MODIS and ETM+ were selected and mosaiced based on pixel geolocation information for non cloudy pixels over the North American continent. The Top-of- Atmosphere (TOA) reflectances were computed for the spectrally matching bands between ETM+ and MODIS over the regions of interest (ROI). The matching pixel pairs were aggregated from a finer to a coarser pixel resolution and the TOA reflectance values covering a wide dynamic range of the sensors were compared and analyzed. Considering the uncertainties of the absolute calibration of the both sensors, radiometric stability was verified for the band pairs. The Railroad Valley Playa, Nada (RVPN) was included in the path of this continuous orbit, which served as a verification point between the shortterm and the long-term trending results from previous studies. This work focuses on monitoring the short-term on-orbit stability of MODIS and the ETM+ RSB. It also provides an assessment of the absolute calibration differences between the two sensors over their wide dynamic ranges.

Choi, Taeyoung↗
NASA NTRS2002

Response of thyroid follicular cells to gamma irradiation compared to proton irradiation: II. The role of connexin 32

The objective of this study was to determine whether connexin 32-type gap junctions contribute to the "contact effect" in follicular thyrocytes and whether the response is influenced by radiation quality. Our previous studies demonstrated that early-passage follicular cultures of Fischer rat thyroid cells express functional connexin 32 gap junctions, with later-passage cultures expressing a truncated nonfunctional form of the protein. This model allowed us to assess the role of connexin 32 in radiation responsiveness without relying solely on chemical manipulation of gap junctions. The survival curves generated after gamma irradiation revealed that early-passage follicular cultures had significantly lower values of alpha (0.04 Gy(-1)) than later-passage cultures (0.11 Gy(-1)) (P < 0.0001, n = 12). As an additional way to determine whether connexin 32 was contributing to the difference in survival, cultures were treated with heptanol, resulting in higher alpha values, with early-passage cultures (0.10 Gy(-1)) nearly equivalent to untreated late-passage cultures (0.11 Gy(-1)) (P > 0.1, n = 9). This strongly suggests that the presence of functional connexin 32-type gap junctions was contributing to radiation resistance in gamma-irradiated thyroid follicles. Survival curves from proton-irradiated cultures had alpha values that were not significantly different whether cells expressed functional connexin 32 (0.10 Gy(-1)), did not express connexin 32 (0.09 Gy(-1)), or were down-regulated (early-passage plus heptanol, 0.09 Gy(-1); late-passage plus heptanol, 0.12 Gy(-1)) (P > 0.1, n = 19). Thus, for proton irradiation, the presence of connexin 32-type gap junctional channels did not influence their radiosensitivity. Collectively, the data support the following conclusions. (1) The lower alpha values from the gamma-ray survival curves of the early-passage cultures suggest greater repair efficiency and/or enhanced resistance to radiation-induced damage, coincident with the expression of connexin 32-type gap junctions. (2) The increased sensitivity of FRTL-5 cells to proton irradiation was independent of their ability to communicate through connexin 32 gap junctions. (3) The fact that the beta components of the survival curves from both gamma rays and proton beams were similar (average 0.022 +/- 0.008 Gy(-2), P > 0.1, n = 39) suggests that at higher doses the loss of viability occurs at a relatively constant rate and is independent of radiation quality and the presence of functional gap junctions.

Non-NASA Center↗
NASA NTRS2000

Physics of Hard Spheres Experiment: Significant and Quantitative Findings Made

Direct examination of atomic interactions is difficult. One powerful approach to visualizing atomic interactions is to study near-index-matched colloidal dispersions of microscopic plastic spheres, which can be probed by visible light. Such spheres interact through hydrodynamic and Brownian forces, but they feel no direct force before an infinite repulsion at contact. Through the microgravity flight of the Physics of Hard Spheres Experiment (PHaSE), researchers have sought a more complete understanding of the entropically driven disorder-order transition in hard-sphere colloidal dispersions. The experiment was conceived by Professors Paul M. Chaikin and William B. Russel of Princeton University. Microgravity was required because, on Earth, index-matched colloidal dispersions often cannot be density matched, resulting in significant settling over the crystallization period. This settling makes them a poor model of the equilibrium atomic system, where the effect of gravity is truly negligible. For this purpose, a customized light-scattering instrument was designed, built, and flown by the NASA Glenn Research Center at Lewis Field on the space shuttle (shuttle missions STS 83 and STS 94). This instrument performed both static and dynamic light scattering, with sample oscillation for determining rheological properties. Scattered light from a 532- nm laser was recorded either by a 10-bit charge-coupled discharge (CCD) camera from a concentric screen covering angles of 0 to 60 or by sensitive avalanche photodiode detectors, which convert the photons into binary data from which two correlators compute autocorrelation functions. The sample cell was driven by a direct-current servomotor to allow sinusoidal oscillation for the measurement of rheological properties. Significant microgravity research findings include the observation of beautiful dendritic crystals, the crystallization of a "glassy phase" sample in microgravity that did not crystallize for over 1 year in 1g (Earth's gravity), and the emergence of face-centered-cubic (FCC) crystals late in the coarsening process (as small crystallites lost particles to the slow ripening of large crystallites). Significant quantitative findings from the microgravity experiments have been developed describing complex interactions among crystallites during the growth process, as concentration fields overlap in the surrounding disordered phase. Time-resolved Bragg scattering under microgravity captures one effect of these interactions quite conclusively for the sample at a volume fraction of 0.528. From the earliest time until the sample is almost fully crystalline, the size and overall crystallinity grow monotonically, but the number of crystallites per unit volume (number density) falls. Apparently nucleation is slower than the loss of crystallites because of the transfer of particles from small to large crystals. Thus, coarsening occurs simultaneously with growth, rather than following the completion of nucleation and growth as is generally assumed. In the same sample, an interesting signature appears in the apparent number density of crystallites and the volume fraction within the crystallites shortly before full crystallinity is reached. A brief upturn in both indicates the creation of more domains of the size of the average crystallite simultaneous with the compression of the crystallites. Only the emergence of dendritic arms offers a reasonable explanation. The arms would be "seen" by the light scattering as separate domains whose smaller radii of curvature would compress the interior phase. In fiscal year 1999, numerous papers, a doctoral dissertation, and the PHaSE final report were produced. Although this flight project has been completed, plans are in place for a follow-on colloid experiment by Chaikin and Russel that employs a light microscope within Glenn's Fluids and Combustion Facility on the International Space Station. PHaSE is providing us with a deeper understanding of the nure of phase transitions. The knowledge derived has added to the understanding of condensed matter. In addition, the burgeoning study of the dynamics of colloidal self-assembly may lead to the development of a range of photonic materials that control the desirable properties of light. Thus, applications of ordered colloidal structures include not only ultrastructure ceramics, but also photonic crystals and photothermal nanosecond light-switching devices. Industries dealing with semiconductors, electro-optics, ceramics, and composites stand to benefit from such advancements.

Doherty, Michael P.↗