
In this paper, we consider numerical approximations of a binary
fluidsurfactant phasefield model coupled with the fluid flow, in which the
system is highly nonlinear that couples the incompressible NavierStokes
equations and two CahnHilliard 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 NavierStokes equation, and a subtle
implicitexplicit treatment for the stress and convective terms. We prove the
wellposedness 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 topk 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 anyk
ranking algorithm: for a given time budget, re turn as many of the topranked
results as possible. Then, given additional time, produce the next lowerranked
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 topranked 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 nontrivial 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
CahnHilliard 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
semiexplicitly, 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 smecticA liquid crystals. The model, derived from the variational
approach of the modified OseenFrank energy, is a highly nonlinear system that
couples the incompressible NavierStokes equations and a constitutive equation
for the layer variable. We develop two linear, secondorder timemarching
schemes based on the "Invariant Energy Quadratization" method for nonlinear
terms in the constitutive equation, the projection method for the NavierStokes
equations, and some subtle implicitexplicit treatments for the convective and
stress terms. Moreover, we prove the wellposedness 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
CahnHilliard type phase equation, we adopt the AllenCahn type phase field
model with the static contact line boundary condition that coupled with
incompressible NavierStokes equations with Navier boundary condition. The
projection method is used to deal with the NavierStokes equa tions and an
auxiliary function is introduced for the nonconvex GinzburgLandau 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 NavierStokes equations with a generalized Navier boundary
condition, and the CahnHilliard 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 NavierStokes equations, the invariant energy quadratization for the
nonlinear gradient terms in the bulk and boundary, and a subtle
implicitexplicit treatment for the stress and convective terms. The
wellposedness of the semidiscretized system and the unconditional energy
stabilities are proved. Various numerical results based on a spectralGalerkin
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 ICARUST600
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 fewGeV 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
CahnHilliard equation with concentration dependent mobility, and the
logarithmic FloryHuggins 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 timediscrete
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
semiexplicitly. 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
fluidsurfactant phase field model, in which the free energy contains a
nonlinear coupling entropy, a GinzburgLandau double well potential, and a
logarithmic FloryHuggins potential. The resulting system consists of two
coupled, nonlinear CahnHilliard 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 semiexplicitly. 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
CahnHilliard phasefield 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
semiexplicitly. Consequently, the resulting numerical schemes lead to a
wellposed 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 fluidsurfactant phase field model that consists two nonlinearly coupled
CahnHilliard equations. The main challenge in solving the system numerically
is how to develop easytoimplement 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 timestepping schemes using the
socalled "Invariant Energy Quadratization" approach for the double well
potentials and a subtle explicitimplicit technique for the nonlinear coupling
potential. Moreover, the resulting linear system is wellposed 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
smecticA 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 NavierStokes equations, and two nonlinear
coupled secondorder elliptic equations. Based on some subtle
explicitimplicit 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 GinzburgLandau 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 NavierStokes equation, and a subtle implicitexplicit 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 twophase hydrodynamics coupled
phasefield model that incorporates the variable densities, viscosities and
moving contact line boundary conditions. The model is a nonlinear, coupled
system that consists of incompressible NavierStokes equations with the
generalized Navier boundary condition, and the CahnHilliard equations with
moving contact line boundary conditions. By some subtle explicitimplicit
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 spectralGalerkin 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 CahnHilliard 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 timediscrete 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 semiexplicitly. 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.

Cyclictocycle variability, CCV, of intakejet flow in an optical engine was
measured using particle image velocimetry (PIV), revealing the possibility of
two different flow patterns. A phasedependent 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 subsampled within a single test;
each data set contained individual cycles showing one flow pattern or the
other. Threedimensional velocity data from a largeeddy 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 cycletocycle variations
in the swirl ratio result in the two different flow patterns in the PIV plane.
Having an unambiguous metric to reveal largescale flow CCV, causes for this
variability were examined within the possible sources present in the available
testing. In particular, variations in intakeport and cylinder pressure,
lateral valve oscillations, and engine RPM were examined as potential causes
for the cycletocycle flow ariations using the phasedependent POD
coefficients. No direct correlation was seen between the intake port pressure,
or the pressure drop across the intake valve, and the incylinder flow pattern.
A correlation was observed between dominant flow pattern and cycletocycle
variations in intake valve horizontal position. RPM values and incylinder flow
patterns did not correlate directly. However, a shift in flow pattern was
observed between early and late cycles in a 2900cycle 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, pseudocaustic 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 $\alpha1$ 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
BigBOSSlike survey and the improvement can reach a factor of $\sim 2$ for
surveys with much denser background populations.