Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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,459 records · Page 81

Bayesian calibration of thermodynamic parameters for geochemical speciation modeling of cementitious materials

Chemical equilibrium modeling of cementitious materials requires aqueous-solid equilibrium constants of the controlling mineral phases (K{sub sp}) and the available concentrations of primary components. Inherent randomness of the input and model parameters, experimental measurement error, the assumptions and approximations required for numerical simulation, and inadequate knowledge of the chemical process contribute to uncertainty in model prediction. A numerical simulation framework is developed in this paper to assess uncertainty in K{sub sp} values used in geochemical speciation models. A Bayesian statistical method is used in combination with an efficient, adaptive Metropolis sampling technique to develop probability density functions for K{sub sp} values. One set of leaching experimental observations is used for calibration and another set is used for comparison to evaluate the applicability of the approach. The estimated probability distributions of K{sub sp} values can be used in Monte Carlo simulation to assess uncertainty in the behavior of aqueous-solid partitioning of constituents in cement-based materials.

36 MATERIALS SCIENCE↗

X-ray imaging crystal spectroscopy for use in plasma transport research

This research describes advancements in the spectral analysis and error propagation techniques associated with x-ray imaging crystal spectroscopy (XICS) that have enabled this diagnostic to be used to accurately constrain particle, momentum, and heat transport studies in a tokamak for the first time. Doppler tomography techniques have been extended to include propagation of statistical uncertainty due to photon noise, the effect of non-uniform instrumental broadening as well as flux surface variations in impurity density. These methods have been deployed as a suite of modeling and analysis tools, written in interactive data language (IDL) and designed for general use on tokamaks. Its application to the Alcator C-Mod XICS is discussed, along with novel spectral and spatial calibration techniques. Example ion temperature and radial electric field profiles from recent I-mode plasmas are shown, and the impact of poloidally asymmetric impurity density and natural line broadening is discussed in the context of the planned ITER x-ray crystal spectrometer.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

THE INTERNAL ROTATION PROFILE OF THE B-TYPE STAR KIC 10526294 FROM FREQUENCY INVERSION OF ITS DIPOLE GRAVITY MODES

The internal angular momentum distribution of a star is the key to determining its evolution. Fortunately, stellar internal rotation can be probed through studies of rotationally split nonradial oscillation modes. In particular, the detection of nonradial gravity modes (g modes) in massive young stars has recently become feasible thanks to the Kepler space mission. Our goal is to derive the internal rotation profile of the Kepler B8V star KIC 10526294 through asteroseismology. We interpret the observed rotational splittings of its dipole g modes using four different approaches based on the best seismic models of the star and their rotational kernels. We show that these kernels can resolve differential rotation within the radiative envelope if a smooth rotational profile is assumed and if the observational errors are small. Based on Kepler data, we find that the rotation rate near the core-envelope boundary is well constrained to 163 ± 89 nHz. The seismic data are consistent with rigid rotation but a profile with counter-rotation within the envelope has a statistical advantage over constant rotation. Our study should be repeated for other massive stars with a variety of stellar parameters in order to determine the physical conditions that control the internal rotation profile of young massive stars, with the aim of improving the input physics of their models.

79 ASTRONOMY AND ASTROPHYSICS↗

Sub-divertor fuel isotopic content detection limit for JET and its impact on ICRF core heating and DTE2 operation

The ability to detect and control fuel isotopic content down to a 1% concentration level is greatly important for the upcoming JET DTE2 campaign, as well as its associated TT and DD phases. A reduction of H minority concentration even from 2% down to 1% is shown in this paper to have significant impact on the effectiveness of ion cyclotron range of frequencies core heating, while the ability to maintain T or D concentration at or below 1% is critical to limiting fusion neutron generation in the DD and TT phases, correspondingly. The sub-divertor measurement of (global) isotopic concentration, based on Penning-activated optical spectroscopy, can deliver minimally this 1% detection for DTE2 as long as light collection from the Penning emission can be optimized and gradual window transmission deterioration can be minimized. Here, this is simulated with a statistical analysis developed to understand the uncertainty sources in the JET DTE1 data, as well as to guide the optimization of an upgraded, fuel-isotopic content (and helium-ash concentration) gas analysis system for the JET divertor in preparation for DTE2. While this random error can be reduced to allow measurement substantially below 1% concentration, analysis also shows a systematic error of up to 1% understood to be due to plasma–surface interactions in the Penning excitation, suggesting that 1% may still be the low-end limit for the sub-divertor measurement, unless a Penning-source conditioning approach is also developed.

fusion fuel isotopic content↗

Bayesian Uncertainty Quantification in Predictions of Flows in Highly Heterogeneous Media and Its Applications to the CO 2 Sequestration

In this proposal, we have worked on Bayesian uncertainty quantification for predictions of fows in highly heterogeneous media. The research in this proposal is broad and includes: prior modeling for heterogeneous permeability fields; effective parametrization of heterogeneous spatial priors; efficient ensemble- level solution techniques; efficient multiscale approximation techniques; study of the regularity of complex posterior distribution and the error estimates due to parameter reduction; efficient sampling techniques; applications to multi-phase ow and transport. We list our publications below and describe some of our main research activities. Our multi-disciplinary team includes experts from the areas of multiscale modeling, multilevel solvers, Bayesian statistics, spatial permeability modeling, and the application domain.

97 MATHEMATICS AND COMPUTING↗

Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations

A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.

16 TIDAL AND WAVE POWER↗

IMPROVED PERFORMANCES IN SUBSONIC FLOWS OF AN SPH SCHEME WITH GRADIENTS ESTIMATED USING AN INTEGRAL APPROACH

In this paper, we present results from a series of hydrodynamical tests aimed at validating the performance of a smoothed particle hydrodynamics (SPH) formulation in which gradients are derived from an integral approach. We specifically investigate the code behavior with subsonic flows, where it is well known that zeroth-order inconsistencies present in standard SPH make it particularly problematic to correctly model the fluid dynamics. In particular, we consider the Gresho–Chan vortex problem, the growth of Kelvin–Helmholtz instabilities, the statistics of driven subsonic turbulence and the cold Keplerian disk problem. We compare simulation results for the different tests with those obtained, for the same initial conditions, using standard SPH. We also compare the results with the corresponding ones obtained previously with other numerical methods, such as codes based on a moving-mesh scheme or Godunov-type Lagrangian meshless methods. We quantify code performances by introducing error norms and spectral properties of the particle distribution, in a way similar to what was done in other works. We find that the new SPH formulation exhibits strongly reduced gradient errors and outperforms standard SPH in all of the tests considered. In fact, in terms of accuracy, we find good agreement between the simulation results of the new scheme and those produced using other recently proposed numerical schemes. These findings suggest that the proposed method can be successfully applied for many astrophysical problems in which the presence of subsonic flows previously limited the use of SPH, with the new scheme now being competitive in these regimes with other numerical methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Flexible and scalable methods for quantifying stochastic variability in the era of massive time-domain astronomical data sets

We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and variability-based classification. We show that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs. We present the likelihood function for light curves sampled from CARMA processes, placing them on a statistically rigorous foundation, and we present a Bayesian method to infer the probability distribution of the PSD given the measured light curve. Because calculation of the likelihood function scales linearly with the number of data points, CARMA modeling scales to current and future massive time-domain data sets. We conclude by applying our CARMA modeling approach to light curves for an X-ray binary, two active galactic nuclei, a long-period variable star, and an RR Lyrae star in order to illustrate their use, applicability, and interpretation.

79 ASTRONOMY AND ASTROPHYSICS↗

Uncertainty bounds for multivariate machine learning predictions on high-strain brittle fracture

Simulation of the crack network evolution on high strain rate impact experiments performed in brittle materials is very compute-intensive. The cost increases even more if multiple simulations are needed to account for the randomness in crack length, location, and orientation, which is inherently found in real-world materials. Constructing a machine learning emulator can make the process faster by orders of magnitude. There has been little work, however, on assessing the error associated with their predictions. Estimating these errors is imperative for meaningful overall uncertainty quantification. In this work, we extend the heteroscedastic uncertainty estimates to bound a multiple output machine learning emulator. Overall, we find that the response prediction is accurate within its predicted errors, but with a somewhat conservative estimate of uncertainty.

36 MATERIALS SCIENCE↗

The Residual Setup Errors of Different IGRT Alignment Procedures for Head and Neck IMRT and the Resulting Dosimetric Impact

Purpose: To assess residual setup errors during head and neck radiation therapy and the resulting consequences for the delivered dose for various patient alignment procedures. Methods and Materials: Megavoltage cone beam computed tomography (MVCBCT) scans from 11 head and neck patients who underwent intensity modulated radiation therapy were used to assess setup errors. Each MVCBCT scan was registered to its reference planning kVCT, with seven different alignment procedures: automatic alignment and manual registration to 6 separate bony landmarks (sphenoid, left/right maxillary sinuses, mandible, cervical 1 [C1]-C2, and C7-thoracic 1 [T1] vertebrae). Shifts in the different alignments were compared with each other to determine whether there were any statistically significant differences. Then, the dose distribution was recalculated on 3 MVCBCT images per patient for every alignment procedure. The resulting dose-volume histograms for targets and organs at risk (OARs) were compared to those from the planning kVCTs. Results: The registration procedures produced statistically significant global differences in patient alignment and actual dose distribution, calling for a need for standardization of patient positioning. Vertically, the automatic, sphenoid, and maxillary sinuses alignments mainly generated posterior shifts and resulted in mean increases in maximal dose to OARs of >3% of the planned dose. The suggested choice of C1-C2 as a reference landmark appears valid, combining both OAR sparing and target coverage. Assuming this choice, relevant margins to apply around volumes of interest at the time of planning to take into account for the relative mobility of other regions are discussed. Conclusions: Use of different alignment procedures for treating head and neck patients produced variations in patient setup and dose distribution. With concern for standardizing practice, C1-C2 reference alignment with relevant margins around planning volumes seems to be a valid option.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

ON THE RELIABILITY OF POLARIZATION ESTIMATION USING ROTATION MEASURE SYNTHESIS

We benchmark the reliability of the rotation measure (RM) synthesis algorithm using the 1005 Centaurus A field sources of Feain et al. The RM synthesis solutions are compared with estimates of the polarization parameters using traditional methods. This analysis provides verification of the reliability of RM synthesis estimates. We show that estimates of the polarization parameters can be made at lower signal-to-noise ratio (S/N) if the range of RMs is bounded, but reliable estimates of individual sources with unusual RMs require unconstrained solutions and higher S/N. We derive from first principles the statistical properties of the polarization amplitude associated with RM synthesis in the presence of noise. The amplitude distribution depends explicitly on the amplitude of the underlying (intrinsic) polarization signal. Hence, it is necessary to model the underlying polarization signal distribution in order to estimate the reliability and errors in polarization parameter estimates. We introduce a Bayesian method to derive the distribution of intrinsic amplitudes based on the distribution of measured amplitudes. The theoretically derived distribution is compared with the empirical data to provide quantitative estimates of the probability that an RM synthesis solution is correct as a function of S/N. We provide quantitative estimates of the probability that any given RM synthesis solution is correct as a function of measured polarized amplitude and the intrinsic polarization amplitude compared to the noise.

79 ASTRONOMY AND ASTROPHYSICS↗

The Dark Energy Survey Year 3 and eBOSS: constraining galaxy intrinsic alignments across luminosity and colour space

ABSTRACT We present direct constraints on galaxy intrinsic alignments (IAs) using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift z ∼ 0.2–1.0, luminous red galaxies from eBOSS at z ∼ 0.8, and also an SDSS-III BOSS CMASS sample at z ∼ 0.5. We measure two-point IA correlations, which we fit using a model that includes lensing, magnification, and photometric redshift error. Fitting on scales 6 Mpc h−1 < rp < 70 Mpc h−1, we make a detection of IAs in each sample, at 5σ–22σ (assuming a simple one-parameter model for IAs). Using these red samples, we measure the IA–luminosity relation. Our results are statistically consistent with previous results, but offer a significant improvement in constraining power, particularly at low luminosity. With this improved precision, we see detectable dependence on colour between broadly defined red samples. It is likely that a more sophisticated approach than a binary red/blue split, which jointly considers colour and luminosity dependence in the IA signal, will be needed in future. We also compare the various signal components at the best-fitting point in parameter space for each sample, and find that magnification and lensing contribute $\sim 2\!-\!18~{{\ \rm per\ cent}}$ of the total signal. As precision continues to improve, it will certainly be necessary to account for these effects in future direct IA measurements. Finally, we make equivalent measurements on a sample of emission-line galaxies from eBOSS at z ∼ 0.8. We constrain the non-linear alignment amplitude to be $A_1=0.07^{+0.32}_{-0.42}$ (|A1| < 0.78 at 95 per cent CL).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental determination of precision, resolution, accuracy and trueness of time-of-flight neutron diffraction strain measurements

We discuss a simple statistical analysis which yields the precision, resolution, accuracy and trueness of diffraction-based lattice strain measurements. The procedure consists of measuring the thermal expansion induced in each component of an ideal, non-reacting, two-component, crystalline powder sample in situ. One component, with a high coefficient of thermal expansion (CTE), serves as an internal thermometer. The quantities of interest are obtained by determining the smallest, statistically significant, thermal lattice strain which can be detected through diffraction analysis in the second, low CTE, component in response to controlled temperature changes. This procedure also provides a robust check of the alignment of the diffraction system and is able to reveal the presence of systematic errors. The application of this technique to a time-of-flight engineering diffractometer/strain scanner is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

BLENDING STUDY FOR SRR SALT DISPOSITION INTEGRATION: TANK 50H SCALE-MODELING AND COMPUTER-MODELING FOR BLENDING PUMP DESIGN, PHASE 2

The Salt Disposition Integration (SDI) portfolio of projects provides the infrastructure within existing Liquid Waste facilities to support the startup and long term operation of the Salt Waste Processing Facility (SWPF). Within SDI, the Blend and Feed Project will equip existing waste tanks in the Tank Farms to serve as Blend Tanks where 300,000-800,000 gallons of salt solution will be blended in 1.3 million gallon tanks and qualified for use as feedstock for SWPF. Blending requires the miscible salt solutions from potentially multiple source tanks per batch to be well mixed without disturbing settled sludge solids that may be present in a Blend Tank. Disturbing solids may be problematic both from a feed quality perspective as well as from a process safety perspective where hydrogen release from the sludge is a potential flammability concern. To develop the necessary technical basis for the design and operation of blending equipment, Savannah River National Laboratory (SRNL) completed scaled blending and transfer pump tests and computational fluid dynamics (CFD) modeling. A 94 inch diameter pilot-scale blending tank, including tank internals such as the blending pump, transfer pump, removable cooling coils, and center column, were used in this research. The test tank represents a 1/10.85 scaled version of an 85 foot diameter, Type IIIA, nuclear waste tank that may be typical of Blend Tanks used in SDI. Specifically, Tank 50 was selected as the tank to be modeled per the SRR, Project Engineering Manager. SRNL blending tests investigated various fixed position, non-rotating, dual nozzle pump designs, including a blending pump model provided by the blend pump vendor, Curtiss Wright (CW). Primary research goals were to assess blending times and to evaluate incipient sludge disturbance for waste tanks. Incipient sludge disturbance was defined by SRR and SRNL as minor blending of settled sludge from the tank bottom into suspension due to blending pump operation, where the sludge level was shown to remain constant. To experimentally model the sludge layer, a very thin, pourable, sludge simulant was conservatively used for all testing. To experimentally model the liquid, supernate layer above the sludge in waste tanks, two salt solution simulants were used, which provided a bounding range of supernate properties. One solution was water (H{sub 2}O + NaOH), and the other was an inhibited, more viscous salt solution. The research performed and data obtained significantly advances the understanding of fluid mechanics, mixing theory and CFD modeling for nuclear waste tanks by benchmarking CFD results to actual experimental data. This research significantly bridges the gap between previous CFD models and actual field experiences in real waste tanks. A finding of the 2009, DOE, Slurry Retrieval, Pipeline Transport and Plugging, and Mixing Workshop was that CFD models were inadequate to assess blending processes in nuclear waste tanks. One recommendation from that Workshop was that a validation, or bench marking program be performed for CFD modeling versus experiment. This research provided experimental data to validate and correct CFD models as they apply to mixing and blending in nuclear waste tanks. Extensive SDI research was a significant step toward bench marking and applying CFD modeling. This research showed that CFD models not only agreed with experiment, but demonstrated that the large variance in actual experimental data accounts for misunderstood discrepancies between CFD models and experiments. Having documented this finding, SRNL was able to provide correction factors to be used with CFD models to statistically bound full scale CFD results. Through the use of pilot scale tests performed for both types of pumps and available engineering literature, SRNL demonstrated how to effectively apply CFD results to salt batch mixing in full scale waste tanks. In other words, CFD models were in error prior to development of experimental correction factors determined during this research, which provided a technique to use CFD models for salt batch mixing and transfer pump operations. This major scientific advance in mixing technology resulted in multi-million dollar cost savings to SRR. New techniques were developed for both experiment and analysis to complete this research. Supporting this success, research findings are summarized in the Conclusions section of this report, and technical recommendations for design and operation are included in this section of the report.

08 HYDROGEN↗

Evaluating Application Resilience with XRay

The rising count and shrinking feature size of transistors within modern computers is making them increasingly vulnerable to various types of soft faults. This problem is especially acute in high-performance computing (HPC) systems used for scientific computing, because these systems include many thousands of compute cores and nodes, all of which may be utilized in a single large-scale run. The increasing vulnerability of HPC applications to errors induced by soft faults is motivating extensive work on techniques to make these applications more resilient to such faults, ranging from generic techniques such as replication or checkpoint/restart to algorithm-specific error detection and tolerance techniques. Effective use of such techniques requires a detailed understanding of how a given application is affected by soft faults to ensure that (i) efforts to improve application resilience are spent in the code regions most vulnerable to faults and (ii) the appropriate resilience technique is applied to each code region. This paper presents XRay, a tool to view the application vulnerability to soft errors, and illustrates how XRay can be used in the context of a representative application. In addition to providing actionable insights into application behavior XRay automatically selects the number of fault injection experiments required to provide an informative view of application behavior, ensuring that the information is statistically well-grounded without performing unnecessary experiments.

97 MATHEMATICS AND COMPUTING↗

Computational and analytical comparison of flux discretizations for the semiconductor device equations beyond Boltzmann statistics

We compare three thermodynamically consistent numerical fluxes known in the literature, appearing in a Voronoï finite volume discretization of the van Roosbroeck system with general charge carrier statistics. Our discussion includes an extension of the Scharfetter–Gummel scheme to non-Boltzmann (e.g. Fermi–Dirac) statistics. It is based on the analytical solution of a two-point boundary value problem obtained by projecting the continuous differential equation onto the interval between neighboring collocation points. Hence, it serves as a reference flux. The exact solution of the boundary value problem can be approximated by computationally cheaper fluxes which modify certain physical quantities. One alternative scheme averages the nonlinear diffusion (caused by the non-Boltzmann nature of the problem), another one modifies the effective density of states. To study the differences between these three schemes, we analyze the Taylor expansions, derive an error estimate, visualize the flux error and show how the schemes perform for a carefully designed p-i-n benchmark simulation. We present strong evidence that the flux discretization based on averaging the nonlinear diffusion has an edge over the scheme based on modifying the effective density of states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interaction of loading pattern and nuclear data uncertainties in reactor core calculations

Along with best-estimate calculations for design and safety analysis, understanding uncertainties is important to determine appropriate design margins. In this framework, nuclear data uncertainties and their propagation to full core calculations are a critical issue. To deal with this task, different error propagation techniques, deterministic and stochastic are currently developed to evaluate the uncertainties in the output quantities. Among these is the sampling based uncertainty and sensitivity software XSUSA which is able to quantify the influence of nuclear data covariance on reactor core calculations. In the present work, this software is used to investigate systematically the uncertainties in the power distributions of two PWR core loadings specified in the OECD UAM-Benchmark suite. With help of a statistical sensitivity analysis, the main contributors to the uncertainty are determined. Using this information a method is studied with which loading patterns of reactor cores can be optimized with regard to minimizing power distribution uncertainties. It is shown that this technique is able to halve the calculation uncertainties of a MOX/UOX core configuration. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Porting Adaptive Ensemble Molecular Dynamics Workflows to the Summit Supercomputer

Molecular dynamics (MD) simulations must take very small (femtosecond) integration steps in simulation-time to avoid numerical errors. Efficient use of parallel programming models and accelerators in state-of-the art MD programs now is pushing Moore’s limit for time-per-MD step. As a result, directly simulating timescales beyond milliseconds will not be attainable directly, even at exascale. However, concepts from statistical physics can be used to combine many parallel simulations to provide information about longer timescales and to adequately sample the simulation space, while preserving details about the dynamics of the system. Implementing such an approach requires a workflow program that allows adaptable steering of task assignments based on extensive statistical analysis of intermediate results. Here we report the implementation of such an adaptable workflow program to drive simulations on the Summit IBM Power System AC922, a pre-exascale supercomputer at the Oak Ridge Leadership Computing Facility (OLCF). We compare to experiences on Titan, Summit’s predecessor, report the performance of the workflow and its components, and describe the porting process. We find that using a workflow program managed by a Mongo database can provide the fault tolerance, scalable performance, task dispatch rate, and reconfigurability required for robust and portable implementation of ensemble simulations such as are used in enhanced-sampling molecular dynamics. This type of workflow generator can also be used to provide adaptive steering of ensemble simulations for other applications in addition to MD.

Ossyra, John R.↗