Mathematical modelling of the oxidation of nickel titanium alloys

A one-dimensional bulk reaction model for the oxidation of nickeltitanium is formulated, with preferential oxidation of titaniumbeing included. The modelling is directed at the better understandingof the dominant mechanisms involved in the oxidation processand their significance for the biocompatibility of the alloy.Two different regimes for the relative diffusivities of oxygenand the metals are investigated. By assuming fast bulk reactions,different asymptotic structures emerge in different parameterregimes and the resulting models take the form of moving boundaryproblems. Different profiles of nickel concentration are obtained:in particular a nickel-rich layer (observed in practice) ispresent below the oxide/metal interface for the case when oxygenand the metals diffuse at comparable rates.     

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 excitation of a closely spaced array by a line source

An infinite row of periodically spaced, identical rigid circularcylinders is excited by an acoustic line source which is parallelto the generators of the cylinders. A method for calculatingthe scattered field accurately and efficiently is presented.When the cylinders are sufficiently close together, Rayleigh–Blochsurface waves that propagate energy to infinity along the arrayare excited. An expression is derived which enables the amplitudesof these surface waves to be computed without requiring thesolution to the full scattering problem.     

Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical analysis for the periodic response of nonsmooth dynamic systems of high-dimension

This paper deals with the numerical solution of high dimensional dynamic systems with nonsmooth characteristics, such as structural elastoplasticity, mechanical clearance and dry friction. As a stable periodic response of the system attracts its transient response, it may be extracted from the transient response somehow. The paper suggests a novel iteration scheme to extrapolate the periodic response from a short time history of its transient response by curve fitting. Compared with the current schemes such as shooting and incremental harmonic balancing, the present scheme makes full use of the information within the transient response and the system characters, then its computation efficiency has increased by an order and the numerical convergence depending on the initial iteration has greatly improved.The subject supported by the Chinese Foundation of Aeronautical Science     

Letter: Some Remarks on Standing Gravitational Waves

An intuitive definition of standing gravitational waves is proposed. Some main classes of exact vacuum solutions are searched for standing gravitational waves, in most cases with a negative result. Only some Einstein–Rosen waves meet the conditions.     

Production of superfine fibers by the electrospinning method

Weakly conducting fluid flows in an electric field in the presence of phase interfaces are investigated theoretically and experimentally. The object of the investigation is an axisymmetric jet. The analysis is carried out within the framework of electrohydrodynamics (EHD) with allowance for surface charge transfer on mobile interfaces. The jet shape is calculated at large distances from the outflow point. The theoretical and experimental data are compared for Newtonian and polymeric fluids.     

Cellulose-based fibers from liquid–crystalline solutions. II. Processing and morphology of acetate/butyrate esters

Fibers were spun from isotropic and anisotropic dimethylacetamide solutions of cellulose esters. Take-up speeds of the dry jet/wet spinning process varied. Water served as the coagulant. The mechanical properties of the fibers increased as spinning progressed from the isotropic to the anisotropic state of the solution. A trade-off in solubility and fiber properties was noted as the butyryl acetyl ratio decreased. Whereas high butyryl content enhances both overall solubility and the formation of liquid–crystalline solutions at lower concentration, it results in lower fiber modulus and strength. Morphology of the fibers depended on the coagulation rate which was influenced by the concentration of the sppinning solution. The level of orientation and crystallinity of the fibers increased somewhat when they were spun from liquid-crystalline solutions. © 1993 John Wiley & Sons, Inc.