Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling analysis”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗
NASA NTRSDate not supplied

ANALYSIS OF THE MSL/MEDLI ENTRY DATA WITH COUPLED CFD AND MATERIAL RESPONSE.

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heat-shield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material. PICA is a lightweight carbon fiber/polymeric resin material that offers out-standing performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on Open-FOAM (PATO) using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA, is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Mars Science Laboratory↗
NASA NTRS1978

Monte Carlo solution of Boltzmann equation for a simple model of highly nonequilibrium diatomic gases - Translational rotational energy relaxation

Theoretical results pertaining to internally excited translational-rotational energy relaxation in a spatially uniform diatomic gas far removed from solid boundaries are obtained by solving the Boltzmann equation by means of the Monte Carlo direct simulation method. The analysis is based on calculations involving three different types of initial conditions: equilibrium, nonequilibrium-equipartition (i.e., equipartition is satisfied, but the distributions are perturbed), and nonequilibrium-nonequipartition (i.e., both equipartition and the distributions are perturbed). Results of monatomic-gas simulations are also included to facilitate comparisons with the coupled translational-rotational relaxation simulations, and some simulations for a normal shock-wave structure are briefly examined. The results show that: (1) single-step transitions are the significant mechanisms of intermodal energy transfer; (2) translational-rotational transitions are coupled most efficiently for low-lying states of rotationally excited molecules and least efficiently for highly rotationally excited molecules; and (3) relaxation occurs via a successive set of distributions that are not Maxwell-Boltzmann (nonlocal Maxwellian).

Yoshikawa, K. K.↗
NASA NTRSDate not supplied

Analysis of MSL/MEDLI Entry Data with Coupled CFD and Material Response

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Thermal Protection Systems↗
DOE OSTI2021

DPC Disposal Thermal Scoping Analysis

This is a progress report on thermal modeling for dual-purpose canister (DPCs) direct disposal that covers several available calculation methods and addresses creep and temperature-dependent properties in a salt repository. Three modeling approaches are demonstrated: A semi-analytical calculation method that uses linear solutions with superposition and imaging, to represent a central waste package in a larger array; A finite difference model of coupled thermal creep, implemented in FLAC2D; and An integrated finite difference thermal-hydrologic modeling approach for repositories in different generic host media, implemented in PFLOTRAN. These approaches are at different levels of maturity, and future work is expected to add refinements and establish the best applications for each.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1980

Topographic coupling of surface and internal Kelvin waves

An analysis is presented for computing the diffraction of barotropic Kelvin waves by a localized topographical irregularity on flat-bottom ocean with an arbitrary vertical stratification. It was shown that all baroclinic Kelvin waves will be generated downstream of the bump, with the first baroclinic mode having the largest amplitude. The Poincare waves predominate in the lowest modes, and are more directionally anisotropic. It was concluded that baroclinic Poincare waves radiating offshore from the bump topography could contribute to the internal wave field in the open ocean and provide an alternative mechanism to dissipate the barotropic tides.

Chao, S.-Y.↗
NASA NTRS2005

Testing the Archivas Cluster (Arc) for Ozone Monitoring Instrument (OMI) Scientific Data Storage

The Ozone Monitoring Instrument (OMI) launched on NASA's Aura Spacecraft, the third of the major platforms of the EOS program on July 15,2004. In addition to the long term archive and distribution of the data from OM1 through the Goddard Earth Science Distributed Active Archive Center (GESDAAC), we are evaluating other archive mechanisms that can archive the data in a more immediately available method where it can be used for futher data production and analysis. In 2004, Archivas, Inc. was selected by NASA s Small Business Innovative Research (SBIR) program for the development of their Archivas Cluster (ArC) product. Arc is an online disk based system utilizing self-management and automation on a Linux cluster. Its goal is to produce a low cost solution coupled with the ease of management. The OM1 project is an application partner of the SBIR program, and has deployed a small cluster (5TB) based on the beta Archwas software. We performed extensive testing of the unit using production OM1 data since launch. In 2005, Archivas, Inc. was funded in SBIR Phase II for further development, which will include testing scalability with the deployment of a larger (35TB) cluster at Goddard. We plan to include Arc in the OM1 Team Leader Computing Facility (TLCF) hosting OM1 data for direct access and analysis by the OMI Science Team. This presentation will include a brief technical description of the Archivas Cluster, a summary of the SBIR Phase I beta testing results, and an overview of the OMI ground data processing architecture including its interaction with the Phase II Archivas Cluster and hosting of OMI data for the scientists.

Tilmes, Curt↗
DOE OSTI2023

Evaluation Of Glass Density to Support the Estimation of Fissile Mass Loadings in Sludge Batch 10 Glasses

Per a directive from the Department of Energy Savannah River Operations Office (DOE-SR) in 2008, the fissile mass loading concentration must remain below 897 g/m 3 in each high-level waste (HLW) glass canister produced by the Defense Waste Processing Facility (DWPF). To support Sludge Batch 5 (SB5) processing, the Savannah River National Laboratory (SRNL) developed a technical basis that facilitates the evaluation of fissile mass loading of the glass product. The calculation is based on the iron (Fe) concentration in the glass as determined by measurements from the Slurry Mix Evaporator acceptability analysis as well as the glass density. In April 2022, a subsequent DOE-SR directive increased the fissile mass loading limit to 2500 g/m 3 beginning with Sludge Batch 11. Thus, the 897 g/m 3 limit still applies to Sludge Batch 10 (SB10) processing. For SB5 through initial Sludge Batch 9 (SB9) processing prior to coupled operation with the Salt Waste Processing Facility (SWPF), SRNL provided DWPF a bounding glass density value that was based on a statistical evaluation of density measurements. To eliminate the need for experimental work, a composition-based density model for HLW glasses was developed at SRNL in 2019. The objective of this report is to present the bounding glass density determined with the composition-based density model for SB10 sludge-only (SO) and coupled processing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1972

Structural evaluation of Marman V-band coupling and flange with conoseal gasket

Results are described of a development test program directed at evaluating the structural capabilities of the Marman V-band coupling and flange with conoseal gasket. The intended end use was for the 75K NERVA flight engine propellant lines. Of major importance in the structural evaluation was the ability to predict stresses throughout the assembly for a variety of loading conditions. Computer finite element analysis was used to predict these stresses but, for the subject configuration, large uncertainties were introduced in modeling the complex geometry and boundary conditions. The purpose of the structural tests was to obtain actual stresses and deflections for correlation with, and updating of the finite element model. Results of the incomplete test program are inconclusive with respect to determining suitability for use on the NERVA engine.

Oates, J. H.↗
NASA NTRS1985

Measurement and Modeling of the Optical Scattering Properties of Crop Canopies

The specular reflection process is shown to be a key aspect of radiation transfer by plant canopies. Polarization measurements are demonstrated as the tool for determining the specular and diffuse portions of the canopy radiance. The magnitude of the specular fraction of the reflectance is significant compared to the magnitude of the diffuse fraction. Therefore, it is necessary to consider specularly reflected light in developing and evaluating light-canopy interaction models for wheat canopies. Models which assume leaves are diffuse reflectors correctly predict only the diffuse fraction of the canopy reflectance factor. The specular reflectance model, when coupled with a diffuse leaf model, would predict both the specular and diffuse portions of the reflectance factor. The specular model predicts and the data analysis confirms that the single variable, angle of incidence of specularly reflected sunlight on the leaf, explains much of variation in the polarization data as a function of view-illumination directions.

Vanderbilt, V. C.↗
NASA NTRS1998

Noncontact Determination of Antisymmetric Plate Wave Velocity in Ceramic Matrix Composites

High-temperature materials are of increasing importance in the development of more efficient engines and components for the aeronautics industry. In particular, ceramic matrix composite (CMC) and metal matrix composite (MMC) structures are under active development for these applications. The acousto-ultrasonic (AU) method has been shown to be useful for assessing mechanical properties in composite structures. In particular, plate wave analysis can characterize composites in terms of their stiffness moduli. It is desirable to monitor changes in mechanical properties that occur during thermomechanical testing and to monitor the health of components whose geometry or position make them hard to reach with conventional ultrasonic probes. In such applications, it would be useful to apply AU without coupling directly to the test surface. For a number of years, lasers have been under investigation as remote ultrasonic input sources and ultrasound detectors. The use of an ultrasonic transducer coupled through an air gap has also been under study. So far at the NASA Lewis Research Center, we have been more successful in using lasers as ultrasonic sources than as output devices. On the other hand, we have been more successful in using an air-coupled piezoelectric transducer as an output device than as an input device. For this reason, we studied the laser in/air-coupled-transducer out combination-using a pulsed NdYAG laser as the ultrasonic source and an air-coupled-transducer as the detector. The present work is focused on one of the AU parameters of interest, the ultrasonic velocity of the antisymmetric plate-wave mode. This easily identified antisymmetric pulse can be used to determine shear and flexure modulus. It was chosen for this initial work because the pulse arrival times are likely to be the most precise. The following schematic illustrates our experimental arrangement for using laser in/air-transducer out on SiC/SiC composite tensile specimens. The NdYAG pulse was directed downward by a 90 infrared prism to the top of the specimen, but at the edge of one end. An energy sensor measured a single pulse at 13 millijoules (mJ) before it passed through the prism, which attenuated 15 percent of its energy. It also provided an output trigger for the waveform time-delay synthesizer.

Kautz, Harold E.↗
NASA NTRS2002

Obtaining Reliable Predictions of Terrestrial Energy Coupling From Real-Time Solar Wind Measurements

Measurements of the interplanetary magnetic field (IMF) from the ACE (Advanced Composition Explorer), Wind, IMP-8 (Interplanetary Monitoring Platform), and Geotail spacecraft have revealed that the IMF variations are contained in phase planes that are tilted with respect to the propagation direction, resulting in continuously variable changes in propagation times between spacecraft, and therefore, to the Earth. Techniques for using 'minimum variance analysis' have been developed in order to be able to measure the phase front tilt angles, and better predict the actual propagation times from the L1 orbit to the Earth, using only the real-time IMF measurements from one spacecraft. The use of empirical models with the IMF measurements at L1 from ACE (or future satellites) for predicting 'space weather' effects has also been demonstrated.

Weimer, Daniel R.↗
DOE OSTI2025

Utilization of the LS-APGD microplasma/orbitrap-FTMS booster system for detection and isotopic analysis of neodymium nanoparticles

Detection and isotopic analysis of particle populations has seen rapid growth across several application areas, including environmental analysis, nuclear forensics, and food safety. The ability to characterize the particles' unique elemental and isotopic fingerprints could provide information related to formation, processing history, and transport. Regarding nuclear forensics, isotopic analysis of particles derived from diverse materials is often used as a tool to trace the origin and processing history. Mass spectrometric-based techniques currently used for particle population analysis often suffer from limited mass resolution, particularly when dealing with real-world samples that are affected by isobaric and polyatomic interferences from the matrix. To address these analytical challenges, we propose a novel method utilizing the liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma ionization source coupled to an ultrahigh resolution Orbitrap mass spectrometer, further enhanced with the FTMS X2T Booster data acquisition and processing unit. The FTMS Booster enables acquisition of extended transient times of up to 3 s, significantly improving mass resolution, thereby reducing or even eliminating the need for prior separation of isobaric or polyatomic interferences. Additionally, the detection of low-abundance isotopes was improved by increasing the signal-to-noise (S/N) ratio. As proof of concept, this study demonstrates the feasibility of the LS-APGD/Orbitrap-FTMS X2T Booster platform for direct analysis using a suspension of well-characterized ∼120 nm neodymium particles. The quality of the isotope ratios values obtained from a few hundred particles were in good agreement with those obtained from homogeneous ionic solutions. These results highlight the potential of the LS-APGD/Orbitrap platform for rapid, accurate, and interference-resilient isotope ratio analysis of particle populations without the need for dissolution and subsequent chemical separations, offering significant advantages for nuclear forensics, safeguards, and environmental applications. The effort here also points to further paths forward, hopefully towards single particle (SP) analysis using microplasma ionization and the ultrahigh resolving power of the Orbitrap mass analyzer.

FTMS X2T booster↗
DOE OSTI2024

Anode Upcycling via Tailored Solvent Treatment

To achieve a truly closed-loop direct recycling process for lithium-ion batteries, all component materials must be recovered. To date, direct recycling method development has primarily focused on the high-value transition-metal cathode materials, while the inherently lower-value graphite has been challenging to recover in a cost-effective manner. However, end-of-life graphite contains a unique engineered value due to the presence of the solid electrolyte interphase (SEI). Growth of the SEI during the cell's active lifetime stabilizes the electronically reactive graphite surface through an irreversible consumption of Li, and thus necessitates both excess lithiation of the cathode and a costly and time-intensive formation procedure during manufacturing. An optimized pre-formed SEI that capitalizes on existing SEI components from end-of-life batteries has the potential to significantly reduce cathode lithiation requirements and eliminate the critical bottleneck of formation cycling during cell remanufacturing. Further, retaining Li at the anode obviates the need for a separate Li leaching and recovery step, improving the overall efficiency of the direct recycling line. In this work, we present a novel approach to "upcycling" spent graphite through use of tailored chemical treatment to remove adverse (i.e., highly resistive and/or poorly passivating) SEI species while retaining beneficially passivating components. We have explored a rational set of solvents spanning a range of polarity, proticity, and molecular size to evaluate structure-property-performance relationships between applied solvent(s), removed and remaining SEI species, and electrochemical response of the resulting graphite product. Further, we have developed and optimized a robust and holistic analysis procedure that couples symmetric-cell electrochemical testing, multi-modal materials characterization, and advanced electrochemical modeling. These analysis results inform a set of correlative metrics for graphite performance relative to both solvent properties and upcycled SEI composition. We demonstrate effective tunability in the residual SEI composition by varying solvent identity and concentration, and report on several promising solvent systems that achieve comparable or performance to pristine graphite.

anode recycling↗
DOE OSTI2024

Fourier Analysis and Loss Modeling for Inductive Wireless Electric Vehicle Charging with Reduced Stray Field

With the growth of electric vehicle (EV) popularity, different charging options to increase user convenience and reduce charging times are being considered and researched. Among these, inductive wireless power transfer (WPT) systems for EVs are being designed to meet specifications such as stray field, power level, efficiency, misalignment tolerance, and ground clearance, which are all heavily influenced by the coil geometry. The proposed Fourier Analysis Method (FAM) is an analytical method to directly design coil geometries to meet stray field and power level requirements through an optimization of Fourier basis function coefficients. The outputs of the optimization are complex, planar coil geometries that meet the power level and stray field constraints with minimized loss factors. Contours of these potentials determine the coil conductor paths and loss models predict the system efficiency and performance over misalignment. The Fourier representation of the geometry is used to conveniently calculate the coupling over misalignment, external proximity effect loss, and ferrite loss. A 6.6 kW prototype WPT system with low stray field and high efficiency is built from the optimization results to validate the models and showcase the usefulness of the FAM design approach.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1996

MSC/NASTRAN DMAP Alter Used for Closed-Form Static Analysis With Inertia Relief and Displacement-Dependent Loads

Solving for the displacements of free-free coupled systems acted upon by static loads is a common task in the aerospace industry. Often, these problems are solved by static analysis with inertia relief. This technique allows for a free-free static analysis by balancing the applied loads with the inertia loads generated by the applied loads. For some engineering applications, the displacements of the free-free coupled system induce additional static loads. Hence, the applied loads are equal to the original loads plus the displacement-dependent loads. A launch vehicle being acted upon by an aerodynamic loading can have such applied loads. The final displacements of such systems are commonly determined with iterative solution techniques. Unfortunately, these techniques can be time consuming and labor intensive. Because the coupled system equations for free-free systems with displacement-dependent loads can be written in closed form, it is advantageous to solve for the displacements in this manner. Implementing closed-form equations in static analysis with inertia relief is analogous to implementing transfer functions in dynamic analysis. An MSC/NASTRAN (MacNeal-Schwendler Corporation/NASA Structural Analysis) DMAP (Direct Matrix Abstraction Program) Alter was used to include displacement-dependent loads in static analysis with inertia relief. It efficiently solved a common aerospace problem that typically has been solved with an iterative technique.

Source record↗
NASA NTRS2018

WCRP Task Team for the Intercomparison of ReAnalyses (TIRA): Motivation and Progress

Reanalyses have proven to be an important resource for weather and climate related research, as well as societal applications at large. New atmospheric reanalyses in various forms are produced every few years. In addition, land and ocean communities are producing disciplinary uncoupled reanalyses, and regional reanalyses are emerging in addition to global reanalyses. Systems of data information are being developed to facilitate the intercomparison of reanalyses among themselves and with observations and climate change experiments (e.g. CREATE-IP and WRIT). Current research and development in reanalysis is directed at (1) extending the length of reanalyzed period and (2) use of coupled Earth system models for climate reanalysis. While the World Climate Research Programme (WCRP) involvement in the reanalyses communities through its Data Advisory Council (WDAC) has been substantial, for example in organizing international conferences on reanalyses. However, these efforts have previously not included an analysis of the differences among reanalyses and their inherent uncertainties which are important questions for both users and developers of reanalyses. Therefore, the WCRP Task Team for the Intercomparison of ReAnalyses (TIRA) task team has been formed as collaborative effort to address reanalyses role in WCRP. This constitutes a logical progression that fills the needs of the community and contributes to the WCRP mission. The primary charge to TIRA is to develop a reanalysis intercomparison project plan that will attain the following objectives.1)To foster understanding and estimation of uncertainties in reanalysis data by intercomparison and other means)To communicate new developments and best practices among the reanalyses producing centers3)To enhance the understanding of data and assimilation issues and their impact on uncertainties, leading to improved reanalyses for climate assessment4)To communicate the strengths and weaknesses of reanalyses, their fitness for purpose, and best practices in the use of reanalysis datasets by the scientific communityThis contribution outlines the need for a task team on reanalyses, their intercomparison, the objectives of the team and progress thus far. Results for some pilot intercomparisons will also be discussed.

Bosilovich, Michael↗