Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional probability”

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,801 records · Page 100

Shock recovery analogs and the origin of mesosiderites

The origin of mesosiderites, which consist of approximately equal-weight proportions of Fe-Ni metal and silicates (gabbros, basalts, orthopyroxenites, dunites), remains an interesting and complex problem in meteoritics. There is general agreement that multiple impact events were probably involved in the formation of these brecciated stony-iron meteorites, but given the heterogeneity of mesosiderites, additional processes have been invoked to explain the unique and intricate textural and compositional makeup of mesosiderites. We conducted a series of shock recovery experiments to test the impact event(s) scenario. The results indicated significant similarities between the shocked analogs and many mesosiderites. We have compared our analogs with a suite of thin sections of Barea mesosiderite. I have conducted a series of flash heating experiments in which equal-weight proportions of gabbro and stainless steel (SS304) powders were compressed into small charges and heated under reducing conditions for short times. These experiments were used to bracket localized, peak postshock temperatures in our analog shots and to compare the mixing relations between the silicate and metal. The shock recovery experiments used porous metal-silicate powder starting mixtures, therefore our experiments are most analogous to an impact scenario where the target is an asteroidal regolith surface composed of a loose mixture of Fe-Ni metal and heterogeneous silicates. Analog experiments may really describe a secondary impact process similar to the late-stage, localized impact melting event. This leaves one of the crucial questions about mesosiderite genesis unanswered, namely what is the source of the Fe-Ni metal that is so intimately distributed in these meteorites?

Rowan, L. R.↗

The response of the high-latitude dayside ionosphere to an abrupt northward transition in the IMF

We examine the response of the high-latitude ionosphere in the prenoon sector to a northward turning of the IMF. The event was observed in the 11-13 UT interval on June 1, 1987, in the course of a multiday SUNDIAL campaign. The transition in the IMF was observed by the IMP-8 satellite which was located upstream of the earth at a distance of 36 Re. The ionospheric response in the 70-80 deg invariant latitude interval was monitored by two radars. The preexisting plasma convection observed by the radars exhibited large velocities (500-1000 m/s) and stable longterm trends, consistent with the inertial rotation of the convection pattern expected of the conditions then prevailing, B(z) less than 0, B(y) greater than 0. The plasma flow rapidly abated in response to the IMF transition. The electron density measurements made by radar in the meridional plane showed that the ionosphere had been rich in structure with the active deposition of ionization by particle precipitation. Subsequently it resembled an inactive, unstructured mid-latitude configuration. There was a dramatic decrease in the amount of backscatter observed by the HF radar. We analyze the times of transition in the various data sets and show that the ionosphere began to show the effects of the IMF transition about 2 min after its probable arrival at the magnetopause boundary.

Ruohoniemi, J. M.↗

Investigation of Flow Instabilities in the Inlet Ducts of DP-1C VTOL Aircraft

An investigation of flow instabilities in the inlet ducts of a two-engine vertical takeoff and landing aircraft DP-1C is described in this report. Recent tests revealed that the engines stall during run ups while the aircraft is operating on the ground. These pop stalls occurred at relatively low power levels, sometimes as low as 60 percent of the engine full speed. Inability to run the engines up to the full speed level is attributed to in-ground effects associated with hot gas ingestion. Such pop stalls were never experienced when the aircraft was tested on a elevated grid platform, which ensured that the aircraft was operating in out-of-the-ground-effect conditions. Based on available information on problems experienced with other vertical takeoff and landing aircraft designs, it was assumed that the engine stalls were caused by partial ingestion of hot gases streaming forward from the main exit nozzle under the aircraft inlets, which are very close to the ground. It was also suggested that the nose wheel undercarriage, located between the inlets, may generate vortices or an unstable wake causing intense mixing of hot exit gases with incoming inlet flow, which would enhance the hot gas ingestion. After running a short three-day series of tests with fully instrumented engine inlets, it is now believed the most probable reason for engine pop stalls are random ingestions of a vortex generated between the two streams moving in opposite directions: outbound hot gas stream from the main nozzle close to the ground and inbound inlet flow above. Originally, the vortex is in a horizontal plane. However, at a certain velocity ratio of these two streams, the vortex attaches either to the ground or the aircraft surface at one end and the other end is swallowed by one of the aircraft inlets. Once the vortex enters the inlet duct, a puff of hot air can be sucked through the vortex core into the engine, which causes a serious inlet flow field distortion followed by an engine stall. Once the engine stalls, the outflow from the inlet pushes the vortex away and the engine resumes normal operation. This hypothesis needs to be verified experimentally; e.g., by extensive smoke flow visualization ahead of the aircraft inlets.

Lepicovsky, Jan↗

Statistical Issues for Uncontrolled Reentry Hazards

A number of statistical tools have been developed over the years for assessing the risk of reentering objects to human populations. These tools make use of the characteristics (e.g., mass, shape, size) of debris that are predicted by aerothermal models to survive reentry. The statistical tools use this information to compute the probability that one or more of the surviving debris might hit a person on the ground and cause one or more casualties. The statistical portion of the analysis relies on a number of assumptions about how the debris footprint and the human population are distributed in latitude and longitude, and how to use that information to arrive at realistic risk numbers. This inevitably involves assumptions that simplify the problem and make it tractable, but it is often difficult to test the accuracy and applicability of these assumptions. This paper looks at a number of these theoretical assumptions, examining the mathematical basis for the hazard calculations, and outlining the conditions under which the simplifying assumptions hold. In addition, this paper will also outline some new tools for assessing ground hazard risk in useful ways. Also, this study is able to make use of a database of known uncontrolled reentry locations measured by the United States Department of Defense. By using data from objects that were in orbit more than 30 days before reentry, sufficient time is allowed for the orbital parameters to be randomized in the way the models are designed to compute. The predicted ground footprint distributions of these objects are based on the theory that their orbits behave basically like simple Kepler orbits. However, there are a number of factors - including the effects of gravitational harmonics, the effects of the Earth's equatorial bulge on the atmosphere, and the rotation of the Earth and atmosphere - that could cause them to diverge from simple Kepler orbit behavior and change the ground footprints. The measured latitude and longitude distributions of these objects provide data that can be directly compared with the predicted distributions, providing a fundamental empirical test of the model assumptions.

Matney, Mark↗

Microstructural and Isotopic Constraints on WL Rim Formation

Coordinated microanalyses of Wark-Lovering (WL) rims are needed to best understand their origin and to decipher their subsequent evolution both in the nebular and parent body settings. Here we present the mineralogy, petrology, microstructures, O isotopic compositions, and Al-Mg systematics of a WL rim on a Type B CAI, ‘‘Big Guy”, from the reduced CV3 chondrite Vigarano [1]. Our SEM and TEM study reveals seven distinct mineral layers in the WL rim that include: (1) gehlenite with rare grossite, (2) hibonite, (3) spinel with minor hibonite and perovskite, (4) zoned melilite (Åk(sub ~0-10)), (5) anorthite, (6) zoned diopside grading outwards from Al,Ti-rich to Al,Tipoor, and (7) forsterite intergrown with diopside. We infer a two-stage history in which WL rim formation was initiated by flash melting and extensive evaporation of the original inclusion edge, followed by subsequent condensation under highly dynamic conditions. The outermost edge of the CAI mantle is mineralogically and texturally distinct compared to the underlying mantle that is composed of coarse, zoned melilite (Åk(sub ~10-60)) grains. The mantle edge contains finegrained gehlenite with hibonite and rare grossite that likely formed by rapid crystallization from a Ca,Al-rich melt produced during a flash vaporization event [2]. These gehlenite and hibonite layers are surrounded by successive layers of spinel, melilite, diopside, and forsterite, indicating their sequential gas-solid reactions onto hibonite. Anorthite occurs as a discontinuous layer that corrodes adjacent melilite and Al-diopside, and appears to have replaced them [3,4], probably even later than the forsterite layer formation. All the WL rim minerals analyzed using the JSC NanoSIMS 50L are 16O-rich (Δ17O ≈ −23‰), indicating their formation in an 16O-rich gas reservoir. Our data are in contrast with many CV CAIs that show heterogeneous Δ 17O values across their WL rims [5]. Our Al-Mg data obtained using the UCLA ims-1290 ion microprobe of the CAI interior and the WL rim define a well-correlated isochron with (26Al/27Al)(sub 0) = 4.94 × 10(exp −5), indicating their synchronous formation ≤5 × 10(exp 4) years after the canonical CAI value. In contrast, no 26Mg excesses are observed in the WL rim anorthite, which suggests its later formation or later isotopic resetting in an 16O-rich gas reservoir, after 26Al had decayed.

Han, J.↗

Development of a new control rod drive mechanism design for the ISU AGN-201M reactor

The Aerojet General Nucleonics (AGN) model 201-Modified, known as the AGN-201M reactor, plays an essential role in the educational and research activities at Idaho State University (ISU). The ISU AGN-201M's original Control Rod Drive Mechanism (CRDM) has been in operation for more than fifty years with no large-scale redesigns. The CRDM is required to eject the fuel rods within one second during a SCRAM event (also known as a 'reactor trip') and adjust the control rods' insertion speed, and keeps the rod insertion sequence correct. The existing control rod drive mechanisms meet these criteria but experience a few concerns due to aging. Concerns include complex maintenance and costly repairs for old electromechanical components, rod position feedback errors, and the impediment of the plate during a SCRAM due to binding of the lead screws of the existing mechanism. During a binding event, the drive mechanism becomes locked, preventing the control rod's magnetic plate from moving in or out under the reactor's normal and emergency operating conditions. Although, the binding has no effect on the ability of the rod to exit the core during the SCRAM. To counteract the concerns and issues with the current CRDM, a new design has been proposed using newer components and a simplified design. The new design utilizes more advanced electric and mechanical components that are commercially available. The new CRDM system is divided into four main aspects: (1) control rod movement design (motor, lead screw, guide rods), (2) control rod ejection (springs, electromagnet), (3) control rod position and feedback (position transducer, microswitches), and (4) material selection and structural analysis. The new design aims to reduce the overall complexity and probability of failure to improve the reactor's overall reliability. With proper material selection and improved structural design, the new drives are lighter with little to no change in structural integrity. The new control rod drive mechanism eliminates binding scenarios by using a single lead screw and implementing additional guide rods. An advanced linear position sensor and microswitches replace the existing and aging synchro system for accurate rod position feedback resulting in better reactivity control. The new design meets the reactor's operational limits by having an average reactivity insertion of 0.065% Δk/k per second, which corresponds to a total control rod insertion time of 19.23 s, while the control rod's ejection time remains less than one second during a SCRAM event. The new design ensures the reactor's long-term viability for educational and research activities by increasing the reliability and safety of operation for years to come.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Bayesian inference of nuclear-matter density from proton scattering

Background: Proton elastic scattering at intermediate energy is widely employed as a tool for determining the matter radius of atomic nuclei. Here, the sensitivity of the approach relies on high-resolution measurements at small scattering angles and low-momentum transfer. Under these conditions, the Glauber multiple scattering theory accurately describes the proton-nucleus elastic cross section. Purpose: Investigate the sensitivity of the Glauber multiple scattering theory to uncertainties associated with input parameters such as the nuclear-matter density distribution and nucleon-nucleon data. Method: A joint Bayesian inference was performed using 12 angular distributions of elastic scattering at different energies on 58 Ni, 90 Zr, and 208 Pb targets. A Metropolis-Hastings algorithm was implemented to make an uncertainty quantification analysis for the input parameters used in the Glauber multiple scattering theory. Results: The experimental cross sections were fitted simultaneously using a joint Bayesian inference approach. Posterior probability density distributions of 42 input parameters were obtained from the analysis. A moderate correlation between the nuclear density parameters and the nucleon-nucleon cross sections was found. This correlation impacts the extraction of the nuclear-matter radius. Conclusions: The present analysis provided a consistent method for extracting the nuclear-matter density distribution of 58 Ni, 90 Zr, and 208 Pb from data across different incident energies. Due to the correlation of the nucleon-nucleon cross sections with the other input parameters, a constrained Bayesian inference using free nucleon-nucleon cross section data was performed. The nuclear-matter radii obtained from the analysis are in good agreement with multiple results reported in the literature.

190 ≤ A ≤ 219↗

Performance of candidate SEPS solar cells as a function of low temperature and low intensity exposure and 1 MeV electron irradiation

The behavior of 144 high-performance shallow-junction silicon solar cells under conditions of low temperature and intensity was examined. The cells represented nine combinations of thickness, base resistivity, front surface texture, and rear surface treatment. At least 16 cells of each type were individually tested both as active elements under light levels between 0.04 and 1.0 solar constant and in the dark as passive rectifiers under external forward bias. One cell type was also irradiated with 1 MeV electrons up to 2.7 x 10 to the 15th e/sq cm. The inferior cells demonstrated high ohmic and/or non-ohmic shunting. No series resistance or Schottky barrier effects were observed. Fluences beyond 10 to the 13th e/sq cm lowered cell current, probably by increasing volume recombination.

Whitaker, A. F.↗

Cell Maintenance Systems

Living human cells require attachment to a suitable surface and special culture conditions in order to grow. These requirements are modified and amplified when cells are taken into a weightless environment. Special handling and maintenance systems are required for routine laboratory procedures conducted in the Orbiter and in the Spacelab. Methods were developed to maintain cells in special incubators designed for the Orbiter middeck, however, electrophoresis and other experiments require cells to be harvested off of the culture substrate before they can be processed or used. The cell transport assembly (CTA) was flown on STS-8, and results show that improvements are required to maintain adequate numbers of cells in this device longer than 48 hours. The life sciences middeck centrifuge probably can be used, but modifications will be required to transfer cells from the CTA and keep the cells sterile. Automated systems such as the Skylab SO-15 flight hardware and crew operated systems are being evaluated for use on the Space Shuttle, Spacelab, and Space Station research modules.

Morrison, D. R.↗

Concentration transient analysis of antimony surface segregation during Si(100) molecular beam epitaxy

Antimony surface segregation during Si(100) molecular beam epitaxy (MBE) was investigated at temperatures T(sub s) = 515 - 800 C using concentration transient analysis (CTA). The dopant surface coverage Theta, bulk fraction gamma, and incorporation probability sigma during MBE were determined from secondary-ion mass spectrometry depth profiles of modulation-doped films. Programmed T(sub s) changes during growth were used to trap the surface-segregated dopant overlayer, producing concentration spikes whose integrated area corresponds to Theta. Thermal antimony doping by coevaporation was found to result in segregation strongly dependent on T(sub s) with Theta(sub Sb) values up to 0.9 monolayers (ML): in films doped with Sb(+) ions accelerated by 100 V, Theta(sub Sb) was less than or equal to 4 x 10(exp -3) ML. Surface segregation of coevaporated antimony was kinematically limited for the film growth conditions in these experiments.

Markert, L. C.↗

Radio-Optical Alignments in a Low Radio Luminosity Sample

We present an optically-based study of the alignment between the radio axes and the optical major axes of eight z approximately 0.7 radio galaxies in a 7C sample. The radio galaxies in this sample are approximately 20-times less radio luminous than 3C galaxies at the same redshift, and are significantly less radio-luminous than any other well-defined samples studied to date. Using Nordic Optical Telescope images taken in good seeing conditions at rest-frame wavelengths just longward of the 4000A break, we find a statistically significant alignment effect in the 7C sample. Furthermore, in two cases where the aligned components are well separated from the host we have been able to confirm spectroscopically that they are indeed at the same redshift as the radio galaxy. However, a quantitative analysis of the alignment in this sample and in a corresponding 3C sample from HST (Hubble Space Telescope) archival data indicates that the percentage of aligned flux may be lower and of smaller spatial scale in the 7C sample. Our study suggests that alignments on the 50-kpc scale are probably closely related to the radio luminosity, whereas those on the 15 kpc scale are not. We discuss these results in the context of popular models for the alignment effect.

Lacy, Mark↗

SafeMAP: Safe Multi-Agent Planning Framework Based on Dynamic Probabilistic Risk Assessment

This paper proposes a risk-aware framework for Safe Multi-Agent Planning (SafeMAP) that unifies disparate models for multi-agent systems in a Markovian process that allows for simultaneous system health monitoring, decision making under uncertainty, and multi-agent system collaboration. As operations beyond low earth orbit mature, there is an increased need for autonomous cyber-physical systems with onboard decision making capabilities. Multi-agent cyber-physical systems in particular offer the potential of increased efficiency, resiliency, and mission capabilities for future applications such as multi-rover terrain operations, distributed satellite operations, and management of smart lunar habitats. SafeMAP utilizes physics-based models of each agent and the relevant components, probability models of the environment and component operational states, and reward models for mission-specific objectives such as scientific task completion or resource consumption. The output of SafeMAP is a set of mission plans that satisfy the mission objective under specified risk/reward constraints. A readable interpretation of each of these generated mission plans is provided as an additional output. SafeMAP has been demonstrated on a simulated case study involving a four-rover system performing surface mapping operations and science tasks. Results of this paper demonstrate SafeMAP’s ability to generate explainable mission plans that satisfy the mission objective while minimizing risk under nominal and off-nominal conditions.

Mohammad Hejase↗

SafeMAP: Safe Multi-Agent Planning framework based on Dynamic Probabilistic Risk Assessment

This paper proposes a risk-aware framework for Safe Multi-Agent Planning (SafeMAP) that unifies disparate models for multi-agent systems in a Markovian process that allows for simultaneous system health monitoring, decision making under uncertainty, and multi-agent system collaboration. As operations beyond low earth orbit mature, there is an increased need for autonomous cyber-physical systems with onboard decision making capabilities. Multi-agent cyber-physical systems in particular offer the potential of increased efficiency, resiliency, and mission capabilities for future applications such as multi-rover terrain operations, distributed satellite operations, and management of smart lunar habitats. SafeMAP utilizes physics-based models of each agent and the relevant components, probability models of the environment and component operational states, and reward models for mission-specific objectives such as scientific task completion or resource consumption. The output of SafeMAP is a set of mission plans that satisfy the mission objective under specified risk/reward constraints. A readable interpretation of each of these generated mission plans is provided as an additional output. SafeMAP has been demonstrated on a simulated case study involving a four-rover system performing surface mapping operations and science tasks. Results of this paper demonstrate SafeMAP’s ability to generate explainable mission plans that satisfy the mission objective while minimizing risk under nominal and off-nominal conditions.

Mohammad Hejase↗

Deep Convective Cloud Calibration Sensitivity Studies in Support of Radiometrically Scaling GEO Imagers With VIIRS

The NASA CERES SYN1deg product provides the scientific community regional hourly TOA and surface broadband fluxes and clouds. For consistent geostationary (GEO) derived fluxes and clouds the GEO imagers are radiometrically scaled to the Aqua-MODIS calibration reference. The CERES project utilizes GEO and MODIS or VIIRS analogous channel coincident, collocated, and co-angled radiance pairs as the primary method to inter-calibrate the GEO imagers. Tropical deep convective clouds (DCC) are bright, near Lambertian, top of atmosphere pseudo invariant Earth targets that do not rely on coincident ray-matched radiance pairs to radiometrically scale sensors to a common calibration reference. The DCC invariant target (DCC-IT) methodology collectively analyzes all tropical DCC identified pixel radiances by way of probability density function (PDF) distributions. Perfectly inter-calibrated sensor pairs should reveal nearly identical PDF distributions given the same DCC identification criterion. The PDF median, mean, mode, and inflection point statistics were tested as a function of DCC identification criterion using SNPP-VIIRS and Himawari-8 AHI 0.65μm channel radiances during January 2019. It was found that the PDF inflection point provided inter-calibration factors within 0.25% that were nearly independent of DCC identification criterion. The PDF median provided inter-calibration factors within 0.25% for the coldest BT and most stringent homogeneity factors. The PDF mean and mode statistics were inadequate under any DCC conditions. It is critical for the DCC pixel radiances to be anisotropically corrected. The DCC-IT methodology will also be tested for other visible and SWIR bands.

DCC↗

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗

Environmentally Adaptive, Multiband Software-Defined-Radar for Monitoring of Item of Interest in a Dynamically Cluttered Room

High-resolution multi-frequency radars when networked together, form a robust and seamless monitoring system for high-valued items. The multi-look-direction radar network investigated by LLNL and DSI, operates at 77 GHz and 5.8GHz, and is capable of detection at centimeter resolution movement of objects of interest. This dual frequency radar (unlike optical methods) can operate under various environments with conditions such a smoke, dust, and other natural or manmade obscurations. In addition, the multiband network radar is based on low-cost FCC approved EM specifications and can be scaled with the dimensions of the room and objects under consideration by user-defined inputs. Furthermore, the underlying radar signal processing based on software-defined radar engine is adaptive to its operational environment, where operators can input the size of the facility and object of interest and the distance between the object and radar. Finally, this multi-frequency radar employs a host of change detection and machine learning algorithms to reach performance levels in terms of high probability of detection and low false alarm rates while the monitoring state of the object can be reported remotely.

42 ENGINEERING↗

Analysis of the Westland Data Set

The "Westland" set of empirical accelerometer helicopter data with seeded and labeled faults is analyzed with the aim of condition monitoring. The autoregressive (AR) coefficients from a simple linear model encapsulate a great deal of information in a relatively few measurements; and it has also been found that augmentation of these by harmonic and other parameters call improve classification significantly. Several techniques have been explored, among these restricted Coulomb energy (RCE) networks, learning vector quantization (LVQ), Gaussian mixture classifiers and decision trees. A problem with these approaches, and in common with many classification paradigms, is that augmentation of the feature dimension can degrade classification ability. Thus, we also introduce the Bayesian data reduction algorithm (BDRA), which imposes a Dirichlet prior oil training data and is thus able to quantify probability of error in all exact manner, such that features may be discarded or coarsened appropriately.

Wen, Fang↗

A Revised Approach to the ULDB Design

The National Aeronautics and Space Administration Balloon Program has experienced problems in the scaling up of the proposed Ultra Long Duration Balloon. Full deployment of the balloon envelope has been the issue for the larger balloons. There are a number of factors that contribute to this phenomenon. Analytical treatments of the deployment issue are currently underway. It has also been acknowledged that the current fabrication approach using foreshortening is costly, labor intensive, and requires significant handling during production thereby increasing the chances of inducing damage to the envelope. Raven Industries has proposed a new design and fabrication approach that should increase the probability of balloon deployment, does not require foreshortening, will reduce the handling, production labor, and reduce the final balloon cost. This paper will present a description of the logic and approach used to develop this innovation. This development consists of a serial set of steps with decision points that build upon the results of the previous steps. The first steps include limited material development and testing. This will be followed by load testing of bi-axial reinforced cylinders to determine the effect of eliminating the foreshortening. This series of tests have the goal of measuring the strain in the material as it is bi-axially loaded in a condition that closely replicated the application in the full-scale balloon. Constant lobe radius pumpkin shaped test structures will be designed and analyzed. This matrix of model tests, in conjunction with the deployment analyses, will help develop a curve that should clearly present the deployment relationship for this kind of design. This will allow the "design space" for this type of balloon to be initially determined. The materials used, analyses, and ground testing results of both cylinders and small pumpkin structures will be presented. Following ground testing, a series of test flights, staged in increments of increasing suspended load and balloon volume, will be conducted. The first small scale test flight has been proposed for early Spring 2004. Results of this test flight of this new design and approach will presented. Two additional domestic test flights from Ft. Sumner, New Mexico, and Palestine, Texas, and one circumglobal test flight from Australia are planned as part of this development. Future plans for both ground testing and test flights will also be presented.

Smith, Michael↗

Refine your search to explore more results.