The dynamics of bubble clouds induced by high-intensity focused ultrasound are investigated in a regime where the cloud size is similar to the ultrasound wavelength. High-speed images show that the cloud is asymmetrical; the bubbles nearest the source grow to a larger radius than the distal ones. Similar structures of bubble clouds are observed in numerical simulations that mimic the laboratory experiment. To elucidate the structure, a parametric study is conducted for plane ultrasound waves with various amplitudes and diffuse clouds with different initial void fractions. Based on an analysis of the kinetic energy of liquid induced by bubble oscillations, a new scaling parameter is introduced to characterize the dynamics. The new parameter generalizes the cloud interaction parameter originally introduced by d'Agostino and Brennen (1989). The dynamic interaction parameter controls the energy localization and consequent anisoptropy of the cloud. Moreover, the amplitude of the far-field, bubble-scattered acoustics is likewise correlated with the proposed parameter. Findings of the present study not only shed light on the physics of cloud cavitation, but may also be of use to quantification of the effects of cavitation on outcomes of ultrasound therapies including HIFU-based lithotripsy.

A cellular automaton that is a generalization of the box-ball system with either many kinds of balls or finite carrier capacity is proposed and studied through two discrete integrable systems: nonautonomous discrete KP lattice and nonautonomous discrete two-dimensional Toda lattice. Applying reduction technique and ultradiscretization procedure to these discrete systems, we derive two types of time evolution equations of the proposed cellular automaton, and particular solutions to the ultradiscrete equations.

A nonautonomous version of the ultradiscrete hungry Toda lattice with a finite lattice boundary condition is derived by applying reduction and ultradiscretization to a nonautonomous two-dimensional discrete Toda lattice. It is shown that the derived ultradiscrete system has a direct connection to the box-ball system with many kinds of balls and finite carrier capacity. Particular solutions to the ultradiscrete system are constructed by using the theory of some sort of discrete biorthogonal polynomials.

We report on two superoutbursts of the AM CVn-type object CR Boo in 2014 April--March and 2015 May--June. A precursor outburst acompanied both of these superoutbursts. During the rising branch of the main superoutburst in 2014, we detected growing superhumps (stage A superhumps) whose period was $0.017669(24)$ d. Assuming that this period reflects the dynamical precession rate at the radius of the 3:1 resonance, we could estimate the mass ratio ($q=M_2/M_1$) of 0.101(4) by using the stage A superhump period and the orbital one of 0.0170290(6) d. This mass ratio is consistent with that expected by the theoretical evolutionary model of AM CVn-type objects. The detection of precursor outbursts and stage A superhumps is the second case in AM CVn-type objects. There are two interpretations of the outbursts of AM CVn-type objects. One is a dwarf nova (DN) outbursts analogy, which is caused by thermal and tidal instabilities. Another is the VY Scl-type variation, which is caused by the variation of the mass-transfer rate of the secondary. This detection of the superhump variations strongly suggests the former interpretation.

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 102 SU UMa-type dwarf novae observed mainly during the 2014-2015 season and characterized these objects. Our project has greatly improved the statistics of the distribution of orbital periods, which is a good approximation of the distribution of cataclysmic variables at the terminal evolutionary stage, and confirmed the presence of a period minimum at a period of 0.053 d and a period spike just above this period. The number density monotonically decreased toward the longer period and there was no strong indication of a period gap. We detected possible negative superhumps in Z Cha. It is possible that normal outbursts are also suppressed by the presence of a disk tilt in this system. There was no indication of enhanced orbital humps just preceding the superoutburst, and this result favors the thermal-tidal disk instability as the origin of superoutbursts. We detected superhumps in three AM CVn-type dwarf novae. Our observations and recent other detections suggest that 8% of objects showing dwarf nova-type outbursts are AM CVn-type objects. AM CVn-type objects and EI Psc-type object may be more abundant than previously recognized. OT J213806, a WZ Sge-type object, exhibited a remarkably different feature between the 2010 and 2014 superoutbursts. Although the 2014 superoutburst was much fainter the plateau phase was shorter than the 2010 one, the course of the rebrightening phase was similar. This object indicates that the O-C diagrams of superhumps can be indeed variable at least in WZ Sge-type objects. Four deeply eclipsing SU UMa-type dwarf novae (ASASSN-13cx, ASASSN-14ag, ASASSN-15bu, NSV 4618) were identified. We studied long-term trends in supercycles in MM Hya and CY UMa and found systematic variations of supercycles of ~20%.

We report on a superoutburst of the AM CVn-type object SDSS J090221.35+381941.9 [J0902; orbital period 0.03355(6) d] in 2014 March-April. The entire outburst consisted of a precursor outburst and the main superoutburst, followed by a short rebrightening. During the rising branch of the main superoutburst, we detected growing superhumps (stage A superhumps) with a period of 0.03409(1) d. During the plateau phase of the superoutburst, superhumps with a shorter period (stage B superhumps) were observed. Using the orbital period and the period of the stage A superhumps, we were able to measure the dynamical precession rate of the accretion disk at the 3:1 resonance, and obtained a mass ratio (q) of 0.041(7). This is the first successful measurement of the mass ratio in an AM CVn-type object using the recently developed stage A superhump method. The value is generally in good agreement with the theoretical evolutionary model. The orbital period of J0902 is the longest among the outbursting AM CVn-type objects, and the borderline between the outbursting systems and systems with stable cool disks appears to be longer than had been supposed.

The spectral transformation technique for symmetric R_{II} polynomials is developed. Use of this technique reveals that the nonautonomous discrete modified KdV (nd-mKdV) lattice is directly connected with the R_{II} chain. Hankel determinant solutions to the semi-infinite nd-mKdV lattice are also presented.

A generalized eigenvalue algorithm for tridiagonal matrix pencils is presented. The algorithm appears as the time evolution equation of a nonautonomous discrete integrable system associated with a polynomial sequence which has some orthogonality on the support set of the zeros of the characteristic polynomial for a tridiagonal matrix pencil. The convergence of the algorithm is discussed by using the solution to the initial value problem for the corresponding discrete integrable system.

A connection between the finite ultradiscrete Toda lattice and the box-ball system is extended to the case where each box has own capacity and a carrier has a capacity parameter depending on time. In order to consider this connection, new carrier rules "size limit for solitons" and "recovery of balls", and a concept "expansion map" are introduced. A particular solution to the extended system of a special case is also presented.