High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Robust handling of non‐linear constraints for GA optimization of aerodynamic shapes

A new approach to the robust handling of non‐linear constraints for GAs (genetic algorithms) optimization is proposed. A specific feature of the approach consists of the change in the conventional search strategy by employing search paths which pass through both feasible and infeasible points (contrary to the traditional approach where only feasible points may be included in a path). The method (driven by full Navier–Stokes computations) was applied to the problem of multiobjective optimization of aerodynamic shapes subject to various geometrical and aerodynamic constraints. The results demonstrated that the method retains high robustness of conventional GAs while keeping CFD computational volume to an acceptable level, which allowed the algorithm to be used in a demanding engineering environment. Copyright © 2004 John Wiley & Sons, Ltd.     

Dynamic surface tension measurements of surfactant solutions using the maximum bubble pressure method – limits of applicability

One of the essential differences in the design of bubble pressure tensiometers consists in the geometry of the measuring capillaries. To reach extremely short adsorption times of milliseconds and below, the so-called deadtime of the capillaries must be of the order of some 10 ms. In particular, for concentrated surfactant solutions, such as micellar solutions, short deadtimes are needed to minimize the initial surfactant load of the generated bubbles. A theoretical model is derived and confirmed by experiments performed for a wide range of experimental conditions, mainly in respect to variations in deadtime and bubble volume.     

On the influence of the surfactant''s polar group on the local and terminal velocities of bubbles

The local and the terminal velocities, the size and the degree of bubbles’ shape deformations were determined as a function of distance from the position of the bubble formation (capillary orifice) in solutions of n-octyltrimethylammonium bromide, n-octyldimethylphosphine oxide, n-octyl-β-D-glucopyranoside and n-octanoic acid.

These surface-active compounds have different polar groups but an identical hydrocarbon chain (C8) in the molecule. The motion of the bubbles was monitored and recorded using a stroboscopic illumination, a CCD camera, and a JVC professional video. The recorded bubble images were analyzed by the image analysis software. The bubbles accelerated rapidly and their shape was deformed immediately after detachment from the capillary. The extent of the bubbles’ shape deformation (ratio of horizontal and vertical diameters) was 1.5 in distilled water and dropped rapidly down to a level of ca. 1.05–1.03 with increasing surfactant concentration. After the acceleration period the bubbles either attained a constant value of the terminal velocity (distilled water and high concentrations of the solutions), or a maximum in the velocity profiles was observed (low concentrations). The values of the terminal velocity diminished drastically with increasing concentration, from the value of 35 cm/s in water down to about 15 cm/s, while the bubble diameter decreased by ca. 10% only. The surfactant adsorption at the surface of the bubbles was evaluated and the minimum adsorption coverages required to immobilize the bubbles’ surface were determined. It was found that this minimum adsorption coverage was ca. 4% for n-octyldimethylphosphine oxide, n-octyl-β-D-glucopyranoside, n-octanoic acid and 25% for n-octyltrimethylammonium bromide. The difference in the adsorption coverage together with the surfactants’ surface activities indicate that it is mainly the adsorption kinetics of the surfactants that governs the fluidity of interfaces of the rising bubbles.     

Elliptic and parabolic problems in unbounded domains

We consider elliptic and parabolic problems in unbounded domains. We give general existence and regularity results in Besov spaces and semi‐explicit representation formulas via operator‐valued fundamental solutions which turn out to be a powerful tool to derive a series of qualitative results about the solutions. We give a sample of possible applications including asymptotic behavior in the large, singular perturbations, exact boundary conditions on artificial boundaries and validity of maximum principles. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Remarks on the Regularity to 3-D Ideal Magnetohydrodynamic Equations

In this paper we are interested in the sufficient conditions which guarantee the regularityof solutions of 3-D ideal magnetohydrodynamic equations in the arbitrary time interval [0,T].Fivesufficient conditions are given.Our results are motivated by two main ideas:one is to control theaccumulation of vorticity alone;the other is to generalize the corresponding geometric conditions of3-D Euler equations to 3-D ideal magnetohydrodynamic equations.     

Delta and singular delta locus for one‐dimensional systems of conservation laws

This work gives a condition for existence of singular and delta shock wave solutions to Riemann problem for 2×2 systems of conservation laws. For a fixed left‐hand side value of Riemann data, the condition obtained in the paper describes a set of possible right‐hand side values. The procedure is similar to the standard one of finding the Hugoniot locus. Fluxes of the considered systems are globally Lipschitz with respect to one of the dependent variables. The association in a Colombeau‐type algebra is used as a solution concept. Copyright © 2004 John Wiley &Sons, Ltd.     

Some heteroclinic solutions of a model of skin pattern formation

In this paper we study travelling wave solutions to a system of four non‐linear partial differential equations, which arise in a tissue interaction model for skin morphogenesis. Under the ‘small‐stress’ assumption we prove the existence and uniqueness (up to a translation) of solutions with the dermis and epidermis cell densities being positive, which are a perturbation of a uniform epidermal cell density. We discuss the problem of the minimal wave‐speed. Copyright © 2004 John Wiley & Sons, Ltd.     

A monotone iterative scheme for a nonlinear second order equation based on a generalized anti–maximum principle

In this paper, a generalized anti–maximum principle for the second order differential operator with potentials is proved. As an application, we will give a monotone iterative scheme for periodic solutions of nonlinear second order equations. Such a scheme involves the L^p norms of the growth, 1 ≤ p ≤ ∞, while the usual one is just the case p = ∞.