Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion processes”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

On the Chemical Mixing Induced by Internal Gravity Waves

Detailed modeling of stellar evolution requires a better understanding of the (magneto)hydrodynamic processes that mix chemical elements and transport angular momentum. Understanding these processes is crucial if we are to accurately interpret observations of chemical abundance anomalies, surface rotation measurements, and asteroseismic data. Here, we use two-dimensional hydrodynamic simulations of the generation and propagation of internal gravity waves in an intermediate-mass star to measure the chemical mixing induced by these waves. We show that such mixing can generally be treated as a diffusive process. We then show that the local diffusion coefficient does not depend on the local fluid velocity, but rather on the wave amplitude. We then use these findings to provide a simple parameterization for this diffusion, which can be incorporated into stellar evolution codes and tested against observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Multiscale models and stochastic simulation methods for computing rare but key binding events in cell biology

The main difficulty in simulating diffusion processes at a molecular level in cell microdomains is due to the multiple scales involving nano- to micrometers. Few to many particles have to be simulated and simultaneously tracked while there are exploring a large portion of the space for binding small targets, such as buffers or active sites. Bridging the small and large spatial scales is achieved by rare events representing Brownian particles finding small targets and characterized by long-time distribution. These rare events are the bottleneck of numerical simulations. A naive stochastic simulation requires running many Brownian particles together, which is computationally greedy and inefficient. Solving the associated partial differential equations is also difficult due to the time dependent boundary conditions, narrow passages and mixed boundary conditions at small windows. We present here two reduced modeling approaches for a fast computation of diffusing fluxes in microdomains. The first approach is based on a Markov mass-action law equations coupled to a Markov chain. The second is a Gillespie's method based on the narrow escape theory for coarse-graining the geometry of the domain into Poissonian rates. The main application concerns diffusion in cellular biology, where we compute as an example the distribution of arrival times of calcium ions to small hidden targets to trigger vesicular release.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Co-Based superalloy morphology evolution: A phase field study based on experimental thermodynamic and kinetic data

Cobalt-based superalloys with two phase $\gamma$/$\gamma '$ microstructures offer great promise as candidates for next-generation high-temperature alloys for applications, such as turbine blades. It is essential to understand the thermodynamic and kinetic factors that influence the microstructural evolution of these alloys in order to optimize the alloy compositions and processing steps with a goal to improve their coarsening, creep and rafting behavior. We are using a continuum phase field approach to study the diffusion process and to predict the equilibrium shapes of Co-Al-W $\gamma '$ precipitates. In order to obtain quantitatively predictive capabilities, we extract chemical free energies for the $\gamma$/$\gamma '$ phases based on CALculation of PHAse Diagrams (CALPHAD) thermodynamic data and diffusion mobilities for Co alloys based on CALPHAD kinetic data. We also use experimental or first-principles data for other quantities, such as misfit strain and interface information, for the parameterization of our model. A particular focus of our study is to understand how different energy balances, misfit strain and kinetics affect the coarsening and rafting behavior of $\gamma '$ precipitates, and the sensitivity of the final precipitate shape to materials parameters. Here, we find that the equilibrium shape of the precipitate results from a delicate competition between chemical, interfacial, and elastic energies, and it is very sensitive to changes in model parameters. Here, we examine how modeling input parameters affect the equilibrium shape of precipitates and relate these parameters to experimentally available values.

36 MATERIALS SCIENCE↗
DOE OSTI2017

High efficiency screen-printed n-type silicon solar cell using co-diffusion of APCVD boron emitter and POCl 3 back surface field

We present the fabrication and analysis of Passivated Emitter and Rear Totally Diffused (PERT) solar cells on n-type silicon using a co-diffusion process. In a single high temperature step, a BSG/SiO x stack deposited by APCVD and a POCl 3 back surface field diffuse into the wafer to form the boron doped emitter and phosphorus doped back surface field. The SiO x layer on top of BSG acts as a masking layer to prevent cross-doping of phosphorus as well as a blocking layer for boron out-diffusion. This resulted in an initial sheet resistance of 76 Ω with good uniformity and a final p + emitter sheet resistance of 97 Ω after boron rich layer removal. Additionally, bulk lifetime was investigated before and after the high temperature step that resulted in an increase from 1.2 ms to 1.5 ms due to a POCl 3 gettering effect. Finally, a peak cell efficiency of 20.3% was achieved and each recombination component in terms of saturation current density was calculated and analyzed to understand the cell for further efficiency enhancement.

14 SOLAR ENERGY↗
DOE OSTI2018

Ultrafast molecular transport on carbon surfaces: The diffusion of ammonia on graphite

We present a combined experimental and theoretical study of the self-diffusion of ammonia on exfoliated graphite. Using neutron time-of-flight spectroscopy we are able to resolve the ultrafast diffusion process of adsorbed ammonia, NH 3 , on graphite. Together with van der Waals corrected density functional theory calculations we show that the diffusion of NH 3 follows a hopping motion on a weakly corrugated potential energy surface with an activation energy of about 4 meV which is particularly low for this type of diffusive motion. The hopping motion includes further a significant number of long jumps and the diffusion constant of ammonia adsorbed on graphite is determined with D = 3.9•10 -8 m 2 /s at 94 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Bayesian uncertainty quantification of tristructural isotropic particle fuel silver release: Decomposing model inadequacy plus experimental noise and parametric uncertainties

Tristructural isotropic (TRISO) particle fuel is one of the most promising fuel concepts enabling high temperature and high burnup reactor operation. One dominant source of radioactivity released from the TRISO particles is silver (Ag), which is subject to a high release fraction and long decay life compared to other fission products. Previous modeling efforts using the fuel performance code BISON indicated nonnegligible uncertainties in modeling the diffusion process of fission products in TRISO compared to the Advanced Gas Reactor experiments. The overall uncertainties observed when modeling the fission product diffusion can result from uncertainties in model parameters, noisy experimental measurements, and deficiencies in the developed models. The three types of underlying uncertainties have not yet been properly quantified in open literature. Here, this paper presents the Bayesian uncertainty quantification (UQ) using massively parallelizable Markov chain Monte Carlo samplers. The uncertainties due to model parameters, model inadequacy, and experimental measurement noise are quantified, with the σ term used to represent the sum of the model inadequacy and measurement noise uncertainties. It is worth noting that this is the first time the σ term is inferred for nuclear fuel experiments, as compared to using prescribed values for uncertainty quantification in previous work. The parallelizable Markov chain Monte Carlo samplers efficiently infer the model parameters and the σ term, giving insight into physical parameters like diffusion coefficients and the combined model discrepancy and measurement noise. A subsequent forward uncertainty quantification (UQ) is also performed based on the calibration results to generate more accurate predictions of the Ag release. The model inadequacy plus experimental noise is the most dominant source of uncertainty compared to the parametric uncertainty. All the UQ analyses presented in this work are based on the second series of the irradiation experiments in the Advanced Gas Reactor program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Using Active Learning to Rapidly Develop Machine Learned Diffusion Coefficients of CO 2 Conversion Reagents in Metal–Organic Frameworks

Here, we used a combined molecular dynamics/active learning (AL) approach to create machine learning models that can predict the diffusion coefficient of epichlorohydrin and chloropropene carbonate, the reactant and product of a common CO 2 cycloaddition reaction, in metal–organic frameworks (MOFs). Nanoporous MOFs are effective catalysts for the cycloaddition of CO 2 to epoxides. The diffusion rates within nanoporous catalysts can control the rate of reaction as the reactants and products must diffuse to the active sites within the MOF and then out of the nanoporous material for reusability. However, the diffusion process is routinely ignored when searching for new materials in catalytic applications. Here we verified improvement during the AL process by consistently tracking metrics on the same groups of MOFs to ensure consistency. Metal identity was found to have little impact on diffusion rates, while structural features like pore limiting diameter act as a threshold where a minimum value is needed for high diffusion rates. We identified the MOFs with the highest epichlorohydrin and chloropropene carbonate diffusion coefficients which can be used for further studies of reaction energetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Extended families of critical and stationary droplets for nonequilibrium phase transitions in spatially discrete bistable systems

Bistable nonequilibrium systems are realized in catalytic reaction-diffusion processes, biological transport and regulation, spatial epidemics, etc. Behavior in spatially continuous formulations, described at the mean-field level by reaction-diffusion type equations (RDEs), often mimics that of classic equilibrium van der Waals type systems. When accounting for noise, similarities include a discontinuous phase transition at some value, p eq , of a control parameter, p , with metastability and hysteresis around p eq . For each p , there is a unique critical droplet of the more stable phase embedded in the less stable or metastable phase which is stationary (neither shrinking nor growing), and with size diverging as p → p eq . Spatially discrete analogs of these mean-field formulations, described by lattice differential equations (LDEs), are more appropriate for some applications, but have received less attention. It is recognized that LDEs can exhibit richer behavior than RDEs, specifically propagation failure for planar interphases separating distinct phases. Herein, we show that this feature, together with an orientation dependence of planar interface propagation also deriving from spatial discreteness, results in the occurrence of entire families of stationary droplets. The extent of these families increases approaching the transition and can be infinite if propagation failure is realized. In addition, there can exist a regime of generic two-phase coexistence where arbitrarily large droplets of either phase always shrink. Such rich behavior is qualitatively distinct from that for classic nucleation in equilibrium and spatially continuous nonequilibrium systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

MOSCATO Development and Integration in Fiscal Year 2023

In FY23, significant developments were made to the MOSCATO code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: Code Restructuring: The MOSCATO code underwent a restructuring process, enabling both Nek5000 and NekRS to share the same interface when using MOSCATO. This integration facilitates smoother interactions between the different modules and streamlines the overall simulation process. Automatic Domain Extension: To allow simulations to propagate over longer time frames, an automatic domain extension capability was added to the metal diffusion module. This enhancement enables more accurate and reliable simulations of diffusion processes, especially in scenarios with extended time requirements. Expanded Corrosion Reactions: The updated MOSCATO package now includes additional reactions that contribute to corrosion. This addition enhances the code’s ability to model corrosion effects more comprehensively, providing a more realistic representation of the systems being studied. Example Case Setup: To assist users and showcase MOSCATO’s capabilities, we developed example cases for various applications. Notably, we created example cases for thermal cells and corrosion loops. These pre-configured cases allow users to get started quickly and effectively with their simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

FVCOM-plume – A Three-dimensional Lagrangian Outfall Plume Dilution and Transport Model for Dynamic Tidal Environments: Model Development

Although numerous steady-state effluent plume dilution models are in use, their application in tidal environments remains a challenge. Three-dimensional dynamic circulation models are also inadequate, often due to the lack of required resolution and simplifying assumptions. To overcome these limitations, we present FVCOM-plume—an outfall plume dilution and transport model that operates within the Finite Volume Coastal Ocean Model (FVCOM) framework. It provides simultaneous inclusion of near-field dilution and far-field plume transport processes. The near-field is based on UM3 model using Lagrangian Control Volume approach to compute buoyant plume trajectory and dilution from multiport diffusers. The far-field uses neutrally buoyant particles with point masses and the random walk method to solve unsteady advection-diffusion processes. A density kernel approach is used to compute concentrations at point locations and analyze far-field plume characteristics. The results demonstrate the ability of FVCOM-plume to simultaneously capture near-field and far-field effluent plume dynamics in tidal environments.

Plume Dilution Models↗
DOE OSTI2024

Site Disorder Drives Cyanide Dynamics and Fast Ion Transport in Li 6 PS 5 CN

Halide argyrodite solid-state electrolytes of the general formula Li 6 PS 5 X exhibit complex static and dynamic disorder that plays a crucial role in ion transport processes. Here, we unravel the rich interplay between site disorder and dynamics in the plastic crystal argyrodite Li 6 PS 5 CN and the impact on ion diffusion processes through a suite of experimental and computational methodologies, including temperature-dependent synchrotron powder X-ray diffraction, AC electrochemical impedance spectroscopy, 7 Li solid-state NMR, and machine learning-assisted molecular dynamics simulations. Sulfide and (pseudo)halide site disorder between the two anion sublattices unilaterally improves long-range lithium diffusion irrespective of the (pseudo)halide identity, which demonstrates the importance of site disorder in dictating bulk ionic conductivity in the argyrodite family. Furthermore, we find that anion site disorder modulates the presence and time scales of cyanide rotational dynamics. Ordered configurations of anions enable fast, quasi-free rotations of cyanides that occur on time scales of 10 11 Hz at T = 300 K. In contrast, we find that cyanide dynamics are slow or frozen in Li 6 PS 5 CN when site disorder between the cyanide and sulfide sublattices is present at T = 300 K. We rationalize the observed differences in cyanide dynamics in the context of elastic dipole interactions between neighboring cyanide anions and local strain induced by the configurations of site disorder that may impact the energetic landscape for cyanide rotational dynamics. Through this study, we find that anion disorder plays a decisive role in dictating the extent and time scales of both lithium ion and cyanide dynamics in Li 6 PS 5 CN.

36 MATERIALS SCIENCE↗
DOE OSTI2017

ZEEMAN DOPPLER MAPS: ALWAYS UNIQUE, NEVER SPURIOUS?

Numerical models of atomic diffusion in magnetic atmospheres of ApBp stars predict abundance structures that differ from the empirical maps derived with (Zeeman) Doppler mapping (ZDM). An in-depth analysis of this apparent disagreement investigates the detectability by means of ZDM of a variety of abundance structures, including (warped) rings predicted by theory, but also complex spot-like structures. Even when spectra of high signal-to-noise ratio are available, it can prove difficult or altogether impossible to correctly recover shapes, positions, and abundances of a mere handful of spots, notwithstanding the use of all four Stokes parameters and an exactly known field geometry; the recovery of (warped) rings can be equally challenging. Inversions of complex abundance maps that are based on just one or two spectral lines usually permit multiple solutions. It turns out that it can by no means be guaranteed that any of the regularization functions in general use for ZDM (maximum entropy or Tikhonov) will lead to a true abundance map instead of some spurious one. Attention is drawn to the need for a study that would elucidate the relation between the stratified, field-dependent abundance structures predicted by diffusion theory on the one hand, and empirical maps obtained by means of “canonical” ZDM, i.e., with mean atmospheres and unstratified abundances, on the other hand. Finally, we point out difficulties arising from the three-dimensional nature of the atomic diffusion process in magnetic ApBp star atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Development of APCVD BSG and POCl 3 Codiffusion Process for Double-Side TOPCon Solar Cell Precursor Fabrication

This paper presents a commercially viable process for fabricating a high-quality double-side tunnel oxide passivating contact (DS-TOPCon) cell precursor using APCVD-deposited boron silicate glass and ex-situ POCl 3 diffusion in a single high-temperature step, eliminating the need for additional masking and diffusion processes. A two-tier temperature profile was developed, involving a pre-annealing at above 900°C in nitrogen (N2) ambient followed by POCl 3 diffusion at 840°C. We investigated the effect of varying pre-annealing temperatures, ranging from 875°C to 950°C, on the passivation quality and metal-Si contact properties of both n-TOPCon and p-TOPCon layers. The resultant DS-TOPCon cell precursor after silicon nitride (SiNX) passivation exhibited an excellent iV OC of close to 730 mV. In addition, a rapid asymmetric poly-Si thinning technique, developed in this work, enabled adjustment of the front n + poly-Si thickness while maintaining the rear p + poly-Si thickness. Two types of DS-TOPCon cell architectures can be fabricated: i) full-area thin (≈40nm) n-TOPCon layer on the front and ii) selective-area thick (≈200nm) n-TOPCon fingers underneath the metal grid. Device simulations suggest that full-area DS-TOPCon cell with 40 nm n + poly-Si and selective-area DS-TOPCon cell with 200 nm n + poly fingers on the front, fabricated from our current DS-TOPCon cell precursor, can achieve cell efficiencies of 22.1 and 23.5%, respectively. Detailed power loss analysis and device simulation reveal that further improvements in material and device parameters have potential to push the cell efficiencies of DS-TOPCon cell structure beyond 25%, making it a promising alternative to fabricate a high-efficiency next-generation solar cells at low cost.

14 SOLAR ENERGY↗
DOE OSTI2017

Theoretical insight into Cobalt subnano-clusters adsorption on α-Al{sub 2}O{sub 3} (0001)

The investigation on the structural stability, nucleation, growth and interaction of cobalt cluster Con(n=2–7) on the α-Al{sub 2}O{sub 3}(0001) surface by using density functional theory methods has been reported. Energetically, the most favorable adsorption sites were identified and the strongest adsorption energy cluster is the tetrahedral Co{sub 4} cluster. On the other hand, the nucleation of Con(n=2–7) clusters on the surface is exothermic and thermodynamically favorable. Moreover, even-odd alternation was found with respect to clusters nucleation as a function of the number of cobalt atoms (for n=1–7). Meanwhile, the Co{sub n} clusters can be adsorbed on the surface stably owing to the charge transfer from Co atoms to Al and O atoms of the Al{sub 2}O{sub 3} substrate. In addition, we establish the crucial importance of monomer, dimer and trimer diffusion on the surface. The diffusion of the monomer cobalt from Al{sup (3)} to O{sup (5)} or O{sup (5)} to Al{sup (4)} site is quite easy on the Al{sub 2}O{sub 3}(0001) surface, whereas the diffusion of the Co{sub 2} dimer is thermodynamically unfavorable by compared with that of the Co adatom and Co{sub 3} trimer. - Graphical abstract: Diffusion process of Co adatom on the α-Al{sub 2}O{sub 3} (0001) surface, Al{sup (3)} site→O{sup (5)} site→Al{sup (4)} site. Potential energy surface for diffusion of a single Co atom from Al{sup (3)} to O{sup (5)} site, and from O{sup (5)} to Al{sup (4)} site on the surface. The activation energy of the two migration processes from Al{sup (3)} to O{sup (5)} and O{sup (5)} to Al{sup (4)} are 0.06 and 0.09 eV, respectively. This implies the monomer is quite mobile on the surface under typical growth conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Quantum Effects on H 2 Diffusion in Zeolite RHO: Inverse Kinetic Isotope Effect for Sieving

We use canonical variational theory (CVT) with small-curvature tunneling (SCT) contributions to investigate quantum effects on the H 2 diffusion process in the pure-silica zeolite RHO. At low temperature we find an inverse kinetic isotopic sieving effect in that the heavier isotopic species diffuses faster than the lighter one. Three quantum effects contribute to this kinetic isotope effect (KIE). The first one is quantum mechanical tunneling; this—on its own—would lead to a normal kinetic isotopic sieving effect, in which lighter diprotium diffuses faster than dideuterium. The second factor, which we find to dominate in the present case, is zero-point energy (ZPE). Deuterium has a lower ZPE, which leads to a smaller effective barrier for tunneling because the transition state has a larger ZPE than the precursor stable state; this results in an inverse KIE. The third factor, the thermal vibrational energy (computed from the quantized vibrational partition function), also favors a normal KIE, but it is outweighed by the ZPE effect. The vibrations of the zeolite host framework are found to play an important role at low temperatures, and our calculations consider up to 7296 degrees of freedom at the equilibrium structure and the saddle point and up to 221 degrees of freedom along the reaction path. In conclusion, the importance of quantum considerations on the dynamics is elucidated by comparison to a purely classical treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Temperature Dependence of Self-diffusion in Cr2O3 from First Principles

Understanding and predicting the dominant diffusion processes in Cr2O3 is essen- tial to its optimization for anti-corrosion coatings, spintronics, and other applications. Despite significant theoretical effort in modeling defect mediated diffusion in Cr2O3 the correlation with experimentally measured diffusivities remains poor partly due to the insufficient accuracy of the theoretical approaches. Here an attempt to resolve these discrepancies is made through high accuracy density functional theory simulations cou- pled with grand canonical formalism of defect thermochemistry. In this approach, point defect formation energies were computed using hybrid exchange correlation functional. This level of theory proved to be essential for achieving the agreement with experi- mental self-diffusion coefficients. The analysis of the resulting self-diffusion coefficients indicate that chromium has higher mobility at low temperatures and high oxygen partial pressures, in particular at standard temperature and pressure conditions. At high vacuum, high temperature conditions, oxygen diffusion becomes dominant. At Cr/Cr2O3 interfaces O vacancies were found to be more mobile than Cr vacancies at all temperatures. Cr diffuses preferentially along the c-axis at low temperatures but switches to basal plane at higher temperatures. O diffusion is primarily bound to basal plane at all temperatures.

Medasani, Bharat K.↗