Theory of weak turbulence of waves on a fluid surface

The nonlinear interaction of waves in a fluid of finite depth is discussed. Forbidden decay processes in the gravitational portion of the spectrum are eliminated from the Hamiltonian by means of a canonical transformation. This provides an opportunity to obtain a kinetic equation which takes into account scattering of capillary waves by gravitational waves, in addition to decays in the subsystem of gravitational waves. The distribution N_k P^1/2h^1/4k^–4 is obtained for capillary waves in shallow water with constant flow of energy P with respect to the spectrum in the space of the wave numbers k. The interaction of the gravitational and capillary turbulence spectra is discussed. An induced distribution of gravitational waves is found which results from their interaction with capillary waves. It is an increasing function of the wave numbers q in the region bounded by the capillary constant k_o, N_q q^9/4 (q < k_o). The coupling of spectra in the gravitational and capillary regions and the conversion from slightly turbulent distributions to universal distributions are discussed.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 97–106, November–December, 1974.     

Study on stress wave propagation in fractured rocks with fractal joint surfaces

This paper presents the experimental and theoretical investigation of property of stress wave propagation in jointed rocks by means of SHPB technique and fractal geometry method. Our aim focuses on the influence of the rough joint surface configuration on stress wave propagation. The comparison of behavior of reflection and transmission waves, deformation and energy dissipation of a rough joint surface characterized by its fractal feature with that of a smooth plane joint has been carried out. It has shown that the rough joint surface distinctly affects the stress wave propagation and energy dissipation in the jointed rocks. The rougher the joint surface was, the more permanent deformation occurred and the more attenuation stress wave took place as well. A nonlinear relationship between the normalized energy dissipation ratio W_J/W_I of the jointed rock and the joint roughness in terms of the fractal dimension has been formulated. It seems that the ratio W_J/W_I, presenting how much energy has been dissipated in the joint, nonlinearly increased with the increment of the fractal dimension D of the jointed surface. The ratio W_J/W_I of a roughly jointed rock, however, tends to be the same as that of a smoothly jointed rock if the fractal dimension is less than a critical value D_c = 2.20. The energy dissipation ratio at the critical point D_c seem to be a constant, not dependent of rock type but fractal joint configuration.     

A modified Lin equation for the energy balance in isotropic turbulence

At sufficiently large Reynolds numbers, turbulence is expected to exhibit scale-invariance in an intermediate (“inertial”) range of wavenumbers, as shown by power law behavior of the energy spectrum and also by a constant rate of energy transfer through wavenumber. However, there is an apparent contradiction between the definition of the energy flux (i.e., the integral of the transfer spectrum) and the observed behavior of the transfer spectrum itself. This is because the transfer spectrum T(k) is invariably found to have a zero-crossing at a single point (at k = k_*), implying that the corresponding energy flux cannot have an extended plateau but must instead have a maximum value at k = k_*. This behavior was formulated as a paradox and resolved by the introduction of filtered/partitioned transfer spectra, which exploited the symmetries of the triadic interactions (J. Phys. A: Math. Theor., 2008). In this paper we consider the more general implications of that procedure for the spectral energy balance equation, also known as the Lin equation. It is argued that the resulting modified Lin equations (and their corresponding Navier–Stokes equations) offer a new starting point for both numerical and theoretical methods, which may lead to a better understanding of the underlying energy transfer processes in turbulence. In particular the filtered partitioned transfer spectra could provide a basis for a hybrid approach to the statistical closure problem, with the different spectra being tackled using different methods.     

Sharp vorticity gradients in two-dimensional turbulence and the energy spectrum

Formation of sharp vorticity gradients in two-dimensional (2D) hydrodynamic turbulence and their influence on the turbulent spectra are considered. The analog of the vortex line representation as a transformation to the curvilinear system of coordinates moving together with the di-vorticity lines is developed and compressibility of this mapping appears as the main reason for the formation of the sharp vorticity gradients at high Reynolds numbers. In the case of strong anisotropy the sharp vorticity gradients can generate spectra which fall off as k ^?3 at large k, which appear to take the same form as the Kraichnan spectrum for the enstrophy cascade. For turbulence with weak anisotropy the k dependence of the spectrum due to the sharp gradients coincides with the Saffman spectrum: E(k) ~ k ^?4. Numerical investigations of decaying turbulence reveal exponential growth of di-vorticity with a spatial distributed along straight lines. Thus, indicating strong anisotropy and accordingly the spectrum is close to the k ^?3-spectrum.     

Direct Evaluation of Cohesive Law in Mode I of Pinus pinaster by Digital Image Correlation

The direct identification of the cohesive law in pure mode I of Pinus pinaster is addressed in this work. The approach couples the double cantilever beam (DCB) test with digital image correlation (DIC). Wooden beam specimens loaded in the radial-longitudinal (RL) fracture propagation system are used. The strain energy release rate in mode I (G _I) is uniquely determined from the load–displacement curve by means of the compliance-based beam method (CBBM). This method relies on the concept of equivalent elastic crack length (a _eq) and therefore does not require the monitoring of crack propagation during test. DIC measurements are processed with two different purposes. Firstly, the physical evidence of a _eq is discussed with regard to actual estimation of the crack length based on post-processing full-field displacement measurements. Secondly, the crack tip opening displacement in mode I (w _I) is determined from the displacements near the initial crack tip. The cohesive law in mode I (σ _I???w _I) is then identified by numerical differentiation of the G _I???w _I relationship. The methodology and accuracy on this reconstruction are addressed. Moreover, the proposed procedure is validated by finite element analyses including cohesive zone modelling. It is concluded that the proposed data reduction scheme is adequate for assessing the cohesive law in pure mode I of P. pinaster.     

Diffusion to a chain of drops (bubbles) at large péclet numbers

The problem of diffusion of a substance, dissolved in a flow, to absorbing drops (bubbles) moving one after another in a viscous incompressible fluid is investigated. An approximate analytic expression is obtained for the differential and integral flows of the substance to the surface of each drop with consideration of the changes of the concentration and velocity fields due to the presence of other drops. A chain of spherical drops of equal radius arranged on the axis of a uniform forward flow is examined. It is shown that if the distance between drops, referred to the radius of the drops, satisfies the inequality 1lP^1/2 (P is the Péclet number), then the integral inflow of the substance to the surface of the second drop of the chain is 2.41 times less than the integral inflow to the first (the drops are enumerated along the flow); the total diffusion flow to the surface of an arbitrary drop with number k is determined by the expression I_k=I₁[k^1/2 – (k–1)^1/2], where I_k is the total flow to the first drop of the chain. The case of diffusion interaction of a solid particle and drop is examined. It is shown that for particles moving one after another with the same velocity in a fluid at rest the presence of a drop before the solid particle leads to a marked decrease of the total diffusion flow of the solid particle [by O(P^1/6) times], whereas the presence of a solid particle before a drop does not affect (in the main approximation with respect to the characteristic diffusion parameter) the total flow of the latter.I _k=I _i[k ^1/2–(k–1)^1/2]Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–69, January–February, 1978.     

Rhythmic precipitate patterns and fractal structure

Liesegang patterns of parallel precipitate bands are obtained when solutions containing co-precipitate ions interdiffuse in a 1D gel matrix.The sparingly soluble salt formed,displays a beautiful stratification of discs of precipitate perpendicular to the 1D tube axis.The Liesegang structures are analyzed from the viewpoint of their fractal nature.Geometric Liesegang patterns are constructed in conformity with the well-known empirical laws such as the time,band spacing and band width laws.The dependence of the band spacing on the initial concentrations of diffusing(outer)and immobile(inner)electrolytes(A0 and B0,respectively)is taken to follow the Matalon-Packter law.Both mathematical fractal dimensions and box-count dimensions are calculated.The fractal dimension is found to increase with increasing A0 and decreasing B0.We also analyze mosaic patterns with random distribution of crystallites,grown under different conditions than the classical Liesegang gel method,and report on their fractal properties.Finally,complex Liesegang patterns wherein the bands are grouped in multiplets are studied,and it is shown that the fractal nature increases with the multiplicity.     

Impact energy absorption by a rate-dependent locking target

The impact by an elastic cylindrical piston on a thin plate-like target resting on a rigid foundation is considered. The relationship between force F acting on the target and displacement x is given by F=kx+q dx/dt provided dx/dt≥0 and 0≤x<d (k, q and d≥0). When x=d locking occurs, and F can assume any value ≥kd without increase in x. The displacement is assumed to be completely irreversible. The motion of the impactor is assumed to be governed by the elementary wave equation and, since the target is thin, wave motion in the target is neglected. The energy W=εFdx and its components W _k=kεx dx (the energy absorbed in a corresponding quasistatic process) and W _q=qε(dx/dt)² dt (the excessive energy because of the rate-dependence) are determined in terms of the impact energy as functions of non-dimensional parameters representing k, q and d. With the aid of diagrams, it is shown under what circumstances locking occurs, and under what circumstances W _k or W _q, or both, are large.     

An improved anisotropy-resolving subgrid-scale model for flows in laminar–turbulent transition region

Some types of mixed subgrid-scale (SGS) models combining an isotropic eddy-viscosity model and a scale-similarity model can be used to effectively improve the accuracy of large eddy simulation (LES) in predicting wall turbulence. Abe (2013) has recently proposed a stabilized mixed model that maintains its computational stability through a unique procedure that prevents the energy transfer between the grid-scale (GS) and SGS components induced by the scale-similarity term. At the same time, since this model can successfully predict the anisotropy of the SGS stress, the predictive performance, particularly at coarse grid resolutions, is remarkably improved in comparison with other mixed models. However, since the stabilized anisotropy-resolving SGS model includes a transport equation of the SGS turbulence energy, k_SGS, containing a production term proportional to the square root of k_SGS, its applicability to flows with both laminar and turbulent regions is not so high. This is because such a production term causes k_SGS to self-reproduce. Consequently, the laminar–turbulent transition region predicted by this model depends on the inflow or initial condition of k_SGS. To resolve these issues, in the present study, the mixed-timescale (MTS) SGS model proposed by Inagaki et al. (2005) is introduced into the stabilized mixed model as the isotropic eddy-viscosity part and the production term in the k_SGS transport equation. In the MTS model, the SGS turbulence energy, k_es, estimated by filtering the instantaneous flow field is used. Since the k_es approaches zero by itself in the laminar flow region, the self-reproduction property brought about by using the conventional k_SGS transport equation model is eliminated in this modified model. Therefore, this modification is expected to enhance the applicability of the model to flows with both laminar and turbulent regions. The model performance is tested in plane channel flows with different Reynolds numbers and in a backward-facing step flow. The results demonstrate that the proposed model successfully predicts a parabolic velocity profile under laminar flow conditions and reduces the dependence on the grid resolution to the same degree as the unmodified model by Abe (2013) for turbulent flow conditions. Moreover, it is shown that the present model is effective at transitional Reynolds numbers. Furthermore, the present model successfully provides accurate results for the backward-facing step flow with various grid resolutions. Thus, the proposed model is considered to be a refined anisotropy-resolving SGS model applicable to laminar, transitional, and turbulent flows.     

Simulating the formation of spiral wave in the neuronal system

Some experimental results have confirmed that a spiral wave could be observed in the cortex of brain. The biological Hodgkin–Huxley neurons are used to construct a regular network with nearest-neighbor connection, artificial line defects are generated to block the traveling wave in the network, and the potential mechanism for formation of spiral wave is investigated. A target wave is generated in a local area by imposing two external forcing currents with diversity (I ₀?I ₁) in different areas of the network. It is confirmed that spiral wave could be induced by the defects even if no specific initial values are used. A single perfect spiral wave can occupy the network when the coupling intensity exceeds certain threshold; otherwise, a group of spiral waves emerges in the network. Certain channel noise can enhance the diversity (I ₀?I ₁) for generating target wave, and then spiral waves are induced by blocking the target wave with defects under no-flux and/or periodic boundary conditions in the network.     

On the nonlinear slewing dynamics and control of the Space Station based Mobile Servicing System

A relatively general Lagrangian formulation for studying the nonlinear dynamics and control of space-craft with interconnected flexible members in a tree-type topology is developed. Versatility of the formulation is illustrated through a dynamical study of the Space Station based two-link Mobile Servicing System (MSS). The performance of the MSS undergoing inplane and out-of-plane slewing maneuvers is compared. Results indicate that, in absence of control, the maneuvers induce undesirable librational motion of the Space Station as well as vibration of the links. Nonlinear control, based on the Feedback Linearization Technique (FLT), appears promising. Quasi-Closed Loop Control (QCLC), a variation of the FLT, is applied to control the libration of the Space Station. Once the attitude of the Space Station is controlled, the performance of the MSS improves significantly. For a 5-minute maneuver of the MSS, the maximum control torque required is only 34.5 Nm.Nomenclature f _i ¹ , f _i,j ¹ fundamental frequency of bodies B _i and B _i,j, respectively - l _c, l _i, l _i,j length of bodies B _c, B _i, and B _i,j, respectively - m _c, m _i, m _i,j mass of bodies B _c, B _i, and B _i,j, respectively - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaerbhv2BYDwAHbacfiGab8xCayaaraqefavySfgDP52BGWuAU9gD% 5bxzaGGbciaa+zgacaWFSaGaa8hiaiqa-fhagaqeaiaa-jhaaaa!4B1F!\[\bar qf, \bar qr\] vector representing flexible and rigid generalized coordinates - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaerbhv2BYDwAHbacfiGaa8hkaiqa-fhagaqeaiaa-jhacaWFPaqe% favySfgDP52BGWuAU9gD5bxzaGGbciaa+rgaaaa!4A18!\[(\bar qr)d\] vector representing the desired rigid generalized coordinates - (I _xx)k, (I _yy)k, (I _zz)k principal inertia of body B _k about X _k, Y _k, and Z _k axes, respectively; ksc, i or i, j - K _p, K _v displacement and velocity gain matrices - N _q total number of generalized coordinates - % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaerbwvMCKfMBHbacfiGab8xuayaaraqefavySfgDP52BGWuAU9gD% 5bxzaGGbciaa+zgaieaacaqFSaGaa0hiaiqa-ffagaqeaGqaciaa8j% haaaa!4AEF!\[\bar Qf, \bar Qr\] control effort vectors for flexible and rigid coordinates, respectively - Q , Q , Q control effort for pitch, roll and yaw degree of freedom, respectively - _k ^y , _k ^z tip deflection of a beam type appendage (B _k) in the Y _k and Z _k directions, respectively.     

Drag of a sphere in an unsteady flow of viscous fluid at small reynolds numbers

An unsteady flow of viscous incompressible fluid past a sphere is investigated. The values of the inertial and unsteady terms in the Navier-Stokes equations are characterized by translational (R) and vibrational (R_k) Reynolds numbers, which are assumed small. The solution is constructed in the form of an expansion with respect to max(R, R _k ^1/2 ) by the method of matched asymptotic expansions. A correction to the Stokes force, correct to o[max(R, R _k ^1/2 )], is calculated. It is shown that the result depends strongly on the ratio R/R _k ^1/2 and goes over into the well-known equations for the cases R 0, R_k 0.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 11–16, January–February, 1988.     

Vorticity dynamics in turbulence growth

The mechanisms of vorticity amplification in the formation of turbulence are investigated by means of direct numerical simulations of the Navier–Stokes equations with different initial conditions and Reynolds numbers. The simulations show good universality of the enstrophy evolution, that occurs in two stages. The first stage is dominated by the effect of vortex stretching, and it finishes with a k ^?3 power-law energy spectrum. The second stage is dominated by the action of viscosity on the small scales, and it finishes with a Kolmogorov k ^?5/3 energy spectrum.     

DNS and LES of turbulent flow in a closed channel featuring a pattern of hemispherical roughness elements

Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) were performed for fully-developed turbulent flow in channels with smooth walls and walls featuring hemispherical roughness elements at shear Reynolds numbers Re_τ = 180 and 400, with the goal of studying the effect of these roughness elements on the wall-layer structure and on the friction factor. The LES and DNS approaches were verified first by comparison with existing DNS databases for smooth walls. Then, a parametric study for the hemispherical roughness elements was conducted, including the effects of shear Reynolds number, normalized roughness height (k⁺ = 10–20) and relative roughness spacing (s⁺/k⁺ = 2–6). The sensitivity study also included the effect of distribution pattern (regular square lattice vs. random pattern) of the roughness elements on the walls. The hemispherical roughness elements generate turbulence, thus increasing the friction factor with respect to the smooth-wall case, and causing a downward shift in the mean velocity profiles. The simulations revealed that the friction factor decreases with increasing Reynolds number and roughness spacing, and increases strongly with increasing roughness height. The effect of random element distribution on friction factor and mean velocities is however weak. In all cases, there is a clear cut between the inner layer near the wall, which is affected by the presence of the roughness elements, and the outer layer, which remains relatively unaffected. The study reveals that the presence of roughness elements of this shape promotes locally the instantaneous flow motion in the lateral direction in the wall layer, causing a transfer of energy from the streamwise Reynolds stress to the lateral component. The study indicates also that the coherent structures developing in the wall layer are rather similar to the smooth case but are lifted up by almost a constant wall-unit shift y⁺ (∼10–15), which, interestingly, corresponds to the relative roughness k⁺ = 10.     

The strain energy criterion of mixed-mode brittle fracture

In this paper, we proposed a new cr tenon or mixea-moae brittle fracture, i.e., the strain energy criterion, which can be stated as (K_Ⅰ/K_Ⅰc)² +(K_Ⅰ/K_Ⅰc)²+(K_Ⅱ/K_Ⅱc)²=1. This criterion is shown to be in good agreement with known experimental data.In this paper, an experimental criterion:(K_Ⅰ/K_Ⅰc)^m+(K_Ⅱ/K_Ⅱc)ⁿ=1, 1≤_n^m≤2.is also proposed.     

Instability of small-amplitude convective flows in a rotating layer with free boundaries

The stability of steady convective flows in a horizontal layer with free boundaries, heated from below and rotating about a vertical axis, is studied in the Boussinesq approximation (Rayleigh-Bénard convection). The flows considered are convective rolls or square cells that are sums of two perpendicular rolls with equal wave numbers k. It is assumed that the Rayleigh number is almost critical in order for convective flows with a wave number k: R = R _c (k) + ε² to arise, the amplitude of the supercritical states being of the order of ε. It is shown that the flows are always unstable relative to perturbations that are the sum of one long-and two short-wave modes corresponding to linear rolls turned through small angles in opposite directions.     

Plasma drift-anisotropy instability and anomalous transport processes

The author expounds a nonlinear theory of the instability of a weakly inhomogeneous plasma with hot ions when a”loss cone” is present in their velocity distribution. Flute-type instabilities (k_Z= 0) are considered, which for a strong enough irregularity may build up even in short traps with”magnetic mirrors.” It is shown that the total ion flux through the magnetic mirrors, which is caused by turbulent diffusion into the loss cone, exceeds by a factor of (2n/R_H ?n)^1/2 the ion flux across the magnetic field as a result of diffusion in coordinate space (here n, ?n, are the density and its gradient, R_H is the ion Larmor radius). The diffusion time of ions into the “loss cone” is of the order τ=10Ω_? ^? (2n/R ^??n)^3/2 (Ω_H is the ion Larmor frequency). A plasma contained in magnetic traps is always in a nonequilibrium thermodynamic state. The nature of the nonequilibrium is connected with the specific geometry of the containing magnetic field. Here we will consider open traps with magnetic mirrors in which the nonequilibrium of the plasma is caused by:

1. 1)
   the curvature of the magnetic field force lines and its associated effective gravity field;     
   
   

Slip motion and stability of a single degree of freedom elastic system with rate and state dependent friction

We consider quasistatic motion and stability of a single degree of freedom elastic system undergoing frictional slip. The system is represented by a block (slider) slipping at speed V and connected by a spring of stiffness k to a point at which motion is enforced at speed V₀ We adopt rate and state dependent frictional constitutive relations for the slider which describe approximately experimental results of Dieterich and Ruina over a range of slip speeds V. In the simplest relation the friction stress depends additively on a term A In V and a state variable θ; the state variable θ evolves, with a characteristic slip distance, to the value ? B In V, where the constants A, B are assumed to satisfy B > A > 0. Limited results are presented based on a similar friction law using two state variables.Linearized stability analysis predicts constant slip rate motion at V₀ to change from stable to unstable with a decrease in the spring stiffness k below a critical value k_cr. At neutral stability oscillations in slip rate are predicted. A nonlinear analysis of slip motions given here uses the Hopf bifurcation technique, direct determination of phase plane trajectories, Liapunov methods and numerical integration of the equations of motion. Small but finite amplitude limit cycles exist for one value of k, if one state variable is used. With two state variables oscillations exist for a small range of k which undergo period doubling and then lead to apparently chaotic motions as k is decreased.Perturbations from steady sliding are imposed by step changes in the imposed load point motion. Three cases are considered: (1) the load point speed V₀ is suddenly increased; (2) the load point is stopped for some time and then moved again at a constant rate; and (3) the load point displacement suddenly jumps and then stops. In all cases, for all values of k:, sufficiently large perturbations lead to instability. Primary conclusions are: (1) ‘stick-slip’ instability is possible in systems for which steady sliding is stable, and (2) physical manifestation of quasistatic oscillations is sensitive to material properties, stiffness, and the nature and magnitude of load perturbations.     

Measurement of liquid film thickness in a micro parallel channel with interferometer and laser focus displacement meter

In the present study, liquid film thicknesses in parallel channels with heights of H = 0.1, 0.3 and 0.5 mm are measured with two different optical methods, i.e., interferometer and laser focus displacement meter. Ethanol is used as a working fluid. Liquid film thicknesses obtained from two optical methods agree very well. At low capillary numbers, dimensionless liquid film thickness is in accordance with Taylor’s law. However, as capillary number increases, dimensionless liquid film thickness becomes larger than Taylor’s law for larger channel heights. It is attributed to the dominant inertial effect at high capillary numbers. Using channel height H for dimensionless liquid film thickness δ₀/H and hydraulic diameter D_h = 2H as the characteristic length for Reynolds and Weber numbers, liquid film thickness in a parallel channel can be predicted well by the circular tube correlation previously proposed by the authors. This is because curvature differences between bubble nose and flat film region are identical in circular tubes and parallel channels.     

Numerical solution of the turbulent flow of an incompressible fluid in curved planar and axially symmetrical channels

The governing equations for axially symmetric flow, where the Reynolds stresses are expressed by scalar turbulent viscosity, are the Reynolds equations. The turbulence model k, ? is used in the well-known form for fully developed turbulent flow.The numerical method, a continuation of the MAC system¹, is adapted so that even for high Reynolds cell numbers precision (δx²) can be achieved for the steady flow. Irregular cells join the rectangular network on the curved surface. Von Neumann's stability condition of the linearised numerical system is investigated. Special problems concerning the numerical solution of the turbulence model equations are stated and a special procedure is worked out to ensure that the fields k, ? do not converge to physically meaningless values. The program for the computer is universal in that the boundary problems can be assigned by input data.As an example, an axially symmetrical diffuser with an area ratio of widening 1.40 is computed. Fields of velocity and pressure at the wall as well as fields v_T and k are assessed. The results are compared with an experiment. The conclusion is that this method is suitable for the problems mentioned in this study as well as for unsteady flow.