Method for direct measurement of the film tension of black foam films at various capillary pressures

A new experimental method to directly measure the film tension of black foam films is developed on the basis of the Laplace equation. The method allows the determination of the tension of curved (spherical) films with various radii and capillary pressures. Measurements with Newton black films from sodium dodecyl sulphate aqueous solution have been carried out. The results show that in the studied range of curvature radii (70 ÷ 360 m) the film tension does not depend on the curvature and the capillary pressure.     

Advancements in subsea riser analysis using quasi-rotations and the Newton–Raphson method

Beam structures undergoing finite deflections and rotations in space have extensive application in the subsea industry particularly for the analysis of holistic systems with larger numbers of mooring and riser components. In using the finite element analysis approach, there is an increasing requirement for large element sizes which preserve accuracy with regard to the coupling of axial, bending and torsion response.The authors outline a method for improving the current state of practice for the analysis of riser systems. The approach draws on the convected coordinates method, Euler–Bernoulli beam theory, the principle of virtual work and the finite element method. Two quasi-rotation measures are developed including a quasi-material rotation definition for rotational deformation relative to the convected axis of a beam and a quasi-space rotation definition to deal with the path dependent nature of rotations in three dimensions.The novel aspect of this work is to relate the rate of change of the quasi-material rotation vector along the beam axis to a linear transformation of the beam axis rate-of-rotation vector through utilising the convected coordinates axes system. In this way, incremental values of quasi-material rotation are directly linked to incremental values of nodal quasi-space rotation and a global Newton–Raphson solution technique for interconnecting beam elements is straightforward to assemble.Furthermore, this leads to accurate definitions of coupled axial, bending and torque response for beams with significant deflection. The approach has particular advantages in the analysis of subsea riser sections. Also, the accuracy of the solution is preserved for a fewer number of elements compared to alternative solutions for computationally sensitive load cases with highly non-linear loading regimes.     

Parallel computation of two-point boundary-value problems via particular solutions

Nonlinear two-point boundary-value problems (TPBVP) can be reduced to the iterative solution of a sequence of linear problems by means of quasilinearization techniques. Therefore, the efficient solution of linear problems is the key to the efficient solution of nonlinear problems.Among the techniques available for solving linear two-point boundary-value problems, the method of particular solutions (MPS) is particularly attractive in that it employs only one differential system, the original nonhomogeneous system, albeit with different initial conditions. This feature of MPS makes it ideally suitable for implementation on parallel computers in that the following requirements are met: the computational effort is subdivided into separate tasks (particular solutions) assigned to the different processors; the tasks have nearly the same size; there is little intercommunication between the tasks.For the TPBVP, the speedup achievable is ofO(n), wheren is the dimension of the state vector, hence relatively modest for the differential systems of interest in trajectory optimization and guidance. This being the case, we transform the TPBVP into a multi-point boundary-value problem (MPBVP) involvingm time subintervals, withm–1 continuity conditions imposed at the interface of contiguous subintervals. For the MPBVP, the speedup achievable is ofO(mn), hence substantially higher than that achievable for the TPBVP. It reduces toO(m) if the parallelism is implemented only in the time domain and not in the state domain.A drawback of the multi-point approach is that it requires the solution of a large linear algebraic system for the constants of the particular solutions. This drawback can be offset by exploiting the particular nature of the interface conditions: if the vector of constants for the first subinterval is known, the vector of constants for the subsequent subintervals can be obtained with linear transformations. Using decomposition techniques together with the discrete version of MPS, the size of the linear algebraic system for the multi-point case becomes the same as that for the two-point case.Numerical tests on the Intel iPSC/860 computer show that substantial speedup can be achieved via parallel algorithms vis-a-vis sequential algorithms. Therefore, the present technique has considerable interest for real-time trajectory optimization and guidance.Dedicated to the Memory of Professor Jan M. SkowronskiThis paper, based on Refs. 1–3, is a much condensed version of the material contained in these references.The technical assistance of the Research Center on Parallel Computation of Rice University, Houston, Texas is gratefully acknowledged.