skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A Gaussian-like immersed-boundary kernel with three continuous derivatives and improved translational invariance

; ;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 316; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-10-06 15:26:41; Journal ID: ISSN 0021-9991
Country of Publication:
United States

Citation Formats

Bao, Yuanxun, Kaye, Jason, and Peskin, Charles S. A Gaussian-like immersed-boundary kernel with three continuous derivatives and improved translational invariance. United States: N. p., 2016. Web. doi:10.1016/
Bao, Yuanxun, Kaye, Jason, & Peskin, Charles S. A Gaussian-like immersed-boundary kernel with three continuous derivatives and improved translational invariance. United States. doi:10.1016/
Bao, Yuanxun, Kaye, Jason, and Peskin, Charles S. Fri . "A Gaussian-like immersed-boundary kernel with three continuous derivatives and improved translational invariance". United States. doi:10.1016/
title = {A Gaussian-like immersed-boundary kernel with three continuous derivatives and improved translational invariance},
author = {Bao, Yuanxun and Kaye, Jason and Peskin, Charles S.},
abstractNote = {},
doi = {10.1016/},
journal = {Journal of Computational Physics},
number = C,
volume = 316,
place = {United States},
year = {Fri Jul 01 00:00:00 EDT 2016},
month = {Fri Jul 01 00:00:00 EDT 2016}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/

Save / Share:
  • Highlights: {yields} A radial basis function based smoothing technique is introduced. {yields} It is more general and easier to implement compared to other techniques. {yields} With this technique, a combined immersed boundary and spectral method is developed. {yields} It is shown that the proposed method works better in terms of reducing the non-physical Gibbs oscillation. - Abstract: It is known that, when the immersed boundary method (IBM) is implemented within spectral-like methods, the Gibbs oscillation seriously deteriorates the calculation of derivatives near the body surface. In this paper, a radial basis function (RBF) based smoothing technique is proposed with themore » intention of eliminating or efficiently reducing the Gibbs oscillation without affecting the flow field outside the body. Based on this technique, a combined IBM/spectral scheme is developed to solve the incompressible Navier-Stokes equations. Numerical simulations of flow through a periodic lattice of cylinders of various cross sections are performed. The results demonstrate that the proposed methodology is able to give accurate and nearly oscillation-free numerical solutions of incompressible viscous flows.« less
  • We derive the necessary boundary conditions required to ensure unique solutions of a new noncompact-kernel integral equation for three-particle scattering. We give an abstract analysis to determine all solutions of this equation. Configuration space boundary conditions are derived to single out the physical solution.
  • We present a complete analysis of the multipartite entanglement of three-mode Gaussian states of continuous-variable systems. We derive standard forms which characterize the covariance matrix of pure and mixed three-mode Gaussian states up to local unitary operations, showing that the local entropies of pure Gaussian states are bound to fulfill a relationship which is stricter than the general Araki-Lieb inequality. Quantum correlations can be quantified by a proper convex roof extension of the squared logarithmic negativity, the continuous-variable tangle, or contangle. We review and elucidate in detail the proof that in multimode Gaussian states the contangle satisfies a monogamy inequalitymore » constraint [G. Adesso and F. Illuminati, New J. Phys8, 15 (2006)]. The residual contangle, emerging from the monogamy inequality, is an entanglement monotone under Gaussian local operations and classical communications and defines a measure of genuine tripartite entanglements. We determine the analytical expression of the residual contangle for arbitrary pure three-mode Gaussian states and study in detail the distribution of quantum correlations in such states. This analysis yields that pure, symmetric states allow for a promiscuous entanglement sharing, having both maximum tripartite entanglement and maximum couplewise entanglement between any pair of modes. We thus name these states GHZ/W states of continuous-variable systems because they are simultaneous continuous-variable counterparts of both the GHZ and the W states of three qubits. We finally consider the effect of decoherence on three-mode Gaussian states, studying the decay of the residual contangle. The GHZ/W states are shown to be maximally robust against losses and thermal noise.« less
  • Dielectrophoresis, a nonlinear electrokinetic transport mechanism, has become popular in many engineering applications including manipulation, characterization and actuation of biomaterials, particles and biological cells. In this paper, we present a hybrid immersed interface–immersed boundary method to study AC dielectrophoresis where an algorithm is developed to solve the complex Poisson equation using a real variable formulation. An immersed interface method is employed to obtain the AC electric field in a fluid media with suspended particles and an immersed boundary method is used for the fluid equations and particle transport. The convergence of the proposed algorithm as well as validation of themore » hybrid scheme with experimental results is presented. In this paper, the Maxwell stress tensor is used to calculate the dielectrophoretic force acting on particles by considering the physical effect of particles in the computational domain. Thus, this study eliminates the approximations used in point dipole methods for calculating dielectrophoretic force. A comparative study between Maxwell stress tensor and point dipole methods for computing dielectrophoretic forces are presented. The hybrid method is used to investigate the physics of dielectrophoresis in microfluidic devices using an AC electric field. The numerical results show that with proper design and appropriate selection of applied potential and frequency, global electric field minima can be obtained to facilitate multiple particle trapping by exploiting the mechanism of negative dielectrophoresis. Our numerical results also show that electrically neutral particles form a chain parallel to the applied electric field irrespective of their initial orientation when an AC electric field is applied. This proposed hybrid numerical scheme will help to better understand dielectrophoresis and to design and optimize microfluidic devices.« less
  • Purpose: To present the validity of our beam modeling with double and triple Gaussian dose kernels for spot scanning proton beams in Nagoya Proton Therapy Center. This study investigates the conformance between the measurements and calculation results in absolute dose with two types of beam kernel. Methods: A dose kernel is one of the important input data required for the treatment planning software. The dose kernel is the 3D dose distribution of an infinitesimal pencil beam of protons in water and consists of integral depth doses and lateral distributions. We have adopted double and triple Gaussian model as lateral distributionmore » in order to take account of the large angle scattering due to nuclear reaction by fitting simulated inwater lateral dose profile for needle proton beam at various depths. The fitted parameters were interpolated as a function of depth in water and were stored as a separate look-up table for the each beam energy. The process of beam modeling is based on the method of MDACC [X.R.Zhu 2013]. Results: From the comparison results between the absolute doses calculated by double Gaussian model and those measured at the center of SOBP, the difference is increased up to 3.5% in the high-energy region because the large angle scattering due to nuclear reaction is not sufficiently considered at intermediate depths in the double Gaussian model. In case of employing triple Gaussian dose kernels, the measured absolute dose at the center of SOBP agrees with calculation within ±1% regardless of the SOBP width and maximum range. Conclusion: We have demonstrated the beam modeling results of dose distribution employing double and triple Gaussian dose kernel. Treatment planning system with the triple Gaussian dose kernel has been successfully verified and applied to the patient treatment with a spot scanning technique in Nagoya Proton Therapy Center.« less