SEARCH · Search NASA

Results for “zone refining”

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Habitable Zones Around Main-Sequence Stars: Dependence on Planetary Mass

The ongoing discoveries of extra-solar planets are unveiling a wide range of terrestrial mass (size) planets around their host stars. In this Letter, we present estimates of habitable zones (HZs) around stars with stellar effective temperatures in the range 2600 K-7200 K, for planetary masses between 0.1M and 5M. Assuming H2O-(inner HZ) and CO2-(outer HZ) dominated atmospheres, and scaling the background N2 atmospheric pressure with the radius of the planet, our results indicate that larger planets have wider HZs than do smaller ones. Specifically, with the assumption that smaller planets will have less dense atmospheres, the inner edge of the HZ (runaway greenhouse limit) moves outward (approx.10% lower than Earth flux) for low mass planets due to larger greenhouse effect arising from the increased H2O column depth. For larger planets, the H2O column depth is smaller, and higher temperatures are needed before water vapor completely dominates the outgoing long-wave radiation. Hence the inner edge moves inward (approx.7% higher than Earth's flux). The outer HZ changes little due to the competing effects of the greenhouse effect and an increase in albedo. New, three-dimensional climate model results from other groups are also summarized, and we argue that further, independent studies are needed to verify their predictions. Combined with our previous work, the results presented here provide refined estimates of HZs around main-sequence stars and provide a step toward a more comprehensive analysis of HZs.

main sequence starts↗

NASA NTRS2003

Mineralogy and Petrology of Amoeboid Olivine Inclusions in CO3 Chondrites: Relationship to Parent-Body Aqueous Alteration

Petrographic and mineralogic studies of amoeboid olivine inclusions (AOIs) in CO3 carbonaceous chondrites reveal that they are sensitive indicators of parent-body aqueous and thermal alteration. As the petrologic subtype increases from 3.0 to 3.8, forsteritic olivine (Fa(sub 0-1)) is systematically converted into ferroan olivine (Fa(sub 60-75)). We infer that the Fe, Si and O entered the assemblage along grain boundaries, forming ferroan olivine that filled fractures and voids. As temperatures increased, Fe(+2) from the new olivine exchanged with Mg(+2) from the original AOI to form diffusive haloes around low-FeO cores. Cations of Mn(+2), Ca(+2) and Cr(+3) were also mobilized. The systematic changes in AOI textures and olivine compositional distributions can be used to refine the classification of CO3 chondrites into subtypes. In subtype 3.0, olivine occurs as small forsterite grains (Fa(sub 0-1)), free of ferroan olivine. In petrologic subtype 3.2, narrow veins of FeO-rich olivine have formed at forsterite grain boundaries. With increasing alteration, these veins thicken to form zones of ferroan olivine at the outside AOI margin and within the AOI interior. By subtype 3.7, there is a fairly broad olivine compositional distribution in the range Fa(sub 63-70), and by subtype 3.8, no forsterite remains and the high-Fa peak has narrowed, Fa(sub 64-67). Even at this stage, there is incomplete equilibration in the chondrite as a whole (e.g., data for coarse olivine grains in Isna (CO3.8) chondrules and lithic clasts show a peak at Fa(sub39)). We infer that the mineral changes in A01 identified in the low petrologic types required aqueous or hydrothermal fluids whereas those in subtypes greater than or equal to 3.3 largely reflect diffusive exchange within and between mineral grains without the aid of fluids.

Chizmadia, Lysa J.↗

NASA NTRS1966

Orbital tube flare for tubing connectors

Marshall engineering personnel have developed and are now refining a new and unique method of mechanically flaring tubing. The system eliminates the use of the conventional split backing die which is used with most tube-flaring equipment. The Marshall-developed tool is adapted for use on standard Leonard 3-CP-HD flaring machines. The stationary split backing die is replaced in function with a die ring, whose i.d.-surface rolls or orbits around the o.d. of the tube flare zone. An i.d.-flaring cone moves axially and forces the metal against the rolling die to form the flare. Major advantages of this system are the elimination of flash, or split die marks, on the back angle of the flare, and improved flare-surface finishes, permitting more uniform sealing characteristics and improvements of the total geometric configuration of the flare as required by the MSFC specification MC-146. The Orbital Adapter system of tooling is currently being used to manufacture flared tube hardware for the Saturn V S-IO space launch vehicle at MSFC, for the Agena space launch vehicle (used for the Gemini project) at Lockheed, Sunnyvale, for LEM checkout equipment for Grumman, and for the Saturn S-IV-B space launch vehicle by Douglas Aircraft Company. Several other aerospace firms are evaluating the system for possible incorporation on their tube flaring requirements. Marshall engineers are now concentrating their efforts on refining and reducing the total cost of the overall system.

Tubing↗

NASA NTRS1998

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗

NASA NTRS1993

Refinements on the pinhole camera measurements of the LDEF attitude

The results from the University of Alabama in Huntsville (UAH) pinhole camera were reanalyzed to include the effects of corotation of the atmosphere with the Earth as well as satellite oscillation. Previous results from the instrument showed that the satellite had stable attitude offsets in yaw of 8.0 deg and 1.0 deg in pitch; these offsets are unchanged by the present analysis. The primary impact zone of oxygen, i.e., the directly exposed spot on a silver detector, had a ratio of major to minor axes equal to 1.05, which was interpreted as being caused by a small oscillation of plus or minus 0.35 degrees (with precision plus or minus 0.15 degrees). The present analysis shows that the observed effect can largely be accounted for by atmospheric corotation, but that an additional oscillation in yaw of the order of a degree cannot be excluded. The sensitivity of the pinhole camera to satellite oscillations is shown to decrease nonlinearly with decreasing magnitude of the oscillation and to vary inversely with the gas temperature.

Peters, Palmer N.↗

NASA NTRS2022

Measurement of Surface Deformation Related to the December 2018 Mt. Etna Eruption Using Time-Series Interferometry and Magma Modeling for Hazard Zone Mapping

Mount Etna has erupted several times since it was first formed. Recently, Mount Etna began erupting again over 24–27 December 2018. Because it erupts frequently, Mount Etna should be observed on a frequent basis. From June 2018 to October 2019, 34 and 56 interferometric synthetic aperture radar (InSAR) images were acquired from the ascending and descending tracks of the Sentinel-1 satellite, respectively. We employed the Stanford Method for Persistent Scatterers (StaMPS) and a refined small baseline subset (SBAS) InSAR method to produce a surface deformation time-series map. In the time-series analysis, the phase signal remained unaltered with time. The Okada model was then applied to the result to generate a modeled interferogram, and the Q-LavHA program was run to generate a lava flow prediction model. A direct comparison of the results showed that Persistent Scatterers Interferometry (PSI)-StaMPS and the refined SBAS technique were comparable in terms of the displacement pattern, with slightly different velocity values obtained for individual points. In particular, a velocity range of −25 to 21 cm/yr was obtained from PSI-StaMPS, whereas a range of −30 to 25 cm/yr was obtained from the refined SBAS method. Upon computation of the vertical and east-west displacement components based on ascending and descending track data using both methods, deformation velocities of 51.5 and 52.5 cm/yr in the westerly direction on the western flank of Mount Etna were obtained from PSI-StaMPS and the refined SBAS method, respectively, whereas on the eastern flank, deformation toward the east was estimated to occur at a velocity of 50.1 or 54.2 cm/yr, respectively. PSIStaMPS estimated a vertical deformation velocity of −5.3 to 18.3 cm/yr, whereas the refined SBAS method produced a velocity range of approximately −7 to 19 cm/yr. The interferogram obtained via Okada modeling showed two fault sources in the 2018 Mount Etna eruption and a total volume change of approximately 12.39 × 10 6 m 3 . From the modeling results, a lava flow prediction model was generated using the Q-LavHA program. The approaches described in this study can be used by government officials, authorities, and other decision-makers to monitor and assess the risk of volcanic activity in the region.

Suci Ramayanti↗

NASA NTRS1980

Numerical calculation of steady inviscid full potential compressible flow about wind turbine blades

The air flow through a propeller-type wind turbine rotor is characterized by three-dimensional rotating cascade effects about the inner portions of the rotor blades and compressibility effects about the tip regions of the blades. In the case of large rotor diameter and/or increased rotor angular speed, the existence of small supersonic zones terminated by weak shocks is possible. An exact nonlinear mathematical model (called a steady Full Potential Equation - FPE) that accounts for the above phenomena has been rederived. An artificially time dependent version of FPE was iteratively solved by a finite volume technique involving an artificial viscosity and a three-level consecutive mesh refinement. The exact boundary conditions were applied by generating a boundary conforming periodic computation mesh.

Dulikravich, D. S.↗

NASA NTRS1994

An outer gap model of high-energy emission from rotation-powered pulsars

We describe a refined calculation of high-energy emission from rotation-powered pulsars based on the outer gap model of Cheng, Ho, & Ruderman (1986 a, b). In this calculation, vacuum gaps form in regions near the speed-of-light cylinder of the pulsar magnetosphere along the boundary between the closed and open field line zones. We have improved upon previous efforts to model the spectra from these pulsars (e.g., Cheng et al. 1986b; Ho 1989) by following the variation in particle production and radiation properties with position in the outer gap. Curvature, synchotron, and inverse-Compton scattering fluxes vary significantly over the gap, and their interactions via photon-photon pair production build up the radiating charge populations at varying rates. We have also incorporated an approximate treatment of the transport particle and photon fluxes between gap emission zones. These effects, along with improved computations of the particle and photon distribution, provide very important modifications of the model gamma-ray flux. In paticular, we attempt to make specific predictions of pulse profile shapes and spectral variations as a function of pulse phase and suggest further extensions to the model which may provide accurate computations of the observed high-energy emissions.

Chiang, James↗

NASA NTRS2016

Multi-Frequency Investigation into Scattering from Vegetation over the Growth Cycle

In this investigation, we aim to collect and use time-series multi-frequency microwave data over winter wheat during entire growth cycle to characterize vegetation dynamics and to quantify its effects on soil moisture retrievals. We plan to incorporate C-band radar and VHF receiver within the existing L-band radarradiometer system called ComRAD (SMAPs ground based simulator). With C-bands ability to sense vegetation details and VHFs root-zone soil moisture within ComRADs footprint, we will be able to test our discrete scatterer vegetation models and parameters at various surface conditions. The purpose of this study is to determine optical depth and effective scattering albedo of vegetation of a given type (i.e. winter wheat) at various stages of growth that are need to refine soil moisture retrieval algorithms being developed for the SMAP mission.

Lang, R. H.↗

NASA NTRS2005

Cyclic Spin Testing of Superalloy Disks With a Dual Grain Microstructure

An aggressive cyclic spin test program was run to verify the reliability of superalloy disks with a dual grain structure, fine grain bore and coarse grain rim, utilizing a disk design with web holes bisecting the grain size transition zone. Results of these tests were compared with conventional disks with uniform grain structures. Analysis of the test results indicated the cyclic performance of disks with a dual grain structure could be estimated to a level of accuracy which does not appear to prohibit the use of this technology in advanced gas turbine engines, although further refinement of lifing methodology is clearly warranted.

Gayda, John↗

NASA NTRS1976

Development and application of operational techniques for the inventory and monitoring of resources and uses for the Texas coastal zone

The author has identified the following significant results. The most significant ADP result was the modification of the DAM package to produce classified printouts, scaled and registered to U.S.G.S., 71/2 minute topographic maps from LARSYS-type classification files. With this modification, all the powerful scaling and registration capabilities of DAM become available for multiclass classification files. The most significant results with respect to image interpretation were the application of mapping techniques to a new, more complex area, and the refinement of an image interpretation procedure which should yield the best results.

Jones, R.↗

NASA NTRS1974

The experimental clean combustor program - Description and status

The experimental clean combustor program is a contract effort with primary objectives including the generation and demonstration of technology for development of advanced commercial CTOL aircraft engines with lower exhaust emissions than current aircraft, and the demonstration of this technology in full scale engines in 1976. The program is being conducted in three phases. These consist of screening of low pollutant combustors, refinement of the best combustors, and engine demonstration of the best combustors. Pollution goals are emission index values of 20 and 4 for carbon monoxide and unburned hydrocarbons respectively at engine idle conditions, and an oxides of nitrogen emission index level of 10 at engine takeoff conditions. Pollution reduction approaches which are being investigated include multiple burning zone combustors, investigations of improved fuel distribution, and preparation and the staging of combustor airflow.

Niedzwiecki, R. W.↗

NASA NTRS2011

Time-Distance Analysis of Deep Solar Convection

Recently it was shown by Hanasoge, Duvall, and DeRosa (2010) that the upper limit to convective flows for spherical harmonic degrees l<? is considerably smaller than the flows predicted by the ASH simulations (Miesch et a7. ref) at the depth r/R=0.95 ' The deep-focusing Lime-distance technique used to develop the upper limit was applied to linear acoustic simulations of a solar interior perturbed by convective flows in order to calibrate the technique. This technique has been applied to other depths in the convection zone and the results will be presented. The deep-focusing technique has considerable sensitivity to the flow ' signals at the desired subsurface location ' However, as shown by Birch {ref}, there is remaining much sensitivity to near-surface signals. Modifications to the technique using multiple bounce signals have been examined in a search for a more refined sensitivity, or kernel function. Initial results are encouraging and results will be presented'

Duvall, T. L., Jr.↗

NASA NTRS1994

Development of microchannel plate x-ray optics

The goal of this research program was to develop a novel technique for focusing x-rays based on the optical system of a lobster's eye. A lobster eye employs many closely packed reflecting surfaces arranged within a spherical or cylindrical shell. These optics have two unique properties: they have unlimited fields of view and can be manufactured via replication of identical structures. Because the angular resolution is given by the ratio of the size of the individual optical elements to the focal length, optical elements with sizes on the order of one hundred microns are required to achieve good angular resolution with a compact telescope. We employed anisotropic etching of single crystal silicon wafers for the fabrication of micron-scale optical elements. This technique, commonly referred to as silicon micromachining, is based on silicon fabrication techniques developed by the microelectronics industry. An anisotropic etchant is a chemical which etches certain silicon crystal planes much more rapidly than others. Using wafers in which the slowly etched crystal planes are aligned perpendicularly to the wafer surface, it is possible to etch a pattern completely through a wafer with very little distortion. Our optics consist of rectangular pores etched completely through group of zone axes (110) oriented silicon wafers. The larger surfaces of the pores (the mirror elements) were aligned with the group of zone axes (111) planes of the crystal perpendicular to the wafer surface. We have succeeded in producing silicon lenses with a geometry suitable for 1-d focusing x-ray optics. These lenses have an aspect ratio (40:1) suitable for x-ray reflection and have very good optical surface alignment. We have developed a number of process refinements which improved the quality of the lens geometry and the repeatability of the etch process. A significant progress was made in obtaining good optical surface quality. The RMS roughness was decreased from 110 A for our initial lenses to 30 A in the final lenses. A further factor of three improvement in surface quality is required for the production of efficient x-ray optics. In addition to the silicon fabrication, an x-ray beam line was constructed at Columbia for testing the optics.

Kaaret, Philip↗

NASA NTRS2017

Europa Mission Configuration Update to Accommodate Maturing Instrument Designs

A mission to Jupiter’s moon, Europa has been of interest to NASA and JPL since the Galileo mission's magnetometer data predicted the presence of a subsurface ocean. The planned Europa mission would be equipped with a suite of instruments to perform both remote and in-situ sensing, the scope of which ranges from gravity science to characterizing the surface composition, with one of the objectives being to confirm the existence of the subsurface ocean. As the selected instruments mature, the challenge has been to select and refine a spacecraft configuration that is flexible enough to accommodate these changes without degrading the scientific capability of the spacecraft. An increase in instrument volume and power have prompted the growth of spacecraft engineering subsystems, which include the solar arrays and the avionics module. These changes have the potential to cause obstructions to the instrument fields of view, stray light keep-out-zones, and thermal radiative fields of view. Many of these criteria are addressed by the implementation of a dedicated instrument platform. Payloads that are unique to the Europa mission require new approaches to instrument accommodation, such as coupling the ice penetrating radar's radiating elements to the solar array. This paper will discuss these accommodation strategies for the Europa spacecraft to generate a new baseline design.

Gentile, Matthew↗

NASA NTRSDate not supplied

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗

NASA NTRS1987

Advances in the computation of transonic separated flows over finite wings

The transonic flow around a low-aspect-ratio infinite wing in a wind tunnel is investigated by means of numerical simulations, with a focus on shock-induced separated flows. A coarse global grid with far-field boundaries matching those of the test section is subdivided into zones: the flow in the clustered zones near the wing is analyzed by solving the Reynolds-averaged Navier-Stokes equations, while that farther from the wing is modeled with the Euler equations. The results are presented graphically, and it is shown that a mushroomlike separated flow with two counterrotating vortices can be simulated when the correct shock strength is imposed (by careful selection of the artificial dissipation, the boundary conditions, the grid refinement, the algebraic turbulence model, and the geometry representation).

Kaynak, Unver↗

NASA NTRS2017

Summary of Data from the Sixth AIAA CFD Drag Prediction Workshop: CRM Cases 2 to 5

Results from the Sixth AIAA CFD Drag Prediction Workshop Common Research Model Cases 2 to 5 are presented. As with past workshops, numerical calculations are performed using industry-relevant geometry, methodology, and test cases. Cases 2 to 5 focused on force/moment and pressure predictions for the NASA Common Research Model wing-body and wing-body-nacelle-pylon configurations, including Case 2 - a grid refinement study and nacelle-pylon drag increment prediction study; Case 3 - an angle-of-attack buffet study; Case 4 - an optional wing-body grid adaption study; and Case 5 - an optional wing-body coupled aero-structural simulation. The Common Research Model geometry differed from previous workshops in that it was deformed to the appropriate static aeroelastic twist and deflection at each specified angle-of-attack. The grid refinement study used a common set of overset and unstructured grids, as well as user created Multiblock structured, unstructured, and Cartesian based grids. For the supplied common grids, six levels of refinement were created resulting in grids ranging from 7x10(exp 6) to 208x10(exp 6) cells. This study (Case 2) showed further reduced scatter from previous workshops, and very good prediction of the nacelle-pylon drag increment. Case 3 studied buffet onset at M=0.85 using the Medium grid (20 to 40x10(exp 6) nodes) from the above described sequence. The prescribed alpha sweep used finely spaced intervals through the zone where wing separation was expected to begin. Although the use of the prescribed aeroelastic twist and deflection at each angle-of-attack greatly improved the wing pressure distribution agreement with test data, many solutions still exhibited premature flow separation. The remaining solutions exhibited a significant spread of lift and pitching moment at each angle-of-attack, much of which can be attributed to excessive aft pressure loading and shock location variation. Four Case 4 grid adaption solutions were submitted. Starting with grids less than 2x10(exp 6) grid points, two solutions showed a rapid convergence to an acceptable solution. Four Case 5 coupled aerostructural solutions were submitted. Both showed good agreement with experimental data. Results from this workshop highlight the continuing need for CFD improvement, particularly for conditions with significant flow separation. These comparisons also suggest the need for improved experimental diagnostics to guide future CFD development.

Tinoco, Edward N.↗