Showing 1–27 of 27 results for author: Roy, T

Shape matters: Understanding the effect of electrode geometry on cell resistance and chemo-mechanical stress

Authors: Tiras Y. Lin, Hanyu Li, Nicholas W. Brady, Nicholas R. Cross, Victoria M. Ehlinger, Thomas Roy, Daniel Tortorelli, Christine Orme, Marcus A. Worsley, Giovanna Bucci

Abstract: Rechargeable batteries that incorporate shaped three-dimensional electrodes have been shown to have increased power and energy densities for a given footprint area when compared to a conventional geometry, i.e., a planar cathode and anode that sandwich an electrolyte. Electrodes can be shaped to enable a higher loading of active material, while keeping the ion transport distance small, however, th… ▽ More Rechargeable batteries that incorporate shaped three-dimensional electrodes have been shown to have increased power and energy densities for a given footprint area when compared to a conventional geometry, i.e., a planar cathode and anode that sandwich an electrolyte. Electrodes can be shaped to enable a higher loading of active material, while keeping the ion transport distance small, however, the relationship between electrical and mechanical performance remains poorly understood. A variety of electrode shapes have been explored, where the electrodes are individually shaped or intertwined with one another. Advances in manufacturing and shape and topology optimization have made such designs a reality. In this paper, we explore sinusoidal half cells and interdigitated full cells. First, we use a simple electrostatics model to understand the cell resistance as a function of shape. We focus on low-temperature conditions, where the electrolyte conductivity decreases and the governing dimensionless parameters change. Next, we use a chemo-mechanics model to examine the stress concentrations that arise due to intercalation-driven volume expansion. We show that shaped electrodes provide a significant reduction in resistance, however, they result in unfavorable stress concentrations. Overall, we find that the fully interdigitated electrodes may provide the best balance with respect to this trade-off. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Topology Optimization for the Full-Cell Design of Porous Electrodes in Electrochemical Energy Storage Devices

Authors: Hanyu Li, Giovanna Bucci, Nicholas W. Brady, Nicholas R. Cross, Victoria M. Ehlinger, Tiras Y. Lin, Miguel Salazar de Troya, Daniel Tortorelli, Marcus A. Worsley, Thomas Roy

Abstract: In this paper, we introduce a density-based topology optimization framework to design porous electrodes for maximum energy storage. We simulate the full cell with a model that incorporates electronic potential, ionic potential, and electrolyte concentration. The system consists of three materials, namely pure liquid electrolyte and the porous solids of the anode and cathode, for which we determine… ▽ More In this paper, we introduce a density-based topology optimization framework to design porous electrodes for maximum energy storage. We simulate the full cell with a model that incorporates electronic potential, ionic potential, and electrolyte concentration. The system consists of three materials, namely pure liquid electrolyte and the porous solids of the anode and cathode, for which we determine the optimal placement. We use separate electronic potentials to model each electrode, which allow interdigitated designs. As the result, a penalization is required to ensure that the anode and cathode do not touch, i.e. causing a short circuit. We compare multiple 2D designs generated for different fixed conditions, e.g. material properties. A 3D design with complex channel and interlocking structure is also created. All optimized designs are far superior to the traditional monolithic electrode design with respect to energy storage metrics. We observe up to 750% increase in energy storage for cases with slow effective ionic diffusion within the porous electrode. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Train Small, Model Big: Scalable Physics Simulators via Reduced Order Modeling and Domain Decomposition

Authors: Seung Whan Chung, Youngsoo Choi, Pratanu Roy, Thomas Moore, Thomas Roy, Tiras Y. Lin, Du Y. Nguyen, Christopher Hahn, Eric B. Duoss, Sarah E. Baker

Abstract: Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the E-pilot, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scale… ▽ More Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the E-pilot, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. This method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about $15 - 40$ times faster with only $\sim$ $1\%$ relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation. △ Less

Submitted 5 December, 2023; originally announced January 2024.

Comments: 40 pages, 12 figures. Submitted to Computer Methods in Applied Mechanics and Engineering

Report number: LLNL-JRNL-857774 MSC Class: 65F55; 65N55 (primary) 76D07 (secondary)

Spatial wavefront shaping with a multipolar-resonant metasurface for structured illumination microscopy

Authors: Tamal Roy, Peter T. Brown, Douglas P. Shepherd, Lisa V. Poulikakos

Abstract: Structured illumination microscopy (SIM) achieves superresolution in fluorescence imaging through patterned illumination and computational image reconstruction, yet current methods require bulky, costly modulation optics and high-precision optical alignment. This work demonstrates how nano-optical metasurfaces, rationally designed to tailor the optical wavefront at sub-wavelength dimensions, hold… ▽ More Structured illumination microscopy (SIM) achieves superresolution in fluorescence imaging through patterned illumination and computational image reconstruction, yet current methods require bulky, costly modulation optics and high-precision optical alignment. This work demonstrates how nano-optical metasurfaces, rationally designed to tailor the optical wavefront at sub-wavelength dimensions, hold great potential as ultrathin, single-surface, all-optical wavefront modulators for SIM. We computationally demonstrate this principle with a multipolar-resonant metasurface composed of silicon nanostructures which generate versatile optical wavefronts in the far field upon variation of the polarization or angle of incident light. Algorithmic optimization is performed to identify the seven most suitable illumination patterns for SIM generated by the metasurface based on three key criteria. We find that multipolar-resonant metasurface SIM (mrm-SIM) achieves resolution comparable to conventional methods by applying the seven optimal metasurface-generated wavefronts to simulated fluorescent objects and reconstructing the objects using proximal gradient descent. The work presented here paves the way for a metasurface-enabled experimental simplification of structured illumination microscopy. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: TR and PTB contributed equally to this work

Journal ref: Opt. Mater. Express 14(5), 1239-1252 (2024)

Interplay of degeneracy and non-degeneracy in fluctuations propagation in coherent feed-forward loop motif

Authors: Tuhin Subhra Roy, Mintu Nandi, Pinaki Chaudhury, Sudip Chattopadhyay, Suman K Banik

Abstract: We present a stochastic framework to decipher fluctuations propagation in classes of coherent feed-forward loops. The systematic contribution of the direct (one-step) and indirect (two-step) pathways is considered to quantify fluctuations of the output node. We also consider both additive and multiplicative integration mechanisms of the two parallel pathways (one-step and two-step). Analytical exp… ▽ More We present a stochastic framework to decipher fluctuations propagation in classes of coherent feed-forward loops. The systematic contribution of the direct (one-step) and indirect (two-step) pathways is considered to quantify fluctuations of the output node. We also consider both additive and multiplicative integration mechanisms of the two parallel pathways (one-step and two-step). Analytical expression of the output node's coefficient of variation shows contributions of intrinsic, one-step, two-step, and cross-interaction in closed form. We observe a diverse range of degeneracy and non-degeneracy in each of the decomposed fluctuations term and their contribution to the overall output fluctuations of each coherent feed-forward loop motif. Analysis of output fluctuations reveals a maximal level of fluctuations of the coherent feed-forward loop motif of type 1. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 20 pages, 4 figures

Journal ref: J. Stat. Mech. 2023 (2023) 093502

Exploration of superconducting multi-mode cavity architectures for quantum computing

Authors: Alessandro Reineri, Silvia Zorzetti, Tanay Roy, Xinyuan You

Abstract: Superconducting radio-frequency (SRF) cavities coupled to transmon circuits have proven to be a promising platform for building high-coherence quantum information processors. An essential aspect of this realization involves designing high quality factor three-dimensional superconducting cavities to extend the lifetime of quantum systems. To increase the computational capability of this architectur… ▽ More Superconducting radio-frequency (SRF) cavities coupled to transmon circuits have proven to be a promising platform for building high-coherence quantum information processors. An essential aspect of this realization involves designing high quality factor three-dimensional superconducting cavities to extend the lifetime of quantum systems. To increase the computational capability of this architecture, we are exploring a multimode approach. This paper presents the design optimization process of a multi-cell SRF cavity to perform quantum computation based on an existing design developed in the scope of particle accelerator technology. We perform parametric electromagnetic simulations to evaluate and optimize the design. In particular, we focus on the analysis of the interaction between a nonlinear superconducting circuit known as the transmon and the cavity. This parametric design optimization is structured to serve as a blueprint for future studies on similar systems. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 2023 Intl Conf on Quantum Computing and Engineering (QCE23)

Report number: FERMILAB-CONF-23-405-SQMS

An Autoencoder-based Snow Drought Index

Authors: Sinan Rasiya Koya, Kanak Kanti Kar, Shivendra Srivastava, Tsegaye Tadesse, Mark Svoboda, Tirthankar Roy

Abstract: In several regions across the globe, snow has a significant impact on hydrology. The amounts of water that infiltrate the ground and flow as runoff are driven by the melting of snow. Therefore, it is crucial to study the magnitude and effect of snowmelt. Snow droughts, resulting from reduced snow storage, can drastically impact the water supplies in basins where snow predominates, such as in the w… ▽ More In several regions across the globe, snow has a significant impact on hydrology. The amounts of water that infiltrate the ground and flow as runoff are driven by the melting of snow. Therefore, it is crucial to study the magnitude and effect of snowmelt. Snow droughts, resulting from reduced snow storage, can drastically impact the water supplies in basins where snow predominates, such as in the western United States. Hence, it is important to detect the time and severity of snow droughts efficiently. We propose Snow Drought Response Index or SnoDRI, a novel indicator that could be used to identify and quantify snow drought occurrences. Our index is calculated using cutting-edge ML algorithms from various snow-related variables. The self-supervised learning of an autoencoder is combined with mutual information in the model. In this study, we use random forests for feature extraction for SnoDRI and assess the importance of each variable. We use reanalysis data (NLDAS-2) from 1981 to 2021 for the Pacific United States to study the efficacy of the new snow drought index. We evaluate the index by confirming the coincidence of its interpretation and the actual snow drought incidents. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Temporal Fusion Transformers for Streamflow Prediction: Value of Combining Attention with Recurrence

Authors: Sinan Rasiya Koya, Tirthankar Roy

Abstract: Over the past few decades, the hydrology community has witnessed notable advancements in streamflow prediction, particularly with the introduction of cutting-edge machine-learning algorithms. Recurrent neural networks, especially Long Short-Term Memory (LSTM) networks, have become popular due to their capacity to create precise forecasts and realistically mimic the system dynamics. Attention-based… ▽ More Over the past few decades, the hydrology community has witnessed notable advancements in streamflow prediction, particularly with the introduction of cutting-edge machine-learning algorithms. Recurrent neural networks, especially Long Short-Term Memory (LSTM) networks, have become popular due to their capacity to create precise forecasts and realistically mimic the system dynamics. Attention-based models, such as Transformers, can learn from the entire data sequence concurrently, a feature that LSTM does not have. This work tests the hypothesis that combining recurrence with attention can improve streamflow prediction. We set up the Temporal Fusion Transformer (TFT) architecture, a model that combines both of these aspects and has never been applied in hydrology before. We compare the performance of LSTM, Transformers, and TFT over 2,610 globally distributed catchments from the recently available Caravan dataset. Our results demonstrate that TFT indeed exceeds the performance benchmark set by the LSTM and Transformers for streamflow prediction. Additionally, being an explainable AI method, TFT helps in gaining insights into the streamflow generation processes. △ Less

Submitted 20 May, 2023; originally announced May 2023.

Differentiable modeling to unify machine learning and physical models and advance Geosciences

Authors: Chaopeng Shen, Alison P. Appling, Pierre Gentine, Toshiyuki Bandai, Hoshin Gupta, Alexandre Tartakovsky, Marco Baity-Jesi, Fabrizio Fenicia, Daniel Kifer, Li Li, Xiaofeng Liu, Wei Ren, Yi Zheng, Ciaran J. Harman, Martyn Clark, Matthew Farthing, Dapeng Feng, Praveen Kumar, Doaa Aboelyazeed, Farshid Rahmani, Hylke E. Beck, Tadd Bindas, Dipankar Dwivedi, Kuai Fang, Marvin Höge , et al. (5 additional authors not shown)

Abstract: Process-Based Modeling (PBM) and Machine Learning (ML) are often perceived as distinct paradigms in the geosciences. Here we present differentiable geoscientific modeling as a powerful pathway toward dissolving the perceived barrier between them and ushering in a paradigm shift. For decades, PBM offered benefits in interpretability and physical consistency but struggled to efficiently leverage lar… ▽ More Process-Based Modeling (PBM) and Machine Learning (ML) are often perceived as distinct paradigms in the geosciences. Here we present differentiable geoscientific modeling as a powerful pathway toward dissolving the perceived barrier between them and ushering in a paradigm shift. For decades, PBM offered benefits in interpretability and physical consistency but struggled to efficiently leverage large datasets. ML methods, especially deep networks, presented strong predictive skills yet lacked the ability to answer specific scientific questions. While various methods have been proposed for ML-physics integration, an important underlying theme -- differentiable modeling -- is not sufficiently recognized. Here we outline the concepts, applicability, and significance of differentiable geoscientific modeling (DG). "Differentiable" refers to accurately and efficiently calculating gradients with respect to model variables, critically enabling the learning of high-dimensional unknown relationships. DG refers to a range of methods connecting varying amounts of prior knowledge to neural networks and training them together, capturing a different scope than physics-guided machine learning and emphasizing first principles. Preliminary evidence suggests DG offers better interpretability and causality than ML, improved generalizability and extrapolation capability, and strong potential for knowledge discovery, while approaching the performance of purely data-driven ML. DG models require less training data while scaling favorably in performance and efficiency with increasing amounts of data. With DG, geoscientists may be better able to frame and investigate questions, test hypotheses, and discover unrecognized linkages. △ Less

Submitted 26 December, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Journal ref: Nat Rev Earth Environ 4, 552-567 (2023)

A scalable DG solver for the electroneutral Nernst-Planck equations

Authors: Thomas Roy, Julian Andrej, Victor A. Beck

Abstract: The robust, scalable simulation of flowing electrochemical systems is increasingly important due to the synergy between intermittent renewable energy and electrochemical technologies such as energy storage and chemical manufacturing. The high Péclet regime of many such applications prevents the use of off-the-shelf discretization methods. In this work, we present a high-order Discontinuous Galerki… ▽ More The robust, scalable simulation of flowing electrochemical systems is increasingly important due to the synergy between intermittent renewable energy and electrochemical technologies such as energy storage and chemical manufacturing. The high Péclet regime of many such applications prevents the use of off-the-shelf discretization methods. In this work, we present a high-order Discontinuous Galerkin scheme for the electroneutral Nernst-Planck equations. The chosen charge conservation formulation allows for the specific treatment of the different physics: upwinding for advection and migration, and interior penalty for diffusion of ionic species as well the electric potential. Similarly, the formulation enables different treatments in the preconditioner: AMG for the potential blocks and ILU-based methods for the advection-dominated concentration blocks. We evaluate the convergence rate of the discretization scheme through numerical tests. Strong scaling results for two preconditioning approaches are shown for a large 3D flow-plate reactor example. △ Less

Submitted 20 December, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: 26 pages, 7 figures

MSC Class: 65N30; 65F08; 78A57 ACM Class: G.1.8; G.1.3

Topology optimization for the design of porous electrodes

Authors: Thomas Roy, Miguel A. Salazar de Troya, Marcus A. Worsley, Victor A. Beck

Abstract: Porous electrodes are an integral part of many electrochemical devices since they have high porosity to maximize electrochemical transport and high surface area to maximize activity. Traditional porous electrode materials are typically homogeneous, stochastic collections of small scale particles and offer few opportunities to engineer higher performance. Fortunately, recent breakthroughs in advanc… ▽ More Porous electrodes are an integral part of many electrochemical devices since they have high porosity to maximize electrochemical transport and high surface area to maximize activity. Traditional porous electrode materials are typically homogeneous, stochastic collections of small scale particles and offer few opportunities to engineer higher performance. Fortunately, recent breakthroughs in advanced and additive manufacturing are yielding new methods to structure and pattern porous electrodes across length scales. These architected electrodes are emerging as a promising new technology to continue to drive improvement; however, it is still unclear which structures to employ and few tools are available to guide their design. In this work we address this gap by applying topology optimization to the design of porous electrodes. We demonstrate our framework on two applications: a porous electrode driving a steady Faradaic reaction and a transiently operated electrode in a supercapacitor. We present computationally designed electrodes that minimize energy losses in a half-cell. For low conductivity materials, the optimization algorithm creates electrode designs with a hierarchy of length scales. Further, the designed electrodes are found to outperform undesigned, homogeneous electrodes. Finally, we present three-dimensional porous electrode designs. We thus establish a topology optimization framework for designing porous electrodes. △ Less

Submitted 19 May, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 21 pages, 14 figures

MSC Class: 35Q93; 49M41; 74P15; 35Q92 ACM Class: G.1.6; G.1.8; J.6.1

An alternate approach to simulate the dynamics of perturbed liquid drops

Authors: Tanu Singla, Tanushree Roy, P. Parmananda, M. Rivera

Abstract: Liquid drops when subjected to external periodic perturbations can execute polygonal oscillations. In this work, a simple model is presented that demonstrates these oscillations and their characteristic properties. The model consists of a spring-mass network such that the masses are analogous to liquid molecules and the springs are to intermolecular forces. Neo-Hookean springs are considered to re… ▽ More Liquid drops when subjected to external periodic perturbations can execute polygonal oscillations. In this work, a simple model is presented that demonstrates these oscillations and their characteristic properties. The model consists of a spring-mass network such that the masses are analogous to liquid molecules and the springs are to intermolecular forces. Neo-Hookean springs are considered to represent these intermolecular forces. The restoring force of a neo-Hookean spring depends nonlinearly on its length such that the force of a compressed spring is much higher than the force of a spring elongated by the same amount. This is equivalent to the incompressibility of liquids, making these springs suitable to simulate the polygonal oscillations. It is shown that this spring-mass network can imitate most of the characteristic features of experimentally reported polygonal oscillations. Additionally, it is shown that the network can execute certain dynamics which so far have not been observed in a perturbed liquid drop. The features of dynamics which are observed in the perturbed network are: polygonal oscillations, rotation of network, numerical relations (rational and irrational) between the frequencies of polygonal oscillations and the forcing signal, and the dependency of the shape of the polygons on the parameters of perturbation. △ Less

Submitted 20 September, 2021; v1 submitted 17 June, 2021; originally announced June 2021.

Magnetized ion-acoustic shock waves in degenerate quantum plasma

Authors: S. Jahan, T. S. Roy, B. E. Sharmin, N. A. Chowdhury, A. Mannan, A. A. Mamun

Abstract: A theoretical investigation has been carried out to examine the ion-acoustic shock waves (IASHWs) in a magnetized degenerate quantum plasma system containing inertialess ultra-relativistically degenerate electrons, and inertial non-relativistic positively charged heavy and light ions. The Burgers' equation is derived by employing reductive perturbation method. It can be seen that under considerati… ▽ More A theoretical investigation has been carried out to examine the ion-acoustic shock waves (IASHWs) in a magnetized degenerate quantum plasma system containing inertialess ultra-relativistically degenerate electrons, and inertial non-relativistic positively charged heavy and light ions. The Burgers' equation is derived by employing reductive perturbation method. It can be seen that under consideration of non-relativistic positively charged heavy and light ions, the plasma model supports only positive electrostatic shock structure. It is also observed that the charge state and number density of the non-relativistic heavy and light ions enhance the amplitude of IASHWs, and the steepness of the shock profile is decreased with ion kinematic viscosity ($ηいーた$). The findings of our present investigation will be helpful in understanding the nonlinear propagation of IASHWs in white dwarfs and neutron stars. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 8 figures, 6 pages

Modulational instability of dust-ion-acoustic waves in pair-ion plasma having non-thermal non-extensive electrons

Authors: M. K. Islam, A. A. Noman, J. Akter, N. A. Chowdhury, A. Mannan, T. S. Roy, M. Salahuddin, A. A. Mamun

Abstract: The modulational instability (MI) criteria of dust-ion-acoustic (DIA) waves (DIAWs) have been investigated in a four-component pair-ion plasma having inertial pair-ions, inertialess non-thermal non-extensive electrons, and immobile negatively charged massive dust grains. A nonlinear Schrödinger equation (NLSE) is derived by using reductive perturbation method. The nonlinear and dispersive coeffici… ▽ More The modulational instability (MI) criteria of dust-ion-acoustic (DIA) waves (DIAWs) have been investigated in a four-component pair-ion plasma having inertial pair-ions, inertialess non-thermal non-extensive electrons, and immobile negatively charged massive dust grains. A nonlinear Schrödinger equation (NLSE) is derived by using reductive perturbation method. The nonlinear and dispersive coefficients of the NLSE can predict the modulationally stable and unstable parametric regimes of DIAWs and associated first and second order DIA rogue waves (DIARWs). The MI growth rate and the configuration of the DIARWs are examined, and it is found that the MI growth rate increases (decreases) with increasing the number density of the negatively charged dust grains in the presence (absence) of the negative ions. It is also observed that the amplitude and width of the DIARWs increase (decrease) with the negative (positive) ion mass. The implications of the results to laboratory and space plasmas are briefly discussed. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 12 pages and 08 figures

First and second-order dust-ion-acoustic rogue waves in non-thermal plasma

Authors: S. Banik, R. K. Shikha, A. A. Noman, N. A. Chowdhury, A. Mannan, T. S. Roy, A. A. Mamun

Abstract: A nonlinear Schrödinger equation (NLSE) has been derived by employing reductive perturbation method for investigating the modulational instability of dust-ion-acoustic waves (DIAWs) in a four-component plasma having stationary negatively charged dust grains, inertial warm ions, and inertialess non-thermal electrons and positrons. It is observed that under consideration, the plasma system supports… ▽ More A nonlinear Schrödinger equation (NLSE) has been derived by employing reductive perturbation method for investigating the modulational instability of dust-ion-acoustic waves (DIAWs) in a four-component plasma having stationary negatively charged dust grains, inertial warm ions, and inertialess non-thermal electrons and positrons. It is observed that under consideration, the plasma system supports both modulationally stable and unstable domains, which are determined by the sign of the dispersive and nonlinear coefficients of NLSE, of the DIAWs. It is also found that the nonlinearity as well as the height and width of the first and second-order rogue waves increases with the non-thermality of electron and positron. The relevancy of our present investigation to the observations in space plasmas is pinpointed. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 5 pages, 5 figures

Conformation Controllable Inelastic Charge Transport and Shot Noise Behavior in Metal-String Single Molecular Devices

Authors: Talem Rebeda Roy, Arijit Sen

Abstract: It is often intriguing experimentally to take stock of how conformational changes in the device configuration may impact the overall charge transport behavior of single-molecule junctions. Based on the allied approach of density functional theory and non-equilibrium Green's function formalism, we explore here the effect of junction heterogeneity on inelastic charge transport in various metal-strin… ▽ More It is often intriguing experimentally to take stock of how conformational changes in the device configuration may impact the overall charge transport behavior of single-molecule junctions. Based on the allied approach of density functional theory and non-equilibrium Green's function formalism, we explore here the effect of junction heterogeneity on inelastic charge transport in various metal-string based single-molecule devices. The constituent active elements being sensitive to the resonant levels, transition metal centers are found to influence stretching, bending, and torsional excitation modes, while rocking and scissoring modes are controlled largely by the axial ligands. For certain molecular conformations and electrode orientations, phonon-assisted quantum interference effect may crop up, leading to the suppression of higher wavenumber vibrational modes. The resulting inelastic spectra are likely to take the shape of dominant Fano resonance or anti-resonance, depending on whether phonons are emitted or absorbed. Such nanoscale quantum interference effect is manifested especially in those metal-string molecular junctions for which the energy gap (between localized and delocalized virtual states) lies well within the optical phonon energies ($Δでるた{E}_{|HOMO-LUMO|} <$ 40 meV). It also turns out that single molecular shot noise can exhibit nearly Poissonian behavior if the inter-channel tunneling through frontier orbitals is accompanied by phonon absorption or emission following a slow relaxation process. Our results thus suggest that charge transport properties across metal-string complexes can be potentially tuned by selective architecture of the metal centers and also, by preferred orientation of nanoscale electrodes in a bid to build up molecular devices with desirable controllability. △ Less

Submitted 21 January, 2020; v1 submitted 20 January, 2020; originally announced January 2020.

Comments: 8 figures

Journal ref: Published in Applied Surface Science, Volume 508, pp 145196 (2020)

A constrained pressure-temperature residual (CPTR) method for non-isothermal multiphase flow in porous media

Authors: Thomas Roy, Tom B Jönsthövel, Christopher Lemon, Andrew J Wathen

Abstract: For both isothermal and thermal petroleum reservoir simulation, the Constrained Pressure Residual (CPR) method is the industry-standard preconditioner. This method is a two-stage process involving the solution of a restricted pressure system. While initially designed for the isothermal case, CPR is also the standard for thermal cases. However, its treatment of the energy conservation equation does… ▽ More For both isothermal and thermal petroleum reservoir simulation, the Constrained Pressure Residual (CPR) method is the industry-standard preconditioner. This method is a two-stage process involving the solution of a restricted pressure system. While initially designed for the isothermal case, CPR is also the standard for thermal cases. However, its treatment of the energy conservation equation does not incorporate heat diffusion, which is often dominant in thermal cases. In this paper, we present an extension of CPR: the Constrained Pressure-Temperature Residual (CPTR) method, where a restricted pressure-temperature system is solved in the first stage. In previous work, we introduced a block preconditioner with an efficient Schur complement approximation for a pressure-temperature system. Here, we extend this method for multiphase flow as the first stage of CPTR. The algorithmic performance of different two-stage preconditioners is evaluated for reservoir simulation test cases. △ Less

Submitted 27 May, 2020; v1 submitted 5 July, 2019; originally announced July 2019.

Comments: 28 pages, 2 figures. Sources/sinks description in arXiv:1902.00095

MSC Class: 65F08; 65F10; 76S05; 65M08; 65M60

ChAKRA : The high resolution charged particle detector array at VECC

Authors: Samir Kundu, T. K. Rana, C. Bhattacharya, K. Banerjee, R. Pandey, Santu Manna J. K. Meena, A. K. Saha, J. K. Sahoo, P. Dhara A. Dey D. Gupta T. K. Ghosh Pratap Roy, G. Mukherjee, R Mandal Saha, S. Roy, S. R. Bajirao, A. Sen, S. Bhattacharya

Abstract: A large 4$πぱい$ array of charged particle detectors has been developed at Variable Energy Cyclotron Centre to facilitate high resolution charged particle reaction and spectroscopy studies by detecting event-by-event the charged reaction products emitted in heavy ion reactions at energy $\sim$ 10-60 MeV/A. The forward part ($θしーた\sim \pm $ $7^{0}$ - $\pm 45^{0}$) of the array consists of 24 highly granul… ▽ More A large 4$πぱい$ array of charged particle detectors has been developed at Variable Energy Cyclotron Centre to facilitate high resolution charged particle reaction and spectroscopy studies by detecting event-by-event the charged reaction products emitted in heavy ion reactions at energy $\sim$ 10-60 MeV/A. The forward part ($θしーた\sim \pm $ $7^{0}$ - $\pm 45^{0}$) of the array consists of 24 highly granular, high resolution charged particle telescopes, each of which is made by three layers [single sided silicon strip($Δでるた$E) + double sided silicon strip (E/$Δでるた$E) + CsI(Tl)(E)]of detectors. The backward part of the array consists of 112 CsI(Tl) detectors which are capable of detecting primarily the light charged particles (Z $\le$ 2) emitted in the angular range of $θしーた\sim \pm $ $45^{0}$ - $\pm 175^{0}$. The extreme forward part of the array ($θしーた\sim \pm $ $3^{0}$ - $\pm 7^{0}$) is made up of 32 slow-fast plastic phoswich detectors that are capable of detecting light (Z $\le$2) and heavy charged particles (3 $\le$ Z $\lesssim$ 20) as well as handling high count rates. The design, construction and characterization of the array has been described. △ Less

Submitted 11 April, 2019; originally announced April 2019.

Comments: Paper has been submitted in NIMA

Reconstruction of Hydraulic Data by Machine Learning

Authors: Corentin J. Lapeyre, Nicolas Cazard, Pamphile T. Roy, Sophie Ricci, Fabrice Zaoui

Abstract: Numerical simulation models associated with hydraulic engineering take a wide array of data into account to produce predictions: rainfall contribution to the drainage basin (characterized by soil nature, infiltration capacity and moisture), current water height in the river, topography, nature and geometry of the river bed, etc. This data is tainted with uncertainties related to an imperfect knowl… ▽ More Numerical simulation models associated with hydraulic engineering take a wide array of data into account to produce predictions: rainfall contribution to the drainage basin (characterized by soil nature, infiltration capacity and moisture), current water height in the river, topography, nature and geometry of the river bed, etc. This data is tainted with uncertainties related to an imperfect knowledge of the field, measurement errors on the physical parameters calibrating the equations of physics, an approximation of the latter, etc. These uncertainties can lead the model to overestimate or underestimate the flow and height of the river. Moreover, complex assimilation models often require numerous evaluations of physical solvers to evaluate these uncertainties, limiting their use for some real-time operational applications. In this study, we explore the possibility of building a predictor for river height at an observation point based on drainage basin time series data. An array of data-driven techniques is assessed for this task, including statistical models, machine learning techniques and deep neural network approaches. These are assessed on several metrics, offering an overview of the possibilities related to hydraulic time-series. An important finding is that for the same hydraulic quantity, the best predictors vary depending on whether the data is produced using a physical model or real observations. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: Submitted to SimHydro 2019

A block preconditioner for non-isothermal flow in porous media

Authors: Thomas Roy, Tom B. Jönsthövel, Christopher Lemon, Andrew J. Wathen

Abstract: In petroleum reservoir simulation, the industry standard preconditioner, the constrained pressure residual method (CPR), is a two-stage process which involves solving a restricted pressure system with Algebraic Multigrid (AMG). Initially designed for isothermal models, this approach is often used in the thermal case. However, it does not have a specific treatment of the additional energy conservat… ▽ More In petroleum reservoir simulation, the industry standard preconditioner, the constrained pressure residual method (CPR), is a two-stage process which involves solving a restricted pressure system with Algebraic Multigrid (AMG). Initially designed for isothermal models, this approach is often used in the thermal case. However, it does not have a specific treatment of the additional energy conservation equation and temperature variable. We seek to develop preconditioners which better capture thermal effects such as heat diffusion. In order to study the effects of both pressure and temperature on fluid and heat flow, we consider a model of non-isothermal single phase flow through porous media. For this model, we develop a block preconditioner with an efficient Schur complement approximation. Both the pressure block and the approximate Schur complement are approximately inverted using an AMG V-cycle. The resulting solver is scalable with respect to problem size and parallelization. △ Less

Submitted 17 June, 2019; v1 submitted 31 January, 2019; originally announced February 2019.

Comments: 35 pages, 3 figures

MSC Class: 65F08; 65F10; 76S05; 65M08; 65M60

Resampling Strategies to Improve Surrogate Model-based Uncertainty Quantification - Application to LES of LS89

Authors: Pamphile Tupui Roy, Luis Miguel Segui, Jean-Christophe Jouhaud, Laurent Gicquel

Abstract: Uncertainty Quantification (UQ) is receiving more and more attention for engineering applications in particular from robust optimization. Indeed, running a computer experiment only provides a limited knowledge in terms of uncertainty and variability of the input parameters. These experiments are often computationally expensive and surrogate models can be constructed to address this issue. The outc… ▽ More Uncertainty Quantification (UQ) is receiving more and more attention for engineering applications in particular from robust optimization. Indeed, running a computer experiment only provides a limited knowledge in terms of uncertainty and variability of the input parameters. These experiments are often computationally expensive and surrogate models can be constructed to address this issue. The outcome of a UQ study is in this case directly correlated to the surrogate's quality. Thus, attention must be devoted to the Design of Experiments (DoE) to retrieve as much information as possible. This work presents two new strategies for parameter space resampling to improve a Gaussian Process (GP) surrogate model. These techniques indeed show an improvement of the predictive quality of the model with high dimensional analytical input functions. Finally, the methods are successfully applied to a turbine blade Large Eddy Simulation application: the aerothermal flow around the LS89 blade cascade. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: Accepted in International Journal for Numerical Methods in Fluids

Dynamic metasurface lens based on MEMS Technology

Authors: Tapashree Roy, Shuyan Zhang, Il Woong Jung, Mariano Troccoli, Federico Capasso, Daniel Lopez

Abstract: In the recent years, metasurfaces, being flat and lightweight, have been designed to replace bulky optical components with various functions. We demonstrate a monolithic Micro-Electro-Mechanical System (MEMS) integrated with a metasurface-based flat lens that focuses light in the mid-infrared spectrum. A two-dimensional scanning MEMS platform controls the angle of the lens along the two orthogonal… ▽ More In the recent years, metasurfaces, being flat and lightweight, have been designed to replace bulky optical components with various functions. We demonstrate a monolithic Micro-Electro-Mechanical System (MEMS) integrated with a metasurface-based flat lens that focuses light in the mid-infrared spectrum. A two-dimensional scanning MEMS platform controls the angle of the lens along the two orthogonal axes (tip-tilt) by +-9 degrees, thus enabling dynamic beam steering. The device can compensate for off-axis incident light and thus correct for aberrations such as coma. We show that for low angular displacements, the integrated lens-on-MEMS system does not affect the mechanical performance of the MEMS actuators and preserves the focused beam profile as well as the measured full width at half maximum. We envision a new class of flat optical devices with active control provided by the combination of metasurfaces and MEMS for a wide range of applications, such as miniaturized MEMS-based microscope systems, LIDAR scanners, and projection systems. △ Less

Submitted 10 December, 2017; originally announced December 2017.

Journal ref: APL Photonics 3, 021302 (2018)

Extremely broadband ultralight thermally emissive metasurfaces

Authors: Ali Naqavi, Samuel P. Loke, Michael D. Kelzenberg, Dennis M. Callahan, Emily C. Warmann, Pilar Espinet-González, Nina Vaidya, Tatiana A. Roy, Jing-Shun Huang, Tatiana G. Vinogradova, Alexander J. Messer, Harry A. Atwater

Abstract: We report the design, fabrication and characterization of ultralight highly emissive metaphotonic structures with record-low mass/area that emit thermal radiation efficiently over a broad spectral (2 to 35 microns) and angular (0-60 degrees) range. The structures comprise one to three pairs of alternating nanometer-scale metallic and dielectric layers, and have measured effective 300 K hemispheric… ▽ More We report the design, fabrication and characterization of ultralight highly emissive metaphotonic structures with record-low mass/area that emit thermal radiation efficiently over a broad spectral (2 to 35 microns) and angular (0-60 degrees) range. The structures comprise one to three pairs of alternating nanometer-scale metallic and dielectric layers, and have measured effective 300 K hemispherical emissivities of 0.7 to 0.9. To our knowledge, these structures, which are all subwavelength in thickness are the lightest reported metasurfaces with comparable infrared emissivity. The superior optical properties, together with their mechanical flexibility, low outgassing, and low areal mass, suggest that these metasurfaces are candidates for thermal management in applications demanding of ultralight flexible structures, including aerospace applications, ultralight photovoltaics, lightweight flexible electronics, and textiles for thermal insulation. △ Less

Submitted 8 October, 2017; originally announced October 2017.

Comparison of Polynomial Chaos and Gaussian Process surrogates for uncertainty quantification and correlation estimation of spatially distributed open-channel steady flows

Authors: Pamphile Tupui Roy, Nabil El Moçayd, Sophie Ricci, Jean-Christophe Jouhaud, Nicole Goutal, Matthias De Lozzo, Mélanie C. Rochoux

Abstract: Data assimilation is widely used to improve flood forecasting capability, especially through parameter inference requiring statistical information on the uncertain input parameters (upstream discharge, friction coefficient) as well as on the variability of the water level and its sensitivity with respect to the inputs. For particle filter or ensemble Kalman filter, stochastically estimating probab… ▽ More Data assimilation is widely used to improve flood forecasting capability, especially through parameter inference requiring statistical information on the uncertain input parameters (upstream discharge, friction coefficient) as well as on the variability of the water level and its sensitivity with respect to the inputs. For particle filter or ensemble Kalman filter, stochastically estimating probability density function and covariance matrices from a Monte Carlo random sampling requires a large ensemble of model evaluations, limiting their use in real-time application. To tackle this issue, fast surrogate models based on Polynomial Chaos and Gaussian Process can be used to represent the spatially distributed water level in place of solving the shallow water equations. This study investigates the use of these surrogates to estimate probability density functions and covariance matrices at a reduced computational cost and without the loss of accuracy, in the perspective of ensemble-based data assimilation. This study focuses on 1-D steady state flow simulated with MASCARET over the Garonne River (South-West France). Results show that both surrogates feature similar performance to the Monte-Carlo random sampling, but for a much smaller computational budget; a few MASCARET simulations (on the order of 10-100) are sufficient to accurately retrieve covariance matrices and probability density functions all along the river, even where the flow dynamic is more complex due to heterogeneous bathymetry. This paves the way for the design of surrogate strategies suitable for representing unsteady open-channel flows in data assimilation. △ Less

Submitted 17 October, 2017; v1 submitted 3 May, 2017; originally announced May 2017.

Luminescence engineering in plasmonic meta-surfaces

Authors: Tapashree Roy, Edward T. F. Rogers, Nikolay I. Zheludev

Abstract: Photoluminescence is a phenomenon of significant interest due to its wide range of technological applications in plasmonics, nanolasers, spasers, lasing spasers, loss compensation and gain in metamaterials, and luminescent media. Nanostructured materials are known to have very different luminescence characteristics to bulk samples or planar films. Here we show that by engineering a nanostructured… ▽ More Photoluminescence is a phenomenon of significant interest due to its wide range of technological applications in plasmonics, nanolasers, spasers, lasing spasers, loss compensation and gain in metamaterials, and luminescent media. Nanostructured materials are known to have very different luminescence characteristics to bulk samples or planar films. Here we show that by engineering a nanostructured meta-surface, we can choose the position of photoluminescence absorption and emission lines of thin gold films. The nanostructuring also aids to strong enhancement of the emission from gold, by a factor of 76 in our experiments. This enhancement is determined by the relative position of the engineered absorption and emission lines to the exciting laser wavelength and the intrinsic properties of the constituent material. These luminescence-engineered materials combined with a resonant material, as in the lasing spaser, or with the power of reconfigurable metamaterials promise huge potential as tunable nanoscale light sources. △ Less

Submitted 10 June, 2016; originally announced June 2016.

Shape Evolution of Sandwitched Droplet in Microconfined Shear Flow

Authors: Kaustav Chaudhury, Debabrata DasGupta, Tamal Roy, Suman Chakraborty

Abstract: Droplets confined in a microfluidic channel often exhibit intriguing shapes, primarily attributable to complex hydrodynamic interactions over small scales. We show that effect of varied substrate wettability conditions may further complicate these interactions remarkably, and often non-trivially. Our studies reveal that the combined influence of substrate wettability and fluidic confinement eventu… ▽ More Droplets confined in a microfluidic channel often exhibit intriguing shapes, primarily attributable to complex hydrodynamic interactions over small scales. We show that effect of varied substrate wettability conditions may further complicate these interactions remarkably, and often non-trivially. Our studies reveal that the combined influence of substrate wettability and fluidic confinement eventually culminates towards influencing the droplet transients, distortion of the local shear flow field, as well as drop stabilization against breakup or detachment, allowing one to develop different regimes of shapes evolution that are fascinatingly distinct from the ones reported in earlier studies on drop breakup in micro-confined shear flows. We further demonstrate that combined consequences of wall effects and interfacial wettability characteristics can be exploited to pattern microfluidic substrates with pre-designed patches, bearing far-ranging scientific and technological consequences in several scientific and technological applications of contemporary relevance. △ Less

Submitted 10 September, 2015; originally announced September 2015.

Comments: Article: 14 pages, 3 figures; Supplementary: 2 pages, 1 figure

Sub-wavelength focusing meta-lens

Authors: Tapashree Roy, Edward T. F. Rogers, Nikolay I. Zheludev

Abstract: We show that planar a plasmonic metamaterial with spatially variable meta-atom parameters can focus transmitted light into sub-wavelength hot-spots located beyond the near-field of the metamaterial. By nano-structuring a gold film we created an array of meta-lenses generating foci of 160 nm (0.2λらむだ) in diameter when illuminated by a wavelength of 800 nm. We attribute the occurrence of sub-wavelength… ▽ More We show that planar a plasmonic metamaterial with spatially variable meta-atom parameters can focus transmitted light into sub-wavelength hot-spots located beyond the near-field of the metamaterial. By nano-structuring a gold film we created an array of meta-lenses generating foci of 160 nm (0.2λらむだ) in diameter when illuminated by a wavelength of 800 nm. We attribute the occurrence of sub-wavelength hotspots beyond the near field to the phenomenon of superoscillation. △ Less

Submitted 7 November, 2012; originally announced November 2012.