Variational quantitative phase-field modeling of nonisothermal sintering process.

Phase-field modeling has become a powerful tool in describing the complex pore-structure evolution and the intricate multiphysics in nonisothermal sintering processes. However, the quantitative validity of conventional variational phase-field models involving diffusive processes is a challenge. Artificial interface effects, like the trapping effects, may originate at the interface when the kinetic properties of two opposing phases are different. On the other hand, models with prescribed antitrapping terms do not necessarily guarantee the thermodynamics variational nature of the model. This issue has been solved for liquid-solid interfaces via the development of the variational quantitative solidification phase-field model. However, there is no related work addressing the interfaces in nonisothermal sintering, where the free surfaces between the solid phase and surrounding pore regions exhibit strong asymmetry of mass and thermal properties. Also, additional challenges arise due to the conserved order parameter describing the free surfaces. In this work, we present a variational and quantitative phase-field model for nonisothermal sintering processes. The model is derived via an extended nondiagonal phase-field model. The model evolution equations have naturally cross-coupling terms between the conserved kinetics (i.e., mass and thermal transfer) and the nonconserved one (grain growth). These terms are shown via asymptotic analysis to be instrumental in ensuring the elimination of interface artifacts, while also examined to not modify the thermodynamic equilibrium condition (characterized by a dihedral angle). Moreover, we demonstrate that the trapping effects and the existence of surface diffusion in conservation laws are direction-dependent. An anisotropic interpolation scheme of the kinetic mobilities that differentiates between the normal and tangential directions along the interface is discussed. Numerically, we demonstrate the importance of the cross-couplings and the anisotropic interpolation by presenting thermal-microstructural evolutions.

Systematic derivation of hydrodynamic equations for viscoelastic phase separation.

We present a detailed derivation of a simple hydrodynamic two-fluid model, which aims at the description of the phase separation of non-entangled polymer solutions, where viscoelastic effects play a role. It is directly based upon the coarse-graining of a well-defined molecular model, such that all degrees of freedom have a clear and unambiguous molecular interpretation. The considerations are based upon a free-energy functional, and the dynamics is split into a conservative and a dissipative part, where the latter satisfies the Onsager relations and the second law of thermodynamics. The model is therefore fully consistent with both equilibrium and non-equilibrium thermodynamics. The derivation proceeds in two steps: firstly, we derive an extended model comprising two scalar and four vector fields, such that inertial dynamics of the macromolecules and of the relative motion of the two fluids is taken into account. In the second step, we eliminate these inertial contributions and, as a replacement, introduce phenomenological dissipative terms, which can be modeled easily by taking into account the principles of non-equilibrium thermodynamics. The final simplified model comprises the momentum conservation equation, which includes both interfacial and elastic stresses, a convection-diffusion equation where interfacial and elastic contributions occur as well, and a suitably convected relaxation equation for the end-to-end vector field. In contrast to the traditional two-scale description that is used to derive rheological equations of motion, we here treat the hydrodynamic and the macromolecular degrees of freedom on the same basis. Nevertheless, the resulting model is fairly similar, though not fully identical, to models that have been discussed previously. Notably, we find a rheological constitutive equation that differs from the standard Oldroyd-B model. Within the framework of kinetic theory, this difference may be traced back to a different underlying statistical-mechanical ensemble that is used for averaging the stress. To what extent the model is able to reproduce the full phenomenology of viscoelastic phase separation is presently an open question, which shall be investigated in the future.

Analysis of a viscoelastic phase separation model.

A new model for viscoelastic phase separation is proposed, based on a systematically derived conservative two-fluid model. Dissipative effects are included by phenomenological viscoelastic terms. By construction, the model is consistent with the second law of thermodynamics. We study well-posedness of the model in two space dimensions, i.e., existence of weak solutions, a weak-strong uniqueness principle, and stability with respect to perturbations, which are proven by means of relative energy estimates. Our numerical simulations based on the new viscoelastic phase separation model are in good agreement with physical experiments. Furthermore, a good qualitative agreement with mesoscopic simulations is observed.

Parametric Sequential Method for MRI-Based Wall Shear Stress Quantification.

Wall shear stress (WSS) has been suggested as a potential biomarker in various cardiovascular diseases and it can be estimated from phase-contrast Magnetic Resonance Imaging (PC-MRI) velocity measurements. We present a parametric sequential method for MRI-based WSS quantification consisting of a geometry identification and a subsequent approximation of the velocity field. This work focuses on its validation, investigating well controlled high-resolution in vitro measurements of turbulent stationary flows and physiological pulsatile flows in phantoms. Initial tests for in vivo 2D PC-MRI data of the ascending aorta of three volunteers demonstrate basic applicability of the method to in vivo.

MR-based wall shear stress measurements in fully developed turbulent flow using the Clauser plot method.

In arterial blood flow wall shear stress (WSS) quantifies the frictional force that flowing blood exerts on a vessel wall. WSS can be directly estimated from phase-contrast (PC) MR velocity measurements and has been suggested as a biomarker in cardio-vascular diseases. We present and investigate the application of the Clauser plot method for estimating WSS in fully developed turbulent stationary flow using PC velocity measurements. The Clauser plot method estimates WSS from the logarithmic region of boundary layer in fully developed turbulent stationary flow. The Clauser plot method was evaluated using 2D PC-MR phantom measurements at 3 T for different in-plane resolutions at various Reynolds numbers. WSS values derived from the Clauser plot were compared to results from Laser Doppler Velocimetry (LDV) measurements and theoretical results calculated using the friction factor formula for smooth pipe flow. For all Reynolds numbers, WSS values derived from the Clauser plot were in good agreement with results from LDV measurements and values using the friction factor formula (relative deviations ∼5%). Furthermore, Clauser plot derived results were almost independent of spatial resolution, in contrast to WSS results obtained with our in-house software tool for MR-based WSS quantification showing relative deviations of more than 100%. In fully developed turbulent flow, the Clauser plot method provides highly consistent WSS independent of the underlying spatial resolution. Therefore, it renders a valuable approach for MR-based WSS estimates in controllable flow settings. Although its direct in vivo applicability is severely limited because of the different flow character, it may serve as helpful approach for validation of MR-based WSS quantification algorithms prior to their clinical application.

High-performance image reconstruction in fluorescence tomography on desktop computers and graphics hardware.

Image reconstruction in fluorescence optical tomography is a three-dimensional nonlinear ill-posed problem governed by a system of partial differential equations. In this paper we demonstrate that a combination of state of the art numerical algorithms and a careful hardware optimized implementation allows to solve this large-scale inverse problem in a few seconds on standard desktop PCs with modern graphics hardware. In particular, we present methods to solve not only the forward but also the non-linear inverse problem by massively parallel programming on graphics processors. A comparison of optimized CPU and GPU implementations shows that the reconstruction can be accelerated by factors of about 15 through the use of the graphics hardware without compromising the accuracy in the reconstructed images.

Electrical forces determine glomerular permeability.

There is ongoing controversy about the mechanisms that determine the characteristics of the glomerular filter. Here, we tested whether flow across the glomerular filter generates extracellular electrical potential differences, which could be an important determinant of glomerular filtration. In micropuncture experiments in Necturus maculosus, we measured a potential difference across the glomerular filtration barrier that was proportional to filtration pressure (-0.045 mV/10 cm H2O). The filtration-dependent potential was generated without temporal delay and was negative within Bowman's space. Perfusion with the cationic polymer protamine abolished the potential difference. We propose a mathematical model that considers the relative contributions of diffusion, convection, and electrophoretic effects on the total flux of albumin across the filter. According to this model, potential differences of -0.02 to -0.05 mV can induce electrophoretic effects that significantly influence the glomerular sieving coefficient of albumin. This model of glomerular filtration has the potential to provide a mechanistic theory, based on experimental data, about the filtration characteristics of the glomerular filtration barrier. It provides a unique approach to the microanatomy of the glomerulus, renal autoregulation, and the pathogenesis of proteinuria.

Nonlinear inversion schemes for fluorescence optical tomography.

Fluorescence optical tomography is a non-invasive imaging modality that employs the absorption and re-emission of light by fluorescent dyes. The aim is to reconstruct the fluorophore distribution in a body from measurements of light intensities at the boundary. Due to the diffusive nature of light propagation in tissue, fluorescence tomography is a nonlinear and severely ill-posed problem, and some sort of regularization is required for a stable solution. In this paper we investigate reconstruction methods based on Tikhonov regularization with nonlinear penalty terms, namely total-variation regularization and a levelset-type method using a nonlinear parameterization of the unknown function. Moreover, we use the full threedimensional nonlinear forward model, which arises from the governing system of partial differential equations. We discuss the numerical realization of the regularization schemes by Newtontype iterations, present some details of the discretization by finite element methods, and outline the efficient implementation of sensitivity systems via adjoint methods. As we will demonstrate in numerical tests, the proposed nonlinear methods provide better reconstructions than standard methods based on linearized forward models and linear penalty terms. We will additionally illustrate, that the careful discretization of the methods derived on the continuous level allows to obtain reliable, mesh independent reconstruction algorithms.