On the theory of surface waves in water generated by moving disturbances

The wave profile generated by an obstacle moving at constant veiocity U over a water surface of infinite extent appears to be stationary with respect to the moving body provided, of course, the motion has been maintained for a long time. When the gravitational and capillary effects are both taken into account, the surface waves so generated may possess a minimum phase velocity c[sub]m characterized by a certain wave length, say [lambda][sub]m (see Ref. 1, p. 459). If the velocity U of the solid body is greater than c[sub]m, then the physically correct solution of this two-dimensional problem requires that the gravity waves (of wave length greater than [lambda][sub]m) should exist only on the downstream side and the capillary waves (of wave length less than [lambda][sub]m) only on the upstream side. If one follows strictly the so-called steady-state formulation so that the time does not appear in the problem, one finds in general that it is not possible to characterize uniquely the mathematical solution with the desired physical properties by imposing only the boundedness conditions at infinity. [Footnote: In the case of a three-dimensional steady-state problem, even the condition that the disturbance should vanish at infinity is not sufficient to characterize the unique solution.] Some stronger radiation conditions are actually necessary. In the linearized treatment of this stationary problem, several methods have been employed, most of which are aimed at obtaining the correct solution by introducing some artificial device, either of a mathematical or physical nature. One of these methods widely used was due to Rayleigh, and was further discussed by Lamb. In the analysis of this problem Rayleigh introduced a "small dissipative force", proportional to the velocity relative to the moving stream. This "law" of friction does not originate from viscosity and is hence physically fictitious, for in the final result this dissipation factor is made to vanish eventually. In the present investigation, Rayleigh's friction coefficient is shown to correspond roughly to a time convergence factor for obtaining the steady-state solution from an initial value problem. (It is not a space-limit factor for fixing the boundary conditions at space infinity, as has usually been assumed in explanation of its effect). Thus, the introduction of Rayleigh's coefficient is only a mathematical device to render the steady-state solution mathematically determinate and physically acceptable. For a physical understanding, however, it is confusing and even misleading; for example, in an unsteady flow case it leads to an incomplete solution, as has been shown by Green. Another approach, purely of a physical nature, was used by Michell in his treatment of the velocity potential for thin ships. To make the problem determinate, he chose the solution which represents the gravity waves propagating only downstream and discarded the part corresponding to the waves traveling upstream. For two-dimensional problems with the capillary effect, this method would mean a superposition of simple waves so as to make the solution physically correct. Some other methods appear to be limited in the necessity of interpreting the principal value of a certain kind of improper integral. In short, as to their physical soundness and mathematical rigor, or even to their merits or demerits, the preference of one method over the others has remained nevertheless a matter of considerable dispute. Only until recently the steady-state problem has been treated by first formulating a corresponding initial value problem. A brief historical sketch of these methods is given in the next section. The purpose of this paper is to try to understand the physical mechanism underlying the steady configuration of the surface wave phenomena and to clarify to a certain extent the background of the artifices adopted for solution of steady-state problems. The point of view to be presented here is that this problem should be formulated first as an initial value problem (for example, the body starts to move with constant velocity at a certain time instant), and then the stationary state is sought by passing to the limit as the time tends to infinity. If at any finite time instant the boundary condition that the disturbance vanishes at infinity (because of the finite wave velocity) is imposed, then the limiting solution as the time tends to infinity is determinate and bears automatically the desired physical properties. Also, from the integral representation of the linearized solution, the asymptotic behavior of the wave form for large time is derived in detail, showing the distribution of the wave trains in space. This asymptotic solution exhibits an interesting picture which reveals how the dispersion* generates two monochromatic wave trains, with the capillary wave in front of, and the gravity wave behind, the surface pressure. *[Footnote: By dispersive medium is meant one in which the wave velocity of a propagating wave depends on the wave length, so that a number of wave trains of different wave lengths tends to form groups, propagating with group velocities which are in general different from the phase velocities of individual wave trains. In case of waves on the water surface, both the gravity and surface tension are responsible for dispersion.] The special cases U< c[sub]m and U = c[sub]m are also discussed. The viscous effect and the effect of superposition are commented upon later. Through this detailed investigation it is found that the dispersive effect, not the viscous effect plays the significant role in producing the final stationary wave configuration

Critical manifold of the kagome-lattice Potts model

Any two-dimensional infinite regular lattice G can be produced by tiling the plane with a finite subgraph B of G; we call B a basis of G. We introduce a two-parameter graph polynomial P_B(q,v) that depends on B and its embedding in G. The algebraic curve P_B(q,v) = 0 is shown to provide an approximation to the critical manifold of the q-state Potts model, with coupling v = exp(K)-1, defined on G. This curve predicts the phase diagram both in the ferromagnetic (v>0) and antiferromagnetic (v<0) regions. For larger bases B the approximations become increasingly accurate, and we conjecture that P_B(q,v) = 0 provides the exact critical manifold in the limit of infinite B. Furthermore, for some lattices G, or for the Ising model (q=2) on any G, P_B(q,v) factorises for any choice of B: the zero set of the recurrent factor then provides the exact critical manifold. In this sense, the computation of P_B(q,v) can be used to detect exact solvability of the Potts model on G. We illustrate the method for the square lattice, where the Potts model has been exactly solved, and the kagome lattice, where it has not. For the square lattice we correctly reproduce the known phase diagram, including the antiferromagnetic transition and the singularities in the Berker-Kadanoff phase. For the kagome lattice, taking the smallest basis with six edges we recover a well-known (but now refuted) conjecture of F.Y. Wu. Larger bases provide successive improvements on this formula, giving a natural extension of Wu's approach. The polynomial predictions are in excellent agreement with numerical computations. For v>0 the accuracy of the predicted critical coupling v_c is of the order 10^{-4} or 10^{-5} for the 6-edge basis, and improves to 10^{-6} or 10^{-7} for the largest basis studied (with 36 edges).Comment: 31 pages, 12 figure

Critical surfaces for general inhomogeneous bond percolation problems

We present a method of general applicability for finding exact or accurate approximations to bond percolation thresholds for a wide class of lattices. To every lattice we sytematically associate a polynomial, the root of which in $[0,1]$ is the conjectured critical point. The method makes the correct prediction for every exactly solved problem, and comparison with numerical results shows that it is very close, but not exact, for many others. We focus primarily on the Archimedean lattices, in which all vertices are equivalent, but this restriction is not crucial. Some results we find are kagome: $p_c=0.524430...$, $(3,12^2): p_c=0.740423...$, $(3^3,4^2): p_c=0.419615...$, $(3,4,6,4):p_c=0.524821...$, $(4,8^2):p_c=0.676835...$, $(3^2,4,3,4)$: $p_c=0.414120...$ . The results are generally within $10^{-5}$ of numerical estimates. For the inhomogeneous checkerboard and bowtie lattices, errors in the formulas (if they are not exact) are less than $10^{-6}$.Comment: Submitted to J. Stat. Mec

Predictive modeling of die filling of the pharmaceutical granules using the flexible neural tree

In this work, a computational intelligence (CI) technique named flexible neural tree (FNT) was developed to predict die filling performance of pharmaceutical granules and to identify significant die filling process variables. FNT resembles feedforward neural network, which creates a tree-like structure by using genetic programming. To improve accuracy, FNT parameters were optimized by using differential evolution algorithm. The performance of the FNT-based CI model was evaluated and compared with other CI techniques: multilayer perceptron, Gaussian process regression, and reduced error pruning tree. The accuracy of the CI model was evaluated experimentally using die filling as a case study. The die filling experiments were performed using a model shoe system and three different grades of microcrystalline cellulose (MCC) powders (MCC PH 101, MCC PH 102, and MCC DG). The feed powders were roll-compacted and milled into granules. The granules were then sieved into samples of various size classes. The mass of granules deposited into the die at different shoe speeds was measured. From these experiments, a dataset consisting true density, mean diameter (d50), granule size, and shoe speed as the inputs and the deposited mass as the output was generated. Cross-validation (CV) methods such as 10FCV and 5x2FCV were applied to develop and to validate the predictive models. It was found that the FNT-based CI model (for both CV methods) performed much better than other CI models. Additionally, it was observed that process variables such as the granule size and the shoe speed had a higher impact on the predictability than that of the powder property such as d50. Furthermore, validation of model prediction with experimental data showed that the die filling behavior of coarse granules could be better predicted than that of fine granules

High resolution, low temperature photoabsorption cross-section of C2H2 with application to Saturn's atmosphere

New laboratory observations of the VUV absorption cross-section of C2H2, obtained under physical conditions approximating stratospheres of the giant planets, were combined with IUE observations of the albedo of Saturn, for which improved data reduction techniques have been used, to produce new models for that atmosphere. When the effects of C2H2 absorption are accounted for, additional absorption by other molecules is required. The best-fitting model also includes absorption by PH3, H2O, C2H6 and CH4. A small residual disagreement near 1600 A suggests that an additional trace species may be required to complete the model

Anomalies in non-stoichiometric uranium dioxide induced by pseudo-phase transition of point defects

A uniform distribution of point defects in an otherwise perfect crystallographic structure usually describes a unique pseudo phase of that state of a non-stoichiometric material. With off-stoichiometric uranium dioxide as a prototype, we show that analogous to a conventional phase transition, these pseudo phases also will transform from one state into another via changing the predominant defect species when external conditions of pressure, temperature, or chemical composition are varied. This exotic transition is numerically observed along shock Hugoniots and isothermal compression curves in UO2 with first-principles calculations. At low temperatures, it leads to anomalies (or quasi-discontinuities) in thermodynamic properties and electronic structures. In particular, the anomaly is pronounced in both shock temperature and the specific heat at constant pressure. With increasing of the temperature, however, it transforms gradually to a smooth cross-over, and becomes less discernible. The underlying physical mechanism and characteristics of this type of transition are encoded in the Gibbs free energy, and are elucidated clearly by analyzing the correlation with the variation of defect populations as a function of pressure and temperature. The opportunities and challenges for a possible experimental observation of this phase change are also discussed.Comment: 11 pages, 5 figure

Modeling and Analysis of Power Processing Systems (MAPPS). Volume 2: Appendices

The computer programs and derivations generated in support of the modeling and design optimization program are presented. Programs for the buck regulator, boost regulator, and buck-boost regulator are described. The computer program for the design optimization calculations is presented. Constraints for the boost and buck-boost converter were derived. Derivations of state-space equations and transfer functions are presented. Computer lists for the converters are presented, and the input parameters justified

Modeling and Analysis of Power Processing Systems (MAPPS). Volume 1: Technical report

Computer aided design and analysis techniques were applied to power processing equipment. Topics covered include: (1) discrete time domain analysis of switching regulators for performance analysis; (2) design optimization of power converters using augmented Lagrangian penalty function technique; (3) investigation of current-injected multiloop controlled switching regulators; and (4) application of optimization for Navy VSTOL energy power system. The generation of the mathematical models and the development and application of computer aided design techniques to solve the different mathematical models are discussed. Recommendations are made for future work that would enhance the application of the computer aided design techniques for power processing systems