Incompressible SPH simulation of wave interaction with porous structure

In this paper an incompressible Smoothed Particle Hydrodynamics (ISPH) method is applied to investigate the flow motion in and around the porous structure. In order to describe in a simple and effective way the flow through the interface between the porous region and pure fluid region within the SPH framework, a heuristic boundary treatment method has been proposed. The ISPH model is first verified against a theoretical model of wave propagation over a porous bed and then further validated by comparing the predicted wave surface profiles and flow velocity fields with the experiment data for a typical case of flow motion around and inside a submerged porous structure. The good agreement has demonstrated that the improved ISPH model developed in this work is capable of modelling wave interaction with porous structures

lassopack: Model selection and prediction with regularized regression in Stata

This article introduces lassopack, a suite of programs for regularized regression in Stata. lassopack implements lasso, square-root lasso, elastic net, ridge regression, adaptive lasso and post-estimation OLS. The methods are suitable for the high-dimensional setting where the number of predictors $p$ may be large and possibly greater than the number of observations, $n$. We offer three different approaches for selecting the penalization (`tuning') parameters: information criteria (implemented in lasso2), $K$-fold cross-validation and $h$-step ahead rolling cross-validation for cross-section, panel and time-series data (cvlasso), and theory-driven (`rigorous') penalization for the lasso and square-root lasso for cross-section and panel data (rlasso). We discuss the theoretical framework and practical considerations for each approach. We also present Monte Carlo results to compare the performance of the penalization approaches.Comment: 52 pages, 6 figures, 6 tables; submitted to Stata Journal; for more information see https://statalasso.github.io

Massive charged particle's tunneling from spherical charged black hole

We generalize the Parikh-Wilczek scheme to the tunneling of a massive charged particle from a general spherical charged black hole. We obtain that the tunneling probability depends on the energy, the mass and the charge of the particle. In particular, the modified Hawking temperature is related to the charge. Only at the leading order approximation can the standard Hawking temperature be reproduced. We take the Reissner-Nordstr\"{o}m black hole as an example to clarify our points of view, and find that the accumulation of Hawking radiation makes it approach an extreme black hole.Comment: 10 pages, no figures; v2: a minor typo corrected; v3: 11 pages, clarification and reference added, final version to be published in EPL; v4: minor modifications to match the published versio

Markov bases and subbases for bounded contingency tables

In this paper we study the computation of Markov bases for contingency tables whose cell entries have an upper bound. In general a Markov basis for unbounded contingency table under a certain model differs from a Markov basis for bounded tables. Rapallo, (2007) applied Lawrence lifting to compute a Markov basis for contingency tables whose cell entries are bounded. However, in the process, one has to compute the universal Gr\"obner basis of the ideal associated with the design matrix for a model which is, in general, larger than any reduced Gr\"obner basis. Thus, this is also infeasible in small- and medium-sized problems. In this paper we focus on bounded two-way contingency tables under independence model and show that if these bounds on cells are positive, i.e., they are not structural zeros, the set of basic moves of all $2 \times 2$ minors connects all tables with given margins. We end this paper with an open problem that if we know the given margins are positive, we want to find the necessary and sufficient condition on the set of structural zeros so that the set of basic moves of all $2 \times 2$ minors connects all incomplete contingency tables with given margins.Comment: 22 pages. It will appear in the Annals of the Institution of Statistical Mathematic

A Centimeter-Sized Quaternary Ti-Zr-Be-Ag Bulk Metallic Glass

A novel centimeter-sized Ti-based bulk metallic glass (BMG) was developed by the addition of Ag in the ternary Ti41Zr25Be34 glassy alloy. By replacing Be with Ag, the glass forming ability (GFA), the yield strength, and the supercooled liquid temperature of the quaternary Ti41Zr25Be34−xAgx (x=2, 4, 6, 8 at.%) glassy alloys have been obviously enhanced. Among the developed Ti-Zr-Be-Ag alloy systems, the Ti41Zr25Be28Ag6 alloy possesses the largest critical diameter (Dmax) of 10 mm, while the yield strength is also enhanced to 1961 MPa, which is much larger than that of Ti41Zr25Be34 (1755 MPa) alloy. The experimental results show that Ag is an effective element for improving the GFA and the yield strength of Ti-Zr-Be glassy alloy

Generalized tt-jj Model

By parameterizing the t-j model we present a new electron correlation model with one free parameter for high-temperature superconductivity. This model is of $SU_{q}(1,2)$ symmetry. The energy spectrums are shown to be modulated by the free parameter in the model. The solution and symmetric structures of the Hilbert space, as well as the Bethe ansatz approach are discussed for special cases.Comment: 13 page, Latex, to appear in J. Phys.

A Generalized Representation Formula for Systems of Tensor Wave Equations

In this paper, we generalize the Kirchhoff-Sobolev parametrix of Klainerman and Rodnianski to systems of tensor wave equations with additional first-order terms. We also present a different derivation, which better highlights that such representation formulas are supported entirely on past null cones. This generalization is a key component for extending Klainerman and Rodnianski's breakdown criterion result for Einstein-vacuum spacetimes to Einstein-Maxwell and Einstein-Yang-Mills spacetimes.Comment: 29 page