• In this paper, we consider numerical approximations of a binary fluid-surfactant phase-field model coupled with the fluid flow, in which the system is highly nonlinear that couples the incompressible Navier-Stokes equations and two Cahn-Hilliard type equations. We develop two, linear and second order time marching schemes for solving this system, by combining the "Invariant Energy Quadratization" approach for the nonlinear potentials, the projection method for the Navier-Stokes equation, and a subtle implicit-explicit treatment for the stress and convective terms. We prove the well-posedness of the linear system and its unconditional energy stability rigorously. Various 2D and 3D numerical experiments are performed to validate the accuracy and energy stability of the proposed schemes.
  • Many problems in areas as diverse as recommendation systems, social network analysis, semantic search, and distributed root cause analysis can be modeled as pattern search on labeled graphs (also called "heterogeneous information networks" or HINs). Given a large graph and a query pattern with node and edge label constraints, a fundamental challenge is to nd the top-k matches ac- cording to a ranking function over edge and node weights. For users, it is di cult to select value k . We therefore propose the novel notion of an any-k ranking algorithm: for a given time budget, re- turn as many of the top-ranked results as possible. Then, given additional time, produce the next lower-ranked results quickly as well. It can be stopped anytime, but may have to continues until all results are returned. This paper focuses on acyclic patterns over arbitrary labeled graphs. We are interested in practical algorithms that effectively exploit (1) properties of heterogeneous networks, in particular selective constraints on labels, and (2) that the users often explore only a fraction of the top-ranked results. Our solution, KARPET, carefully integrates aggressive pruning that leverages the acyclic nature of the query, and incremental guided search. It enables us to prove strong non-trivial time and space guarantees, which is generally considered very hard for this type of graph search problem. Through experimental studies we show that KARPET achieves running times in the order of milliseconds for tree patterns on large networks with millions of nodes and edges.
  • In this paper, we consider numerical approximations for the anisotropic Cahn-Hilliard equation. The main challenge of constructing numerical schemes with unconditional energy stabilities for this model is how to design proper temporal discretizations for the nonlinear terms with the strong anisotropy. We propose two, second order time marching schemes by combining the recently developed SAV approach with the linear stabilization approach, where three linear stabilization terms are added. These terms are shown to be crucial to remove the oscillations caused by the anisotropic coefficients, numerically. The novelty of the proposed schemes is that all nonlinear terms can be treated semi-explicitly, and one only needs to solve three decoupled linear equations with constant coefficients at each time step. We further prove the unconditional energy stabilities rigorously, and present various 2D and 3D numerical simulations to demonstrate the stability and accuracy.
  • In this paper, we consider the numerical approximations for a hydrodynamical model of smectic-A liquid crystals. The model, derived from the variational approach of the modified Oseen-Frank energy, is a highly nonlinear system that couples the incompressible Navier-Stokes equations and a constitutive equation for the layer variable. We develop two linear, second-order time-marching schemes based on the "Invariant Energy Quadratization" method for nonlinear terms in the constitutive equation, the projection method for the Navier-Stokes equations, and some subtle implicit-explicit treatments for the convective and stress terms. Moreover, we prove the well-posedness of the linear system and their unconditionally energy stabilities rigorously. Various numerical experiments are presented to demonstrate the stability and the accuracy of the numerical schemes in simulating the dynamics under shear flow and the magnetic field.
  • In this paper, we propose a novel, thermodynamically consistent phase field model to simulate the deformation and breakup of a ferrodroplet that is immersed in a viscous medium and subject to an applied uniform magnetic field. Instead of using the magnetic body force in the traditional Rosensweig model, the key idea of this model is to propose a new magnetic energy that enables direct effects of the magnetic field on the interface evolution. The model can thereby be derived from the variational principle via minimizing the free energy of the total system. This energy based modeling idea can be easily extended to include more external fields in a convenient and consistent way for more complicated applications involving multiple external fields. We validate the model by performing a series of numerical simulations, including the comparison with analytic solutions, the investigation of the effect of different types of magnetic fields, the dynamical behaviors of the ferrodroplet breakup under the strong magnetic field, magnetic/velocity/pressure field distributions, the magnetic energy density, and the inertial phenomenon, etc.
  • In this paper, we present an efficient energy stable scheme to solve a phase field model incorporating contact line condition. Instead of the usually used Cahn-Hilliard type phase equation, we adopt the Allen-Cahn type phase field model with the static contact line boundary condition that coupled with incompressible Navier-Stokes equations with Navier boundary condition. The projection method is used to deal with the Navier-Stokes equa- tions and an auxiliary function is introduced for the non-convex Ginzburg-Landau bulk potential. We show that the scheme is linear, decoupled and energy stable. Moreover, we prove that fully discrete scheme is also energy stable. An efficient finite element spatial discretization method is implemented to verify the accuracy and efficiency of proposed schemes. Numerical results show that the proposed scheme is very efficient and accurate
  • We consider the numerical approximations for a phase field model consisting of incompressible Navier--Stokes equations with a generalized Navier boundary condition, and the Cahn-Hilliard equation with a dynamic moving contact line boundary condition. A crucial and challenging issue for solving this model numerically is the time marching problem, due to the high order, nonlinear, and coupled properties of the system. We solve this issue by developing two linear, second order accurate, and energy stable schemes based on the projection method for the Navier--Stokes equations, the invariant energy quadratization for the nonlinear gradient terms in the bulk and boundary, and a subtle implicit-explicit treatment for the stress and convective terms. The well-posedness of the semidiscretized system and the unconditional energy stabilities are proved. Various numerical results based on a spectral-Galerkin spatial discretization are presented to verify the accuracy and efficiency of the proposed schemes.
  • The measurement of muon momentum by Multiple Coulomb Scattering is a crucial ingredient to the reconstruction of {\nu}{\mu} CC events in the ICARUS-T600 liquid argon TPC in absence of magnetic field, as in the search for sterile neutrinos at Fermilab where ICARUS will be exposed to ~1 GeV Booster neutrino beam. A sample of ~1000 stopping muons produced by charged current interactions of CNGS {\nu}{\mu} in the surrounding rock at the INFN Gran Sasso underground Laboratory provides an ideal benchmark in the few-GeV range since their momentum can be directly and independently obtained by the calorimetric measurement. Stopping muon momentum in the 0.5- 4.5 GeV/c range has been reconstructed via Multiple Coulomb Scattering with resolution ranging from 10 to 25 % depending on muon energy, track length and uniformity of the electric field in the drift volume.
  • In this paper, we consider the numerical approximations for the fourth order Cahn-Hilliard equation with concentration dependent mobility, and the logarithmic Flory-Huggins potential. One challenge in solving such a diffusive system numerically is how to develop proper temporal discretization for nonlinear terms in order to preserve the energy stability at the time-discrete level. We resolve this issue by developing a set of the first and second order time marching schemes based on a novel, called "Invariant Energy Quadratization" approach. Its novelty is that the proposed scheme is linear and symmetric positive definite because all nonlinear terms are treated semi-explicitly. We further prove all proposed schemes are unconditionally energy stable rigorously. Various 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy and efficiency of the proposed schemes thereafter.
  • In this paper, we consider the numerical solution of a binary fluid-surfactant phase field model, in which the free energy contains a nonlinear coupling entropy, a Ginzburg-Landau double well potential, and a logarithmic Flory-Huggins potential. The resulting system consists of two coupled, nonlinear Cahn-Hilliard type equations. We develop a set of first and second order time marching schemes for this system using the "Invariant Energy Quadratization" approach, in particular, the system is transformed into an equivalent one by introducing appropriate auxiliary variables and all nonlinear terms are then treated semi-explicitly. Both schemes are linear and lead to symmetric positive definite systems at each time step, thus they can be efficiently solved. We further prove that these schemes are unconditionally energy stable in the discrete sense. Various 2D and 3D numerical experiments are performed to validate the accuracy and energy stability of the proposed schemes.
  • How to develop efficient numerical schemes while preserving the energy stability at the discrete level is a challenging issue for the three component Cahn-Hilliard phase-field model. In this paper, we develop first and second order temporal approximation schemes based on the "Invariant Energy Quadratization" approach, where all nonlinear terms are treated semi-explicitly. Consequently, the resulting numerical schemes lead to a well-posed linear system with the symmetric positive definite operator to be solved at each time step. We rigorously prove that the proposed schemes are unconditionally energy stable. Various 2D and 3D numerical simulations are presented to demonstrate the stability and the accuracy of the schemes.
  • In this paper, we consider the numerical approximations for the commonly used binary fluid-surfactant phase field model that consists two nonlinearly coupled Cahn-Hilliard equations. The main challenge in solving the system numerically is how to develop easy-to-implement time stepping schemes while preserving the unconditional energy stability. We solve this issue by developing two linear and decoupled, first order and a second order time-stepping schemes using the so-called "Invariant Energy Quadratization" approach for the double well potentials and a subtle explicit-implicit technique for the nonlinear coupling potential. Moreover, the resulting linear system is well-posed and the linear operator is symmetric positive definite. We rigorously prove the first order scheme is unconditionally energy stable. Various numerical simulations are presented to demonstrate the stability and the accuracy thereafter.
  • In this paper, we consider numerical approximations for the model of smectic-A liquid crystal flows. The model equation, that is derived from the variational approach of the de Gennes free energy, is a highly nonlinear system that couples the incompressible Navier-Stokes equations, and two nonlinear coupled second-order elliptic equations. Based on some subtle explicit--implicit treatments for nonlinear terms, we develop a unconditionally energy stable, linear and decoupled time marching numerical scheme. We also rigorously prove that the proposed scheme obeys the energy dissipation law at the discrete level. Various numerical simulations are presented to demonstrate the accuracy and the stability thereafter.
  • In this paper, we consider numerical approximations of a hydrodynamically coupled phase field diblock copolymer model, in which the free energy contains a kinetic potential, a gradient entropy, a Ginzburg-Landau double well potential, and a long range nonlocal type potential. We develop a set of second order time marching schemes for this system using the "Invariant Energy Quadratization" approach for the double well potential, the projection method for the Navier-Stokes equation, and a subtle implicit-explicit treatment for the stress and convective term. The resulting schemes are linear and lead to symmetric positive definite systems at each time step, thus they can be efficiently solved. We further prove that these schemes are unconditionally energy stable. Various numerical experiments are performed to validate the accuracy and energy stability of the proposed schemes.
  • We consider the numerical approximations of a two-phase hydrodynamics coupled phase-field model that incorporates the variable densities, viscosities and moving contact line boundary conditions. The model is a nonlinear, coupled system that consists of incompressible Navier--Stokes equations with the generalized Navier boundary condition, and the Cahn--Hilliard equations with moving contact line boundary conditions. By some subtle explicit--implicit treatments to nonlinear terms, we develop two efficient, unconditionally energy stable numerical schemes, in particular, a linear decoupled energy stable scheme for the system with static contact line condition, and a nonlinear energy stable scheme for the system with dynamic contact line condition. An efficient spectral-Galerkin spatial discretization is implemented to verify the accuracy and efficiency of proposed schemes. Various numerical results show that the proposed schemes are efficient and accurate.
  • In this paper, we consider the numerical approximations for the fourth order viscous Cahn-Hilliard equation with the hyperbolic relaxation. The main challenge in solving such a diffusive system numerically is how to develop high order temporal discretization for the hyperbolic and nonlinear terms that allows large time step while preserving the unconditional energy stability, i.e., the energy dissipative structure at the time-discrete level. We resolve this issue by developing two second order time marching schemes using the recently developed "Invariant Energy Quadratization" approach where all nonlinear terms are treated semi-explicitly. In each time step, one only needs to solve a linear system that is symmetric positive definite. We rigorously prove all proposed schemes are unconditionally energy stable, which is further verified by time step refinement test numerically. Moreover, a number of 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy and efficiency of the proposed schemes.
  • Cyclic-to-cycle variability, CCV, of intake-jet flow in an optical engine was measured using particle image velocimetry (PIV), revealing the possibility of two different flow patterns. A phase-dependent proper orthogonal decomposition (POD) analysis showed that one or the other flow pattern would appear in the average flow, sampled from test to test or sub-sampled within a single test; each data set contained individual cycles showing one flow pattern or the other. Three-dimensional velocity data from a large-eddy simulation (LES) of the engine showed that the PIV plane was cutting through a region of high shear between the intake jet and another large flow structure. Rotating the measurement plane 10{\deg} revealed one or the other flow structure observed in the PIV measurements. Thus, it was hypothesized that cycle-to-cycle variations in the swirl ratio result in the two different flow patterns in the PIV plane. Having an unambiguous metric to reveal large-scale flow CCV, causes for this variability were examined within the possible sources present in the available testing. In particular, variations in intake-port and cylinder pressure, lateral valve oscillations, and engine RPM were examined as potential causes for the cycle-to-cycle flow ariations using the phase-dependent POD coefficients. No direct correlation was seen between the intake port pressure, or the pressure drop across the intake valve, and the in-cylinder flow pattern. A correlation was observed between dominant flow pattern and cycle-to-cycle variations in intake valve horizontal position. RPM values and in-cylinder flow patterns did not correlate directly. However, a shift in flow pattern was observed between early and late cycles in a 2900-cycle test after an approximately 5 rpm engine speed perturbation.
  • A cosmological preferred direction was reported from the type Ia supernovae (SNe Ia) data in recent years. We use the Union2.1 data to give a simple classification of such studies for the first time. Because the maximum anisotropic direction is independent of isotropic dark energy models, we adopt two cosmological models ($\Lambda$CDM, $w$CDM) for the hemisphere comparison analysis and $\Lambda$CDM model for dipole fit approach. In hemisphere comparison method, the matter density and the equation of state of dark energy are adopted as the diagnostic qualities in the $\Lambda$CDM model and $w$CDM model, respectively. In dipole fit approach, we fit the fluctuation of distance modulus. We find that there is a null signal for the hemisphere comparison method, while a preferred direction ($b=-14.3^\circ \pm 10.1^\circ, l=307.1^\circ \pm 16.2^\circ$) for the dipole fit method. This result indicates that the dipole fit is more sensitive than the hemisphere comparison method.
  • Using analytical method, we study the Singular Isothermal Quadrupole (SIQ) lens system, which is the simplest lens model that can produce four images. In this case, the radial mass distribution is in accord with the profile of the Singular Isothermal Sphere (SIS) lens, and the tangential distribution is given by adding a quadrupole on the monopole component. The basic properties of the SIQ lens have been studied in this paper, including deflection potential, deflection angle, magnification, critical curve, caustic, pseudo-caustic and transition locus. Analytical solutions of the image positions and magnifications for the source on axes are derived. As have been found, naked cusps will appear when the relative intensity $k$ of quadrupole to monopole is larger than 0.6. According to the magnification invariant theory of the SIQ lens, the sum of the signed magnifications of the four images should be equal to unity \citep{dal98}. However, if a source lies in the naked cusp, the summed magnification of the left three images is smaller than the invariant 1. With this simple lens system, we study the situations that a point source infinitely approaches a cusp or a fold. The sum of magnifications of cusp image triplet is usually not equal to 0, and it is usually positive for major cusp while negative for minor cusp. Similarly, the sum of magnifications of fold image pair is usually neither equal to 0. Nevertheless, the cusp and fold relations are still equal to 0, in that the sum values are divided by infinite absolute magnifications by definition.
  • Cosmic magnification has been detected through cross correlation between foreground and background populations (galaxies or quasars). It has been shown that weighing each background object by its $\alpha-1$ can significantly improve the cosmic magnification measurement \citep{Menard02,Scranton05}. Here, $\alpha$ is the logarithmic slope of the luminosity function of background populations. However, we find that this weighting function is optimal only for sparse background populations in which intrinsic clustering is negligible with respect to shot noise. We derive the optimal weighting function for general case including scale independent and scale dependent weights. The optimal weighting function improves the S/N (signal to noise ratio) by $\sim 20%$ for a BigBOSS-like survey and the improvement can reach a factor of $\sim 2$ for surveys with much denser background populations.