Search NASASearch

SEARCH · Search NASA

Results for “high-temperature refractory thermal energy storage”

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.

969 records · Page 32

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

AC and DC Fault Management for Megawatt Electrified Aircraft Electrical Powertrains - Task 2: Power Quality Filtering Using Nanocrystalline Soft Magnetic Inductor

The NASA RTAPS program on AC and DC Fault Management for Megawatt Electrified Power Train is a multi-year joint project with Pratt & Whitney (P&W), Collins Aerospace (CA), and RTX Technology Research Center (RTRC). This research program focuses on the high-voltage distribution issues that present a significant technological obstacle in the adoption of Electrified Aviation Propulsion (EAP) systems. One challenge to the adoption of high-voltage distribution systems with power electronic converters is the need for filter elements to limit the generation and propagation of noise, protect the cable systems from premature aging and prevent against excessive heating within subcomponents due to high-frequency induced currents. While increased distribution voltages aide in reducing the cable mass for a fixed power system, the associated mass with the filtering elements for power electronic converter can grow with increasing distribution voltages – thereby mitigating any benefit associated with increasing the distribution system voltage. To enable high-voltage distribution systems with high system specific power densities, new magnetic materials must be developed. Therefore, the second task of the NASA RTAPS program is associated with the design and application of advanced soft magnetic components for Megawatt class electric propulsion systems, specifically the motor drive system. This report covers the collaborative work between NASA Glenn Research Center (GRC), RTRC, P&W and CA in the development of three types of magnetic components over the span of the three-year program. These critical magnetic components are the DC side EMI filter, which limits the propagation of harmful electromagnetic noise to the rise of the distribution system, and the AC side damping with the dv/dt filter, which limits the fast rise time of the power electronic converter output voltage to limit the degradation on the cable/motor insulation systems. Each of these components are investigated from component level design and are optimized at the system level with a combined modelling and testing effort. In the final experimental evaluation of the NASA developed soft magnetic material with a dv/dt filter, a commercial-off-the-shelf (COTS) magnetic core and the GRC magnetic core are optimized and loaded at 320Arms to evaluate their difference in performance. After a run time of 30 minutes in a MW-class motor driver at RTRC, the NASA GRC cores were found to not only offer a lower temperature rise of nearly 25°𝐶, but also a reduction in measured core loss of 25% (12.75W to 9.5W).

Elecrified Aircraft Propulsion

Thermodynamic Properties of Uranium-Zirconium Carbide Solid Solutions with Oxygen

Uranium-zirconium carbide solid solutions are potential candidates for fuel in space nuclear reactors. To support development of advanced robust fuels for space applications, (U y Zr 1-y )C x O z solid solutions were synthesized via carbothermic reduction with uranium, carbon, and oxygen contents raging from y=0.05 – 0.3, x=0.95 – 0.98, and z = 0.02 – 0.07. These solid solutions were subsequently characterized by X-ray diffraction and elemental combustion analysis. High-temperature oxidative drop solution calorimetry in molten sodium molybdate solvent at 1073 K was employed to determine the formation, oxidation, and mixing enthalpies of (U y Zr 1-y )C x O z solid solutions, and thereby for evaluating their thermodynamic stability. The solid solutions are enthalpically less stable than the UC and ZrC end members. They also become less stable and more reactive towards oxidation with increasing uranium content. Incorporation of oxygen into solid solutions increases their solid solution-stability and improves their oxidation resistance.

Calorimetry

The Operational NASA Anomaly Gas Analyzer and Exploration Follow-Ons

The NASA Anomaly Gas Analyzer (AGA) is the culmination of nearly ten years of advancement of core technology originally developed by Vista Photonics through the Small Business Innovation Research (SBIR) program and expanded using NASA program funding. The AGA is a portable, battery operated, optical gas detection instrument for continuous monitoring of H 2 O, CO 2 , O 2 , CO, NH 3 , HCN, HF, HCl in the spacecraft environment. Real time, in-situ monitored, gas concentrations are communicated at 1 Hz through a built-in display on-orbit and additionally through a serial interface on the ground. The instruments function over the typical range of temperatures and pressures required for the spacecraft environment. A ten-year operational life is targeted. Seven AGA instruments are currently deployed on the International Space Station and two Orion-specific units are manifested for Artemis 2. A variation of the operational AGA is under development for lunar Human Landing Systems. A wider calibrated operating pressure envelope is required in this application. Otherwise, the instruments straightforwardly address lessons learned from the original AGA by improving aspects of thermal control, ruggedness, and ease-of-use. A subset of AGA-developed sensed gases (H 2 O, CO 2 , O 2 ) is the subject of the Primary Constituent Monitor (PCM) development for the Habitation and Logistics Outpost (HALO) destined for lunar orbit. Two PCMs will be installed in HALO as part of the environmental control loop where they will interface with vehicle power and communications. While the basic sensor techniques are unchanged from the AGA, these instruments have been radiation hardened to survive for many years in the harsh target environment. Current development status for these instruments will be presented along with recent developments of a naval submarine variant.

Air Monitoring

Space Transportation System Technology Symposium. Volume 2 - Dynamics and Aeroelasticity

The Space Shuttle, being an hybrid – an airplane and a launch vehicle – represents the greatest challenge that the dynamicist and the aeroelastician have faced. Some specific problem areas related to the Space Shuttle are listed on figure 1. Dynamics and aeroelasticity envelop many disciplines, including aerodynamics, vibration, random processes, structures, fluid flow, mechanics, etc., but, of more importance, they involve the interaction and coupling of many of these various disciplines. Fundamentally, we are concerned with structural integrity and safe flight, i.e., trying to ensure that the vehicle will remain structurally intact as well as function properly in the presence of the many faceted dynamic environment. A new area which may have an impact on our task is the effect of the high temperature environment. In the past, we have been able to successfully decouple the temperature effects from our problem formulation. For the Space Shuttle, this problem must be closely reexamined. As has been pointed out in the opening remarks of the Conference, the Dynamics and Aeroelasticity Technology Group comprises one of several technology groups which are attempting to provide the necessary research to support a successful and safe vehicle. The group has members from most of the NASA Centers as well as from the Air Force groups. We meet periodically to review ongoing work, search for new problem areas; and we are constantly updating and revising our program. The group is organized into three panels as shown on figure 2: a panel on Dynamic Loads and Response, one on Aeroelasticity, and one on Flight Dynamics and Environment. The Conference papers accordingly have been grouped in these same three areas, with each Panel Chairman acting as moderator for his particular session.

Harry L Runyan

3D Woven Mid-Density Carbon Phenolic (3MDCP) Full-Scale Heatshield Development for Mars Sample Return (MSR) Earth Entry System (EES)

The Mars Sample Return (MSR) mission will be returning samples of Martian soil to Earth. The samples will be returned in the Earth Entry System (EES). The EES capsule is protected by a new thermal protection system; the 3-Dimensionally Woven, Mid-Density, Carbon Phenolic (3MDCP) material which is derived from the Heatshield for Extreme Entry Environment Technology (HEEET). The baseline 3MDCP design is a single piece thermal protection system that avoids the manufacturing and certification challenges associated with a tiled configuration. A 3MDCP heatshield begins as a flat woven preform that is formed to the final heatshield shape and then infused with phenolic resin. A NASA Ames lead team, in collaboration with Spirit Textiles (formerly TEAM Inc.), Fiber Materials Inc. (FMI, a Spirit AeroSystems Company), and Kratos SRE, has been developing and demonstrating the processes to fabricate the 3MDCP heatshield and characterize the resulting properties. The process of forming the flat woven preform into the final heatshield shape, a sphere-cone geometry with a 52.5° cone angle, involves local but substantial movement of the yarns in the weave. This results in continuous fibers across the single piece heatshield, albeit with property variations between different regions of the heatshield. This presentation will provide a high-level status of 3MDCP development for the MSR EES heatshield. This will include an overview of the manufacturing processes with an emphasis on the impact of forming on fiber orientation, material properties, and performance. Results comparing preliminary formed and flat materials will be shown. Additionally, the presentation will layout the broader plan for property testing and tie-in to aerothermal performance.

Peter Edward Marshall

Ranger Spacecraft

The thermal-control philosophy of the spacecraft currently under development by the Jet Propulsion Laboratory is design by passive means to maintain all components within the tolerances specified by cognizant engineers. Due to the complexity of the configurations, calculations are) of necessity, fairly generalized and final design is based upon tests in an environmental chamber. The Ranger series spacecraft is designed with a basic structure which is common to all models, with additional hardware to suit the individual mission. This basic structure of Rangers A-1 and A-2 is seen as the hexagonal instrument section, the erectable solar panels, the movable antenna, and the omniantenna. The Ranger A-1 and A-2 configuration is for engineering tests and space-exploration, with the scientific instrumentation isolation requirement dictating the spread-out design. The spacecraft stands 12 feet high, weighs 700 to 800 pounds, and has an internal power of 150 watts. Rangers A-3, A-4, and A-5 are designed to rough land a capsule on the moon. For these, a capsule and retrorocket replace the scientific instruments, occupying the space inside the tower structure. The spacecraft must survive many environments. Chronologically they are: 1) Folded configuration inside an aerodynamic shroud on the pad. 2) Thermal flux from shroud aerodynamically heated during boost phase. 3) Coasting up to 30 minutes attached to Agena stage after booster and shroud are separated. 4) Agena stage burning. 5) Coasting and tumbling after separation from Agena until it passes from earth's shadow. 6) Upon reaching sunlight, panels open and begin sun acquisition. 7) Antenna seeks earth after spacecraft locks onto sun. 8) Space phase- "steady state" with vehicle's vertical axis locked on sun, communicating with earth. The philosophy is to design for the sun-acquired mode, making allowances for the transient conditions.

T. O. Thostesen

Battery charge control: Which approach is best?

The question of the best approach for battery charge control is addressed from two points of view: from the battery point-of-view; and from the spacecraft/mission point-of-view. In terms of performance, more aggressive charging can result in higher discharge voltage and capacity; however, aggressive charging usually implies more overcharge at higher rates and a concomitant higher end of charge temperature (higher stress, shorter life). In terms of life, benign charging can result in lower discharge voltage and capacity (poorer performance). Additionally, benign, or low stress charging, usually implies less overcharge at lower rates with lower end of charge temperature (lower stress, longer life). A discussion of these topics is presented in viewgraph form.

Lurie, Charles

Test Data Analysis of the Thermodynamic Vent System-Augmented Top Spray Injector Liquid Nitrogen Transfer Experiments

Traditionally, a cryogenic tank must be pre-chilled to some “target” temperature before the main vent valve can be closed to attempt a non-vented fill (NVF) of cryogenic liquid propellant. This methodology is particularly attractive for performing in-space transfer of cryogens due to the unknown location of the liquid/vapor interface in microgravity and the high likelihood of venting liquid if the vent valve is opened during transfer. This paper presents in-depth test data analysis of a Thermodynamic Vent System (TVS) augmented injector used for cryogenic tank chilldown and fill experiments of a thin-walled Titanium tank. Eight tests were conducted using liquid nitrogen across a range of inlet conditions and boundary conditions, and three different chilldown/fill methods. For four of the tests, the injector sprays liquid into the tank as normal, but also uses a TVS heat exchanger to cool the metallic injector itself as well as the main incoming liquid stream. Results show that using the TVS augmented injector simplifies transfer operation via enhanced condensation at the injector surface at the cost of sacrificing only a small amount of propellant.

Thermodynamic Vent System

A Panchromatic View of the Restless SN 2009ip Reveals the Explosive Ejection of a Massive Star Envelope

The double explosion of SN 2009ip in 2012 raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN 2009ip during its remarkable rebrightenings. High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the Very Large Array, Swift, Fermi, Hubble Space Telescope, and XMM) constrain SN 2009ip to be a low energy (E (is) approximately 10(exp 50) erg for an ejecta mass approximately 0.5 solar mass) and asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at approximately 5 × 10( exp 14) cm with M approximately 0.1 solar mass, ejected by the precursor outburst approximately 40 days before the major explosion. We interpret the NIR excess of emission as signature of material located further out, the origin of which has to be connected with documented mass-loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic-ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, which later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the massive luminous progenitor star survived. Irrespective of whether the explosion was terminal, SN 2009ip brought to light the existence of new channels for sustained episodic mass loss, the physical origin of which has yet to be identified.

view

Lyotropic Liquid Crystalline Phase Behavior of Boron Nitride Nanotube Aqueous Dispersions

Boron nitride nanotubes (BNNTs) are gaining significant interest due to their outstanding mechanical and thermal properties, as well as their potential to serve as a model nanorod system. Processing BNNT liquid crystalline (LC) dispersions enables precise control over BNNT orientation in macroscopic assemblies, while their low absorption in the visible spectrum facilitates studying LCs at exceptionally high concentrations. Here, we investigate the behavior of BNNTs in aqueous solutions stabilized by the surfactant sodium deoxycholate (SDC), examining the effect of BNNT purity and BNNT−SDC concentrations on lyotropic LC formation. We disperse up to 15wt% BNNT in SDC solutions and use polarized light microscopy to detail the transition from an isotropic state to a biphasic regime, where isotropic and nematic domains coexist due to phase separation, to a single fully nematic phase. Cryogenic electron microscopy provides direct evidence of BNNT alignment within nematic domains. Our results show that enhanced depletion-induced attractions, driven by increased surfactant concentration, lower the threshold concentration of BNNT required to form nematic domains. In contrast, low surfactant concentrations relative to BNNT result in insufficiently coated nanotube surfaces, leading to poor dispersions and BNNT aggregation. Additionally, we fabricate well-aligned BNNT films from aqueous LC dispersions. Our findings advance the understanding of BNNT LCs, offering insight into controlling their orientation and highlighting their potential for high-performance materials.

Joe Khoury

An ESR study of the UV degradation of FEP

Spacecraft in low earth orbit are subjected to significant levels of high energy radiation, including ultraviolet (UV) and visible ultraviolet (VUV) wavelengths. The effects of UV radiation are enhanced over those at the surface of the earth, where the only incident wavelengths are greater than 290 nm. In low earth orbit the incident UV wavelengths extend below 290 nm into the VUV region, where the Lyman alpha-emissions of atomic hydrogen occur at 121 nm. In addition to electromagnetic radiation, in low earth orbit polymer materials may also be subjected to atomic oxygen particle radiation, which will result in direct oxidation of the polymer.

G. A. George

X-Ray Polarization From the Atoll 4U 1735−44 Suggests a Low Inclination

X-ray polarimetry is a new tool capable of probing the geometry of accretion onto weakly magnetized neutron stars. Here we present the first X-ray spectro-polarimetric results from coordinated observations of the atoll source 4U 1735−44, conducted with the Imaging X-ray Polarimetry Explorer (IXPE), NICER, and NuSTAR . Over the 2–8 keV energy range, we obtained a marginal detection of polarization with the polarization degree of 1.4% ± 0.7% and polarization angle of −29◦ ± 14◦, corresponding to a 3σ upper limit on the polarization degree of 3.5%. The best-fit model to describe the spectrum comprises a thermal component associated with the accretion disk, a Comptonized black body component, and a relativistic reflection component. From the reflection model, we infer a disk inclination of ∼40◦. The spectroscopic and polarimetric properties of 4U 1735−44 are consistent with those observed in other atoll sources studied by IXPE, with its low polarization likely due to its low inclination

X-rays: binaries

The Slow Demise of the Long-Lived Sn 2005ip

The Type IIn supernova (SN IIn) 2005ip is one of the most well-studied and long-lasting examples of an SN interacting with its circumstellar environment. The optical light curve plateaued at a nearly constant level for more than five years, suggesting ongoing shock interaction with an extended and clumpy circumstellar medium (CSM). Here, we present continued observations of the SN from ∼1000 to 5000 d post-explosion at all wavelengths, including X-ray, ultraviolet, near-infrared (NIR), and mid-infrared. The UV spectra probe the pre-explosion mass loss and show evidence for CNO processing. From the bolometric light curve, we find that the total radiated energy is in excess of 10 50 erg, the progenitor star’s pre-explosion mass-loss rate was >~ 1 × 10 −2 Mꙩ yr −1 , and the total mass lost shortly before explosion was >~1Mꙩ, though the mass lost could have been considerably larger depending on the efficiency for the conversion of kinetic energy to radiation. The ultraviolet through NIR spectrum is characterized by two high-density components, one with narrow high-ionization lines, and one with broader low-ionization H I, He I, [OI], Mg II, and Fe II lines. The rich Fe II spectrum is strongly affected by Lyα fluorescence, consistent with spectral modelling. Both the Balmer and He I lines indicate a decreasing CSM density during the late interaction period. We find similarities to SN 1988Z, which shows a comparable change in spectrum at around the same time during its very slow decline. These results suggest that, at long last, the shock interaction in SN 2005ip may finally be on the decline.

stars

LAI Assimilation Schedules to Constrain Uncertain Cultivars and Soils in CERES-Maize

In anticipation of large-domain crop model applications where precise local configuration and calibration is not possible, we describe benefits and potential drawbacks of employing a crop pest module to achieve leaf area index (LAI) assimilation into a high performing CERES-Maize crop model configuration at a field experiment site in Perry, Iowa. Simulation experiments explore the use of MODIS satellite-derived LAI to constrain and adjust LAI to counter imprecise cultivar and soil configurations often occurring in the absence of high-quality local information. Simulations using single-day, window, and continuous LAI replacement across 5 cultivars and 2 soil calibration approaches for nine corn rotation years from 2004 to 2020 led to different yield outcomes and reverberations throughout the field environment. Evaluating variance and mean bias, results indicate minimal interventions in early vegetative and grain-filling stages were more beneficial than use of continuous LAI adjustments, as they minimized disruptions to the internal resource balances governing plant stresses and grain production. LAI adjustment was particularly helpful in constraining growth related to uncertain thermal unit requirements and leaf tip appearance rates (P1 and PHINT cultivar parameters, respectively). Findings underscore the need to assimilate additional state variables to ensure internal biophysical coherence. This approach shows promise for applications spanning wider domains with prediction time pressure where detailed configuration, more complex assimilation methods, or recalibration of crop model parameters may not be practical.

Phenology

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES

Applications Technology Satellites Battery and Power System Design

A summary of the ATS battery design which is onboard the Applications Technology Satellite (ATS) is provided. The 15 ampere hour nickel cadmium cells were manufactured by Gulton, 19 series connected cells per battery, and there are two batteries in each spacecraft. The operating design life was two years in a synchronous orbit, and a maximum depth of discharge of 50 percent. The design temperature for the batteries in the spacecraft was 0 to 25 C, and the charge control consisted of 1 volt versus temperature on a constant percentage voltage. Also, C/10 current limit, and a commandable trickle charge rate, using C/20 or C/60. The undervoltage was sent across a 9 cell and a 10 cell group, and it was set at one volt average per group on either group.

F E Ford