The effect of an insoluble surfactant on the skin friction of a bubble

In this paper, we develop a novel moving mesh method suitable for solving axisymmetric free-boundary problems, including the Marangoni effect induced by surfactant or temperature variation. This method employs a body-fitted grid system where the gas-liquid interface is one line of the grid system. We model the surfactant equation of state with a non-linear Langmuir law, and, for simplicity, we limit ourselves to the situation of an insoluble surfactant. We solve complicated dynamic boundary conditions accurately on the gas-liquid interface in the framework of finite-volume methods. Our method is used to study the effect of a surfactant on the skin friction of a bubble in a uniaxial flow. For the limiting case where the surface diffusivity is zero, the effect of a tangential stress generated by the surface tension gradient, allows us to explain a new phenomenon in high concentration regimes: larger surface tension, but also larger deformation. Furthermore, this condition leads to the formation of boundary layers and flow separation at high Reynolds numbers. The influence of these complex flow patterns is examined. © 2005 Elsevier SAS. All rights reserved.

Influence of insoluble surfactant on the deformation and breakup of a bubble or thread in a viscous fluid

The influence of surfactant on the breakup of a prestretched bubble in a quiescent viscous surrounding is studied by a combination of direct numerical simulation and the solution of a long-wave asymptotic model. The direct numerical simulations describe the evolution toward breakup of an inviscid bubble, while the effects of small but non-zero interior viscosity are readily included in the long-wave model for a fluid thread in the Stokes flow limit. The direct numerical simulations use a specific but realizable and representative initial bubble shape to compare the evolution toward breakup of a clean or surfactant-free bubble and a bubble that is coated with insoluble surfactant. A distinguishing feature of the evolution in the presence of surfactant is the interruption of bubble breakup by formation of a slender quasi-steady thread of the interior fluid. This forms because the decrease in surface area causes a decrease in the surface tension and capillary pressure, until at a small but non-zero radius, equilibrium occurs between the capillary pressure and interior fluid pressure. The long-wave asymptotic model, for a thread with periodic boundary conditions, explains the principal mechanism of the slender thread's formation and confirms, for example, the relatively minor role played by the Marangoni stress. The large-time evolution of the slender thread and the precise location of its breakup are, however, influenced by effects such as the Marangoni stress and surface diffusion of surfactant. © 2008 Cambridge University Press.

Mode I crack growth resistance of metallic foams

A Dugdale-type cohesive zone model is used to predict the mode I crack growth resistance (R-curve) of metallic foams, with the fracture process characterized by an idealized traction-separation law that relates the crack surface traction to crack opening displacement. A quadratic yield function, involving the von Mises effective stress and mean stress, is used to account for the plastic compressibility of metallic foams. Finite element calculations are performed for the crack growth resistance under small scale yielding and small scale bridging in plane strain, with K-field boundary conditions. The following effects upon the fracture process are quantified: material hardening, bridging strength, T-stress (the non-singular stress acting parallel to the crack plane), and the shape of yield surface. To study the failure behaviour and notch sensitivity of metallic foams in the presence of large scale yielding, a study is made for panels embedded with either a centre-crack or an open hole and subjected to tensile stressing. For the centre-cracked panel, a transition crack size is predicted for which the fracture response switches from net section yielding to elastic-brittle fracture. Likewise, for a panel containing a centre-hole, a transition hole diameter exists for which the fracture response switches from net section yielding to a local maximum stress criterion at the edge of the hole.

Analysis and simulation of combustion oscillations in a laboratory spray combustor

The unstable combustion that can occur in combustion chambers is a major problem for aeroengines and ground-based industrial gas turbines. Nowadays, CFD provides a flexible, low cost tool to supplement direct measurement. This paper presents simulations of combustion oscillations in a liquid-fuelled experimental rig at the University of Cambridge. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone and to develop inlet and outlet boundary conditions just upstream and downstream of the combustion region enabling the CFD calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in good agreement with experimental measurements. More details about the unstable combustion can be obtained from the simulation results. The approach developed here is expected to provide a powerful tool for the design and operation of stable combustion systems. Copyright © 2009 by ASME.

Microstructural optimization of cellular acoustic materials

A microstructure based acoustic model is introduced, which can be used to optimize the microstructure of cellular materials and thus to obtain their optimal acoustic property. This acoustic model is an unsteady one which is appropriate in the limit of low Reynolds numbers. The model involves three elements. This first involves the propagation of acoustic waves passing the cylinders whose axes are aligned parallel to the direction of propagation. The second model relates to the propagation of acoustic waves passing the cylinders whose axes are aligned perpendicular to the direction of propagation. In both cases the interaction between adjacent cylinders is taken into account by considering the effect of polygonal periodic boundary conditions. As these two models are linear they are combined to give the characteristics of propagation at arbitrary incidence. The third model involves propagation passing spheres in order to represent the joints. Heat transfer is also included. These three models are then used to expand the design space and calculate the optimum cell structure for desired acoustic performance in a number of different applications. Moreover, the application fields are also analyzed.

Investigation of the losses in the stator duct spacers of large A.C. motors using the A -ip 3-D eddy current formulation

This paper describes how the A -if) formulation may be applied to determine the losses in the stator duct spacers of large a.c. motors. The model is described in terms of its geometry and boundary conditions. The novel aspects of the application of the formulation to this problem are explained. These include the modelling of fixed currents sources (the stator windings), the location of the necessary cut surfaces and the determination of their magnetic scalar potential differences, and the implementation of periodic boundary conditions for vector variables. Results are presented showing how the duct spacer losses vary with load, and with the relative permeability of the spacer material. The effects of modelling iron nonlinearity, of both the spacer and the steel laminations, are also illustrated. © 1996 IEEE.

On the domain decomposition and transmission line modelling finite element method for time-domain induction motor analysis

This paper demonstrates how a finite element model which exploits domain decomposition is applied to the analysis of three-phase induction motors. It is shown that a significant gain in cpu time results when compared with standard finite element analysis. Aspects of the application of the method which are particular to induction motors are considered: the means of improving the convergence of the nonlinear finite element equations; the choice of symmetrical sub-domains; the modelling of relative movement; and the inclusion of periodic boundary conditions. © 1999 IEEE.

A discrete dislocation analysis of fatigue crack growth

Cyclic loading of a plane strain mode I crack under small scale yielding is analyzed using discrete dislocation dynamics. The dislocations are all of edge character, and are modeled as line singularities in an elastic solid. At each stage of loading, superposition is used to represent the solution in terms of solutions for edge dislocations in a half-space and a non-singular complementary solution that enforces the boundary conditions, which is obtained from a linear elastic, finite element solution. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated into the formulation through a set of constitutive rules. An irreversible relation between the opening traction and the displacement jump across a cohesive surface ahead of the initial crack tip is also specified, which permits crack growth to emerge naturally. It is found that crack growth can occur under cyclic loading conditions even when the peak stress intensity factor is smaller than the stress intensity required for crack growth under monotonic loading conditions; however below a certain threshold value of ΔKI no crack growth was seen.

Secondary flow near an undulatory surface induced by wall blocking effect

Prandtl's secondary mean motions of the second kind near an undulating surface were explained in terms of turbulent blocking effect and kinematic boundary conditions at the surface, and its order of magnitude was estimated. Isotropic turbulence is distorted by the undulating surface of wavelength λ and amplitude h with a low slope, so that h « λ. The prime mechanism for generating the mean flow is that the far-field Isotropic turbulence is distorted by the non-local blocking effect of the surface to become anisotropic axisymmetric turbulence near the surface with principal axis that is not aligned with the local curvature of the undulation. Then the local analysis can be applied and the mechanism is similar to the mean flow generation mechanism for homogeneous axisymmetric turbulence over a planer surface, i.e. gradients of the Reynolds stress caused by the turbulent blocking effect generate the mean motions. The results from this simple analysis are consistent with previous exact analysis in which the effects of curvature are strictly taken into account. The results also qualitatively agree with flow visualization over an undulating surface in a mixing-box.

Adjoint linearised Euler solver in the frequency domain for jet noise modelling

The paper is devoted to extending the new efficient frequency-domain method of adjoint Green's function calculation to curvilinear multi-block RANS domains for middle and farfield sound computations. Numerical details of the method such as grids, boundary conditions and convergence acceleration are discussed. Two acoustic source models are considered in conjunction with the method and acoustic modelling results are presented for a benchmark low-Reynolds-number jet case.

TOPS (Task optimization in the presence of signal-dependent noise) model

This paper proposes a movement trajectory planning model, which is a maximum task achievement model in which signal-dependent noise is added to the movement command. In the proposed model, two optimization criteria are combined, maximum task achievement and minimum energy consumption. The proposed model has the feature that the end-point boundary conditions for position, velocity, and acceleration need not be prespecified. Consequently, the method can be applied not only to the simple point-to-point movement, but to any task. In the method in this paper, the hand trajectory is derived by a psychophysical experiment and a numerical experiment for the case in which the target is not stationary, but is a moving region. It is shown that the trajectory predicted from the minimum jerk model or the minimum torque change model differs considerably from the results of the psychophysical experiment. But the trajectory predicted from the maximum task achievement model shows good qualitative agreement with the hand trajectory obtained from the psychophysical experiment. © 2004 Wiley Periodicals, Inc.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Vibration of fluid loaded conical shells.

An analytical model is presented to describe the vibration of a truncated conical shell with fluid loading in the low frequency range. The solution for the dynamic response of the shell is presented in the form of a power series. Fluid loading is taken into account by dividing the shell into narrow strips which are considered to be locally cylindrical. Analytical results are presented for different boundary conditions and have been compared with the computational results from a boundary element model. Limitations of the model to the low frequency range are discussed.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Synchronization of modulated traveling baroclinic waves in a periodically forced, rotating fluid annulus.

Frequency entrainment and nonlinear synchronization are commonly observed between simple oscillatory systems, but their occurrence and behavior in continuum fluid systems are much less well understood. Motivated by possible applications to geophysical fluid systems, such as in atmospheric circulation and climate dynamics, we have carried out an experimental study of the interaction of fully developed baroclinic instability in a differentially heated, rotating fluid annulus with an externally imposed periodic modulation of the thermal boundary conditions. In quasiperiodic and chaotic amplitude-modulated traveling wave regimes, the results demonstrate a strong interaction between the natural periodic modulation of the wave amplitude and the externally imposed forcing. This leads to partial or complete phase synchronization. Synchronization effects were observed even with very weak amplitudes of forcing, and were found with both 1:1 and 1:2 frequency ratios between forcing and natural oscillations.