Source author record

David Salac

David Salac 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

11works
6topics
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

11 published item(s)

preprint2020arXiv

Non-equilibrium dynamics of initially spherical vesicles in general flow

Many vesicles have a spherical resting shape and exposure to fluid flows induces an exchange between sub-optical area and visible (systematic) deformation, while the total area is conserved. The dynamics which controls the exchange between sub-optical and visible area depends on membrane properties such as bending rigidity and initial tension. Conversely, observation of these dynamics can be used to determine the membrane properties. The goal of this work is to create a general numerical model which accounts for the exchange between sub-optical and visible area. Unlike prior modeling efforts, the model does not pre-assume a shape type, such as nearly-spherical, or applied flow field, allowing the model to capture a wider variety of flow conditions. Based on implicit interface tracking and using a volume-preserving multiphase Navier-Stokes solver, the model is compared to several experimental results, showing excellent agreement. It is used to explore regimes not possible with previously published models, such as the variable viscosity case, and how these system properties influence experimentally measurable parameters such as deformation parameter. By creating a more generalized framework for the modeling of vesicles with sub-optical area, it will now be possible to make predictions on vesicle material properties from a wider variety of experimental results.

preprint2016arXiv

Energy Variation of Soft Matter Interfaces

The variation of energies associated with soft matter interfaces where surface inhomogeneities are present. These energies include the total bending and splay energy, the variable surface tension energy, a coupling energy between the total curvature and an underlying surface concentration field, the energy due to an external field, and a phase segregation energy. When considering these energies the variation of material properties such a bending rigidity are taken into account, which results in more general variation expressions. These variations can be used to determine the equilibrium interface and concentration configuration or to determine the driving forces for non-equilibrium situations. While the focus of this work are energies associated with multicomponent vesicles, it can easily be extended to any soft matter interface.

preprint2016arXiv

Level Set Jet Schemes for Stiff Advection Equations: The SemiJet Method

Many interfacial phenomena in physical and biological systems are dominated by high order geometric quantities such as curvature. Here a semi-implicit method is combined with a level set jet scheme to handle stiff nonlinear advection problems. The new method offers an improvement over the semi-implicit gradient augmented level set method previously introduced by requiring only one smoothing step when updating the level set jet function while still preserving the underlying methods higher accuracy. Sample results demonstrate that accuracy is not sacrificed while strict time step restrictions can be avoided.

preprint2016arXiv

Vesicles in magnetic fields

Liposome vesicles tend to align with an applied magnetic field. This is due to the directional magnetic susceptibility difference of the lipids which form the membrane of these vesicles. In this work a model of liposome vesicles exposed to magnetic field is presented. Starting from the base energy of a lipid membrane in a magnetic field, the force applied to the surrounding fluids is derived. This force is then used to investigate the dynamics of vesicle in the presence of magnetic fields.

preprint2015arXiv

A General, Mass-Preserving Navier-Stokes Projection Method

The conservation of mass is common issue with multiphase fluid simulations. In this work a novel projection method is presented which conserves mass both locally and globally. The fluid pressure is augmented with a time-varying component which accounts for any global mass change. The resulting system of equations is solved using an efficient Schur-complement method. Using the proposed method four numerical examples are performed: the evolution of a static bubble, the rise of a bubble, the breakup of a thin fluid thread, and the extension of a droplet in shear flow. The method is capable of conserving the mass even in situations with morphological changes such as droplet breakup.

preprint2015arXiv

Dynamics of Three-Dimensional Vesicles in DC Electric fields

A numerical and systematic parameter study of three-dimensional vesicle electrohydrodynamics is presented to investigate the effects of different fluid and membrane properties. The dynamics of vesicles in the presence of DC electric fields is considered, both in the presence and absence of linear shear flow. For suspended vesicles it is shown that the conductivity ratio and viscosity ratio between the interior and exterior fluids, as well as the vesicle membrane capacitance, substantially affect the minimum electric field strength required to induce a full Prolate-Oblate-Prolate transition.In addition, there exists a critical electric field strength above which a vesicle will no longer tumble when exposed to linear shear flow.

preprint2014arXiv

A numerical model for the trans-membrane voltage of vesicles

The Immersed Interface Method is employed to solve the time-varying electric field equations around a three-dimensional vesicle. To achieve second-order accuracy the implicit jump conditions for the electric potential, up to the second normal derivative, are derived. The trans-membrane potential is determined implicitly as part of the algorithm. The method is compared to an analytic solution based on spherical harmonics and verifies the second-order accuracy of the underlying discretization even in the presence of solution discontinuities. A sample result for an elliptic interface is also presented.

preprint2014arXiv

Electrohydrodynamics of Three-Dimensional Vesicles: A Numerical Approach

A three-dimensional numerical model of vesicle electrohydrodynamics in the presence of DC electric fields is presented. The vesicle membrane is modeled as a thin capacitive interface through the use of a semi-implicit level set Jet scheme. The enclosed volume and surface area are conserved both locally and globally by a new Navier-Stokes projection method. The electric field calculations explicitly take into account the capacitive interface by an implicit Immersed Interface Method formulation, which calculates the electric potential field and the trans-membrane potential simultaneously. The results match well with previously published experimental, analytic and two-dimensional computational works.

preprint2013arXiv

The Three-Dimensional Jump Conditions for the Stokes Equations with Discontinuous Viscosity, Singular Forces, and an Incompressible Interface

The three-dimensional jump conditions for the pressure and velocity fields, up to the second normal derivative,across an incompressible/inextensible interface in the Stokes regime are derived herein. The fluid viscosity is only piecewise continuous in the domain while the embedded interface exerts singular forces on the surround fluids. This gives rise to discontinuous solutions in the pressure and velocity field. These jump conditions are required to develop accurate numerical methods, such as the Immersed Interface Method, for the solutions of the Stokes equations in such situations.

preprint2012arXiv

The Semi Implicit Gradient Augmented Level Set Method

Here a semi-implicit formulation of the gradient augmented level set method is presented. By tracking both the level set and it's gradient accurate subgrid information is provided,leading to highly accurate descriptions of a moving interface. The result is a hybrid Lagrangian-Eulerian method that may be easily applied in two or three dimensions. The new approach allows for the investigation of interfaces evolving by mean curvature and by the intrinsic Laplacian of the curvature. In this work the algorithm, convergence and accuracy results are presented. Several numerical experiments in both two and three dimensions demonstrate the stability of the scheme.

preprint2011arXiv

The Augmented Fast Marching Method for Level Set Reinitialization

Including derivative information in the modelling of moving interfaces has been proposed as one method to increase the accuracy of numerical schemes with minimal additional cost. Here a new level set reinitialization technique using the fast marching method is presented. This augmented fast marching method will calculate the signed distance function and up to the second-order derivatives of the signed distance function for arbitrary interfaces. In addition to enforcing the condition $|\nablaϕ|^2=1$, where $ϕ$ is the level set function, the method ensures that $\nabla(|\nablaϕ|)^2=0$ and $\nabla\nabla(|\nablaϕ|)^2=0$ are also satisfied. Results indicate that for both two- and three-dimensional interfaces the resulting level set and curvature field are smooth even for coarse grids. Convergence results show that using first-order upwind derivatives and the augmented fast marching method result in a second-order accurate level set and gradient field and a first-order accurate curvature field.