Ice formation on a smooth or rough cold surface due to the impact of a supercooled water droplet

Ice accretion is considered in the impact of a supercooled water droplet on a smooth or rough solid surface, the roughness accounting for earlier icing. In this theoretical investigation the emphasis and novelty lie in the full nonlinear interplay of the droplet motion and the growth of the ice surface being addressed for relatively small times, over a realistic range of Reynolds numbers, Froude numbers, Weber numbers, Stefan numbers and capillary underheating parameters. The Prandtl number and the kinetic under-heating parameter are taken to be order unity. The ice accretion brings inner layers into play forcibly, affecting the outer flow. (The work includes viscous effects in an isothermal impact without phase change, as a special case, and the differences between impact with and without freezing.) There are four main findings. First, the icing dynamically can accelerate or decelerate the spreading of the droplet whereas roughness on its own tends to decelerate spreading. The interaction between the two and the implications for successive freezings are found to be subtle. Second, a focus on the dominant physical effects reveals a multi-structure within which restricted regions of turbulence are implied. The third main finding is an essentially parabolic shape for a single droplet freezing under certain conditions. Fourth is a connection with a body of experimental and engineering work and with practical findings to the extent that the explicit predictions here for ice-accretion rates are found to agree with the experimental range.

Multiple spacecraft observations of interplanetary shocks: Characteristics of the upstream ULF turbulence

All interplanetary shocks observed by ISEE-3 and either ISEE-1 or ISEE-2 or both in 1978 and 1979 are examined for evidence of upstream waves. In order to characterize the properties of these shocks it is necessary to determine accurate shock normals. An overdetermined set of equations were inverted to obtain shock normals, velocities and error estimates for all these shocks. Tests of the method indicate it is quite reliable. Using these normals the Mach number and angle were between the interplanetary magnetic field and the shock normal for each shock. The upstream waves were separated into two classes: whistler mode precursors which occur at low Mach numbers and upstream turbulence whose amplitude at Mach numbers greater than 1.5 is controlled by the angle of the field to the shock normal. The former waves are right hand circularly polarized and quite monochromatic. The latter waves are more linearly polarized and have a broadband featureless spectrum

Low-dimensional models for turbulent plane Couette flow in a minimal flow unit

We model turbulent plane Couette flow in the minimal flow unit (MFU) – a domain whose spanwise and streamwise extent is just sufficient to maintain turbulence – by expanding the velocity field as a sum of optimal modes calculated via proper orthogonal decomposition from numerical data. Ordinary differential equations are obtained by Galerkin projection of the Navier–Stokes equations onto these modes. We first consider a 6-mode (11-dimensional) model and study the effects of including losses to neglected modes. Ignoring these, the model reproduces turbulent statistics acceptably, but fails to reproduce dynamics; including them, we find a stable periodic orbit that captures the regeneration cycle dynamics and agrees well with direct numerical simulations. However, restriction to as few as six modes artificially constrains the relative magnitudes of streamwise vortices and streaks and so cannot reproduce stability of the laminar state or properly account for bifurcations to turbulence as Reynolds number increases. To address this issue, we develop a second class of models based on ‘uncoupled’ eigenfunctions that allow independence among streamwise and cross-stream velocity components. A 9-mode (31-dimensional) model produces bifurcation diagrams for steady and periodic states in qualitative agreement with numerical Navier–Stokes solutions, while preserving the regeneration cycle dynamics. Together, the models provide empirical evidence that the ‘backbone’ for MFU turbulence is a periodic orbit, and support the roll–streak–breakdown–roll reformation picture of shear-driven turbulence

A functional analysis of change propagation

A thorough understanding of change propagation is fundamental to effective change management during product redesign. A new model of change propagation, as a result of the interaction of form and function is presented and used to develop an analysis method that determines how change is likely to propagate. The analysis produces a Design Structure Matrix, which clearly illustrates change propagation paths and highlights connections that could otherwise be ignored. This provides the user with an in-depth knowledge of product connectivity, which has the potential to support the design process and reduce the product's susceptibility to future change

Some calculations on the ground and lowest-triplet state of the helium isoelectronic sequence with the nucleus in motion

The method described in the preceding paper for the solution of two-electron atoms, which was used to calculate the 1 1S and 2 3S states of helium and heliumlike atoms within the fixed-nucleus approximation, has been applied to the case where all three particles are in relative motion. The solutions in the present case automatically include the effects of the mass polarization term and are compared with the results obtained for the term by using first-order perturbation theory with the fixed-nucleus wave functions. The input data for a particular atom consist of the atomic number, as before, but now the corresponding mass of the nucleus must be given also. Nonrelativistic energies with the nuclear mass included in the calculation have been obtained for the 1 1S and 2 3S states for Z ranging from 1 to 10. The energy with the nucleus in motion can be expressed only to eight significant figures (SF's) given the accuracy with which the relevant physical constants are known at present. All the results given here are computed as if these constants were known to ten SF's so that errors not incurred due to rounding. Convergence of the energies to ten SF's for both the singlet and triplet state was reached with a matrix of size 444 for Z values from 2 to 10. Convergence for the H- ion was a little slower

Effects of liquid slosh on rendezvous dynamics

Mathematical model using Euler-Lagrange equation to investigate liquid slosh effects on rendezvous dynamic

Sharp quadrature error bounds for the nearest-neighbor discretization of the regularized stokeslet boundary integral equation

The method of regularized stokeslets is a powerful numerical method to solve the Stokes flow equations for problems in biological fluid mechanics. A recent variation of this method incorporates a nearest-neighbor discretization to improve accuracy and efficiency while maintaining the ease-of-implementation of the original meshless method. This method contains three sources of numerical error, the regularization error associated from using the regularized form of the boundary integral equations (with parameter $\varepsilon$), and two sources of discretization error associated with the force and quadrature discretizations (with lengthscales $h_f$ and $h_q$). A key issue to address is the quadrature error: initial work has not fully explained observed numerical convergence phenomena. In the present manuscript we construct sharp quadrature error bounds for the nearest-neighbor discretisation, noting that the error for a single evaluation of the kernel depends on the smallest distance ($\delta$) between these discretization sets. The quadrature error bounds are described for two cases: with disjoint sets ($\delta>0$) being close to linear in $h_q$ and insensitive to $\varepsilon$, and contained sets ($\delta=0$) being quadratic in $h_q$ with inverse dependence on $\varepsilon$. The practical implications of these error bounds are discussed with reference to the condition number of the matrix system for the nearest-neighbor method, with the analysis revealing that the condition number is insensitive to $\varepsilon$ for disjoint sets, and grows linearly with $\varepsilon$ for contained sets. Error bounds for the general case ($\delta\geq 0$) are revealed to be proportional to the sum of the errors for each case.Comment: 12 pages, 6 figure