A Physically Consistent Implicit Viscosity Solver for SPH Fluids

In this paper, we present a novel physically consistent implicit solver for the simulation of highly viscous fluids using the Smoothed Particle Hydrodynamics (SPH) formalism. Our method is the result of a theoretical and practical in‐depth analysis of the most recent implicit SPH solvers for viscous materials. Based on our findings, we developed a list of requirements that are vital to produce a realistic motion of a viscous fluid. These essential requirements include momentum conservation, a physically meaningful behavior under temporal and spatial refinement, the absence of ghost forces induced by spurious viscosities and the ability to reproduce complex physical effects that can be observed in nature. On the basis of several theoretical analyses, quantitative academic comparisons and complex visual experiments we show that none of the recent approaches is able to satisfy all requirements. In contrast, our proposed method meets all demands and therefore produces realistic animations in highly complex scenarios. We demonstrate that our solver outperforms former approaches in terms of physical accuracy and memory consumption while it is comparable in terms of computational performance. In addition to the implicit viscosity solver, we present a method to simulate melting objects. Therefore, we generalize the viscosity model to a spatially varying viscosity field and provide an SPH discretization of the heat equation

Turbulent Micropolar SPH Fluids with Foam

In this paper we introduce a novel micropolar material model for the simulation of turbulent inviscid fluids. The governing equations are solved by using the concept of Smoothed Particle Hydrodynamics (SPH). SPH fluid simulations suffer from numerical diffusion which leads to a lower vorticity, a loss in turbulent details and finally in less realistic results. To solve this problem we propose a micropolar fluid model. The micropolar fluid model is a generalization of the classical Navier-Stokes equations, which are typically used in computer graphics to simulate fluids. In contrast to the classical Navier-Stokes model, micropolar fluids have a microstructure and therefore consider the rotational motion of fluid particles. In addition to the linear velocity field these fluids have a field of microrotation which represents existing vortices and provides a source for new ones. Our novel micropolar model can generate realistic turbulences, is linear and angular momentum conserving, can be easily integrated in existing SPH simulation methods and its computational overhead is negligible. Another important visual feature of turbulent liquids is foam. Therefore, we present a post-processing method which considers microrotation in the foam generation. It works completely automatic and requires only one user-defined parameter to control the amount of foam

Response of Sunflower (Helianthus annuus L.) Leaf Surface Defenses to Exogenous Methyl Jasmonate

Helianthus annuus, the common sunflower, produces a complex array of secondary compounds that are secreted into glandular trichomes, specialized structures found on leaf surfaces and anther appendages of flowers. The primary components of these trichome secretions are sesquiterpene lactones (STL), a diverse class of compounds produced abundantly by the plant family Compositae and believed to contribute to plant defense against herbivory. We treated wild and cultivated H. annuus accessions with exogenous methyl jasmonate, a plant hormone that mediates plant defense against insect herbivores and certain classes of fungal pathogens. The wild sunflower produced a higher density of glandular trichomes on its leaves than the cultivar. Comparison of the profiles of glandular trichome extracts obtained by liquid chromatography–mass spectroscopy (LC-MS) showed that wild and cultivated H. annuus were qualitatively similar in surface chemistry, although differing in the relative size and proportion of various compounds detected. Despite observing consistent transcriptional responses to methyl jasmonate treatment, we detected no significant effect on glandular trichome density or LC-MS profile in cultivated or wild sunflower, with wild sunflower exhibiting a declining trend in overall STL production and foliar glandular trichome density of jasmonate-treated plants. These results suggest that glandular trichomes and associated compounds may act as constitutive defenses or require greater levels of stimulus for induction than the observed transcriptional responses to exogenous jasmonate. Reduced defense investment in domesticated lines is consistent with predicted tradeoffs caused by selection for increased yield; future research will focus on the development of genetic resources to explicitly test the ecological roles of glandular trichomes and associated effects on plant growth and fitness

Renal transport kinetics of furosemide in the isolated perfused rat kidney

Direct quantitative data and corresponding theory are provided for the effect of protein binding on the renal transport of furosemide. Drug studies were performed with various combinations of bovine serum albumin and dextran. This resulted in a percent unbound ( fu ) of furosemide ranging from 0.785 to 85.8%. The corrected renal ( CLr/GFR ) and secretion ( CLs/GFR ) clearances of furosemide were observed to increase with percent free, but in a nonproportional manner. Plots of CLr/GFR or CLs/GFR vs. fu appeared to have a prominent y intercept as well as a convex ascending curve. In addition, the excretion ratio [ ER=CLr/ (fu · GFR) ] was reduced from 60.8 to 8.72 as fu increased. Overall, the data were best fitted to a model in which two Michaelis-Menten terms wre used to describe renal tubular transport, and secretion was dependent upon free drug concentrations in the perfusate. The results demonstrate that the renal mechanisms of furosemide excretion are more complex than previously reported and that active secretion may involve two different transport systems over the concentration range studied.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/45036/1/10928_2005_Article_BF01065259.pd

Density maps for improved SPH boundary handling

In this paper, we present the novel concept of density maps for robust handling of static and rigid dynamic boundaries in fluid simulations based on Smoothed Particle Hydrodynamics (SPH). In contrast to the vast majority of existing approaches, we use an implicit discretization for a continuous extension of the density field throughout solid boundaries. Using the novel representation we enhance accuracy and efficiency of density and density gradient evaluations in boundary regions by computationally efficient lookups into our density maps. The map is generated in a preprocessing step and discretizes the density contribution in the boundary's near-field. In consequence of the high regularity of the continuous boundary density field, we use cubic Lagrange polynomials on a narrow-band structure of a regular grid for discretization. This strategy not only removes the necessity to sample boundary surfaces with particles but also decouples the particle size from the number of sample points required to represent the boundary. Moreover, it solves the ever-present problem of particle deficiencies near the boundary. In several comparisons we show that the representation is more accurate than particle samplings, especially for smooth curved boundaries. We further demonstrate that our approach robustly handles scenarios with highly complex boundaries and even outperforms one of the most recent sampling based techniques

Fast Corotated FEM using Operator Splitting

In this paper we present a novel operator splitting approach for corotated FEM simulations. The deformation energy of the corotated linear material model consists of two additive terms. The first term models stretching in the individual spatial directions and the second term describes resistance to volume changes. By formulating the backward Euler time integration scheme as an optimization problem, we show that the first term is invariant to rotations. This allows us to use an operator splitting approach and to solve both terms individually with different numerical methods. The stretching part is solved accurately with an optimization integrator, which can be done very efficiently because the system matrix is constant over time such that its Cholesky factorization can be precomputed. The volume term is solved approximately by using the compliant constraints method and Gauss‐Seidel iterations. Further, we introduce the analytic polar decomposition which allows us to speed up the extraction of the rotational part of the deformation gradient and to recover inverted elements. Finally, this results in an extremely fast and robust simulation method with high visual quality that outperforms standard corotated FEMs by more than two orders of magnitude and even the fast but inaccurate PBD and shape matching methods by more than one order of magnitude without having their typical drawbacks. This enables a very efficient simulation of complex scenes containing more than a million elements