Researcher profile

K. B. Nakshatrala

K. B. Nakshatrala contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

How to pose material design problems for flow through porous media applications?: Sensitivity of dissipation rate to medium's permeability holds the key

Recent studies have advocated using the total dissipation rate under topology optimization to realize material designs involving the flow of fluids through porous media. However, these studies decided how to pose the design problem, such as maximizing the total dissipation rate for some situations while minimizing for others, by solving one-dimensional problems and justifying their choices using numerical experiments. The rigor is lacking -- a bottleneck for further scientific advancements to computational material design. This paper provides the missing theoretical justification. We identify four classes of boundary value problems using the adjoint state method and analytically calculate the sensitivity of the total dissipation rate to the permeability field. For two of those classes in which the flow of fluids is pressure-driven, the sensitivity is positive -- the total dissipation rate increases if the medium's permeability increases. While for the other two classes, in which the flow is velocity-driven, the trend is the opposite. These sensitivities provide rigorous answers to the central question: how to pose a material design problem for flow through porous media applications. The impact of our work is multi-fold. First, this study further elevates the role of the dissipation rate in posing well-posed material design problems using topology optimization. Second, besides the theoretical significance, the results benefit computational scientists and practitioners to realize optimal designs. Third, given their simplicity yet far-reaching impact, both the approach and results possess immense pedagogical value.

preprint2013arXiv

A framework for coupling flow and deformation of the porous solid

In this paper, we consider the flow of an incompressible fluid in a deformable porous solid. We present a mathematical model using the framework offered by the theory of interacting continua. In its most general form, this framework provides a mechanism for capturing multiphase flow, deformation, chemical reactions and thermal processes, as well as interactions between the various physics in a conveniently implemented fashion. To simplify the presentation of the framework, results are presented for a particular model than can be seen as an extension of Darcy's equation (which assumes that the porous solid is rigid) that takes into account elastic deformation of the porous solid. The model also considers the effect of deformation on porosity. We show that using this model one can recover identical results as in the framework proposed by Biot and Terzaghi. Some salient features of the framework are as follows: (a) It is a consistent mixture theory model, and adheres to the laws and principles of continuum thermodynamics, (b) the model is capable of simulating various important phenomena like consolidation and surface subsidence, and (c) the model is amenable to several extensions. We also present numerical coupling algorithms to obtain coupled flow-deformation response. Several representative numerical examples are presented to illustrate the capability of the mathematical model and the performance of the computational framework.

preprint2013arXiv

A numerical methodology for enforcing maximum principles and the non-negative constraint for transient diffusion equations

Transient diffusion equations arise in many branches of engineering and applied sciences (e.g., heat transfer and mass transfer), and are parabolic partial differential equations. It is well-known that, under certain assumptions on the input data, these equations satisfy important mathematical properties like maximum principles and the non-negative constraint, which have implications in mathematical modeling. However, existing numerical formulations for these types of equations do not, in general, satisfy maximum principles and the non-negative constraint. In this paper, we present a methodology for enforcing maximum principles and the non-negative constraint for transient anisotropic diffusion equation. The method of horizontal lines (also known as the Rothe method) is applied in which the time is discretized first. This results in solving steady anisotropic diffusion equation with decay equation at every discrete time level. The proposed methodology for transient anisotropic diffusion equation will satisfy maximum principles and the non-negative constraint on general computational grids, and with no additional restrictions on the time step. We illustrate the performance and accuracy of the proposed formulation using representative numerical examples. We also perform numerical convergence of the proposed methodology. For comparison, we also present the results from the standard single-field semi-discrete formulation and the results from a popular software package, which all will violate maximum principles and the non-negative constraint.

preprint2013arXiv

Modification to Darcy model for high pressure and high velocity applications and associated mixed finite element formulations

The Darcy model is based on a plethora of assumptions. One of the most important assumptions is that the Darcy model assumes the drag coefficient to be constant. However, there is irrefutable experimental evidence that viscosities of organic liquids and carbon-dioxide depend on the pressure. Experiments have also shown that the drag varies nonlinearly with respect to the velocity at high flow rates. In important technological applications like enhanced oil recovery and geological carbon-dioxide sequestration, one encounters both high pressures and high flow rates. It should be emphasized that flow characteristics and pressure variation under varying drag are both quantitatively and qualitatively different from that of constant drag. Motivated by experimental evidence, we consider the drag coefficient to depend on both the pressure and velocity. We consider two major modifications to the Darcy model based on the Barus formula and Forchheimer approximation. The proposed modifications to the Darcy model result in nonlinear partial differential equations, which are not amenable to analytical solutions. To this end, we present mixed finite element formulations based on least-squares formalism and variational multiscale formalism for the resulting governing equations. The proposed modifications to the Darcy model and its associated finite element formulations are used to solve realistic problems with relevance to enhanced oil recovery. We also study the competition between the nonlinear dependence of drag on the velocity and the dependence of viscosity on the pressure. To the best of the authors' knowledge such a systematic study has not been performed.

preprint2010arXiv

A numerical study of fluids with pressure dependent viscosity flowing through a rigid porous medium

In this paper we consider modifications to Darcy's equation wherein the drag coefficient is a function of pressure, which is a realistic model for technological applications like enhanced oil recovery and geological carbon sequestration. We first outline the approximations behind Darcy's equation and the modifications that we propose to Darcy's equation, and derive the governing equations through a systematic approach using mixture theory. We then propose a stabilized mixed finite element formulation for the modified Darcy's equation. To solve the resulting nonlinear equations we present a solution procedure based on the consistent Newton-Raphson method. We solve representative test problems to illustrate the performance of the proposed stabilized formulation. One of the objectives of this paper is also to show that the dependence of viscosity on the pressure can have a significant effect both on the qualitative and quantitative nature of the solution.

preprint2010arXiv

A stabilized mixed formulation for unsteady Brinkman equation based on the method of horizontal lines

In this paper, we present a stabilized mixed formulation for unsteady Brinkman equation. The formulation is systematically derived based on the variational multiscale formalism and the method of horizontal lines. The derivation does not need the assumption that the fine-scale variables do not depend on the time, which is the case with the conventional derivation of multiscale stabilized formulations for transient mixed problems. An expression for the stabilization parameter is obtained in terms of a bubble function, and appropriate bubble functions for various finite elements are also presented. Under the proposed formulation, equal-order interpolation for the velocity and pressure (which is computationally the most convenient) is stable. Representative numerical results are presented to illustrate the performance of the proposed formulation. Spatial and temporal convergence studies are also performed, and the proposed formulation performed well.