Source author record

Huaxiong Huang

Huaxiong Huang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

9 published item(s)

preprint2022arXiv

A Bubble Model for the Gating of Kv Channels

Voltage-gated Kv channels play fundamental roles in many biological processes, such as the generation of the action potential. The gating mechanism of Kv channels is characterized experimentally by single-channel recordings and ensemble properties of the channel currents. In this work, we propose a bubble model coupled with a Poisson-Nernst-Planck (PNP) system to capture the key characteristics, particularly the delay in the opening of channels. The coupled PNP system is solved numerically by a finite-difference method and the solution is compared with an analytical approximation. We hypothesize that the stochastic behaviour of the gating phenomenon is due to randomness of the bubble and channel sizes. The predicted ensemble average of the currents under various applied voltages across the channels is consistent with experimental observations, and the Cole-Moore delay is captured by varying the holding potential.

preprint2022arXiv

Mathematical Model for Chemical Reactions in Electrolyte Applied to Cytochrome $c$ Oxidase: an Electro-osmotic Approach

A mathematical model for chemical reactions in electrolytes is developed using an Energy variational method consistent with classical thermodynamics. Electrostatics and chemical reactions are included in properly defined energetic and dissipative functionals. The energy variation method is generalized to deal with open systems with inputs of charge, mass, and energy. The open systems can transform input energy of one type into output energy of another type. The generalized method is used to analyze the conversion of electrical current into proton flow by cytochrome c oxidase, an important enzyme in mitochondria that helps generate the 'energetic currency of life' ATP. The structure of the oxidase guides flows of current and mass that interact according to the energetic and dissipative functionals of the generalized theory. Kirchhoff's current law provides important coupling when the enzyme is in its natural setting in the mitochondrial membrane. The natural function of the enzyme is the result. Electron flows are converted into proton flows and gradients.

preprint2021arXiv

A phase field model for mass transport with semi-permeable interfaces

In this paper, a thermal-dynamical consistent model for mass transfer across permeable moving interfaces is proposed by using the energy variation method. We consider a restricted diffusion problem where the flux across the interface depends on its conductance and the difference of the concentration on each side. The diffusive interface phase-field framework used here has several advantages over the sharp interface method. First of all, explicit tracking of the interface is no longer necessary. Secondly, the interfacial condition can be incorporated with a variable diffusion coefficient. A detailed asymptotic analysis confirms the diffusive interface model converges to the existing sharp interface model as the interface thickness goes to zero. A decoupled energy stable numerical scheme is developed to solve this system efficiently. Numerical simulations first illustrate the consistency of theoretical results on the sharp interface limit. Then a convergence study and energy decay test are conducted to ensure the efficiency and stability of the numerical scheme. To illustrate the effectiveness of our phase-field approach, several examples are provided, including a study of a two-phase mass transfer problem where drops with deformable interfaces are suspended in a moving fluid.

preprint2019arXiv

An energy stable $C^0$ finite element scheme for a quasi-incompressible phase-field model of moving contact line with variable density

In this paper, we focus on modeling and simulation of two-phase flow with moving contact lines and variable density. A thermodynamically consistent phase-field model with General Navier Boundary Condition is developed based on the concept of quasi-incompressibility and the energy variational method. Then a mass conserving and energy stable C0 finite element scheme is developed to solve the PDE system. Various numerical simulation results show that the proposed schemes are mass conservative, energy stable and the 2nd order for P1 element and 3rd order for P2 element convergence rate in the sense of L2 norm.

preprint2019arXiv

Numerical Method for Parameter Inference of Nonlinear ODEs with Partial Observations

Parameter inference of dynamical systems is a challenging task faced by many researchers and practitioners across various fields. In many applications, it is common that only limited variables are observable. In this paper, we propose a method for parameter inference of a system of nonlinear coupled ODEs with partial observations. Our method combines fast Gaussian process based gradient matching (FGPGM) and deterministic optimization algorithms. By using initial values obtained by Bayesian steps with low sampling numbers, our deterministic optimization algorithm is both accurate and efficient.

preprint2013arXiv

A mathematical model of the metabolic and perfusion effects on cortical spreading depression

Cortical spreading depression (CSD) is a slow-moving ionic and metabolic disturbance that propagates in cortical brain tissue. In addition to massive cellular depolarization, CSD also involves significant changes in perfusion and metabolism -- aspects of CSD that had not been modeled and are important to traumatic brain injury, subarachnoid hemorrhage, stroke, and migraine. In this study, we develop a mathematical model for CSD where we focus on modeling the features essential to understanding the implications of neurovascular coupling during CSD. In our model, the sodium-potassium--ATPase, mainly responsible for ionic homeostasis and active during CSD, operates at a rate that is dependent on the supply of oxygen. The supply of oxygen is determined by modeling blood flow through a lumped vascular tree with an effective local vessel radius that is controlled by the extracellular potassium concentration. We show that during CSD, the metabolic demands of the cortex exceed the physiological limits placed on oxygen delivery, regardless of vascular constriction or dilation. However, vasoconstriction and vasodilation play important roles in the propagation of CSD and its recovery. Our model replicates the qualitative and quantitative behavior of CSD -- vasoconstriction, oxygen depletion, extracellular potassium elevation, prolonged depolarization -- found in experimental studies. We predict faster, longer duration CSD in vivo than in vitro due to the contribution of the vasculature. Our results also help explain some of the variability of CSD between species and even within the same animal. These results have clinical and translational implications, as they allow for more precise in vitro, in vivo, and in silico exploration of a phenomenon broadly relevant to neurological disease.

preprint2012arXiv

A different perspective on retirement income sustainability: the blueprint for a ruin contingent life annuity (RCLA)

The purpose of this article is twofold. First, we motivate the need for a new type of stand-alone retirement income insurance product that would help individuals protect against personal longevity risk and possible "retirement ruin" in an economically efficient manner. We label this product a ruin-contingent life annuity (RCLA), which we elaborate-on and explain with various numerical examples and a basic pricing model. Second, we argue that with the proper perspective a similar product actually exists, albeit not available on a stand-alone basis. Namely, they are fused and embedded within modern variable annuity (VA) policies with guaranteed living income benefit (GLiB) riders. Indeed, the popularity of GLiB riders on VA policies point towards the potential commercial success of such a stand-alone vehicle.

preprint2012arXiv

Optimal retirement consumption with a stochastic force of mortality

We extend the lifecycle model (LCM) of consumption over a random horizon (a.k.a. the Yaari model) to a world in which (i.) the force of mortality obeys a diffusion process as opposed to being deterministic, and (ii.) a consumer can adapt their consumption strategy to new information about their mortality rate (a.k.a. health status) as it becomes available. In particular, we derive the optimal consumption rate and focus on the impact of mortality rate uncertainty vs. simple lifetime uncertainty -- assuming the actuarial survival curves are initially identical -- in the retirement phase where this risk plays a greater role. In addition to deriving and numerically solving the PDE for the optimal consumption rate, our main general result is that when utility preferences are logarithmic the initial consumption rates are identical. But, in a CRRA framework in which the coefficient of relative risk aversion is greater (smaller) than one, the consumption rate is higher (lower) and a stochastic force of mortality does make a difference. That said, numerical experiments indicate that even for non-logarithmic preferences, the stochastic mortality effect is relatively minor from the individual's perspective. Our results should be relevant to researchers interested in calibrating the lifecycle model as well as those who provide normative guidance (a.k.a. financial advice) to retirees.

preprint2008arXiv

The lipid bilayer at the mesoscale: a physical continuum model

We study a continuum model of the lipid bilayer based on minimizing the free energy of a mixture of water and lipid molecules. This paper extends previous work by Blom & Peletier (2004) in the following ways. (a) It formulates a more physical model of the hydrophobic effect to facilitate connections with microscale simulations. (b) It clarifies the meaning of the model parameters. (c) It outlines a method for determining parameter values so that physically-realistic bilayer density profiles can be obtained, for example for use in macroscale simulations. Points (a)-(c) suggest that the model has potential to robustly connect some micro- and macroscale levels of multiscale blood flow simulations. The mathematical modelling in point (a) is based upon a consideration of the underlying physics of inter-molecular forces. The governing equations thus obtained are minimized by gradient flows via a novel numerical approach; this enables point (b). The numerical results are shown to behave physically in terms of the effect of background concentration, in contrast to the earlier model which is shown here to not display the expected behaviour. A "short-tail" approximation of the lipid molecules also gives an analytical tool which yields critical values of some parameters under certain conditions. Point (c) involves the first quantitative comparison of the numerical data with physical experimental results.