Search NASA⌕ Search

Engineering topics

Tokman, Mayya

Publications and source records attributed to Tokman, Mayya.

7 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Efficient Integration of Algebraic Constraint for Exponential Time Integration

This work is a continuation of a previous project where the high efficiency of exponential integration methods for DRMHD systems was demonstrated. Often algebraic constraints must also be enforced for DRMHD systems of interest. Straightforward application of exponential methods to DRMHD equations with constraints leads to prohibitively computationally expensive methods. In this work, we propose new exponential schemes that allow the constraints to be removed from evaluation of exponential matrix functions which drastically reduces computational complexity by eliminating the need to perform computations enforcing constraints in exponential calculations while still preserving a high order of accuracy and allowing for a large time step even when the problem is stiff. This idea is similar to the W-methods and is achieved by carefully designing a method with the desired order of accuracy, even with an incomplete Jacobian matrix used as an argument of exponential-like functions. The constraints are accounted for by including them in the evaluation of the right-hand-side forcing function of the spatially discretized system. We study performance of the new methods on test problems and outline future research directions that this work opens.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Exponential integrators for non-linear diffusion

The goal of this project is to compare the performance of exponential time integrators with traditional methods such as diagonally implicit Runge-Kutta methods in the context of solving the system of reduced magnetohydrodynamics (RMHD). In this report, we present initial results of a proof of concept study that shows that exponential integrators can be an efficient alternative to traditional integration schemes.

42 ENGINEERING↗
DOE OSTI2016

On the performance of exponential integrators for problems in magnetohydrodynamics

Exponential integrators have been introduced as an efficient alternative to explicit and implicit methods for integrating large stiff systems of differential equations. Over the past decades these methods have been studied theoretically and their performance was evaluated using a range of test problems. While the results of these investigations showed that exponential integrators can provide significant computational savings, the research on validating this hypothesis for large scale systems and understanding what classes of problems can particularly benefit from the use of the new techniques is in its initial stages. Resistive magnetohydrodynamic (MHD) modeling is widely used in studying large scale behavior of laboratory and astrophysical plasmas. In many problems numerical solution of MHD equations is a challenging task due to the temporal stiffness of this system in the parameter regimes of interest. In this paper we evaluate the performance of exponential integrators on large MHD problems and compare them to a state-of-the-art implicit time integrator. Both the variable and constant time step exponential methods of EPIRK-type are used to simulate magnetic reconnection and the Kevin–Helmholtz instability in plasma. Performance of these methods, which are part of the EPIC software package, is compared to the variable time step variable order BDF scheme included in the CVODE (part of SUNDIALS) library. We study performance of the methods on parallel architectures and with respect to magnitudes of important parameters such as Reynolds, Lundquist, and Prandtl numbers. We find that the exponential integrators provide superior or equal performance in most circumstances and conclude that further development of exponential methods for MHD problems is warranted and can lead to significant computational advantages for large scale stiff systems of differential equations such as MHD.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Electric Field-Driven Water Dipoles: Nanoscale Architecture of Electroporation

Electroporation is the formation of permeabilizing structures in the cell membrane under the influence of an externally imposed electric field. The resulting increased permeability of the membrane enables a wide range of biological applications, including the delivery of normally excluded substances into cells. While electroporation is used extensively in biology, biotechnology, and medicine, its molecular mechanism is not well understood. This lack of knowledge limits the ability to control and fine-tune the process. In this article we propose a novel molecular mechanism for the electroporation of a lipid bilayer based on energetics analysis. Using molecular dynamics simulations we demonstrate that pore formation is driven by the reorganization of the interfacial water molecules. Our energetics analysis and comparisons of simulations with and without the lipid bilayer show that the process of poration is driven by field-induced reorganization of water dipoles at the water-lipid or water-vacuum interfaces into more energetically favorable configurations, with their molecular dipoles oriented in the external field. Although the contributing role of water in electroporation has been noted previously, here we propose that interfacial water molecules are the main players in the process, its initiators and drivers. The role of the lipid layer, to a first-order approximation, is then reduced to a relatively passive barrier. This new view of electroporation simplifies the study of the problem, and opens up new opportunities in both theoretical modeling of the process and experimental research to better control or to use it in new, innovative ways.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗