Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical validation”

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 757 records · Page 42

Qualitative and Quantitative Evaluation for Representative Human Reliability Analysis Methods

The Korea Institute of Nuclear Safety (KINS) is the regulatory expert organization established by the Korean government to strengthen the nation’s technical capabilities relating to nuclear safety regulation. KINS oversees the technical aspects of nuclear safety regulation, including safety reviews, inspections, education, and safety research—all conducted based on technical knowledge and accumulated regulatory experience. In 2023, KINS requested that Idaho National Laboratory (INL) validates representative human reliability analysis (HRA) methods used throughout the world, thus affording KINS with a basis for determining an HRA method adequate for its domestic regulatory purposes. The present paper mainly examines INL’s efforts in this regard. The resulting INL study covered four representative HRA methods widely used by nuclear utilities and regulatory institutes. These methods were qualitatively evaluated by applying specific evaluation criteria and determining how well each method reflected critical HRA issues. For this assessment, INL benchmarked the Halden International HRA Empirical Study. Using the Halden empirical data, along with information on human failure events (HFEs), the present study employed the selected HRA methods to estimate human error probabilities (HEPs) for the HFEs. It also performed statistical analyses to compare the HEPs predicted via the HRA methods against those from the Halden empirical data.

99 - GENERAL AND MISCELLANEOUS↗

MAS-CCD: New technique for measuring low-level charge content based on the multiple amplifier architecture

Low-noise detectors are a key technology for the next generation of astronomical instruments aimed at spectroscopy of faint objects and the search for exoplanets. In this context, the multiple-amplifier sensing charge-coupled device (MAS-CCD) emerges as a promising technology for future scientific instruments. A critical parameter affecting the performance of these devices is spurious charge, produced by the clocking of the gates. Its measurement is typically challenging with existing methods. In practice, the optimization of this parameter often relies on empirical procedures that require significant time and careful consideration of the trade-off with full-well capacity. In this work, we present a new technique to estimate spurious charge based on covariance analysis of the output amplifiers of the MAS-CCD, which measures the same charge packet in different amplifiers at different times. The method enables fast and precise measurements of spurious charge under operating conditions where conventional approaches are difficult to apply. We develop the theoretical framework of the method and validate the model through simulations. The results demonstrate the feasibility of this approach and suggest that it could serve as a basis for reliable large-scale characterization of sensor performance.

Gamero, Miqueas E. [Buenos Aires, CONICET; Fermila↗

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗

Charting the chemical space of Zintl phases with graph neural networks and bonding insights

A large number of Zintl phases have been discovered by solid-state chemists driven by empirical knowledge, chemical intuition and in some cases, through serendipitous accidents. These discoveries have only scratched the surface, given the vast compositional and structural diversity that Zintl phases can accommodate. The large chemical space of Zintl phases, as well as intermetallic compounds in general, remain under-explored. Here, we use graph neural networks and the upper bound energy minimization approach to efficiently scan a large chemical space of >90 000 hypothetical Zintl phases and accurately discover 1810 new thermodynamically stable phases with 90% precision, as validated with first-principles calculations. We show that our approach is more than 2× more accurate in predicting DFT stability than M3GNet (40% precision) on the same dataset. Using a random forest model and SHAP analysis, we demonstrate the critical role of ionic bonding in the thermodynamic stability of Zintl phases. Our results not only expand the known chemical landscape of Zintl phases but also highlight the efficacy of machine learning frameworks combined with domain knowledge in uncovering chemically meaningful insights across complex intermetallics.

36 MATERIALS SCIENCE↗

Principles in Remote Sensing of Aerosol from MODIS Over Land and Ocean

The well-calibrated spectral radiances measured by MODIS will be processed to retrieve daily aerosol properties that include optical thickness and mass loading over land and optical thickness, the mean particle size of the dominant mode and the ratio between aerosol modes over ocean. In addition, after launch, aerosol single scattering albedo will be calculated as an experimental product. The retrieval process over land is based on a dark target method that identifies appropriate targets in the mid-IR channels and uses an empirical relationship found between the mid-ER and the visible channels to estimate surface reflectance in the visible from the mid-HZ reflectance measured by satellite. The method employs new aerosol models for industrial, smoke and dust aerosol. The process for retrieving aerosol over the ocean makes use of the wide spectral band from 0.55-2.13 microns and a look-up table constructed from combinations of five accumulation modes and five coarse modes. Both the over land and over ocean algorithms have been validated with satellite and airborne radiance measurements. We estimate that MODIS will be able to measure aerosol optical thickness (t) to within 0.05 +/- 0.2t over land and to within 0.05 +/- 0.05t over ocean. Much of the earth's surface is located far from aerosol sources and experience very low aerosol optical thickness. Will the accuracy expected from MODIS retrievals be sufficient to measure the global aerosol direct and indirect forcing? We are attempting to answer this question using global model results and cloud climatology.

Remer, L. A.↗

Improved Treatment of Host-galaxy Correlations in Cosmological Analyses with Type Ia Supernovae

Improving the use of Type Ia supernovae (SNe Ia) as standard candles requires a better approach to incorporate the relationship between SNe Ia and the properties of their host galaxies. Using a spectroscopically confirmed sample of ~1600 SNe Ia, we develop the first empirical model of underlying populations for SNe Ia light-curve properties that includes their dependence on host-galaxy stellar mass; we find a significant correlation between stretch population and stellar mass (99.9% confidence) and a weaker correlation between color and stellar mass (90% confidence). These populations are important inputs to simulations that are used to model selection effects and correct distance biases within the BEAMS with Bias Correction (BBC) framework. Here we improve BBC to also account for SNe Ia-host correlations, and we validate this technique on simulated data samples. Here, we recover the input relationship between SNe Ia luminosity and host-galaxy stellar mass (the mass step, γ) with a bias of 0.004 ±0.001 mag, which is a factor of 5 improvement over previous methods that have a γ bias of ~0.02 ± 0.001 mag. We adapt BBC for a novel dust-based model of intrinsic brightness variations, which results in a greatly reduced mass step for data (γ = 0.017 ± 0.008) and for simulations (γ = 0.006 ± 0.007). Analyzing simulated SNe Ia, the biases on the dark energy equation of state, w, vary from Δw = 0.006(5) to 0.010(5) with our new BBC method; these biases are significantly smaller than the 0.02(5) w bias using previous BBC methods that ignore SNe Ia-host correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

Labyrinth Seal Flutter Analysis and Test Validation in Support of Robust Rocket Engine Design

High energy-density turbomachines, like the SSME turbopumps, utilize labyrinth seals, also referred to as knife-edge seals, to control leakage flow. The pressure drop for such seals is order of magnitude higher than comparable jet engine seals. This is aggravated by the requirement of tight clearances resulting in possible unfavorable fluid-structure interaction of the seal system (seal flutter). To demonstrate these characteristics, a benchmark case of a High Pressure Oxygen Turbopump (HPOTP) outlet Labyrinth seal was studied in detail. First, an analytical assessment of the seal stability was conducted using a Pratt & Whitney legacy seal flutter code. Sensitivity parameters including pressure drop, rotor-to-stator running clearances and cavity volumes were examined and modeling strategies established. Second, a concurrent experimental investigation was undertaken to validate the stability of the seal at the equivalent operating conditions of the pump. Actual pump hardware was used to construct the test rig, also referred to as the (Flutter Rig). The flutter rig did not include rotational effects or temperature. However, the use of Hydrogen gas at high inlet pressure provided good representation of the critical parameters affecting flutter especially the speed of sound. The flutter code predictions showed consistent trends in good agreement with the experimental data. The rig test program produced a stability threshold empirical parameter that separated operation with and without flutter. This empirical parameter was used to establish the seal build clearances to avoid flutter while providing the required cooling flow metering. The calibrated flutter code along with the empirical flutter parameter was used to redesign the baseline seal resulting in a flutter-free robust configuration. Provisions for incorporation of mechanical damping devices were introduced in the redesigned seal to ensure added robustness

El-Aini, Yehia↗

Conservational PDF Equations of Turbulence

Recently we have revisited the traditional probability density function (PDF) equations for the velocity and species in turbulent incompressible flows. They are all unclosed due to the appearance of various conditional means which are modeled empirically. However, we have observed that it is possible to establish a closed velocity PDF equation and a closed joint velocity and species PDF equation through conditions derived from the integral form of the Navier-Stokes equations. Although, in theory, the resulted PDF equations are neither general nor unique, they nevertheless lead to the exact transport equations for the first moment as well as all higher order moments. We refer these PDF equations as the conservational PDF equations. This observation is worth further exploration for its validity and CFD application

Shih, Tsan-Hsing↗

Cure Cycle Design Methodology for Fabricating Reactive Resin Matrix Fiber Reinforced Composites: A Protocol for Producing Void-free Quality Laminates

For the fabrication of resin matrix fiber reinforced composite laminates, a workable cure cycle (i.e., temperature and pressure profiles as a function of processing time) is needed and is critical for achieving void-free laminate consolidation. Design of such a cure cycle is not trivial, especially when dealing with reactive matrix resins. An empirical "trial and error" approach has been used as common practice in the composite industry. Such an approach is not only costly, but also ineffective at establishing the optimal processing conditions for a specific resin/fiber composite system. In this report, a rational "processing science" based approach is established, and a universal cure cycle design protocol is proposed. Following this protocol, a workable and optimal cure cycle can be readily and rationally designed for most reactive resin systems in a cost effective way. This design protocol has been validated through experimental studies of several reactive polyimide composites for a wide spectrum of usage that has been documented in the previous publications.

Hou, Tan-Hung↗

Inference of the optical depth to reionization τ from Planck CMB maps with convolutional neural networks

The optical depth to reionization, τ, is the least constrained parameter of the cosmological Λ cold dark matter (ΛCDM) model. To date, its most precise value is inferred from large-scale polarized cosmic microwave background (CMB) power spectra from the High Frequency Instrument (HFI) aboard the Planck satellite. These maps are known to contain significant contamination by residual non-Gaussian systematic effects, which are hard to model analytically. Therefore, robust constraints on τ are currently obtained through an empirical cross-spectrum likelihood built from simulations. In this paper, we present a likelihood-free inference of τ from polarized Planck HFI maps which, for the first time, is fully based on neural networks (NNs). NNs have the advantage of not requiring an analytical description of the data and can be trained on state-of-the-art simulations, combining the information from multiple channels. By using Gaussian sky simulations and Planck SRoll2 simulations, including CMB, noise, and residual instrumental systematic effects, we trained, tested, and validated NN models considering different setups. We inferred the value of τ directly from Stokes Q and U maps at ~4° pixel resolution, without computing angular power spectra. On Planck data, we obtained τ NN = 0.0579 ± 0.0082, which is compatible with current EE cross-spectrum results but with a ~30% larger uncertainty, which can be assigned to the inherent nonoptimality of our estimator and to the retraining procedure applied to avoid biases. While this paper does not improve on current cosmological constraints on τ, our analysis represents a first robust application of NN-based inference on real data, and highlights its potential as a promising tool for complementary analysis of near-future CMB experiments, also in view of the ongoing challenge to achieve the first detection of primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗

A Modeling Toolkit for Comparing AC and DC Electrical Distribution Efficiency in Buildings

Recently, there has been considerable research interest in the potential for DC distribution systems in buildings instead of the traditional AC distribution systems. Due to the need for performing power conversions between DC and AC electricity, DC distribution may provide electrical efficiency advantages in some systems. To support comparative evaluations of AC-only, DC-only, and hybrid AC/DC distribution systems in buildings, a new modeling toolkit called the Building Electrical Efficiency Analysis Model (BEEAM) was developed and is described in this paper. To account for harmonics in currents or voltages arising from nonlinear devices, the toolkit implements harmonic power flow, along with nonlinear device behavioral descriptions derived from empirical measurements. This paper describes the framework, network equations, device representations, and an implementation of the toolkit in an open source software package, including a component library and graphical interface for creating circuits. Simulations of electrical behavior and device and system efficiencies using the toolkit are compared with experimental measurements of a small office environment in a variety of operating and load configurations. A detailed analysis of uncertainty estimation is also provided. Key findings were that a comparison of predicted versus measured efficiencies and power losses in the validation testbed using the initial toolkit implementation predicted device- and system-level efficiencies with reasonably good accuracy under both balanced and unbalanced AC scenarios. An uncertainty analysis also revealed that the maximum estimated error for system efficiency across all scenarios was 3%, and measured and modeled system efficiency agreed within the experimental uncertainty in approximately half of the scenarios. Based on the correspondence between simulation and measurement, the toolkit is proposed by the authors as a potentially useful tool for comparing efficiency in AC, DC, and hybrid AC/DC distribution systems in buildings.

DC distribution↗

On Electromagnetic Field-to-Wire Coupling Versus Conducted Injection Techniques

Since the inception of conducted injection techniques to model radiated susceptibility/immunity coupling, considerable debate has ensued regarding its validity. This paper affirms the viewpoint of Szentkuti, (1989) builds upon test results of Adams (1992) and Trout (1996), and discusses Perini's theoretical observations (1993, 1995A, 1995B). Analytical and test results are presented which further demonstrate under what specific conditions conducted and radiated techniques can be correlated, and how the work of Adams, Trout, and Perini fits into the general problem of modeling field-to-wire coupling. At frequencies where transmission line and antenna effects are minimal, conducted immunity techniques provide excellent correlation with analytical and empirical predictions of radiated coupling. From a practical standpoint, conducted injection techniques provide realistic coupling at frequencies and amplitude levels that would be uneconomical to achieve with traditional radiated techniques.

Javor, Ken↗

Supersonic Stall Flutter of High Speed Fans

An analytical model is developed for predicting the onset of supersonic stall bending flutter in axial flow compressors. The analysis is based on a modified two dimensional, compressible, unsteady actuator disk theory. It is applied to a rotor blade row by considering a cascade of airfoils whose geometry and dynamic response coincide with those of a rotor blade element at 85 percent of the span height (measured from the hub). The rotor blades are assumed to be unshrouded (i.e., free standing) and to vibrate in their first flexural mode. The effects of shock waves and flow separation are included in the model through quasi-steady, empirical, rotor total-pressure-loss and deviation-angle correlations. The actuator disk model predicts the unsteady aerodynamic force acting on the cascade blading as a function of the steady flow field entering the cascade and the geometry and dynamic response of the cascade. Calculations show that the present model predicts the existence of a bending flutter mode at supersonic inlet Mach numbers. This flutter mode is suppressed by increasing the reduced frequency of the system or by reducing the steady state aerodynamic loading on the cascade. The validity of the model for predicting flutter is demonstrated by correlating the measured flutter boundary of a high speed fan stage with its predicted boundary. This correlation uses a level of damping for the blade row (i.e., the log decrement of the rotor system) that is estimated from the experimental flutter data. The predicted flutter boundary is shown to be in good agreement with the measured boundary.

Adamczyk, J. J.↗

The Charge Transfer Efficiency and Calibration of WFPC2

A new determination of WFPC2 photometric corrections is presented, using HSTphot reduction of the WFPC2 Omega Centauri and NGC 2419 observations from January 1994 through March 2000 and a comparison with ground-based photometry. No evidence is seen for any position-independent photometric offsets (the "long-short anomaly"); all systematic errors appear to be corrected with the CTE and zero point solution. The CTE loss time dependence is determined to be very significant in the Y direction, causing time-independent CTE solutions to be valid only for a small range of times. On average, the present solution produces corrections similar to Whitmore, Heyer, & Casertano, although with an improved functional form that produces less scatter in the residuals and determined with roughly a year of additional data. In addition to the CTE loss characterization, zero point corrections are also determined as functions of chip, gain, filter, and temperature. Of interest, there are chip-to-chip differences of order 0.01 - 0.02 magnitudes relative to the Holtzman et al. calibrations, and the present study provides empirical zero point determinations for the non-standard filters such as the frequently-used F450W, F606W, and F702W.

Dolphin, Andrew E.↗

Validation of High-Resolution CFD Method for Slosh Damping Extraction of Baffled Tanks

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involve the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations and knowledge were mainly carried out by extensive experimental studies. A Volume-Of-Fluid (VOF) based CFD program developed at NASA MSFC was applied to extract slosh damping in a baffled tank from the first principle. First, experimental data using water with subscale smooth wall tank were used as the baseline validation. CFD simulation was demonstrated to be capable of accurately predicting natural frequency and very low damping value from the smooth wall tank at different fill levels. The damping due to a ring baffle at different liquid fill levels from barrel section and into the upper dome was then investigated to understand the slosh damping physics due to the presence of a ring baffle. Based on this study, the Root-Mean-Square error of our CFD simulation in estimating slosh damping was less than 4.8%, and the maximum error was less than 8.5%. Scalability of subscale baffled tank test using water was investigated using the validated CFD tool, and it was found that unlike the smooth wall case, slosh damping with baffle is almost independent of the working fluid and it is reasonable to apply water test data to the full scale LOX tank when the damping from baffle is dominant. On the other hand, for the smooth wall, the damping value must be scaled according to the Reynolds number. Comparison of experimental data, CFD, with the classical and modified Miles equations for upper dome was made, and the limitations of these semi-empirical equations were identified.

Yang, H. Q.↗

Empirical Evaluation of GPS Clock Accuracy for Isochronous Droop-based Inverters

Grid-forming inverters rely upon an internal reference to regulate the grid 's voltage and frequency. Without a reliable reference, issues such as loss of synchronism, frequency instability, and unequal power sharing can occur. One proposed internal phase reference for grid-forming inverters is generated by the controller's clock, but this is affected by nonidealities in the crystal oscillator leading to time drift in the phase reference. To improve the controller's clock accuracy, the Global Positioning System (GPS) can be used to calibrate the internal clock. Previously, GPS-based clock calibration was explored theoretically and various methods were developed for mitigating frequency/phase drift at the power system level. This paper experimentally validates the direct performance of GPS time-calibration at the controller level and justifies its use in maintaining an accurate phase reference for droop-based grid-forming inverters.

clock drift↗

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

A comparison of measured and theoretical predictions for STS ascent and entry sonic booms

Sonic boom measurements have been obtained during the flights of STS-1 through 5. During STS-1, 2, and 4, entry sonic boom measurements were obtained and ascent measurements were made on STS-5. The objectives of this measurement program were (1) to define the sonic boom characteristics of the Space Transportation System (STS), (2) provide a realistic assessment of the validity of xisting theoretical prediction techniques, and (3) establish a level of confidence for predicting future STS configuration sonic boom environments. Detail evaluation and reporting of the results of this program are in progress. This paper will address only the significant results, mainly those data obtained during the entry of STS-1 at Edwards Air Force Base (EAFB), and the ascent of STS-5 from Kennedy Space Center (KSC). The theoretical prediction technique employed in this analysis is the so called Thomas Program. This prediction technique is a semi-empirical method that required definition of the near field signatures, detailed trajectory characteristics, and the prevailing meteorological characteristics as an input. This analytical procedure then extrapolates the near field signatures from the flight altitude to an altitude consistent with each measurement location.

Garcia, F., Jr.↗