Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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

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

At least 1,279 records · Page 71

flat10MIP: an emissions-driven experiment to diagnose the climate response to positive, zero and negative CO2 emissions

Abstract. The proportionality between global mean temperature and cumulative emissions of CO2 predicted in Earth system models (ESMs) is the foundation of carbon budgeting frameworks. Deviations from this behavior could impact estimates of required net-zero timings and negative emissions requirements to meet the Paris Agreement climate targets. However, existing ESM diagnostic experiments do not allow for direct estimation of these deviations as a function of defined emissions pathways. Here, we perform a set of climate model diagnostic experiments for the assessment of transient climate response to cumulative CO2 emissions (TCRE), the Zero Emissions Commitment (ZEC), and climate reversibility metrics in an emissions-driven framework. The emissions-driven experiments provide consistent independent variables simplifying simulation, analysis and interpretation, with emissions rates more comparable to recent levels than existing protocols using model-specific compatible emissions from the CMIP DECK 1pctCO2 experiment, where emissions rates tend to increase during the experiment, such that at the time of CO2 doubling in year 70, emissions are much greater than present-day values. A base experiment, “esm-flat10”, has constant emissions of CO2 of 10 GtC per year (near-present-day values), and initial results show that the TCRE estimated in this experiment is about 0.1 K less than that obtained using 1pctCO2. A subset of ESMs exhibit land carbon sinks that saturate during this experiment. A branch experiment, esm-flat10-zec, illustrates that both positive and negative ZEC effects are less pronounced under esm-flat10 than under 1pctCO2 – the magnitude of ZEC50 in ESMs is, on average, reduced by 30 % compared with 1pctCO2 branch experiments. A final experiment, esm-flat10-cdr, assesses climate reversibility under negative emissions, where we find that peak warming may occur before or after net zero and that the asymmetry in temperature at a given level of cumulative emissions between the positive and negative emissions phases is well described by ZEC in most models. Further, we find that existing probabilistic simple climate model (SCM) ensembles tend to overestimate temperature reversibility compared with ESMs, highlighting the need for additional constraints. We propose a set of climate diagnostic indicators to quantify various aspects of climate reversibility. These experiments were suggested as potential candidates in CMIP7 and have since been adopted as “fast track” simulations.

Sanderson, Benjamin M↗

Evaluation of atmospheric density models and preliminary functional specifications for the Langley Atmospheric Information Retrieval System (LAIRS)

This document presents the results of an extensive survey and comparative evaluation of current atmosphere and wind models for inclusion in the Langley Atmospheric Information Retrieval System (LAIRS). It includes recommended models for use in LAIRS, estimated accuracies for the recommended models, and functional specifications for the development of LAIRS.

Lee, T.↗

WIPP PA Sensitivity Studies in Support of CRA-2014 Completeness Determination and Compliance Decision - 18209

The Waste Isolation Pilot Plant (WIPP) has been developed by the U.S. Department of Energy (DOE) for the geologic (deep underground) disposal of transuranic (TRU) waste. Containment of TRU waste at the WIPP is regulated by the U.S. Environmental Protection Agency (EPA). The DOE demonstrates compliance with the containment requirements by means of performance assessment (PA) calculations. WIPP PA calculations estimate the probability and consequence of potential radionuclide releases from the repository to the accessible environment for a regulatory period of 10,000 years after facility closure. Since TRU waste was first received at the WIPP in 1999, compliance recertification applications (CRAs) have been required by the EPA on five-year cycle. The latest CRA (CRA-2014) was submitted by the DOE to the EPA in 2014. As part of the application, a formal set of PA calculations (CRA-2014 PA) was performed and results were documented in detail. The CRA-2014 PA incorporated several parameter updates based on new data, as well as the implementation of a different model for panel closures, barriers planned for the opening of each waste panel. The results presented as part of the CRA-2014 showed via the required complementary cumulative distribution function (CCDF) release curves that the estimated releases were within compliance limits. As part of their process for reviewing the completeness of the CRA-2014 [1], the EPA formally requested that the DOE perform sensitivity studies based on EPA-requested changes to model parameters. The DOE performed and documented the results of three sensitivity studies described below: CRA14-SEN2 [2], CRA14-SEN3 [3], and CRA14-SEN4 [4]. In each case, the new analysis results were compared to CRA- 2014 PA (CRA14) results to determine the impact of the EPA-requested changes. CRA14-SEN2: the application of EPA-requested modified parameters in the operations and experimental (non-waste) areas of the repository to simulate an accelerated (instantaneous) creep closure, the inclusion of capillary pressure effects on relative permeability, and an increase in initial/residual brine saturation and residual gas saturation. The modifications to the repository model resulted in increased pressures and decreased brine saturations in waste areas and increased pressures and brine saturations in the operations and experimental areas. However, the overall effects on total high-probability mean releases from the repository were entirely insignificant. It was concluded that the modeling assumptions associated with the operations and experimental areas of the repository have an insignificant effect on the prediction of total releases from the repository and/or adequacy of the current (CRA14) model to demonstrate compliance with the regulatory limits. CRA14-SEN3: the application of EPA-requested modified parameters in the panel closure areas and operations and experimental (non-waste) areas of the repository to simulate an accelerated (instantaneous) creep closure, the inclusion of capillary pressure effects on relative permeability, and an increase in initial/residual brine saturation and residual gas saturation. The modifications to the repository model predominantly result in increased pressures and decreased brine saturations in waste areas, increased pressures and brine saturations in the operations and experimental areas, and decreased brine and gas flows across all panel closures. The pressure increases in repository waste regions yielded slightly decreased brine saturations (on average) in those areas. However, the overall effects on total high-probability mean releases from the repository are entirely insignificant. It was concluded that the modeling assumptions made for this sensitivity study which were associated with the panel closures and north end areas of the repository have only a small effect on the prediction of total releases from the repository and/or adequacy of the current (CRA14) model to demonstrate compliance with the regulatory limits. CRA14-SEN4: the application of several changes to the current WIPP PA model (CRA14), including: parameter changes, use of an updated code, and a revised computational grid for the BRAGFLO code. The modifications to the repository model resulted in increased releases in all primary release mechanisms. The impacts of each EPA-requested change to CRA14 were analyzed with respect to each release mechanism in the CRA14-SEN4 study. Overall, total high-probability predicted mean releases from the repository were increased by about 15%. Total low-probability predicted mean releases were increased by about 107%. It was concluded that the EPA-requested changes to the CRA14 result in increases to the predicted total releases from the repository. However, releases calculated in the CRA14-SEN4 analysis remained below regulatory limits, demonstrating continued compliance of the WIPP. In January 2017, the EPA determined that CRA-2014 was complete and then began a technical evaluation of CRA-2014 to determine whether the WIPP continued to comply with the radiation protection standards for disposal. During their evaluation, the EPA utilized the sensitivity study results to inform their ultimate determination in July 2017 that the WIPP continues to comply with regulatory standards. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Chandra cosmos legacy survey: clustering of X-ray-selected AGNs at 2.9 ≤ z ≤ 5.5 using photometric redshift probability distribution functions

We present the measurement of the projected and redshift-space two-point correlation function (2pcf) of the new catalog of Chandra COSMOS-Legacy active galactic nucleus (AGN) at 2.9 ≤ z ≤ 5.5 (〈L{sub bol}〉 ∼ 10{sup 46} erg s{sup −1}) using the generalized clustering estimator based on phot-z probability distribution functions in addition to any available spec-z. We model the projected 2pcf, estimated using π {sub max} = 200 h{sup −1} Mpc with the two-halo term and we derive a bias at z ∼ 3.4 equal to b = 6.6{sub −0.55}{sup +0.60}, which corresponds to a typical mass of the hosting halos of log M {sub h} = 12.83{sub −0.11}{sup +0.12} h{sup −1} M {sub ⊙}. A similar bias is derived using the redshift-space 2pcf, modeled including the typical phot-z error σ {sub z} = 0.052 of our sample at z ≥ 2.9. Once we integrate the projected 2pcf up to π {sub max} = 200 h{sup −1} Mpc, the bias of XMM and Chandra COSMOS at z = 2.8 used in Allevato et al. is consistent with our results at higher redshifts. The results suggest only a slight increase of the bias factor of COSMOS AGNs at z ≳ 3 with the typical hosting halo mass of moderate-luminosity AGNs almost constant with redshift and equal to log M {sub h} = 12.92{sub −0.18}{sup +0.13} at z = 2.8 and log M {sub h} = 12.83{sub −0.11}{sup +0.12} at z ∼ 3.4, respectively. The observed redshift evolution of the bias of COSMOS AGNs implies that moderate-luminosity AGNs still inhabit group-sized halos at z ≳ 3, but slightly less massive than observed in different independent studies using X-ray AGNs at z ≤ 2.

79 ASTRONOMY AND ASTROPHYSICS↗

The Relative Importance of Assumed Infrasound Source Terms and Effects of Atmospheric Models on the Linear Inversion of Infrasound Time Series at the Source Physics Experiment

We invert far-field infrasound data for the equivalent seismoacoustic time-domain moment tensor to assess the effects of variable atmospheric models and source phenomena. The infrasound data were produced by a series of underground chemical explosions that were conducted during the Source Physics Experiment (SPE), which was originally designed to study seismoacoustic signal phenomena. The first goal is to investigate the sensitivity of the inversion to the variability of the estimated atmospheric model. The second goal is to determine the relative contribution of two presumed source mechanisms to the observed infrasonic wavefield. Rather than using actual atmospheric observations to estimate the necessary atmospheric Green’s functions, we build a series of atmospheric models that rely on publicly available, regional-scale atmospheric observations. The atmospheric observations are summarized and interpolated onto a 3D grid to produce a model of sound speed at the time of the experiment. For each of four SPE acoustic datasets that we invert, we produced a suite of three atmospheric models for each chemical explosion event, based on 10 yrs of meteorological data: an average model, which averages the atmospheric conditions for 10 yrs prior to each SPE event, as well as two extrema models. To parameterize the inversion, we assume that the source of infrasonic energy results from the linear combination of explosion-induced surface spall and linear seismic-to-elastic mode conversion at the Earth’s free surface. We find that the inversion yields relatively repeatable results for the estimated spall source. Conversely, the estimated isotropic explosion source is highly variable. This suggests that 1) the majority of the observed acoustic energy is produced by the spall and/or 2) our modeling of the elastic energy, and the subsequent conversion to acoustic energy, is too simplistic.

58 GEOSCIENCES↗

Evaluation of probabilistic photometric redshift estimation approaches for The Rubin Observatory Legacy Survey of Space and Time (LSST)

Many scientific investigations of photometric galaxy surveys require redshift estimates, whose uncertainty properties are best encapsulated by photometric redshift (photo-z) posterior probability density functions (PDFs). A plethora of photo-z PDF estimation methodologies abound, producing discrepant results with no consensus on a preferred approach. We present the results of a comprehensive experiment comparing twelve photo-z algorithms applied to mock data produced forLarge Synoptic Survey Telescope The Rubin Observatory Legacy Survey of Space and Time (lsst) Dark Energy Science Collaboration (desc). By supplying perfect prior information, in the form of the complete template library and a representative training set as inputs to each code, we demonstrate the impact of the assumptions underlying each technique on the output photo-z PDFs. In the absence of a notion of true, unbiased photo-z PDFs, we evaluate and interpret multiple metrics of the ensemble properties of the derived photo-z PDFs as well as traditional reductions to photo-z point estimates. Additionally, we report systematic biases and overall over/under-breadth of the photo-z PDFs of many popular codes, which may indicate avenues for improvement in the algorithms or implementations. Furthermore, we raise attention to the limitations of established metrics for assessing photo-z PDF accuracy; though we identify the conditional density estimate (CDE) loss as a promising metric of photo-z PDF performance in the case where true redshifts are available but true photo-z PDFs are not, we emphasize the need for science-specific performance metrics.

79 ASTRONOMY AND ASTROPHYSICS↗

Long-term modulation of galactic cosmic radiation and its model for space exploration

As the human exploration of space has received new attention in the United States, studies find that exposure to space radiation could adversely impact the mission design. Galactic Cosmic Radiation (GCR), with its very wide range of charges and energies, is particularly important for a mission to Mars, because it imposes a stiff mass penalty for spacecraft shielding. Dose equivalent versus shielding thickness calculations, show a rapid initial drop on exposure with thickness, but an asymptotic behavior at a higher shielding thickness. Uncertainties in the radiobiology are largely unknown. For a fixed radiation risk, this leads to large uncertainties in shielding thickness for small uncertainties in estimated dose. In this paper we investigate the application of steady-state, spherically-symmetric diffusion-convection theory of solar modulation to individual measurements of differential energy spectra from 1954 to 1989 in order to estimate the diffusion coefficient, kappa(r,t), as a function of time. We have correlated the diffusion coefficient to the Climax neutron monitor rates and show that, if the diffusion coefficient can be separated into independent functions of space and time: kappa(r,t)=K(t) k(sub 0) beta Pk(sub 1) (r), where beta is the particle velocity and P the rigidity, then (i) The time dependent quantity 1/K(t), which is proportional to the deceleration potential, phi(r,t), is linearly related to the Climax neutron monitor counting rate. (ii) The coefficients obtained from hydrogen or helium intensity measurements are the same. (iii) There are different correlation functions for odd and even solar cycles. (iv) The correlation function for the Climax neutron monitor counting rate for given time, t, can be used to estimate mean deceleration parameter phi(t) to within +/- 15% with 90% confidence. We have shown that k(r,t) determined from hydrogen and/or helium data, can be used to fit the oxygen and iron differential energy spectra with a root mean square error of about +/- 10%, and essentially independent of the particle charge or energy. We have also examined the ion chamber and C-14 measurements which allow the analysis to be extended from the year 1906 to 1990. Using this model we have defined reference GCR spectra at solar minimum and solar maximum. These can be used for space exploration studies and provide a quantitative estimate of the error in dose due to changes in GCR intensities.

Badhwan, G. D.↗

Chemical Potential of Metal Atoms in Supported Nanoparticles: Dependence upon Particle Size and Support

Metal nanoparticles supported on oxide surfaces form the basis for many industrial catalysts and promise to play an ever-increasing role in future energy and environmental technologies. The chemical potential (μ) of the metal atoms in these particles depends strongly on particle size and support and is an important factor that determines their catalytic properties, including their binding strengths to adsorbed reaction intermediates and their long-term stability against sintering. We present here a method for estimating this chemical potential as a function of particle size for different metal/oxide combinations. We show that this chemical potential for late transition metals is well approximated for a particle of diameter D by μ( D ) = [(3γ m – E adh )(1 + D o / D )](2 V m / D ), where γ m is the surface energy of the bulk metal, E adh is the adhesion energy at the bulk metal/oxide interface, and D o is ~1.5 nm, and V m is the molar volume of the bulk metal. We further show that E adh increases with (1) increasing heat of formation of the most stable oxide of the metal from metal gas atoms plus O 2 (gas) per mole of metal atoms, (2) decreasing enthalpy of reduction of the oxide to its next lower oxidation state plus O 2 (gas), per mole of oxygen atoms, and (3) increasing density of surface oxygen atoms on the oxide surface. The linear scaling of E adh with these properties allows estimations of E adh for a variety of metal/oxide combinations. Using this E adh estimate in the above equation with known values for γ m allows estimates of metal chemical potential versus particle size for late transition metals on various oxide supports. WE conclude this will improve our ability to understand structure–property relations in catalysis and design better catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatial resolution estimation of LANDSAT-4 TM and MSS data

In order to verify that the LANDSAT-4 sensors are operating within specifications, it is useful to estimate the system parameters by analysis of the measured data. One parameter of particular interest is the sensor point-spread function (PSF) which determines the resolution of the system. A method of estimating the PSF is described that utilizes data obtained during scanning of ground elements having identifiable geometric and radiometric structure. These data are then processed in such a manner as to recover either the PSF itself or to estimate the parameters of an assumed functional representation of the PSF.

Mcgillem, C. D.↗

A nonparametric software reliability growth model

Miller and Sofer have presented a nonparametric method for estimating the failure rate of a software program. The method is based on the complete monotonicity property of the failure rate function, and uses a regression approach to obtain estimates of the current software failure rate. This completely monotone software model is extended. It is shown how it can also provide long-range predictions of future reliability growth. Preliminary testing indicates that the method is competitive with parametric approaches, while being more robust.

Miller, Douglas R.↗

Parameter identification for nonlinear aerodynamic systems

Parameter identification for nonlinear aerodynamic systems is examined. It is presumed that the underlying model can be arranged into an input/output (I/O) differential operator equation of a generic form. The algorithm estimation is especially efficient since the equation error can be integrated exactly given any I/O pair to obtain an algebraic function of the parameters. The algorithm for parameter identification was extended to the order determination problem for linear differential system. The degeneracy in a least squares estimate caused by feedback was addressed. A method of frequency analysis for determining the transfer function G(j omega) from transient I/O data was formulated using complex valued Fourier based modulating functions in contrast with the trigonometric modulating functions for the parameter estimation problem. A simulation result of applying the algorithm is given under noise-free conditions for a system with a low pass transfer function.

Pearson, Allan E.↗

Influence of defects on the thermal conductivity of compressed LiF

We report defect formation in LiF, which is used as an observation window in ramp and shock experiments, has significant effects on its transmission properties. Given the extreme conditions of the experiments it is hard to measure the change in transmission directly. Using molecular dynamics, we estimate the change in conductivity as a function of the concentration of likely point and extended defects using a Green-Kubo technique with careful treatment of size effects. With this data, we form a model of the mean behavior and its estimated error; then, we use this model to predict the conductivity of a large sample of defective LiF resulting from a direct simulation of ramp compression as a demonstration of the accuracy of its predictions. Given estimates of defect densities in a LiF window used in an experiment, the model can be used to correct the observations of thermal energy through the window. Also, the methodology we develop is extensible to modeling, with quantified uncertainty, the effects of a variety of defects on the thermal conductivity of solid materials.

36 MATERIALS SCIENCE↗

ASCR Workshop Position Paper: Challenges and Opportunities in High Energy Physics

High energy particle physics and cosmology concern themselves with estimating fundamental parameters of nature, such as the masses and interactions of fundamental particles like the Higgs boson and the rate of expansion of the universe. In doing so, they analyze exabyte-scale datasets, some of the largest in all of science, and face many challenges in subsequent data analysis. These challenges are shared between the two disciplines, but we focus on particle physics to highlight one specific domain. In particle physics, the standard method for estimating parameters involves performing Monte Carlo (MC) integration as a function of both parameters of interest and nuisance parameters using an expensive simulator, counting the number of observed collision events (i.i.d. samples) from an experiment in the corresponding integration domains, and forming a Poisson likelihood function. This likelihood function is then used in a Frequentist manner to construct a maximum likelihood point estimate (MLE) and confidence set for the parameters. To sufficiently populate the high-dimensional integration domains, simulators consume billions of CPU-hours annually and produce hundreds of petabytes of intermediate output data. Several techniques have been developed to: optimize definitions of the integration domains so as to be maximally sensitive to a particular subset of parameters, efficiently estimate the integrals, and build robust surrogate models by interpolating between integral evaluations at different parameter points. One can view this whole endeavor as classical Simulation-Based Inference (SBI).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations

Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations M n+ (M n+ = Li + , Na + , K + , Cu + , Cu 2+ , Co 2+ , Zn 2+ , Ni 2+ , Mg 2+ , Ca 2+ , Ti 2+ , Fe 2+ , Fe 3+ , Co 3+ , Cr 3+ , Ti 4+ and Cr 6+ ) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + M n+ ) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different M n+ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions ( R 2 values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.

Gong, Kai↗

Uncertainty reduction in residual stress measurements by an optimised inverse solution using nonconsecutive polynomials

Many destructive methods for measuring residual stresses such as the slitting method require an inverse analysis to solve the problem. The accuracy of the result as well as an uncertainty component (the model uncertainty) depends on the basis functions used in the inverse solution. The use of a series expansion as the basis functions for the inverse solution was analysed in a previous work for the particular case where functions orders grew consecutively. The present work presents a new estimation of the model uncertainty and a new improved methodology to select the final basis functions for the case where the basis is composed of polynomials. Including nonconsecutive polynomial orders in the basis generates a larger space of possible solutions to be evaluated and allows the possibility to include higher-order polynomials. The paper includes a comparison with two other inverse analyses methodologies applied to synthetically generated data. With the new methodology, the final error is reduced and the uncertainty estimation improved.

36 MATERIALS SCIENCE↗

Navigation for the new millennium: Autonomous navigation for Deep Space 1

The autonomous optical navigation system technology for the Deep Space 1 (DS1) mission is reported on. The DS1 navigation system will be the first to use autonomous navigation in deep space. The systems tasks are to: perform interplanetary cruise orbit determination using images of distant asteroids; control and maintain the orbit of the spacecraft with an ion propulsion system and conventional thrusters, and perform late knowledge updates of target position during close flybys in order to facilitate high quality data return from asteroid MaAuliffe and comet West-Kohoutek-Ikemura. To accomplish these tasks, the following functions are required: picture planning; image processing; dynamical modeling and integration; planetary ephemeris and star catalog handling; orbit determination; data filtering and estimation; maneuver estimation, and spacecraft ephemeris updating. These systems and functions are described and preliminary performance data are presented.

Reidel, J. E.↗

Implicit model following and parameter identification of unstable aircraft

A transformation in the s-plane is described which has utility in implicit model-following optimal control design application and in estimation or parameter identification problems. The objective of the transformation is, for the control problem, to achieve an unstable closed-loop system, and, for the estimation problem, to alleviate algorithm convergence problems that may arise in identifying unstable systems. For the control problem, the transformation is a shift along the real (sigma) axis of the plant and model poles and zeros. This transformation is shown to be equivalent to a modified performance index but offers the advantage of compatibility with existing optimal control solution algorithms. For the estimation problem, the data are multiplied by an exponential function and the assumed measurement and process noise covariances are appropriately modified. Examples of both control and estimation applications are presented.

Lebacqz, J. V.↗

System Identification Applied to Dynamic CFD Simulation and Wind Tunnel Data

Demanding aerodynamic modeling requirements for military and civilian aircraft have provided impetus for researchers to improve computational and experimental techniques. Model validation is a key component for these research endeavors so this study is an initial effort to extend conventional time history comparisons by comparing model parameter estimates and their standard errors using system identification methods. An aerodynamic model of an aircraft performing one-degree-of-freedom roll oscillatory motion about its body axes is developed. The model includes linear aerodynamics and deficiency function parameters characterizing an unsteady effect. For estimation of unknown parameters two techniques, harmonic analysis and two-step linear regression, were applied to roll-oscillatory wind tunnel data and to computational fluid dynamics (CFD) simulated data. The model used for this study is a highly swept wing unmanned aerial combat vehicle. Differences in response prediction, parameters estimates, and standard errors are compared and discussed

Murphy, Patrick C.↗