894 research outputs found
Uniformly high order accurate essentially non-oscillatory schemes 3
In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear
Dewetting-controlled binding of ligands to hydrophobic pockets
We report on a combined atomistic molecular dynamics simulation and implicit
solvent analysis of a generic hydrophobic pocket-ligand (host-guest) system.
The approaching ligand induces complex wetting/dewetting transitions in the
weakly solvated pocket. The transitions lead to bimodal solvent fluctuations
which govern magnitude and range of the pocket-ligand attraction. A recently
developed implicit water model, based on the minimization of a geometric
functional, captures the sensitive aqueous interface response to the
concave-convex pocket-ligand configuration semi-quantitatively
Application of the level-set method to the implicit solvation of nonpolar molecules
A level-set method is developed for numerically capturing the equilibrium
solute-solvent interface that is defined by the recently proposed variational
implicit solvent model (Dzubiella, Swanson, and McCammon, Phys. Rev. Lett. {\bf
104}, 527 (2006) and J. Chem.\Phys. {\bf 124}, 084905 (2006)). In the level-set
method, a possible solute-solvent interface is represented by the zero
level-set (i.e., the zero level surface) of a level-set function and is
eventually evolved into the equilibrium solute-solvent interface. The evolution
law is determined by minimization of a solvation free energy {\it functional}
that couples both the interfacial energy and the van der Waals type
solute-solvent interaction energy. The surface evolution is thus an energy
minimizing process, and the equilibrium solute-solvent interface is an output
of this process. The method is implemented and applied to the solvation of
nonpolar molecules such as two xenon atoms, two parallel paraffin plates,
helical alkane chains, and a single fullerene . The level-set solutions
show good agreement for the solvation energies when compared to available
molecular dynamics simulations. In particular, the method captures solvent
dewetting (nanobubble formation) and quantitatively describes the interaction
in the strongly hydrophobic plate system
Justifying and Explaining Disproportionality 1968—2008: A Critique of Underlying Views of Culture
This is the publisher's version, also found here: http://cec.metapress.com/content/e7835w1374x7g141/?p=0ee3fcac0688484aacb1d3c83c9c469c&pi=2Special education has made considerable advances in research, policy, and practice in its short history. However, students from historically underserved groups continue to be disproportionately identified as requiring special education. Support for color-blind practices and policies can justify racial disproportionality in special education and signal a retrenchment to deficit views about students from historically underserved groups. We respond to these emerging concerns through an analysis of arguments that justify disproportionality. We also identify explanations of the problem and critique the views of culture that underlie these explanations. We conclude with a brief discussion of implications and future directions
Some flows in shape optimization
Geometric flows related to shape optimization problems of Bernoulli type are
investigated. The evolution law is the sum of a curvature term and a nonlocal
term of Hele-Shaw type. We introduce generalized set solutions, the definition
of which is widely inspired by viscosity solutions. The main result is an
inclusion preservation principle for generalized solutions. As a consequence,
we obtain existence, uniqueness and stability of solutions. Asymptotic behavior
for the flow is discussed: we prove that the solutions converge to a
generalized Bernoulli exterior free boundary problem
A Replica Inference Approach to Unsupervised Multi-Scale Image Segmentation
We apply a replica inference based Potts model method to unsupervised image
segmentation on multiple scales. This approach was inspired by the statistical
mechanics problem of "community detection" and its phase diagram. Specifically,
the problem is cast as identifying tightly bound clusters ("communities" or
"solutes") against a background or "solvent". Within our multiresolution
approach, we compute information theory based correlations among multiple
solutions ("replicas") of the same graph over a range of resolutions.
Significant multiresolution structures are identified by replica correlations
as manifest in information theory overlaps. With the aid of these correlations
as well as thermodynamic measures, the phase diagram of the corresponding Potts
model is analyzed both at zero and finite temperatures. Optimal parameters
corresponding to a sensible unsupervised segmentation correspond to the "easy
phase" of the Potts model. Our algorithm is fast and shown to be at least as
accurate as the best algorithms to date and to be especially suited to the
detection of camouflaged images.Comment: 26 pages, 22 figure
Globally optimal 3D image reconstruction and segmentation via energy minimisation techniques
This paper provides an overview of a number of techniques developed within our group to perform 3D reconstruction and image segmentation based of the application of energy minimisation concepts. We begin with classical snake techniques and show how similar energy minimisation concepts can be extended to derive globally optimal segmentation methods. Then we discuss more recent work based on geodesic active contours that can lead to globally optimal segmentations and reconstructions in 2D. Finally we extend the work to 3D by introducing continuous flow globally minimal surfaces. Several applications are discussed to show the wide applicability and suitability of these techniques to several difficult image analysis problems
A new ghost cell/level set method for moving boundary problems:application to tumor growth
In this paper, we present a ghost cell/level set method for the evolution of interfaces whose normal velocity depend upon the solutions of linear and nonlinear quasi-steady reaction-diffusion equations with curvature-dependent boundary conditions. Our technique includes a ghost cell method that accurately discretizes normal derivative jump boundary conditions without smearing jumps in the tangential derivative; a new iterative method for solving linear and nonlinear quasi-steady reaction-diffusion equations; an adaptive discretization to compute the curvature and normal vectors; and a new discrete approximation to the Heaviside function. We present numerical examples that demonstrate better than 1.5-order convergence for problems where traditional ghost cell methods either fail to converge or attain at best sub-linear accuracy. We apply our techniques to a model of tumor growth in complex, heterogeneous tissues that consists of a nonlinear nutrient equation and a pressure equation with geometry-dependent jump boundary conditions. We simulate the growth of glioblastoma (an aggressive brain tumor) into a large, 1 cm square of brain tissue that includes heterogeneous nutrient delivery and varied biomechanical characteristics (white matter, gray matter, cerebrospinal fluid, and bone), and we observe growth morphologies that are highly dependent upon the variations of the tissue characteristics—an effect observed in real tumor growth
- …